Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model intercomparison projects”

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.

At least 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

The fully coupled regionally refined model of E3SM version 2: overview of the atmosphere, land, and river results

Abstract. This paper provides an overview of the United States (US) Department of Energy's (DOE's) Energy Exascale Earth System Model version 2 (E3SMv2) fully coupled regionally refined model (RRM) and documents the overall atmosphere, land, and river results from the Coupled Model Intercomparison Project 6 (CMIP6) DECK (Diagnosis, Evaluation, and Characterization of Klima) and historical simulations – a first-of-its-kind set of climate production simulations using RRM. The North American (NA) RRM (NARRM) is developed as the high-resolution configuration of E3SMv2 with the primary goal of more explicitly addressing DOE's mission needs regarding impacts to the US energy sector facing Earth system changes. The NARRM features finer horizontal resolution grids centered over NA, consisting of 25→100 km atmosphere and land, a 0.125∘ river-routing model, and 14→60 km ocean and sea ice. By design, the computational cost of NARRM is ∼3× of the uniform low-resolution (LR) model at 100 km but only ∼ 10 %–20 % of a globally uniform high-resolution model at 25 km. A novel hybrid time step strategy for the atmosphere is key for NARRM to achieve improved climate simulation fidelity within the high-resolution patch without sacrificing the overall global performance. The global climate, including climatology, time series, sensitivity, and feedback, is confirmed to be largely identical between NARRM and LR as quantified with typical climate metrics. Over the refined NA area, NARRM is generally superior to LR, including for precipitation and clouds over the contiguous US (CONUS), summertime marine stratocumulus clouds off the coast of California, liquid and ice phase clouds near the North Pole region, extratropical cyclones, and spatial variability in land hydrological processes. The improvements over land are related to the better-resolved topography in NARRM, whereas those over ocean are attributable to the improved air–sea interactions with finer grids for both atmosphere and ocean and sea ice. Some features appear insensitive to the resolution change analyzed here, for instance the diurnal propagation of organized mesoscale convective systems over CONUS and the warm-season land–atmosphere coupling at the southern Great Plains. In summary, our study presents a realistically efficient approach to leverage the fully coupled RRM framework for a standard Earth system model release and high-resolution climate production simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Significant impact of forcing uncertainty in a large ensemble of climate model simulations

Forcing due to solar and volcanic variability, on the natural side, and greenhouse gas and aerosol emissions, on the anthropogenic side, are the main inputs to climate models. Reliable climate model simulations of past and future climate change depend crucially upon them. Here we analyze large ensembles of simulations using a comprehensive Earth System Model to quantify uncertainties in global climate change attributable to differences in prescribed forcings. The different forcings considered here are those used in the two most recent phases of the Coupled Model Intercomparison Project (CMIP), namely CMIP5 and CMIP6. We show significant differences in simulated global surface air temperature due to volcanic aerosol forcing in the second half of the 19th century and in the early 21st century. The latter arise from small-to-moderate eruptions incorporated in CMIP6 simulations but not in CMIP5 simulations. We also find significant differences in global surface air temperature and Arctic sea ice area due to anthropogenic aerosol forcing in the second half of the 20th century and early 21st century. These differences are as large as those obtained in different versions of an Earth System Model employing identical forcings. In simulations from 2015 to 2100, we find significant differences in the rates of projected global warming arising from CMIP5 and CMIP6 concentration pathways that differ slightly but are equivalent in terms of their nominal radiative forcing levels in 2100. Our results highlight the influence of assumptions about natural and anthropogenic aerosol loadings on carbon budgets, the likelihood of meeting Paris targets, and the equivalence of future forcing scenarios.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Blocking Simulations in GFDL GCMs for CMIP5 and CMIP6

Abstract The frequency of atmospheric blocking has been largely underestimated by general circulation models (GCMs) participating in the Coupled Model Intercomparison Project (CMIP). Errors in the onset, persistence, barotropicity, geographical preference, seasonality, intensity, and moving speed of global blocking were diagnosed in 10 Geophysical Fluid Dynamics Laboratory (GFDL) GCMs for recent CMIP5 and CMIP6 using a detection approach that combines zonal eddies and the reversal of zonal winds. The blocking frequency, similar at 500 and 250 hPa, is underestimated by 50% in the Atlantic–Europe region during December–February but is overestimated by 60% in the Pacific–North America region during that season and by 70% in the southwest Pacific during July–August. These blocking biases at 500 hPa were investigated in the five CMIP6 models that showed improvements over the CMIP5 versions. The Atlantic–Europe underestimate corresponds to lower instantaneous blocking rates, lower persistent blocking rates, and higher persistent stationary ridge rates; the number of blocks with a duration of 4–5 days is only 40%–65% of that in observations. In contrast, the overestimate consists of excessive blocks with a duration longer than 12 days in the Pacific–North America and up to twice as many 4–6-day events in the southwest Pacific. Simulated December–February blocks up to 12 days in the Pacific–North America region tend to be stronger and to move more slowly than those in observations. Diagnostic sensitivity tests indicated that the zonal mean and zonal eddy components of the mean state play a key role, as replacing each with that of observations substantially reduced many of the outstanding biases in these GCMs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

The Southern hemisphere midlatitude circulation response to rapid adjustments and sea surface temperature driven feedbacks

Rapid adjustments – the response of meteorology to external forcing while sea surface temperatures (SST) and sea ice are held fixed – can affect the midlatitude circulation and contribute to long-term forced circulation responses in climate simulations. This study examines rapid adjustments in the Southern hemisphere (SH) circulation using nine models from the Precipitation Driver and Response Model Intercomparison Project (PDRMIP), which perform fixed SST and coupled ocean experiments for five perturbations: a doubling of carbon dioxide (2xCO2), a tripling of methane (3xCH4), a five-fold increase in sulphate aerosol (5xSO4), a ten-fold increase in black carbon aerosol (10xBC) and a 2% increase in solar constant (2%Sol). In the coupled experiments, the SH eddy-driven jet shifts poleward and strengthens for forcings that produce global warming (and vice versa for 5xSO4), with the strongest response found in austral summer. In austral winter, the responses project more strongly onto a change in jet strength. For 10xBC, which induces strong shortwave absorption, the multi-model mean (MMM) rapid adjustment in DJF jet latitude is ~75% of the change in the coupled simulations. For the other forcings, which induce larger SST changes, the effect of SST-mediated feedbacks on the SH circulation is larger than the rapid adjustment. Nevertheless, for these perturbations the magnitude of the MMM jet shift due to the rapid adjustment is still around 25-33% of that in the coupled experiments. The results demonstrate the need to understand the mechanisms for rapid adjustments in the midlatitude circulation, in addition to the effect of changing SSTs.

Wood, Thomas↗
DOE OSTI2022

On the Correspondence Between Atmosphere-Only and Coupled Simulations for Radiative Feedbacks and Forcing From CO 2

Atmosphere-only experiments are widely used to investigate climate feedbacks simulated in more computationally expensive fully coupled global climate model simulations. We confirm that this remains a valid approach by comparing the radiative feedbacks and forcing between coupled and atmosphere-only simulations for the latest models taking part in the 6th phase of the Coupled Model Intercomparison Project (CMIP6). For global-mean cloud feedbacks, we find a better than previously known correspondence between these experiments, which applies even to the response of individual cloud properties (amount, altitude, and optical depth) and holds even when considering atmosphere-only simulations of only 1 yr duration. For regional cloud feedbacks, the correspondence between the two experiments is generally present at every geographic location except for the tropical Pacific but takes longer experiments to reveal. For the lapse rate and surface albedo feedbacks, the correspondence between the two experiments is weaker due to the non-uniform warming pattern and loss of sea ice in the coupled experiment. For the across-model relationship between 4xCO 2 effective radiative forcing and feedback, we find a different behavior across experiments in CMIP6 than in CMIP5, casting doubt on the physical significance of previous results that highlighted an anti-correlation between the two quantities. Overall, these results confirm the utility of atmosphere-only experiments particularly to study cloud feedbacks, which are the dominant source of inter-model spread in climate sensitivity.

58 GEOSCIENCES↗
DOE OSTI2024

Understanding the Cascade: Removing GCM Biases Improves Dynamically Downscaled Climate Projections

Polarization surrounding bias correction (BC) in creating climate projections arises from its lack of physicality. Here, we perform and analyze 18 dynamical downscaling simulations (with and without BC) to better understand the physical impacts of BC, applied before downscaling, on regional climate output across the western United States. Without BC, downscaled precipitation is systematically and unrealistically wet biased compared to a hierarchy of observationally based datasets over the 1980–2014 period due to cascading mean–state Global Climate Model (GCM) biases: (a) overly strong lower–tropospheric lapse rates (5 K/km), (b) overly cold (2 K) tropospheric temperatures, and (c) anomalous mid–tropospheric cyclonic vorticity advection. With BC, downscaled precipitation (snow) biases are virtually eliminated (halved). Identified GCM biases are common to the broader Coupled Model Intercomparison Project ensemble. Physical effects of BC on the quality of the regionalized projections, pending an evaluation of BC's distortion of the downscaled climate response, may motivate its broader application by dynamical downscalers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Impact of solar geoengineering on temperatures over the Indonesian Maritime Continent

Climate change has been projected to increase the intensity and magnitude of extreme temperature in Indonesia. Solar radiation management (SRM) has been proposed as a strategy to temporarily combat global warming, buying time for negative emissions. Though the global impacts of SRM have been extensively studied in recent years, regional impacts, especially in the tropics, have received much less attention. This paper investigates the potential stratospheric sulfate aerosol geoengineering to modify mean and extreme temperature, as well as the relative humidity and wet bulb temperature (WBT) change over Indonesian Maritime Continent (IMC) based on simulations from three different earth system models. We applied statistical downscaling and corrected the bias of model output to reproduce historical temperatures and relative humidity over IMC. We compared SRM as specified by the Geoengineering Model Intercomparison Project (GeoMIP) G4 experiment – 5 Tg of SO 2 into the lower stratosphere between 2020 and 2069 – concurrent with the RCP4.5 emissions scenario. SRM is able to significantly reduce the temperature means and extremes, and although differences in magnitude of response and spatial pattern occur, there is a generally consistent response. Additionally, the spatial response of changes forced by RCP4.5 scenario and SRM impact is notably heterogeneous in the archipelago, highlighting uncertainties that would be critical in assessing socio-economic consequences of both doing, and not doing SRM. In general, SRM has bigger impacts in reducing temperatures over land than oceans, and the southern monsoon region shows more variability. SRM is also effective at reducing the likelihood of WBT > 27OC events compared to RCP4.5 after some years of SRM deployment as well as during the post-termination period of SRM. Regional downscaling may be an effective tool in obtaining policy-relevant information about local effects of different future scenarios involving SRM.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Model spread in tropical low cloud feedback tied to overturning circulation response to warming

Among models participating in the Coupled Model Intercomparison Project phase 6 (CMIP6), here we show that the magnitude of the tropical low cloud feedback, which contributes considerably to uncertainty in estimates of climate sensitivity, is intimately linked to tropical deep convection and its effects on the tropical atmospheric overturning circulation. First, a reduction in tropical ascent area and an increased frequency of heavy precipitation result in high cloud reduction and upper-tropospheric drying, which increases longwave cooling and reduces subsidence weakening, favoring low cloud reduction (Radiation-Subsidence Pathway). Second, increased longwave cooling decreases tropospheric stability, which also reduces subsidence weakening and low cloudiness (Stability-Subsidence Pathway). In summary, greater high cloud reduction and upper-tropospheric drying (negative longwave feedback) lead to a more positive cloud feedback among CMIP6 models by contributing to a greater reduction in low cloudiness (positive shortwave feedback). Varying strengths of the two pathways contribute considerably to the intermodel spread in climate sensitivity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Evaluation of Global Climate Models for Use in Energy Analysis

The interplay between energy, climate, and weather is becoming more complex due to increasing contributions of renewable energy generation, energy storage, electrified end uses, and the increasing frequency of extreme weather events. Energy system analyses commonly rely on meteorological inputs to estimate renewable energy generation and energy demand; however, these inputs rarely represent the estimated impacts of future climate change. Climate models and publicly available climate change datasets can be used for this purpose, but the selection of inputs from the myriad of available models and datasets is a nuanced and subjective process. In this work, we assess datasets from various global climate models (GCMs) from the Coupled Model Intercomparison Project Phase 6 (CMIP6). We present evaluations of their skills with respect to the historical climate and comparisons of their future projections of climate change for two climate change scenarios. We present the results for different climatic and energy system regions and include interactive figures in the accompanying software repository. Previous work has presented similar GCM evaluations, but none have presented variables and metrics specifically intended for comprehensive energy systems analysis including impacts on energy demand, thermal cooling, hydropower, water availability, solar energy generation, and wind energy generation. We focus on GCM output meteorological variables that directly affect these energy system components including the representation of extreme values that can drive grid resilience events. The objective of this work is not to recommend the best climate model and dataset for a given analysis, but instead to provide a reference to facilitate the selection of climate models and scenarios in subsequent work.

14 SOLAR ENERGY↗
DOE OSTI2022

Warming Pattern over the Northern Hemisphere Midlatitudes in Boreal Summer 1979–2020

Abstract Significant surface air temperature warming during summer 1979–2020 is not uniformly distributed in the northern midlatitudes over land but rather is confined to several longitudinal sectors including Europe, central Siberia and Mongolia, and both coasts of North America. These hot spots are accompanied by a chain of high pressure ridges from an anomalous, circumglobal Rossby wave train in the upper troposphere. From reanalysis data and several baseline experiments from phase 6 of the Coupled Model Intercomparison Project (CMIP6), we find that the circulation trend pattern is associated with fluctuations of the Atlantic multidecadal variability (AMV) and the interdecadal Pacific oscillation. The phase shift of AMV in the 1990s is particularly noteworthy for accelerating warming averaged over the northern midlatitude land. The amplitude of the observed trend in both surface air temperature and the upper-level geopotential height generally falls beyond the range of multidecadal trends simulated by the CMIP6 preindustrial control runs, supporting the likelihood that anthropogenic forcing played a critical role in the observed trend. On the other hand, the fidelity of the simulated low-frequency modes of variability and their teleconnections, especially on multidecadal time scales, is difficult to assess because of the relatively short observational records. Our mechanistic modeling results indicate that synoptic eddy–mean flow interaction is a key to the formation of the anomalous wave train but how the multidecadal modes can modulate the synoptic eddies through atmosphere–ocean and atmosphere–land interactions remains poorly understood. This gap in our knowledge makes it challenging to quantify the roles of the low-frequency modes and external forcings in causing the observed multidecadal trends.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Marine ice sheet experiments with the Community Ice Sheet Model

Ice sheet models differ in their numerical treatment of dynamical processes. Simulations of marine-based ice are sensitive to the choice of Stokes flow approximation and basal friction law and to the treatment of stresses and melt rates near the grounding line. We study the effects of these numerical choices on marine ice sheet dynamics in the Community Ice Sheet Model (CISM). In the framework of the Marine Ice Sheet Model Intercomparison Project 3d (MISMIP3d), we show that a depth-integrated, higher-order solver gives results similar to a 3D (Blatter–Pattyn) solver. We confirm that using a grounding line parameterization to approximate stresses in the grounding zone leads to accurate representation of ice sheet flow with a resolution of ~2 km, as opposed to ~0.5 km without the parameterization. In the MISMIP+ experimental framework, we compare different treatments of sub-shelf melting near the grounding line. In contrast to recent studies arguing that melting should not be applied in partly grounded cells, it is usually beneficial in CISM simulations to apply some melting in these cells. This suggests that the optimal treatment of melting near the grounding line can depend on ice sheet geometry, forcing, or model numerics. In both experimental frameworks, ice flow is sensitive to the choice of basal friction law. To study this sensitivity, we evaluate friction laws that vary the connectivity between the basal hydrological system and the ocean near the grounding line. CISM yields accurate results in steady-state and perturbation experiments at a resolution of ~2 km (arguably 4 km) when the connectivity is low or moderate and ~1 km (arguably 2 km) when the connectivity is strong.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

The GFDL Variable-Resolution Global Chemistry-Climate Model for Research at the Nexus of US Climate and Air Quality Extremes

We present a variable-resolution global chemistry-climate model (AM4VR) developed at NOAA's Geophysical Fluid Dynamics Laboratory (GFDL) for research at the nexus of US climate and air quality extremes. AM4VR has a horizontal resolution of 13 km over the US, allowing it to resolve urban-to-rural chemical regimes, mesoscale convective systems, and land-surface heterogeneity. With the resolution gradually reducing to 100 km over the Indian Ocean, we achieve multi-decadal simulations driven by observed sea surface temperatures at 50% of the computational cost for a 25-km uniform-resolution grid. In contrast with GFDL's AM4.1 contributing to the sixth Coupled Model Intercomparison Project at 100 km resolution, AM4VR features much improved US climate mean patterns and variability. In particular, AM4VR shows improved representation of: precipitation seasonal-to-diurnal cycles and extremes, notably reducing the central US dry-and-warm bias; western US snowpack and summer drought, with implications for wildfires; and the North American monsoon, affecting dust storms. AM4VR exhibits excellent representation of winter precipitation, summer drought, and air pollution meteorology in California with complex terrain, enabling skillful prediction of both extreme summer ozone pollution and winter haze events in the Central Valley. AM4VR also provides vast improvements in the process-level representations of biogenic volatile organic compound emissions, interactive dust emissions from land, and removal of air pollutants by terrestrial ecosystems. We highlight the value of increased model resolution in representing climate–air quality interactions through land-biosphere feedbacks. AM4VR offers a novel opportunity to study global dimensions to US air quality, especially the role of Earth system feedbacks in a changing climate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Estimating Uncertainty in Simulated ENSO Statistics

Abstract Large ensembles of model simulations are frequently used to reduce the impact of internal variability when evaluating climate models and assessing climate change induced trends. However, the optimal number of ensemble members required to distinguish model biases and climate change signals from internal variability varies across models and metrics. Here we analyze the mean, variance and skewness of precipitation and sea surface temperature in the eastern equatorial Pacific region often used to describe the El Niño–Southern Oscillation (ENSO), obtained from large ensembles of Coupled model intercomparison project phase 6 climate simulations. Leveraging established statistical theory, we develop and assess equations to estimate, a priori, the ensemble size or simulation length required to limit sampling‐based uncertainties in ENSO statistics to within a desired tolerance. Our results confirm that the uncertainty of these statistics decreases with the square root of the time series length and/or ensemble size. Moreover, we demonstrate that uncertainties of these statistics are generally comparable when computed using either pre‐industrial control or historical runs. This suggests that pre‐industrial runs can sometimes be used to estimate the expected uncertainty of statistics computed from an existing historical member or ensemble, and the number of simulation years (run duration and/or ensemble size) required to adequately characterize the statistic. This advance allows us to use existing simulations (e.g., control runs that are performed during model development) to design ensembles that can sufficiently limit diagnostic uncertainties arising from simulated internal variability. These results may well be applicable to variables and regions beyond ENSO.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Regime-oriented causal model evaluation of Atlantic–Pacific teleconnections in CMIP6

The climate system and its spatio-temporal changes are strongly affected by modes of long-term internal variability, like the Pacific decadal variability (PDV) and the Atlantic multidecadal variability (AMV). As they alternate between warm and cold phases, the interplay between PDV and AMV varies over decadal to multidecadal timescales. Here, we use a causal discovery method to derive fingerprints in the Atlantic–Pacific interactions and to investigate their phase-dependent changes. Dependent on the phases of PDV and AMV, different regimes with characteristic causal fingerprints are identified in reanalyses in a first step. In a second step, a regime-oriented causal model evaluation is performed to evaluate the ability of models participating in the Coupled Model Intercomparison Project Phase 6 (CMIP6) in representing the observed changing interactions between PDV, AMV and their extra-tropical teleconnections. The causal graphs obtained from reanalyses detect a direct opposite-sign response from AMV to PDV when analyzing the complete 1900–2014 period and during several defined regimes within that period, for example, when AMV is going through its negative (cold) phase. Reanalyses also demonstrate a same-sign response from PDV to AMV during the cold phase of PDV. Historical CMIP6 simulations exhibit varying skill in simulating the observed causal patterns. Generally, large-ensemble (LE) simulations showed better network similarity when PDV and AMV were out of phase compared to other regimes. Also, the two largest ensembles (in terms of number of members) were found to contain realizations with similar causal fingerprints to observations. For most regimes, these same models showed higher network similarity when compared to each other. This work shows how causal discovery on LEs complements the available diagnostics and statistical metrics of climate variability to provide a powerful tool for climate model evaluation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Evaluation of CMIP6 GCMs Over the CONUS for Downscaling Studies

Despite the necessity of Global Climate Models (GCMs) sub-selection in downscaling studies, an objective approach for their selection is currently lacking. Building on the previously established concepts in GCMs evaluation frameworks, we develop a weighted averaging technique to remove the redundancy in the evaluation criteria and rank 37 GCMs from the sixth phase of the Coupled Models Intercomparison Project over the contiguous United States. GCMs are rated based on their average performance across 66 evaluation measures in the historical period (1981–2014) after each metric is weighted between zero and one, depending on its uniqueness. The robustness of the outcome is tested by repeating the process with the empirical orthogonal function analysis in which each GCM is ranked based on its sum of distances from the reference in the principal component space. The two methodologies work in contrasting ways to remove the metrics redundancy but eventually develop similar GCMs rankings. A disparity in GCMs' behavior related to their sensitivity to the size of the evaluation suite is observed, highlighting the need for comprehensive multi-variable GCMs evaluation at varying timescales for determining their skillfulness over a region. The sub-selection goal is to use a representative set of skillful models over the region of interest without substantial overlap in their future climate responses and modeling errors in representing historical climate. Additional analyses of GCMs' independence and spread in their future projections provide the necessary information to objectively select GCMs while keeping all aspects of necessity in view.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

The Effect of an Equatorial Continent on the Tropical Rain Belt. Part 1: Annual Mean Changes in the ITCZ

The TRACMIP (Tropical Rain Belts with an Annual Cycle and Continent Model Intercomparison Project) ensemble includes slab-ocean aquaplanet controls and experiments with a highly idealized tropical continent, characterized by modified aquaplanet grid cells with increased evaporative resistance, increased albedo, reduced heat capacity, and no ocean heat transport (zero Q-flux). In the annual mean, an equatorial cold tongue develops west of the continent and induces dry anomalies and a split in the oceanic intertropical convergence zone (ITCZ). Ocean cooling is initiated by advection of cold, dry air from the winter portion of the continent; warm, humid anomalies in the summer portion are restricted to the continent by anomalous surface convergence. The surface energy budget suggests that ocean cooling persists and intensifies because of a positive feedback between a colder surface, drier and colder air, reduced downwelling longwave (LW) flux, and enhanced net surface LW cooling (LW feedback). A feedback between wind, evaporation, and SST (so-called WES feedback) also contributes to the establishment and maintenance of the cold tongue. Simulations with a gray-radiation model and simulations that diverge from protocol (with negligible winter cooling) confirm the importance of moist-radiative feedbacks and of rectification effects on the seasonal cycle. This mechanism coupling the continental and oceanic climate might be relevant to the double ITCZ bias. The key role of the LW feedback suggests that the study of interactions between monsoons and oceanic ITCZs requires full-physics models and a hierarchy of land models that considers evaporative processes alongside heat capacity as a defining characteristic of land.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Evaluating the Effects of Precipitation and Evapotranspiration on Soil Moisture Variability Within CMIP5 Using SMAP and ERA5 Data

=The effects of precipitation (Pr) and evapotranspiration (ET) on surface soil moisture (SSM) play an essential role in the land-atmosphere system. Here we evaluate multimodel differences of these effects within the Coupled Model Intercomparison Project Phase 5 (CMIP5) compared to Soil Moisture Active Passive (SMAP) products and ECMWF Reanalysis v5 (ERA5) as references in a frequency domain. The variability of SSM, Pr, and ET within three frequency bands (1/7 ~ 1/30 days -1 , 1/30 ~ 1/90 days -1 , and 1/90 ~ 1/365 days -1 ) after normalization is quantified using Fourier transform. We analyze the impact of ET and Pr on SSM variability based on a transfer function assuming that these variables form a linear time-invariant (LTI) system. For the total effects of ET and Pr on SSM variability, the CMIP5 estimations are smaller than the reference data in the two higher frequency bands and are larger than the reference data in the lowest frequency band. Besides, the effects on SSM by Pr and ET are found to be different across the three frequency bands. In each frequency band, the variability of the factor that dominates SSM (i.e., Pr or ET) from CMIP5 is smaller than that from the references. This study identifies the spatiotemporal distribution of differences between CMIP5 models and references (SMAP and ERA5) in simulating ET and Pr effects on SSM within three frequency bands. This study provides insightful information on how soil moisture variability is affected by varying precipitation and evapotranspiration at different time scales within Earth System Models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Overestimated natural biological nitrogen fixation translates to an exaggerated CO 2 fertilization effect in Earth system models

CO 2 fertilization of the terrestrial biosphere is limited by nitrogen. Biological nitrogen fixation (BNF) is the dominant natural nitrogen source to the terrestrial biosphere and can alleviate nitrogen limitation but is poorly constrained in Earth system models (ESMs). Here, in this study, we compare terrestrial BNF from an ensemble of ESMs of the 6th Coupled Model Intercomparison Project to a new global synthesis of observations across natural and agricultural biomes. We find that compared to observations, ESMs underestimate agricultural BNF but overestimate natural BNF in the present day by over 50%. Natural BNF is overestimated in the most productive ecosystems that contribute most to the terrestrial carbon sink (forests and grasslands). ESMs with different BNF representations yield a range of BNF responses to CO 2 enrichment. Some ESMs with phenomenological representations of BNF predict a natural BNF increase in response to a doubling of CO 2 that aligns with a meta-analysis of CO 2 enrichment experiments (31% increase) but fail to account for the substantial carbon cost of BNF. In contrast, ESMs with mechanistic representations of BNF account for its carbon cost as well as its regulation by nitrogen limitation but overestimate the BNF response to a doubling of CO 2 (135% increase). Overall, all current BNF representations in ESMs fall short of fully capturing its response to rising atmospheric CO 2 . Finally, we find a positive correlation between modeled present-day natural BNF and the CO 2 fertilization effect across ESMs, suggesting that overestimated natural BNF translates to an exaggerated CO 2 fertilization effect of approximately 11% in ESMs.

Biological nitrogen fixation↗