scholarly journals A space–time Bayesian hierarchical modeling framework for projection of seasonal maximum streamflow

2022 ◽  
Vol 26 (1) ◽  
pp. 149-166
Author(s):  
Álvaro Ossandón ◽  
Manuela I. Brunner ◽  
Balaji Rajagopalan ◽  
William Kleiber

Abstract. Timely projections of seasonal streamflow extremes can be useful for the early implementation of annual flood risk adaptation strategies. However, predicting seasonal extremes is challenging, particularly under nonstationary conditions and if extremes are correlated in space. The goal of this study is to implement a space–time model for the projection of seasonal streamflow extremes that considers the nonstationarity (interannual variability) and spatiotemporal dependence of high flows. We develop a space–time model to project seasonal streamflow extremes for several lead times up to 2 months, using a Bayesian hierarchical modeling (BHM) framework. This model is based on the assumption that streamflow extremes (3 d maxima) at a set of gauge locations are realizations of a Gaussian elliptical copula and generalized extreme value (GEV) margins with nonstationary parameters. These parameters are modeled as a linear function of suitable covariates describing the previous season selected using the deviance information criterion (DIC). Finally, the copula is used to generate streamflow ensembles, which capture spatiotemporal variability and uncertainty. We apply this modeling framework to predict 3 d maximum streamflow in spring (May–June) at seven gauges in the Upper Colorado River basin (UCRB) with 0- to 2-month lead time. In this basin, almost all extremes that cause severe flooding occur in spring as a result of snowmelt and precipitation. Therefore, we use regional mean snow water equivalent and temperature from the preceding winter season as well as indices of large-scale climate teleconnections – El Niño–Southern Oscillation, Atlantic Multidecadal Oscillation, and Pacific Decadal Oscillation – as potential covariates for 3 d spring maximum streamflow. Our model evaluation, which is based on the comparison of different model versions and the energy skill score, indicates that the model can capture the space–time variability in extreme streamflow well and that model skill increases with decreasing lead time. We also find that the use of climate variables slightly enhances skill relative to using only snow information. Median projections and their uncertainties are consistent with observations, thanks to the representation of spatial dependencies through covariates in the margins and a Gaussian copula. This spatiotemporal modeling framework helps in the planning of seasonal adaptation and preparedness measures as predictions of extreme spring streamflows become available 2 months before actual flood occurrence.

2021 ◽  
Author(s):  
Álvaro Ossandón ◽  
Manuela I Brunner ◽  
Balaji Rajagopalan ◽  
William Kleiber

Abstract. Timely projections of seasonal streamflow extremes can be useful for the early implementation of annual flood risk adaptation strategies. However, predicting seasonal extremes is challenging particularly under non-stationary conditions and if extremes are connected in space. The goal of this study is to implement a space-time model for projection of seasonal streamflow extremes that considers the nonstationarity and spatio-temporal dependence of high flows. We develop a space-time model to project seasonal streamflow extremes for several lead times up to 2 months using a Bayesian Hierarchical Modelling (BHM) framework. This model is based on the assumption that streamflow extremes (3-day maxima) at a set of gauge locations are realizations of a Gaussian elliptical copula and generalized extreme value (GEV) margins with nonstationary parameters. These parameters are modeled as a linear function of suitable covariates from the previous season selected using the deviance information criterion (DIC). Finally, the copula is used to generate streamflow ensembles, which capture spatio-temporal variability and uncertainty. We apply this modelling framework to predict 3-day maximum flow in spring (May-June) at seven gauges in the Upper Colorado River Basin (UCRB) with 0 to 2 months lead time. In this basin, almost all extremes that cause severe flooding occur in spring as a result of snowmelt and precipitation. Therefore, we use regional mean snow water equivalent and temperature from the preceding winter season as well as indices of large-scale climate teleconnections – ENSO, AMO, and PDO – as potential covariates for 3-day maximum flow. Our model evaluation, which is based on the comparison of different model versions and the energy skill score, indicates that the model can capture the space-time variability of extreme flow well and that model skill increases with decreasing lead time. We also find that the use of climate variables slightly enhances skill relative to using only snow information. Median projections and their uncertainties are consistent with observations thanks to the representation of spatial dependencies through covariates in the margins and a Gaussian copula. This spatio-temporal modeling framework helps to plan seasonal adaptation and preparedness measures as predictions of extreme spring flows become available 2 months before actual flood occurrence.


2018 ◽  
Vol 34 (1) ◽  
pp. 15-30 ◽  
Author(s):  
Corey K. Potvin ◽  
Chris Broyles ◽  
Patrick S. Skinner ◽  
Harold E. Brooks ◽  
Erik Rasmussen

Abstract The Storm Prediction Center (SPC) tornado database, generated from NCEI’s Storm Data publication, is indispensable for assessing U.S. tornado risk and investigating tornado–climate connections. Maximizing the value of this database, however, requires accounting for systemically lower reported tornado counts in rural areas owing to a lack of observers. This study uses Bayesian hierarchical modeling to estimate tornado reporting rates and expected tornado counts over the central United States during 1975–2016. Our method addresses a serious solution nonuniqueness issue that may have affected previous studies. The adopted model explains 73% (>90%) of the variance in reported counts at scales of 50 km (>100 km). Population density explains more of the variance in reported tornado counts than other examined geographical covariates, including distance from nearest city, terrain ruggedness index, and road density. The model estimates that approximately 45% of tornadoes within the analysis domain were reported. The estimated tornado reporting rate decreases sharply away from population centers; for example, while >90% of tornadoes that occur within 5 km of a city with population > 100 000 are reported, this rate decreases to <70% at distances of 20–25 km. The method is directly extendable to other events subject to underreporting (e.g., severe hail and wind) and could be used to improve climate studies and tornado and other hazard models for forecasters, planners, and insurance/reinsurance companies, as well as for the development and verification of storm-scale prediction systems.


2020 ◽  
Author(s):  
Claudie Beaulieu ◽  
Matthew Hammond ◽  
Stephanie Henson ◽  
Sujit Sahu

<p>Assessing ongoing changes in marine primary productivity is essential to determine the impacts of climate change on marine ecosystems and fisheries. Satellite ocean color sensors provide detailed coverage of ocean chlorophyll in space and time, now with a combined record length of just over 20 years. Detecting climate change impacts is hindered by the shortness of the record and the long timescale of memory within the ocean such that even the sign of change in ocean chlorophyll is still inconclusive from time-series analysis of satellite data. Here we use a Bayesian hierarchical space-time model to estimate long-term trends in ocean chlorophyll. The main advantage of this approach comes from the principle of ”borrowing strength” from neighboring grid cells in a given region to improve overall detection. We use coupled model simulations from the CMIP5 experiment to form priors to provide a “first guess” on observational trend estimates and their uncertainty that we then update using satellite observations. We compare the results with estimates obtained with the commonly used vague prior, reflecting the case where no independent knowledge is available.  A global average net positive chlorophyll trend is found, with stronger regional trends that are typically positive in high and mid latitudes, and negative at low latitudes outside the Atlantic. The Bayesian hierarchical model used here provides a framework for integrating different sources of data for detecting trends and estimating their uncertainty in studies of global change.</p>


PLoS ONE ◽  
2021 ◽  
Vol 16 (7) ◽  
pp. e0253926
Author(s):  
Xiang Zhang ◽  
Taolin Yuan ◽  
Jaap Keijer ◽  
Vincent C. J. de Boer

Background Mitochondrial dysfunction is involved in many complex diseases. Efficient and accurate evaluation of mitochondrial functionality is crucial for understanding pathology as well as facilitating novel therapeutic developments. As a popular platform, Seahorse extracellular flux (XF) analyzer is widely used for measuring mitochondrial oxygen consumption rate (OCR) in living cells. A hidden feature of Seahorse XF OCR data is that it has a complex data structure, caused by nesting and crossing between measurement cycles, wells and plates. Surprisingly, statistical analysis of Seahorse XF data has not received sufficient attention, and current methods completely ignore the complex data structure, impairing the robustness of statistical inference. Results To rigorously incorporate the complex structure into data analysis, here we developed a Bayesian hierarchical modeling framework, OCRbayes, and demonstrated its applicability based on analysis of published data sets. Conclusions We showed that OCRbayes can analyze Seahorse XF OCR experimental data derived from either single or multiple plates. Moreover, OCRbayes has potential to be used for diagnosing patients with mitochondrial diseases.


2020 ◽  
Author(s):  
WILLIAM KLEIBER ◽  
Álvaro Ossandón ◽  
Balaji Rajagopalan ◽  
Manuela Brunner

2021 ◽  
Author(s):  
Xiang Zhang ◽  
Taolin Yuan ◽  
Jaap Keijer ◽  
Vincent C. J. de Boer

Mitochondrial dysfunction is involved in many complex diseases. Efficient and accurate evaluation of mitochondrial functionality is crucial for understanding pathology as well as facilitating novel therapeutic developments. As a popular platform, Seahorse extracellular flux (XF) analyzer is widely used for measuring mitochondrial oxygen consumption rate (OCR) in living cells. A hidden feature of Seahorse XF OCR data is that it has a complex data structure, caused by nesting and crossing between measurement cycles, wells and plates. Surprisingly, statistical analysis of Seahorse XF data has not received sufficient attention, and current methods completely ignore the complex data structure, impairing the robustness of statistical inference. To rigorously incorporate the complex structure into data analysis, here we developed a Bayesian hierarchical modeling framework, OCRbayes, and demonstrated its applicability based on analysis of published data sets. We showed that OCRbayes can analyze Seahorse XF OCR experimental data derived from either single or multiple plates. Moreover, OCRbayes has potential to be used for diagnosing patients with mitochondrial diseases.


2021 ◽  
Vol 58 (1) ◽  
pp. 42-67 ◽  
Author(s):  
Mads Stehr ◽  
Anders Rønn-Nielsen

AbstractWe consider a space-time random field on ${{\mathbb{R}^d} \times {\mathbb{R}}}$ given as an integral of a kernel function with respect to a Lévy basis with a convolution equivalent Lévy measure. The field obeys causality in time and is thereby not continuous along the time axis. For a large class of such random fields we study the tail behaviour of certain functionals of the field. It turns out that the tail is asymptotically equivalent to the right tail of the underlying Lévy measure. Particular examples are the asymptotic probability that there is a time point and a rotation of a spatial object with fixed radius, in which the field exceeds the level x, and that there is a time interval and a rotation of a spatial object with fixed radius, in which the average of the field exceeds the level x.


Sign in / Sign up

Export Citation Format

Share Document