Search NASA⌕ Search

SEARCH · Search NASA

Results for “model reanalysis”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Stratosphere-Troposphere Coupling of Extreme Stratospheric Wave Activity in CMIP6 Models

Extreme stratospheric wave activity has been suggested to be connected to surface temperature anomalies, but some key processes are not well understood. Using observations, we show that the stratospheric events featuring weaker-than-normal wave activity are associated with increased North American (NA) cold extreme risks before and near the event onset, accompanied by less frequent atmospheric river (AR) events on the west coast of the United States. Strong stratospheric wave events, on the other hand, exhibit a tropospheric weather regime transition. Further, they are preceded by NA warm anomalies and increased AR frequency over the west coast, followed by increased risks of NA cold extremes and north-shifted ARs over the Atlantic. Moreover, these links between the stratosphere and troposphere are attributed to the vertical structure of wave coupling. Weak wave events show a wave structure of westward tilt with increasing altitudes, while strong wave events feature a shift from westward tilt to eastward tilt during their life cycle. This wave phase shift indicates vertical wave coupling and likely regional planetary wave reflection. Further examinations of CMIP6 models show that models with a degraded representation of stratospheric wave structure exhibit biases in the troposphere during strong wave events. Specifically, models with a stratospheric ridge weaker than the reanalysis exhibit a weaker tropospheric signal. Our findings suggest that the vertical coupling of extreme stratospheric wave activity should be evaluated in the model representation of stratosphere-troposphere coupling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Simulating the Trajectory and Biomass Growth of Free-Floating Macroalgal Cultivation Platforms along the U.S. West Coast

Trajectory tracking and macroalgal growth models were coupled to support a novel macroalgae-harvesting concept known as the Nautical Off-shore Macroalgal Autonomous Device (NOMAD). The NOMAD consists of 5 km long carbon-fiber longlines that are seeded and free float southward along the U.S. West Coast for approximately 3 months before harvesting off the California coast, taking advantage of favorable environmental conditions. The trajectory and macroalgal growth models were applied to answer planning questions pertinent to the techno-economic analysis such as identifying the preferred release location, approximate pathway, timing until harvest, and estimated growth. Trajectories were determined with the General NOAA Operational Modeling Environment (GNOME) model, using 11 years of current and wind data, determining probabilities by running nearly 40,000 Monte Carlo simulations varying the start time and location. An accompanying macroalgal growth model was used to estimate the growth of macroalgae based on the trajectory tracks and environmental forcing products, including light, temperature and nutrients. Model results show that NOMAD lines transit south in the months of April to September due to seasonal currents, taking approximately 3 months to reach Southern California. During transit, NOMAD lines are dispersed but typically avoid beaching or passing through marine sanctuaries. NOMAD lines can yield up to 30 kg wet weight per meter of cultivation line.

09 BIOMASS FUELS↗
DOE OSTI2024

Unique impacts of strong and westward-extended western Pacific subtropical high on ozone pollution over eastern China

As a subtropical anticyclonic high-pressure system that typically forms over the northwestern Pacific Ocean in summer, the Western Pacific subtropical high (WPSH) affects meteorological conditions and ozone pollution in China. The relationship between maximum daily 8-h average ozone (MDA8 O 3 ) concentrations and the extremely strong and westward-extended WPSH occurred in 2022 is investigated using observations, reanalysis data and GEOS-Chem model simulations. The intensity of WPSH has a significant positive correlation with MDA8 O3 over southern China during July-August in 2022, with a correlation coefficient of +0.44, but the correlation is negative (–0.40) in northern China. During the strong WPSH days, MDA8 O 3 increased by 16.5µgm -3 (16.4% relative to July-August average) over southern China and decreased by 19.0µgm -3 (14.5%) in northern China compared to the weak WPSH days. The unique dipole pattern in the relationship between ozone levels and the WPSH in 2022 exhibited a contrast to that during 2015–2021. The difference is primarily due to the extremely strong WPSH intensity and its unusual westward expansion in 2022. In this case, an anomalous anticyclone at 500 hPa dominates over southern China, which creates conditions conducive for ozone formation and accumulation. The anticyclone weakened horizontal winds and reduced the dispersion of ozone, alongside a high temperature and low relative humidity, which favored the chemical production of ozone. In contrast, abnormal northerly winds enhanced ozone diffusion in northern China and the low temperature reduced ozone chemical production. Here, this study reveals the mechanism for the significant impact of strong and westward-extended WPSH on ozone concentrations over China, emphasizing the role of the WPSH location in modulating meteorology and ozone levels.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Climatology of Severe Local Storm Environments and Synoptic-Scale Features over North America in ERA5 Reanalysis and CAM6 Simulation

Severe local storm (SLS) activity is known to occur within specific thermodynamic and kinematic environments. These environments are commonly associated with key synoptic-scale features—including southerly Great Plains low-level jets, drylines, elevated mixed layers, and extratropical cyclones—that link the large-scale climate to SLS environments. This work analyzes spatiotemporal distributions of both extreme values of SLS environmental parameters and synoptic-scale features in the ERA5 reanalysis and in the Community Atmosphere Model, version 6 (CAM6), historical simulation during 1980–2014 over North America. Compared to radiosondes, ERA5 successfully reproduces SLS environments, with strong spatiotemporal correlations and low biases, especially over the Great Plains. Both ERA5 and CAM6 reproduce the climatology of SLS environments over the central United States as well as its strong seasonal and diurnal cycles. ERA5 and CAM6 also reproduce the climatological occurrence of the synoptic-scale features, with the distribution pattern similar to that of SLS environments. Compared to ERA5, CAM6 exhibits a high bias in convective available potential energy over the eastern United States primarily due to a high bias in surface moisture and, to a lesser extent, storm-relative helicity due to enhanced low-level winds. Additionally, composite analysis indicates consistent synoptic anomaly patterns favorable for significant SLS environments over much of the eastern half of the United States in both ERA5 and CAM6, though the pattern differs for the southeastern United States. Overall, our results indicate that both ERA5 and CAM6 are capable of reproducing SLS environments as well as the synoptic-scale features and transient events that generate them.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Identifying atmospheric rivers and their poleward latent heat transport with generalizable neural networks: ARCNNv1

Abstract. Atmospheric rivers (ARs) are extreme weather events that can alleviate drought or cause billions of US dollars in flood damage. By transporting significant amounts of latent energy towards the poles, they are crucial to maintaining the climate system's energy balance. Since there is no first-principle definition of an AR grounded in geophysical fluid mechanics, AR identification is currently performed by a multitude of expert-defined, threshold-based algorithms. The variety of AR detection algorithms has introduced uncertainty into the study of ARs, and the thresholds of the algorithms may not generalize to new climate datasets and resolutions. We train convolutional neural networks (CNNs) to detect ARs while representing this uncertainty; we name these models ARCNNs. To detect ARs without requiring new labeled data and labor-intensive AR detection campaigns, we present a semi-supervised learning framework based on image style transfer. This framework generalizes ARCNNs across climate datasets and input fields. Using idealized and realistic numerical models, together with observations, we assess the performance of the ARCNNs. We test the ARCNNs in an idealized simulation of a shallow-water fluid in which nearly all the tracer transport can be attributed to AR-like filamentary structures. In reanalysis and a high-resolution climate model, we use ARCNNs to calculate the contribution of ARs to meridional latent heat transport, and we demonstrate that this quantity varies considerably due to AR detection uncertainty.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Influence of Lake Ice Biases in Reanalysis Data on Downscaled Climate Simulations over the Great Lakes Region

This data package contains observation-based and model-simulated datasets (all provided in NetCDF format) for evaluating how wintertime lake-ice representation affects regional weather and climate over the Laurentian Great Lakes (freshwater lake ecosystem) during the high–ice-cover winter of 2009. The observational component includes: (1) Stage IV gridded precipitation at 4 km, hourly resolution for January–February 2009 over the Great Lakes region (radar–gauge multisensor precipitation analyses); (2) Great Lakes Surface Environmental Analysis (GLSEA) satellite-derived lake-ice coverage at 1.3 km, daily resolution for the 2009 winter months, providing ice coverage over Lakes Superior, Michigan, Huron, Erie, and Ontario; and (3) in situ measurements at the Standard Rock site on Lake Superior from the Great Lakes Evaporation Network (GLEN) at hourly resolution, including near-surface atmospheric variables and sensible and latent heat fluxes (air–lake exchange) at a fixed point location. The modeling component provides corresponding fields from two simulations, both archived at 4 km, hourly resolution: a standalone Weather Research Forecasting model (WRF) run driven by the European Centre for Medium-Range Weather Forecasts (ECMWF) Reanalysis 5 (ERA5), and a two-way coupled model using WRF and the Finite Volume Community Ocean Model (WRF-FVCOM, a 3-D hydrodynamic lake model). These outputs include variables relevant to air–lake interaction and lake-effect processes (e.g., near-surface temperature, humidity, wind, precipitation, and surface turbulent fluxes), enabling direct comparison with the observational datasets. Users can analyze and visualize these NetCDF files with common tools such as Python (e.g., xarray, netCDF4, numpy, pandas), NCO/CDO, Panoply, or ncview; NetCDF variables can also be converted to other formats (e.g., CSV, GeoTIFF) using these utilities.

EARTH SCIENCE > ATMOSPHERE > ATMOSPHERIC TEMPERATU↗
DOE OSTI2023

The Misrepresentation of the Southern African Easterly Jet in Models and Its Implications for Aerosol, Clouds, and Precipitation Distributions

The southern African easterly jet (AEJ-S) is an important midtropospheric feature critical to understanding the tropical convective system over central Africa and the aerosol–cloud interactions over the southeast Atlantic Ocean. However, it remains unclear how well models represent the AEJ-S and its influence on aerosol transport, clouds, and precipitation distribution. Here, in this work, we use ground- and satellite-based observations and reanalysis datasets to assess the representation of AEJ-S in the Coupled Model Intercomparison Project phase 6 (CMIP6) models between September and October during the peak of midtropospheric winds, aerosol transport, clouds, and precipitation. We find that most CMIP6 models have difficulty accurately simulating the strength, position, and spatial distribution of the AEJ-S. Specifically, the AEJ-S is relatively weaker and at a slightly lower altitude in the ensemble of CMIP6 models than represented by observation and reanalysis datasets. To assess the influence of the misrepresented the AEJ-S on CMIP6-simulated aerosol, clouds, and precipitation distributions, we performed composite analyses using models with low and high biases based on the estimates of their midtropospheric easterly wind speed. We find that the misrepresentation of the AEJ-S in CMIP6 models is associated with the overestimation of clouds and precipitation over central Africa, the underestimation of clouds over the southeast Atlantic Ocean, and the limitation of aerosol transport over the continent or the deviation of its spatial distribution from the typical zonal transport over the Atlantic Ocean. Because aerosols, clouds, and precipitation are important components of the regional climate system, we conclude that accurate representation of the AEJ-S is essential over central Africa and the southeast Atlantic Ocean.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A spatial evaluation of Arctic sea ice and regional limitations in CMIP6 historical simulations

The Arctic sea ice response to a warming climate is assessed in a subset of models participating in Phase 6 of the Coupled Model Intercomparison Project (CMIP6), using several metrics in comparison with satellite observations and results from the Pan-Arctic Ice Ocean Modeling and Assimilation System and the Regional Arctic System Model. Our study examines the historical representation of sea ice extent, volume, and thickness using spatial analysis metrics, such as the integrated ice-edge error, Brier score, and Spatial Probability Score. We find that the CMIP6 multi-model mean captures the mean annual cycle and 1979-2014 sea ice trends remarkably well. However, individual models experience a wide range of uncertainty in the spatial distribution of sea ice when compared against satellite measurements and reanalysis data. Our metrics expose common and individual regional model biases, which sea ice temporal analyses alone do not capture. We identify large ice edge and ice thickness errors in Arctic sub-regions, implying possible model specific limitations in or lack of representation of some key physical processes. We postulate that many of them could be related to the oceanic forcing, especially in the marginal and shelf seas, where seasonal sea ice changes are not adequately simulated. We therefore conclude that an individual model’s ability to represent the observed/reanalysis spatial distribution still remains a challenge. We propose the spatial analysis metrics as useful tools to diagnose model limitations, narrow down possible processes affecting them, and guide future model improvements critical to the representation and projections of Arctic climate change.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Insights From Dayflow: A Historical Streamflow Reanalysis Dataset for the Conterminous United States

Abstract Reconstructed historical streamflow time series can supplement limited streamflow gauge observations. However, there are common challenges of typical modeling approaches: process‐based hydrologic models can be data/computation‐intensive, and statistics‐based models can be region/stream‐specific. Here we present a nationally scalable modeling framework integrating the simulated runoff from the Variable Infiltration Capacity (VIC) model with the Routing Application for Parallel computatIon of Discharge (RAPID) routing model leveraging high‐performance computing. We demonstrate an efficient method of assimilating streamflow at US Geological Survey (USGS) streamflow monitoring sites using a simple hierarchical approach in the VIC‐RAPID framework. The result is a reconstructed 36‐year (1980–2015) daily and monthly streamflow dataset (Dayflow) at ∼2.7 million NHDPlusV2 stream reaches in the conterminous US (CONUS). We perform a comprehensive evaluation at 7,526 USGS sites and characterize their error statistics. The results demonstrate that 49% of the USGS sites demonstrate Kling–Gupta Efficiency (KGE) > 0.5 and 58% of the sites show percentage bias within ±20% for the daily naturalized streamflow. Streamflow data assimilation across CONUS shows an overall improvement over naturalized streamflow, notably in the western semiarid‐to‐arid regions. Comparison to other national and global streamflow reanalysis datasets such as the National Water Model and Global Reach‐scale A priori Discharge Estimates for SWOT demonstrates improved KGE, reduced bias, and directions for Dayflow improvements. Investigations of error statistics with key hydrologic, hydroclimatic, and geomorphologic basin characteristics reveal region‐specific patterns which may help improve future framework applications. Overall, Dayflow may enable a better understanding of hydrologic conditions in a changing environment, especially in locations currently not represented by streamflow monitoring networks.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

2023 National Offshore Wind data set (NOW-23)

The 2023 National Offshore Wind data set (NOW-23) is the latest wind resource data set for offshore regions in the United States, which supersedes, for its offshore component, the Wind Integration National Dataset (WIND) Toolkit, which was published about a decade ago and is currently one of the primary resources for stakeholders conducting wind resource assessments in the continental United States. The NOW-23 data set was produced using the Weather Research and Forecasting Model (WRF) version 4.2.1. A regional approach was used: for each offshore region, the WRF setup was selected based on validation against available observations. The WRF model was initialized with the European Centre for Medium Range Weather Forecasts 5 Reanalysis (ERA-5) data set, using a 6-hour refresh rate. The model is configured with an initial horizontal grid spacing of 6 km and an internal nested domain that refined the spatial resolution to 2 km. The model is run with 61 vertical levels, with 12 levels in the lower 300m of the atmosphere, stretching from 5 m to 45 m in height. The MYNN planetary boundary layer and surface layer schemes were used the North Atlantic, Mid Atlantic, Great Lakes, Hawaii, and North Pacific regions. On the other hand, using the YSU planetary boundary layer and MM5 surface layer schemes resulted in a better skill in the South Atlantic, Gulf of Mexico, and South Pacific regions. A more detailed description of the WRF model setup can be found in the WRF namelist files linked at the bottom of this page. For all regions, the NOW-23 data set coverage starts on January 1, 2000. For Hawaii and the North Pacific regions, NOW-23 goes until December 31, 2019. For the South Pacific region, the model goes until 31 December, 2022. For all other regions, the model covers until December 31, 2020. Outputs are available at 5 minute resolution, and for all regions we have also included output files at hourly resolution. The NOW-23 data are provided here as HDF5 files. Examples of how to use the HSDS Service to Access the NOW-23 files are linked below. A list of the variables included in the NOW-23 files is also linked below. No filters have been applied to the raw WRF output.

Array↗
DOE OSTI2025

A Multivariate Space‐Time Dynamic Model for Characterizing the Atmospheric Impacts Following the Mt. Pinatubo Eruption

The June 1991 Mt. Pinatubo eruption resulted in a massive increase of sulfate aerosols in the atmosphere, absorbing radiation and leading to global changes in surface and stratospheric temperatures. A volcanic eruption of this magnitude serves as a natural analog for stratospheric aerosol injection, a proposed solar radiation modification method to combat a warming climate. The impacts of such an event are multifaceted and region-specific. Our goal is to characterize the multivariate and dynamic nature of the atmospheric impacts following the Mt. Pinatubo eruption. We developed a multivariate space-time dynamic linear model to understand the full extent of the spatially- and temporally-varying impacts. Specifically, spatial variation is modeled using a flexible set of basis functions for which the basis coefficients are allowed to vary in time through a vector autoregressive (VAR) structure. This novel model is cast in a Dynamic Linear Model (DLM) framework and estimated via a customized MCMC approach. We demonstrate how the model quantifies the relationships between key atmospheric parameters prior to and following the Mt. Pinatubo eruption with reanalysis data from MERRA-2 and highlight when such a model is advantageous over univariate models.

Dynamic Linear Model↗
DOE OSTI2022

Aridity-Dependent Land Surface Skin Temperature Biases in CMIP5/6

Land surface skin temperature, a critical indicator of climate change, connects the water and energy cycles between the land and the atmosphere. Here, in this study, we evaluate the simulations of land surface skin temperature from the Coupled Model Intercomparison Project Phase 5 (CMIP5) and CMIP6 models with satellite-based datasets and reanalysis. We find systematic cold skin temperature biases over arid regions in CMIP5/CMIP6 simulations. Over arid and semi-arid regions, latent heat biases drive skin temperature biases by evaporative cooling. Over humid regions, surface downward shortwave and albedo biases are relatively more critical. Spatial patterns of biases remain similar in the latest CMIP6 simulations, suggesting systematic biases in land-atmosphere interactions. These biases need to be corrected or considered while using models for future projections.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

E3SMv3 ‘Nephele’ Perturbed Parameter Ensemble Value-Added Product for ENA, SGP, NSA, and Global Means

Processed output from the E3SMv3 Nephele perturbed parameter ensemble (PPE). This product contains global means of present-day (using 2010 aerosol forcings) and preindustrial (using 1850 aerosol forcings) monthly and daily global-mean quantities, as well as hourly mean quantities at the E3SMv3 model grid cell closest to the ENA, SGP, and NSA ARM sites. In this PPE, the E3SMv3 model was run for 27 months in an atmosphere-only configuration with free tropospheric winds nudged to reanalysis. Each ensemble member represents a unique combination of 25 model parameter values (in the microphysics, convective microphysics, and aerosol schemes) that were perturbed randomly and simultaneously. The E3SMv3 default configuration is included as member 0. See https://doi.org/10.22541/essoar.174907165.57104591/v1 (in review) for more information on the Nephele PPE and the perturbed parameters.

Model atmospheric temperature↗
DOE OSTI2023

Implementing and Improving CBMZ-MAM3 Chemistry and Aerosol Modules in the Regional Climate Model WRF-CAM5: An Evaluation over the Western US and Eastern North Pacific

The representation of aerosols in climate-chemistry models is important for air quality and climate change research, but it can require significant computational resources. The objective of this study was to improve the representation of aerosols in climate–chemistry models, specifically in the carbon bond mechanism, version Z (CBMZ), and modal aerosol modules with three lognormal modes (MAM3) in the WRF-CAM5 model. The study aimed to enhance the model’s chemistry capabilities by incorporating biomass burning emissions, establishing a conversion mechanism between volatile organic compounds (VOCs) and secondary organic carbons (SOCs), and evaluating its performance against observational benchmarks. The results of the study demonstrated the effectiveness of the enhanced chemistry capabilities in the WRF-CAM5 model. Six simulations were conducted over the western U.S. and northeastern Pacific region, comparing the model’s performance with observational benchmarks such as reanalysis, ground-based, and satellite data. The findings revealed a significant reduction in root-mean-square errors (RMSE) for surface concentrations of black carbon (BC) and organic carbon (OC). Specifically, the model exhibited a 31% reduction in RMSE for BC concentrations and a 58% reduction in RMSE for OC concentrations. These outcomes underscored the importance of accurate aerosol representation in climate-chemistry models and emphasized the potential for improving simulation accuracy and reducing errors through the incorporation of enhanced chemistry modules in such models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Reducing Southern Ocean Shortwave Radiation Errors in the ERA5 Reanalysis with Machine Learning and 25 Years of Surface Observations

Earth system models struggle to simulate clouds and their radiative effects over the Southern Ocean, partly due to a lack of measurements and targeted cloud microphysics knowledge. We have evaluated biases of downwelling shortwave radiation in the ERA5 climate reanalysis using 25 years (1995–2019) of summertime surface measurements, collected on the Research and Supply Vessel (RSV) Aurora Australis, the Research Vessel (R/V) Investigator, and at Macquarie Island. During October–March daylight hours, the ERA5 simulation of SW down exhibited large errors (mean bias = 54 W m -2 , mean absolute error = 82 W m -2 , root-mean-square error = 132 W m -2 , and R 2 = 0.71). To determine whether we could improve these statistics, we bypassed ERA5’s radiative transfer model for SW down with machine learning–based models using a number of ERA5’s gridscale meteorological variables as predictors. These models were trained and tested with the surface measurements of SW down using a 10-fold shuffle split. An extreme gradient boosting (XGBoost) and a random forest–based model setup had the best performance relative to ERA5, both with a near complete reduction of the mean bias error, a decrease in the mean absolute error and root-mean-square error by 25% ± 3%, and an increase in the R 2 value of 5% ± 1% over the 10 splits. Large improvements occurred at higher latitudes and cyclone cold sectors, where ERA5 performed most poorly. We further interpret our methods using Shapley additive explanations. Our results indicate that data-driven techniques could have an important role in simulating surface radiation fluxes and in improving reanalysis products.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Implementation of an Orographic Drag Scheme Considering Orographic Anisotropy in All Flow Directions in the Earth System Model CAS‐ESM 2.0

Abstract A reasonable representation of orographic anisotropy in earth system models is vital for improving weather and climate modeling. In this study, we implement the orographic drag scheme, including 3‐D orographic anisotropy (3D‐AFD), into the Chinese Academy of Sciences Earth System Model version 2 (CAS‐ESM 2.0). Three groups of simulations named sensitivity run, medium‐range forecast, and seasonal forecast respectively are conducted using the updated CAS‐ESM model and validated against station observation and reanalysis data. These simulations are run using the three schemes—3D‐AFD scheme, the 3D orographic anisotropy scheme for the eight‐direction (3D‐8x), and the 2D isotropic scheme (2D)—to compare their performance in CAS‐ESM 2.0. Sensitivity runs show that the 3D‐AFD provide more continuous transition of calculated drag as function of wind direction than the 3D‐8x, while the drag calculated using the 2D scheme show no variation with change of wind direction. Enhanced drag in the medium range and seasonal forecast using the updated CAS‐ESM alleviates part of the winter wind speed bias over the Tibetan Plateau (TP) and the cold bias over TP and the Siberian polar region. It is shown that the 3D‐AFD scheme alleviates more bias than that of the 3D‐8x scheme (by wind speed reduction of 1 ∼ 2 m/s and temperature of 1 ∼ 2 K) especially in the seasonal forecast. The results suggest that reasonable representation of the orographic anisotropy is important in climate modeling.

Xie, Jinbo↗
DOE OSTI2026

Improving thermodynamic nudging in the E3SM Atmosphere Model version 2 (EAMv2): strategy and hindcast skills on weather systems

Nudging techniques are commonly employed to constrain atmospheric simulations toward observed states, facilitating model evaluation and sensitivity studies. However, if applied improperly – particularly to thermodynamic variables such as temperature and humidity – nudging can distort physical processes and introduce spurious biases, undermining the credibility of the simulations. This study presents an improved nudging implementation that applies vertically modulated tendencies to reduce adverse impacts on model physics. The framework is tested in version 2 of the Energy Exascale Earth System Model (EAMv2) using a suite of hindcast simulations nudged toward ERA5 reanalysis. We systematically evaluate the individual and combined effects of nudging wind, temperature, and humidity fields on the model's ability to represent large-scale atmospheric states and high-impact weather systems. Results show that the revised strategy – particularly when nudging temperature and humidity at selected levels – enhances hindcast skill by improving agreement with ERA5 without degrading the hydrological cycle or precipitation processes. Additional improvements in surface temperature, outgoing longwave radiation, and precipitation biases are achieved through targeted nudging of land-surface variables. The proposed approach strengthens the representation of large-scale conditions relevant to tropical cyclones, atmospheric rivers, and extratropical cyclones in the low-resolution EAMv2. These findings demonstrate that carefully designed thermodynamic nudging, especially of temperature and humidity, improves the realism of constrained simulations and broadens the utility of nudged EAMv2 for atmospheric modeling, machine learning, and high-impact weather research.

Atmospheric river↗
DOE OSTI2025

Development and evaluation of a new 4DEnVar-based weakly coupled ocean data assimilation system in E3SMv2

The development, implementation, and evaluation of a new weakly coupled ocean data assimilation (WCODA) system for the fully coupled Energy Exascale Earth System Model version 2 (E3SMv2) utilizing the four-dimensional ensemble variational (4DEnVar) method are presented in this study. The 4DEnVar method, based on the dimension-reduced projection four-dimensional variational (DRP-4DVar) approach, replaces the adjoint model with the ensemble technique, thereby reducing computational demands. Monthly mean ocean temperature and salinity data from the EN4.2.1 reanalysis are integrated into the ocean component of E3SMv2 from 1950 to 2021 with the goal of providing realistic initial conditions for decadal predictions and predictability studies. The performance of the WCODA system is assessed using various metrics, including the reduction rate of the cost function, root mean square error (RMSE) differences, correlation differences, and model biases. Results indicate that the WCODA system effectively assimilates the reanalysis data into the climate model, consistently achieving negative reduction rates of the cost function and notable improvements in RMSE and correlation across various ocean layers and regions. Significant enhancements are observed in the upper ocean layers across the majority of global ocean regions, particularly in the north Atlantic, north Pacific, and Indian Ocean. Model biases in sea surface temperature and salinity are also substantially reduced. For sea surface temperature, cold biases in the north Pacific and north Atlantic are diminished by about 1–2 °C, and warm biases in the Southern Ocean are corrected by approximately 1.5–2.5 °C. In terms of salinity, improvements are observed with bias reductions of about 0.5–1 psu in the north Atlantic and north Pacific and up to 1.5 psu in parts of the Southern Ocean. The ultimate goal of the WCODA system is to advance the predictive capabilities of E3SM for subseasonal to decadal climate predictions, thereby supporting research on strategic energy-sector policies and planning.

54 ENVIRONMENTAL SCIENCES↗