Extracting common pulse-like signals from multiple ice core time series

International audience To understand the nature and cause of natural climate variability, it is important to possess an accurate estimate of past climate forcings. Direct measurements that are reliable only exist for the past few decades. Therefore knowledge of prior variations has to be established...

Full description

Bibliographic Details
Published in:Computational Statistics & Data Analysis
Main Authors: Gazeaux, Julien, Batista, Deborah, Ammann, Caspar M., Naveau, Philippe, Jégat, Cyrille, Gao, Chaochao
Other Authors: TROPO - LATMOS, Laboratoire Atmosphères, Milieux, Observations Spatiales (LATMOS), Université de Versailles Saint-Quentin-en-Yvelines (UVSQ)-Université Pierre et Marie Curie - Paris 6 (UPMC)-Institut national des sciences de l'Univers (INSU - CNRS)-Centre National de la Recherche Scientifique (CNRS)-Université de Versailles Saint-Quentin-en-Yvelines (UVSQ)-Université Pierre et Marie Curie - Paris 6 (UPMC)-Institut national des sciences de l'Univers (INSU - CNRS)-Centre National de la Recherche Scientifique (CNRS), Department of Mathematical Sciences Denver, University of Colorado Denver, Research Applications Laboratory Boulder (RAL), National Center for Atmospheric Research Boulder (NCAR), Laboratoire des Sciences du Climat et de l'Environnement Gif-sur-Yvette (LSCE), Université de Versailles Saint-Quentin-en-Yvelines (UVSQ)-Institut national des sciences de l'Univers (INSU - CNRS)-Université Paris-Saclay-Centre National de la Recherche Scientifique (CNRS)-Direction de Recherche Fondamentale (CEA) (DRF (CEA)), Commissariat à l'énergie atomique et aux énergies alternatives (CEA)-Commissariat à l'énergie atomique et aux énergies alternatives (CEA), École des Mines de Paris, Department of Environmental Sciences New Brunswick, School of Environmental and Biological Sciences New Brunswick, Rutgers, The State University of New Jersey New Brunswick (RU), Rutgers University System (Rutgers)-Rutgers University System (Rutgers)-Rutgers, The State University of New Jersey New Brunswick (RU), Rutgers University System (Rutgers)-Rutgers University System (Rutgers)
Format: Article in Journal/Newspaper
Language:English
Published: HAL CCSD 2013
Subjects:
Online Access:https://hal.science/hal-00922271
https://doi.org/10.1016/j.csda.2012.01.024
Description
Summary:International audience To understand the nature and cause of natural climate variability, it is important to possess an accurate estimate of past climate forcings. Direct measurements that are reliable only exist for the past few decades. Therefore knowledge of prior variations has to be established based on indirect information derived from natural archives. The challenge has always been to find a strict objective method that can identify volcanic events and offer sound amplitude estimates in these noisy records. An automatic procedure is introduced here to estimate the magnitude of strong, but short-lived, volcanic signals from a suite of polar ice core series. Rather than treating records from individual ice cores separately and then averaging their respective magnitudes, our extraction algorithm jointly handles multiple time series to identify their common, but hidden, volcanic pulses. The statistical procedure is based on a multivariate multi-state space model. Exploiting the joint fluctuations, it provides an accurate estimator of the timing, peak magnitude and duration of individual pulse-like deposition events within a set of different series. This ensures a more effective separation of the real signals from spurious noise that can occur in any individual time series, and thus a higher sensitivity to identify smaller scale events. At the same time, it provides a measure of confidence through the posterior probability for each pulse-like event, indicating how well a pulse can be recognized against the background noise. The flexibility and robustness of our approach, as well as important underlying assumptions and remaining limitations, are discussed by applying our method to first simulated and then real world ice core time series.