skip to main content


Title: Comparison of Local, Regional, and Scaling Models for Rainfall Intensity–Duration–Frequency Analysis
Abstract Intensity–duration–frequency (IDF) analyses of rainfall extremes provide critical information to mitigate, manage, and adapt to urban flooding. The accuracy and uncertainty of IDF analyses depend on the availability of historical rainfall records, which are more accessible at daily resolution and, quite often, are very sparse in developing countries. In this work, we quantify performances of different IDF models as a function of the number of available high-resolution (Nτ) and daily (N24h) rain gauges. For this aim, we apply a cross-validation framework that is based on Monte Carlo bootstrapping experiments on records of 223 high-resolution gauges in central Arizona. We test five IDF models based on (two) local, (one) regional, and (two) scaling frequency analyses of annual rainfall maxima from 30-min to 24-h durations with the generalized extreme value (GEV) distribution. All models exhibit similar performances in simulating observed quantiles associated with return periods up to 30 years. When Nτ > 10, local and regional models have the best accuracy; bias correcting the GEV shape parameter for record length is recommended to estimate quantiles for large return periods. The uncertainty of all models, evaluated via Monte Carlo experiments, is very large when Nτ ≤ 5; however, if N24h ≥ 10 additional daily gauges are available, the uncertainty is greatly reduced and accuracy is increased by applying simple scaling models, which infer estimates on subdaily rainfall statistics from information at daily scale. For all models, performances depend on the ability to capture the elevation control on their parameters. Although our work is site specific, its results provide insights to conduct future IDF analyses, especially in regions with sparse data.  more » « less
Award ID(s):
1831475 1928724
NSF-PAR ID:
10194664
Author(s) / Creator(s):
Date Published:
Journal Name:
Journal of Applied Meteorology and Climatology
Volume:
59
Issue:
9
ISSN:
1558-8424
Page Range / eLocation ID:
1519 to 1536
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Urgency of Precipitation Intensity-Duration-Frequency (IDF) estimation using the most recent data has grown significantly due to recent intense precipitation and cloud burst circumstances impacting infrastructure caused by climate change. Given the continually available digitized up-to-date, long-term, and fine resolution precipitation dataset from the United States Department of Agriculture Forest Service’s (USDAFS) Experimental Forests and Ranges (EF) rain gauge stations, it is both important and relevant to develop precipitation IDF from onsite dataset (Onsite-IDF) that incorporates the most recent time period, aiding in the design, and planning of forest road-stream crossing structures (RSCS) in headwaters to maintain resilient forest ecosystems. Here we developed Onsite-IDFs for hourly and sub-hourly duration, and 25-yr, 50-yr, and 100-yr design return intervals (RIs) from annual maxima series (AMS) of precipitation intensities (PIs) modeled by applying Generalized Extreme Value (GEV) analysis and L-moment based parameter estimation methodology at six USDAFS EFs and compared them with precipitation IDFs obtained from the National Oceanic and Atmospheric Administration Atlas 14 (NOAA-Atlas14). A regional frequency analysis (RFA) was performed for EFs where data from multiple precipitation gauges are available. NOAA’s station-based precipitation IDFs were estimated for comparison using RFA (NOAA-RFA) at one of the EFs where NOAA-Atlas14 precipitation IDFs are unavailable. Onsite-IDFs were then evaluated against the PIs from NOAA-Atlas14 and NOAA-RFA by comparing their relative differences and storm frequencies. Results show considerable relative differences between the Onsite- and NOAA-Atlas14 (or NOAA-RFA) IDFs at these EFs, some of which are strongly dependent on the storm durations and elevation of precipitation gauges, particularly in steep, forested sites of H. J. Andrews (HJA) and Coweeta Hydrological Laboratory (CHL) EFs. At the higher elevation gauge of HJA EF, NOAA-RFA based precipitation IDFs underestimate PI of 25-yr, 50-yr, and 100-yr RIs by considerable amounts for 12-h and 24-h duration storm events relative to the Onsite-IDFs. At the low-gradient Santee (SAN) EF, the PIs of 3- to 24-h storm events with 100-yr frequency (or RI) from NOAA-Atlas14 gauges are found to be equivalent to PIs of more frequent storm events (25–50-yr RI) as estimated from the onsite dataset. Our results recommend use of the Onsite-IDF estimates for the estimation of design storm peak discharge rates at the higher elevation catchments of HJA, CHL, and SAN EF locations, particularly for longer duration events, where NOAA-based precipitation IDFs underestimate the PIs relative to the Onsite-IDFs. This underscores the importance of long-term high resolution EF data for new applications including ecological restorations and indicates that planning and design teams should use as much local data as possible or account for potential PI inconsistencies or underestimations if local data are unavailable.

     
    more » « less
  2. Abstract

    The estimation of the frequency of intense rainfall events is a crucial step for quantifying their impact on human societies and on the environment. This process is hindered by large gaps in ground observational networks at the global scale, such that extensive areas remain ungauged. The increasing availability of satellite‐based rainfall estimates, while providing data with unprecedented resolution and global coverage, also introduces new challenges: the scale disparity between gridded and rain‐gauge precipitation products on the one hand, and the short length of the available satellite records on the other. Here we propose a statistical framework for the estimation of rainfall extremes that is specifically designed to simultaneously address these two key issues, providing a new way of estimating extreme rainfall magnitudes from space. A downscaling procedure is here introduced to recover the spatial correlation and the probability density function of daily rainfall at the point (gauge) scale from coarse‐scale satellite estimates. The results are then combined with a recent statistical model of extremes (the Metastatistical Extreme Value distribution), which optimizes the use of the information obtained from relatively short satellite observational time series. The methodology is tested using data from the Tropical Rainfall Measuring Mission Multisatellite Precipitation Analysis over the Little Washita River, Oklahoma. We find that our approach satisfactorily reproduces downscaled daily rainfall probability density functions and can significantly improve the Tropical Rainfall Measuring Mission Multisatellite Precipitation Analysis‐based estimation of quantiles with return times larger than the length of the available data set (19 years here), which are especially important for several water‐related applications.

     
    more » « less
  3. Extreme storm surges can overwhelm many coastal flooding protection measures in place and cause severe damages to private communities, public infrastructure, and natural ecosystems. In the US Mid-Atlantic, a highly developed and commercially active region, coastal flooding is one of the most significant natural hazards and a year-round threat from both tropical and extra-tropical cyclones. Mean sea levels and high-tide flood frequency has increased significantly in recent years, and major storms are projected to increase into the foreseeable future. We estimate extreme surges using hourly water level data and harmonic analysis for 1980–2019 at 12 NOAA tide gauges in and around the Delaware and Chesapeake Bays. Return levels (RLs) are computed for 1.1, 3, 5, 10, 25, 50, and 100-year return periods using stationary extreme value analysis on detrended skew surges. Two traditional approaches are investigated, Block Maxima fit to General Extreme Value distribution and Points-Over-Threshold fit to Generalized Pareto distribution, although with two important enhancements. First, the GEV r -largest order statistics distribution is used; a modified version of the GEV distribution that allows for multiple maximum values per year. Second, a systematic procedure is used to select the optimum value for r (for the BM/GEVr approach) and the threshold (for the POT/GP approach) at each tide gauge separately. RLs have similar magnitudes and spatial patterns from both methods, with BM/GEVr resulting in generally larger 100-year and smaller 1.1-year RLs. Maximum values are found at the Lewes (Delaware Bay) and Sewells Point (Chesapeake Bay) tide gauges, both located in the southwest region of their respective bays. Minimum values are found toward the central bay regions. In the Delaware Bay, the POT/GP approach is consistent and results in narrower uncertainty bands whereas the results are mixed for the Chesapeake. Results from this study aim to increase reliability of projections of extreme water levels due to extreme storms and ultimately help in long-term planning of mitigation and implementation of adaptation measures. 
    more » « less
  4. Abstract

    This study proposed a framework to evaluate multivariate return periods of hurricanes using event‐based frequency analysis techniques. The applicability of the proposed framework was demonstrated using point‐based and spatial analyses on a recent catastrophic event, Hurricane Ian. Univariate, bivariate, and trivariate frequency analyses were performed by applying generalized extreme value distribution and copula on annual maximum series of flood volume, peak discharge, total rainfall depth, maximum wind speed, wave height and storm surge. As a result of point‐based analyses, return periods of Hurricane Ian was investigated by using our framework; univariate return periods were estimated from 39.2 to 60.2 years, bivariate from 824.1 to 1,592.6 years, and trivariate from 332.1 to 1,722.9 years for the Daytona‐St. Augustine Basin. In the Florida Bay‐Florida Keys Basin, univariate return periods were calculated from 7.5 to 32.9 years, bivariate from 36.5 to 114.9 years, and trivariate from 25.0 to 214.8 years. Using the spatial analyses, we were able to generate the return period map of Hurricane Ian across Florida. Based on bivariate frequency analyses, 18% of hydrologic unit code 8 (HUC8) basins had an average return period of more than 30 years. Sources of uncertainty, due to the scarcity of analysis data, stationarity assumption and impact of other weather systems such as strong frontal passages, were also discussed. Despite these limitations, our framework and results will be valuable in disaster response and recovery.

     
    more » « less
  5. Abstract We address the challenge, due to sparse observational records, of investigating long-term changes in the storm surge climate globally. We use two centennial and three satellite-era daily storm surge time series from the Global Storm Surge Reconstructions (GSSR) database and assess trends in the magnitude and frequency of extreme storm surge events at 320 tide gauges across the globe from 1930, 1950, and 1980 to present. Before calculating trends, we perform change point analysis to identify and remove data where inhomogeneities in atmospheric reanalysis products could lead to spurious trends in the storm surge data. Even after removing unreliable data, the database still extends existing storm surge records by several decades for most of the tide gauges. Storm surges derived from the centennial 20CR and ERA-20C atmospheric reanalyses show consistently significant positive trends along the southern North Sea and the Kattegat Bay regions during the periods from 1930 and 1950 onwards and negative trends since 1980 period. When comparing all five storm surge reconstructions and observations for the overlapping 1980–2010 period we find overall good agreement, but distinct differences along some coastlines, such as the Bay of Biscay and Australia. We also assess changes in the frequency of extreme surges and find that the number of annual exceedances above the 95th percentile has increased since 1930 and 1950 in several regions such as Western Europe, Kattegat Bay, and the US East Coast. 
    more » « less