the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
GPS displacement dataset for study of elastic surface mass variations
Donald Argus
Felix Landerer
David Wiese
Matthias Ellmer
Abstract. Quantification of uncertainty in surface mass change signals derived from GPS measurements poses challenges, especially when dealing with large datasets with continental or global coverage. This study recommends a method for preparing and processing GPS measurements for use in a future joint solution with GRACE(-FO)-like observations. We assess the structure and quantify the uncertainty of vertical land displacement derived from 3045 GPS stations distributed across the continental US. Monthly means of daily positions are available for 15 years. We list the required corrections to isolate surface mass signals in GPS estimates and screen the data using GRACE(-FO) as external validation. Evaluation of GPS timeseries is a critical step, which identifies a) corrections that were missed; b) sites that contain non-elastic signals (e.g., close to aquifers); and c) sites affected by background modelling errors (e.g., errors in the glacial isostatic model). Finally, we quantify uncertainty of GPS vertical land displacement (VLD) estimates through stochastic modeling and quantification of spatially correlated errors. Our aim is to assign weights to GPS estimates of VLD, which will be used in a joint solution with GRACE(-FO). We prescribe white, colored and spatially correlated noise. To quantify spatially correlated noise, we build on the common mode imaging approach adding a geophysical constraint (i.e., surface hydrology) to derive an error estimate for the surface mass signal. We study the uncertainty derived using each technique and find that three techniques exhibit an average noise level between 2–3 mm: white noise, flicker noise, and RMS of residuals about a seasonality and trend fit. Prescribing random walk noise increases the error level such that half of the stations have noise > 4 mm, which is systematic with the noise level derived through modeling of spatial correlated noise. The new data set is suitable for use in a future joint solution with GRACE(-FO)-like observations.
- Preprint
(2446 KB) - Metadata XML
-
Supplement
(185 KB) - BibTeX
- EndNote
Athina Peidou et al.
Status: final response (author comments only)
-
RC1: 'Comment on essd-2023-183', Anonymous Referee #1, 29 Aug 2023
This paper presents some procedures that help in making vertical land displacements useful as complementary data to GRACE for inverting hydrological loads, including noise analysis, using discrepancies between GPS and GRACE to detect missed GPS offsets, and an improvement of the CMC Imaging method of Kreemer and Blewitt (2021). The data produced is ultimately very useful for researchers interested in hydrologic loading. Where will the VLD time-series be available?
Possible issue with Correlation presentation:
The calculation of correlation between GPS and GRACE (per watershed) is not clear to me, given that you have X GPS time-series and Y GRACE mascon time-series. I suspect there is a missing step, that doesn’t seem to be explained: how GRACE data is translated to VLD at the station level. Please explain that. Furthermore, I don’t see how a poor correlation can be attributed to a single station (as the case for St. Lawrence) ; this is assuming that the authors didn’t mean to say that there was a missed GPS offset in all of the time-series, which seems unlikely.
Minor points:
Line 85. Kreemer and Blewitt did not introduce the term CMN, but rather used the term CMC (common mode components), which was first introcuded by Tian and Shen (2016)
Line 118. Originally is not used properly in this contextLine 118. “layout” should read “lay out”
line 126-127: “We overcome CMC’s limitation of include spatially correlated hydrology signals...”. This sentence is grammatically not correct, nor is the context clear: which limitation? That includes both noise and unmodeled signal?
Line 192: I think P_lm is missing in front of “are the associated Legendre polynomials”
Line 220: “𝐴 and 𝐵 being the amplitude and phase”. That is not how formula 2 is written, which is erroneous, should be A*cos(2*pi*t+B)
Line 302-303: “We identified the need for antenna offset corrections (in the case of Great Lakes)”. Before the St. Lawrence watershed was named in this context, now Great Lakes. Later (line 458) Lake Michigan is mentioned. Which is it?
Citation: https://doi.org/10.5194/essd-2023-183-RC1 -
RC2: 'Comment on essd-2023-183', Anonymous Referee #2, 04 Sep 2023
The article is a valuable compilation of methods and approaches used in the geodetic literature to analyze time series. However, it contains several shortcomings that limit its understanding or use of the dataset. The main limitation is that the authors do not make the created dataset and its uncertainties available (as they promised in the abstract and introduction), which greatly prevents me from assessing its potential. Some of the results are presented laconically (as the uncertainties; I presume these are the uncertainties of trend, or…?), and others lack adequate explanations (where is the RW in the data coming from? Perhaps it is the result of inadequate series length? or other effects?). The authors are keen to streamline the procedure for qualifying the data for further analysis, but sometimes speeding up preprocessing or classification can lead to erroneous conclusions.
Major shortcomings include:
- I believe that a key requirement for publishing an article in ESSD is that the proposed dataset be available. Unfortunately, I found nothing about this either in the article or in the supplementary materials. Do the authors plan to make available GPS VLDs that can be taken directly for comparison? Their uncertainties were also described, but are not available. Will these be the uncertainties of individual observations or the uncertainties of trends?
- It would also be useful to make available a list of stations that the authors considered being those responding elastically to the applied load, their coordinates, and a list of stations excluded from comparisons as those responding poroelastically.
- The authors do not explain all the abbreviations used in the article.
- Please show some of the GRACE and GLDAS time series against the GPS time series so that readers can get a general idea of how the time series agree with each other. Not all users of the dataset need to be geodesists.
- It is not clear how the authors converted the TWS available for JPL mascons to displacements. Equation 1 does not describe the entire procedure.
- It is also unclear how the authors interpolated the gap between GRACE and GRACE-FO. This is because they present trend estimates for the period 2018-2021.5, which includes a gap of more than a year in GRACE observations.
- Equation 2 lacks time dependency y(t).
- It is not clear what “ccs” means in the CMN algorithm.
- It is not clear with which noise model the authors perform the analysis for CMC and CMCHF.
- Table 1: What μ means here? Why the authors do not compare the amplitudes of noise or the percentage contribution of different noises to the analyzed combinations? This might be interpretable.
- Figure 6: Is this the uncertainty of trend, or…?
- The authors mentioned a case of unlogged offset. It is presented in the supplementary materials, but only for one case. I think the authors should approach the topic more descriptively and present other cases in which unlogged offsets were also corrected manually.
- Lines 229-230: The authors mention describing interesting cases in Supplements, but they are missing there.
- Several sentences in the text are missing a noun or verb, the authors should carefully review the entire text.
Citation: https://doi.org/10.5194/essd-2023-183-RC2
Athina Peidou et al.
Data sets
GPS displacements and uncertainties to study elastic surface mass variations in North America Athina Peidou, Donald Argus, Matthias Ellmer, Felix Landerer, David Wiese https://zenodo.org/record/8184285
Athina Peidou et al.
Viewed
HTML | XML | Total | Supplement | BibTeX | EndNote | |
---|---|---|---|---|---|---|
298 | 76 | 16 | 390 | 27 | 7 | 10 |
- HTML: 298
- PDF: 76
- XML: 16
- Total: 390
- Supplement: 27
- BibTeX: 7
- EndNote: 10
Viewed (geographical distribution)
Country | # | Views | % |
---|
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1