Journal
JOURNAL OF GEOPHYSICAL RESEARCH-SOLID EARTH
Volume 121, Issue 3, Pages 2054-2068Publisher
AMER GEOPHYSICAL UNION
DOI: 10.1002/2015JB012552
Keywords
geodesy; Theil-Sen; automation; time series; breakdown; uncertainty
Categories
Funding
- NASA ACCESS subaward [S14-NNX14AJ52A-S1]
- NASA Sea Level Rise subaward [1551941]
- NASA ESI grant [NNX12AK26G]
- NSF EarthScope grant [EAR-1252210]
- USGS NEHRP grant [G15AC00078]
- CNES/TOSCA grant
- NASA [43465, NNX12AK26G] Funding Source: Federal RePORTER
- Directorate For Geosciences
- Division Of Earth Sciences [1252210] Funding Source: National Science Foundation
Ask authors/readers for more resources
Automatic estimation of velocities from GPS coordinate time series is becoming required to cope with the exponentially increasing flood of available data, but problems detectable to the human eye are often overlooked. This motivates us to find an automatic and accurate estimator of trend that is resistant to common problems such as step discontinuities, outliers, seasonality, skewness, and heteroscedasticity. Developed here, Median Interannual Difference Adjusted for Skewness (MIDAS) is a variant of the Theil-Sen median trend estimator, for which the ordinary version is the median of slopes v(ij)=(x(j)-x(i))/(t(j)-t(i)) computed between all data pairs i>j. For normally distributed data, Theil-Sen and least squares trend estimates are statistically identical, but unlike least squares, Theil-Sen is resistant to undetected data problems. To mitigate both seasonality and step discontinuities, MIDAS selects data pairs separated by 1year. This condition is relaxed for time series with gaps so that all data are used. Slopes from data pairs spanning a step function produce one-sided outliers that can bias the median. To reduce bias, MIDAS removes outliers and recomputes the median. MIDAS also computes a robust and realistic estimate of trend uncertainty. Statistical tests using GPS data in the rigid North American plate interior show 0.23mm/yr root-mean-square (RMS) accuracy in horizontal velocity. In blind tests using synthetic data, MIDAS velocities have an RMS accuracy of 0.33mm/yr horizontal, 1.1mm/yr up, with a 5th percentile range smaller than all 20 automatic estimators tested. Considering its general nature, MIDAS has the potential for broader application in the geosciences.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available