Search NASA⌕ Search

SEARCH · Search NASA

Results for “circulation model”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Bimodality in simulated precipitation frequency distributions and its relationship with convective parameterizations

Abstract Bimodality in precipitation frequency distributions is often evident in atmospheric models, but rarely in observations. This study i) proposes a metric to objectively quantify the bimodality in precipitation distributions, ii) evaluates model simulations contributed to the Coupled Model Intercomparison Project (CMIP) phase 5 (CMIP5), phase 6 (CMIP6), and the DYnamics of the Atmospheric general circulation Modeled On Non-hydrostatic Domains (DYAMOND) project by comparing them to satellite-based and reanalysis precipitation products, and iii) investigates possible origins of bimodal precipitation distributions. Our results reveal that about 83% (20 out of 24) of CMIP5 and 70% (21 out of 30) of CMIP6 models used in this study exhibit bimodal distributions. The few DYAMOND models that use a deep convective parameterization also show bimodal distributions, while most DYAMOND models do not. Predictably, the bimodality originates from the separation of precipitation process between resolved grid-scale and parameterized subgrid-scale. However, in a larger number of models bimodality arises from the parameterized subgrid-scale convective precipitation alone.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Future projections in the climatology of global low-level jets from CORDEX-CORE simulations

The potential changes in the strength and location of five low-level jets (LLJs) located within four Coordinated Regional Climate Downscaling Experiment (CORDEX) domains are examined for present and future climate conditions using an ensemble of simulations conducted with the RegCM4 regional model at a 25 km horizontal grid spacing. Lateral and lower boundary forcing fields are from three General Circulation Models (GCMs), and we analyse a historical period (1995–2014) along with two future periods (2041–2060 and 2080–2099) under the Representative Concentration Pathways 2.6 and 8.5. The RegCM4, as driven by the GCMs, is capable of capturing most of the observed climatological features of the LLJs, both in terms of spatial location and seasonal evolution. Analysis of the influence of global warming on the LLJs shows a consistent strengthening of the jets and a shift in their location under both warming scenarios. The Monsoon and West African westerly LLJs exhibit a northward shift, while the Caribbean and South American LLJs present a westward expansion. The use of an ensemble of high-resolution simulations is found to provide a key element for a robust assessment of changes in LLJs associated with future global warming scenarios.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Non-Synchronous Rotation on Europa Driven by Ocean Currents

It has been suggested that the ice shell of Jupiter's moon Europa may drift non-synchronously due to tidal torques. Here we argue that torques applied by the underlying ocean are also important and can result in non-synchronous rotation. The resulting spin rate can be slightly slower than the synchronous angular rate that would have kept the same point of the ice shell facing Jupiter. We develop an ice shell rotation model, driven by ocean stress calculated using a high-resolution state-of-the-art ocean general circulation model, and take into account the viscoelastic deformation of the ice shell. We use the ice shell model results together with observed limits on the ice shell drift speed to constrain ice shell parameters such as effective viscosity, which is currently uncertain by at least four orders of magnitude. Our results suggest, at best, sluggish ice shell convection. Depending on the relaxation time scale of the ice shell and on the ocean currents, the ice shell may exhibit negligible drift, constant drift, or oscillatory drift superimposed on random fluctuations. The expected rotation rate exceeds ~30 m/yr; future spacecraft observations can be used to test these predictions and yield insight into the properties of the ice shell and underlying ocean.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

Dissipation Scaled Internal Wave Drag in a Global Heterogeneously Coupled Internal/External Mode Total Water Level Model

This study showcases a global, heterogeneously coupled total water level system wherein salinity and temperature outputs from a coarser-resolution (~12 km) ocean general circulation model are used to calculate density-driven terms within a global, higher-resolution (~2.5 km) depth-averaged total water level model. We demonstrate that the inclusion of baroclinic forcing in the barotropic model requires modification of the internal wave drag term to prevent excess degradation of tidal results compared to the barotropic model. By scaling the internal tide dissipation by an easy to calculate dissipation ratio, the resulting heterogeneously coupled model has complex root mean square errors (RMSE) of 2.27 cm in the deep ocean and 12.16 cm in shallow waters for the M 2 tidal constituent. While this represents a 10%–20% deterioration as compared to the barotropic model, the improvements in total water level prediction more than offset this degradation. Global median RMSE compared to observations of total water levels, 30-day sea levels, and non-tidal residuals improve by 1.86 (18.5%), 2.55 (42.5%), and 0.36 (5.3%) cm respectively. The drastic improvement in model performance highlights the importance of including density-driven effects within global hydrodynamic models and will help to improve the results of both hindcasts and forecasts in modeling extreme and nuisance flooding. With only an 11% increase in model run time compared to the fully barotropic total water level model, this approach paves the way for high resolution coastal water level and flood models to be used alongside climate models, improving operational forecasting of total water levels.

Blakely, Coleman Peter [University of Notre Dame, ↗
DOE OSTI2021

A clustering-based approach to ocean model–data comparison around Antarctica

The Antarctic Continental Shelf seas (ACSS) are a critical, rapidly changing element of the Earth system. Analyses of global-scale general circulation model (GCM) simulations, including those available through the Coupled Model Intercomparison Project, Phase 6 (CMIP6), can help reveal the origins of observed changes and predict the future evolution of the ACSS. However, an evaluation of ACSS hydrography in GCMs is vital: previous CMIP ensembles exhibit substantial mean-state biases (reflecting, for example, misplaced water masses) with a wide inter-model spread. Because the ACSS are also a sparely sampled region, grid-point-based model assessments are of limited value. Our goal is to demonstrate the utility of clustering tools for identifying hydrographic regimes that are common to different source fields (model or data), while allowing for biases in other metrics (e.g., water mass core properties) and shifts in region boundaries. We apply K-means clustering to hydrographic metrics based on the stratification from one GCM (Community Earth System Model version 2; CESM2) and one observation-based product (World Ocean Atlas 2018; WOA), focusing on the Amundsen, Bellingshausen and Ross seas. When applied to WOA temperature and salinity profiles, clustering identifies “primary” and “mixed” regimes that have physically interpretable bases. For example, meltwater-freshened coastal currents in the Amundsen Sea and a region of high-salinity shelf water formation in the southwestern Ross Sea emerge naturally from the algorithm. Both regions also exhibit clearly differentiated inner- and outer-shelf regimes. The same analysis applied to CESM2 demonstrates that, although mean-state model biases in water mass T–S characteristics can be substantial, using a clustering approach highlights that the relative differences between regimes and the locations where each regime dominates are well represented in the model. CESM2 is generally fresher and warmer than WOA and has a limited fresh-water-enriched coastal regimes. Given the sparsity of observations of the ACSS, this technique is a promising tool for the evaluation of a larger model ensemble (e.g., CMIP6) on a circum-Antarctic basis.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Future projections in tropical cyclone activity over multiple CORDEX domains from RegCM4 CORDEX-CORE simulations

The characteristics of tropical cyclone (TC) activity over five TC basins lying within four Coordinated Regional Downscaling Experiment (CORDEX) domains are examined for present and future climate conditions using a new ensemble of simulations completed as part of the CORDEX-CORE initiative with the regional climate model RegCM4. The simulations are conducted at a 25 km horizontal grid spacing and are driven by three CMIP5 general circulation models (GCMs) under two Representative Concentration Pathways (RCP2.6 and RCP8.5). The RegCM4 captures most features of the observed TC climatology, except for the TC intensity, which is thus statistically adjusted using a bias correction procedure to account for the effect of the coarse model resolution. The RegCM4 exhibits an improved simulation of several TC statistics compared to the driving GCMs, over most basins analyzed. In future climate conditions we find significant increases in TC frequency over the North Indian Ocean, the Northwest Pacific and Eastern Pacific regions, which are consistent with an increase in mid-tropospheric relative humidity. The North Atlantic and Australasia regions show a decrease in TC frequency, mostly associated with an increase in wind shear. We also find a consistent increase in future storm rainfall rates associated with TCs and in the frequency of the most intense TCs over most domains. Finally, our study shows robust responses often, but not always, in line with previous studies, still implying the presence of significant uncertainties in the projection of TC characteristics, which need to be addressed using large ensembles of simulations with high-resolution models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Simulations With the Marine Biogeochemistry Library (MARBL)

Abstract The Marine Biogeochemistry Library (MARBL) is a prognostic ocean biogeochemistry model that simulates marine ecosystem dynamics and the coupled cycles of carbon, nitrogen, phosphorus, iron, silicon, and oxygen. MARBL is a component of the Community Earth System Model (CESM); it supports flexible ecosystem configuration of multiple phytoplankton and zooplankton functional types; it is also portable, designed to interface with multiple ocean circulation models. Here, we present scientific documentation of MARBL, describe its configuration in CESM2 experiments included in the Coupled Model Intercomparison Project version 6 (CMIP6), and evaluate its performance against a number of observational data sets. The model simulates present‐day air‐sea CO 2 flux and many aspects of the carbon cycle in good agreement with observations. However, the simulated integrated uptake of anthropogenic CO 2 is weak, which we link to poor thermocline ventilation, a feature evident in simulated chlorofluorocarbon distributions. This also contributes to larger‐than‐observed oxygen minimum zones. Moreover, radiocarbon distributions show that the simulated circulation in the deep North Pacific is extremely sluggish, yielding extensive oxygen depletion and nutrient trapping at depth. Surface macronutrient biases are generally positive at low latitudes and negative at high latitudes. CESM2 simulates globally integrated net primary production (NPP) of 48 Pg C yr −1 and particulate export flux at 100 m of 7.1 Pg C yr −1 . The impacts of climate change include an increase in globally integrated NPP, but substantial declines in the North Atlantic. Particulate export is projected to decline globally, attributable to decreasing export efficiency associated with changes in phytoplankton community composition.

Meteorology & Atmospheric Sciences↗
DOE OSTI2021

Investigating the Causes and Impacts of Convective Aggregation in a High Resolution Atmospheric GCM

Abstract A ∼50 km resolution atmospheric general circulation model (GCM) is used to investigate the impact of radiative interactions on spatial organization of convection, the model's mean state, and extreme precipitation events in the presence of realistic boundary conditions. Mechanism‐denial experiments are performed in which synoptic‐scale feedbacks between radiation and dynamics are suppressed by overwriting the model‐generated atmospheric radiative cooling rates with its monthly varying climatological values. When synoptic‐scale radiative interactions are disabled, the annual mean circulation and precipitation remain almost unchanged, however tropical convection becomes less aggregated, with an increase in cloud fraction and relative humidity in the free troposphere but a decrease in both variables in the boundary layer. Changes in cloud fraction and relative humidity in the boundary layer exhibit more sensitivity to the presence of radiative interactions than variations in the degree of aggregation. The less aggregated state is associated with a decrease in the frequency of extreme precipitation events, coincident with a decrease in the dynamical contribution to the magnitude of extreme precipitation. At regional scales, the spatial contrast in radiative cooling between dry and moist regions diminishes when radiative interactions are suppressed, reducing the upgradient transport of energy, degree of aggregation, and frequency of extreme precipitation events. However, the mean width of the tropical rain belt remains almost unaffected when radiative interactions are disabled. These results offer insights into how radiation‐circulation coupling affects the spatial organization of convection, distributions of clouds and humidity, and weather extremes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Explicit Representation of Orographic Anisotropy for All Directions Improves Nanling Mountain Rainfall Simulation

Climate models exhibit significant rainfall bias in mountainous regions. One reason is the insufficient representation of orographic anisotropy in these models. In this study, we implement the orographic drag scheme with 3-D orographic anisotropy (all flow directions (AFD)) into a general circulation model and investigate its impact on Nanling rainfall simulation where models have systematic dry bias in the summer. It is shown that the AFD alleviated the Nanling mountain rainfall bias by over 60%. This is through an anomalous “lower-convergence-higher-divergence” deceleration pattern of the flow windward of the Nanling Mountains that enhanced vertical motion and low-level moisture convergence. The results suggest the importance of explicit orographic anisotropy representation in rainfall simulation in mountainous regions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Object-Based Evaluation of Tropical Precipitation Systems in DYAMOND Simulations over the Maritime Continent

The present study uses an object-based evaluation metric to examine the precipitation bias over the Maritime Continent in the global cloud-resolving models. We specifically focus on the difference between the models that directly resolve convection and those using convection parameterization. The 40-day hindcast experiments of the DYnamics of the Atmospheric general circulation Modeled On Non-hydrostatic Domain (DYAMOND) intercomparison project are evaluated against the high-resolution satellite rainfall products. The hindcast of the Central Weather Bureau Global Forecast System (CWBGFS) under the DYAMOND protocol is also included. The results indicate that most models simulate insufficient numbers of large precipitation system [object-based precipitation system (OPS), > 370 km in scale], indicating weaker convection organization. The observation indicates that the maximum precipitation within the OPS intensifies with increasing object size. All of the models capture this positive relationship, but most of them overestimate the sensitivity. Most of the models overestimate both the frequency and intensity of small OPS (< 160 km), except for the models with convection parameterization [i.e., CWBGFS, European Centre for Medium-Range Weather Forecasts Integrated Forecasting System (IFS)-9 km]. Although most of the models can reproduce the observed peak time of diurnal precipitation over the land area in the Maritime Continent, the simulated fractional contribution of different sizes of OPS to the total precipitation varies from model to model, and their peak times do not follow the observed ones with delayed peak times as the size of OPS increases from small, mid-size, to large categories. Most of the models reasonably capture the mean diurnal cycle peak time, but only the models with convection parameterization and Model for Prediction Across Scales (MPAS) can represent the diurnal evolution of fractional contribution from different OPSs. The implications of the current results to the upscale processes of the tropical convection systems in the global models are also discussed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

System code evaluation of near-term accident tolerant claddings during pressurized water reactor station blackout accidents

Following the Fukushima Daiichi nuclear accident in 2011, researches on Accident-tolerant fuels (ATFs) are currently of high interest in not only the nuclear industry but also governmental and international organizations. In this work, a quantitative evaluation of the performance of monolithic FeCrAl cladding and Cr-coated Zircaloy cladding has been performed for Pressurized Water Reactor (PWR) Station Blackout (SBO) accidents. A generic PWR model has been built in system thermal-hydraulics code TRACE based on the Surry Nuclear Power Station with counter-current natural circulation modelling capability for hotleg and steam generator U-tube components during the accidents. The base model results are then compared to MELCOR and RELAP simulations to verify the system component implementation in TRACE. Two PWR SBO scenarios were investigated, including: short-term SBO and long-term SBO with early reactor coolant pump (RCP) seal failure. These scenarios are defined to be very similar to the accidents studied in the State-of-the-Art Reactor Consequence Analysis (SOARCA) project. TRACE code is modified to reflect the oxidation kinetics of FeCrAl and Cr-coating. Larson-Miller creep rupture model is also implemented in TARCE using its built-in control systems to simulate the creep rupture of hotlegs. Additionally, the comparison between the TRACE models with and without the counter-current flow modeling resulted in significant difference when comparing ATF cladding to Zircaloy for short term SBO, while it marginal impacted the performance during long term SBO with RCP seal failure. For short term SBO, both ATF cladding underwent hot leg creep rupture ~20 min after Zircaloy cladding. While Zircaloy and Cr-coated cladding had generated significant amount of hydrogen gas (>10 kg) before the creep rupture event, FeCrAl cladding had only generated <0.5 kg of hydrogen gas. For long term SBO with RCP seal failure, significant hydrogen generation and fuel melting was predicted before hot leg creep rupture for the ATF cladding while providing only 10–20 min additional coping time compared to Zircaloy.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Assessing the Representation of Synoptic Variability Associated With California Extreme Precipitation in CMIP6 Models

Abstract Days of extreme precipitation over California are evaluated in Coupled Model Intercomparison Project Phase 6 (CMIP6) models and the ERA‐Interim reanalysis. In the current climate, the model spread in composited precipitation on extreme precipitation days is closely related to the magnitude of composited integrated vapor transport (IVT) across models, a proxy for the intensity of atmospheric rivers. Most models underestimate the magnitude of IVT associated with extreme precipitation, according to ERA‐Interim. This is due mostly to the contribution of moisture, which almost all models overestimate, while the contribution of lower‐tropospheric wind speed is generally closer to the reanalyses. Moreover, most models underestimate the variance in the latitude of maxima of numerous variables among days of extreme California precipitation. That is, in the general circulation models there is a lack of diversity in the latitude of the disturbances bringing winter precipitation to California. In the future climate, most models project a decrease in the frequency of southward‐displaced disturbances among California extreme precipitation days. Hence, the greatest increases in extreme precipitation are over northern California. However, the historical underestimate of the latitudinal variance of disturbances calls into question the reliability of these projections. This bias should be especially considered for dynamical downscaling efforts over the region.

Norris, Jesse↗
DOE OSTI2023

Parameterization and Explicit Modeling of Cloud Microphysics: Approaches, Challenges, and Future Directions

Cloud microphysical processes occur at the smallest end of scales among cloud-related processes and thus must be parameterized not only in large-scale global circulation models (GCMs) but also in various higher-resolution limited-area models such as cloud-resolving models (CRMs) and large-eddy simulation (LES) models. Instead of giving a comprehensive review of existing microphysical parameterizations that have been developed over the years, this study concentrates purposely on several topics that we believe are understudied but hold great potential for further advancing bulk microphysics parameterizations: multi-moment bulk microphysics parameterizations and the role of the spectral shape of hydrometeor size distributions; discrete vs “continuous” representation of hydrometeor types; turbulence-microphysics interactions including turbulent entrainment-mixing processes and stochastic condensation; theoretical foundations for the mathematical expressions used to describe hydrometeor size distributions and hydrometeor morphology; and approaches for developing bulk microphysics parameterizations. Also presented are the spectral bin scheme and particle-based scheme (especially, super-droplet method) for representing explicit microphysics. Their advantages and disadvantages are elucidated for constructing cloud models with detailed microphysics that are essential to developing processes understanding and bulk microphysics parameterizations. Particle-resolved direct numerical simulation (DNS) models are described as an emerging technique to investigate turbulence-microphysics interactions at the most fundamental level by tracking individual particles and resolving the smallest turbulent eddies in turbulent clouds. Outstanding challenges and future research directions are explored as well.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

An Integrated Hydroclimatic Assessment of Future Reservoir and Hydropower Operations in the U.S.

The engineering of rivers by dams is a formative feature of human-nature systems and the interconnectivity of water, energy, and the climate. Sufficient and broad-based representations of dams in large-scale hydrological models prove essential to mapping their extensive regulation of river flow and biogeochemistry and gauging climate-linked provisions, including freshwater supply and hydropower. We present an integrated modeling framework to investigate future streamflow and hydropower generation in the Contiguous U.S. (1990–2075), leveraging an ensemble of six downscaled and bias-corrected General Circulation Models (GCMs) from the high-end SSP585 scenario of the CMIP6. To achieve this, we develop a reservoir operations and parameterization scheme for 1,384 dams in a high-resolution river network, including simulated hydropower generation for 326 dams. For the GCM ensemble mean, we simulate a widespread increase in regulated streamflow into the late-century (11% annual and 17% in winter for the dam median) with region-specific changes in summer streamflow that feature prominent declines in the Northwest (−7%). Mediation by reservoirs is shown to dampen intra-annual streamflow changes, delivering additional summer releases that partially mitigate declining flows. Total hydropower generation is projected to increase modestly (+3%), with boosted generation in the winter (+9%) and spring (+5%) offsetting declined summer generation (−3.4%), suggesting strong adaptation potential for hydropower in the future energy portfolio. Further analysis reveals that the choice of GCM, particularly in western regions, has significant bearing on projected streamflow and hydropower changes.

13 HYDRO ENERGY↗
DOE OSTI2024

Global Corn Heat Stress: Mean and SD of Degree Days Above 29°C based on NEX-GDDP-CMIP6 Climate Projections

Description This global dataset provides the estimated mean and standard deviation (SD) of corn heat stress (degree days above 29°C) for a set of climate models in NEX-GDDP-CMIP6 at 0.25-degree resolution. The NEX-GDDP-CMIP6 dataset is comprised of global downscaled climate scenarios derived from the General Circulation Model (GCM) runs conducted under the Coupled Model Intercomparison Project Phase 6 (CMIP6). The current dataset includes: Long-Term Average Degree Days Above 29°C- Historical Long-Term Average Degree Days Above 29°C- SSP245 Long-Term Standard Deviation of Degree Days Above 29°C- Historical Long-Term Standard Deviation of Degree Days Above 29°C- SSP245 The mean and SD are calculated over 1985-2014 for the historical period and over 2035-2064 for future projections. A full description of methods, including growing season, daily temperature distribution, and statistical coefficients, can be found in Haqiqi (2024). The source climate data are obtained from https://ds.nccs.nasa.gov/thredds2/catalog/catalog.html and are described in Thrasher et al (2022). The codes used to create this dataset are available at https://github.com/ihaqiqi/dd29c_nex_cmip6. Acknowledgments This work was supported by the US Department of Energy, Office of Science, Biological and Environmental Research Program, Earth and Environmental Systems Modeling, MultiSector Dynamics under Cooperative Agreement DE-SC0022141. The data processing, computation, and storage were completed on Purdue Anvil supercomputer and cyberinfrastructure supported by the National Science Foundation HDR award # 2118329: "NSF Institute for Geospatial Understanding through an Integrative Discovery Environment (I-GUIDE)". References Haqiqi. I. (2024). Trade can buffer climate-induced risks and volatilities in crop supply. Environmental Research: Food Systems. https://doi.org/10.1088/2976-601X/ad7d12 Thrasher, B., Wang, W., Michaelis, A., Melton, F., Lee, T., & Nemani, R. (2022). NASA global daily downscaled projections, CMIP6. Scientific Data, 9(1), 262. https://doi.org/10.1038/s41597-022-01393-4

Climate Change↗
DOE OSTI2024

Why Does Atmospheric Radiative Heating Weaken Midlatitude Cyclones?

Abstract Recent work has indicated that atmospheric radiative heating reduces the kinetic energy of large‐scale eddies in the midlatitudes. However, a physical mechanism that connects radiation to the midlatitude eddy kinetic energy is still uncertain. Using a high‐resolution general circulation model we perform an experiment in which the radiative cooling profile at each model time step is overwritten with the climatological mean, computed from a control simulation. This approach separates the mean and transient effects of radiative heating on the extratropical circulation. We find that, when radiative heating is fixed, the globally‐averaged eddy kinetic energy is enhanced by ∼6%. We show that thermal radiation dampens temperature anomalies near the surface and tropopause in low‐pressure systems, destroying eddy available potential energy and eddy kinetic energy. We identify this as a possible mechanism by which atmospheric radiative heating weakens midlatitude cyclones.

Mischell, Eric↗
DOE OSTI2021

Tropical Cirrus in Global Storm-Resolving Models. Part II: Cirrus Life Cycle and Top-of-Atmosphere Radiative Fluxes

Cirrus clouds of various thicknesses and radiative characteristics extend over much of the tropics, especially around deep convection. They are difficult to observe due to their high altitude and sometimes small optical depths. They are also difficult to simulate in conventional global climate models, which have coarse grid spacings and simplified parameterizations of deep convection and cirrus formation. We investigate the representation of tropical cirrus in global storm-resolving models (GSRMs), which have higher spatial resolution and explicit convection and could more accurately represent cirrus cloud processes. This study uses GSRMs from the DYnamics of the Atmospheric general circulation Modeled On Non-hydrostatic Domains (DYAMOND) project. The aggregate life cycle of tropical cirrus is analyzed using joint albedo and outgoing longwave radiation (OLR) histograms to assess the fidelity of models in capturing the observed cirrus cloud populations over representative tropical ocean and land regions. The proportions of optically-thick deep convection, anvils, and cirrus vary across models and are portrayed in the vertical distribution of cloud cover and top-of-atmosphere radiative fluxes. Model differences in cirrus populations, likely driven by subgrid processes such as ice microphysics, dominate over regional differences between convectively-active tropical land and ocean locations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Hierarchical Testing of a Hybrid Machine Learning‐Physics Global Atmosphere Model

Machine learning (ML)-based models have demonstrated high skill and computational efficiency, often outperforming conventional physics-based models in weather and subseasonal predictions. While prior studies have assessed their fidelity in capturing synoptic-scale atmospheric dynamics, their performance across timescales and under out-of-distribution forcing, such as +3K or +4K uniform-warming forcings, and the sources of biases remain elusive, to establish the model's reliability for Earth science. Here, we design three sets of experiments targeting synoptic-scale phenomena, interannual variability, and out-of-distribution uniform-warming forcings. We evaluate the Neural General Circulation Model (NeuralGCM), a hybrid model integrating a dynamical core with ML-based component, against observations and physics-based Earth system models (ESMs). At the synoptic scale, NeuralGCM captures the evolution and propagation of extratropical cyclones with performance comparable to ESMs. At the interannual scale, when forced by El Niño-Southern Oscillation sea surface temperature (SST) anomalies, NeuralGCM successfully reproduces associated teleconnection patterns but exhibits deficiencies in capturing nonlinear response. Under out-of-distribution uniform-warming forcings, NeuralGCM simulates similar responses in global-average temperature and precipitation and reproduces large-scale tropospheric circulation features similar to those in ESMs. Notable weaknesses include overestimating the tracks and spatial extent of extratropical cyclones, biases in the teleconnected wave train triggered by tropical SST anomalies, and differences in upper-level warming and stratospheric circulation responses to SST warming compared to physics-based ESMs. The causes of these weaknesses were explored. Despite the noted weaknesses, NeuralGCM reproduces responses across experiments reasonably and performs comparably to ESMs. By integrating a dynamical core with ML, NeuralGCM shows potential for developing ML-based ESMs.

global warming↗