Estimating the state of a geophysical system with sparse observations: time delay methods to achieve accurate initial states for prediction
Abstract
The problem of forecasting the behavior of a complex dynamical system through analysis of observational timeseries data becomes difficult when the system expresses chaotic behavior and the measurements are sparse, in both space and/or time. Despite the fact that this situation is quite typical across many fields, including numerical weather prediction, the issue of whether the available observations are "sufficient" for generating successful forecasts is still not well understood. An analysis by Whartenby et al. (2013) found that in the context of the nonlinear shallow water equations on a β plane, standard nudging techniques require observing approximately 70 % of the full set of state variables. Here we examine the same system using a method introduced by Rey et al. (2014a), which generalizes standard nudging methods to utilize time delayed measurements. Here, we show that in certain circumstances, it provides a sizable reduction in the number of observations required to construct accurate estimates and highquality predictions. In particular, we find that this estimate of 70 % can be reduced to about 33 % using time delays, and even further if Lagrangian drifter locations are also used as measurements.
 Authors:
 Univ. of California, San Diego, CA (United States). Dept. of Physics
 Univ. of California, San Diego, CA (United States). Dept. of Physics; Univ. of California, San Diego, CA (United States). Scripps Inst. of Oceanography, Marine Physical Lab.
 Publication Date:
 Research Org.:
 Krell Inst., Ames, IA (United States)
 Sponsoring Org.:
 USDOE; National Science Foundation (NSF); US Department of the Navy, Office of Naval Research (ONR)
 OSTI Identifier:
 1366514
 Grant/Contract Number:
 FG0297ER25308; PHY0961153
 Resource Type:
 Journal Article: Accepted Manuscript
 Journal Name:
 Nonlinear Processes in Geophysics (Online)
 Additional Journal Information:
 Journal Name: Nonlinear Processes in Geophysics (Online); Journal Volume: 24; Journal Issue: 1; Journal ID: ISSN 16077946
 Publisher:
 European Geosciences Union  Copernicus
 Country of Publication:
 United States
 Language:
 English
 Subject:
 58 GEOSCIENCES
Citation Formats
An, Zhe, Rey, Daniel, Ye, Jingxin, and Abarbanel, Henry D. I. Estimating the state of a geophysical system with sparse observations: time delay methods to achieve accurate initial states for prediction. United States: N. p., 2017.
Web. doi:10.5194/npg2492017.
An, Zhe, Rey, Daniel, Ye, Jingxin, & Abarbanel, Henry D. I. Estimating the state of a geophysical system with sparse observations: time delay methods to achieve accurate initial states for prediction. United States. doi:10.5194/npg2492017.
An, Zhe, Rey, Daniel, Ye, Jingxin, and Abarbanel, Henry D. I. Mon .
"Estimating the state of a geophysical system with sparse observations: time delay methods to achieve accurate initial states for prediction". United States.
doi:10.5194/npg2492017. https://www.osti.gov/servlets/purl/1366514.
@article{osti_1366514,
title = {Estimating the state of a geophysical system with sparse observations: time delay methods to achieve accurate initial states for prediction},
author = {An, Zhe and Rey, Daniel and Ye, Jingxin and Abarbanel, Henry D. I.},
abstractNote = {The problem of forecasting the behavior of a complex dynamical system through analysis of observational timeseries data becomes difficult when the system expresses chaotic behavior and the measurements are sparse, in both space and/or time. Despite the fact that this situation is quite typical across many fields, including numerical weather prediction, the issue of whether the available observations are "sufficient" for generating successful forecasts is still not well understood. An analysis by Whartenby et al. (2013) found that in the context of the nonlinear shallow water equations on a β plane, standard nudging techniques require observing approximately 70 % of the full set of state variables. Here we examine the same system using a method introduced by Rey et al. (2014a), which generalizes standard nudging methods to utilize time delayed measurements. Here, we show that in certain circumstances, it provides a sizable reduction in the number of observations required to construct accurate estimates and highquality predictions. In particular, we find that this estimate of 70 % can be reduced to about 33 % using time delays, and even further if Lagrangian drifter locations are also used as measurements.},
doi = {10.5194/npg2492017},
journal = {Nonlinear Processes in Geophysics (Online)},
number = 1,
volume = 24,
place = {United States},
year = {Mon Jan 16 00:00:00 EST 2017},
month = {Mon Jan 16 00:00:00 EST 2017}
}
Web of Science

Previously, the authors reported an equation of state (EOS) modeling approach that successfully calculated the PVTX properties of supercritical fluid mixtures. The model is based on a corresponding states assumption applied to a highly accurate EOS for the reference CH{sub 4} system. The CH{sub 4} EOS was parameterized from 273 to 723 K and 1 to 3,000 bar by using experimental PVT data. Molecular dynamics simulated PVT data were used to extend the parameterization in the CH{sub 4} system to 2,000 K and 20 kbar. Mixing in the H{sub 2}OCO{sub 2}CH{sub 4}N{sub 2} system was successfully described by using amore »

New spectral methods of estimating delay time in reactor noise investigations
In experiments using neutron noise methods to estimate delay time in reactor noise investigations, the random processes of two detectors independent of the physical nature of the noise field under investigation are represented by an expression containing a local component, which conveys information on the delay, and a global component, which identically affects both detectors. In cases where the global component is dominant, the abscissa of the mutual correlation function global maximum is equal to zero, and the spectral density function is far from a linear relationship, which makes traditional methods unsuitable for obtaining the delay time. In this paper,more » 
A statespace Bayesian framework for estimating biogeochemical transformations using timelapse geophysical data
We develop a statespace Bayesian framework to combine timelapse geophysical data with other types of information for quantitative estimation of biogeochemical parameters during bioremediation. We consider characteristics of endproducts of biogeochemical transformations as state vectors, which evolve under constraints of local environments through evolution equations, and consider timelapse geophysical data as available observations, which could be linked to the state vectors through petrophysical models. We estimate the state vectors and their associated unknown parameters over time using Markov chain Monte Carlo sampling methods. To demonstrate the use of the statespace approach, we apply it to complex resistivity data collected duringmore » 
Initial time singularities and admissible initial states for a system of coupled scalar fields
We discuss the problem of initial states for a system of coupled scalar fields out of equilibrium in the oneloop approximation. The fields consist of classical background fields, taken constant in space, and quantum fluctuations. If the initial state is the adiabatic vacuum, i.e., the ground state of a Fock space of particle excitations that diagonalize the mass matrix, the energymomentum tensor is infinite at t=0, its most singular part behaves as 1/t. When the system is coupled to gravity this presents a problem that we solve by a Bogoliubov transformation of the naive initial state. As a side resultmore » 
Resolving seismic anisotropy: Sparse matrix methods for geophysical inverse problems
Two techniques for the singular value decomposition (SVD) of large sparse matrices are directly applicable to geophysical inverse problems: subspace iteration and the Lanczos method. Both methods require very little incore storage and efficiently compute a set of singular values and singular vectors. A comparison of the singular value and vector estimates of these iterative approaches with the results of a conventional incore SVD algorithm demonstrates their accuracy. Hence, it is possible to conduct an assessment of inversion results for larger sparse inverse problems such as those arising in seismic tomography. As an example, the authors examine the resolution matrixmore »