Contribution to IAG JSG 03.

IAG JSG03 contribution

By C.C.Tscherning, University of Copenhagen.
IAG JSG 0.3 objectives are found here.

Use of LSC and RPM for anomalous potential prediction.

First are described initial experiments reported at EGU2014.
Secondly, experiments inspired by the comparison results a number of new tests with noise added to the data and using different reference fields are in progress.
The various simulated test data sets have been used for the computation of values of the anomalous potential and a comparison of values generated using EGM2008 to the maximal degree. An European and a Pacific test area are used. Results reported below have been obtained using the programs empcov, covfit and geocol19.
Covariance functions:
Europe-EGM to 240 . Pacific -240
Europe-EGM to 36 . Pacific -36

From ground gravity and EGM2008 to deg. 240.

All "grob" and "fine" Tr data used in the test-areas with EGM2008 coefficients to maximal degree 240 as observations.

Europe

Anomalous potential differences, Simulated EGM2008 values - LSC predicted.
COMPARISON OF PREDICTIONS AND OBSERVATIONS
_DATA TYPE = 11, m**2/s**2.
NUMBER: 8712
0 OBSERVATIONS PREDICTIONS DIFFERENCE ERROR ESTIMATES
MEAND -0.01 0.02 -0.03 0.80
ST.DEVI. 4.08 4.19 0.82 0.38
MAX 19.93 22.15 10.43 2.58
MIN -16.24 -15.62 -3.67 0.35
File with differences and error-estimate. Format: id#, lat., long., h, sim-model, difference, error-estimate. Units m**2/s**2.

Pacific.

Anomalous potential differences, Simulated EGM2008 values - LSC predicted.
COMPARISON OF PREDICTIONS AND OBSERVATIONS
0DATA TYPE = 11,
NUMBER: 8640
0 OBSERVATIONS PREDICTIONS DIFFERENCE ERROR ESTIMATES
MEAND -0.127288 -0.071543 -0.055745 0.788378
ST.DEVI. 4.817569 4.786428 0.243873 0.292099
MAX 22.203400 22.426302 1.766949 1.529044
MIN -28.710900 -27.744831 -3.059059 0.108776

File with differences and error-estimate.

From airborne gravity and EGM2008 to deg. 240.

Europe.

COMPARISON OF PREDICTIONS AND OBSERVATIONS
0DATA TYPE = 11
NUMBER: 8712
0 OBSERVATIONS PREDICTIONS DIFFERENCE ERROR ESTIMATES
MEAND -0.01 0.10 -0.11 2.97
ST.DEVI. 4.08 3.55 2.73 1.73
MAX 19.93 13.74 16.87 5.02
MIN -16.24 -14.33 -15.45 0.77
File with differences and error-estimate.

Pacific.

COMPARISON OF PREDICTIONS AND OBSERVATIONS
DATA TYPE = 11
NUMBER: 8640
0 OBSERVATIONS PREDICTIONS DIFFERENCE ERROR ESTIMATES
MEAND -0.13 -0.05 -0.08 1.65
ST.DEVI. 4.82 4.60 1.92 1.31
MAX 22.20 23.52 13.43 4.80
MIN -28.71 -27.81 -15.65 0.38
File with differences and error-estimate.

From GOCE Tzz and EGM2008 to deg. 36.

Note: here is predicted both T and Tr !

Europe.

COMPARISON OF PREDICTIONS AND OBSERVATIONS
DATA TYPE = 11 (T, m**2/s**2.
NUMBER: 8712
0 OBSERVATIONS PREDICTIONS DIFFERENCE ERROR ESTIMATES
MEAND___ -0.114069 0.145354 -0.259423 6.795504
ST.DEVI._ 32.548924 31.441345 6.032833 0.105642
MAX_____ 69.859200 63.356536 31.348700 6.897625
MIN_____ -92.096800 -75.522027 -27.790924 6.333314
File with differences and error-estimate.
COMPARISON OF PREDICTIONS AND OBSERVATIONS
0DATA TYPE = 12, Tr, mgal.
0 OBSERVATIONS PREDICTIONS DIFFERENCE ERROR ESTIMATES
MEAND 3.56 2.29 1.27 31.82
ST.DEVI. 49.13 37.10 29.44 0.86
MAX 175.51 110.54 126.27 32.42
MIN -196.53 -78.82 -170.48 27.97

From GOCE Tzz and EGM2008 to deg. 60.

COMPARISON OF PREDICTIONS AND OBSERVATIONS
DATA TYPE = 11 (T, m**2/s**2.
NUMBER: 8712
0 OBSERVATIONS PREDICTIONS DIFFERENCE ERROR ESTIMATES
MEAND 378.84 378.89 -0.05 7.09
ST.DEVI. 96.65 96.57 5.14 0.33
MAX 531.56 526.40 33.01 8.91
MIN 1.80 6.09 -32.87 6.19
COMPARISON OF PREDICTIONS AND OBSERVATIONS
0DATA TYPE = 12, Tr, mgal.
0 OBSERVATIONS PREDICTIONS DIFFERENCE ERROR ESTIMATES
MEAND -21.58 -21.81 0.23 30.97
ST.DEVI. 46.35 38.09 25.53 1.04
MAX 227.61 165.87 144.21 35.00
MIN -209.41 -123.33 -161.60 25.60

From GOCE Tzz and EGM2008 to deg.120.

COMPARISON OF PREDICTIONS AND OBSERVATIONS
DATA TYPE = 11 (T, m**2/s**2.
NUMBER: 8712
0 OBSERVATIONS PREDICTIONS DIFFERENCE ERROR ESTIMATES
MEAND 378.835383 378.911519 -0.076136 6.946962
ST.DEVI. 96.646546 96.519195 5.111822 0.202579
MAX 531.564000 526.292504 29.413262 7.837187
MIN 1.802000 5.491850 -34.383029 6.124164
0DATA TYPE = 12, Tr, mgal.
NUMBER: 2501
0 OBSERVATIONS PREDICTIONS DIFFERENCE ERROR ESTIMATES
MEAND -21.58 -21.84 0.26 31.21
ST.DEVI. 46.35 38.06 25.45 0.90
MAX 227.61 166.39 141.01 33.58
MIN -209.41 -123.39 -149.43 25.71

Pacific.

Note: here is predicted both T and Tr !
COMPARISON OF PREDICTIONS AND OBSERVATIONS
DATA TYPE = 11 (T, m**2/s**2.
NUMBER: 8640
0 OBSERVATIONS PREDICTIONS DIFFERENCE ERROR ESTIMATES
MEAND 4.62 4.80 -0.18 8.41
ST.DEVI. 49.07 48.52 6.66 0.44
MAX 115.52 105.67 26.89 8.70<
MIN -106.43 -96.42 -33.56 7.02
File with differences and error-estimate.
COMPARISON OF PREDICTIONS AND OBSERVATIONS
DATA TYPE = 12 (Tr, mgal).
NUMBER: 8833
0 OBSERVATIONS PREDICTIONS DIFFERENCE ERROR ESTIMATES
MEAND -2.22 -3.60 1.39 36.61
ST.DEVI. 65.25 52.66 38.24 2.71
MAX 191.28 117.51 232.00 38.18
MIN -286.54 -124.70 -198.30 28.12

Prediction of dT from Tr and EGM2008 to 36, Europe.

NUMBER: 514
0 OBSERVATIONS PREDICTIONS DIFFERENCE ERROR ESTIMATES
MEAND 2.800407 2.816829 -0.016422 0.138451
ST.DEVI. 10.916129 10.920700 0.172240 0.018781
MAX 23.037000 22.873253 0.278124 0.178281
MIN -18.875500 -18.472153 -0.488373 0.089580

Prediction from dT (GRACE) and EGM2008 to 36, Europe.

Note: here is predicted both T and Tr !
DATA TYPE = 11, T (m**2/s**2).
NUMBER: 8712
0 OBSERVATIONS PREDICTIONS DIFFERENCE ERROR ESTIMATES
MEAND -0.11 1.64 -1.76 28.75
ST.DEVI. 32.55 18.40 19.15 0.58
MAX 69.86 36.95 46.97 29.90
MIN -92.10 -39.03 -65.46 27.27
DATA TYPE = 12, Tr (mgal).
NUMBER: 8833
0 OBSERVATIONS PREDICTIONS DIFFERENCE ERROR ESTIMATES
MEAND 3.56 -0.26 3.82 45.65
ST.DEVI. 49.13 16.18 42.24 0.69
MAX 175.51 32.59 148.70 46.57
MIN -196.53 -31.75 -185.46 42.36

Prediction from dT (GRACE) and EGM2008 to 36, Pacific.

Note: here is predicted both T and Tr !
DATA TYPE = 11, T (m**2/s**2).
NUMBER: 8640
0 OBSERVATIONS PREDICTIONS DIFFERENCE ERROR ESTIMATES
MEAND 4.62 4.43 0.19 24.64
ST.DEVI. 49.07 38.92 18.57 0.84
MAX 115.52 81.98 51.56 26.46
MIN -106.43 -73.85 -53.38 22.09

DATA TYPE = 12, Tr, mgal.
NUMBER: 8833
0 OBSERVATIONS PREDICTIONS DIFFERENCE ERROR ESTIMATES
MEAND -5.87 -7.31 1.44 52.00
ST.DEVI. 70.38 44.45 48.09 2.55
MAX 169.32 60.57 245.41 54.17
MIN -310.24 -129.16 -237.78 43.63

Results from reduced point mass modelling compared to LSC.

Report .

Experiments with the removal of various degrees of EGM2008.

Statistics and summary results.

Experiments with new covariance functions and noisy data.

Europe, Tr fein+grob as obs., N=240. T_val predicted.

0 OBSERVATIONS PREDICTIONS DIFFERENCE ERROR ESTIMATES
MEAND -0.01 0.03 -0.04 0.82
ST.DEVI. 4.08 4.19 0.83 0.38
MAX 19.93 22.22 10.58 2.61
MIN -16.24 -15.58 -3.78 0.36

Europe, Tr fein+grob+Tr_flugA as obs., N=240. T_val predicted.

0 OBSERVATIONS PREDICTIONS DIFFERENCE ERROR ESTIMATES
MEAND -0.01 0.06 -0.07 0.78
ST.DEVI. 4.08 4.21 0.85 0.38
MAX 19.93 22.20 10.60 2.61
MIN -16.24 -15.59 -3.75 0.33

Europe, Tr fein+grob+Tr_flugA+Tr_flugC as obs., N=240. T_val predicted.

0 OBSERVATIONS PREDICTIONS DIFFERENCE ERROR ESTIMATES
MEAND -0.01 0.05 -0.06 0.68
ST.DEVI. 4.08 4.21 0.85 0.36
MAX 19.93 22.18 10.63 2.61
MIN -16.24 -15.61 -3.73 0.28

Pacific, Tr fein+grob as obs., N=240. T_val predicted.

0 OBSERVATIONS PREDICTIONS DIFFERENCE ERROR ESTIMATES
MEAND -0.127288 -0.073991 -0.053297 0.584574
ST.DEVI. 4.817569 4.795121 0.222599 0.211449
MAX 22.203400 22.395545 1.879677 1.277941
MIN -28.710900 -27.920164 -2.925266 0.091494
:

Pacific, Tr fein+grob+Tr_flugP as obs., N=240. T_val predicted.

0 OBSERVATIONS PREDICTIONS DIFFERENCE ERROR ESTIMATES
MEAND -0.127288 -0.074644 -0.052644 0.467215
ST.DEVI. 4.817569 4.806413 0.169535 0.200361
MAX 22.203400 22.319976 0.912890 1.275459
MIN -28.710900 -27.960953 -1.992999 0.073658

Europe, Tr fein+grob as obs., N=120. T_val predicted.

0 OBSERVATIONS PREDICTIONS DIFFERENCE ERROR ESTIMATES
MEAND -0.67 -0.66 -0.01 0.59
ST.DEVI. 10.90 11.12 0.83 0.36
MAX 36.04 37.40 10.46 2.23
MIN -41.60 -41.31 -5.37 0.22

Europe, Tr fein+grob+Tr_flugA as obs., N=120. T_val predicted.

0 OBSERVATIONS PREDICTIONS DIFFERENCE ERROR ESTIMATES
MEAND -0.67 -0.65 -0.02 0.54
ST.DEVI. 10.90 11.09 0.75 0.34
MAX 36.04 37.37 8.33 2.22
MIN -41.60 -41.30 -5.06 0.20

Europe, Tr fein+grob+Tr_flugA+Tr_flugC as obs., N=120. T_val predicted.

0 OBSERVATIONS PREDICTIONS DIFFERENCE ERROR ESTIMATES
MEAND -0.67 -0.62 -0.05 0.45
ST.DEVI. 10.90 11.06 0.73 0.32
MAX 36.04 37.38 8.36 2.22
MIN -41.60 -41.31 -5.02 0.17

Pacific, Tr fein+grob as obs., N=120. T_val predicted.

0 OBSERVATIONS PREDICTIONS DIFFERENCE ERROR ESTIMATES
MEAND -0.20 -0.09 -0.10 0.80
ST.DEVI. 10.16 10.10 0.55 0.35
MAX 40.77 41.06 4.56 2.66
MIN -40.01 -40.22 -5.40 0.17

Pacific, Tr fein+grob+Tr_flugP as obs., N=120. T_val predicted.

0 OBSERVATIONS PREDICTIONS DIFFERENCE ERROR ESTIMATES
MEAND -0.20 -0.08 -0.11 0.60
ST.DEVI. 10.16 10.13 0.37 0.33
MAX 40.77 41.13 2.58 2.65
MIN -40.01 -39.92 -4.98 0.12

Europe, Tr fein+grob as obs., N=60. T_val predicted.

0 OBSERVATIONS PREDICTIONS DIFFERENCE ERROR ESTIMATES
MEAND -3.03 -2.98 -0.05 0.80
ST.DEVI. 22.11 22.55 1.54 0.42
MAX 48.64 49.81 16.64 2.62
MIN -72.47 -69.05 -5.87 0.34

Europe, Tr fein+grob+Tr_flugA as obs., N=60. T_val predicted.

0 OBSERVATIONS PREDICTIONS DIFFERENCE ERROR ESTIMATES
MEAND -3.03 -2.81 -0.22 0.73
ST.DEVI. 22.11 22.38 1.20 0.40
MAX 48.64 49.70 10.71 2.59
MIN -72.47 -68.61 -5.99 0.31

Europe, Tr fein+grob+Tr_flugA+Tr_flugC as obs., N=60. T_val predicted.

0 OBSERVATIONS PREDICTIONS DIFFERENCE ERROR ESTIMATES
MEAND -3.03 -2.82 -0.21 0.60
ST.DEVI. 22.11 22.44 1.20 0.37
MAX 48.64 49.65 10.74 2.58
MIN -72.47 -69.75 -5.95 0.26

Pacific, Tr fein+grob as obs., N=60. T_val predicted.

0 OBSERVATIONS PREDICTIONS DIFFERENCE ERROR ESTIMATES
MEAND 0.71 1.42 -0.70 0.60
ST.DEVI. 32.84 33.55 2.35 0.47
MAX 79.77 81.80 7.55 2.64
MIN -75.09 -75.00 -24.16 0.16

Pacific, Tr fein+grob+Tr_flugP as obs., N=60. T_val predicted.

0 OBSERVATIONS PREDICTIONS DIFFERENCE ERROR ESTIMATES
MEAND 0.71 1.34 -0.62 0.43
ST.DEVI. 32.84 33.24 1.46 0.38
MAX 79.77 80.50 1.98 2.63
MIN -75.09 -74.92 -13.16 0.14

Europe, Tr fein+grob as obs., N=36. T_val predicted.

0 OBSERVATIONS PREDICTIONS DIFFERENCE ERROR ESTIMATES
MEAND -0.11 -0.37 0.26 1.06
ST.DEVI. 32.55 32.76 1.91 0.60
MAX 69.86 69.46 20.02 3.58
MIN -92.10 -88.87 -5.25 0.35

Europe, Tr fein+grob+Tr_flugA as obs., N=36. T_val predicted.

0 OBSERVATIONS PREDICTIONS DIFFERENCE ERROR ESTIMATES
MEAND -0.11 -0.30 0.18 0.99
ST.DEVI. 32.55 32.63 1.71 0.58
MAX 69.86 69.43 17.03 3.57
MIN -92.10 -88.39 -5.59 0.32

Europe, Tr fein+grob+Tr_flugA+Tr_flugC as obs., N=36. T_val predicted.

0 OBSERVATIONS PREDICTIONS DIFFERENCE ERROR ESTIMATES
MEAND -0.67 -0.62 -0.05 0.45
ST.DEVI. 10.90 11.06 0.73 0.32
MAX 36.04 37.38 8.36 2.22
MIN -41.60 -41.31 -5.02 0.17

Pacific, Tr fein+grob as obs., N=36. T_val predicted.

0 OBSERVATIONS PREDICTIONS DIFFERENCE ERROR ESTIMATES
MEAND 4.62 6.21 -1.60 1.04
ST.DEVI. 49.07 50.30 3.19 0.66
MAX 115.52 121.92 3.41 3.60
MIN -106.43 -106.00 -28.70 0.27

Pacific, Tr fein+grob+Tr_flugP as obs., N=36. T_val predicted.

0 OBSERVATIONS PREDICTIONS DIFFERENCE ERROR ESTIMATES
MEAND 4.62 6.11 -1.50 0.80
ST.DEVI. 49.07 50.04 2.55 0.58
MAX 115.52 118.94 3.11 3.59
MIN -106.43 -105.83 -22.53 0.22
Last update 2014-08-12.