The spatiotemporal variability of latent heat flux (LE) and water vapor mixing ratio (rv) variability are not well understood due to the scale-dependent and nonlinear atmospheric energy balance responses to land surface heterogeneity. Airborne in situ and profiling Raman lidar measurements with the wavelet technique are utilized to investigate scale-dependent relationships among LE, vertical velocity (w) variance (s2w), and rv variance (s2wv) over a heterogeneous surface in the Chequamegon Heterogeneous Ecosystem Energy-balance Study Enabled by a High-density Extensive Array of Detectors 2019 (CHEESEHEAD19) field campaign. Our findings reveal distinct scale distributions of LE, s2w, and s2wv at 100 m height, with a majority scale range of 120m-4km in LE, 32m-2km in s2w, and 200 m – 8 km in s2wv. The scales are classified into three scale ranges, the turbulent scale (8m–200m), large-eddy scale (200m–2km), and mesoscale (2 km–8km) to evaluate scale-resolved LE contributed by s2w and s2wv. In the large-eddy scale in Planetary Boundary Layer (PBL), 69-75% of total LE comes from 31-51% of the total sw and 39-59% of the total s2wv. Variations exist in LE, s2w, and s2wv, with a range of 1.7-11.1% of total values in monthly-mean variation, and 0.6–7.8% of total values in flight legs from July to September. These results confirm the dominant role of the large-eddy scale in the PBL in the vertical moisture transport from the surface to the PBL. This analysis complements published scale-dependent LE variations, which lack detailed scale-dependent vertical velocity and moisture information.

The Chequamegon Heterogeneous Ecosystem Energy-balance Study Enabled by a High-density Extensive Array of Detectors 2019 (CHEESEHEAD19) is an ongoing National Science Foundation project based on an intensive field campaign that occurred from June-October 2019. The purpose of the study is to examine how the atmospheric boundary layer responds to spatial heterogeneity in surface energy fluxes. One of the main objectives is to test whether lack of energy balance closure measured by eddy covariance (EC) towers is related to mesoscale atmospheric processes. Finally, the project evaluates data-driven methods for scaling surface energy fluxes, with the aim to improve model-data comparison and integration. To address these questions, an extensive suite of ground, tower, profiling, and airborne instrumentation was deployed over a 10×10 km domain of a heterogeneous forest ecosystem in the Chequamegon-Nicolet National Forest in northern Wisconsin USA, centered on the existing Park Falls 447-m tower that anchors an Ameriflux/NOAA supersite (US-PFa / WLEF). The project deployed one of the world’s highest-density networks of above-canopy EC measurements of surface energy fluxes. This tower EC network was coupled with spatial measurements of EC fluxes from aircraft, maps of leaf and canopy properties derived from airborne spectroscopy, ground-based measurements of plant productivity, phenology, and physiology, and atmospheric profiles of wind, water vapor, and temperature using radar, sodar, lidar, microwave radiometers, infrared interferometers, and radiosondes. These observations are being used with large eddy simulation and scaling experiments to better understand sub-mesoscale processes and improve formulations of sub-grid scale processes in numerical weather and climate models.