Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “regional climate simulation”

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 OSTI2022

Numerical Errors in Ice Microphysics Parameterizations and their Effects on Simulated Regional Climate

The major characteristics of ice microphysics in Weather Research and Forecasting (WRF) Double-Moment 6-class (WDM6) bulk-type cloud microphysics originate from the diagnosed ice number concentration, which is a function of the cloud-ice mixing ratio. In this study, we correct numerical errors in ice microphysics processes of the WDM6, in which the cloud-ice shape is assumed as single bullets and examine the impact on regional climate simulations. Here, by rederiving the relationships between cloud microphysics characteristics, including the one linking the cloud-ice mixing ratio and number concentration, we remove numerical errors intrinsic to the description of cloud-ice characteristics in the original WDM6 microphysics scheme. The revised WDM6 is tested using a WRF framework for regional climate simulations over the East Asian region. We find that our correction to the WDM6 improves the model’s performance in capturing the observed distribution of the monsoon rain band. A reduction in cloud ice is significant in the revised WDM6, which strengthens the Western North Pacific High. By conducting the additional sensitivity experiment in which the characteristics of cloud-ice shape are revised as the one for the column type, our study also finds out that the impacts of the existing numerical errors on the simulated monsoon is as large as the ones of the changes in cloud-ice shape.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

How well are hazards associated with derechos reproduced in regional climate simulations?

Abstract. A 15-member ensemble of convection-permitting regional simulations of the fast-moving and destructive derecho of 29–30 June 2012 that impacted the northeastern urban corridor of the USA is presented. This event generated 1100 reports of damaging winds, generated significant wind gusts over an extensive area of up to 500 000 km2, caused several fatalities, and resulted in widespread loss of electrical power. Extreme events such as this are increasingly being used within pseudo-global-warming experiments to examine the sensitivity of historical, societally important events to global climate non-stationarity and how they may evolve as a result of changing thermodynamic and dynamic contexts. As such it is important to examine the fidelity with which such events are described in hindcast experiments. The regional simulations presented herein are performed using the Weather Research and Forecasting (WRF) model. The resulting ensemble is used to explore simulation fidelity relative to observations for wind gust magnitudes, spatial scales of convection (as is manifest in high composite reflectivity, cREF), and both rainfall and hail production as a function of model configuration (microphysics parameterization, lateral boundary conditions (LBCs), start date, use of nudging, compiler choice, damping, and number of vertical levels). We also examine the degree to which each ensemble member differs with respect to key mesoscale drivers of convective systems (e.g., convective available potential energy and vertical wind shear) and critical manifestations of deep convection, e.g., vertical velocities, cold-pool generation, and how those properties relate to the correct characterization of the associated atmospheric hazards (wind gusts and hail). Use of a double-moment, seven-class scheme with number concentrations for all species (including hail and graupel) results in the greatest fidelity of model-simulated wind gusts and convective structure to the observations of this event. All ensemble members, however, fail to capture the intensity of the event in terms of the spatial extent of convection and the production of high near-surface wind gusts. We further show very high sensitivity to the LBCs employed and specifically that simulation fidelity is higher for simulations nested within ERA-Interim compared to ERA5. Excess convective available potential energy (CAPE) in all ensemble members after the derecho passage leads to excess production of convective cells, wind gusts, cREF > 40 dBZ, and precipitation during a frontal passage on the subsequent day. This event proved very challenging to forecast in real time and to reproduce in the 15-member hindcast simulation ensemble presented here. Future work could examine if simulations with other initial and lateral boundary conditions can achieve greater fidelity.

Shepherd, Tristan (ORCID:0000000186276419)↗
DOE OSTI2026

Evaluation of a high-resolution regional climate simulation for surface and hub-height wind climatology over North America

Assessing the availability of key wind resources requires augmenting observations to support the implementation of wind energy infrastructure. However, observations are limited, necessitating the development of high-resolution, long-term gridded datasets. This study presents a robust, dynamically downscaled climatological dataset, offering 20 years of hourly wind data at a 4 km spatial resolution across North America, and evaluates its performance against observations, including meteorological towers and automated surface-observing system (ASOS) stations, as well as coarse-resolution reanalysis data (the European Centre for Medium-Range Weather Forecasts (ECMWF) reanalysis version 5 (ERA5)). Results demonstrate that the downscaled high-resolution wind data outperform ERA5 in regions of complex terrain and coastal areas, with improved overlap coefficients for wind data distributions and reduced root mean square errors (RMSEs) for hub-height and near-surface diurnal wind patterns. The downscaled simulation also captures the synoptic drivers of seasonal wind direction patterns reasonably well, indicated by high wind rose similarity indices. This study also provides an analysis of interannual variability, utilizing the dataset's full 20-year period, and model uncertainty, generated by varying model initial conditions and physics parameterizations across 1-year ensemble members, which are key considerations for wind resource assessment in wind farm development.

17 WIND ENERGY↗
DOE OSTI2022

Evaluating seasonal and regional distribution of snowfall in regional climate model simulations in the Arctic

In this study, we investigate how the regional climate model HIRHAM5 reproduces the spatial and temporal distribution of Arctic snowfall when compared to CloudSat satellite observations during the examined period of 2007–2010. For this purpose, both approaches, i.e., the assessments of the surface snowfall rate (observation-to-model) and the radar reflectivity factor profiles (model-to-observation), are carried out considering spatial and temporal sampling differences. The HIRHAM5 model, which is constrained in its synoptic representation by nudging to ERA-Interim, represents the snowfall in the Arctic region well in comparison to CloudSat products. The spatial distribution of the snowfall patterns is similar in both identifying the southeastern coast of Greenland and the North Atlantic corridor as regions gaining more than twice as much snowfall as the Arctic average, defined here for latitudes between 66 and 81°N. Excellent agreement (difference less than 1%) in the Arctic-averaged annual snowfall rate between HIRHAM5 and CloudSat is found, whereas ERA-Interim reanalysis shows an underestimation of 45% and significant deficits in the representation of the snowfall rate distribution. From the spatial analysis, it can be seen that the largest differences in the mean annual snowfall rates are an overestimation near the coastlines of Greenland and other regions with large orographic variations as well as an underestimation in the northern North Atlantic Ocean. To a large extent, the differences can be explained by clutter contamination, blind zone or higher resolution of CloudSat measurements, but clearly HIRHAM5 overestimates the orographic-driven precipitation. The underestimation of HIRHAM5 within the North Atlantic corridor south of Svalbard is likely connected to a poor description of the marine cold air outbreaks which could be identified by separating snowfall into different circulation weather type regimes. By simulating the radar reflectivity factor profiles from HIRHAM5 utilizing the Passive and Active Microwave TRAnsfer (PAMTRA) forward-modeling operator, the contribution of individual hydrometeor types can be assessed. Looking at a latitude band at 72–73°N, snow can be identified as the hydrometeor type dominating radar reflectivity factor values across all seasons. The largest differences between the observed and simulated reflectivity factor values are related to the contribution of cloud ice particles, which is underestimated in the model, most likely due to the small sizes of the particles. The model-to-observation approach offers a promising diagnostic when improving cloud schemes, as illustrated by comparison of different schemes available for HIRHAM5.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Assessment of CMIP5 and CMIP6 AMIP Simulated Clouds and Surface Shortwave Radiation Using ARM Observations over Different Climate Regions

The simulations of clouds and surface radiation from 10 CMIP6 models and their CMIP5 predecessors are compared to the ARM ground-based observations over different climate regions. Compared to the ARM radar-lidar derived total cloud fractions (CFT) and cloud fraction vertical distributions over the six selected sites, both CMIP5 and CMIP6 significantly underestimated CFT and low-level CF over the Northern Hemispheric midlatitude sites (SGPC1 and ENAC1), although the biases are generally smaller in CMIP6. Over the tropical oceanic site (TWPC2), 5 out of 10 CMIP6 models better simulated low-level CF than their CMIP5 predecessors. CMIP6 simulations generally agreed well with the ARM observations in CFT and cloud fraction vertical distributions over the tropical continental (MAOM1) and coastal (TWPC3) sites but missed the transitions between dry and wet seasons, similar to CMIP5 simulations. The improvements in downwelling shortwave fluxes (SWdn) at the surface from the majority of CMIP6 compared to CMIP5 primarily resulted from the improved cloud fraction simulations, especially over the SGPC1, ENAC1, and TWPC3 sites. By contrast, both CMIP5 and CMIP6 models exhibited diverse performances of clouds and shortwave radiation over the Arctic site (NSAC1), where CMIP6 models produced more clouds than CMIP5 models, especially for the low-level clouds. The comparisons between observations and CMIP5 and CMIP6 simulations provide valuable quantitative assessments of the accuracy of mean states and variabilities in the model simulations and shed light on general directions to improve climate models in different regions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Leveraging regional mesh refinement to simulate future climate projections for California using the Simplified Convection-Permitting E3SM Atmosphere Model Version 0

Abstract. The spatial heterogeneity related to complex topography in California demands high-resolution (< 5 km) modeling, but global convection-permitting climate models are computationally too expensive to run multi-decadal simulations. We developed a 3.25 km California climate modeling framework by leveraging regional mesh refinement (CARRM) using the U.S. Department of Energy (DOE)'s global Simple Cloud-Resolving E3SM Atmosphere Model (SCREAM) version 0. Four 5-year time periods (2015–2020, 2029–2034, 2044–2049, and 2094–2099) were simulated by nudging CARRM outside California to 1° coupled simulation of E3SMv1 under the Shared Socioeconomic Pathways (SSP)5-8.5 future scenario. The 3.25 km grid spacing adds considerable value to the prediction of the California climate changes, including more realistic high temperatures in the Central Valley and much improved spatial distributions of precipitation and snowpack in the Sierra Nevada and coastal stratocumulus. Under the SSP5-8.5 scenario, CARRM simulation predicts widespread warming of 6–10 °C over most of California, a 38 % increase in statewide average 30 d winter–spring precipitation, a near-complete loss of the alpine snowpack, and a sharp reduction in shortwave cloud radiative forcing associated with marine stratocumulus by the end of the 21st century. We note a climatological wet precipitation bias for the CARRM and discuss possible reasons. We conclude that SCREAM RRM is a technically feasible and scientifically valid tool for climate simulations in regions of interest, providing an excellent bridge to global convection-permitting simulations.

58 GEOSCIENCES↗
DOE OSTI2024

Evaluation of precipitation across the contiguous United States, Alaska, and Puerto Rico in multi-decadal convection-permitting simulations

Abstract This study is an early effort to generate a multi-decadal convection-permitting regional climate dataset that covers nearly the entire North American continent. We assessed a 20 year dynamically downscaled regional climate simulation at a 4 km spatial resolution with explicit convection across the contiguous United States (CONUS), Alaska, and Puerto Rico. Specifically, we evaluated the model’s performance in representing mean, 95th percentile, and extreme precipitation across regions. Our findings indicate that when compared with ERA5 reanalysis, the forcing data, convection-permitting simulation improves representations of seasonal, 95th percentile, and extreme precipitation over a large portion of the CONUS, Alaska, and Puerto Rico, particularly in areas where precipitation is heaviest. The simulation adds value over its forcing data (ERA5) in up to 53% of all grid cells in the CONUS, 68.8% in Alaska, and 84.0% in Puerto Rico. It is important to note that, however, despite improvements, model errors in Puerto Rico remain large. Similar improvements are observed in extreme indices, including consecutive dry days, maximum 5 days precipitation, and extreme precipitation. Analysis of the diurnal cycle of mean hourly precipitation suggests that representations of convective processes—including onset, dissipation, suppression, downstream propagation, and local circulation—improved overall.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Characteristics of Precipitation and Mesoscale Convective Systems Over the Peruvian Central Andes in Multi 5‐Year Convection‐Permitting Simulations

Abstract Using the Weather Research and Forecasting model with two planetary boundary layer schemes, ACM2 and MYNN, convection‐permitting model (CPM) regional climate simulations were conducted for a 6‐year period, including a one‐year spin‐up period, at a 15‐km grid spacing covering entire South America and a nested convection‐permitting 3‐km grid spacing covering the Peruvian central Andes region. These two CPM simulations along with a 4‐km simulation covering South America produced by National Center for Atmospheric Research (NCAR), three gridded precipitation products, and rain gauge data in Peru and Brazil, are used to document the characteristics of precipitation and mesoscale convective systems (MCSs) in the Peruvian central Andes region. Results show that all km‐scale simulations generally capture the spatiotemporal patterns of precipitation and MCSs at both seasonal and diurnal scales, although biases exist in aspects such as precipitation intensity and MCS frequency, size, propagation speed, and associated precipitation intensity. The 3‐km simulation using MYNN scheme generally outperforms the other simulations in capturing seasonal and diurnal precipitation over the mountain, while both it and the 4‐km simulation demonstrate superior performance in the western Amazon Basin, based on the comparison to the gridded precipitation products and gauge data. Dynamic factors, primarily low‐level jet and terrain‐induced uplift, are the key drivers for precipitation and MCS genesis along the east slope of the Andes, while thermodynamic factors control the precipitation and MCS activity in the western Amazon Basin and over elevated mountainous regions. The study suggests model improvements and better model configurations for future regional climate projections.

Meteorology & Atmospheric Sciences↗
DOE OSTI2023

Use-Inspired, Process-Oriented GCM Selection: Prioritizing Models for Regional Dynamical Downscaling

Dynamical downscaling is a crucial process for providing regional climate information for broad uses, using coarser-resolution global models to drive higher-resolution regional climate simulations. The pool of global climate models (GCMs) providing the fields needed for dynamical downscaling has increased from the previous generations of the Coupled Model Intercomparison Project (CMIP). However, with limited computational resources, the need for prioritizing the GCMs for subsequent downscaling studies remains. GCM selection for dynamical downscaling should focus on evaluating processes relevant for providing boundary conditions to the regional models and be inspired by regional uses such as the response of extremes to changes in the boundary conditions. This leads to the need for metrics representing processes of relevance to diverse stakeholders and subregions of a domain. Procedures to account for metric redundancy and the statistical distinguishability of GCM rankings are required. Further, procedures for selecting realizations from ensembles of top-performing GCM simulations can be used to span the range of climate change signals in multiple ways. As a result, distinct weighting of metrics and prioritization of particular realizations may depend on user needs. We provide high-level guidelines for such region-specific evaluations and address how CMIP7 might enable dynamical downscaling of a representative sample of high-quality models across representative shared socioeconomic pathways (SSPs).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Effects of Lower Troposphere Vertical Mixing on Simulated Clouds and Precipitation Over the Amazon During the Wet Season

Planetary boundary layer (PBL) schemes parameterize unresolved turbulent mixing within the PBL and free troposphere (FT). Previous studies reported that precipitation simulation over the Amazon in South America is quite sensitive to PBL schemes and the exact relationship between the turbulent mixing and precipitation processes is, however, not disentangled. In this study, regional climate simulations over the Amazon in January–February 2019 are examined at process level to understand the precipitation sensitivity to PBL scheme. The focus is on two PBL schemes, the Yonsei University (YSU) scheme, and the asymmetric convective model v2 (ACM2) scheme, which show the largest difference in the simulated precipitation. During daytime, while the FT clouds simulated by YSU dissipate, clouds simulated by ACM2 maintain because of enhanced moisture supply due to the enhanced vertical moisture relay transport process: (a) vertical mixing within PBL transports surface moisture to the PBL top, and (b) FT mixing feeds the moisture into the FT cloud deck. Due to the thick cloud deck over Amazon simulated by ACM2, surface radiative heating is reduced and consequently the convective available potential energy is reduced. As a result, precipitation is weaker from ACM2. Two key parameters dictating the vertical mixing are identified, p , an exponent determining boundary layer mixing and λ , a scale dictating FT mixing. Sensitivity simulations with altered p , λ , and other treatments within YSU and ACM2 confirm the precipitation sensitivity. The FT mixing in the presence of clouds appears most critical to explain the sensitivity between YSU and ACM2.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Daytime urban heat stress in North America reduced by irrigation

There is considerable uncertainty regarding the impact of irrigation on heat stress, partly stemming from the choice of heat stress index. Moreover, existing simulations are at scales that cannot appropriately resolve population centres or clouds and thus the potential for human impacts. Using multi-year convection-permitting and urban-resolving regional climate simulations, we demonstrate that irrigation alleviates summertime heat stress across more than 1,600 urban clusters in North America. This holds true for most physiologically relevant heat stress indices. The impact of irrigation varies by climate zone, with more notable irrigation signals seen for arid urban clusters that are situated near heavily irrigated fields. Through a component attribution framework, we show that irrigation-induced changes in wet-bulb temperature, often used as a moist heat stress proxy in the geosciences, exhibit an opposite sign to the corresponding changes in wet bulb globe temperature—a more complete index for assessing both indoor and outdoor heat risk—across climate zones. In contrast, the local changes in both wet-bulb and wet bulb globe temperature due to urbanization have the same sign. Here, our results demonstrate a complex relationship between irrigation and heat stress, highlighting the importance of using appropriate heat stress indices when assessing the potential for population-scale human impacts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A Simple Lagrangian Parcel Model for the Initiation of Summer-time Mesoscale Convective Systems over the Central United States

Mesoscale convective systems (MCSs) account for more than 50% of summertime precipitation over the central United States and have a significant impact on local weather and hydrologic cycle. It is hypothesized that the inadequate treatment of MCSs is responsible for the long-standing warm and dry bias over the central United States in coarse-resolution general circulation model (GCM) simulations. In particular, a better understanding of MCS initiation is still lacking. Here a single-column Lagrangian parcel model is first developed to simulate the basic features of a rising parcel. This simple model demonstrates the collective effects of boundary layer moistening and dynamical lifting in triggering convective initiation and reproduces successfully its early afternoon peak with surface equivalent potential temperature as a controlling factor. It also predicts that convection is harder to trigger in the future climate under global warming, consistent with the results from convection-permitting regional climate simulations. Then, a multicolumn model that includes an array of single-column models aligned in the east–west direction and incorporates idealized cold pool interaction mechanisms is developed. The multicolumn model captures readily the cold pool–induced upscale growth feature in MCS genesis from initially scattered convection that is organized into a mesoscale cluster in a few hours. It also highlights the crucial role of lifting effects due to cold pool collision and spreading, subsidence effect, and gust front propagation speed in controlling the final size of mesoscale clusters and cold pool regions. This simple model should be useful for understanding fundamental mechanisms of MCS initiation and providing guidance for improving MCS simulations in GCMs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Evaluating Contemporary and Future-Scenario Substantial-Precipitation Events in the Missouri River Basin Using Object-Oriented Analysis

The Missouri River Basin is the largest single river basin in the United States, and, as such, it plays an important role in natural ecosystems as well as the country’s economy, through agriculture, hydroelectric power generation, and transportation. Episodes of heavy precipitation can have a substantial negative impact on all these aspects of the basin, so understanding how well these episodes are simulated and projected to change in the future climate is important. We analyzed contemporary and projected mid-century behavior of heavy-precipitation episodes using an object-oriented analysis to diagnose short-term (≥5-day) and extended-period (≥30-day) events with substantial precipitation, using PRISM gridded, observed precipitation and RegCM4 regional-climate simulations that used outputs from two different GCMs for boundary conditions. The simulations were produced for the North American portion of the CORDEX program. A 25 km grid was used for the simulations and for aggregated PRISM precipitation. Overall, the simulated contemporary-climate events compared favorably with the PRISM events’ frequency and duration. The simulated event areas tended to be larger than the areas in the PRISM events, suggesting that the effective resolution of the simulations is greater than 25 km. Event areas and durations change little going from contemporary to scenario climate. The short-term events increase in frequency by an amount commensurate with the increase in mean precipitation simulated for the basis. However, the extended-term events showed little change in frequency, despite the average precipitation increase. Roughly half the extended-period events overlapped with at least one short-term event in both the observations and the simulations. Extended-period events that overlap a short-term event generally have larger areas and longer durations compared to their counterparts with no overlapping short-term events. Understanding the climate dynamics yielding the two types of extended-period events could be useful for assessing future changes in the Missouri River Basin’s heavy precipitation events and their impact.

Fisel, Brandon J.↗
DOE OSTI2024

Summer Convective Precipitation Changes Over the Great Lakes Region Under a Warming Scenario

To understand future summer precipitation changes over the Great Lakes Region (GLR), we performed an ensemble of regional climate simulations through the Pseudo-Global Warming (PGW) approach. We found that different types of convective precipitation respond differently to the PGW signal. Isolated deep convection (IDC), usually concentrated in the southern domain, shows an increase in precipitation to the north of the GLR. Mesoscale convective systems (MCSs), usually concentrated upwind of the GLR, shift to the downwind region with increased precipitation. Thermodynamic variables such as convective available potential energy (CAPE) and convective inhibition energy (CIN) are found to increase across almost the entire studied domain, creating a potential environment more favorable for stronger convection systems and less favorable for weaker ones. Meanwhile, changes in the lifting condensation level (LCL) and level of free convection (LFC) show a strong correlation with variations in convective precipitation, highlighting the significance of these thermodynamic factors in controlling precipitation over the domain. Our results indicate that the decrease in LCL and LCF in areas with increased convective precipitation is mainly due to increased atmospheric moisture. In response to the prescribed warming perturbation, MCSs occur more frequently downwind, while localized IDCs exhibit more intense rain rates, longer durations, and larger rainfall area.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Antecedent Hydrometeorological Conditions of Wildfire Occurrence in the Western U.S. in a Changing Climate

Abstract Wildfires have significant hydrological and ecological impacts in the western U.S. Using a high‐resolution regional climate simulation and wildfire observations for 1984–2018, this study investigates the antecedent hydrometeorological conditions (AHCs) of wildfires in the western U.S. During the warm season (April‐September), the wildfire AHCs feature diverse surface pressure (PS), soil moisture, and longwave/shortwave radiation (LW/SW) conditions. K‐means clustering classifies wildfires into four types with distinct AHCs: low‐PS‐type and high‐PS‐type with lower and higher PS anomalies, respectively, LW‐type featuring intense LW but weak SW anomalies, and wet‐soil‐type with wet soil anomalies. Each fire cluster represents 22%–27% of all the wildfires, featuring different combinations of climate and vegetation conditions and their diverse relations to regional hydrometeorological conditions, with wet‐soil‐type fires often exhibiting opposite correlations with AHCs compared to those of the other three types. In five major Köppen climate zones over the western U.S., clustering‐based predictions improve the seasonal wildfire prediction accuracy ( R 2 ) by 10% compared to prediction without classification. Such improvement comes from separating the opposite relationships between wet‐soil‐type fires and their seasonal AHCs from the other three types, along with separating LW‐type fires, which include most of the lightning‐ignited fires that occur more randomly. Increases in wildfire occurrence during 1984–2018 are dominated by the increases in the LW‐type fires, while the wet‐soil‐type fires have decreased, consistent with the long‐term drying in the western U.S.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Central American mountains inhibit eastern North Pacific seasonal tropical cyclone activity

Abstract The eastern North Pacific (ENP) has the highest density of tropical cyclones (TCs) on earth, and yet the controls on TCs, from individual events to seasonal totals, remain poorly understood. One effect that has not been fully considered is the unique geography of the Central American mountains. Although observational studies suggest these mountains can readily fuel individual TCs through dynamical processes, here we show that these mountains indeed play the opposite role on the seasonal timescale, hindering seasonal ENP TC activity by up to 35%. We found that these mountains significantly interrupt the abundant moisture transport from the Caribbean Sea to the ENP, limiting deep convection over the open ocean area where TCs preferentially occur. This study advances our fundamental understanding of ENP TC genesis mechanisms across the weather-to-climate timescales, and also highlights the importance of topography representation in improving the ENP regional climate simulations, as well as TC seasonal predictions and future projections.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Response of U.S. West Coast Mountain Snowpack to Local Sea Surface Temperature Perturbations: Insights from Numerical Modeling and Machine Learning

Sea surface temperature (SST) significantly modulates the precipitation and temperature over land, with important consequences on land surface processes such as snowpack. Compared to the impact of remote SST, the effect of nearshore/local SST is less well understood. In this study, the impact of local SST on the mountain snowpack of the U.S. West Coast is investigated using two 6-km regional climate simulations driven by the same lateral boundary conditions but with time-varying versus time-invariant and warmer local SSTs during 2003–15. Results show that local SST warming leads to warmer winter with more precipitation over the mountains. Meanwhile, the removal of SST temporal variability results in reduced temperature variability but increased precipitation variability. As a result, winter snow accumulation decreases by ~200 mm per season in the Cascade Mountains in the north but increases by ~100 mm per season in the Sierra Nevada in the south. Such a dipole response results from the competing effects of precipitation and temperature change at different elevations and are amplified by the enhanced atmospheric river moisture transport. To further delineate the relative contributions of different meteorological factors to the snowpack response, two neural network models were developed to predict the snow behaviors at seasonal and monthly scales. These models reveal the dominant influence of the total amount and the average temperature of precipitation on the snowpack response. Furthermore, these findings highlight the sensitivity of mountain snowpack to local SST in the western United States and underscore the importance of local SST and atmospheric rives to accurate snowpack estimations for water management.

54 ENVIRONMENTAL SCIENCES↗