Search

Western U.S. (WUS) rainfall and snowpack vary greatly on interannual and decadal timescales. This combined with their importance to water resources makes future projections of these variables highly societally relevant. Previous studies have shown that precipitation events in the WUS are influenced by the timing, positioning, and duration of extreme integrated water vapor transport (IVT) events (e.g., atmospheric rivers) along the coast. We investigate end-of-21st-century projections of WUS precipitation and IVT in a collection of regional climate models (RCMs) from the North American Coordinated Regional Downscaling Experiment (NA-CORDEX). Several of the NA-CORDEX RCMs project a decrease in cool season precipitation at high elevation (e.g., across the Sierra Nevada) with a corresponding increase in the Great Basin of the U.S. We explore the causes of this terrain-related precipitation change in a subset of the NA-CORDEX RCMs through an examination of IVT-events. Projected changes in frequency and duration of IVT-events depend on the event's extremity: By the end of the century extreme IVT-events increase in frequency whereas moderate IVT-events decrease in frequency. Furthermore, in the future, total precipitation across the WUS generally increases during extreme IVT-events, whereas total precipitation from moderate IVT-events decreases across higher elevations. Thus, we argue that the mean cool season precipitation decreases at high elevations and increases in the Great Basin are largely determined by changes in moderate IVT-events which are projected to be less frequent and bring less high-elevation precipitation.

Accurate estimation of the potential “upper limit” for extreme precipitation is critical for dam safety and water resources management, as dam failures pose significant risks to life and property. Methods used to estimate the theoretical upper limit of precipitation are often outdated and in need of updating. The rarity of extreme events means that old storms with limited observational data are often used to define the upper bound of precipitation. Observations of many important old storms are limited in spatial and temporal coverage, and sometimes of dubious quality. This reduces confidence in flood hazard assessments used in dam safety evaluations and leads to unknown or uncertain societal risk. This paper describes a method for generating and applying ensembles of high-resolution, state-of-the-art numerical model simulations of historical past extreme precipitation events to meet contemporary stakeholder needs. The method was designed as part of a research-to-application-focused partnership project to update state dam safety rules in Colorado and New Mexico. The results demonstrated multiple stakeholder and user benefits that were applied directly into storm analyses utilized for extreme rainfall estimation, and diagnostics were developed and ultimately used to update Colorado state dam safety rules, officially passed in January 2020. We discuss how what started as a prototype research foray to meet a specific user need may ultimately inform wider adoption of numerical simulations for water resources risk assessment, and how the historical event downscaling method performed offers near-term, implementable improvements to current dam safety flood risk estimates that can better serve society today.

Accurate quantitative precipitation estimates (QPEs) at high spatial and temporal resolution are difficult to obtain in regions of complex terrain due to the large spatial heterogeneity of orographically enhanced precipitation, sparsity of gauges, precipitation phase variations, and terrain effects that impact the quality of remotely sensed estimates. The large uncertainty of QPE in these regions also makes the evaluation of high-resolution quantitative precipitation forecasts (QPFs) challenging, as it can be difficult to choose a reference QPE that is reliable at both high and low elevations. In this paper we demonstrate a methodology to combine information from multiple high-resolution hourly QPE products to evaluate QPFs from NOAA’s High-Resolution Rapid Refresh (HRRR) model in a region of Northern California. The methodology uses the quantiles of monthly QPE distributions to determine a range of hourly precipitation that correspond to “good,” “possible,” “underestimated,” or “overestimated” QPFs. In this manuscript, we illustrate the use of the methodology to evaluate QPFs for seven atmospheric river events that occurred during the 2016–17 wet season in Northern California. Because the presence of frozen precipitation is often not captured by traditional QPE products, we evaluate QPFs both for all precipitation, and with likely frozen precipitation excluded. The methodology is shown to provide useful information to evaluate model performance while taking into account the uncertainty of available QPE at various temporal and spatial scales. The potential of the technique to evaluate changes between model versions is also shown.

Snow influences land–atmosphere interactions in snow-dominated areas, and snow melt contributes to basin streamflows. However, estimating snowpack properties such as the snow depth (SD) and snow water equivalent (SWE) from land surface model simulations remains a challenge. This is, in part, due to uncertainties in the atmospheric forcing variables, which propagate into hydrological model predictions. This study implements the Weather Research and Forecasting (WRF)-Hydro framework with the Noah-Multiparameterization (Noah-MP) land surface model in the NOAA’s National Water Model (NWM) version 2.0 configuration to estimate snow in a single column and subsequently the streamflow across the Aroostook River’s sub-basins in Maine for water years (WY) 2014–2016. This study evaluates how differences between two atmospheric forcing datasets, the North American Land Data Assimilation version 2 (NLDAS-2) and in situ (Station), translate into differences in the simulation of snow. NLDAS-2 was used as the meteorological forcing in the retrospective NWM 2.0 simulations. The results from the single-column study showed that differences in the simulated SWE and SD were linked to differences in the 2 m air temperature (T2m), which influenced the precipitation partitioning of rain and snow, as parameterized in Noah-MP. The negative mean bias of −0.7 K (during the accumulation period) in T2m for NLDAS-2, compared to the Station forcing, was a major factor that contributed to the positive mean bias of +52 mm on average in the peak SWE in the NLDAS-2-forced Noah-MP simulation during the study period. The higher T2m values at the Station led to higher sensible heat fluxes towards the snowpack, which led to a higher amount of net energy at the snow’s surface and melt events during the accumulation season in Station-forced Noah-MP simulations. The results from the retrospective NWM version 2.0′s simulation in the basin showed that the streamflow estimates were closer to the United States Geological Survey gage observations at the two larger sub-basins (NSE = 0.9), which were mostly forested, compared to the two smaller sub-basins (NSE ≥ 0.4), which had more agricultural land-use. This study also showed that the spring snowmelt timing was captured quite well by the timing of the decline in the simulated SWE and SD, providing an early indication of melt in most sub-basins. The simulated fractional snow cover area (fSCA) however provided less information about the changes in snow or onset of snowmelt as it was mostly binary (full snow cover in winter), which differed from the more realistic fSCA values shown by the Moderate Resolution Imaging Spectroradiometer.