Search for: All records

Abstract Existing stochastic rainfall generators (SRGs) are typically limited to relatively small domains due to spatial stationarity assumptions, hindering their usefulness for flood studies in large basins. This study proposes StormLab, an SRG that simulates precipitation events at 6‐hr and 0.03° resolution in the Mississippi River Basin (MRB). The model focuses on winter and spring storms caused by water vapor transport from the Gulf of Mexico—the key flood‐generating storm type in the basin. The model generates anisotropic spatiotemporal noise fields that replicate local precipitation structures from observed data. The noise is transformed into precipitation through parametric distributions conditioned on large‐scale atmospheric fields from a climate model, reflecting spatial and temporal nonstationarity. StormLab can produce multiple realizations that reflect the uncertainty in fine‐scale precipitation arising from a specific large‐scale atmospheric environment. Model parameters were fitted monthly from December–May, based on storms identified from 1979 to 2021 ERA5 reanalysis data and Analysis of Record for Calibration (AORC) precipitation. StormLab then generated 1,000 synthetic years of precipitation events based on 10 CESM2 ensemble simulations. Empirical return levels of simulated annual maxima agree well with AORC data and show an overall increase in 1‐ to 500‐year events in the future period (2022–2050). To our knowledge, this is the first SRG simulating nonstationary, anisotropic high‐resolution precipitation over continental‐scale river basins, demonstrating the value of conditioning such stochastic models on large‐scale atmospheric variables. StormLab provides a wide range of extreme precipitation scenarios for design floods in the MRB and can be further extended to other large river basins. 

Abstract Storm direction modulates a hydrograph's magnitude and duration, thus having a potentially large effect on local flood risk. However, how changes in the preferential storm direction affect the probability distribution of peak flows remains unknown. We address this question with a novel Monte Carlo approach where stochastically transposed storms drive hydrologic simulations over medium and mesoscale watersheds in the Midwestern United States. Systematic rotations of these watersheds are used to emulate changes in the preferential storm direction. We found that the peak flow distribution impacts are scale‐dependent, with larger changes observed in the mesoscale watershed than in the medium‐scale watershed. We attribute this to the high diversity of storm patterns and the storms' scale relative to watershed size. This study highlights the potential of the proposed stochastic framework to address fundamental questions about hydrologic extremes when our ability to observe these events in nature is hindered by technical constraints and short time records. 

Abstract Estimating the probabilities of rare floods in mountainous watersheds is challenging due to the hydrometeorological complexity of seasonally varying snowmelt and soil moisture dynamics, as well as spatiotemporal variability in extreme precipitation. Design storm methods and statistical flood frequency analyses often overlook these complexities and how they shape the probabilities of rare floods. This study presents a process‐based approach that combines gridded precipitation, stochastic storm transposition (SST), and physics‐based distributed rainfall‐runoff modeling to simulate flood peak and volume distributions up to the 10,000‐year recurrence interval and to provide insights into the hydrometeorological drivers of those events. The approach is applied to a small mountainous watershed in the Colorado Front Range in the United States. We show that storm transposition in the Front Range can be justified under existing definitions of regional precipitation homogeneity. The process‐based results show close agreement with a statistically based mixture distribution that considers underlying flood drivers. We further demonstrate that antecedent conditions and snowmelt drive frequent peak discharges and rarer flood volumes, while the upper tail of the flood peak distribution appears to be controlled by heavy rainfall and rain‐on‐snow. In particular, we highlight the important role of early fall extreme rainfall in controlling rare flood peaks (but not volumes), despite only one such event having been observed in recent decades. Notwithstanding issues related to the accuracy of gridded precipitation datasets, these findings highlight the potential of SST and process‐based modeling to help understand the relationships between flood drivers and flood frequencies. 

Abstract Computational advances have made atmospheric modeling at convection‐permitting (≤4 km) grid spacings increasingly feasible. These simulations hold great promise in the projection of climate change impacts including rainfall and flood extremes. The relatively short model runs that are currently feasible, however, inhibit the assessment of the upper tail of rainfall and flood quantiles using conventional statistical methods. Stochastic storm transposition (SST) and process‐based flood frequency analysis are two approaches that together can help to mitigate this limitation. SST generates large numbers of extreme rainfall scenarios by temporal resampling and geospatial transposition of rainfall fields from relatively short data sets. Coupling SST with process‐based flood frequency analysis enables exploration of flood behavior at a range of spatial and temporal scales. We apply these approaches with outputs of 13‐year simulations of regional climate to examine changes in extreme rainfall and flood quantiles up to the 500‐year recurrence interval in a medium‐sized watershed in the Midwestern United States. Intensification of extreme precipitation across a range of spatial and temporal scales is identified in future climate; changes in flood magnitudes depend on watershed area, with small watersheds exhibiting the greatest increases due to their limited capacity to attenuate flood peaks. Flood seasonality and snowmelt are predicted to be earlier in the year under projected warming, while the most extreme floods continue to occur in early summer. Findings highlight both the potential and limitations of convection‐resolving climate models to help understand possible changes in rainfall and flood frequency across watershed scales. 

Abstract Although prior studies have evaluated the role of sampling errors associated with local and regional methods to estimate peak flow quantiles, the investigation of epistemic errors is more difficult because the underlying properties of the random variable have been prescribed using ad‐hoc characterizations of the regional distributions of peak flows. This study addresses this challenge using representations of regional peak flow distributions derived from a combined framework of stochastic storm transposition, radar rainfall observations, and distributed hydrologic modeling. The authors evaluated four commonly used peak flow quantile estimation methods using synthetic peak flows at 5,000 sites in the Turkey River watershed in Iowa, USA. They first used at‐site flood frequency analysis using the Pearson Type III distribution with L‐moments. The authors then pooled regional information using (1) the index flood method, (2) the quantile regression technique, and (3) the parameter regression. This approach allowed quantification of error components stemming from epistemic assumptions, parameter estimation method, sample size, and, in the regional approaches, the number ofpooledsites. The results demonstrate that the inability to capture the spatial variability of the skewness of the peak flows dominates epistemic error for regional methods. We concluded that, in the study basin, this variability could be partially explained by river network structure and the predominant orientation of the watershed. The general approach used in this study is promising in that it brings new tools and sources of data to the study of the old hydrologic problem of flood frequency analysis. 

Abstract. Conventional rainfall frequency analysis faces several limitations. These include difficulty incorporating relevant atmospheric variables beyond precipitation and limited ability to depict the frequency of rainfall over large areas that is relevant for flooding. This study proposes a storm-based model of extreme precipitation frequency based on the atmospheric water balance equation. We developed a storm tracking and regional characterization (STARCH) method to identify precipitation systems in space and time from hourly ERA5 precipitation fields over the contiguous United States from 1951 to 2020. Extreme “storm catalogs” were created by selecting annual maximum storms with specific areas and durations over a chosen region. The annual maximum storm precipitation was then modeled via multivariate distributions of atmospheric water balance components using vine copula models. We applied this approach to estimate precipitation average recurrence intervals for storm areas from 5000 to 100 000 km2 and durations from 2 to 72 h in the Mississippi Basin and its five major subbasins. The estimated precipitation distributions show a good fit to the reference data from the original storm catalogs and are close to the estimates from conventional univariate GEV distributions. Our approach explicitly represents the contributions of water balance components in extreme precipitation. Of these, water vapor flux convergence is the main contributor, while precipitable water and a mass residual term can also be important, particularly for short durations and small storm footprints. We also found that ERA5 shows relatively good water balance closure for extreme storms, with a mass residual on average 10 % of precipitation. The approach can incorporate nonstationarities in water balance components and their dependence structures and can benefit from further advancements in reanalysis products and storm tracking techniques.