Bayesian modelling of piecewise trends and discontinuities to improve the estimation of coastal vertical land motion

Abstract

AbstractOne of the major sources of uncertainty affecting vertical land motion (VLM) estimations are discontinuities and trend changes. Trend changes are most commonly caused by seismic deformation, but can also stem from long-term (decadal to multidecadal) surface loading changes or from local origins. Although these issues have been extensively addressed for Global Navigation Satellite System (GNSS) data, there is limited knowledge of how such events can be directly detected and mitigated in VLM, derived from altimetry and tide-gauge differences (SATTG). In this study, we present a novel Bayesian approach to automatically and simultaneously detect such events, together with the statistics commonly estimated to characterize motion signatures. Next to GNSS time series, for the first time, we directly estimate discontinuities and trend changes in VLM data inferred from SATTG. We show that, compared to estimating a single linear trend, accounting for such variable velocities significantly increases the agreement of SATTG with GNSS values (on average by 0.36 mm/year) at 339 globally distributed station pairs. The Bayesian change point detection is applied to 606 SATTG and 381 GNSS time series. Observed VLM, which is identified as linear (i.e. where no significant trend changes are detected), has a substantially higher consistency with large-scale VLM effects of glacial isostatic adjustment (GIA) and contemporary mass redistribution (CMR). The standard deviation of SATTG (and GNSS) trend differences with respect to GIA+CMR trends is by 38% (and 48%) lower for time series with constant velocity compared to variable velocities. Given that in more than a third of the SATTG time series variable velocities are detected, the results underpin the importance to account for such features, in particular to avoid extrapolation biases of coastal VLM and its influence on relative sea-level-change determination. The Bayesian approach uncovers the potential for a better characterization of SATTG VLM changes on much longer periods and is widely applicable to other geophysical time series.