Global-Krigger: A Global Kriging Interpolation Toolbox With Paleoclimatology Examples

dc.contributor.authorCosentino, N. J.
dc.contributor.authorOpazo, N. E.
dc.contributor.authorLambert, F.
dc.contributor.authorOsses, A.
dc.contributor.authorvan 't Wout, E.
dc.date.accessioned2025-01-20T20:10:44Z
dc.date.available2025-01-20T20:10:44Z
dc.date.issued2023
dc.description.abstractMany applications in Earth sciences require spatial prediction, that is, obtaining a continuous scalar field from a set of discrete scalar data points on the Earth's surface. Such applications include model data comparisons and derivation of continuous scalar fields as input for Earth system models. The advantage of kriging as an interpolation method is that it provides predictions with confidence intervals for data sets of irregularly distributed points in space. However, the theory of kriging for non-Euclidean domains such as oblate spheroids (e.g., the Earth's surface) is poorly developed, and existing kriging algorithms for global interpolation oftentimes cannot guarantee the validity of their predictions. Here, we present Global-Krigger, a new kriging interpolation algorithm adapted for local to global applications that (a) incorporates a numerical check to guarantee that the necessary conditions for the kriging system of linear equations are met, and (b) derives a combined uncertainty field due both to spatial variations in data density and measurement error. The robustness of the method is demonstrated by cross-validating predictions against reanalysis fields of traditional climatological scalar variables. We also show an example application in paleoclimatology for Holocene mineral dust deposition fluxes. The toolbox includes a user-friendly graphical user interface that guides users through a range of choices during data pre-interpolation analysis, kriging, and post-processing.Plain Language Summary Many climatic variables such as temperature and precipitation are measured by monitoring stations on the Earth's surface at irregularly distributed locations. Similarly, paleoclimatic proxy variables are available at specific locations where paleoclimatic archives such as marine sediment cores have been retrieved. Such data sets can be used to guide Earth system simulations, but this works best when they are available as a regular grid. Kriging is an interpolation method that produces regular grids from irregularly distributed measurements. However, while kriging theory is well-developed for measurements on a plane surface, its application on curved surfaces is problematic. To improve this, we have developed Global-Krigger, a user-friendly kriging interpolation toolbox for the surface of the Earth. It incorporates new tests to check the validity of kriging results on curved surfaces. The toolbox also provides a quantification of the interpolation uncertainty, taking into account both the uncertainty of the measurements and the spatial distribution of the measurement sites. We show that our algorithm gives robust results by applying it to a sub-sample of known modern climatology fields, and also provides a paleoclimatic application for Holocene mineral dust deposition fluxes. Finally, we develop a graphical user interface to guide users in the process of kriging.
dc.fuente.origenWOS
dc.identifier.doi10.1029/2022GC010821
dc.identifier.eissn1525-2027
dc.identifier.urihttps://doi.org/10.1029/2022GC010821
dc.identifier.urihttps://repositorio.uc.cl/handle/11534/92092
dc.identifier.wosidWOS:001003163400001
dc.issue.numero6
dc.language.isoen
dc.revistaGeochemistry geophysics geosystems
dc.rightsacceso restringido
dc.subjectdust
dc.subjectpaleoclimate
dc.subjectMATLAB
dc.subjectmodel-data comparison
dc.subjectgridded data
dc.subject.ods02 Zero Hunger
dc.subject.odspa02 Hambre cero
dc.titleGlobal-Krigger: A Global Kriging Interpolation Toolbox With Paleoclimatology Examples
dc.typeartículo
dc.volumen24
sipa.indexWOS
sipa.trazabilidadWOS;2025-01-12
Files