skip to main content


Title: A Nonstationary Standardized Precipitation Index (NSPI) Using Bayesian Splines
Abstract

The standardized precipitation index (SPI) measures meteorological drought relative to historical climatology by normalizing accumulated precipitation. Longer record lengths improve parameter estimates, but these longer records may include signals of anthropogenic climate change and multidecadal natural climate fluctuations. Historically, climate nonstationarity has either been ignored or incorporated into the SPI using a quasi-stationary reference period, such as the WMO 30-yr period. This study introduces and evaluates a novel nonstationary SPI model based on Bayesian splines, designed to both improve parameter estimates for stationary climates and to explicitly incorporate nonstationarity. Using synthetically generated precipitation, this study directly compares the proposed Bayesian SPI model with existing SPI approaches based on maximum likelihood estimation for stationary and nonstationary climates. The proposed model not only reproduced the performance of existing SPI models but improved upon them in several key areas: reducing parameter uncertainty and noise, simultaneously modeling the likelihood of zero and positive precipitation, and capturing nonlinear trends and seasonal shifts across all parameters. Further, the fully Bayesian approach ensures all parameters have uncertainty estimates, including zero precipitation likelihood. The study notes that the zero precipitation parameter is too sensitive and could be improved in future iterations. The study concludes with an application of the proposed Bayesian nonstationary SPI model for nine gauges across a range of hydroclimate zones in the United States. Results of this experiment show that the model is stable and reproduces nonstationary patterns identified in prior studies, while also indicating new findings, particularly for the shape and zero precipitation parameters.

Significance Statement

We typically measure how bad a drought is by comparing it with the historical record. With long-term changes in climate or other factors, however, a typical drought today may not have been typical in the recent past. The purpose of this study is to build a model that measures drought relative to a changing climate. Our results confirm that the model is accurate and captures previously noted climate change patterns—a drier western United States, a wetter eastern United States, earlier summer weather, and more extreme wet seasons. This is significant because this model can improve drought measurement and identify recent changes in drought.

 
more » « less
Award ID(s):
2002539
NSF-PAR ID:
10374418
Author(s) / Creator(s):
 ;  
Publisher / Repository:
American Meteorological Society
Date Published:
Journal Name:
Journal of Applied Meteorology and Climatology
Volume:
61
Issue:
7
ISSN:
1558-8424
Page Range / eLocation ID:
p. 761-779
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Analyzing gradual trends in meteorological drought has become increasingly important as anthropogenic climate change and natural climate variability interact to complicate measurement of drought severity. Complex seasonality and long-term trends pose a limitation in understanding spatial trends in nonstationary changes of meteorological drought in the United States. This study seeks to address this issue by simultaneously analyzing recurring seasonal patterns (stationary component) and long-term drought trends (nonstationary component), with a unique focus on nonlinear trends and common regional patterns. We analyzed 696 instrumental precipitation gauges with long historical records in the continental United States, using a novel spline-based model to disaggregate a 3-month meteorological drought index (SPI) into its seasonal and long-term components. The disaggregated components for each gauge were then clustered into subregions with similar seasonality and groupings with similar long-term trends using a two-step process. Our results identify clearly defined regions based on precipitation seasonality, while long-term trends are not spatially coherent with the seasonality. Instead, these findings support prior findings of an increasingly drier western United States and an increasingly wetter eastern United States over the last century, but with more nuanced spatial and temporal patterns. The new clustering analysis based on nonstationary meteorological drought trends can contribute to informing and adapting current water management strategies to long-term drought trends. Significance Statement This study considered 656 precipitation gauges across the continental United States to find regions with similar precipitation seasonality and then to group records with similar long-term climate trends. The study focused on 3-month average precipitation, a key indicator for drought monitoring. We identified eight regions across the United States with similar precipitation seasonality. From 1920 to the present, we found continuous drying trends throughout the western United States, continuously wetter trends in the northern plains, and an overall wetter trend interrupted by a midcentury dry period (1930–50) for much of the central Plains and Midwest. This study’s use of splines, or fitted curves, allowed these nonlinear patterns, which we believe better capture the nuances and intensification of climate change effects on precipitation. 
    more » « less
  2. Abstract

    Flood‐frequency curves, critical for water infrastructure design, are typically developed based on a stationary climate assumption. However, climate changes are expected to violate this assumption. Here, we propose a new, climate‐informed methodology for estimating flood‐frequency curves under non‐stationary future climate conditions. The methodology develops an asynchronous, semiparametric local‐likelihood regression (ASLLR) model that relates moments of annual maximum flood to climate variables using the generalized linear model. We estimate the first two marginal moments (MM) – the mean and variance – of the underlying log‐Pearson Type‐3 distribution from the ASLLR with the monthly rainfall and temperature as predictors. The proposed methodology, ASLLR‐MM, is applied to 40 U.S. Geological Survey streamgages covering 18 water resources regions across the conterminous United States. A correction based on the aridity index was applied on the estimated variance, after which the ASLLR‐MM approach was evaluated with both historical (1951–2005) and projected (2006–2035, under RCP4.5 and RCP8.5) monthly precipitation and temperature from eight Global Circulation Models (GCMs) consisting of 39 ensemble members. The estimated flood‐frequency quantiles resulting from the ASLLR‐MM and GCM members compare well with the flood‐frequency quantiles estimated using the historical period of observed climate and flood information for humid basins, whereas the uncertainty in model estimates is higher in arid basins. Considering additional atmospheric and land‐surface conditions and a multi‐level model structure that includes other basins in a region could further improve the model performance in arid basins.

     
    more » « less
  3. Abstract

    Satellite precipitation products, as all quantitative estimates, come with some inherent degree of uncertainty. To associate a quantitative value of the uncertainty to each individual estimate, error modeling is necessary. Most of the error models proposed so far compute the uncertainty as a function of precipitation intensity only, and only at one specific spatiotemporal scale. We propose a spectral error model that accounts for the neighboring space–time dynamics of precipitation into the uncertainty quantification. Systematic distortions of the precipitation signal and random errors are characterized distinctively in every frequency–wavenumber band in the Fourier domain, to accurately characterize error across scales. The systematic distortions are represented as a deterministic space–time linear filtering term. The random errors are represented as a nonstationary additive noise. The spectral error model is applied to the IMERG multisatellite precipitation product, and its parameters are estimated empirically through a system identification approach using the GV-MRMS gauge–radar measurements as reference (“truth”) over the eastern United States. The filtering term is found to be essentially low-pass (attenuating the fine-scale variability). While traditional error models attribute most of the error variance to random errors, it is found here that the systematic filtering term explains 48% of the error variance at the native resolution of IMERG. This fact confirms that, at high resolution, filtering effects in satellite precipitation products cannot be ignored, and that the error cannot be represented as a purely random additive or multiplicative term. An important consequence is that precipitation estimates derived from different sources shall not be expected to automatically have statistically independent errors.

    Significance Statement

    Satellite precipitation products are nowadays widely used for climate and environmental research, water management, risk analysis, and decision support at the local, regional, and global scales. For all these applications, knowledge about the accuracy of the products is critical for their usability. However, products are not systematically provided with a quantitative measure of the uncertainty associated with each individual estimate. Various parametric error models have been proposed for uncertainty quantification, mostly assuming that the uncertainty is only a function of the precipitation intensity at the pixel and time of interest. By projecting satellite precipitation fields and their retrieval errors into the Fourier frequency–wavenumber domain, we show that we can explicitly take into account the neighboring space–time multiscale dynamics of precipitation and compute a scale-dependent uncertainty.

     
    more » « less
  4. Abstract

    The selection of calibration and validation time periods in hydrologic modelling is often done arbitrarily. Nonstationarity can lead to an optimal parameter set for one period which may not accurately simulate another. However, there is still much to be learned about the responses of hydrologic models to nonstationary conditions. We investigated how the selection of calibration and validation periods can influence water balance simulations. We calibrated Soil and Water Assessment Tool hydrologic models with observed streamflow for three United States watersheds (St. Joseph River of Indiana/Michigan, Escambia River of Florida/Alabama, and Cottonwood Creek of California), using time period splits for calibration/validation. We found that the choice of calibration period (with different patterns of observed streamflow, precipitation, and air temperature) influenced the parameter sets, leading to dissimilar simulations of water balance components. In the Cottonwood Creek watershed, simulations of 50‐year mean January streamflow varied by 32%, because of lower winter precipitation and air temperature in earlier calibration periods on calibrated parameters, which impaired the ability for models calibrated to earlier periods to simulate later periods. Peaks of actual evapotranspiration for this watershed also shifted from April to May due to different parameter values depending on the calibration period's winter air temperatures. In the St. Joseph and Escambia River watersheds, adjustments of the runoff curve number parameter could vary by 10.7% and 20.8%, respectively, while 50‐year mean monthly surface runoff simulations could vary by 23%–37% and 169%–209%, depending on the observed streamflow and precipitation of the chosen calibration period. It is imperative that calibration and validation time periods are chosen selectively instead of arbitrarily, for instance using change point detection methods, and that the calibration periods are appropriate for the goals of the study, considering possible broad effects of nonstationary time series on water balance simulations. It is also crucial that the hydrologic modelling community improves existing calibration and validation practices to better include nonstationary processes.

     
    more » « less
  5. Abstract

    Soil moisture is an important driver of growth in boreal Alaska, but estimating soil hydraulic parameters can be challenging in this data‐sparse region. Parameter estimation is further complicated in regions with rapidly warming climate, where there is a need to minimize model error dependence on interannual climate variations. To better identify soil hydraulic parameters and quantify energy and water balance and soil moisture dynamics, we applied the physically based, one‐dimensional ecohydrological Simultaneous Heat and Water (SHAW) model, loosely coupled with the Geophysical Institute of Permafrost Laboratory (GIPL) model, to an upland deciduous forest stand in interior Alaska over a 13‐year period. Using a Generalized Likelihood Uncertainty Estimation parameterisation, SHAW reproduced interannual and vertical spatial variability of soil moisture during a five‐year validation period quite well, with root mean squared error (RMSE) of volumetric water content at 0.5 m as low as 0.020 cm3/cm3. Many parameter sets reproduced reasonable soil moisture dynamics, suggesting considerable equifinality. Model performance generally declined in the eight‐year validation period, indicating some overfitting and demonstrating the importance of interannual variability in model evaluation. We compared the performance of parameter sets selected based on traditional performance measures such as the RMSE that minimize error in soil moisture simulation, with one that is designed to minimize the dependence of model error on interannual climate variability using a new diagnostic approach we call CSMP, which stands for Climate Sensitivity of Model Performance. Use of the CSMP approach moderately decreases traditional model performance but may be more suitable for climate change applications, for which it is important that model error is independent from climate variability. These findings illustrate (1) that the SHAW model, coupled with GIPL, can adequately simulate soil moisture dynamics in this boreal deciduous region, (2) the importance of interannual variability in model parameterisation, and (3) a novel objective function for parameter selection to improve applicability in non‐stationary climates.

     
    more » « less