ALBERT

Description: Global meteorological analyses from the European Center for Medium Range Weather Forecasts are employed to compute the atmospheric excitation psi of the polar motion for the 9-year period of 1980-1988. Both the matter component psi(matter) and the motion component psi (motion) are computed, the former with and without the oceanic inverted barometer (IB) effect. It is found that psi(motion) contributes significantly to the total excitation psi overall and nonnegligibly to the annual signal in psi, or the annual wobble excitation in particular. The results for the annual wobble excitation, in terms of the prograde component psi(t) and the retrogade component phsi(-) for January 1, are within the (rather large) range of previous estimates. The IB effect has a small impact on psi(+), whereas its impact on psi(-) is considerable.

Description: Temporal variations in the low-degree zonal harmonics of the earth's gravitational field have recently been observed by satellite laser ranging. A host of geophysical processes contribute to these variations. The present paper studies quantitatively a prime contributor, atmospheric mass redistribution, using ECMWF global surface pressure data for the period of 1980-1988. The annual and semiannual amplitudes and phases of the zonal J(l) coefficient with degree l = 2-6 with and without the oceanic inverted-barometer (IB) effect are computed to obtain the predicted effects on the orbit nodal residuals of Lageos and Starlette. These predicted values are then compared with observations. It is found that the atmospheric influence, combined with the hydrological influence agree well with the Lageos observation for the annual term. The corresponding match appears poorer for Starlette.

Description: The GEOS-1 general circulation model has been used to compute atmospheric torques on the oceans and solid Earth for the period 1980-1995. The time series for the various torque components have been analyzed by means of Fourier transform techniques. It was determined that the wind stress torque over land is more powerful than the wind stress torque over water by 55%, 42%, and 80% for the x, y, and z components respectively. This is mainly the result of power in the high frequency range. The pressure torques due to polar flattening, equatorial ellipticity, marine geoid, and continental orography were computed. The orographic or "mountain torque" components are more powerful than their wind stress counterparts (land plus ocean) by 231% (x), 191% (y), and 77% (z). The marine pressure torques due to geoidal undulations are much smaller than the orographic ones, as expected. They are only 3% (x), 4% (y), and 5% (z) of the corresponding mountain torques. The geoidal pressure torques are approximately equal in magnitude to those produced by the equatorial ellipticity of the Earth. The pressure torque due to polar flattening makes the largest contributions to the atmospheric torque budget. It has no zonal component, only equatorial ones. Most of the power of the latter, between 68% and 69%, is found in modes with periods under 15 days. The single most powerful mode has a period of 361 days. The gravitational torque ranks second in power only to the polar flattening pressure torque. Unlike the former, it does produce a zonal component, albeit much smaller (1%) than the equatorial ones. The gravitational and pressure torques have opposite signs, therefore, the gravitational torque nullifies 42% of the total pressure torque. Zonally, however, the gravitational torque amounts to only 6% of the total pressure torque. The power budget for the total atmospheric torque yields 7595 and 7120 Hadleys for the equatorial components and 966 Hadleys for the zonal. The x-component exhibits a large mean value (1811 H), mainly the result of polar flattening pressure torque acting on the ocean surfaces. Atmospheric torque modes with periods of 408, 440, and 476 days appear in the spectrum of the equatorial components.

Description: We conducted numerical simulations of the wind-forcing of the sea level variations in the North Sea using a barotropic ocean model with realistic geography, bathymetry, and boundary conditions, to examine the forcing of the 14-month "pole tide" which is known to be strong along the Denmark- Netherlands coast. The simulation input is the monthly-mean surface wind stress field from the National Centers for Environmental Prediction (NCEP) reanalysis for the 40-year period 1958-1997. The output sea level response was then compared with 10 coastal tide gauge records from the Permanent Service for Mean Sea Level (PSMSL). Besides the strong seasonal variations, several prominent quasi-periodicities exist at around 7 years, 3 years, 14 months, 9 months, and 6.5 months. Correlation and spectral analyses show remarkable agreement between the model output and the observations, particularly in the 14-month, or Chandler period band. The latter indicates that the enhanced pole tide found in the North Sea along the Denmark-Netherlands coast is actually the coastal setup response to wind stress forcing with a periodicity of 14 months. We find no need to invoke a geophysical explanation involving resonance-enhancement of pole tide in the North Sea to explain the observations.

Description: Satellite-based altimetric data taken by GOES-3, SEASAT, and GEOSAT over the Aral Sea, the Black Sea, and the Caspian Sea are analyzed and a least squares collocation technique is used to predict the geoid undulations on a 0.25x0.25 deg. grid and to transform these geoid undulations to free air gravity anomalies. Rapp's 180x180 geopotential model is used as the reference surface for the collocation procedure. The result of geoid to gravity transformation is, however, sensitive to the information content of the reference geopotential model used. For example, considerable detailed surface gravity data were incorporated into the reference model over the Black Sea, resulting in a reference model with significant information content at short wavelengths. Thus, estimation of short wavelength gravity anomalies from gridded geoid heights is generally reliable over regions such as the Black Sea, using the conventional collocation technique with local empirical covariance functions. Over regions such as the Caspian Sea, where detailed surface data are generally not incorporated into the reference model, unconventional techniques are needed to obtain reliable gravity anomalies. Based on the predicted gravity anomalies over these inland seas, speculative tectonic structures are identified and geophysical processes are inferred.

Description: Satellite based altimetric data taken by GEOS-3, SEASAT, and GEOSAT over the Aral Sea are analyzed and a least squares collocation technique is used to predict the geoid undulations on a 0.25 by 0.25 degree grid and to transform these geoid undulations to free air gravity anomalies. Rapp's 180 by 180 geopotential model is taken as the empirical model. The collocation procedure is performed with a set of local residual empirical covariance functions. For comparison, Rapp's global covariance functions and Jordan's self-consistent theoretical covariance functions based on Jordan's formulation and on locally derived parameters are also used to grid geoid undulations and to predict gravity anomalies. The sensitivity of the collocation results to the choice of covariance functions is discussed.

Description: Recent Satellite Laser Ranging derived long wavelength gravity time series analysis has focused to a large extent on the effects of the recent large changes in the Earth s 52, and the potential causes. However, it is difficult to determine whether there are corresponding signals in the shorter wavelength zonals from the existing SLR-derived time variable gravity results, although it appears that geophysical fluid transport is being observed. For example, the recovered J3 time series shows remarkable agreement with NCEP-derived estimates of atmospheric gravity variations. Likewise, some of the non-zonal spherical harmonic coefficient series have significant interannual signal that appears to be related to mass transport. The non-zonal degree 2 terms show reasonable correlation with atmospheric signals, as well as climatic effects such as El Nino Southern Oscillation. While the formal uncertainty of these terms is significantly higher than that for J2, it is also clear that there is useful signal to be extracted. Consequently, the SLR time series is being reprocessed to improve the time variable gravity field recovery. We will present recent updates on the J2 evolution, as well as a look at other components of the interannual variations of the gravity field, complete through degree 4, and possible geophysical and climatic causes.

Description: Via mass transport, the interannual meteorological oscillations (ENSO, NAO, etc.) produce slight variations in the global gravitational field, which can be observed by satellite laser ranging (SLR) and anticipated space missions such as GRACE. The meteorological oscillation modes are identified using empirical orthogonal function and principal component decomposition using de-seasoned surface fields produced from the 40-years of NCEP reanalysis. This decompositional analysis requires an accounting for issues associated with area-weighting and non-zero mean. The atmospheric contributions to gravitational variation, mode by mode and in combinations, are then computed via regional integration. To address the corresponding oceanic contributions, we examine the output of the numerical POCM ocean circulation model (1992-1998) in two ways: ocean bottom pressure field and ocean surface topography minus steric effect. The resulting lowest-degree Stokes coefficients for both NCEP reanalysis and POCM are then compared with space geodetic SLR observations, to identify the importance of each meteorological oscillations in gravitational variation signals.

Description: The oblateness of the Earth's gravity field, 52, has long been observed to undergo a slight decrease due to post-glacial rebound of the mantle. Sometime around 1998 this trend reversed quite suddenly. This reversal persisted until 2001, at which point the atmosphere-corrected time series appears to have reversed yet again. Presently, the time series appears to be returning to the value that would nominally have been reached had the anomaly not occurred. This anomaly signifies a large interannual change in global mass distribution whose J2 effect overshadows that of the post-glacial rebound over such timescales. A number of possible causes have been considered, with oceanic mass redistribution as the leading candidate although other effects, such as glacial melting and core effects may be contributing. The amount by which J2 returns to it's nominal value provides a valuable constraint on the separation of the causes, and will be considered. We will present our latest Satellite Laser Ranging and DORIS Doppler derived time series for J2, and various other low-degree harmonic terms, as well as our investigations into the causes. In addition, we will show the comparison of the J2 results with those derived from CHAMP, as computed at NASA GSFC, and the recently released GRACE gravity model.

Description: The oblateness of the Earth's gravity field, J2, has long been observed to undergo a slight decrease due to post-glacial rebound of the mantle. Sometime around 1998 this trend reversed quite suddenly. This reversal persisted until 2001, at which point the atmosphere-corrected time series appears to have reversed yet again. Presently, the time series appears to be returning to the value that would nominally have been reached had the anomaly not occurred. This anomaly signifies a large interannual change in global mass distribution whose J2 effect overshadows that of the post-glacial rebound over such timescales. A number of possible causes have been considered, with oceanic mass redistribution as the leading candidate although other effects, such as glacial melting and core effects may be contributing.