Global isostatic geoid anomalies for plate and boundary layer models of the lithosphere

Commonly used one dimensional geoid models predict that the isostatic geoid anomaly over old ocean basins for the boundary layer thermal model of the lithosphere is a factor of two greater than that for the plate model. Calculations presented, using the spherical analogues of the plate and boundary layer thermal models, show that for the actual global distribution of plate ages, one dimensional models are not accurate and a spherical, fully three dimensional treatment is necessary. The maximum difference in geoid heights predicted for the two models is only about two meters. The thermal structure of old lithosphere is unlikely to be resolvable using global geoid anomalies. Stripping the effects of plate aging and a hypothetical uniform, 35 km, isostatically-compensated continental crust from the observed geoid emphasizes that the largest-amplitude geoid anomaly is the geoid low of almost 120 m over West Antarctica, a factor of two greater than the low of 60 m over Ceylon

Improving Geoidal Height Estimates from Global Geopotential Model Using Regression Model and GPS Data

Conventionally, for most application, position of a point is often referred to the geoid as the reference surface. Thus there is an important need for the knowledge of the geoid undulation in the area where positioning tasks is performed, This requirement is made more apparent with the advent of high precision using GPS where the resulting ellipsoid height must be converted to orthometric height. An ideal solution is to use a precise gravimetric solution where the geoidal height at each GPS point is computed and applied. Unfortunately, at the moment there is no such solution available in Malaysia. However. efforts are currently being made to develop a precise gravimetric geoid, For the time being, an alternative method would have to be use and the global geopotential model is one of them. [n order to increase the accuracy of computed geoid height from the geopotential model. a regression model is used in conjunction with the GPS data. The resulting accuracy estimates of the geoid height determination increases from around 60 cm 'to about 10 cm leve1

Geoid anomalies and fracture zones in the Pacific Ocean

The high degree and order geoid field in the Pacific is a superposition of fracture zone anomalies and hot-spot swell anomalies. A two-dimensional spectral analysis of this field reveals a very strong north-south wavenumber contribution with a dominant wavelength of about 2000 km, a much smaller contribution from east-west wavenumbers, and negligible contributions from other directions. One dimensional profiles were taken in order to appreciate the magnitudes of the north-south and east-west components. A calculated geoid anomaly using an idealized fracture zone model contains just about the same amount of power in the 2350 km band wavelength as does the north-south profile of the SEASAT geoid field. In an attempt to correlate plate age with geoid anomalies, a digitized age map of the Pacific was used to generate a synthetic geoid, which was subtracted from SEASAT. This procedure produces a residual geoid in which the fracture zone anomalies appear to be diminished, if not removed

Shape of the ocean surface and implications for the Earth's interior: GEOS-3 results

A new set of 1 deg x 1 deg mean free air anomalies was used to construct a gravimetric geoid by Stokes' formula for the Indian Ocean. Utilizing such 1 deg x 1 deg geoid comparisons were made with GEOS-3 radar altimeter estimates of geoid height. Most commonly there were constant offsets and long wavelength discrepancies between the two data sets; there were many probable causes including radial orbit error, scale errors in the geoid, or bias errors in altitude determination. Across the Aleutian Trench the 1 deg x 1 deg gravimetric geoids did not measure the entire depth of the geoid anomaly due to averaging over 1 deg squares and subsequent aliasing of the data. After adjustment of GEOS-3 data to eliminate long wavelength discrepancies, agreement between the altimeter geoid and gravimetric geoid was between 1.7 and 2.7 meters in rms errors. For purposes of geological interpretation, techniques were developed to directly compute the geoid anomaly over models of density within the Earth. In observing the results from satellite altimetry it was possible to identify geoid anomalies over different geologic features in the ocean. Examples and significant results are reported

The geoid spectrum from altimetry

Satellite altimetry information from the world's major oceans was analyzed to arrive at a geoid power spectrum. Using the equivalent of about 7 revolutions of data (mostly from GEOS-3) the power spectrum of the sea surface generally follows the expected values from Kaula's rule applied to the geoid. Analysis of overlapping altimetry arcs (and oceanographic data) shows that the surface spectrum is dominated by the geoid to about 500 cycles (40 km half wavelength) but that sea state departures are significant starting at about 250 cycles (80 km). Estimates of geopotential variances from a derived (smooth) geoid spectrum show significantly less power than Kaula's rule to about 60 cycles, but somewhat more from there to about 400 cycles. At less than 40 km half wavelength, the total power in the marine geoid may be negligible

Analysis of altimetry over inland seas

Satellite-based altimetric data taken by GEOS-3 and SEASAT over the Black Sea and Caspian Sea are analyzed and a least squares collocation technique is used to predict the geoid undulation on a .25-degree by .25-degree grid and to transform these geoid undulations to free air gravity anomalies. This project entailed processing satellite altimeter data over inland seas for recovery of area mean gravity information. Gravity information in this area of the world is not readily available, so the possibility of obtaining it from the processing of altimeter observations is attractive. The principal objective was to complete and extend analyses done in a previous study, verify those results, and document the results and techniques. A secondary objective was to improve the algorithms and results, if possible. The approach used involved editing geoid height data to remove overland data; evaluating geoid height differences at crossover points; removing orbit errors from geoid heights using crossover differences; gridding geoid height data at .25-degree by .25-degree intervals; and estimating the gravity anomalies from gridded geoid heights using the collocation technique

Detailed gravimetric geoid confirmation of short wavelength features of sea surface topography detected by the Skylab S-193 altimeter in the Atlantic Ocean

A detailed gravimetric geoid was computed for the Northwest Atlantic Ocean and Caribbean Sea area in support of the calibration and evaluation of the GEOS-C altimeter. This geoid, computed on a 15 ft. x 15 ft. grid was based upon a combination of surface gravity data with the GSFC GEM-6 satellite derived gravity data. A comparison of this gravimetric geoid with 10 passes of SKYLAB altimeter data is presented. The agreement of the two data types is quite good with the differences generally less than 2 meters. Sea surface manifestations of numerous short wavelength (approximately 100 km) oceanographic features are now indicated in the gravimetric geoid and are also confirmed by the altimetry data

Properties of the lithosphere and asthenosphere deduced from geoid observations

Data from the GEOS-3 and SEASAT Satellites provided a very accurate geoid map over the oceans. Broad bathymetric features in the oceans such as oceanic swells and plateaus are fully compensated. It is shown that the geoid anomalies due to the density structures of the lithosphere are proportional to the first moment of the density distribution. The deepening of the ocean basins is attributed to thermal isostasy. The thickness of the oceanic lithosphere increases with age due to the loss of heat to the sea floor. Bathymetry and the geoid provide constraints on the extent of this heat loss. Offsets in the geoid across major fracture zones can also be used to constrain this problem. Geoid bathymetry correlations show that the Hawaiian and Bermuda swells and the Cape Verde Rise are probably due to lithospheric thinning

Calibration and evaluation of Skylab altimetry for geodetic determination of the geoid

The author has identified the following significant results. The Skylab altimeter experiment has proven the capability of the altimeter for measurement of sea surface topography. The geometric determination of the geoid/mean sea level from satellite altimetry is a new approach having significant applications in many disciplines including geodesy and oceanography. A generalized least squares collocation technique was developed for determination of the geoid from altimetry data. The technique solves for the altimetry geoid and determines one bias term for the combined effect of sea state, orbit, tides, geoid, and instrument error using sparse ground truth data. The influence of errors in orbit and a priori geoid values are discussed. Although the Skylab altimeter instrument accuracy is about plus or minus 1m, significant results were obtained in identification of large geoidal features such as over the Puerto Rico trench. Comparison of the results of several passes shows that good agreement exists between the general slopes of the altimeter geoid and the ground truth, and that the altimeter appears to be capable of providing more details than are now available with best known geoids