Search NASA⌕ Search

SEARCH · Search NASA

Results for “high latitude”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Arctic River Delta Morphologic Variability and Implications for Riverine Fluxes to the Coast

Arctic riverine fluxes are anticipated to increase as the Arctic warms and have a large impact on the Arctic ocean. Deltas modify the spatial and temporal distributions of riverine fluxes, but no thorough studies have been conducted to analyze Arctic delta morphologies to determine their influence on land–ocean fluxes. We performed an analysis of six high–latitude deltas (Colville, Kolyma, Lena, Mackenzie, Yenisei, and Yukon) to characterize delta morphologies and determine the influence of morphology on the distribution of fluxes to the coast. All six deltas deliver material to the coast at discrete locations across small areas despite differences in delta shoreline length. Large Arctic deltas exhibit large variability in channel width, which we hypothesize is due to a feedback with ice cover and retreat that favors the growth of large channels over geologic timescales. Spatial variability in island sizes suggests variability in channel activity, island nourishment, and susceptibility to drowning by sea level rise. Potential lake storage is highest on the Mackenzie delta, thus providing a means for reducing nutrient and sediment loading of the coastal ocean. Here, the connected lakes are also prevalent on the Colville and Yukon deltas, suggesting that these deltas can filter riverine fluxes even when the deltas are not flooded. Differences in Arctic delta morphologies can be explained by varying levels of riverine and marine influence, antecedent topography, and local channel dynamics. Ice cover also plays a large role in controlling Arctic delta morphologies and dynamics that has not been previously represented in interpretations of existing delta metrics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Stable Methane Isotopologues From Northern Lakes Suggest That Ebullition Is Dominated by Sub‐Lake Scale Processes

Abstract Stable isotopes have emerged as popular study targets when investigating emission of methane (CH 4 ) from lakes. Yet little is known on how isotopic patterns conform to variations in emission magnitudes—a highly relevant question. Here, we present a large multiyear data set on stable isotopes of CH 4 ebullition (bubbling) from three small adjacent subarctic lakes. The δ 13 C‐CH 4 and δD‐CH 4 range from −78.4‰ to −53.1‰ and from −369.8‰ to −218.8‰, respectively, and vary greatly among the lakes. The signatures suggest dominant hydrogenotrophic methanogenesis, particularly in the deep zones, but there are also signals of seemingly acetoclastic production in some high fluxing shallow areas, possibly fueled by in situ vegetation, but in‐sediment anaerobic CH 4 oxidation cannot be ruled out as an alternative cause. The observed patterns, however, are not consistent across the lakes. Neither do they correspond to the spatiotemporal variations in the measured bubble CH 4 fluxes. Patterns of acetoclastic and hydrogenotrophic production plus oxidation demonstrate that gains and losses of sediment CH 4 are dominated by sub‐lake scale processes. The δD‐CH 4 in the bubbles was significantly different depending on measurement month, likely due to evaporation effects. On a larger scale, our isotopic data, combined with those from other lakes, show a significant difference in bubble δD‐CH 4 between postglacial and thermokarst lakes, an important result for emission inventories. Although this characteristic theoretically assists in source partitioning studies, most hypothetical future shifts in δD‐CH 4 due to high‐latitude lake area or production pathway are too small to lead to atmospheric changes detectable with current technology.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Validation and Sensitivity Analysis of a 1-D Lake Model Across Global Lakes

Lakes have important influence on weather and climate from local to global scales. However, their prediction using numerical models is notoriously difficult because lakes are highly heterogeneous across the globe, but observations are sparse. In this study, we assessed the performance of a 1-D lake model in simulating the thermal structures of 58 lakes with diverse morphometric and geographic characteristics by following the phase 2a local lake protocol of the Inter-sectoral Impact Model Intercomparison Project (ISIMIP2a). After calibration, the root-mean-square errors (RMSE) were below 2 °C for 70% and 75% of the lakes for epilimnion and full-profile temperature simulations, with an average of 1.71 °C and 1.43 °C, respectively. The model performance mainly depended on lake shape rather than location, supporting the possibility of grouping model parameters by lake shape for global applications. Furthermore, through machine-learning based parameter sensitivity tests, we identified turbulent heat fluxes, wind-driven mixing and water transparency as the major processes controlling lake thermal and mixing regimes. Snow density was also important for modeling the ice phenology of high-latitude lakes. The relative influence of the key processes and the corresponding parameters mainly depended on lake latitude and depth. Turbulent heat fluxes showed a decreasing importance in affecting epilimnion temperature with increasing latitude. Wind-driven mixing was less influential to lake stratification for deeper lakes while the impact of light extinction, on the contrary, showed a positive correlation with depth. Our findings may guide improvements in 1-D lake model parameterizations to achieve higher fidelity in simulating global lake thermal dynamics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

CMIP6 Models Underestimate the Holton‐Tan Effect

Abstract The teleconnection between the Quasi‐Biennial Oscillation (QBO) and the Arctic polar vortex is investigated using Coupled Model Intercomparison Project 6 (CMIP6) models. Output from 14 CMIP6 models is compared with reanalysis, three experiments with prescribed QBOs, one of which has no free polar stratospheric variability, and transient experiments in which a QBO is prescribed in runs previously devoid of a QBO. Each CMIP6 model underestimates the Holton‐Tan effect (HTE), the weakening of the polar vortex expected with QBO easterlies in the tropical lower stratosphere. To establish why, potential vorticity maps are used to investigate longitudinal variations in the teleconnection. Prescribing easterly QBO in the transient experiments promotes more high‐latitude planetary wave breaking by influencing the mid‐latitude stratospheric circulation, particularly over Asia. CMIP6 models that better simulate this response over Asia better simulate the HTE. These models also have stronger 10 hPa QBO westerlies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

The Global Fingerprint of Modern Ice-Mass Loss on 3-D Crustal Motion

Crustal motion generated by rapid ice-mass loss from Earth's glaciers and ice sheets has previously been considered in Global Navigational Satellite System (GNSS) analyses and numerical models across regions of ice retreat. However, the fingerprint of ice-mass loss is not limited to glaciated areas, but is characterized by a global pattern of 3-D crustal deformation. We compute “far-field” vertical and horizontal deformation rates that occurred in response to early 21st century mass flux from the Greenland and Antarctic Ice Sheets, global glaciers and ice caps, and associated ocean loading. We demonstrate that mass changes in the Greenland Ice Sheet and high latitude glacier systems each generated average crustal motion of 0.1–0.4 mm/yr across much of the Northern Hemisphere, with significant year-to-year variability in magnitude and direction. Horizontal motions associated with ice-mass loss exceed vertical rates in many far-field areas, and both should be considered in future analysis of GNSS measurements.

58 GEOSCIENCES↗
DOE OSTI2022

An Assessment of Nonhydrostatic and Hydrostatic Dynamical Cores at Seasonal Time Scales in the Energy Exascale Earth System Model (E3SM)

Abstract In global atmospheric modeling, the differences between nonhydrostatic (NH) and hydrostatic (H) dynamical cores are negligible in dry simulations when grid spacing is larger than 10 km. However, recent studies suggest that those differences can be significant at far coarser resolution when moisture is included. To better understand how NH and H differences manifest in global fields, we perform and analyze an ensemble of 28 and 13 km seasonal simulations with the NH and H dynamical cores in the Energy Exascale Earth System Model global atmosphere model, where the differences between H and NH configurations are minimized. A set of idealized rising bubble experiments is also conducted to further investigate the differences. Although NH and H differences are not significant in global statistics and zonal averages, significant differences in precipitation amount and patterns are observed in parts of the tropics. The most prominent differences emerge near India and the Western Pacific in the boreal summer, and the central‐southern Indian Ocean and Pacific in the boreal winter. Tropical differences influence surrounding regions through modification of the regional circulation and can propagate to the extratropics, leading to significant temperature and geopotential differences over the middle to high latitudes. While the dry bubble experiments show negligible deviation between H and NH dynamics until grid spacing is below 6.25 km, precipitation amount and vertical velocity are different in the moist case even at 25 km resolution.

58 GEOSCIENCES↗
DOE OSTI2022

Are Terrestrial Biosphere Models Fit for Simulating the Global Land Carbon Sink?

The Global Carbon Project estimates that the terrestrial biosphere has absorbed about one-third of anthropogenic CO 2 emissions during the 1959–2019 period. This sink-estimate is produced by an ensemble of terrestrial biosphere models and is consistent with the land uptake inferred from the residual of emissions and ocean uptake. The purpose of our study is to understand how well terrestrial biosphere models reproduce the processes that drive the terrestrial carbon sink. One challenge is to decide what level of agreement between model output and observation-based reference data is adequate considering that reference data are prone to uncertainties. To define such a level of agreement, we compute benchmark scores that quantify the similarity between independently derived reference data sets using multiple statistical metrics. Models are considered to perform well if their model scores reach benchmark scores. Our results show that reference data can differ considerably, causing benchmark scores to be low. Model scores are often of similar magnitude as benchmark scores, implying that model performance is reasonable given how different reference data are. While model performance is encouraging, ample potential for improvements remains, including a reduction in a positive leaf area index bias, improved representations of processes that govern soil organic carbon in high latitudes, and an assessment of causes that drive the inter-model spread of gross primary productivity in boreal regions and humid tropics. The success of future model development will increasingly depend on our capacity to reduce and account for observational uncertainties.

58 GEOSCIENCES↗
DOE OSTI2021

Retreat and Regrowth of the Greenland Ice Sheet During the Last Interglacial as Simulated by the CESM2‐CISM2 Coupled Climate–Ice Sheet Model

Abstract During the Last Interglacial, approximately 129 to 116 ka (thousand years ago), the Arctic summer climate was warmer than the present, and the Greenland Ice Sheet retreated to a smaller extent than its current state. Previous model‐derived and geological reconstruction estimates of the sea‐level contribution of the Greenland Ice Sheet during the Last Interglacial vary widely. Here, we conduct a transient climate simulation from 127 to 119 ka using the Community Earth System Model (CESM2), which includes a dynamic ice sheet component (the Community Ice Sheet Model, CISM2) that is interactively coupled to the atmosphere, land, ocean, and sea ice components. Vegetation distribution is updated every 500 years based on biomes simulated using a monthly climatology to force the BIOME4 equilibrium vegetation model. Results show a substantial retreat of the Greenland Ice Sheet, reaching a minimum extent at 121.9 ka, equivalent to a 3.0 m rise in sea level relative to the present day, followed by gradual regrowth. In contrast, a companion simulation employing static vegetation based on pre‐industrial conditions shows a much smaller ice‐sheet retreat, highlighting the importance of the changes in high‐latitude vegetation distribution for amplifying the ice‐sheet response.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

The Intrinsic 150‐Day Periodicity of the Southern Hemisphere Extratropical Large‐Scale Atmospheric Circulation

Abstract The variability of the Southern Hemisphere (SH) extratropical large‐scale circulation is dominated by the Southern Annular Mode (SAM), whose timescale is extensively used as a key metric in evaluating state‐of‐the‐art climate models. Past observational and theoretical studies suggest that the SAM lacks any internally generated (intrinsic) periodicity. Here, we show, using observations and a climate model hierarchy, that the SAM has an intrinsic 150‐day periodicity. This periodicity is robustly detectable in the power spectra and principal oscillation patterns (aka dynamical mode decomposition) of the zonal‐mean circulation, and in hemispheric‐scale precipitation and ocean surface wind stress. The 150‐day period is consistent with the predictions of a new reduced‐order model for the SAM, which suggests that this periodicity is associated with a complex interaction of turbulent eddies and zonal wind anomalies, as the latter propagate from low to high latitudes. These findings present a rare example of periodic oscillations arising from the internal dynamics of the extratropical turbulent circulations. Based on these findings, we further propose a new metric for evaluating climate models, and show that some of the previously reported shortcomings and improvements in simulating SAM's variability connect to the models' ability in reproducing this periodicity. We argue that this periodicity should be considered in evaluating climate models and understanding the past, current, and projected Southern Hemisphere climate variability.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Adaptation Strategies Strongly Reduce the Future Impacts of Climate Change on Simulated Crop Yields

Abstract Simulations of crop yield due to climate change vary widely between models, locations, species, management strategies, and Representative Concentration Pathways (RCPs). To understand how climate and adaptation affects yield change, we developed a meta‐model based on 8703 site‐level process‐model simulations of yield with different future adaptation strategies and climate scenarios for maize, rice, wheat and soybean. We tested 10 statistical models, including some machine learning models, to predict the percentage change in projected future yield relative to the baseline period (2000–2010) as a function of explanatory variables related to adaptation strategy and climate change. We used the best model to produce global maps of yield change for the RCP4.5 scenario and identify the most influential variables affecting yield change using Shapley additive explanations. For most locations, adaptation was the most influential factor determining the projected yield change for maize, rice and wheat. Without adaptation under RCP4.5, all crops are expected to experience average global yield losses of 6%–21%. Adaptation alleviates this average projected loss by 1–13 percentage points. Maize was most responsive to adaptive practices with a projected mean yield loss of −21% [range across locations: −63%, +3.7%] without adaptation and −7.5% [range: −46%, +13%] with adaptation. For maize and rice, irrigation method and cultivar choice were the adaptation types predicted to most prevent large yield losses, respectively. When adaptation practices are applied, some areas are predicted to experience yield gains, especially at northern high latitudes. These results reveal the critical importance of implementing adequate adaptation strategies to mitigate the impact of climate change on crop yields.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Atmospheric Radiative and Oceanic Biological Productivity Responses to Increasing Anthropogenic Combustion‐Iron Emission in the 1850–2010 Period

Abstract Anthropogenic emission is an important component of the present‐day iron cycle yet its long‐term impacts on climate are poorly understood. Iron mineralogy strongly affects its radiative and oceanic interactions and was unrepresented in previous studies. We perform simulations using a mineralogy‐based inventory and an atmospheric transport model and estimate the 1850–2010 global mean direct radiative forcing (DRF) to be +0.02 to +0.10 W/m 2 . We estimate that the CO 2 sequestration of 0.2–13 ppmv over the last 150 years due to enhanced phytoplankton productivity by anthropogenic iron deposition causes an avoided CO 2 forcing of −0.002 to −0.16 W/m 2 . While globally small, these impacts can be higher in specific regions; the anthropogenic DRF is +0.5 W/m 2 over areas with more coal combustion and metal smelting, and anthropogenic soluble iron sustains >10% of marine net primary productivity in the high‐latitude North Pacific Ocean, a region vulnerable to stratification due to climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Orographic Cirrus and Its Radiative Forcing in NCAR CAM6

Cirrus clouds play an important role in the Earth's radiative energy budget, thereby affecting the climate state and climate change. Orographic gravity wave (OGW)-induced sub-grid scale vertical velocity (i.e., cooling rate) is not resolved by large-scale models and its impact on ice formation in cirrus clouds is not well quantified. In this study, one sub-grid scale OGW scheme (e.g., McFarlane) is used in the Community Atmosphere Model version 6 (CAM6) to generate vertical velocity variance (σw) for cirrus formation. Results from the default model and simulations with the OGW-induced σw are evaluated against the DOE ARM Small Particles in Cirrus (SPARTICUS) campaign observations. The OGW based on the McFarlane scheme increases the sub-grid scale σw over mountains compared to the default model and improves the model agreement with the SPARTICUS observations. Larger σ w due to OGWs can trigger more frequent homogeneous nucleation in orographic cirrus and generates a higher number concentration of ice crystals observed during the SPARTICUS campaign. Moreover, our evaluation of the model simulations against satellite observations indicates that the McFarlane scheme generates high in-cloud ice number concentrations (>200 L -1 ) in the upper troposphere over mountains and high plateaus at mid- and high-latitudes of the winter hemisphere as shown in the observations. More ice crystals with smaller sizes absorb more infrared radiation (+0.523 ± 0.125 W m -2 ). The net radiative cloud forcing change at the top of the atmosphere is +0.330 W m -2 due to the orographic cirrus.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Near‐Surface Hydrology and Soil Properties Drive Heterogeneity in Permafrost Distribution, Vegetation Dynamics, and Carbon Cycling in a Sub‐Arctic Watershed

Abstract Discontinuous permafrost environments exhibit strong spatial heterogeneity at scales too small to be driven by weather forcing or captured by Earth System Models. Here we analyze effects of observed spatial heterogeneity in soil and vegetation properties, hydrology, and thermal dynamics on ecosystem carbon dynamics in a watershed on the Seward Peninsula in Alaska. We apply a Morris global sensitivity analysis to a process‐rich, successfully tested terrestrial ecosystem model (TEM), ecosys , varying soil properties, boundary conditions, and weather forcing. We show that landscape heterogeneity strongly impacts soil temperatures and vegetation composition. Snow depth, O‐horizon thickness, and near‐surface water content, which vary at scales of O(m), control the soil thermal regime more than an air temperature gradient corresponding to a 140 km north–south distance. High shrub productivity is simulated only in talik (perennially unfrozen) soils with high nitrogen availability. Through these effects on plant and permafrost dynamics, landscape heterogeneity impacts ecosystem productivity. Simulations with near‐surface taliks have higher microbial respiration (by 78.0 gC m −2 yr −1 ) and higher net primary productivity (by 104.9 gC m −2 yr −1 ) compared to runs with near‐surface permafrost, and simulations with high shrub productivity have outlying values of net carbon uptake. We explored the prediction uncertainty associated with ignoring observed landscape heterogeneity, and found that watershed net carbon uptake is 60% larger when heterogeneity is accounted for. Our results highlight the complexity inherent in discontinuous permafrost environments and demonstrate that missing representation of subgrid heterogeneity in TEMs could bias predictions of high‐latitude carbon budget.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

An Ensemble of Neural Networks for Moist Physics Processes, Its Generalizability and Stable Integration

Abstract With the recent advances in data science, machine learning has been increasingly applied to convection and cloud parameterizations in global climate models (GCMs). This study extends the work of Han et al. (2020, https://doi.org/10.1029/2020MS002076 ) and uses an ensemble of 32‐layer deep convolutional residual neural networks, referred to as ResCu‐en, to emulate convection and cloud processes simulated by a superparameterized GCM, SPCAM. ResCu‐en predicts GCM grid‐scale temperature and moisture tendencies, and cloud liquid and ice water contents from moist physics processes. The surface rainfall is derived from the column‐integrated moisture tendency. The prediction uncertainty inherent in deep learning algorithms in emulating the moist physics is reduced by ensemble averaging. Results in 1‐year independent offline validation show that ResCu‐en has high prediction accuracy for all output variables, both in the current climate and in a warmer climate with +4K sea surface temperature. The analysis of different neural net configurations shows that the success to generalize in a warmer climate is attributed to convective memory and the 1‐dimensional convolution layers incorporated into ResCu‐en. We further implement a member of ResCu‐en into CAM5 with real world geography and run the neural‐network‐enabled CAM5 (NCAM) for 5 years without encountering any numerical integration instability. The simulation generally captures the global distribution of the mean precipitation, with a better simulation of precipitation intensity and diurnal cycle. However, there are large biases in temperature and moisture in high latitudes. These results highlight the importance of convective memory and demonstrate the potential for machine learning to enhance climate modeling.

Meteorology & Atmospheric Sciences↗
DOE OSTI2024

Ground Heat Flux Reconstruction Using Bayesian Uncertainty Quantification Machinery and Surrogate Modeling

Ground heat flux (G 0 ) is a key component of the land-surface energy balance of high-latitude regions. Despite its crucial role in controlling permafrost degradation due to global warming, G 0 is sparsely measured and not well represented in the outputs of global scale model simulation. In this study, an analytical heat transfer model is tested to reconstruct G 0 across seasons using soil temperature series from field measurements, Global Climate Model, and climate reanalysis outputs. The probability density functions of ground heat flux and of model parameters are inferred using available G 0 data (measured or modeled) for snow-free period as a reference. When observed G 0 is not available, a numerical model is applied using estimates of surface heat flux (dependent on parameters) as the top boundary condition. These estimates (and thus the corresponding parameters) are verified by comparing the distributions of simulated and measured soil temperature at several depths. Aided by state-of-the-art uncertainty quantification methods, the developed G 0 reconstruction approach provides novel means for assessing the probabilistic structure of the ground heat flux for regional permafrost change studies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Different Strategies of Stratospheric Aerosol Injection Would Significantly Affect Climate Extreme Mitigation

Abstract Stratospheric aerosol injection (SAI) has been proposed as a potential supplement to mitigate some climate impacts of anthropogenic warming. Using Community Earth System Model ensemble simulation results, we analyze the response of temperature and precipitation extremes to two different SAI strategies: one injects SO 2 at the equator to stabilize global mean temperature and the other injects SO 2 at multiple locations to stabilize global mean temperature as well as the interhemispheric and equator‐to‐pole temperature gradients. Our analysis shows that in the late 21st century, compared with the present‐day climate, both equatorial and multi‐location injection lead to reduced hot extremes in the tropics, corresponding to overcooling of the mean climate state. In mid‐to‐high latitude regions, in comparison to the present‐day climate, substantial decreases in cold extremes are observed under both equatorial and multi‐location injection, corresponding to residual winter warming of the mean climate state. Both equatorial and multi‐location injection reduce precipitation extremes in the tropics below the present‐day level, associated with the decrease in mean precipitation. Overall, for most regions, temperature and precipitation extremes show reduced change in response to multi‐location injection than to equatorial injection, corresponding to reduced mean climate change for multi‐location injection. In comparison with equatorial injection, in response to multi‐location injection, most land regions experience fewer years with significant change in cold extremes from the present‐day level, and most tropical regions experience fewer years with significant change in hot extremes. The design of SAI strategies to mitigate anthropogenic climate extremes merits further study.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Larger Cloud Liquid Water Enhances Both Aerosol Indirect Forcing and Cloud Radiative Feedback in Two Earth System Models

Previous studies have noticed that the Coupled Model Intercomparison Project Phase 6 (CMIP6) models with a stronger cooling from aerosol-cloud interactions (ACI) also have an enhanced warming from positive cloud feedback, and these two opposing effects are counter-balanced in simulations of the historical period. However, reasons for this anti-correlation are less explored. In this study, we perturb the cloud ice microphysical processes to obtain cloud liquid of varying amounts in two Earth System Models (ESMs). We find that the model simulations with a larger liquid water path (LWP) tend to have a stronger cooling from ACI and a stronger positive cloud feedback. More liquid clouds in the mean-state present more opportunities for anthropogenic aerosol perturbations and also weaken the negative cloud feedback at middle to high latitudes. This work, from a cloud state perspective, emphasizes the influence of the mean-state LWP on effective radiative forcing due to ACI (ERF ACI ).

54 ENVIRONMENTAL SCIENCES↗
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↗