New aspects of the equilibrium pole tide

Summary. We have developed a new spherical harmonic algorithm for the calculation of the loading and self-gravitating equilibrium pole tide. Based on a suggestion of Dahlen, this approach minimizes the distortions in tide height caused by an incomplete representation of the ocean function. With slight modification our approach easily could be used to compute self-gravitating and loading luni-solar tides as well.
Using our algorithm we have compared the static pole tide with tide observations at a variety of locations around the world. We find statistically significant evidence for pole tide enhancements in mid-ocean as well as the shallow seas.
We have also re-investigated the effect of the static tide on the Chandler wobble period. The difference between the wobble period of an oceanless, elastic earth with a fluid core (Smith & Dahlen) and the period of an earth minus static oceans yields a 7.4-day discrepancy. We conclude from tide observations that much of the discrepancy can probably be accounted for by non-equilibrium pole tide behaviour in the deep oceans.     

Modelling the pole tide and its effect on the Earth's rotation

Summary. The pole tide is the response of the ocean to incremental centrifugal forces associated with the Chandler wobble. The tide has a potentially important effect on the period and damping of the wobble, but it is at present not well constrained by observations. Here, we construct both analytical and numerical models for the pole tide. The analytical models consider the tide first in a global ocean and then in an enclosed basin on a beta-plane. The results are found to approach equilibrium linearly with decreasing frequency and inversely with increasing basin depth. The numerical models solve Laplace's tidal equations over the world's oceans using realistic continental boundaries and bottom topography. The results indicate that the effects of the non-equilibrium portion of the deep ocean tide on the Chandler wobble period and damping are negligible.     

Wind stress forcing of the North Sea 'pole tide'

We conduct numerical simulations of the wind forcing of sea level variations in the North Sea using a barotropic ocean model with realistic geography and bathymetry to examine the forcing of the 14 month 'pole tide', which is known to be anomalously large 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 ocean model output sea level response is then compared with 10 coastal tide gauge records from the Permanent Service for Mean Sea Level (PSMSL) over the same period of time. Besides the strong seasonal variations, several prominent quasi-periodicities exist near 7 years, 3 years, 14 months, 9 months and 6.5 months. Correlations 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 set-up response to wind stress forcing with a periodicity of around 14 months. We find no need to invoke a geophysical explanation involving resonance enhancement of the pole tide in the North Sea to explain the observations.     

Estimation of the pole tide gravimetric factor at the chandler period through wavelet filtering

Wavelet analysis for filtering is used to improve estimation of gravity variations induced by Chandler wobble. This method eliminate noise in superconducting gravimeter (SG) records with bandpass filters derived from Daubechies wavelet. The SG records at four European stations (Brussels, Membach, Strasbourg and Vienna) are analysed in this study. First, the earth tidal constituents are removed from the observed data by using synthetic tides, then the gravity residuals are filtered into a narrow period band of 256–512 d by a wavelet bandpass filter. These data are submitted to three regression analysis methods for estimating the gravimetric factor of the Chandler wobble. After processing by wavelet filtering, SG records can provide amplitude factors δ and phase lags κ of the Chandler wobble with much smaller mean square deviation (MSD) than these provided by former studies. It is mainly because the wavelet method can effectively eliminate instrumental drift and provide smoothed data series for the regression analysis.     

Random walk of the Earth's pole related to the Chandler wobble excitation

Summary. The equation governing the polar motion shows that the polar secular drift and the Chandler wobble amplitude are related to each other. In particular, a drift of the mean pole position comes out as a consequence of the maintenance of the Chandler wobble by possible step perturbations of the Earth's inertia tensor.
The minimum excitation functions necessary to explain the Chandler wobble amplitude variations for the period 1901–84 are derived from the Chandler term, with the hypothesis that the excitations follow a uniform random distribution in time. It is shown that they have the statistical properties of the steps of a two-dimensional random walk. These functions are then used to derive, from a statistical simulation, a lower limit of the secular drift which may result from the excitation of the Chandler wobble.
The drift generated by the random walk is of the same order of magnitude as the observed secular drift for the period 1901–84, but their time dependence is different. This indicates that the observed secular drift cannot be explained as the consequence of an excitation of the Chandler wobble by random steps of the Earth's inertia tensor. However, the possible contribution of the Chandler wobble excitation to the polar drift has to be taken into account when other mechanisms, such as lithospheric rebound related to deglaciation, are proposed.     

The influence of earthquakes on the Chandler wobble during 1977–1983

Summary. A direct calculation is made of the effect on the Chandler wobble of 1287 earthquakes that occurred during 1977–1983. The hypocentral parameters (location and origin time) and the moment tensor representation of the best point source for each earthquake as determined by the 'centroidmoment tensor' technique were used to calculate the change in the Chandler wobble's excitation function by assuming this change is due solely to the static deformation field generated by that earthquake. The resulting theoretical earthquake excitation function is compared with the 'observed' excitation function that is obtained by deconvolving a Chandler wobble time series derived from LAGEOS polar motion data. Since only 7 years of data are available for analysis it is not possible to resolve the Chandler band and determine whether or not the theoretical earthquake excitation function derived here is coherent and in phase with the 'observed' excitation function in that band. However, since the power spectrum of the earthquake excitation function is about 56 dB less than that of the 'observed' excitation function at frequencies near the Chandler frequency, it is concluded that earthquakes, via their static deformation field, have had a negligible influence on the Chandler wobble during 1977–1983. However, fault creep or any type of aseismic slip that occurs on a time-scale much less than the period of the Chandler wobble could have an important (and still unmodelled) effect on the Chandler wobble.     

Rift kinematics preserved in deep-time erosional landscape below the northern North Sea

Our understanding of continental rifting is, in large parts, derived from the stratigraphic record. This record is, however, incomplete as it does not often capture the geomorphic and erosional signal of rifting. New 3D seismic reflection data reveal a Late Permian-Early Triassic landscape incised into the pre-rift basement of the northern North Sea. This landscape, which covers at least 542 km², preserves a drainage system bound by two major tectonic faults. A quantitative geomorphic analysis of the drainage system reveals 68 catchments, with channel steepness and knickpoint analysis of catchment-hosted palaeo-rivers showing that the landscape preserved a >2 Myr long period of transient tectonics. We interpret that this landscape records a punctuated uplift of the footwall of a major rift-related normal fault (Vette Fault) at the onset of rifting. The landscape was preserved by a combination of relatively rapid subsidence in the hangingwall of a younger fault (Øygarden Fault) and burial by post-incision sediments. As such, we show how and why erosional landscapes are preserved in the stratigraphic record, and how they can help us understand the tectono-stratigraphic evolution of ancient continental rifts.     

Estimation of period and Q of the Chandler wobble

The period P and Q -value of the Chandler wobble are two fundamental functional of the Earth's internal physical properties and global geodynamics. We revisit the problem of the estimation of P and Q , using 10.8 yr of modern polar motion as well as contemporary atmospheric angular momentum (AAM) data. We make full use of the knowledge that AAM is a major broad-band excitation source for the polar motion. We devise two optimization criteria under the assumption that, after removal of coherent seasonal and long-period signals, the non-AAM excitation is uncorrelated with the AAM. The procedures lead to optimal estimates for P and Q. Our best estimates, judging from comprehensive sets of Monte Carlo simulations, are P = 433.7 ± 1.8 (1σ) days, Q =49 with a la range of (35, 100). In the process we also obtain (as a by-product) an estimate of roughly 0.8 for a 'mixing ratio' of the inverted-barometer (IB) effect in the AAM pressure term, indicating that the ocean behaves nearly as IB in polar motion excitation on temporal scales from months to years     

Joint traveltime inversion of wide-angle seismic data and a deep reflection profile from the central North Sea

We report results from the Seismic Wide-Angle and Broadband Survey carried out over the Mid North Sea High. This paper focuses on integrating the information from a conventional deep multichannel reflection profile and a coincident wide-angle profile obtained by recording the same shots on a set of ocean bottom hydrophones (OBH). To achieve this integration, a new traveltime inversion scheme was developed (reported elsewhere) that was used to invert traveltime information from both the wide-angle OBH records and the reflection profile simultaneously. Results from the inversion were evaluated by producing synthetic seismograms from the final inversion model and comparing them with the observed wide-angle data, and an excellent match was obtained. It was possible to fine-tune velocities in less well-resolved parts of the model by considering the critical distance for the Moho reflection. The seismic velocity model was checked for compatibility with the gravity field, and used to migrate and depth-convert the reflection profile. The unreflective upper crust is characterized by a high velocity gradient, whilst the highly reflective lower crust is associated with a low velocity gradient. At the base of the crust there are several subhorizontal reflectors, a few kilometres apart in depth, and correlatable laterally for several tens of kilometres. These reflectors are interpreted as representing a strike section through northward-dipping reflectors at the base of the crust, identified on orthogonal profiles by Freeman et al. (1988) as being slivers of subducted and imbricated oceanic crust, relics of the mid-Palaeozoic Iapetus Ocean.