Search NASA⌕ Search

SEARCH · Search NASA

Results for “Atmospheric forcing”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Long‐Term Large‐Scale Atmospheric Forcing Data From Three‐Dimensional Constrained Variational Analysis for the ARM SGP Site

Here, this study presents a long‐term three‐dimensional large‐scale forcing data set (VARANAL3D) derived from the three‐dimensional constrained variational analysis (3DCVA) method at the Atmospheric Radiation Measurement (ARM) program Southern Great Plains (SGP) site from 2004 to 2018. Building on the same input data sets as the conventional continuous forcing data set (VARANAL), VARANAL3D maintains overall consistency in domain‐averaged fields while introducing spatial variability, offering critical insights into the influence of mesoscale synoptic systems on cloud‐related processes. Evaluations are conducted across four cloud and precipitation regimes: Clear‐sky, Shallow‐clouds, Afternoon‐precipitation, and Nocturnal‐precipitation, presenting high consistency of the domain‐mean forcing data sets while emphasizing the role of subdomain forcing variability particularly in precipitating regimes. Single column model (SCM) simulations demonstrate that subdomain VARANAL3D forcing improves cloud and precipitation representation, with the ensemble outperforming domain‐mean forcing in three cloudy and precipitating regimes. Overall, these results highlight VARANAL3D's value for investigating the impacts of spatial variability of large‐scale forcing on atmospheric processes. The VARANAL3D data set provides new opportunities for evaluating model physics, advancing the development of scale‐aware parameterizations and deepening our understanding of cloud and precipitation dynamics.

Environmental sciences↗
DOE OSTI2022

Uncertainty in land carbon budget simulated by terrestrial biosphere models: the role of atmospheric forcing

Global estimates of the land carbon sink are often based on simulations by terrestrial biosphere models (TBMs). The use of a large number of models that differ in their underlying hypotheses, structure and parameters is one way to assess the uncertainty in the historical land carbon sink. Here we show that the atmospheric forcing datasets used to drive these TBMs represent a significant source of uncertainty that is currently not systematically accounted for in land carbon cycle evaluations. We present results from three TBMs each forced with three different historical atmospheric forcing reconstructions over the period 1850–2015. We perform an analysis of variance to quantify the relative uncertainty in carbon fluxes arising from the models themselves, atmospheric forcing, and model-forcing interactions. We find that atmospheric forcing in this set of simulations plays a dominant role on uncertainties in global gross primary productivity (GPP) (75% of variability) and autotrophic respiration (90%), and a significant but reduced role on net primary productivity and heterotrophic respiration (30%). Atmospheric forcing is the dominant driver (52%) of variability for the net ecosystem exchange flux, defined as the difference between GPP and respiration (both autotrophic and heterotrophic respiration). In contrast, for wildfire-driven carbon emissions model uncertainties dominate and, as a result, model uncertainties dominate for net ecosystem productivity. At regional scales, the contribution of atmospheric forcing to uncertainty shows a very heterogeneous pattern and is smaller on average than at the global scale. We find that this difference in the relative importance of forcing uncertainty between global and regional scales is related to large differences in regional model flux estimates, which partially offset each other when integrated globally, while the flux differences driven by forcing are mainly consistent across the world and therefore add up to a larger fractional contribution to global uncertainty.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Systematic Evaluation of Atmospheric Forcing, Surface Datasets, and Mesh Effects on Kilometer-Scale Land Surface and River Modeling

Earth system models are advancing toward kilometer-scale resolution to capture local climate impacts and extremes. High-resolution land and river modeling depends on multiple factors, including mesh, surface datasets, and atmospheric forcing, but their relative effects at kilometer scales remain unquantified. We evaluated five Energy Exascale Earth System Model land and river configurations over the Mid-Atlantic region using two mesh (1/8° structured versus variable-resolution unstructured mesh), two surface datasets (default versus newly developed), and three atmospheric forcings (NLDAS2, MSWX, GSWP). Evaluation against satellite, reanalysis, and in situ benchmarks across water, energy, and carbon cycles quantifies how these factors affect model performance. Forcing selection produces the largest bias reductions (12-99% across variables), followed by surface datasets (7-75%) and mesh (up to 21%). Forcing effects vary by variable, with MSWX reducing biases for snow water equivalent, evapotranspiration, albedo, temperature, and gross primary productivity, GSWP for snow cover and runoff, and NLDAS for soil moisture and streamflow. The use of newly developed surface datasets improves gross primary productivity (58% bias reduction) and evapotranspiration but increase soil moisture and albedo biases due to current modeling limitations. Variable-resolution unstructured mesh improves the simulation of small-basin streamflow through better capturing drainage networks, though mesh minimally affects other land variables. These findings provide important guidance for high-resolution modeling development and actionable science.

Land and River modeling↗
DOE OSTI2024

Impacts of Topography-Based Subgrid Scheme and Downscaling of Atmospheric Forcing on Modeling Land Surface Processes in the Conterminous US

The effects of small-scale topography-induced land surface heterogeneity are not well represented in current Earth System Models (ESMs). In this study, a new topography-based subgrid structure referred to as topographic units (TGU) designed to better capture subgrid topographic effects, and methods to downscale atmospheric forcing to the land TGUs have been implemented in the Energy Exascale Earth System Model (E3SM) Land Model (ELM). Effects of the subgrid scheme and downscaling methods on ELM simulated land surface processes are evaluated over the conterminous United States (CONUS). For this purpose, ELM simulations are performed using two configurations without (NoD ELM) and with (D ELM) downscaling, both using TGUs derived for the 0.5-degree grids and the same land surface parameters. Simulations using the two ELM configurations are compared over the CONUS domain, regional levels, and at observational sites (e.g., SNOTEL). The CONUS-level results suggest that D ELM simulates more snowfall and snow water equivalent (SWE), higher runoff, and less ET during spring and summer. Regional-level results suggest more pronounced impacts of downscaling over regions dominated by higher elevation TGUs and regions with maximum precipitation occurring during cool seasons. Results at the SNOTEL sites suggest that D ELM has superior capability of reproducing the observed SWE at 83% of the sites, with more pronounced performance over topographically heterogeneous TGUs with their maximum precipitation occurring during cool seasons. The results highlight the importance of improving representation of small-scale surface heterogeneity in ESMs and motivate future research to understand their effects on land-atmosphere interactions, streamflow, and water resources management over mountainous regions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Improving Simulations of Vegetation Dynamics over the Tibetan Plateau: Role of Atmospheric Forcing Data and Spatial Resolution

The efficacy of vegetation dynamics simulations in offline land surface models (LSMs) largely depends on the quality and spatial resolution of meteorological forcing data. In this study, the Princeton Global Meteorological Forcing Data (PMFD) and the high spatial resolution and upscaled China Meteorological Forcing Data (CMFD) were used to drive the Simplified Simple Biosphere model version 4/Top-down Representation of Interactive Foliage and Flora Including Dynamics (SSiB4/TRIFFID) and investigate how meteorological forcing datasets with different spatial resolutions affect simulations over the Tibetan Plateau (TP), a region with complex topography and sparse observations. By comparing the monthly Leaf Area Index (LAI) and Gross Primary Production (GPP) against observations, we found that SSiB4/TRIFFID driven by upscaled CMFD improved the performance in simulating the spatial distributions of LAI and GPP over the TP, reducing RMSEs by 24.3% and 20.5%, respectively. The multi-year averaged GPP decreased from 364.68 gC m –2 yr –1 to 241.21 gC m –2 yr –1 with the percentage bias dropping from 50.2% to –1.7%. When using the high spatial resolution CMFD, the RMSEs of the spatial distributions of LAI and GPP simulations were further reduced by 7.5% and 9.5%, respectively. This study highlights the importance of more realistic and high-resolution forcing data in simulating vegetation growth and carbon exchange between the atmosphere and biosphere over the TP.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Ocean fronts and eddies force atmospheric rivers and heavy precipitation in western North America

Atmospheric rivers (ARs) are responsible for over 90% of poleward water vapor transport in the mid-latitudes and can produce extreme precipitation when making landfall. However, weather and climate models still have difficulty simulating and predicting landfalling ARs and associated extreme precipitation, highlighting the need to better understand AR dynamics. Here, using high-resolution climate models and observations, we demonstrate that mesoscale sea-surface temperature (SST) anomalies along the Kuroshio Extension can exert a remote influence on landfalling ARs and related heavy precipitation along the west coast of North America. Inclusion of mesoscale SST forcing in the simulations results in approximately a 40% increase in landfalling ARs and up to a 30% increase in heavy precipitation in mountainous regions and this remote impact occurs on two-week time scales. The asymmetrical response of the atmosphere to warm vs. cold mesoscale SSTs over the eddy-rich Kuroshio Extension region is proposed as a forcing mechanism that results in a net increase of moisture flux above the planetary boundary layer, prompting AR genesis via enhancing moisture transport into extratropical cyclones in the presence of mesoscale SST forcing.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Strong Red Noise Ocean Forcing on Atlantic Multidecadal Variability Assessed from Surface Heat Flux: Theory and Application

The role of ocean forcing on Atlantic multidecadal variability (AMV) is assessed from the (downward) heat flux–SST relation in the framework of a new stochastic climate theory forced by red noise ocean forcing. Previous studies suggested that atmospheric forcing drives SST variability from monthly to interannual time scales, with a positive heat flux–SST correlation, while heat flux induced by ocean processes can drive SST variability at decadal and longer time scales, with a negative heat flux–SST correlation. Here, first, we develop a theory to show how the sign of heat flux–SST correlation is affected by atmospheric and oceanic forcing with time scale. In particular, a red noise ocean forcing is necessary for the sign reversal of heat flux–SST correlation. Furthermore, this sign reversal can be detected equivalently in three approaches: the low-pass correlation at lag zero, the unfiltered correlation at long (heat flux) lead, and the real part of the heat flux–SST coherence. Second, we develop a new scheme in combination with the theory to assess the magnitude and time scale of the red noise ocean forcing for AMV in the GFDL SPEAR model (Seamless System for Prediction and Earth System Research) and observations. In both the model and observations, the ocean forcing on AMV is in general comparable with the atmospheric forcing, with a 90% probability greater than the atmospheric forcing in observations. In contrast to the white noise atmospheric forcing, the ocean forcing has a persistence time comparable or longer than a year, much longer than the SST persistence of ~3 months. This slow ocean forcing is associated implicitly with slow subsurface ocean dynamics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

More Realistic Intermediate Depth Dry Firn Densification in the Energy Exascale Earth System Model (E3SM)

Abstract Earth system models account for seasonal snow cover, but many do not accommodate the deeper snowpack on ice sheets (aka firn) that slowly transforms to ice under accumulating snowfall. To accommodate and resolve firn depths of up to 60 m in the Energy Exascale Earth System Model's land surface model (ELM), we add 11 layers to its snowpack and evaluate three dry snow compaction equations in multi‐century simulations. After comparing results from ELM simulations (forced with atmospheric reanalysis) with empirical data, we find that implementing into ELM a two‐stage firn densification model produces more accurate dry firn densities at intermediate depths of 20–60 m. Compared to modeling firn using the equations in the (12 layer) Community Land Model (version 5), switching to the two‐stage firn densification model (with 16 layers) significantly decreases root‐mean‐square errors in upper 60 m dry firn densities by an average of 41 kg m −3 (31%). Simulations with three different firn density parameterizations show that the two‐stage firn densification model should be used for applications that prioritize accurate upper 60 m firn air content (FAC) in regions where the mean annual surface temperature is greater than roughly −31°C. Because snow metamorphism, firn density, and FAC are major components in modeling ice sheet surface albedo, melt water retention, and climatic mass balance, these developments advance broader efforts to simulate the response of land ice to atmospheric forcing in Earth system models.

58 GEOSCIENCES↗
DOE OSTI2026

CLM5 Simulations of Soil Moisture and Gross Primary Productivity for CONUS at 0.125 degrees

This dataset provides 0.125-degree gridded simulations of soil moisture and gross primary productivity (GPP) for the Contiguous United States (CONUS), generated using the Community Land Model version 5 (CLM5) with the biogeochemistry module enabled. The data covers a historical baseline (1980-2015) and mid-century future projections (2020-2055). Future projections are organized into two sets of scenarios to distinguish the impacts of different drivers: (1) Atmospheric Only (ATM): These scenarios apply future atmospheric forcings while holding land use and land cover (LULC) at historical baseline levels. The atmospheric forcings represent moderately versus severely hotter/drier atmospheric conditions (dynamically downscaled perturbed thermodynamics simulations based on CMIP6 SSP245 and SSP585 warming signals), each with cooler versus hotter Earth System Model temperature sensitivity instantiations. These scenarios are identified in the folder names as rcp45_cooler_near, rcp45_hotter_near, rcp85_cooler_near, and rcp85_hotter_near. (2) Coupled Atmospheric and Land-Use (LAND+ATM): These scenarios apply future atmospheric forcing together with future LULC by pairing atmospheric pathways with lower versus higher population/economic growth scenarios representing Shared Socioeconomic Pathways 3 and 5 (SSP3 and SSP5). These scenarios are identified in the file names as ssp3_rcp45_cooler_near, ssp3_rcp45_hotter_near, ssp5_rcp85_cooler_near, and ssp5_rcp85_hotter_near. Please refer to the "README_first.md" file for detailed information on file structure, variables, units, and data formats.

drought↗
DOE OSTI2026

CLM5 Simulations of Soil Moisture and Grain Carbon for CONUS at 0.125 degrees

This dataset provides 0.125° gridded simulations of soil moisture and crop grain carbon for the Contiguous United States (CONUS), generated using the Community Land Model version 5 (CLM5) with the biogeochemistry module enabled. The data covers a historical baseline (1980–2015) and mid-century future projections (2020–2055). Future projections are organized into two sets of scenarios to distinguish the impacts of different drivers: (1) Atmospheric Only (ATM-only): These scenarios apply future atmospheric forcings while holding land use and land cover change (LULCC) at historical baseline levels. The atmospheric forcings represent moderately versus severely hotter/drier atmospheric conditions (dynamically downscaled perturbed thermodynamics simulations based on CMIP6 SSP245 and SSP585 warming signals), each with cooler versus hotter Earth System Model temperature sensitivity instantiations. These scenarios are identified in the file names as atm45cooler, atm45hotter, atm85cooler, and atm85hotter. (2) Coupled Atmospheric and Land-Use (LAND+ATM): These scenarios apply future atmospheric forcing together with future LULCC by pairing atmospheric pathways with lower versus higher population/economic growth scenarios representing Shared Socioeconomic Pathways 3 and 5 (SSP3 and SSP5). These scenarios are identified in the file names as atm45cooler_ssp3, atm45hotter_ssp3, atm45cooler_ssp5, atm45hotter_ssp5, atm85cooler_ssp3, atm85hotter_ssp3, atm85cooler_ssp5, and atm85hotter_ssp5. Please refer to the README file for detailed information on file structure, variables, units, and data formats.

Yao, Lili [Pacific Northwest National Laboratory] ↗
DOE OSTI2021

Demons in the North Atlantic: Variability of Deep Ocean Ventilation

Abstract Translation of atmospheric forcing variability into the ocean interior via ocean ventilation is an important aspect of transient climate change. On a seasonal timescale in the subtropics, this translation is mediated by a so‐called “Demon” that prevents access to all except late‐winter mixed‐layer water. Here, we use an eddy‐permitting numerical circulation model to investigate a similar process operating on longer (interannual) timescales in the subpolar North Atlantic. We find that variations in atmospheric forcing are mediated in their translation to the ocean interior, with year‐to‐year changes in the late‐winter mixed layer depth being the critical factor. The signature of persistent strong atmospheric forcing driving deep mixed layers is preferentially ventilated to the interior when the forcing is ceased. Susceptibility to this effect depends on the location and density of subduction—with the rate at which newly ventilated water escapes its region of subduction being the crucial factor.

MacGilchrist, G. A.↗
DOE OSTI2022

Rain-Induced Stratification of the Equatorial Indian Ocean and Its Potential Feedback to the Atmosphere

Surface freshening through precipitation can act to stably stratify the upper ocean, forming a rain layer (RL). RLs inhibit subsurface vertical mixing, isolating deeper ocean layers from the atmosphere. This process has been studied using observations and idealized simulations. The present ocean modeling study builds upon this body of work by incorporating spatially resolved and realistic atmospheric forcing. Fine-scale observations of the upper ocean collected during the Dynamics of the Madden-Julian Oscillation field campaign are used to verify the General Ocean Turbulence Model (GOTM). Spatiotemporal characteristics of equatorial Indian Ocean RLs are then investigated by forcing a 2D array of GOTM columns with realistic and well-resolved output from an existing regional atmospheric simulation. RL influence on the ocean-atmosphere system is evaluated through analysis of RL-induced modification to surface fluxes and sea surface temperature (SST). This analysis demonstrates that RLs cool the ocean surface on time scales longer than the associated precipitation event. A second simulation with identical atmospheric forcing to that in the first, but with rainfall set to zero, is performed to investigate the role of rain temperature and salinity stratification in maintaining cold SST anomalies within RLs. Approximately one third, or 0.1°C, of the SST reduction within RLs can be attributed to rain effects, while the remainder is attributed to changes in atmospheric temperature and humidity. The prolonged RL-induced SST anomalies enhance SST gradients that have been shown to favor the initiation of atmospheric convection. These findings encourage continued research of RL feedbacks to the atmosphere.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Quantifying sensitivity in numerical weather prediction-modeled offshore wind speeds through an ensemble modeling approach

A decade of research has shown that numerical weather prediction (NWP)-modeled wind speeds can be highly sensitive to the inputs and setups within the NWP model. For wind resource characterization applications, this sensitivity is often addressed by constructing a range of setups and selecting the one that best validates against observations. However, this approach is not possible in areas that lack high-quality hub height observations, especially offshore wind areas. In such cases, techniques to quantify and disseminate confidence in NWP-modeled wind speeds in the absence of observations are needed. We address this need in the present study and propose best practices for quantifying the spread in NWP-modeled wind speeds. We implement an ensemble approach in which we consider 24 different setups to the Weather Research and Forecasting (WRF) model. We construct the ensemble by considering variations in WRF version, WRF namelist, atmospheric forcing, and sea surface temperature (SST) forcing. Our analysis finds that the standard deviation produces more consistent estimates compared to the interquartile range and tends to be the more conservative estimator for ensemble variability. We further find that model spread increases closer to the surface and on shorter time scales. In conclusion, we explore methods to attribute total ensemble variability to the different ensemble components (e.g., atmospheric forcing and SST product) and find that contributions by components also vary depending on time scale. We anticipate that the methods and results presented in this paper will provide a reasonable foundation for further research into ensemble-based wind resource data sets.

17 WIND ENERGY↗
DOE OSTI2021

Disentangling the Effects of Vapor Pressure Deficit and Soil Water Availability on Canopy Conductance in a Seasonal Tropical Forest During the 2015 El Niño Drought

In land models used in Earth system models, vegetation water stress is heuristically related to soil water matric potential, and root water uptake in soil layers is simply proportional to root density to meet the transpiration demand. The efficiency of roots utilizing deeper soil water under water stress cannot be well represented, but it is of critical importance for transpiration regulation and carbon flux. To improve modeling of tropical forest response to drought, a plant hydrodynamic scheme (HYDRO), adapted from a trait-based model in the Functionally Assembled Terrestrial Ecosystem Simulator (FATES), has been implemented in the Energy Exascale Earth System Model (E3SM) Land Model (ELM) to represent the hydraulic dynamics and root water uptake strategy. The adapted model conforms to the big-leaf photosynthesis model in ELM. Vegetation water stress in the model is a function of prognostic leaf water potential and is used to regulate stomatal conductance. Point simulations were performed at two tropical forest sites in the Amazon and Panama. Model sensitivity to the numerical solver, plant hydraulic traits, soil properties, and atmospheric forcing has been evaluated. Our results showed that, even though water stress factors calculated at our study sites using HYDRO are comparable to that simulated by the heuristic water stress model, soil moisture prediction is improved with soil water redistributions through roots. Simulated results using HYDRO is not as sensitive to plant hydraulic traits measured for different tree species as to the soil properties and atmospheric forcing. Simplified numerical solver may cause inaccuracy as to where roots uptake water. This study suggests that plant hydrodynamic representation is critical for improved simulations of soil moisture under water stresses.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Data and Scripts Associated with "Modeling Ecohydrological Responses of Vegetation to Urban Microclimates Using the E3SM Land Model"

This dataset supports the study of vegetation ecohydrological responses to urban microclimates using the land component of the Energy Exascale Earth System Model (ELM) at four urban sites in Knoxville, Tennessee, USA. It includes the model inputs, simulation outputs, and associated scripts for running ELM simulations and analyzing the resulting data. The Model_Inputs folder includes static surface data, satellite-derived phenology (i.e., leaf area index), and atmospheric forcing data used to drive ELM simulations. Detailed descriptions of these datasets are provided in Section 2.3.2 of the associated manuscript. The Model_Outputs folder contains simulation results for the baseline, treatment, and ensemble experiments. Outputs from the baseline and treatment simulations are provided as raw ELM NetCDF files. Because the raw outputs from the 4,000-member ensemble are prohibitively large, the ensemble results are provided as summarized CSV files, which also serve as the source data for Figure 5 of the associated manuscript. The Scripts folder contains three components: E3SM, the core codebase of the Energy Exascale Earth System Model (E3SM); elm-olmt, the Offline Land Model Testbed (OLMT) used to perform the simulations; and knoxville_elm, which contains the analysis scripts used to process model outputs and generate the figures and results presented in the associated manuscript. Additional information is provided in Scripts_readme.txt within the Scripts directory.

Lu, Xiaoman [ORNL] (ORCID:0000000306698780)↗
DOE OSTI2026

High-Resolution WRF-Based Downscaling of Earth System Model Projections for Energy Applications across CONUS

Evaluating energy resources under future scenarios requires meteorological information that adequately resolves regional-scale variability and is suitable for regional energy system studies. Although Earth system model (ESM) outputs provide essential large-scale context, their coarse resolution and inherent biases limit direct use in energy system applications. This work presents a high-resolution dynamical downscaling framework using the Weather Research and Forecasting (WRF) model to generate energy-relevant regional fields for future scenarios over the contiguous United States (CONUS). The framework first identifies an optimal WRF configuration through sensitivity experiments, then evaluates raw and bias-corrected ESM initial and boundary conditions, with soil moisture and soil temperature bias correction implemented as an integral component of the bias-corrected ESM atmospheric forcing prior to WRF dynamical downscaling to improve land-atmosphere interactions. Simulations performed at 4-km resolution show that uncorrected ESM forcing leads to systematically dry and cold soil states, which propagate into elevated near-surface air temperature and solar irradiance biases, particularly during summer for the period 2000-2014. Incorporating bias-corrected atmospheric forcing together with soil state bias correction substantially reduces these errors and improves the representation of surface energy processes in WRF simulations. The results highlight the importance of bias-aware initialization strategies in high-resolution dynamical downscaling for future energy system analysis and planning.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Northern Hemisphere Snow Drought in Earth System Model Simulations and ERA5‐Land Data in 1980–2014

Abstract Low snow levels over the past few decades and predictions of a low‐to‐no snow future have spurred research into snow droughts, which pose a threat to water security and management. Systematic data‐model comparisons of snow drought have been lacking, hindering our understanding of the drivers of snow drought in the past. To address this gap, we analyzed snow drought events using standardized snow water equivalent index derived from monthly results of four numerical experiments using the E3SM Land Model (ELM) and ERA5‐Land data during the period of 1980–2014. Additionally, we compared snow drought duration calculated from models with those from the ERA5‐Land data during selected El Niño‐Southern Oscillation (ENSO) years. The numerical experiments were conducted with ELM driven by two prescribed atmospheric forcings, and with the coupled land‐atmosphere configuration of E3SM with and without plant hydraulics scheme feedback. Analysis reveals that 20%–30% of snow droughts occur due to factors other than above‐normal temperature and low snowfall, such as low soil moisture, warm soil temperature, and low relative humidity, etc., especially in high latitudes (50° North). Furthermore, our study highlights the exacerbating effect of ENSO events on snow drought conditions in various regions, despite some discrepancies between model and ERA5‐Land results. We also identified limitations of the coupled land‐atmosphere models in our current configuration in capturing the spatial patterns of snow droughts. This study underscores the challenge of predicting and mitigating snow drought and the need for a comprehensive understanding of the factors contributing to snow drought.

54 ENVIRONMENTAL SCIENCES↗