An official website of the United States government
Abstract We develop a space‐time Bayesian hierarchical modeling (BHM) framework for two flood risk attributes—seasonal daily maximum flow and the number of events that exceed a threshold during a season (NEETM)—at a suite of gauge locations on a river network. The model uses generalized extreme value (GEV) and Poisson distributions as marginals for these flood attributes with non‐stationary parameters. The rate parameters of the Poisson distribution and location, scale, and shape parameters of the GEV are modeled as linear functions of suitable covariates. Gaussian copulas are applied to capture the spatial dependence. The best covariates are selected using the Watanabe‐Akaike information criterion (WAIC). The modeling framework results in the posterior distribution of the flood attributes at all the gauges and various lead times. We demonstrate the utility of this modeling framework to forecast the flood risk attributes during the summer peak monsoon season (July‐August) at five gauges in the Narmada River basin (NRB) of West‐Central India for several lead times (0–3 months). As potential covariates, we consider climate indices such as El Niño–Southern Oscillation (ENSO), the Indian Ocean Dipole (IOD), and the Pacific Warm Pool Region (PWPR) from antecedent seasons, which have shown strong teleconnections with the Indian monsoon. We also include new indices related to the East Pacific and West Indian Ocean regions depending on the lead times. We show useful long lead skill from this modeling approach which has a strong potential to enable robust risk‐based flood mitigation and adaptation strategies 3 months before flood occurrences.
more »
« less
This content will become publicly available on December 1, 2025
Spatial and Temporal Bayesian Hierarchical Model Over Large Domains With Application to Holocene Sea Surface Temperature Reconstruction in the Equatorial Pacific
We present a novel space‐time Bayesian hierarchical model (BHM) to reconstruct annual Sea Surface Temperature (SST) over a large domain based on SST at limited proxy (i.e., sediment core) locations. The model is tested in the equatorial Pacific. The BHM leverages Principal Component Analysis to identify dominant space‐time modes of contemporary variability of the SST field at the proxy locations and employs these modes in a Gaussian process framework to estimate SSTs across the entire domain. The BHM allows us to model the mean field and covariance, varying in space and time in the process layers of the hierarchy. Using the Markov Chain Monte Carlo (MCMC) method and suitable priors on the model parameters, posterior distributions of the model parameters and, consequently, posterior distributions of the SST fields and the attendant uncertainties are obtained for any desired year. The BHM is calibrated and validated in the contemporary period (1854–2014) and subsequently applied to reconstruct SST fields during the Holocene (0–10 ka). Results are consistent with prior inferences of La Niña‐like conditions during the Holocene. This modeling framework opens exciting prospects for modeling and reconstruction of other fields, such as precipitation, drought indices, and vegetation.
more »
« less
- Award ID(s):
- 2310487
- PAR ID:
- 10623888
- Publisher / Repository:
- American Geophysical Union
- Date Published:
- Journal Name:
- Paleoceanography and Paleoclimatology
- Volume:
- 39
- Issue:
- 12
- ISSN:
- 2572-4517
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Model‐based projections of hydroclimate in western North America (wNA) remain uncertain and depend on how Pacific sea surface temperature (SST) will evolve in the future. However, whether climate models can accurately capture Pacific SST changes and its relationship with wNA hydroclimate in the future remains elusive. Here, we use a synthesis of proxy records and idealized model simulations to elucidate the spatiotemporal evolution and the forcings that drive wNA hydroclimate and Pacific SST during the Holocene (past ∼11,000 years), when the boundary conditions are different from the present. We find that wNA hydroclimate and Pacific SST co‐evolved during the Holocene, where wNA became wetter while the eastern equatorial Pacific and the north Pacific became warmer toward the present. We attribute changes in wNA hydroclimate to precession and carbon dioxide changes, but we are unable to attribute Pacific SST changes unambiguously to any forcing. Our analysis offers a framework to understand the relationship between wNA hydroclimate and Pacific SST and provides an empirical assessment of how these two regions are related over time.more » « less
-
Yamashita, Y.; Kano, M. (Ed.)Bayesian hybrid models (BHMs) fuse physics-based insights with machine learning constructs to correct for systematic bias. In this paper, we demonstrate a scalable computational strategy to embed BHMs in an equation-oriented modelling environment. Thus, this paper generalizes stochastic programming, which traditionally focuses on aleatoric uncertainty (as characterized by a probability distribution for uncertainty model parameters) to also consider epistemic uncertainty, i.e., mode-form uncertainty or systematic bias as modelled by the Gaussian process in the BHM. As an illustrative example, we consider ballistic firing using a BHM that includes a simplified glass-box (i.e., equation-oriented) model that neglects air resistance and a Gaussian process model to account for systematic bias (i.e., epistemic or model-form uncertainty) induced from the model simplification. The gravity parameter and the GP hypermeters are inferred from data in a Bayesian framework, yielding a posterior distribution. A novel single-stage stochastic program formulation using the posterior samples and Gaussian quadrature rules is proposed to compute the optimal decisions (e.g., firing angle and velocity) that minimize the expected value of an objective (e.g., distance from a stationary target). PySMO is used to generate expressions for the GP prediction mean and uncertainty in Pyomo, enabling efficient optimization with gradient-based solvers such as Ipopt. A scaling study characterizes the solver time and number of iterations for up to 2,000 samples from the posterior.more » « less
-
Abstract We present a Bayesian hierarchical space‐time stochastic weather generator (BayGEN) to generate daily precipitation and minimum and maximum temperatures. BayGEN employs a hierarchical framework with data, process, and parameter layers. In the data layer, precipitation occurrence at each site is modeled using probit regression using a spatially distributed latent Gaussian process; precipitation amounts are modeled as gamma random variables; and minimum and maximum temperatures are modeled as realizations from Gaussian processes. The latent Gaussian process that drives the precipitation occurrence process is modeled in the process layer. In the parameter layer, the model parameters of the data and process layers are modeled as spatially distributed Gaussian processes, consequently enabling the simulation of daily weather at arbitrary (unobserved) locations or on a regular grid. All model parameters are endowed with weakly informative prior distributions. The No‐U Turn sampler, an adaptive form of Hamiltonian Monte Carlo, is used to maximize the model likelihood function and obtain posterior samples of each parameter. Posterior samples of the model parameters propagate uncertainty to the weather simulations, an important feature that makes BayGEN unique compared to traditional weather generators. We demonstrate the utility of BayGEN with application to daily weather generation in a basin of the Argentine Pampas. Furthermore, we evaluate the implications of crop yield by driving a crop simulation model with weather simulations from BayGEN and an equivalent non‐Bayesian weather generator.more » « less
-
The latitudinal temperature gradient is a fundamental state parameter of the climate system tied to the dynamics of heat transport and radiative transfer. Thus, it is a primary target for temperature proxy reconstructions and global climate models. However, reconstructing the latitudinal temperature gradient in past climates remains challenging due to the scarcity of appropriate proxy records and large proxy–model disagreements. Here, we develop methods leveraging an extensive compilation of planktonic foraminifera δ 18 O to reconstruct a continuous record of the latitudinal sea-surface temperature (SST) gradient over the last 95 million years (My). We find that latitudinal SST gradients ranged from 26.5 to 15.3 °C over a mean global SST range of 15.3 to 32.5 °C, with the highest gradients during the coldest intervals of time. From this relationship, we calculate a polar amplification factor (PAF; the ratio of change in >60° S SST to change in global mean SST) of 1.44 ± 0.15. Our results are closer to model predictions than previous proxy-based estimates, primarily because δ 18 O-based high-latitude SST estimates more closely track benthic temperatures, yielding higher gradients. The consistent covariance of δ 18 O values in low- and high-latitude planktonic foraminifera and in benthic foraminifera, across numerous climate states, suggests a fundamental constraint on multiple aspects of the climate system, linking deep-sea temperatures, the latitudinal SST gradient, and global mean SSTs across large changes in atmospheric CO 2 , continental configuration, oceanic gateways, and the extent of continental ice sheets. This implies an important underlying, internally driven predictability of the climate system in vastly different background states.more » « less
