Next Article in Journal
Constructing a New Inter-Calibration Method for DMSP-OLS and NPP-VIIRS Nighttime Light
Previous Article in Journal
Modeling and Prediction of Regular Ionospheric Variations and Deterministic Anomalies
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Letter

Seismic Impact of Large Earthquakes on Estimating Global Mean Ocean Mass Change from GRACE

1
School of Communication and Information Engineering, Shanghai University, Shanghai 200444, China
2
Shanghai Astronomical Observatory, Chinese Academy of Sciences, Shanghai 200030, China
3
School of Astronomy and Space Science, University of Chinese Academy of Sciences, Beijing 100049, China
4
Center for Space Research, University of Texas at Austin, TX 78759, USA
5
Shanghai Key Laboratory of Space Navigation and Positing Techniques, Shanghai 200030, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2020, 12(6), 935; https://doi.org/10.3390/rs12060935
Submission received: 10 January 2020 / Revised: 5 March 2020 / Accepted: 11 March 2020 / Published: 13 March 2020
(This article belongs to the Section Remote Sensing in Geology, Geomorphology and Hydrology)

Abstract

:
We analyze the impact of large earthquakes on the estimation of the global mean ocean mass (GMOM) change rate over the 13-year period (January 2003 to December 2015) using the Gravity Recovery and Climate Experiment (GRACE) Release-06 (RL06) monthly gravity solutions released by the Center for Space Research (CSR). We take into account the effects of the December 2004 Mw9.1 and April 2012 Mw8.6 Sumatra earthquakes, the March 2011 Mw9.0 Tohoku-Oki earthquake, and the February 2010 Mw8.8 Chile earthquake. After removing the co- and post-seismic effects of these earthquakes in the oceanic areas by least squares fitting, we estimate the GMOM rate from GRACE monthly observations. Results show that GRACE-observed GMOM rate before the seismic correction is 2.12 ± 0.30 mm/year, while after correction the rate is 2.05 ± 0.30 mm/year. Even though the −0.07 ± 0.02 mm/year seismic influence on GRACE GMOM rate is small on a global scale, it is a systematic bias and should be considered for improved quantification and understanding of the global sea level change.

1. Introduction

Ocean mass variation is caused by water exchange between the ocean and other components of the climate system, including polar ice sheets, mountain glaciers, terrestrial water storage (TWS), and atmospheric water vapor [1]. Due to the small contributions of atmospheric water vapor, most ocean mass change studies focused on contributions of polar ice sheets, mountain glaciers, and TWS change. For example, the Ice Sheet Mass Balance Inter-comparison Exercise (IMBIE) team [2] shows that the contribution of the Antarctica ice sheet melting to the mean sea level rise is 7.6 ± 3.9 mm during the period 1992–2017. In another recent study the IMBIE team [3] indicates that the Greenland ice sheet melting contributes 10.6 ± 0.9 mm to the mean sea level rise for the period 1992–2018. Wouters et al. [4] demonstrates that the contribution of global glaciers excluding Greenland and Antarctica melting to the mean sea level rise is 8.0 ± 1.3 mm during the period 2002–2016. Reager et al. [5] indicates that climate change has resulted in additional storage of TWS, slowing down the ocean mass change by 0.71 ± 0.20 mm/year during the period 2002–2014.
Satellite gravimetry, such as the Gravity Recovery and Climate Experiment (GRACE) and GRACE Follow-On (GFO), provides a unique and large-scale measurement to quantify the global mean ocean mass (GMOM) changes. However, precise quantification of the GMOM change using GRACE/GFO is still challenging (e.g., [6,7]). A number of studies have addressed some of these challenges including the land–ocean leakage errors (e.g., [8,9,10,11]), the significant low-degree term errors (e.g., [12,13,14,15]), and the large glacial isostatic adjustment (GIA) model uncertainties (e.g., [16,17,18,19]). Besides, the gravity effect linked with large offshore earthquakes should also affect the estimation of the GMOM change, and has not been taken into account in previous studies.
During the GRACE mission period (2002–2017) several large earthquakes caused observable gravity change signals, which were detected by GRACE [20,21]. The co- and post-seismic effects have been well investigated with GRACE measurements combined with other observations such as Global Navigation Satellite System (GNSS) displacements, for the events including the 2004 Sumatra Mw9.1 [22], 2011 Tohoku-Oki Mw9.0 [23,24], 2010 Chile Mw8.8 [25], 2005 Nias (Sumatra) Mw8.6 [21,22], 2012 Sumatra Mw8.6 [26], 2007 Bengkulu (Sumatra) Mw8.5 [20,21], 2006 Kuril Islands Mw8.3 [27,28], and 2013 Okhotsk Sea Mw8.3 [29] earthquakes. Among these large earthquakes some have induced noticeable mass change signals, which are clearly captured in the global mass change (expressed in equivalent water height, EWH) trends from GRACE (Figure 1a). Figure 1a shows that the seismic signals within the three regions marked by the red circles are apparent, revealing clear positive–negative spatial patterns corresponding to mass increase and loss induced by the earthquakes. Figure 1b shows the example of the GRACE mass change time series at the point A (2.5°N, 93.5°E) from January 2003 to December 2015, which is located at the oceanic area within the Sumatra earthquake region. The time series exhibit a clear jump in December 2004 (i.e., the occurrence month of the megathrust earthquake) as well as gradual increases thereafter, which are clearly caused by co- and post-seismic deformations, rather than sea water mass changes. Therefore, we need to correct these seismic effects to avoid potential systematic bias in GRACE GMOM rate estimation.
The main purpose of this study is to estimate the effect of large offshore earthquakes on the GMOM change rate 2003–2015 observed by GRACE. The earthquakes we select include the December 2004 and April 2012 Sumatra earthquakes, March 2011 Tohoku-Oki earthquake, and February 2010 Chile earthquake, whose effects are visible in Figure 1a. We correct the seismic effect using least squares fitting considering both co- and post-seismic terms. Then we compare the GMOM change rates before and after the correction, in order to show how significantly these large earthquakes affect the determination of GRACE GMOM rate.

2. Data and Method

We used the Center for Space Research (CSR) Release-06 (RL06) monthly gravity field solutions to estimate the ocean mass change rate during the period January 2003 to December 2015. The solutions contain spherical harmonic (SH) coefficients up to degree and order 60 [30]. We applied the 300 km Gaussian smoothing to suppress the spatial noise of the GRACE observations [31]. Solid earth deformation effect due to the glacial isostatic adjustment (GIA) was removed using the model from A et al. [16], thereafter referred as the A13 GIA model. We replace the GRACE C20 terms with the Satellite Laser Ranging (SLR) observations [32]. Since the degree-1 terms, representing the geocenter motion were absent in the GRACE monthly solutions, we used the degree-1 coefficients provided by the GRACE project (GRACE Technical Note 13) [33] to include the geocenter motion effect, which are derived from GRACE observations using an improved method [12,34]. Then we calculated the global gridded (1° × 1°) surface mass change fields (with respect to the mean field 2003.01–2015.12) from the 142 monthly solutions. By least squares fitting of the time series at each grid point of these monthly solutions, we obtained the global mass change rate distribution.
During the least squares fitting, the seismic effects of the three regions (as circled by red curves in Figure 1) need to be corrected. The correction is carried out by removing the co- and post-seismic terms with least squares fitting of the GRACE observation time series (e.g., [35,36]). When the time series at the point includes the seismic effect of only one earthquake, the fitting equation can be expressed as:
y s e i s m i c _ e f f e c t ( t ) = { C 1 t < t e q C 2 + p · ( 1 e ( t t e q ) / τ ) t t e q
where C1 and C2 are the constant terms before and after the earthquake. p is the total post-seismic gravity change reached at the end of relaxation. τ is the e-folding relaxation time of the post-seismic changes. Here we chose τ = 7 months in the Sumatra and Chile earthquake regions fitting. We also chose τ = 4 months in the Tohoku-Oki earthquake region fitting (the selection of the parameter τ on the time series fitting will be discussed later in Section 4). teq is the earthquake occurrence time. Figure 2a shows an example of the time series fitting at a selected point B (36.5°N, 145.5°E) in the Tohoku-Oki earthquake region. When the time series at the point includes the seismic effects of two earthquakes, the fitting equation can be expressed as:
y s e i s m i c _ e f f e c t ( t ) = { C 1 t t e q 1 C 2 + p 1 · ( 1 e ( t t e q 1 ) / τ 1 ) t e q 1 < t t e q 2 C 3 + p 2 · ( 1 e ( t t e q 2 ) / τ 2 ) t > t e q 2
where C1 is the pre-seismic constant of the first earthquake, C2 is the post-seismic constant of the first earthquake, and C3 is the post-seismic constant of the second earthquake. teq1 and teq2 represent the occurrence times of the first and second earthquakes, respectively. τ1 and τ2 are the e-folding relaxation time parameters of the post-seismic change corresponding to the first and second earthquakes, respectively. Note here it is implicitly assumed that the post-seismic deformation of the first event is fully relaxed before the second event occurs. This assumption should make sense for the case of Sumatra earthquake region because the interval between the two events is more than 7 years, during which most post-seismic deformation of the 2004 Sumatra earthquake has been relaxed. Figure 2b shows an example of the time series fitting at a selected point A (2.5°N, 93.5°E) in the Sumatra earthquake region.
Then we removed the seismic effect from the GRACE observation and obtained the residual time series (GRACE observation minus the co- and post-seismic terms in Equation (1) or Equation (2)) for each grid point over the region. The corrected mass change rate can be figured out from the residual time series with least squares fitting, expressed as:
y r e s i d u a l ( t ) = A a n c · c o s 2 π t + A a n s · s i n 2 π t + A s e c · c o s 4 π t + A s e s · s i n 4 π t + A t r · t + A 0
where A a n c ,   A a n s ,   A s e c   and   A s e s represent the periodic terms including the annual and semi-annual variations. A t r represents the corrected linear trend and A 0 is the constant term. By removing the seismic terms from the GRACE time series, we corrected the biases from the impact of seismic deformation.
After correcting the seismic effect, we used the global EWH rates to calculate the mean ocean mass change rate. During the calculation, we need to correct the signal leakage effect due to limited spatial resolution of GRACE observations. Since the GRACE monthly solutions are truncated to the SH degree 60 and the spatial smoothing is also applied to suppress noise, significant mass-change signals leak from land to ocean in the coastal regions (e.g., the Gulf of Alaska, and the surrounding regions of Greenland and Antarctica in Figure 1a). In this study we used the forward modeling (FM) to correct the leakage in the continent-ocean boundary areas. The FM method uses iterations to correct the signal leakage between land and ocean [37]. We calculated the GMOM change rates before and after the seismic correction, and evaluated the bias due to the effects of the large earthquakes.

3. Results

3.1. Seismic Effect Correction

First, we corrected the seismic effects in the Sumatra earthquake region. Since we considered two earthquakes in this region (the December 2004 Mw9.1 and April 2012 Mw8.6 earthquakes), we adopted Equation (2) (see Section 2) to correct the seismic effect. We noted that the seismic signals of the March 2005 MW 8.6 Nias earthquake were also included in GRACE observations over this region. Nevertheless, its occurrence time was close to the December 2004, so that we implicitly treated the two earthquakes as one event in the time series fitting. The correction region ranged from −13°N to 17°N in latitude and from 81°E to 111°E in longitude. The results are shown in Figure 3. The corrected spatial distribution (Figure 3b) shows that most seismic effects were removed in this region, and the amplitude of the rate from residual variations was basically consistent with that from the surrounding oceanic areas. Figure 3c,d shows the examples of time series fitting at two selected points D (3.5°N, 94.5°E) and E (10.5°N, 98.5°E), locating in the positive and negative mass-change areas, respectively. The linear rates at points D and E were 3.58 cm/year and −1.97 cm/year before the correction, and were 0.28 cm/year and −0.19 cm/year after the correction, respectively. Both examples show significant changes before and after the correction, since the rate values were dramatically reduced. After seismic correction, the linear rates (0.28 cm/year and −0.19 cm/year) were comparable to the amplitudes of the rates at the other grid points in the surrounding oceanic areas, in which all the grid values were within the range −0.3 to +0.3 cm/year. Figure 3b shows that the corrected linear rates at the grid points surrounding the epicenter were consistently small (as those at the grid points in further oceanic areas), suggesting that the seismic impact was well corrected.
Second, we corrected the seismic effects in the Tohoku-Oki earthquake region. We adopted Equation (1) (see Section 2) to correct the seismic effect in this region. The correction region ranges from 30°N to 42°N in latitude and from 132°E to 152°E in longitude. Figure 4a,b shows the EWH rate spatial distributions before and after the seismic correction. Similar to that in the Sumatra earthquake region, the time series fitting examples at both the points F (36.5°N, 144.5°E) and G (39.5°N, 138.5°E) in Figure 4c,d, also shows significant changes before and after seismic correction (1.82 vs. 0.09 cm/year for the point F, and −1.59 vs. 0.09 cm/year for the point G). The corrected linear rates surrounding the epicenter were small (within the range −0.1 to +0.1 cm/year, Figure 4b), which indicates that the seismic impact was well corrected for the Tohoku-Oki earthquake region. It should be noted that a discontinuity feature was visible at the correction region boundaries, which will be discussed later in Section 4.2.
Finally, we corrected the seismic effects in the Chile earthquake region. Similar to that in the Tohoku-Oki earthquake region, we used Equation (1) (see Section 2) to correct the seismic effect in this region. The correction region ranged from −41°N to −31°N in latitude and from 277°E to 297°E in longitude. Results show that the changes before and after seismic correction were noticeable, not only in the EWH rate distribution (Figure 5a vs. Figure 5b), but also in the time series fitting examples at the points H (−37.5°N, 283.5°E) and I (−35.5°N, 291.5°E; Figure 5c vs. Figure 5d). The fitted linear rates at points H and I were 0.89 and −1.78 cm/year before the correction, and 0.09 and 0.01 cm/year after the correction, respectively. In the areas surrounding the epicenter, the corrected linear rates (Figure 5b) were within the range −0.1 to +0.1 cm/year, indicating a good seismic impact correction for the Chile earthquake region.

3.2. Impact of Seismic Effect on the GMOM Change Rate

After correcting the seismic effects in the three regions, we calculated the impact of seismic effect on GRACE GMOM rate. We first estimated the GRACE GMOM rate without seismic effect correction. Based on FM with 100 iterations, the recovered GMOM change rates turn out to be convergent (see the convergence curves with respect to the iterations in Figure 6). Then we took the convergent value at 100 iterations as the leakage-corrected GMOM change rate, which was 2.12 mm/year (shown by the red curve in Figure 6).
We then estimated the GRACE GMOM rate with seismic effect corrections in three cases. When only correcting the seismic effect in the Sumatra earthquake region, the estimated GRACE GMOM rate was 2.08 mm/year (see the blue curve in Figure 6). The spatial distribution of the global EWH rate after correcting the seismic effect in Sumatra region is shown in Figure 7a. When correcting the seismic effects in the Sumatra and Tohoku-Oki earthquake regions, the GMOM rate was 2.07 mm/year (see the magenta curve in Figure 6). The global EWH rate distribution after correcting the seismic effects in these two regions is shown in Figure 7b. When we corrected the seismic effects in all the three regions (considering the Sumatra, Tohoku-Oki and Chile earthquakes), the estimated GMOM rate was 2.05 mm/year (see the black curve in Figure 6). The seismic-corrected global EWH rate distribution is shown in Figure 7c. The corresponding GMOM change rates are summarized in Table 1.
As a result, the influence of these four large earthquakes on the estimation of GMOM change rate was −0.07 mm/year. Among them, the bias caused by the Sumatra earthquake region (including the December 2004 and April 2012 Sumatra earthquakes) accounted for −0.04 mm/year. The bias caused by the March 2011 Tohoku-Oki earthquake was −0.01 mm/year. Additionally, the bias caused by the February 2010 Chile earthquake accounted for −0.02 mm/year.

4. Discussion

4.1. Uncertainty Estimate of GMOM Change Rate

To estimate the uncertainty of the GMOM change rate obtained in this study, we followed the method used in Chen et al. [9,14] based on least squares fitting. We applied the forward modeling (FM) leakage correction (with 100 iterations, as done for the global trend rate field in Section 3.2) to the monthly global EWH change fields (142 months in total, as stated in Section 2), and obtained the GMOM change time series from January 2003 to December 2015. Then we fitted the time series with least squares by taking account of the constant term, linear trend, annual, and semi-annual terms (see Equation (3) in Section 2). We took the 2 sigma of the fitted trend term as the uncertainty of the GMOM rate, which was ± 0.11 mm/year (see Line 2, Column 2 in Table 2). In addition, by correcting the seismic effects in the three earthquake regions (with Equations (1) and (2)) from the monthly global EWH change fields and then applying the FM leakage correction, we obtained the seismic-corrected GMOM change time series. The corresponding fitted GMOM rate with seismic correction and its uncertainty are listed in Table 2 (Line 2, Column 3).
We also calculated the GMOM rates with the Geoforschungszentrum (GFZ) RL06 [38] and Jet Propulsion Laboratory (JPL) RL06 [39] GRACE monthly solutions from January 2003 to December 2015. All the data processing steps were the same as those done for the CSR RL06 solutions, except that the degree-1 coefficients corresponding to the GFZ and JPL solutions [32] were used, respectively. The GFZ and JPL GMOM rates without and with seismic correction are listed in Table 2 (Lines 3 and 4), together with their uncertainties. The differences between each other among the GMOM rates from CSR, GFZ, and JPL solutions reached up to 0.11 mm/year (i.e., the difference between the GFZ and JPL rates before seismic correction), which was right within the estimated uncertainty range. Our estimated uncertainty ± 0.11 mm/year was close to the one provided in Chen et al. [14], in which they obtained a ± 0.14 mm/year uncertainty level (see their Table 1) by fitting the GMOM time series from January 2005 to December 2015. Our uncertainty estimate was slightly smaller because we used a longer time span (2003–2015) of GRACE data in the time series fitting.
It should be noted that the uncertainty estimated from the formal errors by least squares fitting might only represent the contribution from observation noise of the time series. It tends to underestimate the true uncertainty level since other impacts also contribute to the GMOM rate errors. One impact is from the uncertainty of the spatial filtering applied to the GRACE solutions. A recent study in Chen et al. [9] estimates this impact and indicates that different spatial filtering choices on GRACE data may lead to a difference up to 0.2 mm/year in the GMOM rate estimation (see their comparisons among the three cases of applying decorrelation and 300 km Gaussian smoothing, decorrelation only and no filters, in Table 5). Here we took ± 0.20 mm/year as the resulting uncertainty for the GMOM rate from the variety of different filtering choices.
Another impact results from the uncertainty of GIA models used to correct the GIA effect in GRACE observations. Currently there are several GIA models available for the correction in the GMOM rate estimation. Different GIA models may lead to notable differences in the correction. Chen et al. (2018, Table 6) [9] compares the GIA contributions to the GMOM rate from three commonly used models, including the Paulson07 [40], A13 (as used in this study), and ICE6G [41] GIA models. The comparison indicates the GIA contribution difference reaches up to 0.2 mm/year. Thus we took ±0.20 mm/year as the uncertainty corresponding to GIA model errors for the GMOM rate estimation in our study.
Considering the error sources from the observation noise (±0.11 mm/year), the spatial filtering uncertainty (±0.20 mm/year) and the GIA model uncertainty (±0.20 mm/year), we took the combined error estimate ±0.30 mm/year (from the square root of the sum of squares for the three) as a more reasonable uncertainty level for the GMOM rate in our study. As a result, the CSR RL06 GMOM rates should be 2.12 ± 0.30 mm/year and 2.05 ± 0.30 mm/year before and after seismic effect correction, respectively. Here we note that the GRACE’s low-degree term errors (e.g., the degree-1 and C20 terms) might also have notable contribution to the GMOM rate uncertainty. Therefore, the true uncertainty level could be even larger.

4.2. Uncertainty Estimate of Seismic Impact on the GMOM Change Rate

For the seismic effect correction, theoretically we could use predictions from earthquake dislocation models to remove the co- and post-seismic terms in the GRACE observations. However, the uncertainty of current fault slip models and Earth models necessary for the model predictions are fairly large, which leads to large errors for the seismic modeling. Different model inputs may lead to distinctly different co- and post-seismic predictions even for the same earthquake. Therefore, in this study we estimated the influence of seismic terms on the GMOM change by least squares fitting of GRACE observations, rather than using the dislocation model predictions.
Nevertheless, the time series fitting method has its own uncertainty on the estimation of seismic impact on the GMOM change rate. The choice of the post-seismic relaxation time parameter τ affects the fitting results, which inevitably leads to errors for the evaluation of the GMOM change rate. The value of τ varies for earthquakes in different regions. For example, de Linage et al. [35] shows that with a non-linear optimization method, the range of the estimated post-seismic parameter τ for the December 2004 Mw9.1 earthquake is 0.5–0.75 year (i.e., 6–9 months) in the region surrounding the epicenter. Li et al. [42] applies a residual minimization standard of the observations and indicates that the value of post-seismic parameter τ for the March 2011 Mw9.0 Tohoku-Oki earthquake is around 0.32 year (i.e., 4 months). Therefore, in this study we used τ = 7 and τ = 4 months to fit the post-seismic terms for the Sumatra and Tohoku-Oki earthquake regions, respectively. Since the post-seismic signal was relatively weaker for the February 2010 Mw8.8 Chile earthquake (compared to the co-seismic jumps and the fluctuations related with other effects and observation noise, see Figure 5c,d), the value of parameter τ should only have a fairly slight impact on the fitting result of post-seismic terms. Here, we empirically chose τ = 7 months to correct the seismic effect for the Chile earthquake region.
We approximately estimated the impact of fitting parameter τ on our GMOM rate by comparing the results with different τ values. Considering that the contribution from the Sumatra earthquake region was dominant among the seismic effect corrections (see Table 1), we took the Sumatra case as an example. We did the time series fitting with the values τ = 6, 7, 8, and 9 months and calculated the GMOM rates for the four cases (following the processing steps in Section 3.2). We found the maximum change from τ = 6 to 9 months was 0.0002 mm/year for the contribution to the GMOM rate. Corresponding results are listed in Table S1 of the Supplementary Material. Thus it should be reasonable to infer that the τ value choice influence on the GMOM rate in this study was at a level of 0.001 mm/year, and should be far less than the range of ± 0.01 mm/year.
The range of the spatial window selected to remove the seismic effects also contributes to the uncertainty of estimating the seismic impact on the GMOM rate. Currently we directly selected a rectangular area surrounding the epicenter as the correction region. However, in fact, an ideal choice would be choosing the correction area based on the spatial characteristics of seismic signals in the near field of the earthquake. The simple rectangular treatment leads to a visible discontinuity feature in the spatial distribution after the seismic correction in the Sumatra, Tohoku-Oki, and Chile earthquake regions (see Figure 3b, Figure 4b, and Figure 5b in Section 3.1). Here we varied the spatial window ranges in the estimation to investigate its influence on the results. Taking also the Sumatra earthquake region as an example, we compared the corrected GMOM rates in three cases with the ranges of 20° × 20°, 30° × 30° and 40° × 40° (see Table S2 in the Supplementary Material). The differences between each other among the three cases were within 0.01 mm/year.
Another factor to impact the estimated result is the limited number of selected earthquakes during the considered time period (2003–2015). In this study our consideration for seismic correction included five events, namely, the Sumatra 2004, 2012, and Nias 2005 earthquakes, the Tohoku-Oki 2011 earthquake, and the Chile 2010 earthquake. During the period 2003–2015 several other large earthquakes (not included in this study) also caused gravity change effects observable by GRACE, such as the Bengkulu 2007 Mw8.5, Kuril Islands 2006 Mw8.3, and Okhotsk Sea 2013 Mw8.3 events. However, for the linear mass-change rates over the period 2003–2015, their contributions are small and could be treated as not significant. This can be verified by examining the global EWH trend rates in Figure 1a. The positive-negative patterns associated with seismic effects were noticeable only in three regions of the global trend map (as marked by the red circles), while the effects of the above excluded earthquake were not visible. Thus we could infer that their corresponding contributions to the GMOM rate (2003–2015) were well below that from any of the three cases included in this study (i.e., the Sumatra, Tohoku-Oki, and Chile earthquake regions, among which the smallest contribution was 0.01 mm/year, see Table 1).
The GRACE monthly solutions from different releasing centers may also lead to discrepancy of the results. According to Table 2, the discrepancies reached up to 0.02 mm/year among the estimates from CSR RL6, GFZ RL06, and JPL RL06 solutions (with −0.07, −0.08, and −0.09 mm/year for CSR, GFZ, and JPL, respectively). The 0.02 mm/year uncertainty related with different data processing centers was well above those contributions from the error sources discussed above. Therefore, we took ± 0.02 mm/year as a relatively reasonable uncertainty estimate for the seismic impact on the GMOM rate 2003–2015 in our study, which came out as −0.07 ± 0.02 mm/year.

4.3. Comparison with Estimates from Altimery and Argo

In addition to the GRACE estimation, the GMOM change rate could also be obtained from the Altimetry and Argo measurements (by removing the Argo steric contribution from the Altimetry sea level change observation). Due to the spatial coverage limitation of the Argo floats during its early observation period, the sea level change under a global scale is usually investigated with Argo data after 2005. Recently, Chen et al. [14] obtained a GMOM rate of 2.68 ± 0.15 mm/year for the period 2005.01–2015.12, with 3.79 mm/year (from Altimetry) minus 1.11 mm/year (from Argo). In an earlier study Dieng et al. [10] estimated the GMOM rate (2004.01–2015.12) as 2.35 ± 0.14 mm/year, with 3.49 mm/year (Altimetry) minus 1.14 mm/year (Argo). By comparing the above results one may find that the GMOM rate differs notably for different time spans, which is likely because the data period currently available is not long enough to precisely quantify the long-term trend rate. Moreover, the long-term trend quantification is also contaminated by the strong interannual variations that are difficult to efficiently remove in the current observations of limited data length.
For the period 2003.01–2015.12 considered in this study, the GMOM rate by Altimetry minus Argo is difficult to obtain reliably due to the Argo data availability problem over the earlier years of the period. If one assumes that the GMOM rate 2003.01–2015.12 by Altimetry minus Argo is close to that for the period 2004.01–2015.12, neither of our GRACE estimated GMOM rates before and after seismic correction (2.12 ± 0.30 mm/year and 2.05 ± 0.30 mm/year) agrees with the trend rate from Altimetry and Argo (2.35 ± 0.14 mm/year). The GRACE estimate of GMOM rate is noticeably smaller than that obtained from Altimetry and Argo, which is the same situation as in the previous studies Chen et al. [14] and Dieng et al. [10] (mentioned above). The discrepancy might attribute to the low-degree term errors as well as data filtering uncertainties associated with GRACE observations, model uncertainties in the GIA effect correction, Argo measurement uncertainties (especially for the deep ocean), etc. However, when considering the large uncertainty range of the GMOM rate, our GRACE result after seismic correction (2.05 ± 0.30 mm/year) can overlap with the Altimetry–Argo result (2.35 ± 0.14 mm/year) within the error range.
In an earlier study in 2019, Uebbing et al. [19] estimated the GMOM change rate with GRACE observations, and obtained the rate 1.43 mm/year for the period 2002.06–2014.05 and 1.75 mm/year for the period 2002.06–2016.05. Even though the period 2002.06–2016.05 is close to the one considered in this study (2003.01–2015.12), their estimated GMOM rate is distinctly smaller than ours. Considering that we used the same A13 model for the GIA effect correction as used in Uebbing et al. [19], the discrepancy might result from the reasons as follows. Besides the data period difference that may contribute to the discrepancy, the GRACE data version difference may contribute notably to the discrepancy. In the calculation we used the RL06 solutions, while Uebbing et al. [19] use the RL05 solutions. Moreover, the applied leakage correction methods are different (i.e., the FM method in our study vs. the buffer zone method in Uebbing et al. [19]). Note that the FM correction method usually leads to a relatively larger GMOM rate (up to 0.2–0.3 mm/yr) than the buffer zone method, as indicated by Chen et al. [9,14].

4.4. Seismic Impact for Shorter Time Spans and Regional Scales

For the GMOM change rate estimation, the seismic correction may increase to a large extent when choosing a shorter time period that covers the occurrence time of the earthquake. As an example, we selected the period from January 2003 to December 2007, during which the contribution of the 2004 Sumatra earthquake was apparently dominant. The difference between the 2003–2007 GMOM rates before and after seismic correction reached up to −0.25 mm/year (see Figure S1 in the Supplementary Material), which indicates the seismic effect was significantly larger than the −0.07 mm/year correction for the 2003–2015 period.
The spatial scale of the study area also affects the results. For a global scale, the seismic effect was −0.07 mm/year (2003–2015), while for regional scales, the effect might be significantly different. For instance taking the Indian Ocean as the study area, the mass change signals in this region caused by the 2004 and 2012 Sumatra earthquakes will affect the regional ocean the mass rate more dramatically. As a result, the corresponding seismic effect on the Indian Ocean mean mass change rate could be much larger than that on the GMOM rate.

5. Conclusions

We studied the impact of the seismic effect caused by the December 2004, April 2012 Sumatra earthquakes, the March 2011 Tohoku-Oki earthquake, and the February 2010 Chile earthquake on the GMOM rate estimated by GRACE. We used least squares to fit the co- and post-seismic terms, and applied the FM recovery method to correct the land-ocean leakage effect. Results show that GRACE GMOM rates before and after seismic correction were 2.12 ± 0.30 mm/year and 2.05 ± 0.30 mm/year. The uncertainty level ± 0.30 mm/year was estimated by considering the contributions from the GRACE observation noise, the uncertainty due to various spatial filtering choices and the GIA model errors. The bias caused by the four large earthquakes on GMOM rate was −0.07 ± 0.02 mm/year, with the uncertainty estimated by investigating various error sources. Although the seismic effect on the GMOM rate turned out under the uncertainty level of the current GRACE GMOM rate estimation at a global scale (−0.07 vs. ± 0.30 mm/year), it is a systematic bias. It should be considered for improved quantification and understanding of the global sea level change. For regional scales, this bias could be more significant. Moreover, the seismic effect would affect the GMOM rate more dramatically when focusing on shorter time spans. Consequently, it is necessary to consider the seismic effect correction for more precise estimation of GMOM change by satellite gravimetry observations from GRACE, GRACE Follow-On, and especially for the future satellite gravimetry missions with expected higher resolutions and precisions.

Supplementary Materials

The following are available online at https://www.mdpi.com/2072-4292/12/6/935/s1. Table S1: The 2003.01–2015.12 GMOM rates with seismic effect over the Sumatra earthquake region corrected using different τ values in the time series fitting, Table S2: The 2003.01–2015.12 GMOM rates with seismic effect over the Sumatra earthquake region corrected by considering different spatial window ranges, Figure S1: Seismic correction of the GMOM change rate for a shorter time span 2003–2007.

Author Contributions

Conceptualization, J.L.; Methodology, L.T., J.L. and J.C.; Software, L.T. and J.L.; Validation, S.-Y.W., R.W., and X.H.; Writing—Original Draft Preparation, L.T. and J.L.; Writing—Review and Editing, J.L., J.C. and S.-Y.W. J.L. conceived the ideas of the research and wrote the Discussion and Conclusions sections of the manuscript. L.T. did the calculations and wrote the Method and Results sections of the manuscript. J.C. and S.-Y.W. contributed to the preparation of the manuscript through review and comments. R.W., and X.H. contributed to the validation of the methods and results of the research. All authors contributed to the discussion of the results, as well as the revisions of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key Research and Development Program of China [Grant No. 2016YFB0501405], the Natural Science Foundation of China [Grant No. 11873075], the Natural Science Foundation of Shanghai [Grant Nos. 19ZR1466900, 17ZR1435600], and the Open Fund of Stay Key Laboratory of Geodesy and Earth’s Dynamics [Grant No. SKLGED2018-2-2-E]. The APC was funded by [Grant No. 2016YFB0501405] and [Grant No. 11873075].

Acknowledgments

We thank the three anomalous reviewers for providing constructive comments and suggestions to improve our work. The GRACE CSR, GFZ, and JPL RL06 monthly solutions used in this study are available through GFZ’s ftp link (ftp://rz-vm152.gfz-potsdam.de/grace/Level-2/). The authors L.T., J.L., S.Y.W., and X.H. are grateful to the membership support from the Laboratory of Astronomy and Space Technology Applications at Shanghai Astronomical Observatory, Chinese Academy of Sciences. This work made use of the High Performance Computing Resource in the Core Facility for Advanced Research Computing at Shanghai Astronomical Observatory, Chinese Academy of Sciences.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. WCRP Global Sea Level Budget Group. Global sea-level budget 1993-present. Earth Syst. Sci. Data 2018, 10, 1551–1590. [Google Scholar] [CrossRef] [Green Version]
  2. The IMBIE Team. Mass balance of Greenland Ice Sheet from 1992 to 2018. Nature 2019. [Google Scholar] [CrossRef]
  3. The IMBIE Team. Mass balance of Antarctica Ice Sheet from 1992 to 2017. Nature 2019. [Google Scholar] [CrossRef] [Green Version]
  4. Wouters, B.; Gardner, A.S.; Moholdt, G. Global Glacier Mass Loss during the GRACE satellite Mission (2002–2016). Front. Earth Sci. 2019, 7, 96. [Google Scholar] [CrossRef] [Green Version]
  5. Reager, J.T.; Gardner, A.S.; Famiglietti, J.S.; Wiese, D.N.; Eicker, A.; Lo, M.-H. A decade of sea level rise slowed by climate-driven hydrology. Science 2016, 351, 699–703. [Google Scholar] [CrossRef] [Green Version]
  6. Dieng, H.B.; Cazenave, A.; von Schuckmann, K.; Ablain, M.; Meyssignac, B. Sea level budget over 2005–2013: Missing contributions and data errors. Ocean Science 2015, 11, 789–802. [Google Scholar] [CrossRef] [Green Version]
  7. Cazenave, A.; Palanisamy, H.; Ablain, M. Contemporary sea level changes from satellite altimetry: What have we learned? What are the new challenges? Adv. Space Res. 2018, 62, 1639–1653. [Google Scholar] [CrossRef]
  8. Chen, J.L.; Wilson, C.R.; Tapley, B.D. Contribution of ice sheet and mountain glacier melt to recent sea level rise. Nat. Geosci. 2013, 6, 549–552. [Google Scholar] [CrossRef]
  9. Chen, J.L.; Tapley, B.D.; Save, H.; Tamisiea, M.E.; Bettadpur, S.; Ries, J. Quantification of ocean mass change using gravity recovery and climate experiment, satellite altimeter, and Argo floats observations. J. Geophys. Res. Solid Earth 2018, 123, 10221–10225. [Google Scholar] [CrossRef]
  10. Dieng, H.B.; Cazenave, A.; Meyssignac, B.; Ablain, M. New estimate of the current rate of sea level rise from a sea level budget approach. Geophys. Res. Lett. 2017, 44, 3744–3751. [Google Scholar] [CrossRef]
  11. Xi, H.; Zhang, Z.; Lu, Y.; Li, Y. Mass sea level variation in the South China Sea from GRACE, altimetry and model and the connection with ENSO. Adv. Space Res. 2019, 64, 117–128. [Google Scholar] [CrossRef]
  12. Sun, Y.; Riva, R.; Ditmar, P. Optimizing estimates of annual variations and trends in geocenter motion and J2 from a combination of GRACE data and geophysical models. J. Geophys. Res. Solid Earth 2016, 121, 8352–8370. [Google Scholar] [CrossRef] [Green Version]
  13. Jeon, T.; Seo, K.W.; Youm, K.; Chen, J.L.; Wilson, C.R. Global sea level change signatures observed by GRACE satellite gravimetry. Sci. Rep. 2018, 8, 13519. [Google Scholar] [CrossRef] [PubMed]
  14. Chen, J.L.; Tapley, B.D.; Seo, K.-W.; Wilson, C.R.; Ries, J. Improved Quantification of Global Mean Ocean Mass Change Using GRACE Satellite Gravimetry Measurements. Geophys. Res. Lett. 2019, 46, 13984–13991. [Google Scholar] [CrossRef]
  15. Kim, J.-S.; Seo, K.-W.; Jeon, T.; Chen, J.L.; Wilson, C.R. Missing hydrological contribution to sea level rise. Geophys. Res. Lett. 2019, 46. [Google Scholar] [CrossRef] [Green Version]
  16. Wahr, J.; Zhong, S. Computations of the viscoelastic response of a 3-D compressible Earth to surface loading: An application to Glacial Isostatic Adjustment in Antarctica and Canada. Geophys. J. Int. 2013, 192, 557–572. [Google Scholar] [CrossRef]
  17. Peltier, R.W.; Argus, D.F.; Drummond, R. Comment on “An assessment of the ICE-6G_C (VM5a) glacial isostatic adjustment model” by Purcell et al. J. Geophys. Res. Solid Earth 2018, 123, 2019–2028. [Google Scholar] [CrossRef]
  18. Caron, L.; Ivins, E.R.; Larour, E.; Adhikari, S.; Nilsson, J.; Blewitt, G. GIA Model Statistics for GRACE Hydrology, Cryosphere, and Ocean Science. Geophys. Res. Lett. 2018, 45, 2203–2212. [Google Scholar] [CrossRef]
  19. Uebbing, B.; Kushce, J.; Rietbroek, R.; Landerer, F.W. Processing Choices Affect Ocean Mass Estimates from GRACE. J. Geophys. Res. Oceans 2019, 124, 1029–2044. [Google Scholar] [CrossRef]
  20. Han, S.C.; Riva, R.; Sauber, J.; Okal, E. Source parameter inversion for recent great earthquakes from a decade-long observation of global gravity fields. J. Geophys. Res. Solid Earth 2013, 118, 1240–1267. [Google Scholar] [CrossRef] [Green Version]
  21. Chao, B.F.; Liau, J.R. Gravity changes due to large earthquakes detected in GRACE satellite data via empirical orthogonal function analysis. J. Geophys. Res. Solid Earth 2019, 124, 3024–3035. [Google Scholar] [CrossRef]
  22. Han, S.C.; Shum, C.K.; Bevis, M.; Ji, C.; Kuo, C. Crustal dilatation observed by GRACE after the 2004 Sumatra-Andaman earthquake. Science 2006, 313, 658–661. [Google Scholar] [CrossRef] [Green Version]
  23. Han, S.C.; Sauber, J.; Riva, R. Contribution of satellite gravimetry to understanding seismic source processes of the 2011 Tohoku-Oki earthquake. Geophys. Res. Lett. 2011, 38, L24312. [Google Scholar] [CrossRef]
  24. Matsuo, K.; Heki, K. Coseismic gravity changes of the 2011 Tohoku-Oki earthquake from satellite gravimetry. Geophys. Res. Lett. 2011, 38, L00G12. [Google Scholar] [CrossRef] [Green Version]
  25. Han, S.C.; Sauber, J.; Luthcke, S. Regional gravity decrease after the 2010 Maule (Chile) earthquake indicates large-scale mass redistribution. Geophys. Res. Lett. 2010, 37, L23307. [Google Scholar] [CrossRef]
  26. Han, S.C.; Sauber, J.; Pollitz, F. Coseismic compression/dilatation and viscoelastic uplift/subsidence following the 2012 Indian Ocean earthquakes quantified from satellite gravity observations. Geophys. Res. Lett. 2015, 42. [Google Scholar] [CrossRef]
  27. Han, S.C.; Sauber, J.; Pollitz, F. Postseismic gravity change after the 2006–2007 great earthquake doublet and constraints on the asthenosphere structure in the central Kuril Islands. Geophys. Res. Lett. 2016, 43, 3169–3177. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Rahimi, A.; Li, J.; Naeeni, M.R.; Shahrisvand, M.; Fatolazadeh, F. On the extraction of co-seismic signal for the Kuril Island earthquakes using GRACE observations. Geophys. J. Int. 2018, 215, 346–362. [Google Scholar] [CrossRef]
  29. Tanaka, Y.-S.; Heki, K.; Matsuo, K.; Shestakov, N.V. Crustal subsidence observed by GRACE after the 2013 Okhotsk deep-focus earthquake. Geophys. Res. Lett. 2015, 42, 3204–3209. [Google Scholar] [CrossRef] [Green Version]
  30. Bettadpur, S. CSR Level-2 Processing Standards Document for Product Release 06, GRACE 327-742, Revision 5.0, The GRACE Project, Center for Space Research, University of Texas at Austin. Available online: https://podaac-tools.jpl.nasa.gov/drive/files/allData/grace/docs/L2-CSR006_ProcStd_v5.0.pdf (accessed on 18 April 2018).
  31. Wahr, J.; Molenaar, M.; Bryan, F. Time variability of the Earth’s gravity field: Hydrological and oceanic effects and their possible detection using GRACE. J. Geophys. Res. 1998, 103, 30205–30229. [Google Scholar] [CrossRef]
  32. Cheng, M.K.; Ries, J.R. Monthly Estimates of C20 from 5 SLR Satellites Based on GRACE RL06 Models, GRACE Technical Note 11, The GRACE Project, Center for Space Research, University of Texas at Austin. 2019. Available online: https://podaac-tools.jpl.nasa.gov/drive/files/allData/grace/docs/TN-11_C20_SLR.txt (accessed on 22 October 2019).
  33. Landerer, F. Monthly Estimates of Degree-1 (Geocenter) Gravity Coefficients, Generated from GRACE (04-2002–06/2017) and GRACE – FO (06/2018 onward) RL06 Solutions, GRACE Technical Note 13, The GRACE Project, NASA Jet Propulsion Laboratory. 2019. Available online: https://podaac-tools.jpl.nasa.gov/drive/files/allData/grace/docs/TN-13_GEOC_CSR_RL06.txt (accessed on 17 October 2019).
  34. Swenson, S.; Chambers, D.; Wahr, J. Estimating geocenter variations from a combination of GRACE and ocean model output. J. Geophys. Res. 2008, 113, B08410. [Google Scholar] [CrossRef] [Green Version]
  35. de Linage, C.; Rivera, L.; Hinderer, J.; Boy, J.-P.; Rogister, Y.; Lambotte, S.; Biancale, R. Separation of coseismic and postseismic gravity changes for the 2004 Sumatra–Andaman earthquake from 4.6 yr of GRACE observations and modelling of the coseismic change by normal-modes summation. Geophys. J. Int. 2009, 176, 695–714. [Google Scholar] [CrossRef] [Green Version]
  36. Li, J.; Chen, J.L.; Zhang, Z. Seismologic applications of GRACE time-variable gravity measurements. Earthq. Sci. 2014, 27, 229–245. [Google Scholar] [CrossRef] [Green Version]
  37. Chen, J.L.; Wilson, C.R.; Li, J.; Zhang, Z. Reducing leakage error in GRACE-observed long-term ice mass change: A case study in West Antarctica. J. Geod. 2015, 89, 925–940. [Google Scholar] [CrossRef]
  38. Dahle, C.; Flechtner, F.; Murböck, M.; Michalak, G.; Neumayer, H.; Abrykosov, O.; Reinhold, A.; König, R. GFZ Level-2 Processing Standards Document for Level-2 Product Release 06; GRACE 327-743, (Rev. 1.0), (Scientific Technical Report STR—Data; 18/04); GFZ German Research Centre for Geosciences: Potsdam, Germany, 2018. [Google Scholar] [CrossRef]
  39. Yuan, D.-N. JPL Level-2 Processing Standards Document for Product Release 06, GRACE 327-744, Revision 6.0, Jet Propulsion Laboratory, California Institute of Technology. Available online: https://podaac-ftp.jpl.nasa.gov/allData/grace/docs/L2-JPL_ProcStds_v6.0.pdf (accessed on 1 June 2018).
  40. Paulson, A.; Zhong, S.; Wahr, J. Inference of mantle viscosity from GRACE and relative sea level data. Geophys. J. Int. 2007, 171, 497–508. [Google Scholar] [CrossRef] [Green Version]
  41. Peltier, W.R.; Argus, D.F.; Drummond, R. Space geodesy constrains ice-age terminal deglaciation: The global ICE-6G_C (VM5a) model. J. Geophys. Res. Solid Earth 2015, 120, 450–487. [Google Scholar] [CrossRef] [Green Version]
  42. Li, J.; Chen, J.L.; Wilson, C.R. Topographic effects on coseismic gravity change for the 2011 Tohoku-Oki earthquake and comparison with GRACE. J. Geophys. Res. Solid Earth 2016, 121, 5509–5537. [Google Scholar] [CrossRef]
Figure 1. (a) The spatial distribution of equivalent water height (EWH) trend rates 2003–2015 from the Gravity Recovery and Climate Experiment (GRACE) Release-06 (RL06) solutions, with the glacial isostatic adjustment (GIA) effect removed and (b) the time series at point A (2.5°N, 93.5°E) from GRACE, and the linear trend fitted by least squares fitting.
Figure 1. (a) The spatial distribution of equivalent water height (EWH) trend rates 2003–2015 from the Gravity Recovery and Climate Experiment (GRACE) Release-06 (RL06) solutions, with the glacial isostatic adjustment (GIA) effect removed and (b) the time series at point A (2.5°N, 93.5°E) from GRACE, and the linear trend fitted by least squares fitting.
Remotesensing 12 00935 g001
Figure 2. The time series fitting examples at two selected points, with the co- and post-seismic terms taken into account. (a) The time series fitting includes the effect of only one earthquake (i.e., the March 2011 Tohoku-Oki earthquake) and (b) the fitting includes the effect of two earthquakes (i.e., the December 2004 and April 2012 Sumatra earthquakes). The red curve represents the original time series observed by GRACE. The blue curves represent the fitted co- and post-seismic terms. The vertical dashed lines denote the occurrence months of the earthquakes.
Figure 2. The time series fitting examples at two selected points, with the co- and post-seismic terms taken into account. (a) The time series fitting includes the effect of only one earthquake (i.e., the March 2011 Tohoku-Oki earthquake) and (b) the fitting includes the effect of two earthquakes (i.e., the December 2004 and April 2012 Sumatra earthquakes). The red curve represents the original time series observed by GRACE. The blue curves represent the fitted co- and post-seismic terms. The vertical dashed lines denote the occurrence months of the earthquakes.
Remotesensing 12 00935 g002
Figure 3. (a) The EWH rate distribution in the Sumatra earthquake region, (b) the rate distribution after the seismic effect correction, (c) the time series at the point D (3.5°N, 94.5°E), and (d) the time series at the point E (10.5°N, 98.5°E).
Figure 3. (a) The EWH rate distribution in the Sumatra earthquake region, (b) the rate distribution after the seismic effect correction, (c) the time series at the point D (3.5°N, 94.5°E), and (d) the time series at the point E (10.5°N, 98.5°E).
Remotesensing 12 00935 g003
Figure 4. (a) The EWH rate distribution in the Tohoku-Oki earthquake region, (b) the rate distribution after the seismic effect correction, (c) the time series at the point F (36.5°N, 144.5°E), and (d) the time series at the point G (39.5°N, 138.5°E).
Figure 4. (a) The EWH rate distribution in the Tohoku-Oki earthquake region, (b) the rate distribution after the seismic effect correction, (c) the time series at the point F (36.5°N, 144.5°E), and (d) the time series at the point G (39.5°N, 138.5°E).
Remotesensing 12 00935 g004
Figure 5. (a) The equivalent water height rate distribution in the Chile earthquake region, (b) the rate distribution after the seismic effect correction; (c) the time series at the point H (−37.5°N, 283.5°E); and (d) the time series at the point I (−35.5°N, 291.5°E).
Figure 5. (a) The equivalent water height rate distribution in the Chile earthquake region, (b) the rate distribution after the seismic effect correction; (c) the time series at the point H (−37.5°N, 283.5°E); and (d) the time series at the point I (−35.5°N, 291.5°E).
Remotesensing 12 00935 g005
Figure 6. The global mean ocean mass (GMOM) change rates using forward modeling recovery after 100 iterations. The red curve represents the GMOM change rate without seismic correction. The blue curve denotes the GMOM change rate with correction for the Sumatra earthquake region. The magenta curve represents the GMOM change rate with the corrections for the Sumatra and Tohoku-Oki earthquake regions. The black curve denotes the GMOM change rate with the corrections for the three regions.
Figure 6. The global mean ocean mass (GMOM) change rates using forward modeling recovery after 100 iterations. The red curve represents the GMOM change rate without seismic correction. The blue curve denotes the GMOM change rate with correction for the Sumatra earthquake region. The magenta curve represents the GMOM change rate with the corrections for the Sumatra and Tohoku-Oki earthquake regions. The black curve denotes the GMOM change rate with the corrections for the three regions.
Remotesensing 12 00935 g006
Figure 7. (a) The global EWH change rate distribution corrected only for the Sumatra earthquake region; (b) the global rate distribution corrected for the Sumatra and Tohoku-Oki earthquake regions; and (c) the global rate distribution corrected for the Sumatra, Tohoku-Oki, and Chile earthquake regions.
Figure 7. (a) The global EWH change rate distribution corrected only for the Sumatra earthquake region; (b) the global rate distribution corrected for the Sumatra and Tohoku-Oki earthquake regions; and (c) the global rate distribution corrected for the Sumatra, Tohoku-Oki, and Chile earthquake regions.
Remotesensing 12 00935 g007
Table 1. Gravity Recovery and Climate Experiment (GRACE) global mean ocean mass (GMOM) change rates from forward modeling (FM) recovery for the period 2003–2015 without and with the seismic corrections.
Table 1. Gravity Recovery and Climate Experiment (GRACE) global mean ocean mass (GMOM) change rates from forward modeling (FM) recovery for the period 2003–2015 without and with the seismic corrections.
Without CorrectionWith Correction for Sumatra RegionWith Correction for Sumatra + Tohoku-Oki RegionsWith Correction for the Three Regions
GMOM rate (mm/yr)2.12 ± 0.11 *2.08 ± 0.12 *2.07 ± 0.12 *2.05 ± 0.12 *
* Note: the uncertainties are from 2 sigma of least-squares fitting as described subsequently in Section 4.1.
Table 2. The GMOM rates from Center for Space Research (CSR), Geoforschungszentrum (GFZ), and Jet Propulsion Laboratory (JPL) RL06 solutions for the period January 2003 to December 2015. The uncertainties are from 2 sigma of the least-squares fitting.
Table 2. The GMOM rates from Center for Space Research (CSR), Geoforschungszentrum (GFZ), and Jet Propulsion Laboratory (JPL) RL06 solutions for the period January 2003 to December 2015. The uncertainties are from 2 sigma of the least-squares fitting.
Without Correction (mm/year)With Correction (mm/year)
CSR RL062.12 ± 0.112.05 ± 0.12
GFZ RL062.05 ± 0.111.97 ± 0.11
JPL RL062.16 ± 0.112.07 ± 0.11

Share and Cite

MDPI and ACS Style

Tang, L.; Li, J.; Chen, J.; Wang, S.-Y.; Wang, R.; Hu, X. Seismic Impact of Large Earthquakes on Estimating Global Mean Ocean Mass Change from GRACE. Remote Sens. 2020, 12, 935. https://doi.org/10.3390/rs12060935

AMA Style

Tang L, Li J, Chen J, Wang S-Y, Wang R, Hu X. Seismic Impact of Large Earthquakes on Estimating Global Mean Ocean Mass Change from GRACE. Remote Sensing. 2020; 12(6):935. https://doi.org/10.3390/rs12060935

Chicago/Turabian Style

Tang, Lu, Jin Li, Jianli Chen, Song-Yun Wang, Rui Wang, and Xiaogong Hu. 2020. "Seismic Impact of Large Earthquakes on Estimating Global Mean Ocean Mass Change from GRACE" Remote Sensing 12, no. 6: 935. https://doi.org/10.3390/rs12060935

APA Style

Tang, L., Li, J., Chen, J., Wang, S. -Y., Wang, R., & Hu, X. (2020). Seismic Impact of Large Earthquakes on Estimating Global Mean Ocean Mass Change from GRACE. Remote Sensing, 12(6), 935. https://doi.org/10.3390/rs12060935

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop