Statistical significance of trends in Zenith Wet Delay from re-processed GPS solutions

Long series of Zenith Wet Delay (ZWD) obtained as part of a homogeneous re-processing of Global Positioning System solutions constitute a reliable set of data to be assimilated into climate models. The correct stochastic properties, i.e. the noise model of these data, have to be identified to assess the real value of ZWD trend uncertainties since assuming an inappropriate noise model may lead to over- or underestimated error bounds leading to statistically insignificant trends. We present the ZWD time series for 1995–2017 for 120 selected globally distributed stations. The deterministic model in the form of a trend and significant seasonal signals were removed prior to the noise analysis. We examined different stochastic models and compared them to widely assumed white noise (WN). A combination of the autoregressive process of first-order plus WN (AR(1) + WN) was proven to be the preferred stochastic representation of the ZWD time series over the generally assumed white-noise-only approach. We found that for 103 out of 120 considered stations, the AR(1) process contributed to the AR(1) + WN model in more than 50% with noise amplitudes between 9 and 68 mm. As soon as the AR(1) + WN model was employed, 43 trend estimates became statistically insignificant, compared to 5 insignificant trend estimates for a white-noise-only model. We also found that the ZWD trend uncertainty may be underestimated by 5–14 times with median value of 8 using the white-noise-only assumption. Therefore, we recommend that AR(1) + WN model is employed before tropospheric trends are to be determined with the greatest reliability.