Abstract:We develop a new computational approach for the speedy and accurate recovery of unknown spatially varying parameters of (stochastic) PDE models based on the sequentially arriving noisy observations of the evolving system variables. Two essential ingredients in our strategy for the identification of the target spatial fields are given by the Fourier diagonalization (FD) method and the hierarchical Bayesian model. We first apply the FD framework to perform an effectively parallelized computation in quantifying the propagating uncertainties of the dynamic variables in the Fourier space, and to facilitate a drastic saving of computational resources required for the Bayesian estimation of the associated parameters via sequential data assimilation. Yet our case study reveals that a very poor performance of the FD scheme occurs when the fast and slow variables coexist in the system evolution. The key observation is that, estimating the parameters in association with some slow variables, the FD-based Bayesian solver significantly underperforms compared to the remaining cases of faster variables. Due to this highly non-uniform discrepancy in the accuracy across distinct Fourier modes, the approximation of the parameter field cannot be so desired if the consequence is represented in the physical space. As one effort to circumvent this problem, we provide a systematic approach for a radical improvement of the naive FD technique. To do this, we make use of the hierarchical Bayesian model; it refers to the process of gradually enriching the knowledge of the unknown spatial fields from coarse to medium and to fine resolutions by using the Bayesian inference obtained at coarser levels to provide prior information for the estimation at finer levels. One major contribution of this paper is the development of a variant of the classical application of the multi-resolution approach through the close integration with the FD method, leading to the emergence of a new Bayesian paradigm for the data-driven parametric identification. Numerical experiments are performed to corroborate our demonstration concerning the efficacy and effectiveness of the proposed algorithm in obtaining a good degree of accuracy together with a significantly reduced computational cost. (C) 2020 Elsevier Inc. All rights reserved.

High-Dimensional Density Estimation via SCA: An Example in the Modelling of Hurricane Tracks

Hierarchical Density Decompositions for Abnormal Event Diagnosis in Serially Correlated Non-Gaussian Systems

Stochastic Simulation of Tropical Cyclones for Risk Assessment at One Go: A Multivariate Functional PCA Approach

Generalized early warning signals in multivariate and gridded data with an application to tropical cyclones

Steps Towards Modeling Community Resilience under Climate Change: Hazard Model Development

Tipping in a Low-Dimensional Model of a Tropical Cyclone

A Bayesian Hierarchical Model Framework to Quantify Uncertainty of Tropical Cyclone Precipitation Forecasts

Dimensionality Reduction for Onboard Modeling of Uncertain Atmospheres

Spatio-temporal methods for estimating subsurface ocean thermal response to tropical cyclones

Data-driven dimensionality reduction and causal inference for spatiotemporal climate fields

Nonparametric Density Estimation for Data Scattered on Irregular Spatial Domains: A Likelihood-Based Approach Using Bivariate Penalized Spline Smoothing

Modeling the Temporal Correlation in Hurricane Frequency for Damage Assessment of Residential Structures Subjected to Climate Change

Estimating intraevent and interevent variability and spatial correlation of tropical cyclone wind fields and their use for the risk assessment of a portfolio of structures

Local likelihood estimation of complex tail dependence structures, applied to U.S. precipitation extremes

Collective Spectral Density Estimation and Clustering for Spatially-Correlated Data

Hierarchical Sparse Observation Models and Informative Prior for Bayesian Inference of Spatially Varying Parameters.

On the generation of high‐resolution probabilistic design events capturing the joint occurrence of rainfall and storm surge in coastal basins

Improving Piecewise Linear Snow Density Models through Hierarchical Spatial and Orthogonal Functional Smoothing

Kernel density estimation applied to tropical cyclones genesis in northwestern pacific

Nonparametric Conditional Density Estimation In A Deep Learning Framework For Short-Term Forecasting

Distributed Inference for Spatial Extremes Modeling in High Dimensions