Search NASA⌕ Search

SEARCH · Search NASA

Results for “mountainous regions”

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 OSTI2026

Evaluating ecosystem water use efficiency under drought stress: a case study of the Helan Mountain region, northwest China

Context Water use efficiency (WUE) is a fundamental ecological indicator links carbon assimilation and water loss in terrestrial ecosystems. Understanding its responses to drought stress is essential for adaptive ecosystem management, particularly in climate-sensitive mountain landscapes. Objectives This study aimed to investigate drought-driven variations in WUE across major vegetation types in the Helan Mountain region of Northwest China. Specifically, we sought to identify dominant ecological drivers of WUE variability and to disentangle their relative importance and causal pathways. Methods We quantified WUE using the Moderate Resolution Imaging Spectroradiometer (MODIS) products and the Drought Severity Index (DSI) data from 2001 to 2020. To examine WUE – drought relationships across contrasting vegetation types, we employed a spatially explicit analytical framework integrating Random Forest (RF) modeling, partial correlation analysis, and structural equation modeling (SEM). Results Regional WUE exhibited relatively stable interannual dynamics, yet pronounced spatial heterogeneity that was strongly modulated by drought conditions. Vegetation properties, particularly Leaf Area Index (LAI) and Normalized Difference Vegetation Index (NDVI), emerged as the dominant determinants of WUE, with NDVI alone explaining over 20% of its spatial variance in forest and grassland during non-drought periods. SEM analyses revealed that climate forcing influenced WUE mainly through indirect pathways mediated by soil moisture availability and vegetation structural dynamics, rather than through direct climatic controls. Among all regulating factors, LAI acted as the central control node governing ecosystem carbon–water coupling. In contrast, short-term climatic stress, especially atmospheric demand and drought duration, exerted weak or negative direct effects on WUE. Ecosystem-specific responses were observed, with croplands mainly regulated by soil water availability, whereas forests and grasslands showed more sensitive to atmospheric drought stress. Together, these results reveal a hierarchical control framework where soil–vegetation interactions mediate climate impacts on WUE, driving strong spatial heterogeneity in drought responses across mountain landscapes. Conclusions Our findings highlight the pivotal role of indirect drought effects mediated by vegetation and soil processes in shaping ecosystem WUE. The identified soil–vegetation–climate regulatory hierarchy provides mechanistic insight into landscape–scale drought sensitivity and supports integrated modeling approaches for evaluating ecosystem resilience and sustainable management in arid mountain regions.

China↗
DOE OSTI2024

Characterization of Most Promising Sequestration Formations in the Rocky Mountain Region

The project Characterization of Most Promising Sequestration Formations in the Rocky Mountain Region is one of 9 site characterization projects that were implemented as part of ARRA (American Recovery and Reinvestment Act). Data from this project was used to improve resolution of data in NATCARB in the area of study. Data related to this study has already been incorporated in NATCARB Atlas. The Rocky Mountain Carbon Capture and Storage (RMCCS) project investigated multiple geologic formations and characterized a local site on the Colorado Plateau for future CCS opportunities. The RMCCS project focused on the Cretaceous Dakota, Jurassic Entrada, and Pennsylvanian Weber Sandstones, the three largest regional formations. All formations in this project are potential CO2 storage resources for future power plants, natural gas processing plants, cement plants, and oil shale development projects. The area adjacent to Craig, Colorado, (Sand Wash Basin) was the area selected for detailed geologic characterization on the RMCCS project. The basin was selected in part because the geology can be extrapolated to other sites on the Colorado Plateau. Field mapping and seismic surveys were conducted to identify and evaluate the basin's structural configuration. A 9,745-foot deep characterization well was drilled to collect 131 feet of core and a suite of geophysical well log data. Petrophysical tests on samples of core were used to calibrate geophysical log data, which can be used to obtain storage resource estimates and evaluate associated uncertainty as well as simulate the hydrologic behavior of injected CO2. A detailed analysis of the primary formations (Dakota, Entrada and Weber sandstones) yielded a more accurate CO2 storage resource assessment for these formations within the Colorado Plateau; RMCCS estimates indicate a total CO2 storage resource of more than 38,000 million metric tons. The characterization of the Sand Wash Basin (2-D seismic surveys, multiple well logs and lithological, petrophysical and geochemical analyses) allowed for a detailed 3-D model to be constructed. The model served as the framework for analyses ranging from CO2 storage resource, injectivity, and subsurface flow to uncertainty estimates to evaluation of risk.

2-D seismic↗
DOE OSTI2021

Weather radar in Nepal: opportunities and challenges in a mountainous region

Extreme rainfall is one of the major causes of natural hazards in the central Himalayan region, including Nepal. The performance of strategies to manage hazards and related risks relies on the accuracy of quantitative hydrometeorological prediction. Rain gauges have traditionally been used to measure the rainfall amount. However, point measurements with limited gauge coverage cannot accurately represent spatial precipitation variability in complex topography. Weather radar have shown potential for useful information on accurate areal rainfall estimates. The Department of Hydrology and Meteorology (DHM) in Nepal installed their first weather radar in 2019 in the western region of the country. Two more radars will be added to the planned radar network in the near future, in the country's central and eastern regions, respectively. We highlight both the opportunities and challenges with radar installation and observation in the mountainous regions. Radar-rainfall estimates across the Himalayas can be useful to inform decision-making in a broad range of infrastructure sectors, including water, energy, construction, transportation, and agriculture.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Rayleigh Wave Propagation in the Bighorn Mountains Region, Wyoming

Short-period Rayleigh waves, Rg , provide strong constraints on the depth of shallow seismic events and are of interest for monitoring small explosions. Characterizing the seismic sources that generate Rg requires an understanding of how shallow crustal structure affects Rayleigh wave propagation. Here, in support of these efforts, this study utilizes observed waveforms from small shallow explosions recorded on temporary seismic network deployments in the Bighorn region, Wyoming. We study regional near-surface geology by measuring changes in surface-wave amplitude and polarization during propagation through basins, foothills, and mountains. We develop additional insight by carrying out surface-wave eigenfunction analyses and numerical-wave simulations, which together reproduce many characteristics seen in the observed waveforms. Our results show how sedimentary basins in the Bighorn region allow for amplified prograde-polarized higher-mode and retrograde-polarized fundamental-mode Rayleigh waves, whereas adjacent mountains only support retrograde motion. These different modes provide distinct constraints on the Earth structure and source characteristics, potentially enabling targeted inversions in future studies. Finally, our findings provide insight into Rg propagation through complex near-surface geology, improving our understanding of shallow propagation and source effects that are relevant to explosion monitoring efforts.

58 GEOSCIENCES↗
DOE OSTI2023

Seasonal forecasting skill for the High Mountain Asia region in the Goddard Earth Observing System

Seasonal variability of the global hydrologic cycle directly impacts human activities, including hazard assessment and mitigation, agricultural decisions, and water resources management. This is particularly true across the High Mountain Asia (HMA) region, where availability of water resources can change depending on local seasonality of the hydrologic cycle. Forecasting the atmospheric states and surface conditions, including hydrometeorologically relevant variables, at subseasonal-to-seasonal (S2S) lead times of weeks to months is an area of active research and development. NASA's Goddard Earth Observing System (GEOS) S2S prediction system has been developed with this research goal in mind. Here, we benchmark the forecast skill of GEOS-S2S (version 2) hydrometeorological forecasts at 1–3-month lead times in the HMA region, including a portion of the Indian subcontinent, during the retrospective forecast period, 1981–2016. To assess forecast skill, we evaluate 2 m air temperature, total precipitation, fractional snow cover, snow water equivalent, surface soil moisture, and terrestrial water storage forecasts against the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2) and independent reanalysis data, satellite observations, and data fusion products. Anomaly correlation is highest when the forecasts are evaluated against MERRA-2 and particularly in variables with long memory in the climate system, likely due to the similar initial conditions and model architecture used in GEOS-S2S and MERRA-2. When compared to MERRA-2, results for the 1-month forecast skill range from an anomaly correlation of R anom =0.18 for precipitation to R anom =0.62 for soil moisture. Anomaly correlations are consistently lower when forecasts are evaluated against independent observations; results for the 1-month forecast skill range from R anom =0.13 for snow water equivalent to R anom =0.24 for fractional snow cover. We find that, generally, hydrometeorological forecast skill is dependent on the forecast lead time, the memory of the variable within the physical system, and the validation dataset used. Overall, these results benchmark the GEOS-S2S system's ability to forecast HMA hydrometeorology.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Understanding and Predictability of Integrated Mountain Hydroclimate (Workshop Report)

Mountainous systems cover approximately 23% of Earth’s land and are distributed across all continents. They can capture and store atmospheric moisture that is then cycled through the terrestrial surface and subsurface system, released to downstream communities, and cycled back to the atmosphere. Mountain hydroclimate—characterized by steep gradients, geological, ecological, and biogeochemical diversity—is influenced by topographic forcing and elevated warming and susceptible to large subseasonal to multidecadal variability and rapid changes. Terrestrial hydrological and biogeochemical cycles also experience cascading effects from global warming impacts, such as multidecadal declines in mountain snowpack, longer growing seasons, and increased frequency and severity of extreme events like droughts and wildfires. However, little is known about the effects of these impacts and their feedbacks on climate systems and surface-subsurface compartments. Also unknown are the full implications of changing hydroclimate and extreme events on hydro biogeochemical cycles across atmosphere, terrestrial, and human systems in mountain regions and beyond. This knowledge gap is critical, given human reliance on mountain systems for stable water supply and quality. Mountain systems’ increasing vulnerability to climate change and human perturbations motivates the need to improve understanding of integrated mountain hydroclimate (IMHC) systems and their feedbacks and impacts on humans across scales. However, due to large heterogeneity and strong gradients, coupled natural-human processes in mountain regions present significant challenges for observations, modeling, predictions, and projections. Motivated by gaps in mountain hydroclimate understanding, observations, and modeling and the need for credible projections of future changes, the U.S. Department of Energy’s (DOE) Biological and Environmental Research (BER) program organized a virtual workshop on “Understanding and Predictability of Integrated Mountain Hydroclimate.” Sponsored by BER’s Earth and Environmental Systems Sciences Division (EESSD), the workshop aimed to inform and catalyze EESSD’s growing interests in enhancing predictive understanding of IMHC. Organizers structured the workshop to identify (1) knowledge gaps, (2) observational and modeling challenges, (3) short-term (1 to 3 years) and long-term (10 years and beyond) research opportunities, and (4) strategies for fostering collaboration and coordination. To address the outstanding challenges of IMHC, the workshop included two sessions organized by disciplinary, cross-disciplinary, and crosscutting science topics. The disciplinary and cross-disciplinary topics focused on essential IMHC elements: atmosphere, terrestrial, and human systems and their interactions. Breakout sessions on disciplinary and cross-disciplinary topics facilitated identification of crosscutting topics and central emerging themes. Session 1 focused on connecting existing DOE investments to accelerate progress related to scientific challenges in understanding mountain hydroclimate. In Session 2, participants further explored key Session 1 takeaways through the lens of multiagency collaborations and coordination.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Identifying snowfall elevation patterns by assimilating satellite-based snow depth retrievals

Precipitation in mountain regions is highly variable and poorly measured, posing important challenges to water resource management. Traditional methods to estimate precipitation include in-situ gauges, Doppler weather radars, satellite radars and radiometers, numerical modeling and reanalysis products. Each of these methods is unable to adequately capture complex orographic precipitation. Here, we propose a novel approach to characterize orographic snowfall over mountain regions. We use a particle batch smoother to leverage satellite information from Sentinel-1 derived snow depth retrievals and to correct various gridded precipitation products. This novel approach is tested using a simple snow model for an alpine basin located in Trentino Alto Adige, Italy. Here, we quantify the precipitation biases across the basin and found that the assimilation method (i) corrects for snowfall biases and uncertainties, (ii) leads to cumulative snowfall elevation patterns that are consistent across precipitation products, and (iii) results in overall improved basin-wide snow variables (snow depth and snow cover area) and basin streamflow estimates.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Explicit Representation of Orographic Anisotropy for All Directions Improves Nanling Mountain Rainfall Simulation

Climate models exhibit significant rainfall bias in mountainous regions. One reason is the insufficient representation of orographic anisotropy in these models. In this study, we implement the orographic drag scheme with 3-D orographic anisotropy (all flow directions (AFD)) into a general circulation model and investigate its impact on Nanling rainfall simulation where models have systematic dry bias in the summer. It is shown that the AFD alleviated the Nanling mountain rainfall bias by over 60%. This is through an anomalous “lower-convergence-higher-divergence” deceleration pattern of the flow windward of the Nanling Mountains that enhanced vertical motion and low-level moisture convergence. The results suggest the importance of explicit orographic anisotropy representation in rainfall simulation in mountainous regions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Impact of Orographic Drag Schemes on East Asia Rainfall

Current generation of climate models often has significant biases in mountainous regions where the gradient of elevation is steep, and the terrain is complex. Potential reasons for these biases include under‐representation of orographic drag process in climate models. In this study, we assess the impact of orographic drag on East Asia rainfall by comparing the impact of a new orographic drag scheme that considers 3D orographic anisotropy (3D‐oro) with a 2‐D scheme in a general circulation model. Two sets of simulations (medium‐range and seasonal forecast) are carried out for the comparison and validation against observation. It is shown that through local/remote forcing of the drag in the mountainous regions, the 3D‐oro alleviates part of the excessive rainfall in west Tibetan Plateau and parts of insufficient rainfall in Southeast China by about 25%∼50% in the January/winter forecasts; it alleviates about 25%∼50% of the rainfall bias in part of south Tibetan Plateau and of East Asia in the July/summer forecasts. The results suggest the importance of improved orographic drag process and its impact in climate modeling for those regions that are prone to significant impact of hydroclimate events.

orographic anisotropy↗
DOE OSTI2021

Seismic source characterization of the Arabian Peninsula and Zagros Mountains from regional moment tensor and coda envelopes

Reliable estimates of moment magnitude and source mechanism for seismic events in the Middle East can be challenging due to a small number of openly available stations, the complex tectonic setting, and regions of high attenuation. Access to high-quality waveform data from well-calibrated regional seismic stations is fundamental in producing robust and stable estimates of earthquake source parameters, particularly when measurements of absolute ground motion amplitudes are required. Earthquake source mechanism and moment magnitude are invaluable information in the assessment of seismic hazard, plate motions, and the characterization of faults and regional stress field. The expansion of the Kingdom of Saudi Arabia’s national seismic network in the past several decades provides an opportunity to develop the capabilities of routine focal mechanism and moment magnitude estimations. In this study, we use time-domain full waveform moment tensor inversion and coda envelope–derived amplitude measurements to solve for earthquake source mechanism, moment magnitudes, and their source type. Furthermore, we compared the moment magnitudes calculated from the two methods and publicly available earthquake catalogs and discuss the implications of the obtained source parameters.

58 GEOSCIENCES↗
DOE OSTI2023

Editorial: Advances in observations and modeling of snow, forest-snow processes and snow hydrology

Seasonal snowpack plays a pivotal role in the regional water cycle. In Western mountainous regions, where precipitation predominantly occurs as snow during winter and is limited in the summer months when water demand peaks, snowpack serves as a critical water source essential for sustaining water supply, ecosystems, agriculture, and industries. Additionally, snowpack dynamically influences extreme events, such as floods, droughts, and wildfires, through intricate land-atmosphere feedback. Nevertheless, significant knowledge gaps persist within this field due to data limitations and the complexity of the physics involved in modeling snow and snow hydrology (i.e., snow acts as a critical upstream hydrological process). Here, this research topic presents a collection of manuscripts that address some of these important research gaps.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Population Aging and Heat Exposure in the 21st Century: Which U.S. Regions Are at Greatest Risk and Why?

Abstract Background and Objectives The co-occurring trends of population aging and climate change mean that rising numbers of U.S. older adults are at risk of intensifying heat exposure. We estimate county-level variations in older populations’ heat exposure in the early (1995–2014) and mid (2050) 21st century. We identify the extent to which rising exposures are attributable to climate change versus population aging. Research Design and Methods We estimate older adults’ heat exposure in 3,109 counties in the 48 contiguous U.S. states. Analyses use NASA NEX Global Daily Downscaled Product (NEX-GDDP-CMIP6) climate data and county-level projections for the size and distribution of the U.S. age 69+ population. Results Population aging and rising temperatures are documented throughout the United States, with particular “hotspots” in the Deep South, Florida, and parts of the rural Midwest. Increases in heat exposure by 2050 will be especially steep in historically colder regions with large older populations in New England, the upper Midwest, and rural Mountain regions. Rising temperatures are driving exposure in historically colder regions, whereas population aging is driving exposure in historically warm southern regions. Discussion and Implications Interventions to address the impacts of temperature extremes on older adult well-being should consider the geographic distribution and drivers of this exposure. In historically cooler areas where climate change is driving exposures, investments in warning systems may be productive, whereas investments in health care and social services infrastructures are essential in historically hot regions where exposures are driven by population aging.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Low flow characteristics for regulated and unregulated streams in North Carolina and prediction using climate signals

In this work, low flow statistics of regulated and unregulated streams in the state of North Carolina were updated with streamflow data through 2019. About 22% of the streams considered show a significant downward trend, but considerable low-frequency variability confounding trends. The relationship between regional index time series of groundwater-depth low flow and the Atlantic Multidecadal Oscillation (AMO) was examined to determine whether low-frequency climate modes can account for the long-term pattern in low flows. Consequently, a significant correlation was found between AMO and groundwater-depth low flow, such that positive AMO is associated with lower groundwater-depth low flow and vice versa, particularly for the Piedmont region. Predictive equations for annual low flows at the ecoregion level shows that springtime average streamflow and AMO were selected as the primary predictors of low flow for coastal and Piedmont regions, whereas springtime average streamflow and the November–December–January average Oceanic Niño Index were used as the primary predictors for the mountain region. The relative root mean square error (RMSE) of the disaggregated predictions to US Geologic Survey gauge locations was <23% at 79% of the stations, between 24% and 43% at 10% of the stations, and greater than 44% at 1% of the stations. The remaining 10% of stations showed large RMSEs. This latter percentage is characterized by smaller drainage basins and intermittent flows, suggesting the prediction models are not applicable to drainage basins smaller than roughly 20 square km and intermittent streams.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Evaluating probabilistic deep learning methods for uncertainty quantification of temperature downscaling

Deep learning (DL) has emerged as a promising tool for downscaling coarse-resolution climate data to high-resolution outputs, enabling improved regional climate predictions. A critical aspect of DL-based downscaling is the incorporation of uncertainty quantification (UQ), which enhances the interpretability and reliability of predictions—key factors for climate risk assessment and decision-making. This study develops a DL model to downscale 2 m temperature across the contiguous United States using reanalysis datasets. We systematically evaluate three epistemic UQ methods—deep ensembles (DEns), Monte Carlo dropout (MCD), and Flipout—based on their probabilistic accuracy, downscaling performance, sensitivity to geographical features, and computational efficiency. Results indicate that MCD generally outperforms Flipout and DEns in terms of calibration and downscaling accuracy. However, DEns demonstrate lower calibration errors in coastal regions, indicating its higher confidence within these areas. Flipout, in contrast, is more sensitive to elevation gradients and exhibits higher calibration errors in mountainous regions. Hence, the choice of UQ method for this task depends on the specific requirements of the application. For applications that prioritize overall calibration, downscaling accuracy, and computational efficiency, MCD is a strong candidate. These findings highlight the importance of selecting UQ methods based on application-specific requirements, such as geographical context and computational constraints. By addressing the trade-offs between UQ methods, this study provides actionable insights for improving the reliability, scalability, and utility of DL-based downscaling in climate science.

Environmental sciences↗
DOE OSTI2022

Representativeness assessment of the pan-Arctic eddy covariance site network and optimized future enhancements

Abstract. Large changes in the Arctic carbon balance are expected as warming linked to climate change threatens to destabilize ancient permafrost carbon stocks. The eddy covariance (EC) method is an established technique to quantify net losses and gains of carbon between the biosphere and atmosphere at high spatiotemporal resolution. Over the past decades, a growing network of terrestrial EC tower sites has been established across the Arctic, but a comprehensive assessment of the network's representativeness within the heterogeneous Arctic region is still lacking. This creates additional uncertainties when integrating flux data across sites, for example when upscaling fluxes to constrain pan-Arctic carbon budgets and changes therein. This study provides an inventory of Arctic (here > = 60∘ N) EC sites, which has also been made available online (https://cosima.nceas.ucsb.edu/carbon-flux-sites/, last access: 25 January 2022). Our database currently comprises 120 EC sites, but only 83 are listed as active, and just 25 of these active sites remain operational throughout the winter. To map the representativeness of this EC network, we evaluated the similarity between environmental conditions observed at the tower locations and those within the larger Arctic study domain based on 18 bioclimatic and edaphic variables. This allows us to assess a general level of similarity between ecosystem conditions within the domain, while not necessarily reflecting changes in greenhouse gas flux rates directly. We define two metrics based on this representativeness score: one that measures whether a location is represented by an EC tower with similar characteristics (ER1) and a second for which we assess if a minimum level of representation for statistically rigorous extrapolation is met (ER4). We find that while half of the domain is represented by at least one tower, only a third has enough towers in similar locations to allow reliable extrapolation. When we consider methane measurements or year-round (including wintertime) measurements, the values drop to about 1/5 and 1/10 of the domain, respectively. With the majority of sites located in Fennoscandia and Alaska, these regions were assigned the highest level of network representativeness, while large parts of Siberia and patches of Canada were classified as underrepresented. Across the Arctic, mountainous regions were particularly poorly represented by the current EC observation network. We tested three different strategies to identify new site locations or upgrades of existing sites that optimally enhance the representativeness of the current EC network. While 15 new sites can improve the representativeness of the pan-Arctic network by 20 %, upgrading as few as 10 existing sites to capture methane fluxes or remain active during wintertime can improve their respective ER1 network coverage by 28 % to 33 %. This targeted network improvement could be shown to be clearly superior to an unguided selection of new sites, therefore leading to substantial improvements in network coverage based on relatively small investments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Impacts of long-range-transported mineral dust on summertime convective cloud and precipitation: a case study over the Taiwan region

As one of the most abundant atmospheric aerosols and effective ice nuclei, mineral dust affects clouds and precipitation in the Earth system. Here numerical experiments are carried out to investigate the impacts of dust aerosols on summertime convective clouds and precipitation over the mountainous region of Taiwan by acting as ice-nucleating particles. We run the Weather Research and Forecasting model (WRF) with the Morrison two-moment and spectral-bin microphysics (SBM) schemes at 3 km resolution, using dust number concentrations from a global chemical transport model (GEOS-Chem-APM). The case study indicates that the long-range-transported mineral dust, with relatively low number concentrations, can notably affect the properties of convective clouds (ice and liquid water contents, cloud top height, and cloud coverage) and precipitation (spatial pattern and intensity). The effects of dust are evident during strong convective periods, with significantly increased ice water contents in the mixed-phase regime via the enhanced heterogeneous freezing. With both the Morrison and SBM schemes, we see the invigoration effects of dust aerosols on the convective intensity through enhanced condensation and deposition latent heating. The low-altitude dust particles are uplifted to the freezing level by updrafts, which, in turn, enhance the convective cloud development through immersion freezing and convective invigoration. Compared to the Morrison scheme, the SBM scheme predicts more realistic precipitation and different invigoration effects of dust. The differences are partially attributed to the saturation adjustment approach utilized in the bulk scheme, which leads to a stronger enhancement of condensation at midlatitudes to low altitudes and a weaker deposition increase at the upper level.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Evaluating 3 decades of precipitation in the Upper Colorado River basin from a high-resolution regional climate model

Abstract. Convection-permitting regional climate models (RCMs) have recently become tractable for applications at multi-decadal timescales. These types of models have tremendous utility for water resource studies, but better characterization of precipitation biases is needed, particularly for water-resource-critical mountain regions, where precipitation is highly variable in space, observations are sparse, and the societal water need is great. This study examines 34 years (1987–2020) of RCM precipitation from the Weather Research and Forecasting model (WRF; v3.8.1), using the Climate Forecast System Reanalysis (CFS; CFSv2) initial and lateral boundary conditions and a 1 km × 1 km innermost grid spacing. The RCM is centered over the Upper Colorado River basin, with a focus on the high-elevation, 750 km2 East River watershed (ERW), where a variety of high-impact scientific activities are currently ongoing. Precipitation is compared against point observations (Natural Resources Conservation Service Snow Telemetry or SNOTEL), gridded climate datasets (Newman, Livneh, and PRISM), and Bayesian reconstructions of watershed mean precipitation conditioned on streamflow and high-resolution snow remote-sensing products. We find that the cool-season precipitation percent error between WRF and 23 SNOTEL gauges has a low overall bias (x^ = 0.25 %, s = 13.63 %) and that WRF has a higher percent error during the warm season (x^ = 10.37 %, s = 12.79 %). Warm-season bias manifests as a high number of low-precipitation days, though the low-resolution or SNOTEL gauges limit some of the conclusions that can be drawn. Regional comparisons between WRF precipitation accumulation and three different gridded datasets show differences on the order of ± 20 %, particularly at the highest elevations and in keeping with findings from other studies. We find that WRF agrees slightly better with the Bayesian reconstruction of precipitation in the ERW compared to the gridded precipitation datasets, particularly when changing SNOTEL densities are taken into account. The conclusions are that the RCM reasonably captures orographic precipitation in this region and demonstrates that leveraging additional hydrologic information (streamflow and snow remote-sensing data) improves the ability to characterize biases in RCM precipitation fields. Error characteristics reported in this study are essential for leveraging the RCM model outputs for studies of past and future climates and water resource applications. The methods developed in this study can be applied to other watersheds and model configurations. Hourly 1 km × 1 km precipitation and other meteorological outputs from this dataset are publicly available and suitable for a wide variety of applications.

Geology↗