A Lagrangian dispersion model for predicting CO2sources, sinks, and fluxes in a uniform loblolly pine (Pinus taeda L.) stand

A canopy Lagrangian turbulent scalar transport model for predicting scalar fluxes, sources, and sinks within a forested canopy was tested using CO2concentration and flux measurements. The model formulation is based on the localized near‐field theory (LNF) proposed byRaupach[1989a, b]. Using the measured mean CO2concentration profile, the vertical velocity variance profile, and the Lagrangian integral timescale profile within and above a forested canopy, the proposed model predicted the CO2flux and source (or sink) profiles. The model testing was carried out using eddy correlation measurements at 9 m in a uniform 13 m tallPinus taeda L.(loblolly pine) stand at the Blackwood division of the Duke Forest near Durham, North Carolina. The tree height and spacing are relatively uniform throughout. The measured vertical profile leaf area index (LAI) was characterized by three peaks, with a maximum LAI occurring at 6.5 m, in qualitative agreement with the LNF source‐sink predicted profile. The LNF CO2flux predictions were in better agreement with eddy correlation measurements (coefficient of determinationr2= 0.58; and standard error of estimate equal to 0.16 mg kg−1m s−1) thanKtheory. The model reproduced the mean diurnal CO2flux, suggesting better performance over longer averaging time periods. Two key simplifications to the LNF formulation were considered, namely, the near‐Gaussian approximation to the vertical velocity and the absence of longitudinal advection. It was found that both of these assumptions were violated throughout the day, but the resulting CO2flux error at 9 m was not strongly related to these approximations. In contrast to the forward LNF approach utilized by other studies, this investigation demonstrated that the inverse LNF approach is sensitive to near‐field corrections.