Spatiotemporal modeling of European paleoclimate using doubly sparse Gaussian processes

Seth D. Axen, Alexandra Gessner, Christian Sommer, Nils Weitzel, Álvaro Tejero-Cantero

arXiv preprint · 15 November 2022

In one sentence. Doubly sparse spatiotemporal Gaussian process builds a probabilistic European climate-field reconstruction from the Last Glacial Maximum to mid-Holocene, combining simulations with pollen proxies.

Abstract

Paleoclimatology -- the study of past climate -- is relevant beyond climate science itself, such as in archaeology and anthropology for understanding past human dispersal. Information about the Earth's paleoclimate comes from simulations of physical and biogeochemical processes and from proxy records found in naturally occurring archives. Climate-field reconstructions (CFRs) combine these data into a statistical spatial or spatiotemporal model. To date, there exists no consensus spatiotemporal paleoclimate model that is continuous in space and time, produces predictions with uncertainty, and can include data from various sources. A Gaussian process (GP) model would have these desired properties; however, GPs scale unfavorably with data of the magnitude typical for building CFRs. We propose to build on recent advances in sparse spatiotemporal GPs that reduce the computational burden by combining variational methods based on inducing variables with the state-space formulation of GPs. We successfully employ such a doubly sparse GP to construct a probabilistic model of European paleoclimate from the Last Glacial Maximum (LGM) to the mid-Holocene (MH) that synthesizes paleoclimate simulations and fossilized pollen proxy data.

Labels

Labels and summary are generated by Claude and reviewed weekly. Wrong? Edit this entry.