Filtering Turbulent Sparsely Observed Geophysical Flows

Harlim, John; Majda, Andrew J.
April 2010
Monthly Weather Review;Apr2010, Vol. 138 Issue 4, p1050
Academic Journal
Filtering sparsely turbulent signals from nature is a central problem of contemporary data assimilation. Here, sparsely observed turbulent signals from nature are generated by solutions of two-layer quasigeostrophic models with turbulent cascades from baroclinic instability in two separate regimes with varying Rossby radius mimicking the atmosphere and the ocean. In the “atmospheric” case, large-scale turbulent fluctuations are dominated by barotropic zonal jets with non-Gaussian statistics while the “oceanic” case has large-scale blocking regime transitions with barotropic zonal jets and large-scale Rossby waves. Recently introduced, cheap radical linear stochastic filtering algorithms utilizing mean stochastic models (MSM1, MSM2) that have judicious model errors are developed here as well as a very recent cheap stochastic parameterization extended Kalman filter (SPEKF), which includes stochastic parameterization of additive and multiplicative bias corrections “on the fly.” These cheap stochastic reduced filters as well as a local least squares ensemble adjustment Kalman filter (LLS-EAKF) are compared on the test bed with 36 sparse regularly spaced observations for their skill in recovering turbulent spectra, spatial pattern correlations, and RMS errors. Of these four algorithms, the cheap SPEKF algorithm has the superior overall skill on the stringent test bed, comparable to LLS-EAKF in the atmospheric regime with and without model error and far superior to LLS-EAKF in the ocean regime. LLS-EAKF has special difficulty and high computational cost in the ocean regime with small Rossby radius, which creates stiffness in the perfect dynamics. The even cheaper mean stochastic model, MSM1, has high skill, which is comparable to SPEKF for the oceanic case while MSM2 has significantly worse filtering performance than MSM1 with the same inexpensive computational cost. This is interesting because MSM1 is based on a simple new regression strategy while MSM2 relies on the conventional regression strategy used in stochastic models for shear turbulence.


Related Articles

  • Evolutionary Estimation of Macro-Level Diffusion Models Using Genetic Algorithms: An Alternative to Nonlinear Least Squares. Venkatesan, Rajkumar; Krishnan, Trichy V.; Kumar, K. // Marketing Science;Summer2004, Vol. 23 Issue 3, p451 

    In this paper, we provide theoretical arguments and empirical evidence for how Genetic Algorithms (GA) can be used for efficient estimation of macro-level diffusion models. Using simulations we find that GA and Sequential Search-Based-Nonlinear Least Squares (SSB-NLS) provide comparable...

  • A Simple Algorithm for the Restoration of Clipped Speech Signal. Dahimene, Abdelhakim; Noureddine, Mohamed; Azrar, Aarab // Informatica (03505596);Aug2008, Vol. 32 Issue 3, p183 

    This paper deals with the problem of peak clipped speech. Our basic assumption is that the clipped speech is voiced and can be linearly predicted with a high accuracy. The coefficients of linear prediction are computed using two different algorithms: a least square direct method and a recursive...

  • Moving long baseline positioning algorithm with uncertain sound speed. Yan, Weisheng; Chen, Wei; Cui, Rongxin // Journal of Mechanical Science & Technology;Sep2015, Vol. 29 Issue 9, p3995 

    This paper presents a Moving long baseline (MLBL) positioning algorithm for the underwater target considering the uncertain underwater sound speed. First, the positioning and the sound speed models are established. To tackle the uncertain sound speed, a Uncertain least squares (ULS) positioning...

  • Conditional quantiles with varying Gaussians. Xiang, Dao-Hong // Advances in Computational Mathematics;May2013, Vol. 38 Issue 4, p723 

    In this paper we study conditional quantile regression by learning algorithms generated from Tikhonov regularization schemes associated with pinball loss and varying Gaussian kernels. Our main goal is to provide convergence rates for the algorithm and illustrate differences between the...

  • New Method of Spatial Extrapolation of Meteorological Fields on the Mesoscale Level Using a Kalman Filter Algorithm for a Four-Dimensional Dynamic�Stochastic Model. Komarov, V. S.; Lavrinenko, A. V.; Kreminskii, A. V.; Lomakina, N. Ya.; Popov, Yu. B.; Popova, A. I. // Journal of Atmospheric & Oceanic Technology;Feb2007, Vol. 24 Issue 2, p182 

    A new method and an algorithm of spatial extrapolation of mesometeorological fields to a territory uncovered with observations are suggested. The algorithm uses a linear Kalman filter for a four-dimensional dynamic�stochastic model of space�time variations of the atmospheric...

  • Scalable Implementations of Ensemble Filter Algorithms for Data Assimilation. Anderson, Jeffrey L.; Collins, Nancy // Journal of Atmospheric & Oceanic Technology;Aug2007, Vol. 24 Issue 8, p1452 

    A variant of a least squares ensemble (Kalman) filter that is suitable for implementation on parallel architectures is presented. This parallel ensemble filter produces results that are identical to those from sequential algorithms already described in the literature when forward observation...

  • Constrained Least Square Estimation Algorithm for Multisensor Bearings-Only Passive Target Tracking. Zhiwu XUAN; Kai NIE // Applied Mechanics & Materials;2014, Issue 602-605, p1842 

    The sensor needs to maneuver to get better observability in Bearings-Only passive target tracking with single sensor which makes the observation time longer. Multisensor Bearings-Only passive target tracking can solve the problem using exchange data. So the constrained Least Square...

  • New Matrix-weight IMM Algorithm Based on Coefficients Filters. Sha Junchen; Xing Jianping; Gao Liang; Liang Haozhe // Advances in Information Sciences & Service Sciences;Aug2012, Vol. 4 Issue 15, p111 

    In this paper, a fast interacting multiple model algorithm weighted by matrix (FMIMM) is presented based on the weighted least squares (WLS) method, which employs the steady state of the second and third order Kalman filters. And a CA model is introduced into the proposed algorithm to improve...

  • Research on Embedded System with Implementation of a Moving Object Tracking Algorithm based on Improved Meanshift on DM6437. Liu Yanfei; Li Qi; Fang Hao; Xu Huachun // Advanced Materials Research;2014, Vol. 1003, p207 

    To improve the real-time performance of the meanshift algorithm in the embedded system, an improved meanshift algorithm for tracking moving target is proposed in this paper. In order to reduce the influence of background pixel in a target model, the target model is build by using the target...


Read the Article


Sorry, but this item is not currently available from your library.

Try another library?
Sign out of this library

Other Topics