![Page 1: Time-variable gravity field determination from GRACE-FO ... · Epoch-wise covariance information from kinematic PPP: mainly 2/rev; important information on observation quality Ambiguities](https://reader033.vdocument.in/reader033/viewer/2022050602/5fa9284b6feee711ab6c66eb/html5/thumbnails/1.jpg)
Time-variable gravity fielddetermination from GRACE-FO data
at AIUB
Martin Lasser1, Ulrich Meyer1, Daniel Arnold1, Adrian Jaggi1
1Astronomical Institute, University of Bern, Switzerland
GRACE & GRACE-FO Science Team MeetingOct. 8-10, 2019
Pasadena, California
Astronomical Institute, University of Bern AIUB
![Page 2: Time-variable gravity field determination from GRACE-FO ... · Epoch-wise covariance information from kinematic PPP: mainly 2/rev; important information on observation quality Ambiguities](https://reader033.vdocument.in/reader033/viewer/2022050602/5fa9284b6feee711ab6c66eb/html5/thumbnails/2.jpg)
M.Lasser,U.Meyer,D.Arnold,A.Ja
ggi:
StatusofGFO
data
processingatAIU
BGRACE&
GRACE-FO
Scien
ceTeam
Meeting,Oct.8-10,2019,Pasaden
a,California
Content
Kinematic orbits of GRACE-FO satellites
Stochastic behaviour of kinematic positions
Empirical modelling
Inter-satellite link
Impact on gravity field determination
Slide 2 of 16 Astronomical Institute, University of Bern AIUB
![Page 3: Time-variable gravity field determination from GRACE-FO ... · Epoch-wise covariance information from kinematic PPP: mainly 2/rev; important information on observation quality Ambiguities](https://reader033.vdocument.in/reader033/viewer/2022050602/5fa9284b6feee711ab6c66eb/html5/thumbnails/3.jpg)
M.Lasser,U.Meyer,D.Arnold,A.Ja
ggi:
StatusofGFO
data
processingatAIU
BGRACE&
GRACE-FO
Scien
ceTeam
Meeting,Oct.8-10,2019,Pasaden
a,California
Observations
Credit:DLR
GPS Inter-satellite link
absolute relative3-dim 1-dimlower precision very high precision
Slide 3 of 16 Astronomical Institute, University of Bern AIUB
![Page 4: Time-variable gravity field determination from GRACE-FO ... · Epoch-wise covariance information from kinematic PPP: mainly 2/rev; important information on observation quality Ambiguities](https://reader033.vdocument.in/reader033/viewer/2022050602/5fa9284b6feee711ab6c66eb/html5/thumbnails/4.jpg)
M.Lasser,U.Meyer,D.Arnold,A.Ja
ggi:
StatusofGFO
data
processingatAIU
BGRACE&
GRACE-FO
Scien
ceTeam
Meeting,Oct.8-10,2019,Pasaden
a,California
Kinematic orbits of GRACE-FO satellites
LEO precise orbit determination from GPS data
• 10 s sampling (full sampling used in processing)
• undifferenced PPP solution
• fixing zero-difference ambiguities to integers
→ CODE (Center of Orbit Determination in Europe) phase biases
→ CODE ambiguity-fixed clocks
Slide 4 of 16 Astronomical Institute, University of Bern AIUB
![Page 5: Time-variable gravity field determination from GRACE-FO ... · Epoch-wise covariance information from kinematic PPP: mainly 2/rev; important information on observation quality Ambiguities](https://reader033.vdocument.in/reader033/viewer/2022050602/5fa9284b6feee711ab6c66eb/html5/thumbnails/5.jpg)
M.Lasser,U.Meyer,D.Arnold,A.Ja
ggi:
StatusofGFO
data
processingatAIU
BGRACE&
GRACE-FO
Scien
ceTeam
Meeting,Oct.8-10,2019,Pasaden
a,California
Integer-fixing of carrier phase ambiguities
SLR residuals for 2018-11 (in comparison with float)SV type mean [mm] std. dev [mm]GRACE-C float 3.1 18.6GRACE-D float -4.5 16.8GRACE-C fixed 4.3 10.6GRACE-D fixed 1.0 12.0
K-band range validation of kinematic positions
→ Individual arcs feature less scattering
Slide 5 of 16 Astronomical Institute, University of Bern AIUB
![Page 6: Time-variable gravity field determination from GRACE-FO ... · Epoch-wise covariance information from kinematic PPP: mainly 2/rev; important information on observation quality Ambiguities](https://reader033.vdocument.in/reader033/viewer/2022050602/5fa9284b6feee711ab6c66eb/html5/thumbnails/6.jpg)
M.Lasser,U.Meyer,D.Arnold,A.Ja
ggi:
StatusofGFO
data
processingatAIU
BGRACE&
GRACE-FO
Scien
ceTeam
Meeting,Oct.8-10,2019,Pasaden
a,California
Integer-fixing of carrier phase ambiguities
SLR residuals for 2018-11 (in comparison with float)SV type mean [mm] std. dev [mm]GRACE-C float 3.1 18.6GRACE-D float -4.5 16.8GRACE-C fixed 4.3 10.6GRACE-D fixed 1.0 12.0
K-band range validation of kinematic positions
→ Individual arcs feature less scattering
Slide 5 of 16 Astronomical Institute, University of Bern AIUB
![Page 7: Time-variable gravity field determination from GRACE-FO ... · Epoch-wise covariance information from kinematic PPP: mainly 2/rev; important information on observation quality Ambiguities](https://reader033.vdocument.in/reader033/viewer/2022050602/5fa9284b6feee711ab6c66eb/html5/thumbnails/7.jpg)
M.Lasser,U.Meyer,D.Arnold,A.Ja
ggi:
StatusofGFO
data
processingatAIU
BGRACE&
GRACE-FO
Scien
ceTeam
Meeting,Oct.8-10,2019,Pasaden
a,California
Stochastic behaviour of kinematic positions
• SST-hl tracking geometry varies between poles and the equator
• Epoch-wise covariance information from kinematic PPP:mainly 2/rev; important information on observation quality
• Ambiguities feature the correlation between epochs
→ Full covariance matrix from PPP (at least in principle)
→ Integer-fixed: almost no correlations
Slide 6 of 16 Astronomical Institute, University of Bern AIUB
![Page 8: Time-variable gravity field determination from GRACE-FO ... · Epoch-wise covariance information from kinematic PPP: mainly 2/rev; important information on observation quality Ambiguities](https://reader033.vdocument.in/reader033/viewer/2022050602/5fa9284b6feee711ab6c66eb/html5/thumbnails/8.jpg)
M.Lasser,U.Meyer,D.Arnold,A.Ja
ggi:
StatusofGFO
data
processingatAIU
BGRACE&
GRACE-FO
Scien
ceTeam
Meeting,Oct.8-10,2019,Pasaden
a,California
Stochastic behaviour of kinematic positions
• SST-hl tracking geometry varies between poles and the equator
• Epoch-wise covariance information from kinematic PPP:mainly 2/rev; important information on observation quality
• Ambiguities feature the correlation between epochs
→ Full covariance matrix from PPP (at least in principle)
→ Integer-fixed: almost no correlations
Slide 6 of 16 Astronomical Institute, University of Bern AIUB
![Page 9: Time-variable gravity field determination from GRACE-FO ... · Epoch-wise covariance information from kinematic PPP: mainly 2/rev; important information on observation quality Ambiguities](https://reader033.vdocument.in/reader033/viewer/2022050602/5fa9284b6feee711ab6c66eb/html5/thumbnails/9.jpg)
M.Lasser,U.Meyer,D.Arnold,A.Ja
ggi:
StatusofGFO
data
processingatAIU
BGRACE&
GRACE-FO
Scien
ceTeam
Meeting,Oct.8-10,2019,Pasaden
a,California
Stochastic behaviour of kinematic positions
• SST-hl tracking geometry varies between poles and the equator
• Epoch-wise covariance information from kinematic PPP:mainly 2/rev; important information on observation quality
• Ambiguities feature the correlation between epochs
→ Full covariance matrix from PPP (at least in principle)
→ Integer-fixed: almost no correlations
Slide 6 of 16 Astronomical Institute, University of Bern AIUB
![Page 10: Time-variable gravity field determination from GRACE-FO ... · Epoch-wise covariance information from kinematic PPP: mainly 2/rev; important information on observation quality Ambiguities](https://reader033.vdocument.in/reader033/viewer/2022050602/5fa9284b6feee711ab6c66eb/html5/thumbnails/10.jpg)
M.Lasser,U.Meyer,D.Arnold,A.Ja
ggi:
StatusofGFO
data
processingatAIU
BGRACE&
GRACE-FO
Scien
ceTeam
Meeting,Oct.8-10,2019,Pasaden
a,California
Empirical modelling
Non linear least squares: x = x0 + (ATPA)−1ATPδl
→ The observations enter via δl = o− c
• Residuals feature observation noise (o)
• And mismodelings from the computed component (c)
→ c contains the full a priori force model(gravity field, tides, AOD, accelerometer, ...)
Slide 7 of 16 Astronomical Institute, University of Bern AIUB
![Page 11: Time-variable gravity field determination from GRACE-FO ... · Epoch-wise covariance information from kinematic PPP: mainly 2/rev; important information on observation quality Ambiguities](https://reader033.vdocument.in/reader033/viewer/2022050602/5fa9284b6feee711ab6c66eb/html5/thumbnails/11.jpg)
M.Lasser,U.Meyer,D.Arnold,A.Ja
ggi:
StatusofGFO
data
processingatAIU
BGRACE&
GRACE-FO
Scien
ceTeam
Meeting,Oct.8-10,2019,Pasaden
a,California
Empirical modelling
Non linear least squares: x = x0 + (ATPA)−1ATPδl
→ The observations enter via δl = o− c
• Residuals feature observation noise (o)
• And mismodelings from the computed component (c)
→ c contains the full a priori force model(gravity field, tides, AOD, accelerometer, ...)
Slide 7 of 16 Astronomical Institute, University of Bern AIUB
![Page 12: Time-variable gravity field determination from GRACE-FO ... · Epoch-wise covariance information from kinematic PPP: mainly 2/rev; important information on observation quality Ambiguities](https://reader033.vdocument.in/reader033/viewer/2022050602/5fa9284b6feee711ab6c66eb/html5/thumbnails/12.jpg)
M.Lasser,U.Meyer,D.Arnold,A.Ja
ggi:
StatusofGFO
data
processingatAIU
BGRACE&
GRACE-FO
Scien
ceTeam
Meeting,Oct.8-10,2019,Pasaden
a,California
Empirical modelling
Non linear least squares: x = x0 + (ATPA)−1ATPδl
→ The observations enter via δl = o− c
• Residuals feature observation noise (o)
• And mismodelings from the computed component (c)
→ c contains the full a priori force model(gravity field, tides, AOD, accelerometer, ...)
Slide 7 of 16 Astronomical Institute, University of Bern AIUB
![Page 13: Time-variable gravity field determination from GRACE-FO ... · Epoch-wise covariance information from kinematic PPP: mainly 2/rev; important information on observation quality Ambiguities](https://reader033.vdocument.in/reader033/viewer/2022050602/5fa9284b6feee711ab6c66eb/html5/thumbnails/13.jpg)
M.Lasser,U.Meyer,D.Arnold,A.Ja
ggi:
StatusofGFO
data
processingatAIU
BGRACE&
GRACE-FO
Scien
ceTeam
Meeting,Oct.8-10,2019,Pasaden
a,California
Empirical modelling
Non linear least squares: x = x0 + (ATPA)−1ATPδl
→ The observations enter via δl = o− c
Slide 8 of 16 Astronomical Institute, University of Bern AIUB
![Page 14: Time-variable gravity field determination from GRACE-FO ... · Epoch-wise covariance information from kinematic PPP: mainly 2/rev; important information on observation quality Ambiguities](https://reader033.vdocument.in/reader033/viewer/2022050602/5fa9284b6feee711ab6c66eb/html5/thumbnails/14.jpg)
M.Lasser,U.Meyer,D.Arnold,A.Ja
ggi:
StatusofGFO
data
processingatAIU
BGRACE&
GRACE-FO
Scien
ceTeam
Meeting,Oct.8-10,2019,Pasaden
a,California
Empirical covariance estimation
• Set up (constrained) piecewise constant accelerations(strategy used at AIUB)
In case of a stationary process
• Estimate a mean covariance function from the residuals e.
cov(∆tk) =1
n
n∑i=0
e(ti)e(ti + ∆tk)
• Variance/covariance matrix has a symmetric Toeplitz-structureand the individual elements only depend on the distance in time.
→ Mean covariance function for each month
→ Length of correlation 100 min
Slide 9 of 16 Astronomical Institute, University of Bern AIUB
![Page 15: Time-variable gravity field determination from GRACE-FO ... · Epoch-wise covariance information from kinematic PPP: mainly 2/rev; important information on observation quality Ambiguities](https://reader033.vdocument.in/reader033/viewer/2022050602/5fa9284b6feee711ab6c66eb/html5/thumbnails/15.jpg)
M.Lasser,U.Meyer,D.Arnold,A.Ja
ggi:
StatusofGFO
data
processingatAIU
BGRACE&
GRACE-FO
Scien
ceTeam
Meeting,Oct.8-10,2019,Pasaden
a,California
Inter-satellite link
The inter-satellite link is much more precise
→ dominates the resulting gravity field signal
→ 5 s sampling for K-band, 2 s sampling for LRI
• Assign one factor to weight K-band
→ extended to weight based on its arc-wise residual RMS
• Set up constrained PCA to model stochastic behaviour(mainly a 2/rev signal)
Slide 10 of 16 Astronomical Institute, University of Bern AIUB
![Page 16: Time-variable gravity field determination from GRACE-FO ... · Epoch-wise covariance information from kinematic PPP: mainly 2/rev; important information on observation quality Ambiguities](https://reader033.vdocument.in/reader033/viewer/2022050602/5fa9284b6feee711ab6c66eb/html5/thumbnails/16.jpg)
M.Lasser,U.Meyer,D.Arnold,A.Ja
ggi:
StatusofGFO
data
processingatAIU
BGRACE&
GRACE-FO
Scien
ceTeam
Meeting,Oct.8-10,2019,Pasaden
a,California
Inter-satellite link
The inter-satellite link is much more precise
→ dominates the resulting gravity field signal
→ 5 s sampling for K-band, 2 s sampling for LRI
• Assign one factor to weight K-band
→ extended to weight based on its arc-wise residual RMS
• Set up constrained PCA to model stochastic behaviour(mainly a 2/rev signal)
Slide 10 of 16 Astronomical Institute, University of Bern AIUB
![Page 17: Time-variable gravity field determination from GRACE-FO ... · Epoch-wise covariance information from kinematic PPP: mainly 2/rev; important information on observation quality Ambiguities](https://reader033.vdocument.in/reader033/viewer/2022050602/5fa9284b6feee711ab6c66eb/html5/thumbnails/17.jpg)
M.Lasser,U.Meyer,D.Arnold,A.Ja
ggi:
StatusofGFO
data
processingatAIU
BGRACE&
GRACE-FO
Scien
ceTeam
Meeting,Oct.8-10,2019,Pasaden
a,California
Inter-satellite link
• Largedistance intime betweenreferenceepoch ofbackgroundfield andobservations
Slide 11 of 16 Astronomical Institute, University of Bern AIUB
![Page 18: Time-variable gravity field determination from GRACE-FO ... · Epoch-wise covariance information from kinematic PPP: mainly 2/rev; important information on observation quality Ambiguities](https://reader033.vdocument.in/reader033/viewer/2022050602/5fa9284b6feee711ab6c66eb/html5/thumbnails/18.jpg)
M.Lasser,U.Meyer,D.Arnold,A.Ja
ggi:
StatusofGFO
data
processingatAIU
BGRACE&
GRACE-FO
Scien
ceTeam
Meeting,Oct.8-10,2019,Pasaden
a,California
Inter-satellite link
• Largedistance intime betweenreferenceepoch ofbackgroundfield andobservations
Slide 11 of 16 Astronomical Institute, University of Bern AIUB
![Page 19: Time-variable gravity field determination from GRACE-FO ... · Epoch-wise covariance information from kinematic PPP: mainly 2/rev; important information on observation quality Ambiguities](https://reader033.vdocument.in/reader033/viewer/2022050602/5fa9284b6feee711ab6c66eb/html5/thumbnails/19.jpg)
M.Lasser,U.Meyer,D.Arnold,A.Ja
ggi:
StatusofGFO
data
processingatAIU
BGRACE&
GRACE-FO
Scien
ceTeam
Meeting,Oct.8-10,2019,Pasaden
a,California
Impact on gravity field determination
• Formal errors aretoo optimistic
Slide 12 of 16 Astronomical Institute, University of Bern AIUB
![Page 20: Time-variable gravity field determination from GRACE-FO ... · Epoch-wise covariance information from kinematic PPP: mainly 2/rev; important information on observation quality Ambiguities](https://reader033.vdocument.in/reader033/viewer/2022050602/5fa9284b6feee711ab6c66eb/html5/thumbnails/20.jpg)
M.Lasser,U.Meyer,D.Arnold,A.Ja
ggi:
StatusofGFO
data
processingatAIU
BGRACE&
GRACE-FO
Scien
ceTeam
Meeting,Oct.8-10,2019,Pasaden
a,California
Impact on gravity field determination
• Slight improvementdue to bettercharacterisation ofrange-rate data
• Noise inbackground forcefield modelledthrough PCA only
Slide 13 of 16 Astronomical Institute, University of Bern AIUB
![Page 21: Time-variable gravity field determination from GRACE-FO ... · Epoch-wise covariance information from kinematic PPP: mainly 2/rev; important information on observation quality Ambiguities](https://reader033.vdocument.in/reader033/viewer/2022050602/5fa9284b6feee711ab6c66eb/html5/thumbnails/21.jpg)
M.Lasser,U.Meyer,D.Arnold,A.Ja
ggi:
StatusofGFO
data
processingatAIU
BGRACE&
GRACE-FO
Scien
ceTeam
Meeting,Oct.8-10,2019,Pasaden
a,California
Impact on gravity field determination
• Noise inbackground forcefield mainlycharacterised due toGPS contribution
Slide 14 of 16 Astronomical Institute, University of Bern AIUB
![Page 22: Time-variable gravity field determination from GRACE-FO ... · Epoch-wise covariance information from kinematic PPP: mainly 2/rev; important information on observation quality Ambiguities](https://reader033.vdocument.in/reader033/viewer/2022050602/5fa9284b6feee711ab6c66eb/html5/thumbnails/22.jpg)
M.Lasser,U.Meyer,D.Arnold,A.Ja
ggi:
StatusofGFO
data
processingatAIU
BGRACE&
GRACE-FO
Scien
ceTeam
Meeting,Oct.8-10,2019,Pasaden
a,California
Impact on gravity field determination
• Formal errors inhigher degreesrealistic
• Low degrees requirefurther inspection
Slide 15 of 16 Astronomical Institute, University of Bern AIUB
![Page 23: Time-variable gravity field determination from GRACE-FO ... · Epoch-wise covariance information from kinematic PPP: mainly 2/rev; important information on observation quality Ambiguities](https://reader033.vdocument.in/reader033/viewer/2022050602/5fa9284b6feee711ab6c66eb/html5/thumbnails/23.jpg)
M.Lasser,U.Meyer,D.Arnold,A.Ja
ggi:
StatusofGFO
data
processingatAIU
BGRACE&
GRACE-FO
Scien
ceTeam
Meeting,Oct.8-10,2019,Pasaden
a,California
Summary
• Integer-fixing of carrier phase ambiguities closer to K-band
• Empirical noise modelling via GPS improves formal errors
→ Extend the procedure to K-band/Laser link
• Processing strategy tested on GRACE-FO
→ (At some point) reprocess GRACE time series
Thank you for your attention.
Slide 16 of 16 Astronomical Institute, University of Bern AIUB