Gravity potential spherical harmonic series

Time variable gravity fields, reflecting variations of mass distribution in the system Earth is one of the key parameters to understand the changing Earth. Mass variations are caused either by redistribution of mass in, on or above the Earth's surface or by geophysical processes in the Earth's interior. The first set of observations of monthly variations of the Earth gravity field was provided by the US/German GRACE satellite mission beginning in 2002. This mission is still providing valuable information to the science community. However, as GRACE has outlived its expected lifetime, the geoscience community is currently seeking successor missions in order to maintain the long time series of climate change that was begun by GRACE. Several studies on science requirements and technical feasibility have been conducted in the recent years. These studies required a realistic model of the time variable gravity field in order to perform simulation studies on sensitivity of satellites and their instrumentation. This was the primary reason for the European Space Agency (ESA) to initiate a study on ''Monitoring and Modelling individual Sources of Mass Distribution and Transport in the Earth System by Means of Satellites''. The goal of this interdisciplinary study was to create as realistic as possible simulated time variable gravity fields based on coupled geophysical models, which could be used in the simulation processes in a controlled environment. For this purpose global atmosphere, ocean, continental hydrology and ice models were used. The coupling was performed by using consistent forcing throughout the models and by including water flow between the different domains of the Earth system. In addition gravity field changes due to solid Earth processes like continuous glacial isostatic adjustment (GIA) and a sudden earthquake with co-seismic and post-seismic signals were modelled. All individual model results were combined and converted to gravity field spherical harmonic series, which is the quantity commonly used to describe the Earth's global gravity field. The result of this study is a twelve-year time-series of 6-hourly time variable gravity field spherical harmonics up to degree and order 180 corresponding to a global spatial resolution of 1 degree in latitude and longitude. In this paper, we outline the input data sets and the process of combining these data sets into a coherent model of temporal gravity field changes. The resulting time series was used in some follow-on studies and is available to anybody interested.

Observations on the propagation of internal tidal waves on the continental slope off Sierra Leone

During R/V Meteor-cruise no. 30 4 moorings with 17 current meters were placed on the continental slope of Sierra Leone at depths between 81 and 1058 meters. The observation period started on March 8, 1973, 16.55 hours GMT and lasted 19 days for moorings M30_068MOOR, M30_069MOOR, M30_070MOOR on the slope and 9 days for M30_067MOOR on the shelf. One current meter recorded at location M30_067MOOR for 22 days. Hydrographic data were collected at 32 stations by means of the "Kieler Multi-Meeressonde". Harmonic analysis is applied to the first 15 days of the time series to determine the M2 and S2 tides. By vertically averaging of the Fourier coefficients the field of motion is separated into its barotropic and its baroclinic component. The expected error generated by white Gaussian noise is estimated. To estimate the influence of the particular vertical distribution of the current meters, the barotropic M2 tide is calculated by ommitting and interchanging time series of different moorings. It is shown that only the data of moorings M30_069MOOR, M30_070MOOR and M30_067MOOR can be used. The results for the barotropic M2 tide agree well with the previous publications of other authors. On the slope at a depth of 1000 m there is a free barotropic wave under the influence of the Coriolis-force propagating along the slope with an amplitude of 3.4 cm S**-1. On the shelf, the maximum current is substantially greater (5.8 cm s**-1) and the direction of propagation is perpendicular to the slope. As for the continental slope a separation into different baroclinic modes using vertical eigenmodes is not reasonable, an interpretation of the total baroclinic wave field is tried by means of the method of characteristis. Assuming the continental slope to generate several linear waves, which superpose, baroclinic tidal ellipses are calculated. The scattering of the direction of the major axes M30_069MOOR is in contrast to M30_070MOOR, where they are bundled within an angle of 60°. This is presumably caused by the different character of the bottom topography in the vicinity of the two moorings. A detailed discussion of M30_069MOOR is renounced since the accuracy of the bathymetric chart is not sufficient to prove any relation between waves and topography. The bundeling of the major axes at M30_070MOOR can be explained by the longslope changes of the slope, which cause an energy transfer from the longslope barotropic component to the downslope baroclinic component. The maximum amplitude is found at a depth of 245 m where it is expected from the characteristics originating at the shelf edge. Because of the dominating barotropic tide high coherence is found between most of the current meters. To show the influence of the baroclinic tidal waves, the effect of the mean current is considered. There are two periods nearly opposite longshore mean current. For 128 hours during each of these periods, starting on March 11, 05.00, and March 21, 08.30, the coherences and energy spectra are calculated. The changes in the slope of the characteristics are found in agreement with the changes of energy and coherence. Because of the short periods of nearly constant mean current, some of the calculated differences of energy and coherence are not statistically significant. For the M2 tide a calculation of the ratios of vertically integrated total baroclinic energy and vertically integrated barotropic kinetic energy is carried out. Taking into account both components (along and perpendicular to the slope) the obtained values are 0.75 and 0.98 at the slope and 0.38 at the shelf. If each component is considered separately, the ratios are 0.39 and 1.16 parallel to the slope and 5.1 and 15.85 for the component perpendicular to it. Taking the energy transfer from the longslope component to the doenslope component into account, a simple model yields an energy-ratio of 2.6. Considering the limited application of the theory to the real conditions, the obtained are in agreement with the values calculated by Sandstroem.

Atmospheric, oceanic and hydrological polar motion excitation mechanisms based on a combination of geometric and gravimetric space observations

The goal of our study is to determine accurate time series of geophysical Earth rotation excitations to learn more about global dynamic processes in the Earth system. For this purpose, we developed an adjustment model which allows to combine precise observations from space geodetic observation systems, such as Satellite Laser Ranging (SLR), Global Navigation Satellite Systems (GNSS), Very Long Baseline Interferometry (VLBI), Doppler Orbit determination and Radiopositioning Integrated on Satellite (DORIS), satellite altimetry and satellite gravimetry in order to separate geophysical excitation mechanisms of Earth rotation. Three polar motion time series are applied to derive the polar motion excitation functions (integral effect). Furthermore we use five time variable gravity field solutions from Gravity Recovery and Climate Experiment (GRACE) to determine not only the integral mass effect but also the oceanic and hydrological mass effects by applying suitable filter techniques and a land-ocean mask. For comparison the integral mass effect is also derived from degree 2 potential coefficients that are estimated from SLR observations. The oceanic mass effect is also determined from sea level anomalies observed by satellite altimetry by reducing the steric sea level anomalies derived from temperature and salinity fields of the oceans. Due to the combination of all geodetic estimated excitations the weaknesses of the individual processing strategies can be reduced and the technique-specific strengths can be accounted for. The formal errors of the adjusted geodetic solutions are smaller than the RMS differences of the geophysical model solutions. The improved excitation time series can be used to improve the geophysical modeling.