An official website of the United States government
Abstract Bayesian hierarchical models allow ecologists to account for uncertainty and make inference at multiple scales. However, hierarchical models are often computationally intensive to fit, especially with large datasets, and researchers face trade‐offs between capturing ecological complexity in statistical models and implementing these models.We present a recursive Bayesian computing (RB) method that can be used to fit Bayesian models efficiently in sequential MCMC stages to ease computation and streamline hierarchical inference. We also introduce transformation‐assisted RB (TARB) to create unsupervised MCMC algorithms and improve interpretability of parameters. We demonstrate TARB by fitting a hierarchical animal movement model to obtain inference about individual‐ and population‐level migratory characteristics.Our recursive procedure reduced computation time for fitting our hierarchical movement model by half compared to fitting the model with a single MCMC algorithm. We obtained the same inference fitting our model using TARB as we obtained fitting the model with a single algorithm.For complex ecological statistical models, like those for animal movement, multi‐species systems, or large spatial and temporal scales, the computational demands of fitting models with conventional computing techniques can limit model specification, thus hindering scientific discovery. Transformation‐assisted RB is one of the most accessible methods for reducing these limitations, enabling us to implement new statistical models and advance our understanding of complex ecological phenomena.
more »
« less
Homogenized gridded dataset for drought and hydrometeorological modeling for the continental United States
Abstract We present a novel data set for drought in the continental US (CONUS) built to enable computationally efficient spatio-temporal statistical and probabilistic models of drought. We converted drought data obtained from the widely-used US Drought Monitor (USDM) from its native geo-referenced polygon format to a 0.5 degree regular grid. We merged known environmental drivers of drought, including those obtained from the North American Land Data Assimilation System (NLDAS-2), US Geological Survey (USGS) streamflow data, and National Oceanic and Atmospheric Administration (NOAA) teleconnections data. The resulting data set permits statistical and probabilistic modeling of drought with explicit spatial and/or temporal dependence. Such models could be used to forecast drought at short-range, seasonal to sub-seasonal, and inter-annual timescales with uncertainty, extending the reach and value of the current US Drought Outlook from the National Weather Service Climate Prediction Center. This novel data product provides the first common gridded dataset that includes critical variables used to inform hydrological and meteorological drought.
more »
« less
- Award ID(s):
- 2151881
- PAR ID:
- 10500160
- Publisher / Repository:
- Nature Publishing Group
- Date Published:
- Journal Name:
- Scientific Data
- Volume:
- 11
- Issue:
- 1
- ISSN:
- 2052-4463
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Deterministic earthquake prediction remains elusive, but time‐dependent probabilistic seismicity forecasting seems within reach thanks to the development of physics‐based models relating seismicity to stress changes. Difficulties include constraining the earthquake nucleation model and fault initial stress state. Here, we analyze induced earthquakes from the Groningen gas field, where production is strongly seasonal, and seismicity began 3 decades after production started. We use the seismicity response to stress variations to constrain the earthquake nucleation process and calibrate models for time‐dependent forecasting of induced earthquakes. Remarkable agreements of modeled and observed seismicity are obtained when we consider (a) the initial strength excess, (b) the finite duration of earthquake nucleation, and (c) the seasonal variations of gas production. We propose a novel metric to quantify the nucleation model's ability to capture the damped amplitude and the phase of the seismicity response to short‐timescale (seasonal) stress variations which allows further tightening the model's parameters.more » « less
-
Abstract Drought poses a major threat to agricultural production and food security. This study evaluates the changes in drought-induced crop yield loss risk for six crops (alfalfa, barley, corn, soybean, spring wheat, and winter wheat) between 1971–2000 and 1991–2020 across the contiguous US. Using a copula-based probabilistic framework, our results reveal a spatially heterogeneous change in yield risk to meteorological droughts, which varies with crop types. Regional analyses identify the largest temporal decline in yield risk in the Southeast and Upper Midwest, while the Northwest and South show an increase in risk. Among the considered anthropogenic and climatic drivers of crop productivity, changes in climatic variables such as high temperatures (e.g., killing degree days), vapor pressure deficit and precipitation show significantly stronger associations with changes in yield risk than irrigated area and nitrogen fertilizer application. Among the counties that observe drier drought events, only 55% exhibit an increase in crop yield loss risk due to drier droughts. The rest 45% show a decrease in yield loss risk due to mediation of favorable climatic and anthropogenic factors. Alarmingly, more than half (for barley and spring wheat), and one-third (for alfalfa, corn, soybean and winter wheat) of that the risk increasing regions have outsized influence on destabilizing national crop production. The findings provide valuable insights for policymakers, agricultural stakeholders, and decision-makers in terms of the potential ways and locations to be prioritized for enhancing local and national agricultural resilience and ensuring food security.more » « less
-
Abstract The US Drought Monitor is the leading drought monitoring tool in the United States. Updated weekly and freely distributed, it records the drought conditions as geo-referenced polygons showing one of six ordered levels. These levels are determined by a mixture of quantitative environmental measurements and local expert opinion across the entire United States. At present, forecasts of the Drought Monitor only convey the expected direction of drought development (i.e. worsen, persist, subside) and do not communicate any uncertainty. This limits the utility of forecasts. In this paper, we describe a Bayesian spatio-temporal ordinal hierarchical model for use in modelling and projecting drought conditions. The model is flexible, scalable, and interpretable. By viewing drought data as areal rather than point-referenced, we reduce the cost of sampling from the posterior by avoiding dense matrix inversion. Draws from the posterior predictive distribution produce future forecasts of actual drought levels—rather than only the direction of drought development—and all sources of uncertainty are propagated into the posterior. Spatial random effects and an autoregressive model structure capture spatial and temporal dependence, and help ensure smoothness in forecasts over space and time. The result is a framework for modelling and forecasting drought levels and capturing forecast uncertainty.more » « less
-
Abstract As droughts have widespread social and ecological impacts, it is critical to develop long‐term adaptation and mitigation strategies to reduce drought vulnerability. Climate models are important in quantifying drought changes. Here, we assess the ability of 285 CMIP6 historical simulations, from 17 models, to reproduce drought duration and severity in three observational data sets using the Standardized Precipitation Index (SPI). We used summary statistics beyond the mean and standard deviation, and devised a novel probabilistic framework, based on the Hellinger distance, to quantify the difference between observed and simulated drought characteristics. Results show that many simulations have less thanerror in reproducing the observed drought summary statistics. The hypothesis that simulations and observations are described by the same distribution cannot be rejected for more thanof the grids based on ourdistance framework. No single model stood out as demonstrating consistently better performance over large regions of the globe. The variance in drought statistics among the simulations is higher in the tropics compared to other latitudinal zones. Though the models capture the characteristics of dry spells well, there is considerable bias in low precipitation values. Good model performance in terms of SPI does not imply good performance in simulating low precipitation. Our study emphasizes the need to probabilistically evaluate climate model simulations in order to both pinpoint model weaknesses and identify a subset of best‐performing models that are useful for impact assessments.more » « less
