Inference of abrupt changes in noisy geochemical records using transdimensional changepoint models

Abstract

International audienceWe present a method to quantify abrupt changes (or changepoints) in data series, represented as a function of depth or time. These changes are often the result of climatic or environmental variations and can be manifested inmultiple datasets as different responses, but all datasets can have the same changepoint locations/timings. The method we present uses transdimensional Markov chain Monte Carlo to infer probability distributions on the number and locations (in depth or time) of changepoints, the mean values between changepoints and, if required, the noise variance associated with each dataset being considered. This latter point is important as we generally will have limited information on the noise, such as estimates only of measurement uncertainty, and in most cases it is not practical to make repeat sampling/measurement to assess other contributions to the variation in the data.Wedescribe themain features of the approach (and describe themathematical formulation in supplementary material), and demonstrate its validity using synthetic datasets, with known changepoint structure (number and locations of changepoints) and distribution of noise variance for each dataset.We show that when using multiple data, we expect to achieve better resolution of the changepoint structure than when we use each dataset individually. This is conditional on the validity of the assumption of common changepoints between different datasets.We then apply themethod to two sets of real geochemical data, both from peat cores, taken from NE Australia and eastern Tibet. Under the assumption that changes occur at the same time for all datasets, we recover solutions consistent with those previously inferred qualitatively from independent data and interpretations. However, our approach provides a quantitative estimate of the relative probability of the inferred changepoints, allowing an objective assessment of the significance of each change

Transdimensional Love-wave tomography of the British Isles and shear-velocity structure of the East Irish Sea Basin from ambient-noise interferometry

“…In fact, the concept of transdimensionality is not limited to traveltime tomography but can be adapted to a number of different inverse problems including regression (Gallagher et al 2011), inversion of controlled source electromagnetic data (Ray et al 2014), inversion of surface-wave dispersion data (Young et al 2013), and joint inversion of surface-wave dispersion and receiver function data (Bodin et al 2012b). …”

Crustal and uppermost mantle structure variation beneath La Réunion hotspot track

“…Malinverno 2002;Sambridge et al 2006;Hopcroft et al 2009;Dettmer et al 2010;Piana Agostinetti & Malinverno 2010;Gallagher et al 2011;Bodin et al 2012a,b;Iaffaldano et al 2012). TB inversion of RFs is a fully non-linear inversion approach cast in a Bayesian (i.e.…”

Taming uncertainty in geophysical inversion

“…Further afield this approach has been applied in Hydrology, Minsley (2011), signal processing, Kolb and Lekić (2014), environmental monitoring for exploration, Warner et al (2015), inversion of geoid data for viscosity, Rudolph et al (2016), various problems in Geochemistry including estimating thermal properties of the Lithosphere, Gallagher et al (2009Gallagher et al ( , 2011, plate tectonic reconstructions, Iaffaldano et al (2012), and a variety of Earth Science regression problems, similar to that in FIGURE 1, such as paleo sea-level reconstruction, Lambeck et al (2014), and even to estimating the variable rotation of the Earth's solid inner core, Tkalčić et al (2013) and paleomagnetic time series analysis, Ingham et al (2014). The versatility and power of the methodology have given it wide appeal and are leading to increasing numbers of novel applications.…”