This study investigates the sensitivity of the GNSSseg segmentation method to change in: GNSS data processing method, length of time series (17 to 25 years), auxiliary data used in the integrated water vapor (IWV) conversion, and reference time series used in the segmentation (ERA-Interim versus ERA5). Two GNSS data sets (IGS repro1 and CODE REPRO2015), representative of the first and second IGS reprocessing, were compared. Significant differences were found in the number and positions of detected change-points due to different a priori ZHD models, antenna/radome calibrations, and mapping functions. The more recent models used in the CODE solution have reduced noise and allow the segmentation to detect smaller offsets. Similarly, the more recent reanalysis ERA5 has reduced representativeness errors, improved quality compared to ERA-Interim, and achieves higher sensitivity of the segmentation. Only 45–50% of the detected change-points are similar between the two GNSS data sets or between the two reanalyses, compared to 70–80% when the length of the time series or the auxiliary data are changed. About 35% of the change-points are validated with respect to metadata. The uncertainty in the homogenized trends is estimated to be around 0.01–0.02 kg m−2 year−1.

<p>In recent years, the detection and correction of the non-natural irregularities in the long climatic records, so-called homogenization, has been studied. This work is motivated by the problem of identification of origins of the breakpoints in the segmentation of difference series (difference between a candidate series and a reference series). Several segmentation methods have been developed for the difference series, but many of them assume that the reference series is homogenous. However, the homogeneity of the reference series, in reality, is uncertain and unproven. In our study, we applied the segmentation method GNSSseg (Quarello et al., 2020) on the difference between the Integrated water vapour estimates of the CODE REPRO2015 GNSS data set and the ERA5 reanalysis. About 36.5% of change points can be validated from the GPS metadata, and the origins of the remaining 64.5% are questionable (Nguyen et al., 2021). The ambiguity can be leveraged when there is at least one nearby GPS station with respect to which the candidate series can be compared. The proposed method uses weighted t-tests combining the candidate GPS and ERA series and their homologues (denoted GPS' and ERA') from each nearby station. If sufficient consistency emerges from the six tests for all the nearby stations, a decision can be made whether the breakpoint detected in the candidate GPS-ERA series is due to GPS or, alternatively, to ERA. For each quadruplet (GPS, ERA, GPS', ERA'), six t-tests are performed, and the outcomes are combined. In a set of 81 globally distributed GNSS time series spanning more than 25 years, 56 series have at least one nearby station, where 171 breakpoints are detected in segmentation, in which 136 breakpoints are attributed to the GPS. Among those, 94 breakpoints have consistent results between all the nearby stations. GPS-related breakpoints are used for the correction of the mean shift in the difference series. The impact of the breakpoint correction on the GNSS IWV trend estimates is then evaluated. </p><p>Quarello A, Bock O, & Lebarbier E. (2020). A new segmentation method for the homogenisation of GNSS-derived IWV time-series. <em>arXiv: Methodology</em>.</p><p>Nguyen KN, Quarello A, Bock O, Lebarbier E. Sensitivity of Change-Point Detection and Trend Estimates to GNSS IWV Time Series Properties. <em>Atmosphere</em>. 2021; 12(9):1102. https://doi.org/10.3390/atmos12091102</p>