Publication

It is frequently observed in field experiments that the eddy covariance heat fluxes are systematically underestimated as compared to the available energy. The flux imbalance problem is investigated using the NCAR's large-eddy simulation (LES) model imbedded with an online scheme to calculate Reynolds-averaged fluxes. A top-down and a bottom-up tracer are implemented into the LES model to quantify the influence of entrainment and bottom-up diffusion processes on flux imbalance. The results show that the flux imbalance follows a set of universal functions that capture the exponential decreasing dependence on u(*)/w(*), where u(*) and w(*) are friction velocity and the convective velocity scale, respectively, and an elliptic relationship to z/z(i) , where z(i) is the mixing-layer height. The source location in the boundary layer is an important factor controlling the imbalance magnitude and its horizontal and vertical distributions. The flux imbalance of heat and the bottom-up tracer is tightly related to turbulent coherent structures, whereas for the top-down diffusion, such relations are weak to nonexistent. Our results are broadly consistent with previous studies on the flux imbalance problem, suggesting that the published results are robust and are not artefacts of numerical schemes.