Dynamical downscaling remains a powerful tool for studying regional climate processes, and the genesis of high‐resolution historical and future climate data. This technique is particularly important over areas of complex terrain, such as the western United States (WUS), where global models are especially limited in representing regional climate. After identifying a suite of WRF options that best simulate snow and precipitation for an average water year (2010) over the WUS, we evaluate the performance of the dynamically downscaled European Centre for Medium‐range Weather Forecasting's fifth Reanalysis (ERA5) from 1980 to 2020 on 45‐km, 9‐km, and two 3‐km grids. We find that by decreasing the horizontal grid spacing within WRF, improvements to Sierra Nevada and Northern Rocky Mountain snow, Santa Ana and Diablo winds, and coastal meteorology occur. For landfalling atmospheric rivers (ARs), the downscaled reanalysis simulates greater upstream integrated vapor transport (IVT) than ERA5. However, WRF skillfully simulates the positioning of the IVT and the timing and magnitude of AR precipitation. This potential IVT bias, in conjunction with increasing resolution, leads to a wet precipitation bias across the Sierra Nevada in the 3‐km experiment. This conclusion is supported by streamflow analysis, although we note that the bias in the 3‐km experiment can also be explained by in situ undercatch issues. Meanwhile, the 9‐km experiment is more biased than the 3‐km experiment across the Northern Rocky Mountains compared to in situ measured SWE and precipitation, indicating a geographic sensitivity to biases.

The Colorado River Basin is an important natural resource for the semi‐arid southwestern United States (US), where it provides water to more than 40 million people. While nearly 1.5°C of anthropogenic warming has occurred across this region from the 1880s to 2021, climate models show little agreement in the precipitation change during the same historical period, with no trend in the mean of the latest (sixth) generation of Global Climate Models. As such, here we focus on how the CO₂ increase and associated anthropogenic warming over the historical period has impacted runoff across the Colorado Basin. We find that the Colorado Basin's runoff over the historical period has decreased by 8.1% per degree Celsius of warming (°C⁻¹). However, the magnitude of this sensitivity is reduced to 6.8% °C⁻¹ when considering vegetation response to historical CO₂. For present‐day conditions, this translates to runoff reductions of 10.3% due to anthropogenic increases in both temperature and CO₂ since 1880. We demonstrate that Colorado Basin's natural flow has been decreased by roughly the storage of Lake Mead during the 2000–2021 megadrought due to this long term anthropogenic influence, suggesting the basin's first shortage in 2021 would likely not have occurred without anthropogenic warming. We further show warming has led to disproportionate aridification in snowpack regions, causing runoff to decline at double the rate relative to non‐snowpack regions. Thus, despite only making up ∼30% of the basin's drainage area, 86% of runoff decreases in the Colorado Basin is driven by water loss in snowpack regions.

Predicting future climate change over a region of complex terrain, such as the western United States (US), remains challenging due to the low resolution of global climate models (GCMs). Yet the climate extremes of recent years in this region, such as floods, wildfires, and drought, are likely to intensify further as climate warms, underscoring the need for high-quality and high-resolution predictions. Here, we present an ensemble of dynamically downscaled simulations over the western US from 1980–2100 at 9 km grid spacing, driven by 16 latest-generation GCMs. This dataset is titled the Western US Dynamically Downscaled Dataset (WUS-D3). We describe the challenges of producing WUS-D3, including GCM selection and technical issues, and we evaluate the simulations' realism by comparing historical results to temperature and precipitation observations. The future downscaled climate change signals are shaped in physically credible ways by the regional model's more realistic coastlines and topography. (1) The mean warming signals are heavily influenced by more realistic snowpack. (2) Mean precipitation changes are often consistent with wetting on the windward side of mountain complexes, as warmer, moister air masses are uplifted orographically during precipitation events. (3) There are large fractional precipitation increases on the lee side of mountain complexes, leading to potentially significant changes in water resources and ecology in these arid landscapes. (4) Increases in precipitation extremes are generally larger than in the GCMs, driven by locally intensified atmospheric updrafts tied to sharper, more realistic gradients in topography. (5) Changes in temperature extremes are different from what is expected by a shift in mean temperature and are shaped by local atmospheric dynamics and land surface feedbacks. Because of its high resolution, comprehensiveness, and representation of relevant physical processes, this dataset presents a unique opportunity to evaluate societally relevant future changes in western US climate.

In this study, we calibrate a regional climate model’s (RCM) underlying land surface model (LSM). In addition to providing a realistic representation of runoff across the hydroclimatically diverse western United States, this is done to take advantage of the RCM’s ability to physically resolve meteorological forcing data in ungauged regions, and to prepare the calibrated hydrologic model for tight coupling, or the ability to represent land surface–atmosphere interactions, with the RCM. Specifically, we use a 9-km resolution meteorological forcing dataset across the western United States, from the fifth generation ECMWF Reanalysis (ERA5) downscaled by the Weather Research Forecasting (WRF) regional climate model, as an offline forcing for Noah-Multiparameterization (Noah-MP). We detail the steps involved in producing an LSM capable of accurately representing runoff, including physical parameterization selection, parameter calibration, and regionalization to ungauged basins. Based on our model evaluation from 1954 to 2021 for 586 basins with daily natural streamflow, the streamflow bias is reduced from 24.2% to 4.4%, and the median daily Nash–Sutcliffe efficiency (NSE) is improved from 0.12 to 0.36. When validating against basins with monthly natural streamflow data, we obtain a similar reduction in bias and a median monthly NSE improvement from 0.18 to 0.56. In this study, we also discover the optimal setup when using a donor-basin method to regionalize parameters to ungauged basins, which can vary by 0.06 NSE for unique designs of this regionalization method. Significance Statement: This study provides useful guidance for improving a land surface model to accurately represent runoff across a spatially extensive and hydroclimatically diverse region (the western United States). The land surface model is updated to represent runoff more accurately at gauged basins, and then additionally updated for basins without observational data using a mathematical approach called the donor-basin method. We make use of a regional climate model’s reanalysis-derived meteorological data and its underlying land surface model to achieve realistic runoff. The calibrated land surface model can thus be tightly coupled in subsequent studies in a manner that should more accurately reflect runoff conditions. Findings from this study will serve as a useful reference for the atmospheric (and hydrologic) modeling communities and their ability to represent large-scale hydrology accurately.

Streamflow predictions are critical for managing water resources and for environmental conservation, especially in the water-short Western United States. Land surface models (LSMs), such as the variable infiltration capacity (VIC) model and the Noah LSM with multiparameterization options (Noah-MP), play an essential role in providing comprehensive runoff predictions across the region. Virtually all LSMs require parameter estimation (calibration) to optimize their predictive capabilities. Here, we focus on the calibration of VIC and Noah-MP models at a $\mathrm{1}/\mathrm{16}$ ° latitude–longitude resolution across the Western United States. We first performed global optimal calibration of parameters for both models for 263 river basins in the region. We find that the calibration significantly improves the models' performance, with the median daily streamflow Kling–Gupta efficiency (KGE) increasing from 0.37 to 0.70 for VIC, and from 0.22 to 0.54 for Noah-MP. In general, post-calibration model performance is higher for watersheds with relatively high precipitation and runoff ratios, and at lower elevations. At a second stage, we regionalize the river basin calibrations using the donor-basin method, which establishes transfer relationships for hydrologically similar basins, via which we extend our calibration parameters to 4816 hydrologic unit code (HUC)-10 basins across the region. Using the regionalized parameters, we show that the models' capabilities to simulate high and low flow conditions are substantially improved following calibration and regionalization. The refined parameter sets we developed are intended to support regional hydrological studies and hydrological assessments of climate change impacts.