🤖 AI Summary
This study addresses two key challenges in ecological–climatic transition zone analysis: (1) quantifying spatial boundary localization uncertainty under noisy gridded data, and (2) lacking a rigorous statistical framework for temporal dynamic inference. We propose a unified modeling and inference framework integrating heteroscedastic Gaussian process regression with a scale-normalized maximum absolute deviation (MAD)-based global envelope test (GET). This is the first method to jointly achieve probabilistic reconstruction and formal significance testing of boundary curves. Applied to Köppen–Trewartha climate classification data across the Sahel–Sahara transition zone, empirical results reveal overall boundary stability at the decadal scale (1960–1989), while precisely detecting a localized boundary shift induced by the extreme 1983–1984 drought. The framework combines theoretical rigor—grounded in spatial statistics and functional hypothesis testing—with practical interpretability, offering a novel paradigm for dynamic monitoring of ecological–climatic transitions.
📝 Abstract
Spatial boundaries, such as ecological transitions or climatic regime interfaces, capture steep environmental gradients, and shifts in their structure can signal emerging environmental changes. Quantifying uncertainty in spatial boundary locations and formally testing for temporal shifts remains challenging, especially when boundaries are derived from noisy, gridded environmental data. We present a unified framework that combines heteroskedastic Gaussian process (GP) regression with a scaled Maximum Absolute Difference (MAD) Global Envelope Test (GET) to estimate spatial boundary curves and assess whether they evolve over time. The heteroskedastic GP provides a flexible probabilistic reconstruction of boundary lines, capturing spatially varying mean structure and location specific variability, while the test offers a rigorous hypothesis testing tool for detecting departures from expected boundary behaviors. Simulation studies show that the proposed method achieves the correct size under the null and high power for detecting local boundary shifts. Applying our framework to the Sahel Sahara transition zone, using annual Koppen Trewartha climate classifications from 1960 to 1989, we find no statistically significant decade scale changes in the arid and semi arid or semi arid and non arid interfaces. However, the method successfully identifies localized boundary shifts during the extreme drought years of 1983 and 1984, consistent with climate studies documenting regional anomalies in these interfaces during that period.