the estimation of 550 km x 550 km mean gravity …€¦ · the estimation of 550 km x 550 km mean...

26
THE ESTIMATION OF 550 km X 550 km MEAN GRAVITY ANOMALIES M. R. WILLIAMSON and E. M. GAPOSCHKIN (NASA-CR-142297) THE ESTIMATION OF 550 km X N75-1909- 550 km MEAN GRAVITY ANOMALIES (Smithsonian Astrophysical Observatory) 26 p HC $3.75 CSCL 08N Unclas G3/89 13601 Smithsonian Astrophysical Observatory SPECIAL REPORT 363 https://ntrs.nasa.gov/search.jsp?R=19750011023 2020-04-30T17:53:08+00:00Z

Upload: others

Post on 30-Apr-2020

18 views

Category:

Documents


0 download

TRANSCRIPT

THE ESTIMATION OF550 km X 550 km MEANGRAVITY ANOMALIES

M. R. WILLIAMSON and E. M. GAPOSCHKIN

(NASA-CR-142297) THE ESTIMATION OF 550 km X N75-1909-550 km MEAN GRAVITY ANOMALIES (SmithsonianAstrophysical Observatory) 26 p HC $3.75

CSCL 08N UnclasG3/89 13601

Smithsonian Astrophysical ObservatorySPECIAL REPORT 363

https://ntrs.nasa.gov/search.jsp?R=19750011023 2020-04-30T17:53:08+00:00Z

Research in Space ScienceSAO Special Report No. 363

THE ESTIMATION OF 550 km X 550 km

MEAN GRAVITY ANOMALIES

M. R. Williamson and E. M. Gaposchkin

February 20, 1975

Smithsonian InstitutionAstrophysical Observatory

Cambridge, Massachusetts 02138

408-085

TABLE OF CONTENTS

Page

ABSTRACT................................... ... v

1 THE BLOCK COVARIANCE METHOD .................... 1

2 THE 10 X1 0 DATA ......... .......................... 3

3 THE BLOCK GRAVITY ANOMALY ESTIMATES ................ 5

4 REFERENCES... ......................... .. . . 19

FIGURES

1 Distribution of 10 X 10 mean surface-gravity data .............. 4

2 The block covariance function of unit gravity anomalies ........ 7

TABLES

1 Comparison of 10 X 1° mean gravity anomalies withDMAAC (1973) ........... .... .................... 3

2 The block covariance function of unit gravity anomalies ........ 8

3 The block covariance matrix ....................... * 9

4 Estimated block gravity anomalies referred to the 1967 system... 10

iii

ABSTRACT

The calculation of 550 km X 550 km mean gravity anomalies from 1° X 1" mean

free-air gravimetry data is discussed. The block estimate procedure developed by

Kaula is used. Estimates for 1452 of the 1654 blocks are obtained.

RESUME

Dans cet article on discute du calcul des anomalies de la gravite

moyenne sur 550 km x 550 km a partir de donnges de gravim6trie moyenne

sans air sur 1o x 10. On utilise le proced6 d'estimation de blocs d6ve-

lopp' par Kaula. On obtien des estimations pour 1452 des 1654 blocs.

KOHCHIEKT

06cyHaeT F BbuIYFtLeHMe cpe ,HMx rPpaBTaOHHbx aHOMa.EnH

550KM X 550KM, 1CXOZF X3 cpezHiX zaHHblX rpaBMToMeTpM CBB60HC

aTMoc epb 10 X 10. YnoTpe6nFreTc npouegypa oueHK 6 no.a pa3pa-

6OTaHHaH Kayna. HoinyMeHbl OeHKM nAJI 1452 x3 1654 6jOKOB.

THE ESTIMATION OF 550 km X 550 km

MEAN GRAVITY ANOMALIES

M. R. Williamson and E. M. Gaposchkin

1. THE BLOCK COVARIANCE METHOD

The objective is to obtain estimates of mean gravity anomalies for 550 km X 550 km

regions from an incomplete set of 10 X 10 mean free-air gravity anomalies. The esti-

mate procedure is based on the covariance analysis developed by Wiener (1966),

Kolmogoroff, and Kaula (1967). A description and evaluation of the procedure are

given in Gaposchkin (1973).

To obtain the mean gravity anomaly for a 550 km X 550 km block, estimates of all

110 km X 110 km unit mean gravity anomalies in the block are needed. For units

where measurements of 1 X 10 mean gravity anomalies are available, the unit mean

anomalies are taken to be the average of the 10 X 10 mean anomalies within the unit.

The other unit mean anomalies g are calculated from the measured unit anomalies f

in the same block by the equation

N

gi Kik fjk

Kjk= K(f, Tjk)

K is the intra-block covariance function of the unit anomalies f, and Tjk is the distance

between the jth and kth units. The covariance function is estimated from

This work was supported in part by Grant NGR 09-015-002 from the NationalAeronautics and Space Administration.

1

K(f, T7)= 1 E f fk ,jk jk

where the sum includes Njk pairs of measurements with

AT ATT-- < Tjk < T + -

The block anomalies are obtained by averaging the unit anomalies. The mean anomalies

for blocks that include no measured 10 X 10 mean anomalies cannot be estimated by this

procedure.

2

2. THE 10 X 10 DATA

We obtained four sets of 1l X 1l mean free-air gravity anomalies. A set of 29, 209

measured 10 X 10 means was obtained from the Defense Mapping Agency Aerospace

Center (DMAAC) (1973). A set of 1454 10 X 10 means for Australia was obtained from

Mather (1970). Two sets of 10 X 10 means were obtained from Talwani: 4250 10 X 10

means for North America and the North Atlantic (Talwani, Poppe, and Rabinowitz,

1972), and 3944 10 X 10 means for the Indian Ocean (Kahle and Talwani, 1973). The

sets of data were combined, the Mather and the Talwani data were used for regions

when available. The combined set has measured values for 31, 654 of the 64, 800 10 X 10

areas. Figure 1 is a map showing the distribution of the data.

The DMAAC data were compared with the other sets of data and with a previous

compilation by the Aeronautical Chart and Information Center (ACIC) (1971) at common

points. Table 1 shows the results of these comparisons.

Table 1. Comparison of 1° X 10 mean gravity anomalies with DMAAC (1973).

Number of Mean difference RMSData points compared (mgal) (mgal)

Australia 1364 1. 64 24. 16(Mather, 1970)

North America and 3613 -0.18 15.29the North Atlantic(Talwani et al., .1972)

Indian Ocean 2226 -1. 66 23.09(Kahle and Talwani,1973)

Worldwide 19, 164 -0. 23 16. 99

(ACIC, 1971)

3

z-~ _____ --~-- ---. ----------------- 41 g~P~ - ~----r

F-~~---~- --~--~2~2=-~--~ -- ----

i;F_5=-~ ~ ___-Bp "-'- .I ____ ~ _ g t--- ___~

._-1

_ I~_ I-~_,~ ~rlT~__ I..r~~I __

_-~-~--' ----- _ :_1- _ r-_;~~_-~_~ -"-

IP - _-~:1

Ir-~ I . _L_=:r =, ~~------~1~

_ I_ _ ~___

'----:_ ..- =-~-1~_

_

-_-_f

Figure i. Distribution of 1" X 1" mean surface-gravity data.

3. THE BLOCK GRAVITY ANOMALY ESTIMATES

The estimates were obtained with the procedure described in Section 1 by using

the composite data set described in Section 2.

The reference ellipsoid used in the calculation is the Geodetic Reference System

of 1967 (International Union of Geodesy and Geophysics, 1967). A reference ellipsoid

is determined by the geocentric gravitational constant GM, the rotation rate of the

earth w, the semimajor axis of the earth a, and the dynamical form factor J 2 . The

1967 system is defined by the ellipsoid

20 3 -2GM = 3. 98603 X 10 cm sec ,

-5 -1w = 7. 292 115 1467 X 10- 5 rad sec - 1

a = 6. 378160 X 108 cm ,

J2 = 0.001 0827

From these quantities, the flattening f is determined to be

f= 1/298. 247 167 427

A formula for the reference gravity y as a function of latitude 4 is

1l967 = 978031. 85 (1 + 0.005 278 895 sin2 4 + 0.000 023 462 sin4 ) mgal

This formula has an accuracy of 0. 004 mgal. Another formula, which has an accuracy

of 0. 1 mgal, is

y 67 = 978031. 846 (1 + 0. 005 302 86 sin2 O + 5.82 X 10 - 6 sin2 2 ) .

The ACIC data were compiled using this reference system and the International Gravity

Standardization Network, 1971 (Morelli and Gantar, 1974). The other data, however,

are referred to the 1930 international gravity system and the Potsdam Network. To

5

convert these data to the 1967 system, a correction A must be added to the anomalies.

The 1930 international gravity system is defined by the reference gravity

Y1930 = 978049.0 (1 + 0.005 2884 sin 2 + 5. 9X 10- 6 sin2 2 ) mgal

This reference gravity corresponds to an ellipsoid defined by (Heiskanen and Moritz,1967)

20 3 -2GM = 3. 986 329 X 10 cm sec - 2

W = 0.729 211 510 X 10- 4 rad sec -1

a = 6. 378 388 X 108 cm

f = 1/297 .

With the inclusion of the Potsdam correction of 14 mgal, the correction A is

A= 1 9 3 0 -Y 19 6 7 - 14

To an accuracy of 0. 1 mgal

A= 3. 14 - 13. 58 sin2 + 0.02 sin2 24 mgal

This formula was used to convert the 10 X 10 mean anomalies to the 1967 system. Ifthe anomalies are to be combined with satellite data, the effect of the atmosphere mustbe included. The surface gravity given by the anomalies is increased by 0. 87 mgal toobtain the total gravitational force that influences satellite motion. This correction hasnot been made to the anomalies given in this report.

The covariance function is given in Figure 2 and Table 2. The covariance matrixis given in Table 3. The composite set of 31, 654 10 X 10 mean anomalies provided23, 777 of the 41, 350 unit mean anomalies. From these unit mean anomalies, the esti-mate procedure provided 1452 of the 1654 block anomalies. The block anomalies referredto the 1967 system are given in Table 4. A punched card deck can be provided onrequest.

6

408-085

1000 I I

800

N 600

E

400

200 -

I I I I I II 2 3 4 5 6 7

DEGREES

Figure 2. The block covariance function of unit gravity anomalies.

We compared these block anomalies with 1183 block anomalies calculated previously

(Gaposchkin, 1973). The mean difference is 0. 6 mgal and the RMVIS difference is

10 mgal.

7

Table 2. The block covariance function of unit gravity anomalies.

Average angular Covariance functiondistance (mgal2 )

00 948

0.29 1614

0.94 607

1.21 457

1.79 377

2.18 316

2.81 279

3.16 252

3.69 243

4.20 229

4.77 182

5.22 161

5.67 223

6.20 118

6.69 201

8

Table 3. The block covariance matrix.

1 2 3 4 5 6 7 8 9 10 11 12 13

1 9. 4826E+02 5.7175E+02 3.4442E+02 2.6393E+02 2.3447E+02 5.7175E+02 4.2864E+02 3. 1291E+02 2.5177E+02 2.3121E+02 3.4442E+02 3. 1291E02 2.7692E02

2 5.7175E+02 9.4826E+02 5.7175E+02 3.444E+02 2.639344E+02 4.2863E+02 57175E+02 4.2863E+02 3.1291E+02 2. 5177E+02 3. 1291E+02 3.4442E02 3.1291E+02

3 3.4442E+02 5.7175E+02 9.4826E+02 5.7175E+02 3.4442E+02 3. 1291E+02 4.2863E+02 5. 7175E+02 4.2863E+02 3. 1291E+02 2.7692E+02 3. 1291E+02 3. 4442E+02

4 2.6393E+02 3.4442E+02 5.7175E+02 9.4826E+02 5.7175E+02 2.5177E+02 3.1291E+02 4.2863E+02 5.7175E+02 4.2863E+02 2.4411E+02 2.7692E+02 3. 1291E+02

5 2.3447E+02 2.6393E+02 3.4442E+02 5.7175E+02 9.4826E+02 2.3121E+02 2.5177E+02 3. 1291E+02 4.2863E+02 5.7175E+02 2.0662E+02 2.4411E+02 2.7692E+02

6 5.7175E+02 4.2863E+02 3. 1291E+02 2. 5177E+02 2.3121E+02 9.4826E+02 5.7175E+02 3.4442E+02 2.6393E+02 2. 3447E+02 5.7175E+02 4.2863E+02 3. 1291E+02

7 4.2863E+02 5.7175E+02 4.286.3E+02 3. 1291E+02 2.5177E+02 5.7175E+02 9.4826E+02 5.7175E+02 3.4442E+02 2.6393E+02 4.2863E+02 5.7175E+02 4.2863E+02

8 3.1291E+02 4.2863E+02 5.7175E+02 4.2863E+02 3.1291E+02 3.4442E+02 5.7175E+02 9.4826E+02 5.7175E+02 3.4442E+02 3. 1291E+02 4.2863E+02 5.7175E+02

9 2.5177E+02 3.1291E+02 4.2863E+02 5.7175E+02 4.2863E+02 2.6393E+02 3.4442E+02 5.7175E+02 9.4826E+02 5.7175E+02 2.5177E+02 3.1291E+02 4.2863E+02

10 2.3121E+02 2.5177E+02 3. 1291E+02 4.2863E+02 5.7175E+02 2.3447E+02 2.6393E+02 3.4442E+02 5.7175E+02 9.4826E+02 2.3121E+02 2. 5177E+02 3. 1291E+02

11 3.4442E+02 3.1291E+02 2.7692E+02 2.4411E+02 2.0662E+02 5.7175E+02 4.2863E+02 3. 1291E+02 2.5177E+02 2.3121E+02 9.4826E+02 5.7175E+02 3.4442E+02

12 3.1291E+02 3.4442E+02 3. 1291E+02 2.7692E+02 2.4411E+02 4.2863E+02 5.7175E+02 4.2863E+02 3. 1291E+02 2.5177E+02 5.7175E+02 9.4826E+02 5.7175E+02

13 2.7692E+02 3. 1291E+02 3.4442E+02 3. 1291E+02 2.7692E+02 3. 1291E+02 4.2863E+02 5.7175E+02 4.2863E+02 3. 1291E+02 3.4442E+02 5.7175E+02 9.4826E+02

14 2.4411E+02 2.7692E+02 3. 1291E+02 3.4442E+02 3.1291E+02 2.5177E+02 3. 1291E+02 4.2863E+02 5.7175E+02 4.2863E+02 2.6393E+02 3.4442E+02 5.7175E+02

15 2.0662E+02 2.4411E+02 2.7692E+02 3.1291E+02 3.4442E+02 2.3121E+02 2.5177E+02 3. 1291E+02 4.2863E+02 5.7175E+02 2. 447E+02 2.6393E+02 3.4442E+02

16 2.6393E+02 2.5177E+02 2.4411E+02 2.2541E+02 1.7123E+02 3.4442E+02 3. 1291E+02 2.7692E+02 2.4411E+02 2.0662E+02 5.7175E+02 4.2863E+02 3.1291E+02

17 2.5177E+02 2.6393E+02 2.5177E+02 2.4411E+02 2.2541E+02 3.1291E+02 3.4442E+02 3.1291E+02 2.7692E+02 2.4411E+02 4.2863E+02 5.7175E+02 4.2863E+02

18 2.4411E+02 2.5177E+02 2.6393E+02 2.5177E+02 2.4411E+02 2.7692E+02 3.1291E+02 3.4442E+02 3.1291E+02 2.7692E+02 3. 1291E+02 4.2863E+02 5.7175E+02

19 2.2541E+02 2.4411E+02 2. 5177E+02 2.6393E+02 2.5177E+02 2.4411E+02 2.7692E+02 3. 1291E+02 3.4442E+02 3. 1291E+02 2.5177E+02 3. 1291E+02 4.2863E+02

20 1.7123E+02 2.2541E+02 2.4411E+02 2.5177E+02 2.6393E+02 2.0662E+02 2.4411E+02 2.7692E+02 3.1291E+02 3.4442E+02 2.3121E+02 2.5177E+02 3.1291E+02

21 2.3447E+02 2.3121E+02 2.0662E+02 1.7123E+02 2.2093E+02 2.6393E+02 2.5177E+02 2.4411E+02 2.2541E+02 1.7123E+02 3.4442E+02 3. 1291E+02 2.7692E+02

22 2.3121E+02 2.3447E+02 2.3121E+02 2.0662E+02 1.7123E+02 2.5177E+02 2.6393E+02 2.5177E+02 2.4411E+02 2.2541E+02 3. 1291E+02 3.4442E+02 3.1291E+02

23 2.0662E+02 2.3121E+02 2.3447E+02 2.3121E+02 2.0662E+02 2.4411E+02 2.5177E+02 2.6393E+02 2.5177E+02 2.4411E+02 2.7692E+02 3. 1291E+02 3.4442E+02

en 24 1.7123E+02 2.0662E+02 2.3121E+02 2.3447E+02 2.3121E+02 2.2541E+02 2.4411E+02 2.5177E+02 2.6393E+02 2.5177E+02 2.4411E+02 2.7692E+02 3.1291E+02

25 2.2093E+02 1.7123E+02 2.0662E+02 2.3121E+02 2.3447E+02 1.7123E+02 2.2541E+02 2.4411E+02 2.5177E+02 2.6393E+02 2.0662E+02 2.4411E+02 2.7692E+02

14 15 16 17 18 19 20 21 22 23 24 25

1 2.4411E+02 2.0662E+02 2.6393E+02 2.5177E+02 2.4411E+02 2.2541E+02 1.7123E+02 2.3447E+02 2.3121E+02 2.0662E+02 1.7123E+02 2.2093E+02

2 2.7692E+02 2.4411E+02 2.5177E+02 2.6393E+02 2.5177E+02 2.4411E+02 2.2541E+02 2.3121E+02 2.3447E+02 2.3121E+02 2.0662E+02 1.7123E+02

3 3.1291E+02 2.7692E+02 2.4411E+02 2.5177E+02 2.6393E+02 2.5177E+02 2.4411E+02 2.0662E+02 2.3121E+02 2.3447E+02 2.3121E+02 2.0662E+02

4 3.4442E+02 3.1291E+02 2.2541E+02 2.4411E+02 2.5177E+02 2.6393E+02 2.5177E+02 1.7123E+02 2.0662E+02 2.3121E+02 2.3447E+02 2.3121E+02

5 3. 1291E+02 3.4442E+02 1.7123E+02 2.2541E+02 2.4411E+02 2.5177$+02 2.6393E+02 2.2093E+02 1.7123E+02 2.0662E+02 2.3121E+02 2.3447E+02

6 2.5177E+02 2.3121E+02 3.4442E+02 3.1291E+02 2.7692E+02 2.4411E+02 2.0662E+02 2.6393E+02 2. 5177E+02 2.4411E+02 2.2541E+02 1.7123E+02

7 3.1291E+02 2.5177E+02 3.1291E+02 3.4442E+02 3.1291E+02 2.7692E+02 2.4411E+02 2.5177E+02 2.6393E+02 2.5177E+02 2.4411E+02 2.2541E+02

8 4.2863E+02 3.1291E+02 2.7692E+02 3.1291E+02 3.4442E+02 3.1291E+02 2.7692E+02 2.4411E+02 2.5177E+02 2.6393E+02 2.5177E+02 2.4411E+02

9 5.7175E+02 4.2863E+02 2.4411E+02 2.7692E+02 3.1291E+02 3.4442E+02 3.1291E+02 2.2541E+02 2.4411E+02 2.5177E+02 2.6393E+02 2.5177E+02

10 4.2863E+02 5.7175E+02 2.0662E+02 2.4411E+02 2.7692E+02 3.1291E+02 3.4442E+02 1.7123E+02 2.2541E+02 2.4411E+02 2.5177E+02 2.6393E+02

11 2.6393E+02 2.3447E+02 5.7175E+02 4.2863E+02 3.1291E+02 2.5177E+02 2.3121E+02 3.4442E+02 3.1291E+02 2.7692E+02 2.4411E+02 2.0662E+02

12 3.4442E+02 2.6393E+02 4.2863E+02 5.7175E+02 4.2863E+02 3. 1291E+02 2.5177E+02 3.1291E+02 3.4442E+02 3.1291E+02 2.7692E+02 2.4411E+02

13 5.7175E+02 3.4442E+02 3.1291E+02 4.2863E+02 5.7175E+02 4.2863E+02 3.1291E+02 2.7692E+02 3.1291E+02 3.4442E+02 3.1291E+02 2.7692E+02

14 9. 4826E+02 5.7175E+02 2.5177E+02 3. 1291,+02 4.2863E+02 5.7175E+02 4.2863E+02 2.4411E+02 2.7692E+02 3.1291E+02 3.4442E+02 3.1291E+02

15 5.7175E+02 9.4826E+02 2.3121E+02 2.5177E+02 3.1291E+02 4.2863E+02 5.7175E+02 2.0662E+02 2.4411E+02 2.7692E+02 3.1291E+02 3.4442E+02

16 2.5177E+02 2.3121E+02 9.4826E+02 5.7175E+02 3.4442E+02 2.6393E+02 2.3447E+02 5.7175E+02 4.2863E+02 3.1291E+02 2.5177E+02 2.3121E+02

17 3.1291E+02 2.5177E+02 5.7175E+02 9.4826E+02 5.7175E+02 3.4442E+02 2.6393E+02 4.2863E+02 5.7175E+02 4.2863E+02 3.1291E+02 2.5177E+02

18 4.2863E+02 3.1291E+02 3.4442E+02 5.7175E+02 9.4826E+02 5.7175E+02 3.4442E+02 3.1291E+02 4.2863E+02 5.7175E+02 4.2863E+02 3.1291E+02

19 5.7175E+02 4.2863E+02 2.6393E+02 3.4442E+02 5.7175E+02 9.4826E+02 5.7175E+02 2. 5177E+02 3.1291E+02 4.2863E+02 5.7175E+02 4.2863E+02

20 4.2863E+02 5.7175E+02 2.3447E+02 2.6393E+02 3.4442E+02 5.7175E+02 9.4826E+02 2.3121E+02 2.5177E+02 3.1291E+02 4.2863E+02 5.7175E+02

21 2.4411E+02 2.0662E+02 5.7175E+02 4.2863E+02 3.1291E+02 2.5177E+02 2.3121E+02 9.4826E+02 5.7175E+02 3.4442E+02 2.6393E+02 2.3447E+02

22 2.7692E+02 2.4411E+02 4.2863E+02 5.7175E+02 4.2863E+02 3.1291E+02 2.5177E+02 5.7175E+02 9.4826E+02 5.7175E+02 3.4442E+02 2.6393E+02

23 3.1291E+02 2.7692E+02 3.1291E+02 4.2863E+02 5.7175E+02 4.2863E+02 3.1291E+02 3.4442E+02 5.7175E+02 9.4826E+02 5.7175E+02 3.4442E+02

24 3.4442E+02 3.1291E+02 2.5177E+02 3.1291E+02 4.2863E+02 5.7175E+02 4.2863E+02 2.6393E+02 3.4442E+02 5.7175E+02 9.4826E+02 5.7175E+02

25 3.1291E+02 3.4442E+02 2.3121E+02 2.5177E+02 3.1291E+02 4.2863E+02 5.7175E+02 2.3447E+02 2.6393E+02 3.4442E+02 5.7175E+02 9.4826E+02

a. 0~. 0 4n. o- N00 N T 0-. -t. aN 4 0 ~ ~ l a

.4N.MN ~ ~ I .NM .4 . *... N.CJ t . 55 SN: .NN .. N

4.44444 4 ,4 ttros 44,O444044 4 4 44 4 4

4 4 4 4*4 t 44 4 .4 .t

0 a)NNNNN wNNN 0 NmMM

o2 2u:t!2a :5 m4 ""

o4-4

Nl 4 r 4 4 9N 2 r- W~~4N4 4 0rU- 0 w04 0,4 Nm %W NU LlflWf

o0 -4 N0NN0 N N NM 4-o - NNNNN0NNNN ..N ..NN gN N.-.N--g. .NN .N .

o 0.M ~ .4 . 40 0 MIO0 M0MN N y m em N-

0-to 4 0-4 M-tMC~aCN4 N'N No a. 'aMM C 00 M MQM.-

- 0N.l. M -I NM N N-N -N M -* .. N- *

lrlnn ~ ~ ~ ~ ~ ul wlSnninunV~n n n11n nnCbM ;*,

-14 4 Q 4 0 10 4 404 IC 0- 0 4 4 z0NQMa4 a 4 0 M 4zj.0 @ O z 4 4MN-.-0

NM4 400N @ 0N'umm flo ' 9 n -49 0-4 zm m2P 4 ~ 0 r mCNwM P r l0

4w rC -W m ~ 4,C~C a,~A~U 10 ? phnAU0

m ! 'r .M4 mOM 0S m4 l.-M 0N NZ~O -. 4 . 0 4 4' A 0 -

0 .O ,4 0.0a ".aM M N-M4. 0MC NM NaM. 0 N 4NIM4MNNE- l - -MO-

r:- r0n C wMa 0 4~ 4- ~ - ~ . N '' 4 N N.0 ~ NM-mzC -4

ll II i ri'C ,1 IliI Ii1 LnIin I l l wn i nlniIn ,n Ili I IB m li

4PCC. rOO 4-V. 9-CU4 O

NNJJNNN N NNNNNNNNNJJNNJNNJJJJJJNN OwJNN m mm mm mm mw m mm mm mm w mNJmm m

..................................

VV O IV V I 0V,-JO.O J-*N C .J V.JV.OVO J ,ZNJ.)V.Z Z.JZ 4r' O Z; W 'w

S.S 1S 14 .- 4 -4 -4 .S Z -&4 1r 4 -4 '9 'r 4 41 4 WC J, 4O-N CC .S- r

O O O O O O O O C O C O . J.O.O N.C.C.C O.C.C.O.O C. ... .

NN NN~ JN N NJJJJ NNN N N NNNNN J N. JJJNNNNNNNNN

V.V.V.~~NJC W...t. . . .h . . . . . U I 1V .~ . . .V.... . ..... .I. .

NNNNNJNNJNNJNJJNJNNJJNJNNJNNNJJNJNNJNNJNNJNJJNJJNJJNJNNJNNJNNJNJJNJJNJJNJNNJ -J-

oz OV 0 -0 ' V JJ'U'.'V -W

t.IWWWWW wJAJ mmNTJNJJ~rNNN b~aIA

. . r 9 WN N . .o .

0 1 1O0 0r 4 Ontrotr rt t 0 C V.C o NV.CVw . C.9- rUCtt w r " rV.CItri

wm IV IWf~ cc Mccl F

NJ c c -I' I w. Ii INNN wN.NNc~ I-.- vIAN0

z

NNNNNNN NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN~ @

WWWWWWWWWNNNNNNNNNNNNNN

004 ON- N0

00000U 000000000 000000000 0 0 0 000 4a 000000000 t 0000

NNNNNNNNNNNNNNNNN~~~NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN a

WNS.,.' N. I~ *-.I. I ................ ...... NI . .... .u ..

0 '-J 0. 0 W4'W4'U a, 0 0.' - I 0' 4' 4'J 0 ON0o-0 0' 4' 0 4' o 0.-. 0 ' 0-U 40 WNUO' 0o l 0 0t. 0 -*U0 . '0 .JNOW WN

0.0.00W.N N 0O' '0O O'N N.I O0' 'O' 0NW 'N' 00N00004'O NO'O NOtn7'n '

NNNNNNNNNN NN NNNNNNN N- NNNNN0NNNNN N NNN

WW UWW N NN N NN N NNNrW

04'OW-J NODta'0'ON J 0'04 .W N O' t.' f ON -O O0' 04' OW -.4 O' .. 0'O - 0 0'0U)- O'O U. '04

o'00 0O' 0000 000 00'0000NOWOOoN00O'0 0O' oO00000'0000'000,040O'W0

NONNNO~NNNONNNNNONNNNNNO~NNNNONONNNNONNNNN0NONNN

- *-.NN~~NNNNN lNN r~crrrr II 1a- NINNI 'N II

4 - o , ;77Z 0 , 4 -JNN - 4W .O 04 0 ''.W. Ww~ ' tn 40WO O1 4' 4'-OWnn ~ t W N 't

?:"; :1,:- np f :D f f.

-' 00 N 04'OO4'0 N 00 4 'W40 O 0.04 N 00 N O N n4

04' 0000- 0000 4C0-400000 .0 4N00 0' 00-N 00-Nt' O 00-4 t0000400 0000 U' W0 0 4' 0 0, O'p

NNNNNNN""~" NNNNNN---NN--------NNNNNNNNNNNNN----NN 0

NNNNNNNNN NNNNN WWWZ4WWN 0111cc 1 cIm c C

OOO O't n4 4 '

W NO C.0 00U .U.4'4'WWD N N0 C 0 ' ' ~WN U.U.4'4'W 00

........ ... "..;.f r..............................orr,,

NNNNNNNNNNNNNNNNNNNNNNNNNNNNNN rrrNNNNNNNNNNNNNN rNrNN0

o o o o o o o o o 0 0 0 0 0 0 0 OO, 0 N IN I

'00 0' 0'00 00 0'00 00 0'0 .000. '.000 0'.00 0 0' 0 c

00 0000000000000000000000 000000000000000000000000000000 NONON5

am. a a aim a a. cam. sal a ~~NNNN

NNN NN NNNN - - N -.- NNNN NNNNN .NNN *-'NNN NNN

D~0U'LU.U'L~YUrJJ * UtW 0'OCU0'4'4NU'WN .0JOOU'%)4'U'~ ~ ~ YUUC N CrJrU 00C

o0- a01 0 U 0 0Q 0

4WWO.W1,DWW0- W0 4; O'WZW4 OWOW WWW0 WWW 4WM WO WW 00 -. 0-40-NN .O 0.4CDN WW 0),0WW0WWWNW0WWW0WWWW0 o0W100WW W4 W 000 NO0-00000000.000

muma * m u liii .1 IS44

f f fr f f tf f f;D tl t

Nl~UJ W a II - -1 11 WNNNNiNNN-N -I~iN

u!4o oo rt n 40 440t 0,N ONtr w oo U'OtflU, O X4U1M OU-trOZ~t .

000 w 0 0 0 000000w00W-41"000OO~o O~ 0wo N 00 WWOW 00OW 00000000WW 0OW 4

oooooo o~o OOO OO OOO OOO OOO OOO O .P')NiJI ~iU4.

IIIi alal0l .-I r f Nil f f f fi- .- WW ffI O F. .DNJ~iN)i I i

~ S .0 5 U N itN O tN NW N N 0 - "

. 0 00-J i 4 0. 0..-i4iJ W0 S.45 J.- 5

4= .i0 O 0 W W S 0 0 0 N i N . 0 0 0 0N * -N - W"' L

OW 0. fW 00 W0 00 W i0 WOf rv0W4Pf0000000 0000000 !T ;fVIp !: el:P

0 0N.- , 0N 4-mJ4 4 !i: Ni ii0 ~ NN )i N NJ N )

MCC- W 9- W W wNN i~N NNci i -- - -carc c Z 4

WaW W UJ W W W WJU~ WW WODW~N~JN4JNJJN WO MW WWWWW WO DW J

N~O~ 0 ~a~J w0 0 ; 1 40WWV 4 - WW 0

=WW W2~ W W=W =0WW -4W~ t WW W W W~ WWW W0 WW 4tW

mmmmtli m mamum mD

t*zmt l NN~JNNIN ~ aJi 1

WWW~WWDW WWW NW W U W WW WW NW WW WW WW WW WWODW W W WW C t W',

IjO O O O O~'0 .O.O N I IO O~ IW J'J~0 I~J- I. O* I 'III

N -1NNNNNJNNN NJJNNJJN IJNNNNNJ- - --

Tf 12 vrva v f; ;

mtigmummz mzN N0 0 41 0NNJJNNJJNN N N 4-4JNWNWNN-NJNJ MJNNJNJNJ"ZNNWNJNJ ZJN NJNJ 0JNJNJNJ4JV=NJNNNNJJJJJNNN -jS M-4

000 U) OU) w a ww w wza U)U U U ) JNJN N J JNJN N J JNJN

C aNJNL .~

N~I WWI ) I ) J~ 1 -, 1a -- ~ I S I NJ*- I "I I NJN

z

WWWNNNNNNNNNNNMNWtUN NODN NMNNMNCWNNWWWWWONNNNNNNNNNNNNNMOD M 0.

0000 050 000 000 0 S 000 000 000 000 00000000t000000000000t -

-4 -4- r 4 t 000 Ut Ut4 Ut1 -400 4- 4'4'4.0.N0- 1 14-4 r 414 *-'-t- *'*.0 4Ut 40000-4000000000N400000000000 4

------------- UUYYN--- a NNNNN NNNNN NNNN~rO~~~D U E Z

0 00 M0 0 N0 W 0.0 W 0.00 M00 W0 SN4- SC WOt SO 0- M-UtO t ttS-. t tt M 0t tt t 0 0 0t t0 M tt0

NNNNNNNNNNNNNNNNNNNNNNJN -4 09~~f f ff f.~LPE f

0000 C00O-.00400.WN.-000000 0000.W .0 .0N M0410--; w 0t)Wr1 7" 7 ;;::W- ! N N - N .. ." " 11

.....................................tbf aI

a. . . . . . . . . . . . . . . . . . . . . . . . . . . . . .rb 0000000Z000040000000000 00 000000000000000000000 0 -4 .

aaa OW*-UttNN.-.-. am----------------aNNNax.NNN .a aNNN aa am-ma

Ut0-4ONNNN WN N OO NN S O W OO O 00 O OO 00 . NO..J N0.-W000000000.-.Ut,-.Or tr t0UtNOWOSO0NNN NOSOO 40.jO .-. j~ -00

NNNN t~.-. - - - - - - - - - NNNNNNN N N N N N NN - -. NN .- 0

...... ~~~... ..

~~~NNNNNNNNNNNNNNN rcrrrrrrr WWWW

tOO00Ut~~t55WWNN...-.00 O--OUSWN..-0t O-4 Ut 5*W .. t. UtU'SSWWN' C @0~0O OQQE E W m h

.."""~~PfPPP~~~~m " ''P ~:..........

I~~~~ ~~~NN N N~:= ~NNN NN N N N N N NNN NN NN NNNN N N N N N N NNNNlNNNNNNNNNNNNN NNNNN @~ Ui

00000000000000000000000000000000000000-0000-00000 SN

...... t , ... ~"."... ~~~~ ......

t-.O CWS WCm O SCOWW45-Urm~IUtOWOOtCONON-00'.. CS N N - W S W N4 ~f0 0 -00-UtOOO tOtONOS~O.~~WOUW 0 5Ut00 00 00 0 0 0 0 0 0 S4 N 5 0 -J Ot.

-- - - - - - - - - - - - - -NNNNNNNNNNNNNNNNNNN-.--NN N N -£CSW@0Ut0'Ut0'.-StOC Y~tUY~tU~UU ~~UUt~ttf t STU -tS.-UtuJt rt UtLW

ee.Oeo o oeo oo 0 010 01010oo z 1 lzOQ 0 010 oeoo z ww owwwwow 4,

!I i lW-8 I I I lS iII a ii l -- o111

N ~ ~ ~ ~ ~ ~ 4 N~ Q N NNNNNNNNN N N N N: N N N4 N. N N N N Nl o N &

NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN

LMUUJNNNNNN WtoWWJrrrrrrtoWN, oX

w ot 4S 0on a, t o Swo&oo 4ooozoo o o too t oo

• "N" NN" NN"N"NNN" NN NNNNNNNNNNNNNNNNNN NNNNNNNNNN 0." %" ""QOQOQOOOQOQQQQ 0QQQOQOQ0QQQOQOQQQOQOQWlJWWt.YWWOWU 0

O0 0 0 0 0 0 00 0

m0 0 0 0 0

O0

0000000C0J0004 WWWONWWW0UW

S. ............. .......... ... .

NWJJ~DVNWNOfrm00N00QtOOWN0NWQW(WII NO0 *.NW0NN~r 50

;44 V. .~ . . . . .

NN N NNNNNNN .. NNNNNNNNNNN. . NN. .. N ZNN W6~~EN WO W40X0W a~toototot~WOato~ttootto WOOW~tW- P

NNNNNNNNNNNNNNNNNNNNNNNNNJNNNNNNNN NNNNNNNNNNNNN ;

IUWt.tWtW WWWW WNN NN~CNNN N r tcrrrLc t.ZUJUJ

litilililllil C'COOlllOillllNlrrOQimillllillWilir ilNN

ca)oooo oooooooo o ooooo 0oo00ooo oo oo ooo oo 'ooo

w 0 00." "WNWNNWNNNNW ONWNOW NNWNWNmNWNNNee NWNONWNWNNWNNWNWNNWN DO° O

all Ia a II alaS afP~T amII *INN~.P?* *-J1 ~ 5~r~~ Irr I I1 W: I~pp .-. NNN I.-.

NtO~Q)WQtPt00I S.J~Otr X.Jt)'SS to. =ED*L.XJ100'.C90O Wt W-NNt Eo

OO0NYO O WXWf ~CN0 Nm00W m'~~WJ~oWfOD.- WtoOW0.,* OWmQW*P oO.*~~ N .-w00000 0N O ~Oa-Jz

- -- NNrNN - - .. NNNNN,-N -.. NNNNNNN~r 0 ~

o~~~toNO& W NWWWWWWWIt0W *.OWtoWW tW.SNC-tStW ttt.-t toWW WW

a WW tWWW WW tN N N N NN NN N N NN N NN rr

totop tWN" OCCTWTh"tototwNNO ~ -4toottaNPtoCWototo WNN 0000

. .. .. .. .. .. .. .. ..

NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN~NNNNNNNNNN 0>~"~~"N""NW~

,,,,,r""""""".".. . . . . . . . . . a

41 4- 4 4 M~ 1 9- 41 41 41 4 2. ' WV M 4' 41 4' 41 42' 2 .' 4'W 4'S~ 4214 4- 4- M' 4- 4.- 4- 4- 2'2'1 414 *W 4- 2' 4- S- 4- g- 4- 2' 4- W S-4 Z' 2 4'II~~" N.1 I 4'-NN W N..NN I--'l NNN

c t wN N s - -w oN N N- N N rN NwN N N N.

0'~~~~ X o- o -t 4rto' '0 O- X t o ' ot o2 t t 't 't 't N00 4 o0 otot ot ot ot

ass s,, 55 t u g g s g as~s g IS alas saga ii ,,,, , as g a g ,,,

4.4.4.4.4.4.4 WW WW;4r:OW W W W WW W W W W J(JUUW W WW W~~Jb-4 -- 4 J N -- P 4- 4 - 4 - J - 4- 4- 4- 4- J- 4- 4- 4- 4 - J- 4 -

uo . 0a o o"11 Ilo 011 uO 10 -o11 01u., Io "1. "nr ,Du u ,

'0W-4j- N.-4 0W OWW.W W4.WW000 W D W-W W,0, WW WW-0WWDWWN0 0W0WOW W. . N .koW W W

5,4 4 4 1S 4N 19 -1 4-a.- . 4- 4-4 -4 4 -4 - -t -4- S -V -- t r 11 41PN rp-g - t5 1- P1

us as Isles sag ass ssss s as",, ts sag ss t os osass s

WWWWNN WW- - - - --W W 2& 0120,20-m- - -WWW&0NN0010N- - - - --41A-0 -41OD41r 0 - 4C4w- 41 4 OU4.-a W FJ C W-J ww .w ww w J aN-w 00QQ Z~-O Ww O ' 00- 0. O C cm. Q p

01 -4 - Z;, - .0 N4.0 0 D00 4 o.0-wJ00 -1 .04. -4ww4 w74 - 0 0

mom -J..-44...- 0 D 01 4J -W - Jw-JO W4-449W-Jww0J0. J WUaE S- mrNDO NUJ" su-44 ED

-O - 1 4- - 1-4- 4 -N 'P -400-4 0 mmm .fl.0w-NP 4. Wo4.E4N4p Mm WV W 40NmJ. 0-W 0 0

-w tW - -1 -4- " -

,uPE0 4.~N0mJO NW 4urn-w 4 410 0 t O U E -w w w w 1D 0OOU,-4N.-'Nw

w~5~ 5 5 ~ g 5s1g~s s 5 5 5 1 a sa g g a sg g a a w us- . 9 40w0 ON(-N0O 4 1 0 - -000 mwwo00000ww 0 0 0 0 00p~ mcE 'Elzmm DUxwJ ED

N N N D N N N N PPO-.J--.J.J.. N z

r~~~~~~runr ~ ~ ~ ~ ~ ~ ~ ~ ~ 9 Droy9lU l f~inUtl ~ ''lflft ,r~r ~ t

all~~~~~~ ~~~ I5 sit Ia~ ssss. sc tssta ass t

000 0 00 !6~J 0 ,-Oc 0 0 0 p0 z

00.- rrJV,, JJ ~ '~J.~

00 00 tO.N0t~ N0 O W. ,tN 0 U 0 0 Jt N 41O i~J' 41 S Ca,0O' )

j0-Jflw w w w.J0 -0 U1 bJO00"0U'O O-4 t-O 1tJO

00 00 0 00 00 00 0 4'1. LU0~C~U 0, 0 0 0v

000 ~ ~ ~ oa & c CC'C. 00 00000 0 wo* U 0O O O

IJJWW ~ ~ 1 -.----J-J.i. v N C J I a. N Ot

4. REFERENCES

AERONAUTICAL CHART AND INFORMATION CENTER

1971. 10 X 10 Mean Free-Air Gravity Anomalies. ACIC Reference Publ. No. 29,

August, 324 pp.

DEFENSE MAPPING AGENCY AEROSPACE CENTER

1973. 10 X 10 Mean Free-Air Gravity Anomalies. DMAAC Reference Publ. No.

73-0002, December, 100 pp.

GAPOSCHKIN, E. M., Editor

1973. 1973 Smithsonian Standard Earth (III). Smithsonian Astrophys. Obs. Spec.

Rep. No. 353, 388 pp.

HEISKANEN, W. A., and MORITZ, H.

1967. Physical Geodesy. W. H. Freeman and Company, San Francisco, pp. 79-

80.

INTERNATIONAL UNION OF GEODESY AND GEOPHYSICS

1967. Geodetic Reference System 1967. International Association of Geodesy

Spec. Publ. No. 3, Bureau Central de l'Association Internationale de

Geodesie, 116 pp.

KAHLE, H. G., and TALWANI, M.

1973. Gravimetric Indian Ocean geoid. Zs. f. Geophys., vol. 39, pp. 167-187.

KAULA, W. M.

1967. Theory of statistical analysis of data distributed over a sphere. Rev.

Geophys., vol. 5, pp. 38-107.

MATHER, R.

1970. The Australian geodetic datum in earth space. UNISURV Rep. No. 19,

Univ. New South Wales, p. 80.

MORELLI, C., and GANTAR, C.

1974. The International Gravity Standardization Net 1971. International Union

of Geodesy and Geophysics, 200 pp.

TALWANI, M., POPPE, H. R., and RABINOWITZ, P. D.

1972. Gravimetrically determined geoid in the western North Atlantic. In Sea

Surface Topography from Space, vol. II, ed. by J. Apel, NOAA Tech.

Rep. ERL 228 - AOML 7-2, pp. 23-1 to 23-33.

19

WIENER, N.

1966. Extrapolation, Interpolation, and Smoothing of Stationary Time Serieswith Engineering Applications. The MIT Press, Cambridge, Mass.,160 pp.

20

BIOGRAPHICAL NOTES

M. R. WILLIAMSON graduated in applied mathematics from Brown University in

1963. She received an M. S. in 1965 and a Ph.D. in 1970 in physics from Tufts

University.

She has been with the Smithsonian Astrophysical Observatory since 1971 as a

mathematician in the Analytical Satellite Geophysics Group. She has developed pro-

grams for analyzing surface-gravity data and for studying the effects of solar radiation

pressure on satellite motion.

E. M. GAPOSCHKIN graduated in electrical engineering from Tufts University in

1957. He received a Degree of Numerical Analysis in 1959 from Cambridge University

in England and a Ph. D. in geophysics from Harvard University in 1969.

Dr. Gaposchkin joined the staff at Smithsonian Astrophysical Observatory in 1959,

where he has been programer and Division Chief of the Computations Division and

mathematician in the Research and Analysis Department. He has helped develop the

basic computer program used in all analyses of satellite motion.

Since 1968, Dr. Gaposchkin has been principal scientist of the Analytical Satellite

Geophysics Group. His main interests include satellite geodesy and geophysics and

applied mathematics. Dr. Gaposchkin spent the last year at Groupe de Recherches de

Geodesie Spatiale at the Observatoire de Meudon, Universite' de Paris.

NO TICE

This series of Special Reports was instituted under the supervisionof Dr. F. L. Whipple, Director of the Astrophysical Observatory of theSmithsonian Institution, shortly after the launching of the first artificialearth satellite on October 4, 1957. Contributions come from the Staffof the Observatory.

First issued to ensure the immediate dissemination of data for satel-lite tracking, the reports have continued to provide a rapid distributionof catalogs of satellite observations, orbital information, and prelimi-nary results of data analyses prior to formal publication in the appro-priate journals. The Reports are also used extensively for the rapidpublication of preliminary or special results in other fields of astro-physics.

The Reports are regularly distributed to all institutions partici-pating in the U. S. space research program and to individual scientistswho request them from the Publications Division, Distribution Section,Smithsonian Astrophysical Observatory, Cambridge, Massachusetts02138.