Engineering topics

Rhoades, Alan

Publications and source records attributed to Rhoades, Alan.

17 records

DOE OSTI2026

High-resolution mountain topography can inform global snow vulnerability estimates

Snow is changing globally. Computationally intensive snow reanalysis products and downscaled climate model projections allow for the estimation of historical and projected changes in snow over ∼4–10 km resolutions, but these resolutions are coarse relative to the scales needed for water supply and flood planning. Fine-scale digital elevation models (DEMs) are widely available but are underutilized to make first-order assessments of snow vulnerability. Here, we leverage DEMs at a 7.5 arc s (∼250 m) resolution, combining these with historical freezing level height estimates from ERA-5 to derive estimates of changes in the snow-receiving area (SRA) and its variability across global mountain ranges. Results show estimated SRA declines in 29% (1.9 million km2) of the global mountain area from 1982–2020; 66% of the mountainous areas had no change over the historical period. At +1.5 °C of warming relative to the pre-industrial control, global mountain SRA would decline by 9.5% (1.0 million km2) relative to recent conditions. This loss would be approximately doubled with +2 °C of warming. In a +4 °C warming scenario, an additional 34% (3.6 million km2) of SRA would be lost beyond the +2 °C case. Across individual mountain ranges, SRA losses can occur nonlinearly with warming, with some locations that have historically had relatively minor SRA losses at risk of substantially larger losses in warmer climates. Analysis using coarser-resolution DEMs can underestimate or overestimate SRA and its rate of loss, with the largest impacts in relatively warm, low-elevation mountain ranges. Results of this work provide estimates of projected loss in SRA at policy-relevant warming levels; inform the resolutions needed for process-based snow modeling; identify snow vulnerability hotspots; and provide a new integrated approach to snow vulnerability assessment that is achievable at global scales and highlights potential nonlinearities from recent trends to a variety of future warming scenarios.

climate, mountains↗

DOE OSTI2025

Climate models show colorado drying sooner and with greater certainty east of the continental divide

Many studies have examined the aridity of the Colorado River Basin and the possible impacts of climate change which could further strain already over-allocated water resources in the region. Fewer studies have examined the multiple Colorado Rocky Mountain headwater regions specifically. This is especially true of areas East of the Continental Divide, despite water originating there being critical to cities and agriculture in Eastern Colorado and further downstream. This paper explores and compares drying trends in the Eastern and Western Colorado Rocky Mountains using single-model initial-condition large ensembles from ten global climate models. The use of multiple models allows us to identify signals that are consistent across different physics parameterizations, model grids, and other model intricacies. The large ensembles also allow us to quantify the time of emergence of these climate change signals--that is, when did (or when will) the long term change due to anthropogenic greenhouse gasses exceed the internal variability of the climate system. Consistent with previous studies, we find evidence of drying on both sides of the Continental Divide. That drying is more pronounced, occurs sooner, and is more consistent across global climate models in the East, however, highlighting the region’s importance despite generally receiving less attention than the West.

Rugg, Allyson↗

DOE OSTI2025

Understanding drivers and uncertainty in projected African precipitation

We investigate the drivers of projected summer precipitation changes and their uncertainties across Africa in the second half of the 21st century under the SSP2-4.5 and SSP5-8.5 scenarios using CMIP6 models. Our results reveal distinct regional precipitation changes, particularly under SSP5-8.5, with robust increases of 75%, 24%, and 17% over the Sahara, South Eastern Africa, and South Central Africa, respectively, and a decline of up to 5% over West Southern Africa (WSAF). In most regions, precipitation increases are driven by enhanced vertical thermodynamic processes associated with temperature-induced moisture increases and enhanced moisture convergence. In contrast, the WSAF decrease is associated with vertical dynamic processes driven by a weakening of the Hadley circulation’s ascending branch. Model uncertainty accounts for over 85% of total projection uncertainty across all regions and is largely due to subgrid-scale parameterizations. Overall, this study enhances our understanding of climate change impacts on African precipitation.

Taguela, Thierry N↗

DOE OSTI2024

Are Atmospheric Models Too Cold in the Mountains? The State of Science and Insights from the SAIL Field Campaign

Mountains play an outsized role in water resource availability, and the amount and timing of water they provide depend strongly on temperature. To that end, we ask the question: How well are atmospheric models capturing mountain temperatures? We synthesize results showing that high-resolution, regionally relevant climate models produce 2-m air temperature (T2m) measurements colder than what is observed (a “cold bias”), particularly in snow-covered midlatitude mountain ranges during winter. We find common cold biases in 44 studies across global mountain ranges, including single-model and multimodel ensembles. We explore the factors driving these biases and examine the physical mechanisms, data limitations, and observational uncertainties behind T2m. Our analysis suggests that the biases are genuine and not due to observation sparsity or resolution mismatches. Cold biases occur primarily on mountain peaks and ridges, whereas valleys are often warm biased. Our literature review suggests that increasing model resolution does not clearly mitigate the bias. By analyzing data from the Surface Atmosphere Integrated Field Laboratory (SAIL) field campaign in the Colorado Rocky Mountains, we test various hypotheses related to cold biases and find that local wind circulations, longwave (LW) radiation, and surface-layer parameterizations contribute to the T2m biases in this particular location. We conclude by emphasizing the value of coordinated model evaluation and development efforts in heavily instrumented mountain locations for addressing the root cause(s) of T2m biases and improving predictive understanding of mountain climates.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Are Atmospheric Models Too Cold in the Mountains? The State of Science and Insights from the SAIL Field Campaign

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Using Temporal Deep Learning Models to Estimate Daily Snow Water Equivalent Over the Rocky Mountains

Abstract In this study we construct and compare three different deep learning (DL) models for estimating daily snow water equivalent (SWE) from high‐resolution gridded meteorological fields over the Rocky Mountain region. To train the DL models, Snow Telemetry (SNOTEL) station‐based SWE observations are used as the prediction target. All DL models produce higher median Nash‐Sutcliffe Efficiency (NSE) values than a conceptual SWE model and interpolated gridded data sets, although mean squared errors also tend to be higher. Sensitivity of the SWE prediction to the model's input variables is analyzed using an explainable artificial intelligence (XAI) method, yielding insight into the physical relationships learned by the models. This method reveals the dominant role precipitation and temperature play in snowpack dynamics. In applying our models to estimate SWE throughout the Rocky Mountains, an extrapolation problem arises since the statistical properties of SWE (e.g., annual maximum) and geographical properties of individual grid points (e.g., elevation) differ from the training data. This problem is solved by normalizing the SWE with its historical maximum value to alleviate extrapolation for all tested DL models. Our work shows that the DL models are promising tools for estimating SWE, and sufficiently capture relevant physical relationships to make them useful for spatial and temporal extrapolation of SWE values.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Surface Atmosphere Integrated Field Laboratory (SAIL) (Field Campaign Report)

Mountains are the natural water towers of the world, effectively turning water vapor into readily available fresh water through precipitation, snowpack, and runoff. They contribute disproportionately to precipitation over land, but are under-observed, leading to large gaps in the scientific understanding of convection, extreme precipitation and weather, and interactions between atmospheric circulation, radiation, and land-surface conditions. The mountain hydrometeorology community has repeatedly called for integrated atmospheric and land observations of water and energy budgets in complex terrain that span these scales to establish benchmarks against which scale-dependent models can be further developed.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

LSTM-Based Data Integration to Improve Snow Water Equivalent Prediction and Diagnose Error Sources

Accurate prediction of snow water equivalent (SWE) can be valuable for water resource managers. Recently, deep learning methods such as long short-term memory (LSTM) have exhibited high accuracy in simulating hydrologic variables and can integrate lagged observations to improve prediction, but their benefits were not clear for SWE simulations. Here we tested an LSTM network with data integration (DI) for SWE in the western United States to integrate 30-day-lagged or 7-day-lagged observations of either SWE or satellite-observed snow cover fraction (SCF) to improve future predictions. SCF proved beneficial only for shallow-snow sites during snowmelt, while lagged SWE integration significantly improved prediction accuracy for both shallow- and deep-snow sites. The median Nash–Sutcliffe model efficiency coefficient (NSE) in temporal testing improved from 0.92 to 0.97 with 30-day-lagged SWE integration, and root-mean-square error (RMSE) and the difference between estimated and observed peak SWE values d max were reduced by 41% and 57%, respectively. DI effectively mitigated accumulated model and forcing errors that would otherwise be persistent. Moreover, by applying DI to different observations (30-day-lagged, 7-day-lagged), we revealed the spatial distribution of errors with different persistent lengths. For example, integrating 30-day-lagged SWE was ineffective for ephemeral snow sites in the southwestern United States, but significantly reduced monthly-scale biases for regions with stable seasonal snowpack such as high-elevation sites in California. These biases are likely attributable to large interannual variability in snowfall or site-specific snow redistribution patterns that can accumulate to impactful levels over time for nonephemeral sites. These results set up benchmark levels and provide guidance for future model improvement strategies.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Dataset for 'Ombadi et al. (2023). A warming-induced reduction in snow fraction amplifies rainfall extremes, Nature'

This package contains the main codes, sample input data and main result files to reproduce the analysis and results presented in the article: “Ombadi et al. (2023), A warming-induced reduction in snow fraction amplifies rainfall extremes, Nature”. The folder consists of the following: (1) “Raw data”: a folder that contains sample input data which is used in some of the codes for demonstration purposes. It also contains data that was not pre-processed such as Elevation data; (2) “Results”: this folder contains files of the main results presented in the paper including: “Annual-Max-Series”, “Change-rainfall-extremes”, “Change-snow-fraction”, “Warming levels_By scenario_model_year” and “Masks”. Description of these folders is detailed in the "Readme.rtf" file; (3) Python jupyter notebooks (Extract_Annual Max Series (AMS).ipynb, Elevation-dependent amplification of rainfall extremes.ipynb, Sensitivity_to_global_warming.ipynb) demonstrate the main steps of analysis. Further description of those notebooks is provided in the "Readme.rtf" file; (4) R code for extreme value analysis (Extreme_Value_Analysis.R). The sample and pre-processed dataset in "Raw data" is obtained from publicly available repositories of CMIP6 and ERA5 datasets; see Methods for more detail. This research was supported by Office of Science, Office of Biological and Environmental Research of the US Department of Energy under contract no. DE-AC02-05CH11231 for the CASCADE Scientific Focus (funded by the Regional and Global Model Analysis Program area within the Earth and Environmental Systems Modeling Program) and the iNAIADS Early Career Research Project (funded by the Environmental Systems Science program).

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

The Role of Atmospheric Rivers on Groundwater: Lessons Learned From an Extreme Wet Year

Abstract In the coastal regions of the western United States, atmospheric rivers (ARs) are associated with the largest precipitation generating storms and contribute up to half of annual precipitation, but the impact of ARs on the integrated hydrologic cycle, specifically on groundwater storage and hydrodynamics, is largely unknown. To better explore the hydrologic behavior of AR versus non‐AR event precipitation, we present a novel combination of two water tracking methods (one in the atmosphere and one in the subsurface) to explicitly track the full lifecycle of water parcels generated by ARs. Simulations of northern California's Cosumnes River watershed during the record wet 2017 water year are performed via the coupling of a high‐resolution regional climate model and a land surface‐groundwater model accounting for lateral groundwater flow. Despite ARs contributing more precipitation than non‐AR storms, we find less AR water is preferentially stored in aquifers by year end. Fractionally, ARs result in 300% less snow derived groundwater‐recharged compared to non‐AR precipitation. Rain‐on‐snow (RoS) plays an important role in AR‐driven discharge, where over 50% of total discharge from ARs snow is from RoS events. Finally, despite record‐breaking annual precipitation, simulated groundwater depletion occurs by year end due to estimates of groundwater pumping activities. The results from these simulations serve as a partial analogue of future hydrologic conditions where ARs are expected to intensify and provide a greater fraction of annual precipitation due to climate change.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Amplified drought trends in Nepal increase the potential for Himalayan wildfires

Here we report in the spring 2021, Nepal underwent a record wildfire season in which active fires were detected at a rate 10 times greater than the 2002–2020 average. Prior to these major wildfire events, the country experienced a prolonged precipitation deficit and extreme drought during the post-monsoon period (starting in October 2020). An analysis using observational, reanalysis, and climate model ensemble data indicates that both climate variability and climate change-induced severe drought conditions were at play. Further analysis of climate model outputs suggests the likely reoccurrence of drought conditions, thus favoring active wildfire seasons in Nepal throughout the twenty-first century. While the inter-model uncertainty is large and direct modeling of wildfire spread and suppression has not been completed, the demonstrated relationship between a drought index (the standardized precipitation and evapotranspiration index) and subsequent fire activity may offer actionable opportunities for forest managers to employ the monitoring and projection of climate anomalies at sub-seasonal to decadal timescales to inform their management strategies for Nepal’s wildlands.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Data used for Figure 3 of the Nature Reviews Earth and Environment (NREE) paper: "A low-to-no snow future and its impacts on water resources in the western United States

In order to synthesize western United States snowpack projections, this dataset contains the results of 18 peer-review journal articles over three periods of interest (2025-2049, 2050-2074, and 2075-2099) and over 4 mountain ranges (Cascades, Sierra Nevada, Rockies, Wasatch/Uinta) in addition to western-US wide projections. Distinction is made by model type (Earth System Models, bias-corrected statistically downscaled Earth System Models, and regional climate models). RCP4.5 and RCP8.5 emission scenarios are considered. Percent snow water equivalent (SWE) loss considers 1 April, peak SWE and/or seasonal SWE. Heterogeneity in projected snowpack changes exists across mountain ranges and for different modeling approaches, but generally indicate agreement in decreases by the end of the century.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

A mechanistic understanding of oxygen isotopic changes in the Western United States at the Last Glacial Maximum

At the Last Glacial Maximum (LGM), records suggest drier conditions in the northwest United States and wetter conditions in the southwest United States relative to present-day as well as widespread changes in the isotopic composition of water. However, the mechanisms responsible for these changes remain ambiguous. Here, we explore differences in western United States hydroclimate between the LGM and preindustrial with a water isotope tracer enabled Earth System Model. We then use proxy forward models to compare simulated and recorded δ 18 O in speleothems. We find that the pattern of hydroclimate response in the western United States at the LGM relates to a combination of 1) increased frequency and southward shifted wintertime extratropical cyclones in the North Pacific, 2) greater rainout of moisture as it moves over the continent, and 3) reduced evaporation in the cooler LGM climate. The simulated lower δ 18 O of precipitation at the LGM relates predominantly to an increase in cool season moisture removal efficiency, with a secondary contribution from relatively more cool season precipitation. Both surface temperatures and North American ice sheets contribute to these hydroclimate changes at the LGM. Comparisons between δ 18 O from proxy forward models and speleothem records in the western United States show general agreement at the LGM, with increasing depletion moving towards the continental interior. Finally, this study highlights the similarities and differences between hydrologic and δ 18 O changes at the LGM and emphasizes the utility of model-proxy comparison for interpretation.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Evaluating Variable‐Resolution CESM Over China and Western United States for Use in Water‐Energy Nexus and Impacts Modeling

Abstract Climate models are critical tools to study earth system processes and are often further downscaled to provide refined‐resolution data sets for regional water‐energy impact analyses. In this study, the Variable‐Resolution Community Earth System Model (VR‐CESM), a new dynamical downscaling method, is employed to generate 1/8° simulations for China and the western United States over 1970–2006. Precipitation, temperature, snowpack, radiation, and wind speed are compared with two 1° global climate models, regional climate model ensemble, and reference data sets to evaluate their fidelity. Simulated precipitation skill is dependent upon the reference data sets used, a slight (<+10%) and spatially consistent wet bias in China during summer and the western United States during winter. Precipitation bias in VR‐CESM is the smallest among models and demonstrates the added value provided by resolution refinement. However, temperature biases are generally greater than precipitation due to more complicated interactions among processes. VR‐CESM simulates +2°C warm bias in the low‐elevation regions in China but −2°C cold bias in mountains. However, DJF cold bias (−2.23°C to−3.37°C) and JJA warm bias (+1.13°C to +1.87°C) exist in the western United States. The modeling skills of precipitation and temperature at various topographies are similar in the two evaluated regions. The cold bias was likely affected by downward shortwave radiation deficit (−15%) and slow bias of surface wind speeds (−16% to −29%). Despite the wet and cold bias, snow water equivalent and snow cover are underestimated in mountains. We conclude that VR‐CESM outperforms other models evaluated and provides the best estimation of downscaled climate model data for water‐energy studies.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Machine learning to generate gridded extreme precipitation data sets for global land areas with limited in situ measurements

There is a strong need for gridded observational data sets that describe the climatology of extreme precipitation across the globe; such data sets are critical for quantifying precipitation extremes and their corresponding influence on large perturbations on surface and groundwater systems and flooding. Machine learning (ML) methods can be used to generate pseudo in situ measurements of the climatology of extreme precipitation for land regions with poor sampling by matching relevant geographic (e.g., orography) and atmospheric (e.g., surface temperature, precipitation, and pressure) variables for densely sampled regions and those with limited geographic sampling.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Surface Atmosphere Integrated Field Laboratory (SAIL) Science Plan

Mountains are the natural water towers of the world, effectively turning water vapor into readily available fresh water through precipitation, snowpack, and runoff. Unfortunately, Earth system models (ESMs) have persistently been unable to predict the timing and availability of water resources from mountains because the source(s) of model error are difficult to isolate in complex terrain with limited atmospheric or land-surface observations. Further complications arise from the gross scale mismatch between ESM grid box sizes and the relevant scales of mountainous hydrological processes. The mountain hydrometeorology community has repeatedly called for integrated atmospheric and land observations of water and energy budgets in complex terrain that span these scales to establish benchmarks against which scale-dependent models can be further developed.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Sources of Subseasonal-To-Seasonal Predictability of Atmospheric Rivers and Precipitation in the Western United States

Atmospheric rivers (ARs) account for a large portion of winter precipitation in the western US. To evaluate the sources of AR and precipitation predictability at subseasonal-to-seasonal timescales, we examine the relationships between two climate modes, the El Niño/Southern Oscillation (ENSO) and Madden-Julian oscillation (MJO), and winter hydroclimate in the western US. Our analysis uses a large ensemble of the Weather Research and Forecast (WRF) model simulations from 1981 to 2017, facilitating the assessment of uncertainty in climate mode-AR relationship due to climate variability and short length of observations. Over the North Pacific basin, we find ENSO-related latitudinal shifts of ARs, but ENSO has little effect on basin-wide averaged AR frequency. Over the western US, there is some uncertainty in the ENSO-AR connection due to the impacts of ENSO indices and datasets. However, extreme El Niño events defined by the ENSO longitude index are consistently linked to increased landfalling AR activity. The MJO can lead to significantly enhanced or suppressed landfalling AR activity in the large-ensembles, depending on the phase of MJO and time lag, while observations are too short in length to robustly show this signal. ENSO substantially modulates the MJO-AR relationship, triggering variegated responses of landfalling ARs and AR precipitation in La Niña and El Niño years. Our findings highlight the need to evaluate concurrent effects of different climate modes on ARs and precipitation, and may shed light on a path toward more accurate subseasonal-to-seasonal prediction of ARs and precipitation over the western US.

54 ENVIRONMENTAL SCIENCES↗