While yearly budgets of CO2 and evapotranspiration (ET) above forests can be readily obtained from eddy-covariance measurements, the quantification of their respective soil (respiration and evaporation) and canopy (photosynthesis and transpiration) components remains an elusive yet critical research objective. To this end, methods capable of reliably partitioning the measured ET and F_c fluxes into their respective soil and plant sources and sinks are highly valuable. In this work, we investigate four partitioning methods (two new, and two existing) that are based on analysis of conventional high frequency eddy-covariance (EC) data. The physical validity of the assumptions of all four methods, as well as their performance under different scenarios, are tested with the aid of large eddy simulations, which are used to replicate eddy-covariance field experiments. Our results indicate that canopies with large, exposed soil patches increase the mixing and correlation of scalars; this negatively impacts the performance of the partitioning methods, all of which require some degree of uncorrelatedness between CO2 and water vapor. In addition, best performance for all partitioning methods were found when all four flux components are non-negligible, and measurements are collected close to the canopy top. Methods relying on the water-use efficiency (W) perform better when W is known a priori, but are shown to be very sensitive to uncertainties in this input variable especially when canopy fluxes dominate. We conclude by showing how the correlation coefficient between CO2 and water vapor can be used to infer the reliability of different W parameterizations.

Unsteady land-sea breezes (LSBs) resulting from time-varying surface thermal contrasts Δθ(t) are explored in the presence of a constant synoptic pressure forcing, Mg, when the latter is oriented from sea to land (α=0°), versus land to sea (α=180°). Large eddy simulations reveal the development of four distinctive regimes depending on the joint interaction between (Mg, α) and Δθ(t) in modulating the fine-scale dynamics. Time lags, computed as the shifts that maximize correlation coefficients of the dynamics between transient and the corresponding steady state scenarios at Δθ=Δθmax, are found to be significant and to extend 2 hours longer for α=0° compared to α=180°. These diurnal dynamics result in non-equilibrium flows that behave differently over the two patches for both α’s. Turbulence is found to be out of equilibrium with the mean flow, and the mean itself is found to be out of equilibrium with the thermal forcing. The sea surface heat flux is consistently more sensitive than its land counterpart to the time-varying external forcing Δθ(t), and more so for synoptic forcing from land-to-sea (α=180°). Hence, although the land reaches equilibrium faster, the sea patch is found to exert a stronger control on the final turbulence-mean flow equilibrium response. Finally, vertical velocity profile at the shore and shore-normal velocity transects at the first grid level are shown to encode the multiscale regimes of the LSBs evolution, and can thus be used to identify these regimes using k-means clustering.