hierarchical spatio-temporal modeling

This paper addresses the challenge of jointly modeling spatially shared variation and localized temporal dynamics in areal spatiotemporal data. We propose a hierarchical Bayesian model grounded in spatially correlated Gaussian processes. Innovatively, we model time-varying parameters as a stochastic process with a conditional autoregressive (CAR) prior and explicitly incorporate dependence between temporal variability and temporal range, enabling more flexible and interpretable characterization of spatiotemporal evolution. Spatial information is propagated through shared variance components, allowing the model to capture cross-regional common patterns while preserving region-specific temporal trajectories. Bayesian inference is performed via a MALA–MH–Gibbs hybrid MCMC algorithm. Evaluated on real-world datasets—malaria incidence in Mozambique and food insecurity prevalence in Cameroon—the model achieves significantly higher predictive accuracy than conventional approaches, demonstrating clear utility for public health surveillance and policy decision-making.

Bayesian inference for high-dimensional spatio-temporal ordinal data—such as drought severity classifications—is computationally prohibitive, poorly parallelizable, and often necessitates oversimplified models. To address these challenges, we propose a two-stage decoupled MCMC algorithm: in Stage I, temporal ordinal regression models are fitted in parallel across spatial units; in Stage II, posterior samples are reweighted and spatially calibrated to integrate spatial dependence, enabling full spatio-temporal joint inference. Our approach preserves model fidelity while overcoming scalability bottlenecks inherent in conventional MCMC. It represents the first scalable, high-fidelity Bayesian framework for large-scale ordinal spatio-temporal modeling. Evaluated on U.S. drought data comprising hundreds of thousands of observations, the method achieves a tenfold speedup over standard single-stage MCMC, yields highly consistent posterior distributions, and significantly outperforms existing simplified alternatives—demonstrating practical viability for real-world, large-scale applications.

Modeling large-scale, non-Gaussian, noisy, and incomplete spatial data remains challenging due to limitations of conventional copula models in capturing complex spatial heterogeneity and non-Gaussian dependence structures. Method: We propose a Bayesian hierarchical model integrating vine copulas with spatial random effects. Crucially, we design a novel vine copula structure explicitly embedding spatial dependence, enabling low-rank latent process representations and computationally efficient Bayesian inference. Contribution/Results: Our method overcomes key bottlenecks in traditional copula-based spatial modeling. In both parameter estimation and spatial prediction tasks, it significantly outperforms benchmarks—including fixed-rank kriging (FRK)—in accuracy, convergence speed, and robustness. Applied to atmospheric methane concentration mapping over Australia’s Bowen Basin using Sentinel-5P satellite remote sensing data, the approach delivers superior predictive performance under missingness and noise. This work establishes a new paradigm for non-Gaussian spatial statistical modeling.

To address computational challenges in modeling non-Gaussian spatiotemporal data—specifically, the intractability of analytically integrating random effects and weak parameter identifiability impeding MCMC convergence—this paper proposes a Bayesian inference framework based on predictive stacking. Innovatively integrating generalized conjugate multivariate distribution theory with an extended Diaconis–Ylvisaker conjugate prior family, the method enables exact sampling from conditional posteriors and facilitates Bayesian model integration across disparate model configurations. It circumvents the computational bottlenecks inherent in conventional generalized linear mixed models, where closed-form integration of random effects is infeasible. Simulation studies demonstrate substantially improved MCMC convergence speed and parameter estimation accuracy relative to standard MCMC approaches. The framework is further validated on large-scale spatiotemporal count data from the North American Breeding Bird Survey, confirming its robustness and scalability.

This study addresses the challenges in modeling global and local effects within inhomogeneous pairwise interaction Gibbs point processes and the lack of effective methods for testing complete spatial randomness (CSR). To overcome these limitations, the authors propose a hierarchical Bayesian framework that, for the first time, integrates basis function expansions with Bayesian hierarchical modeling to flexibly characterize both the intensity and interaction functions. Building on posterior inference, they develop a Bayesian testing procedure specifically designed for CSR assessment. The approach enables efficient inference via Markov chain Monte Carlo (MCMC) and demonstrates strong empirical performance: when applied to water strider distribution and forest fire data, it successfully uncovers complex spatial dependence structures and provides reliable CSR tests, substantially enhancing the flexibility and inferential power of Gibbs point process modeling.

This work proposes a dynamic Bayesian framework endowed with a Markovian dependency structure to address the challenges of computational inefficiency and limited cross-temporal information sharing in high-dimensional multivariate spatiotemporal modeling. By integrating matrix-variate Gaussian distributions, dynamic linear models, and Bayesian predictive stacking—augmented with an adaptive Markov transition mechanism—the approach enables efficient online forward filtering and backward smoothing within a sequence-parallel hybrid architecture. The proposed method achieves exact inference while substantially enhancing scalability and dynamic adaptability, making it well-suited for large-scale, multivariate, and dynamically evolving spatiotemporal data streams requiring efficient online learning.

Traditional geostatistical methods rely on second-order moments and Gaussian assumptions, which are inadequate for capturing non-Gaussian spatial dependence. This work addresses this limitation by leveraging Sklar’s theorem to introduce a copula-based framework that decouples marginal distributions from the spatial dependence structure. The authors systematically develop spatial copula models applicable at both fixed point sets and process levels, emphasizing Kolmogorov consistency to clarify distinctions between these two modeling paradigms. The framework is further extended to spatio-temporal settings and supports flexible marginal specifications. By integrating spatial statistics, copula theory, and stochastic processes, this study establishes a unified approach for modeling non-Gaussian spatial dependence, elucidating the relationships, strengths, and limitations of existing methodologies, thereby advancing both theoretical understanding and practical applications in the field.

This study addresses the unresolved question of when spatial and non-spatial random effects yield equivalent posterior inference for regression coefficients in Bayesian regression with multilevel areal data. Within a hierarchical Bayesian framework assuming Gaussian responses and employing a Leroux conditional autoregressive (CAR) prior, the authors formally derive—for the first time—a closed-form sample size threshold $m^*$ that determines whether spatial modeling is necessary. This threshold admits a clear interpretation, revealing that the difference in posterior variances converges to zero at an $O(m^{-1})$ rate. Simulation studies confirm that $m^*$ accurately identifies the tipping point in modeling complexity; notably, spatial modeling remains essential regardless of sample size whenever covariates exhibit no within-area variation.

This work addresses the computational burden, streaming nature, and privacy constraints inherent in Bayesian inference for large-scale, complex spatiotemporal data by proposing a distributed recursive Bayesian inference framework based on Integrated Nested Laplace Approximation (INLA). The approach synergistically integrates distributed computing, federated learning, and recursive updating mechanisms to substantially reduce computational complexity while preserving inferential accuracy. Implemented within the R-INLA framework, the method enables scalable, adaptive, and privacy-preserving Bayesian analysis. Empirical evaluations across multiple case studies demonstrate that the proposed framework achieves performance comparable to centralized full-data inference, even under streaming data conditions and stringent privacy requirements.