Diagnosing Land-Atmosphere Coupling Errors in E3SM v3 Using Diurnal Metrics and ARM Observations
Explore the source record for details and available documents.
SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
This dataset contains land model benchmarking results for the Energy Exascale Earth System Model version 2.1 (E3SMv2.1), including outputs from both coupled biogeochemistry simulations and stand-alone land model simulations. These results are compared against several emission-driven historical simulations from the Coupled Model Intercomparison Project Phase 6 (CMIP6). Benchmarking was conducted using the International Land Model Benchmarking (ILAMB) package, version 2.7 (ILAMBv2.7). CMIP6 model outputs were sourced from the Earth System Grid Federation (ESGF), while the E3SMv2.1 results were derived from raw model outputs. These outputs underwent processing steps such as time serialization, conservative regridding, and data standardization to ensure comparability. For spatial interpolation, the Earth System Modeling Framework (ESMF) tool, ESMF_RegridWeightGen, was employed to generate regridding weights, enabling the transformation of E3SM’s native cubed-sphere grid to a regular latitude-longitude grid.
Emissions-driven (prognostic CO 2 ) simulations are essential for representing two-way carbon-climate feedback in Earth System Models. We present an emissions-driven land–atmosphere coupled biogeochemistry (BGC) configuration (BGCLNDATM_progCO2) in version 2.1 of the Energy Exascale Earth System Model (E3SMv2.1). This is the first E3SM configuration that performs land-atmosphere emission-hindcasts. Here, we document its implementation, evaluate the model's performance against observations and other models, and propose a structured evaluation protocol for such emissions-driven simulations. We conducted transient historical simulations (1850–2014) with BGCLNDATM_progCO2 and compare them to reference simulations—a land-atmosphere coupled simulation without BGC and a standalone land simulation with BGC, both using prescribed CO 2 concentrations—and to observations. BGCLNDATM_progCO2 overestimates atmospheric CO 2 concentrations by 11–23 ppm yet stays within the 40-ppm spread CMIP6 emission-driven models and retains physical climate properties comparable to the reference runs. The CO 2 biases are partly attributed to underrepresented oceanic CO 2 uptake and inadequate representations of some terrestrial processes. In general, introducing prognostic CO 2 did not change physical climate metrics at the global scale but had larger regional effects, particularly over land where spatially heterogeneous CO 2 and prognostic leaf area index influenced surface energy balance. Finally, we propose a general evaluation protocol including spin-up assessment, atmospheric CO 2 benchmarking, physical climate evaluation, and land biogeochemical analysis to support scientific rigor and facilitate inter-model comparisons. The new configuration lays the groundwork for future enhancements, including improved terrestrial biogeochemical processes, integrated marine biogeochemistry, and additional human–Earth system interactions. These developments advance E3SM toward fully coupled emissions-driven simulations, enabling more accurate carbon–climate feedback projections and informing mitigation policy by providing physically consistent carbon-budget metrics for mitigation scenarios.
Accurate representation of groundwater table depth (GWTD) is crucial for simulating hydrological cycling in Earth system models (ESM). Nevertheless, there is a notable gap in the literature regarding the validation of GWTD simulations in ESMs and their subsequent impact on downstream hydrological components. This study explores the calibration of parameterization of global GWTD using machine learning within the Energy Exascale Earth System Model (E3SM) Land Model (ELM). Despite achieving significant gains in simulating GWTD through calibration, offline ELM simulations unexpectedly show that these improvements do not translate to substantial enhancements in model performance for other key hydrological variables, including soil moisture (SM), runoff, groundwater contribution to runoff or base flow index (BFI), and evapotranspiration and its partitioning. The performance in SM and runoff was even degraded in some regions, while BFI was mostly overestimated. Although there is significant improvement in GWTD within the critical range of 1–5 m, where groundwater traditionally influences land surface energy fluxes, these improvements occurred mostly in humid areas where the impact of GWTD on surface processes is minimal. Although the impacts of model calibration are generally small in offline ELM simulations, coupled land-atmosphere simulations exhibit much stronger responses to GWTD calibration, highlighting the role of land-atmosphere feedbacks in Earth system modeling. These findings underscore the need for integrated calibration strategies that simultaneously optimize multiple hydrological variables. However, if a single-variable approach is necessary, it is crucial to establish clear priorities for calibration, identifying the most critical variables that have the greatest impact on overall model performance.
Tropical land generally experiences the hottest period (spring) in a year just before the onset of wet season. Previous studies suggested that in a warming climate, the wet season would come later, but its origin is debated and its impact on temperature remains unknown. Here, we find that the warming of hot season would be amplified under global warming, and refer to it as “hot-season-gets-hotter” phenomenon. The amplified hot season warming is closely tied to the amplified warming of hot temperature percentiles. The hot-season-gets-hotter phenomenon is mainly due to the rainfall delay and most evident in the Amazon, where spring is warming by almost 1?K more than the annual mean and the 99th percentile temperatures are warming ~30% more than the mean by the end of 21st century in a high emission scenario. Comparing experiments with and without land-atmosphere coupling, it is further found that the rainfall delay is initially driven by the enhanced effective atmospheric heat capacity and then substantially amplified by positive soil moisture-atmosphere feedback. In the satellite period, observations consistently show that the hot-season-gets-hotter phenomenon has already emerged along with the rainfall delay in the Amazon. Intensified hot and dry spring climate can enhance risks of drought, heatwaves and wildfires, threatening the Amazon forest and habitats in the tropics.
Multiple land-atmosphere coupling processes influence the spatial distribution of boundary-layer energy, moisture, aerosol, or cloud condensation nuclei (CCN) properties. Accurately quantifying the spatial and temporal distributions of CCN concentrations is essential for investigating aerosol-cloud interactions and reducing uncertainties in aerosol radiative forcing estimates of Earth's system.
The U.S. Department of Energy’s Atmospheric Radiation Measurement (ARM) User Facility supports atmospheric science through an integrated network of fixed and mobile observatories. These facilities collect continuous and campaign-based observations of atmospheric properties, with the goal of improving the representation of clouds, aerosols, precipitation, and radiation in Earth system models. The Bankhead National Forest (BNF) site, established as an ARM Mobile Facility (AMF) on 1 October 2024, is situated in a forested region of northern Alabama. Its strategic location in a southeastern U.S. environment characterized by complex terrain, diverse land cover, and frequent convective storms provides a valuable opportunity to examine coupled land-atmosphere processes under natural variability.
Assessing energy resources under future scenarios requires high-resolution meteorological information that is physically consistent and suitable for regional-scale analysis. While Earth system model (ESM) projections provide valuable large-scale information, their coarse resolution and systematic biases limit direct applicability for energy system modeling and planning. In this study, we develop a high-resolution dynamical downscaling framework based on the Weather Research and Forecasting (WRF) model to translate global-scale ESM data into energy-relevant regional projections over the contiguous United States (CONUS). The framework identifies an optimized WRF configuration through numerical experiments and evaluates raw and bias-corrected ESM initial and boundary conditions, with soil moisture (SM) and soil temperature (ST) bias correction implemented as an integral part of the bias-corrected ESM forcing to improve land-atmosphere coupling prior to WRF dynamical downscaling. Using an optimized WRF configuration at 4-km resolution, we show that raw ESM forcing introduces systematic dry and cold soil biases that propagate into pronounced warm biases in near-surface air temperature and positive biases in solar irradiance, particularly during summer. Applying bias-corrected atmospheric forcing together with bias-corrected SM and ST substantially reduces these downstream biases and improves the surface energy balance and near-surface atmospheric fields. These results demonstrate that bias-aware treatment of initial conditions is critical for producing high-resolution downscaled projections suitable for energy system modeling and planning applications.
The accurate representation of interactions between clouds and planetary boundary layer (PBL) is a persistent challenge in climate models, critical for simulating surface energy budget. The emergence of kilometer-grid-scale global storm resolving models (GSRMs) offers the potential for enhanced details of PBL processes in these complex interactions. This study evaluates the representation of PBL-coupled and decoupled clouds in nine GSRM simulations against extensive ground-based observations by the Department of Energy Atmospheric Radiation Measurement (ARM) program, across six sites encompassing diverse regimes such as marine and continental environments in tropics and midlatitude. By differentiating coupling based on the relative positions between cloud bases and PBL tops, our analysis focuses on the simulation of PBL height, cloud frequency, position and vertical extent. The GSRMs generally exhibit commendable agreements with observed cloud structures and PBL diurnal cycles across different ARM sites. In contrast to the relatively consistent representation of decoupled clouds, discrepancies exist between the simulated and the observed coupled clouds, particularly in areas of intense convection, for example, over tropical rainforests and mountainous regions. These biases are probably associated with the models' tendency to underestimate the boundary layer humidity and the frequency of coupled clouds within different ranges of PBL heights. This study underscores the importance for continuous improvements in the representation of boundary layer and convection within these global kilometer-grid-scale models.
The generation of atmospheric aerosols is a fundamental process in aerosol-climate interactions, where particles act as nuclei for cloud formation and ultimately influence Earth’s radiative balance. Understanding aerosol precursors, chemical compositions, and formation mechanisms, particularly for secondary organic aerosols (SOAs), is essential for assessing their overall impact on climate. Aerosol-climate interactions are especially important in the Southeastern United States (SE-US), where the atmosphere is strongly influenced by both biogenic and anthropogenic emissions, as well as intense land-atmosphere convective coupling. Current mechanistic representations of SOA formation under ambient conditions are largely extrapolated from laboratory-scale oxidation experiments involving individual volatile organic compounds (VOCs), which are often assumed to behave independently when present in mixtures. However, recent empirical evidence demonstrates that coexisting VOC precursors can undergo previously unexplored molecular-scale interactions that significantly alter SOA physicochemical properties, including mass loading, yield, and water uptake.
Land-atmosphere interactions play a critical role in the evolution and formation of low-level clouds. The different states of coupling between low-level clouds and the surface are uncertain, primarily over continental regions, where complex thermodynamics complicates their investigation. This study uses observations from the Atmospheric Radiation Measurement User Facility to explore cloud-surface coupling and perform a climatological analysis of this interaction in five countries across three continents. The results reveal consistent coupling thresholds and average percentages across the five sites, with coupled clouds accounting for 66% of the cases and decoupled clouds for 34%. Thermodynamic and dynamic evaluations show distinct differences between coupled and decoupled clouds. Coupled clouds are characterized by humid environments, in which vertical motions connect the surface and lower atmosphere to the cloud base, conditions that favor the formation of boundary layer clouds. Decoupled clouds prefer to occur in a drier and colder environment with vertical motions inside the boundary layer being detached from the cloud base, under which boundary layer clouds are hard to form. Coupled clouds peak during warmer hours and seasons, and vice versa for decoupled clouds. This study underscores the complexity of cloud-land-surface interactions and paves the way for further investigations into cloud formation and evolution under different atmospheric environments.
The central United States frequently exhibits warm and dry biases in simulations of summertime conditions, a persistent feature that remains unresolved. While previous studies linked these biases to misrepresented surface energy exchanges, the role of belowground processes remains poorly understood. Here, we demonstrate that inadequate representation of root water uptake in land surface models contributes to this bias. Using both offline Noah-MP and coupled WRF-Noah-MP simulations with static and dynamic root water uptake schemes, we show that the inclusion of dynamic root processes reduces 2-m air temperature biases and enhances precipitation, primarily by increasing the rain rate of convective systems. Offline and coupled simulations further reveal that the cooling effects and precipitation increases are amplified through positive land-atmosphere feedback, active only in the coupled model. These findings highlight an important role of root in modulating land-atmosphere interactions and underscore the need to refine root-zone processes to improve regional atmospheric simulations.
Accurate representation of urban properties and processes at higher resolutions in global modeling systems is essential for advancing our ability to capture the complexities of urban systems and informing effective resilience strategies. However, the prescription of coarse global-scale urban properties in most state-of-the-art Earth system models (ESMs) is limiting their potential for capturing urban signals as they advance toward kilometer-scale simulation capabilities. To bridge this gap in inadequate urban property representation and to advance urban-resolving Earth system modeling, this work integrates the newly-developed global 1 km-resolution facet-level urban surface property data set, U-Surf, into the land component of Community Earth System Model (CESM)—Community Terrestrial System Model (CTSM). The land-only CTSM simulations are validated against satellite measurements, ground-based urban weather stations, flux tower observations, and reanalysis data. Results demonstrate that the enhanced urban properties allow improved simulations of urban meteorology and surface energy fluxes compared to the default coarse-resolution categorical urban canopy parameters. Spatial scaling analysis reveals regime-dependent information loss during resolution aggregation, as well as substantial scale-dependent variations in urban surface energy flux representation. Furthermore, these findings have critical implications for coupled Earth system modeling when including the effect of land-atmosphere interaction. This work establishes a foundation for future urban-resolving kilometer-scale ESM development, which will enable systematic intra- and inter-city comparisons that inform urban adaptation strategies across diverse global urban environments.
Soil moisture (SM) is a key regulator of ecosystem biogeophysics, influencing plant water relations and land-atmosphere energy exchanges. We evaluate the representation of SM in 16 Earth System Models from the Coupled Model Intercomparison Project Phase 6 (CMIP6) using the International Land Model Benchmarking (ILAMB) framework, focusing on surface (0–5, 0–10 cm) and rootzone (0–100 cm) depths, as well as key ecohydrological variables like gross primary productivity (GPP), leaf area index (LAI), and evapotranspiration (ET), and their coupling. Models are benchmarked against multiple observational and assimilated datasets to assess both state variables and cross-variable relationships. Surface SM is generally well represented (r > 0.87), while rootzone SM variability is systematically overestimated (normalized standard deviation > 1). ET shows strong agreement with observations (r > 0.9), whereas GPP and LAI exhibit larger inter-model spread. Skill in individual variables does not guarantee realistic SM–ecohydrology coupling, which varies strongly across models and depends on the reference dataset. Köppen-based regional analyses reveal strong regime dependence, with several models performing well in Tropical and Temperate regions but degrading in Continental (high-latitude) zones. Across both global and regional benchmarks, models cluster by land surface framework, indicating that structural choices in soil hydrology and soil–plant coupling exert a first-order control on performance. These results provide process-relevant benchmarks and suggest that improving the representation of vertical soil structure, rooting depth distributions, and soil–plant hydraulic coupling will be central to advancing soil moisture realism in next-generation Earth system models.
Previous studies have identified the oversimplified root system representation as a key factor leading to inaccuracies in vegetation‐atmosphere feedbacks. In this study, a dynamic root water uptake scheme in the Noah‐MP land surface model has been coupled to the Weather Research and Forecasting (WRF) model to investigate its impact on the surface hydroclimate variables and land‐atmosphere interactions. To evaluate the impact of the dynamic root, two coupled simulations were conducted, one with the dynamic root water uptake scheme (DynRt) and one with the static root water uptake scheme (StcRt), which is based on the default root representation in Noah‐MP, with slight modifications, primarily in vegetation‐related parameters. Both DynRt and StcRt simulations were conducted with a small ensemble of three members to account for variations in physical parameterizations, initial and boundary forcing and model setup. When compared with reference data sets, the DynRt simulations show improved results than the StcRt simulations, reducing biases in the simulated leaf area index, surface energy fluxes, soil moisture and precipitation. Two different mechanisms through which roots affect land‐atmosphere coupling have been identified. Over the transitional climate zone between the dry and wet climate, the dynamic root scheme affects surface climate and land‐atmosphere coupling mainly through changes in soil moisture through hydraulic redistribution by plant root system. Over the energy‐limited mesic zone, the dynamic root affects regional land‐atmosphere coupling mainly through changes in carbon allocation. This work highlights the importance of dynamic root representation in improving vegetation‐atmosphere simulations by enhancing predictions of water, energy, and carbon fluxes.
Land-atmospheric (L-A) feedbacks have historically been studied using models whose structure and parameterizations influence outcomes and insights. The representation of L-A feedbacks based on observations alone remains an ongoing challenge for understanding boundary layer development and precipitation. To address this gap, we use ground-based passive remote sensing and in-situ observations to present an analysis of the atmosphere during 103 interstorm soil moisture drydown events spanning nine warm seasons (2016–2024) in the U.S. Southern Great Plains region. By separating events based on local L-A coupling signals and characterizing the profiles of atmospheric heat and moisture to surface energy flux behavior, we investigate the physical mechanisms linking land surface processes to boundary layer development. We find that during interstorm drydowns, the atmospheric column follows a consistent pattern: moisture increases within the boundary layer, peaks near its top, and declines rapidly above, while warming occurs through the depth. Drydowns that shift toward evaporation produce stronger and deeper thermodynamic responses than cases dominated by sensible heating, which are weaker and shallower. Additionally, moisture is accumulated faster within the boundary layer during shorter drydowns, with longer drydowns representing slower, moisture-limited growth. Drydowns with wetter initial soil moisture will sustain stronger moistening within and above the boundary layer, accelerating buoyancy growth and convective potential toward the next storm. These results provide observational evidence linking surface flux evolution to boundary layer thermodynamics and offer a process-level benchmark for evaluating coupled L-A representations in models and demonstrating the influence of soil moisture on short-term weather forecasting skill.
Soil moisture is coupled with vegetation and atmosphere, influencing global cycling of water, carbon, and energy. However, it remains unclear how soil moisture-atmosphere interactions affect land-atmosphere carbon and water exchanges simultaneously. Using Earth system model experiments, we show widespread carbon-water trade-offs between net ecosystem production and precipitation-minus-evapotranspiration driven by soil moisture dynamics. Soil moisture positively controls net ecosystem production and negatively affects precipitation-minus-evapotranspiration, through direct soil water stress and indirect soil moisture-atmosphere feedbacks. While soil moisture variability magnifies the interannual variability of net ecosystem production, it moderates that of precipitation-minus-evapotranspiration over land. These opposing effects lead to a pronounced carbon-water trade-off, which originates from the interplay between carbon acquisition through photosynthesis and water extraction through evapotranspiration. This trade-off is projected to intensify in a warming and drying future, as soil moisture increasingly regulates carbon and water exchanges, posing a serious challenge to sustaining both terrestrial carbon sink and water supply.
Eddy covariance (EC), the gold standard for measuring ecosystem scale gas and heat exchanges, has transformed our understanding of the breathing of the biosphere, and thus global change biology. Despite numerous methodological improvements and insights gained from the technique, the community faces persistent challenges that have been present since the first EC measurements. Here, we review the theoretical developments underpinning EC. We present theoretical developments in four important areas that have relevance to EC measurements of the net ecosystem exchanges (NEE) of gases and heat from a single tower: (i) measuring the total vertical flux density, (ii) flux attenuation, (iii) coordinate rotations, and (iv) energy balance closure. Persistent problems with EC measurements, such as the inability to close the energy budget, led us to identify two priorities for revisiting the theory underlying: (i) sensible heat flux calculations, and (ii) constraining the mean vertical wind velocity. We present a framework for improved calculation of sensible heat flux derived from first principles of fluid mechanics and thermodynamics that considers coupled heat and mass transfer so that conservation of both is obeyed. These refinements are motivated by the need for unbiased measurements of energy and mass transfer between the land surface and atmosphere for ecosystem research and to validate satellite observations and land surface models.