A Bayesian Framework for Storm Tracking Using a Hidden-State Representation

Scharenbroich, Lucas; Magnusdottir, Gudrun; Smyth, Padhraic; Stern, Hal; Chia-chi Wang
June 2010
Monthly Weather Review;Jun2010, Vol. 138 Issue 6, p2132
Academic Journal
A probabilistic tracking model is introduced that identifies storm tracks from feature vectors that are extracted from meteorological analysis data. The model assumes that the genesis and lysis times of each track are unknown and estimates their values along with the track’s position and storm intensity over time. A hidden-state dynamics model (Kalman filter) characterizes the temporal evolution of the storms. The model uses a Bayesian methodology for estimating the unknown lifetimes (genesis–lysis pairs) and tracks of the storms. Prior distributions are placed over the unknown parameters and their posterior distributions are estimated using a Markov Chain Monte Carlo (MCMC) sampling algorithm. The posterior distributions are used to identify and report the most likely storm tracks in the data. This approach provides a unified probabilistic framework that accounts for uncertainty in storm timing (genesis and lysis), storm location and intensity, and the feature detection process. Thus, issues such as missing observations can be accommodated in a statistical manner without human intervention. The model is applied to the field of relative vorticity at the 975-hPa level of analysis from the National Centers for Environmental Prediction Global Forecast System during May–October 2000–02, in the tropical east Pacific. Storm tracks in the National Hurricane Center best-track data (HURDAT) for the same period are used to assess the performance of the storm identification and tracking model.


Related Articles

  • Numerical Bayesian inference with arbitrary prior. Tsionas, Efthymios G. // Statistical Papers;2000, Vol. 41 Issue 4, p437 

    Presents a study that examined how advances in Markov Chain Monte Carlo integration can facilitate the routine Bayesian analysis of the linear model when the prior distribution is completely user dependent. Background on the problem of specifying prior information in applying econometric...

  • Improved Extend Kalman Particle Filter Based On Markov Chain Monte Carlo for Nonlinear State Estimation. Huajian Wang // International Proceedings of Computer Science & Information Tech;2012, Vol. 50, p157 

    Considering the problem of poor tracking accuracy and particle degradation in the standard particle filter algorithm, a new improved extend kalman particle filter algorithm based on markov chain monte carlo(MCMC) is discussed. The algorithm uses the Extended Kalman filter to generate the...

  • Testing a Molecular Clock without an Outgroup: Derivations of Induced Priors on Branch-Length Restrictions in a Bayesian Framework. Suchard, Marc A.; Weiss, Robert E.; Sinsheimer, Janet S. // Systematic Biology;Feb2003, Vol. 52 Issue 1, p48 

    Proposes a Bayesian method for testing molecular clock hypotheses for use with aligned sequence data from multiple taxa. Use of a nonreversible nucleotide substitution model; Use of reversible jump Markov chain Monte Carlo to sample from the posterior distribution of the phylogenetic model...

  • Improving Marginal Likelihood Estimation for Bayesian Phylogenetic Model Selection. WANGANG XIE; LEWIS, PAUL O.; YU FAN; LYNN KUO; MING-HUI CHEN // Systematic Biology;Mar2011, Vol. 60 Issue 2, p150 

    The marginal likelihood is commonly used for comparing different evolutionary models in Bayesian phylogenetics and is the central quantity used in computing Bayes Factors for comparing model fit. A popular method for estimating marginal likelihoods, the harmonic mean (HM) method, can be easily...

  • Bayesian Mixture Model Analysis for Detecting Differentially Expressed Genes. Zhenyu Jia; Shizhong Xu // International Journal of Plant Genomics;2008, Vol. 2008, p1 

    Control-treatment design is widely used in microarray gene expression experiments. The purpose of such a design is to detect genes that express differentially between the control and the treatment. Many statistical procedures have been developed to detect differentially expressed genes, but all...

  • A recipe for Bayesian network driven stock assessment. Hammond, T. R. // Canadian Journal of Fisheries & Aquatic Sciences;Sep2004, Vol. 61 Issue 9, p1647 

    Markov Chain Monte Carlo (MCMC), the most widely used algorithm in Bayesian statistics, can fail to converge. Although convergence is tested by various diagnostics, these can only reveal failure, never success. To avoid these difficulties, this paper suggests a recipe for using Bayesian network...

  • R/qtlbim: QTL with Bayesian Interval Mapping in experimental crosses. Brian S. Yandell; Tapan Mehta; Samprit Banerjee; Daniel Shriner; Ramprasad Venkataraman; Jee Young Moon; W. Whipple Neely; Hao Wu; Randy von Smith; Nengjun Yi // Bioinformatics;Mar2007, Vol. 23 Issue 5, p641 

    Summary: R/qtlbim is an extensible, interactive environment for the Bayesian Interval Mapping of QTL, built on top of R/qtl (Broman et al., 2003), providing Bayesian analysis of multiple interacting quantitative trait loci (QTL) models for continuous, binary and ordinal traits in experimental...

  • Adaptive approximate Bayesian computation. Beaumont, Mark A.; Cornuet, Jean-Marie; Marin, Jean-Michel; Robert, Christian P. // Biometrika;Dec2009, Vol. 96 Issue 4, p983 

    Sequential techniques can enhance the efficiency of the approximate Bayesian computation algorithm, as in Sisson et al.’s (2007) partial rejection control version. While this method is based upon the theoretical works of Del Moral et al. (2006), the application to approximate Bayesian...

  • Evidence of a Major Gene From Bayesian Segregation Analyses of Liability to Osteochondral Diseases in Pigs. Kadarmideen, Haja N.; Janss, Luc L. G. // Genetics;Nov2005, Vol. 171 Issue 3, p1195 

    Bayesian segregation analyses were used to investigate the mode of inheritance oil osteochondral lesions (osteochondrosis, OC) in pigs. Data consisted of 1163 animals with OC and their pedigrees included 2891 animals. Mixed-inheritance threshold models (MITM) and several variants of MITM, in...


Read the Article


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

Try another library?
Sign out of this library

Other Topics