The observing system design of multi-disciplinary field measurements involves a variety of considerations on logistics, safety, and science objectives. Typically, this is done based on investigator intuition and designs of prior field measurements. However, there is potential for considerable increase in efficiency, safety, and scientific success by integrating numerical simulations in the design process. Here, we present a novel approach to observing system simulation experiments that aids surface-atmosphere synthesis at the interface of meso- and microscale meteorology. We used this approach to optimize the Chequamegon Heterogeneous Ecosystem Energy-balance Study Enabled by a High-density Extensive Array of Detectors 2019 (CHEESEHEAD19). During pre-field simulation experiments, we considered the placement of 20 eddy-covariance flux towers, operations for 72 hours of low-altitude flux aircraft measurements, and integration of various remote sensing data products. High-resolution Large Eddy Simulations generated a super-sample of virtual ground, airborne, and satellite observations to explore two specific design hypotheses. We then analyzed these virtual observations through Environmental Response Functions to yield an optimal aircraft flight strategy for augmenting a stratified random flux tower network in combination with satellite retrievals. We demonstrate how this novel approach doubled CHEESEHEAD19’s ability to explore energy balance closure and spatial patterning science objectives while substantially simplifying logistics. Owing to its extensibility, the approach lends itself to optimize observing system designs also for natural climate solutions, emission inventory validation, urban air quality, industry leak detection and multi-species applications, among other use cases.

The Earth’s surface is heterogeneous at multiple scales owing to spatial variability in various properties. The atmospheric responses to these heterogeneities through fluxes of energy, water, carbon and other scalars are scale-dependent and non-linear. Although these exchanges can be measured using the eddy covariance technique, widely used tower-based measurement approaches suffer from spectral losses in lower frequencies when using typical averaging times. However, spatially resolved measurements such as airborne eddy covariance measurements can detect such larger scale (meso-{$\beta$}, $\gamma$) transport. To evaluate the prevalence and magnitude of these flux contributions we applied wavelet analysis to airborne flux measurements over a heterogeneous mid-latitude forested landscape, interspersed with open water bodies and wetlands. The measurements were made during the Chequamegon Heterogeneous Ecosystem Energy-balance Study Enabled by a High-density Extensive Array of Detectors (CHEESEHEAD19) intensive field campaign. We ask, how do spatial scales of surface-atmosphere fluxes vary over heterogeneous surfaces across the day and across seasons? Measured fluxes were separated into smaller-scale turbulent and larger-scale mesoscale contributions. We found significant mesoscale contributions to H and LE fluxes through summer to autumn which wouldn’t be resolved in single point tower measurements through traditional time-domain half-hourly Reynolds decomposition. We report scale-resolved flux transitions associated with seasonal and diurnal changes of the heterogeneous study domain. This study adds to our understanding of surface atmospheric interactions over unstructured heterogeneities and can help inform multi-scale model-data integration of weather and climate models at a sub-grid scale.

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.