Fourier contamination of a finite length of data in tidal analysis

Summary. The contamination effect when a discrete Fourier analysis is applied to a short length of data in order to estimate the main diurnal (O₁) and semi-diurnal (M₂) components of the solid body tide is estimated and it is shown that a moderate length of record (30 days) has distinct advantages over larger record lengths of less than 60 days or so.     

On the nature of abnormal quiet days in Sq(H)

Summary. It is well known that the time of occurrence of the minimum in the horizontal component of the Earth's magnetic field ( H _min) on quiet days at a mid-latitude station on the poleward side of the Sq focus shows considerable variability from day to day. This variability has previously been discussed in terms of the incidence of so-called 'abnormal quiet days'(AQD), arbitrarily defined for the station Hartland as quiet days when H _min occurred outside an interval of ± 2½ hr centred on the most common time of 1130 LT. AQDs have some interesting properties, which have been documented, but their precise nature and cause have not been elucidated. In this paper we report the results of a study of AQDs as identified at Hartland using a chain of Northern hemisphere stations situated approximately along the 0° longitude meridian and extending on both sides of the Sq focus. It is found that there are two effects on AQDs: (i) a northward component field varying in LT is superposed at all latitudes throughout the day, so reducing the amplitude of the normal Sq(H) variation at stations poleward of the focus and increasing it on the equatorward side, (ii) a southward perturbation field, of most probable magnitude 6.0 nT for the period studied, is imposed for about 4 hr at a fixed UT at all latitudes, so constituting an 'AQD event' which can lead to the occurrence of H _min at an anomalous local time for a station poleward of the focus. It is shown that the AQD event may be of large geographical extent and that it is related to the interplanetary magnetic field. All the main properties of AQD occurrences are explained, and it is suggested that much of the day-to-day variability in the amplitude and phase of the normal Sq(H) variation probably also arises from the occurrence of AQD events at times close to the diurnal turning points.     

The change in Sq(H) amplitude on Abnormal Quiet Days

Summary. Analysis of geomagnetic data has shown that the superposed northward magnetic field, which reduces the S _q( H ) amplitude at northern mid-latitude stations on Abnormal Quiet Days, and increases the amplitude at stations on the equatorward side of the S _q focus, builds up in amplitude over four to five days before the AQD occurs, and subsides over a similar period after the AQD. It is inferred indirectly that the azimuthal component B_y of the interplanetary magnetic field varies similarly. Data for the opposite meridian show that the imposed field reverses to a southward direction at lower latitudes. The inferred currents to account for these fields are believed to flow in the ionosphere, but to arise from magnetospheric electric fields induced by the solar wind-transported IMF.     

Three-dimensional VP and VP /VS models of the upper crust in the Friuli area (northeastern Italy)

3-D images of P velocity and P - to S -velocity ratio have been produced for the upper crust of the Friuli area (northeastern Italy) using local earthquake tomography. The data consist of 2565 P and 930 S arrival times of high quality. The best-fitting V _P and V _P / V _S 1-D models were computed before the 3-D inversion. V _P was measured on two rock samples representative of the investigated upper layers of the Friuli crust. The tomographic V _P model was used for modelling the gravity anomalies, by converting the velocity values into densities along three vertical cross-sections. The computed gravity anomalies were optimized with respect to the observed gravity anomalies. The crust investigated is characterized by sharp lateral and deep V _P and V _P / V _S anomalies that are associated with the complex geological structure. High V _P / V _S values are associated with highly fractured zones related to the main faulting pattern. The relocated seismicity is generally associated with sharp variations in the V _P / V _S anomalies. The V _P images show a high-velocity body below 6 km depth in the central part of the Friuli area, marked also by strong V _P / V _S heterogeneities, and this is interpreted as a tectonic wedge. Comparison with the distribution of earthquakes supports the hypothesis that the tectonic wedge controls most of the seismicity and can be considered to be the main seismogenic zone in the Friuli area.     

On the linear relation between mb and Ms for discrimination between explosions and earthquakes

Summary. The statistical capability of the m _b: M _s discriminant for the discrimination of earthquake and explosion populations is examined by application of discriminant functions to a group of 83 explosions and 72 earthquakes in Eurasia. Equations are derived for the probability that an event is an earthquake or an explosion. The positive sign of DIS in the decision index equation, DIS _i = 34.3383 – 11.9569 mb _t + 7.1161 M _si , indicates that the event i is an earthquake. Its negative sign indicates that event i is an explosion. The probability of correct classification for an event, P _i , is related to its DIS _i value, by P _i = [1-exp (DIS _i )]⁻¹, where a large, positive DIS indicates a high probability that an event is an earthquake and a large, negative DIS indicates a high probability that an event is an explosion. The discrimination line M _s = 1.680 m _b– 4.825, or m _b= 0.595 M _s+ 2.872 very successfully separates the explosion population from the earthquake population. The points on this line have an equal chance of being an earthquake or an explosion; moreover, for any event, the distance parallel to the M _s-axis from the point representing that event in the m _b: M _s plane to this line is a measure of the probability for the correct classification of that event.     

Period and Qw of The Chandler Wobble

A splined ILS/IPMS data set (1900–1973) from the most homogeneous values available has been analysed by the maximum entropy method of Burg. Principal conclusions are: (1) the spectral character of the Chandler wobble is a single broad peak, (2) the period is 432·95 ± 1·02 mean solar days and, (3) the Q _w is 36±10. Measurements indicate that Q _w is non- stationary in time.     

Seismic image of a CO2 reservoir beneath a seismically active volcano

Mammoth Mountain is a seismically active volcano 200 000 to 50 000 years old, situated on the southwestern rim of Long Valley caldera, California. Since 1989 it has shown evidence of unrest in the form of earthquake swarms (Hill et al. 1990), volcanic 'long-period' earthquakes (Pitt & Hill 1994), increased output of magmatic ³He (Sorey et al. 1993) and the emission of about 500 tonnes day ⁻¹ of CO₂ (Farrar et al. 1995; Hill 1996; M. Sorey, personal communication, 1997), which has killed trees and poses a threat to human safety. Local-earthquake tomography shows that in mid-1989 areas of subsequent tree-kill were underlain by extensive regions where the ratio of the compressional and shear elastic-wave speeds V_P/V_S was about 9 per cent lower than in the surrounding rocks. Theory (Mavko & Mukerji 1995), experiment (Ito, DeVilbiss & Nur 1979), and experience at other geothermal/volcanic areas (Julian et al. 1996) and at petroleum reservoirs (Harris et al. 1996) indicate that V_P/V_S is sensitive to pore-fluid compressibility, through its effect on V_P . The observed V_P/V_S anomaly is probably caused directly by CO₂, and seismic V_P/V_S tomography is thus a promising tool for monitoring gas concentration and movement in volcanoes, which may in turn be related to volcanic activity.     

Regional variations of heat flow differences with depth in Alberta, Canada

Summary. Differences between estimated average heat flow values for the Mesozoic and Cenozoic formations ( Q ₁) and estimated average heat flow values for the Palaeozoic formations below the erosional unconformity ( Q ₂) are calculated for the Alberta part of the western Canadian sedimentary basin. Significant heat flow differences exist for these two intervals and the map of Δ Q = Q ₁– Q ₂ shows that Q ₂ is generally greater than Q ₁ in the western and south-western part of Alberta, while in the northern part of the province Q ₂ is generally less than Q ₁. The regional variations of Δ Q are large, with standard deviation of 26 mW m⁻² and average value –13.5 mW m⁻². A regional trend of Δ Q correlates with topographic relief and the hydraulic head variations in the basin. It is shown that there is a heat flow increase with depth in water recharge areas and a decrease in heat flow with depth in the low topographic elevation water discharge areas when comparing the average heat flow in Mesozoic + Cenozoic and Palaeozoic formations.     

Further evidence for oceanic excitation of polar motion

While the role of the atmosphere in driving variations in polar motion is well established, the importance of the oceans has been recognized only recently. Further evidence for the role of the oceans in the excitation of polar motion is presented. To estimate the equatorial excitation functions, χ ₁ and χ ₂ , for the ocean, we use velocity and mass fields from a constant-density ocean model, driven by observed surface wind stresses and atmospheric pressure, for the period 1993–1995; comparison with similar functions derived from a more complex density-stratified ocean model indicates the effectiveness of the simple constant-density modelling approach. Corresponding atmospheric excitation functions are computed from NCEP/NCAR re-analyses. Results indicate significant improvements in the agreement with the observed polar motion excitation when the simulated oceanic effects are added to atmospheric excitation. Correlations between the polar motion and the geophysical signals at periods of 15–150 days increase from 0.53 to 0.80 and from 0.75 to 0.88 for χ ₁ and χ ₂ , respectively. The oceanic signals are particularly important for seasonal variations in χ ₁ (correlation increases from 0.28 to 0.85 when oceanic excitation is included). A positive impact of the oceans on more rapid polar motion is also observed, up to periods as short as 5 days. The sensitivity of the results to different forcing fields and different amounts of friction in the oceans is also discussed.     

On the connection between the critical value of the "vertical" electrical anisotropy and earthquake prediction

Summary. It is shown that at any orientation of a dipole of direct current at the surface of a uniform anisotropic halfspace there exists a unique critical value η_cr of the "vertical" electric anisotropy η. For η > η_cr not all field lines are closed, and three straight field lines pass through the centre of the dipole. Also there exists a unique value η_eq > η_cr, where equatorial η is equal to zero. If η > η_eq then equatorial ρ_app is negative since the electric vector is turned more than through 90° at the equator of dipole, and it is shown that η_cr= 1.79, η_eq= 1.93 if the axis of dipole is directed symmetrically with respect to the axes of anisotropy. Routine observations of variations of ρ_app with the aim of earthquake prediction should be carried out taking into account the possibility of great fluctuations on η far from the epicentre of an earthquake.     

Stacking for the frequencies and Qs of 0S0 and 1S0

Summary. Using nine IDA records for the Indonesian earthquake of 1977 August 19, we have formed an optimal linear combination of the records and have measured the frequency and Q of ₀ S ₀ and ₁ S ₀. The frequency was measured using the moment ratio method. The attenuation was measured by the minimum width method and by the time-lapse method. The frequency and attenuation were measured simultaneously by varying them to obtain a best fit to the data. A 2000-hr stack, the sum of nine individual records, for ₀ S ₀ gave a frequency of 0.814664 mHz±4 ppm. The values for the Q of ₀ S ₀ for the three different methods of measurement were 5600,5833 and 5700, respectively. The error in the estimates of Q ^-1 is about 5 per cent for the minimum power method. For ₁ S ₀ a 300-hr stack yielded a frequency of 1.63151 mHz±30 ppm. The values of Q for this mode were 1960, 1800 and 1850, respectively, with an error in Q ^-1 of about 12 per cent for the minimum power method.     

A new derivation of the ratio Q(wobble)/Qμ (mantle)

Summary An extension of the Love-Larmor theory to a low-loss unelastic earth model, leads to the surprisingly simple approximation
   
where τ_s= 447.4 sidereal day is the static wobble period, τ_R= 306 sidereal day is the rigid-earth wobble period and τ_w= 433 sidereal day is the observed Chandler period. Q _W, Q _μ are the respective average Q values of the wobble and the Earth's mantle at τ_W. The known numerical factor F is only slightly dependent on the Earth structure.     

Transient and impulse responses of a one-dimensional linearly attenuating medium—II. A parametric study

Summary. We investigate one-dimensional waves in a standard linear solid for geophysically relevant ranges of the parameters. The critical parameters are shown to be T*= t_u/Q_m where t _u is the travel time and Q_m the quality factor in the absorption band, and τ⁻¹ _m , the high-frequency cut-off of the relaxation spectrum. The visual onset time, rise time, peak time, and peak amplitude are studied as functions of T* and τ _m. For very small τ _m , this model is shown to be very similar to previously proposed attenuation models. As τ _m grows past a critical value which depends on T* , the character of the attenuated pulse changes. Seismological implications of this model may be inferred by comparing body wave travel times with a'one second'earth model derived from long-period observations and corrected for attenuation effects assuming a frequency independent Q over the seismic band. From such a comparison we speculate that there may be a gap in the relaxation spectrum of the Earth's mantle for relaxation times shorter than about one second. However, observational constraints from the attenuation of body waves suggest that such a gap might in fact occur at higher frequencies. Such a hypothesis would imply a frequency dependence of Q in the Earth's mantle for short-period body waves.     

Fluxes of CO2 and CH4 in high latitude wetlands: measuring, modelling and predicting response to climate change

This review covers selected aspects of recent international efforts to measure and model greenhouse gas emission from northern wetlands, to identify the environmental factors that control gas emission, and to investigate wetlands'responses (particularly with respect to gas emission) to global change. Both bottom-up and top-to-bottom approaches, based respectively on local observations plus inventory of gas fluxes and inverse modelling of global circulation, agree on the size of the high latitude (>60_°N) contribution to global methane, which should be about 13% or 70 Tg/year. It has been shown that winter and spring fluxes are an essential part in the annual budget of CH₄ and especially CO₂ exchange (varying from 5 to 50%). Soil micro-organisms were shown to be able to respire during winter even at-16_°C. In comparison to aerobically respiring organisms, anaerobic methanogenic bacteria were less active in frozen soil, although they are subjected to significant stimulation by soil freeze-thaw cycles. The absence of immediate coupling of methanogenesis with plant photosynthesis implies that substrates for methane formation are derived from peat decomposition rather than from root exudation.     

Lithospheric flexure due to prograding sediment loads: implications for the origin of offlap/onlap patterns in sedimentary basins

Abstract Simple elastic plate models have been used to determine the stratigraphic patterns that result from prograding sediment loads. The predicted patterns, which include coastal offlap/onlap and downlap in a basinward direction, are generally similar to observations of stratal geometry from Cenozoic sequences of the U.S. Atlantic and Gulf Coast margins. Coastal offlap is a feature of all models in which the water depth and elastic thickness of the lithosphere, T _e (which is a measure of the long-term strength of the lithosphere), are held constant, and is caused by a seaward shift in the sediment load and its compensation as progradation proceeds. The coastal offlap pattern is reduced if sediments prograde into a subsiding basin, since subsidence causes an increase in the accommodation space and loading landward of a prograding wedge. The stratal geometry that results is complex, however, and depends on the sediment supply, the amount of subsidence, and T _e. If the sediment supply to a subsiding basin proceeds in distinct 'pulses' (due, say, to different tectonic events in a source region) then it is possible to determine the relationship between stratal geometry and T _e. Coastal offlap and downlap are features of most models where the lithosphere either has a constant T _e slowly increases T_e with time, or changes T _e laterally; however, in the case where sediments prograde onto lithosphere that rapidly increases T _e with rime, the offlap can be replaced by onlap. Lithospheric flexure due to prograding sediment loads is capable of producing a wide variety of stratal geometries and may therefore be an important factor to take into account when evaluating the relative role of tectonics and eustatic sea-level changes in controlling the stratigraphic record.     

Tidal gravity and ocean tide loading in Europe

b
The results are presented from tidal gravity measurements at five sites in Europe using LaCoste and Romberg ET gravimeters. Improvements that we have made to the accuracies of these gravimeters are discussed. It is shown that the 'standard' calibration of the International Center for Earth Tides, used for worldwide tidal gravity profiles, is 1.2 per cent too high. The M₂ and O₁ observations are compared with model calculations of the Earth's body tide and ocean tide loading and it is shown that there is a very significant improvement in the agreement between observations and models compared to that obtained with previous tidal gravity measurements. For O₁, where the ocean tide loading and attraction in central Europe is only 0.4 per cent of the body tide, our measurements verify that the Dehant-Wahr anelastic body tide model gravimetric factor is accurate to 0.2 per cent. It is also shown that the effects of lateral heterogeneities in Earth structure on tidal gravity are too small to explain the large anomalies in previously published tidal gravity amplitudes. The observations clearly show the importance of conserving tidal mass in the Schwiderski ocean tide model. For sites in central Europe, the M₂ and O₁ observations and the models are in agreement at the 0.1 μgal (10⁻⁹ m s⁻²) level and tidal corrections to this accuracy can now be made to absolute gravity measurements.     

Multiple scattering and energy transfer of seismic waves – separation of scattering effect from intrinsic attenuation – I. Theoretical modelling

Summary. In order to separate the scattering effect from the intrinsic attenuation, we need a multiple scattering model for seismic wave propagation in random heterogeneous media. In this paper, we apply radiative transfer theory to seismic wave propagation and formulate in the frequency domain the energy density distribution in space for a point source. We consider the cases of isotropic scattering and strong forward scattering. Some numerical examples are shown. It is seen that the energy density–distance curves have quite different shapes depending on the values of medium seismic albedo B ₀=η_s/(η_s+η_a) where η_s is the scattering coefficient and η_a is the absorption coefficient of the medium. For a high albedo ( B > 0.5) medium, the energy–distance curve is of arch shape and the position of the peak is a function of the extinction coefficient of the medium η_e=η_s+η_a. Therefore it is possible to separate the scattering effect and the absorption based on the measured energy density distribution curves.     

The detection of single-domain greigite (Fe3S4) using rotational remanent magnetization (RRM) and the effective gyro field (Bg): mineral magnetic and palaeomagnetic applications

The intensity of rotational remanent magnetization (RRM) acquired by single-domain greigite at a rotation frequency of 5 rps was combined with measurements of anhysteretic remanent magnetization (ARM) to calculate the effective biasing field ( Bg ) that produced the RRM. Samples of single-domain greigite had Bg values between -137 and -84 μT, and a MDF_RRM of c. 80 mT. By contrast, a suite of natural and synthetic ferrimagnetic iron oxide samples, including single-domain magnetite and y Fe₃O₄ tape particles, acquired Bg values between -3 and -14 μT, and MDF_RRM ranged between 43 and 68 mT (when RRM was acquired). Multidomain magnetite did not acquire a RRM. Bg values at 5 rps were calculated from previously published data for magnetite particles of different grain sizes, which revealed a minimum Bg value of -24 μT and a MDF_RRM of 57 mT for the finest fraction (0.2-0.8 μm in diameter). In a geological example, measurements of Bg and MDF_RRM were used to detect the presence of greigite in a 4 m long Late Weichselian sediment core. Variations in inclination, declination and the intensity of the natural remanent magnetization (NRM) correlate with changes in magnetic mineralogy.     

A study of geomagnetic storms

Summary. An attempt is made to find interplanetary magnetic field and solar-wind parameters which control the development of geomagnetic storms. For this purpose, the interplanetary energy flux is estimated in terms of the Poynting flux ( E × B /4π), and its time variations are compared with the rate of energy dissipation in terms of the ring-current particle injection u _i( t ), Joule dissipation in the ionosphere u_j ( t ) and auroral particle injection u_p ( t ) for 15 major geomagnetic storms.
It is shown that the growth of geomagnetic storms, namely the time variations of the rate of the total energy dissipation, u ( t ) = u _i( t ) + u _j( t ) + u _p( t ), is closely related to the Poynting flux by the following relation:
where l ₀≅ 7 R _E and θ' is a measure of the angle between the interplanetary magnetic field vector and the magnetospheric field vector at the front of the magnetosphere in the equatorial plane. Further, it is shown that within a factor of 2 for each storm period.
A large increase of u ( t ) is associated with substorm activity. Thus, the energy flux ɛ( t ) entering the magnetosphere is dissipated through magneto-spheric substorm processes within the magnetosphere, and their accumulated effects can be understood as geomagnetic storm phenomena.     

Inefficient propagation of Sn under south-eastem China

Summary. Records of 21 earthquakes with ray paths crossing south-eastern China have been examined to establish the nature of propagation of the seismic wave S_n . S_n is seen at stations within 1500km of earthquakes, but not beyond 1500km. No single localized region of attenuation can account for all the missing or weak S_n waves. The inefficient propagation of S_n through south-eastern China is surprising in light of efficient S_n propagation in other continental platforms which, like south-eastern China, appear to have been stable throughout the Cenozoic.