skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


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
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 Observations of gravitational waves emitted by merging compact binaries have provided tantalizing hints about stellar astrophysics, cosmology, and fundamental physics. However, the physical parameters describing the systems (mass, spin, distance) used to extract these inferences about the Universe are subject to large uncertainties. The most widely used method of performing these analyses requires performing many Monte Carlo integrals to marginalize over the uncertainty in the properties of the individual binaries and the survey selection bias. These Monte Carlo integrals are subject to fundamental statistical uncertainties. Previous treatments of this statistical uncertainty have focused on ensuring that the precision of the inferred inference is unaffected; however, these works have neglected the question of whether sufficient accuracy can also be achieved. In this work, we provide a practical exploration of the impact of uncertainty in our analyses and provide a suggested framework for verifying that astrophysical inferences made with the gravitational-wave transient catalogue are accurate. Applying our framework to models used by the LIGO–Virgo–KAGRA collaboration and in the wider literature, we find that Monte Carlo uncertainty in estimating the survey selection bias is the limiting factor in our ability to probe narrow population models and this will rapidly grow more problematic as the size of the observed population increases. 
    more » « less
  3. 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
  4. Abstract Computational models of the cardiovascular system are increasingly used for the diagnosis, treatment, and prevention of cardiovascular disease. Before being used for translational applications, the predictive abilities of these models need to be thoroughly demonstrated through verification, validation, and uncertainty quantification. When results depend on multiple uncertain inputs, sensitivity analysis is typically the first step required to separate relevant from unimportant inputs, and is key to determine an initial reduction on the problem dimensionality that will significantly affect the cost of all downstream analysis tasks. For computationally expensive models with numerous uncertain inputs, sample‐based sensitivity analysis may become impractical due to the substantial number of model evaluations it typically necessitates. To overcome this limitation, we consider recently proposed Multifidelity Monte Carlo estimators for Sobol’ sensitivity indices, and demonstrate their applicability to an idealized model of the common carotid artery. Variance reduction is achieved combining a small number of three‐dimensional fluid–structure interaction simulations with affordable one‐ and zero‐dimensional reduced‐order models. These multifidelity Monte Carlo estimators are compared with traditional Monte Carlo and polynomial chaos expansion estimates. Specifically, we show consistent sensitivity ranks for both bi‐ (1D/0D) and tri‐fidelity (3D/1D/0D) estimators, and superior variance reduction compared to traditional single‐fidelity Monte Carlo estimators for the same computational budget. As the computational burden of Monte Carlo estimators for Sobol’ indices is significantly affected by the problem dimensionality, polynomial chaos expansion is found to have lower computational cost for idealized models with smooth stochastic response. 
    more » « less
  5. Long‐term trends in equatorial African rainfall have proven difficult to determine because of a dearth in ground‐measured rainfall data. Multiple, satellite‐based products now provide daily rainfall estimates from 1983 to the present at relatively fine spatial resolutions, but in order to assess trends in rainfall, they must be validated alongside ground‐based measurements. The purpose of this paper is twofold: (a) to assess the accuracy of four rainfall products covering the past several decades in western Uganda; and (b) to ascertain recent, multi‐decadal trends in annual rainfall for the region. The four products are African Rainfall Climatology Version 2 (ARC2), Climate Hazards Group InfraRed Precipitation with Stations (CHIRPS), Precipitation Estimation from Remotely Sensed Information using Artificial Neural Networks–Climate Data Record (PERSIANN‐CDR), and TAMSAT African Rainfall Climatology And Timeseries (TARCAT). The bias and accuracy of 10‐day, monthly, and seasonal rainfall totals of the four products were assessed using approximately 10 years of data from 10 rain gauges. The homogeneity of the products over multiple time periods was assessed using change‐point analysis. The accuracy of the four products increased with an increase in temporal scale, and CHIRPS was the only product that could be considered sufficiently accurate at estimating seasonal rainfall totals throughout most of the region. TARCAT tended to underestimate totals, and ARC2 and PERSIANN were in general the least accurate products. Only annual rainfall estimates from CHIRPS and TARCAT were significantly correlated with ground‐measured rainfall totals. TARCAT was the most homogeneous product, while ARC2, CHIRPS, and PERSIANN had significant negative change points that caused a drying bias over the 1983–2016 period. After adjusting the satellite‐based rainfall estimates based on the timing and magnitude of the change points, annual rainfall totals derived from all four products indicated that western Uganda experienced significantly increasing rainfall from 1983 to 2016. 
    more » « less