第27次南极科考度夏期间长城站地震观测

在第27次中国南极科学考察度夏期间,针对长城站原地震台受干扰源（如发电栋、冷库）的影响,对台站进行了重新选址工作,并完成了新台站的基建工程（摆坑、电缆布设等）,架设了一套宽频带地震仪。对新台站记录的数据进行分析表明新台站运行稳定,背景噪音小,记录的地震事件信噪比高,这为下一步开展相关地震学方面的研究将提供高质量的数据。尽管本次科考地震观测时间短,但是新台站不仅观测到了远震记录,而且观测到了长城站附近地区的近震和疑似冰震的记录。单台定位结果显示4次近震分布基本与南设得兰群岛走向平行,且震源深度由东北向西南逐渐变深。     

On reduction of long-period horizontal seismic noise using local barometric pressure

Tilt from atmospheric loading has long been known to be the major source of long-period horizontal seismic noise. We try to quantify these effects for seismic data from the Black Forest Observatory (BFO), which is known to be a very quiet station. Experimental transfer functions between local barometric pressure and horizontal seismic noise are estimated for two long time-series by standard methods. Two simple analytical physical models are developed: the local deformation model (LDM) and the acoustic-gravity wave model (TWM). Subsequently these models, with only two free parameters are fit using least squares to the observed seismic noise for time-series of widely differing lengths. The results are variable, sometimes rather dramatic variance reductions are obtained and sometimes the reduction is hardly significant. The method produces the best results when barometrically induced noise is high. The resulting admittances for the LDM are compared to finite element calculations. Since the methods are simple and can result in conspicuous reductions in noise we provide one more reason for installing barometers at even the best broad-band seismic stations.     

Determination of Lg-wave attenuation using single-station seismograms: a case study for Zimbabwe

In this paper we show how the quality factor Q may be calculated using a single seismograph station and a number of events recorded on analogue seismograms. We followed Nuttli's (1973) method and extended it to one seismograph station. Using the single station Bulawayo (BUL), we determined a mean Q value of 650 for Zimbabwe. Furthermore, we considered different propagation paths over Zimbabwe as the seismic waves travelled to BUL and found a low Q value of 350 for the Deka fault zone. the Q value of 650 obtained in this study agrees well within error with that of 603 reported by Chow et al. (1980) using the multistation, multi-event method, and lies within the range of Q values (454–759) reported by Xie & Mitchell (1990) using a 'back-project'method to image large-scale lateral variations of Lg coda Q . the results obtained are important as the Q value constitutes part of the input data in seismic hazard calculations. the method may be used to determine Q in regions where there may be only one station with reliable analogue seismogram data.     

Effects of the free surface on shear wavetrains

Summary. The behaviour of shear-waves is of great importance in identifying and investigating seismic anisotropy in the Earth. However, shear wavetrains recorded at the Earth's surface do not always reflect the motion at depth, introducing practical problems of interpretation. Shear wavetrains incident on the surface of an isotropic half-space at angles less than critical (about 35°) are broadly preserved, but at greater angles substantial distortions can occur. For stations situated close to the source, as in local earthquake studies, the local SP phase, a radially polarized precursor to S , may occur. The behaviour at the surface of an anisotropic half-space is further complicated by the divergence of phase and energy propagation vectors. All of these complications suggest that detailed seismogram modelling is essential to any study of shear wave propagation in the Earth, and in particular to investigations of anisotropy-induced shear-wave splitting.     

一次夏季强沙尘暴形成机理的综合分析

利用沙尘暴发生前夕酒泉基准站的探空资料、沙尘暴发生前后地面自动站观测资料和有关台站的常规观测资料对2004年7月12日出现在我国内蒙古西部和河西走廊西端的一次沙尘暴天气过程,从大气层结稳定度、地形因素对沙尘暴影响的机理以及下垫面性质等方面进行了分析,并对天气现象发生时地面气象要素和地温的变化进行了研究。结果表明：①大气层结的不稳定起了激发作用;②该地区特殊的地形条件所造成的“狭管效应”起到了加强作用;③沙尘暴爆发前风速迅速增大,为起沙提供了动力条件;④气压和相对湿度在沙尘暴发生前夕和沙尘暴过程的中后期有明显扰动, 沙尘暴结束后温度、气压、相对湿度和地面等气象要素有显著性突变。     

Traveltime measurements from noise correlation: stability and detection of instrumental time-shifts

We test the feasibility of using Green's functions extracted from records of ambient seismic noise to monitor temporal changes in the Earth crust properties by repeated measurements at regional distances. We use about 11 yr of continuous recordings to extract surface waves between three pairs of stations in California. The correlations are computed in a moving 1-month window and we analyse the temporal evolution of measured interstation traveltimes. The comparison of the arrival times in the positive and negative correlation time of Rayleigh and Love waves allows us to separate time-shifts associated with any form of physical change in the medium, those resulting from clock drift or other instrumental errors, and those due to change in the localization of the noise sources. This separation is based on the principle of time symmetry. When possible, we perform our analysis in two different period bands: 5–10 and 10–20 s. The results indicate that significant instrumental time errors (0.5 s) are present in the data. These time-shifts can be measured and tested by closure relation and finally corrected independently of any velocity model. The traveltime series show a periodic oscillation that we interpret as the signature of the seasonal variation of the region of origin of the seismic noise. Between 1999 and 2005, the final arrival time fluctuations have a variance of the order of 0.01 s. This allows us to measure interstation traveltimes with errors smaller than 0.3 per cent of the interstation traveltime and smaller than 1 per cent of the used wave period. This level of accuracy was not sufficient to detect clear physical variation of crustal velocity during the considered 11 yr between the three stations in California. Such changes may be more easily detectable when considering pairs of stations more closely located to each other and in the vicinity of tectonically active faults or volcanoes.     

A Preliminary Investigation into the Relationship between Long-Period Seismic Noise and Local Fluctuations in the Atmospheric Pressure Field

Summary The displacement response of an elastic half space to a plane pressure wave is examined in order to establish the conditions under which sources of this type can contribute significantly to the long-period seismic noise field. The study is restricted to pressure waves which propagate at velocities well below the seismic wave velocities characteristic of the half space. The numerical studies indicate that pressure waves with amplitudes of 100 μbar or more can contribute significantly to the long-period vertical background noise observed at the surface, provided that the detectors are located on sections of alluvial fill or poorly to moderately indurated sandstones and shales whose thicknesses are greater than about a kilometre. These same waves can also create significant tilt noise on long-period horizontal seismographs located at or near the surface, regardless of the rock type. The seismic disturbances created by pressure waves decay rapidly away from the surface. Therefore, it appears that it may be possible to eliminate the effects of atmospherically generated noise by placing the detectors at moderate depths.     

Seismo-acoustic analysis of the Buncefield oil depot explosion in the UK, 2005 December 11

A massive vapour cloud explosion occurred at the Buncefield fuel depot near Hemel Hempstead, UK, in the morning of 2005 December 11. The explosion was the result of an overflow from one of the storage tanks with the release of over 300 tons of petrol and generating a vapour cloud that spread over an area of 80 000 m², before being ignited. Considerable damage was caused in the vicinity of the explosion and a total of 43 people were injured. The explosion was detected by seismograph stations in the UK and the Netherlands and by infrasound arrays in the Netherlands. We analysed the seismic recordings to determine the origin time of 06:01:31.45 ±0.5 s (UTC) from P -wave arrival times. Uncertainties in determination of origin time from acoustic arrival times alone were less than 10 s. Amplitudes of P -, Lg and primary acoustic waves were measured to derive decay relationships as function of distance. From the seismic amplitudes we estimated a yield of 2–10 tons equivalent to a buried explosion. Most seismic stations recorded primary and secondary acoustic waves. We used atmospheric ray tracing to identify the various travel paths, which depend on temperature and wind speed as function of altitude, leading to directional variation. Refracted waves were observed from the troposphere, stratosphere and thermosphere with a good match between observed and calculated traveltimes. The various wave types were also identified through array processing, which provides backazimuth and slowness, of recordings from an infrasound array in the Netherlands. The amplitude of stratospheric refracted acoustic waves recorded by the array microbarometers was used to estimate a yield of 21.6 (±5) tons TNT equivalent. We have demonstrated through joint seismo-acoustic analysis of the explosion that both the seismic velocity model and the atmospheric model are sufficient to explain the observed traveltimes.     

Coseismic and other short-term strain changes recorded with Sacks-Evertson strainmeters in a deep mine, South Africa

Summary. During 1977 March and April, three Sacks-Evertson borehole dilatometers were installed at the ends of boreholes drilled into the sidewall of an experimental tunnel at a depth of 3.1 km in the ERPM gold mine near Johannesburg. In the following year coseismic strain changes ranging from 5 ± 10⁻¹⁰ to values exceeding 5 ± 10⁻⁶ were recorded for hundreds of mine tremors in the magnitude range -1 to 3.7 and at hypocentral distances of 50 m to about 2 km. Hypocentral coordinates and magnitudes were determined from seismograms recorded from an underground array of geophones. Amplitudes and polarities of the coseismic strain steps are generally in excellent agreement with theoretical expectations based on point-source dislocation theory; specifically, the strain steps are proportional to the seismic moment divided by the cube of hypocentral distance. At a strain level of 5 ± 10⁻⁹ or greater the tremors do not appear to be preceded by any short-term indications of instability even for tremors producing coseismic steps greater than 5 ± 10⁻⁶ and for which the strainmeters were within a source radius of the hypocentre. Continuous strain changes observed at the times when the mine excavation, at a distance of about 100 m, is extended are in good agreement with calculated changes based on the theory of elasticity. A similar calculation is consistent with post-seismic strain changes observed to follow some of the closer tremors. These post-seismic strains show a logarithmic dependence on time following the tremor and appear to be due to the interaction of a tremor with the adjacent mine excavation rather than to deformation within the actual seismic source region.     

Could the energy near the FCN and the FICN be explained by luni-solar or atmospheric forcing?

The observed time-series of precession/nutation show residuals with respect to an empirical model based on the rigid Earth theoretical nutations and a frequency dependent transfer function with resonances to the Earth's normal modes. These residuals display energy mainly in the frequency domain around 430 and 500 days in the inertial frame. In this frequency band, the energy is possibly related to two normalmode frequencies: the free core nutation (FCN) and the free inner core nutation (FICN). In this paper, we examine the possibility of obtaining this energy from the resonance effect induced by a luni-solar (or planetary) forcing, or by an atmospheric forcing at a frequency very close to these Earth free nutations. The amplification factor due to the resonance is computed from an analytical formula expressed in the case of a simplified three-layer ellipsoidal rotating earth (with an elastic inner core, a liquid outer core and an elastic mantle), as well as the empirical formula based on the analysis of VLBI observations. For the tidal forcing, the theoretical results do not show any resonance at the level of precision we have examined but it is still possible to find one frequency near the FCN or FICN frequencies which could be excited. In contrast, for the atmospheric pressure the level of energy needed could be obtained from the diurnal pressure, depending on the noise level of the Earth's global pressure. We also show that the combination of three waves can explain the observed decrease of energy with time. While the tidal potential amplitudes are too small, a pressure noise level of 0.5 Pa would be sufficient to excite these waves.     

The Structure of Atmospheric Turbulence and its Application to the Design of Pipe Arrays

Summary Turbulent boundary layers at the surface of the Earth limit the detection of infrasonic waves with periods greater than 1 s. Pipe arrays designed to improve the signal-to-noise ratios of infrasonic waves usually assume that the background noise due to this turbulent boundary layer is incoherent between the array inlets. The power at various points on a surface was measured; coherences between these points were determined and they were found to be significant in the period range 1–100 s. Such coherent noise must be considered when pipe arrays are designed.     

Broad-band calibration of marine seismic sources used by R/V Polarstern for academic research in polar regions

Air guns and air-gun arrays of different volumes are used for scientific seismic surveys with R/V Polarstern in polar regions. To assess the potential risk of these research activities on marine mammal populations, knowledge of the sound pressure field of the seismic sources is essential. Therefore, a broad-band (0–80 kHz) calibration study was conducted at the Heggernes Acoustic Range, Norway. A GI (2.4 l), a G (8.5 l) and a Bolt gun (32.8 l) were deployed as single sources, 3 GI (7.4 l), 3 G (25.6 l) and 8 VLF™ Prakla-Seismos air guns (24.0 l) as arrays. Each configuration was fired along a line of 3–4 km length running between two hydrophone chains with receivers in 35, 100, 198 and 263 m depth. Peak-to-peak, zero-to-peak, rms and sound exposure levels (SEL) were analysed as functions of range. They show the typical dipole-like directivity of marine seismic sources with amplitude cancellation close to the sea surface, higher amplitudes in greater depths, and sound pressure levels which continuously decrease with range. Levels recorded during the approach are lower than during the departure indicating a shadowing effect of Polarsterns's hull. Backcalculated zero-to-peak source levels range from 224–240 dB re 1 μPa @ 1 m. Spectral source levels are highest below 100 Hz and amount to 182–194 dB re 1 μPa Hz^–1. They drop off continuously with range and frequency. At 1 kHz they are ∼30 dB, at 80 kHz ∼60 dB lower than the peak level. Above 1 kHz amplitude spectra are dominated by Polarstern's self-noise. From the rms and sound exposure levels of the deepest hydrophone radii for different thresholds are derived. For a 180 dB rms-level threshold radii maximally vary between 200 and 600 m, for a 186 dB SEL threshold between 50 and 300 m.     

Air-Coupled Seismic Waves at Long Range from Apollo Launchings

Summary Microphones and seismographs were co-located in arrays on Skidaway Island, Georgia, for the launchings of Apollo 13 and 14, 374 km to the south. Simultaneous acoustic and seismic waves were recorded for both events at times appropriate to the arrival of the acoustic waves from the source. Significant comparisons of the true signals are (1) the acoustic signal is relatively broadband compared to the nearly monochromatic seismic signal; (2) the seismic signal is much more continuous than the more pulse-like acoustic signal; (3) ground loading from the pressure variations of the acoustic waves is shown to be too small to account for the seismic waves; (4) the measured phase velocities of both acoustic and seismic waves across the local instrument arrays differ by less than 6 per cent and possibly 3 per cent if experimental error is included. It is concluded that the seismic waves are generated by resonant coupling to the acoustic waves along some 10 km of path on Skidaway Island. The thickness of unconsolidated sediment on the island is appropriate to a resonant ground wave frequency of 3.5 to 4 Hz, as observed. Under appropriate conditions, ground wave observations may prove more effective means of detecting certain aspects of acoustic signals in view of the filtering of wind noise and amplification through resonance.     

Precise measurement of seismic travel times from an air gun source

Summary. Seismic travel times for extrema, zero-crossings, or entire body waves need to be determined precisely to one part in 10³ or better in several varieties of seismic studies employing an impulsive artificial source. Examples are crosshole surveys which delineate rock crack distribution separating the holes and monitoring of crustal seismic travel times in earthquake precursor studies. A timing resolution of one part in 10³ has been achieved previously using digitally recorded seismic data. These methods, however, do not use interpolation between digitized data points as a method to increase the timing resolution. We report travel-time determinations based on interpolation between digitized points which achieve a precision of two parts in 10⁴, a five-fold improvement over the existing methods. In addition, the effects of seismic noise on travel-time measurement have been compared for the extremum location, the unnormalized correlation, and the normalized correlation method. The following conclusions are drawn from this comparison: (1) the normalized correlation method provides an 18–55 per cent improvement in the standard deviation of the mean over the extremum location method, and (2) results as accurate as those by the normalized correlation can be obtained by the unnormalized correlation if a complete up-and-down swing of the waveform is used as the master trace and if the master trace is close to being sinusoidal. The advantage of the unnormalized correlation over the normalized correlation is speed; the unnormalized correlation is faster by a factor of 28 in computing time.     

Some remarks about the degree-one deformation of the Earth

The degree-one deformation of the Earth (and the induced discrepancy between the figure centre and the mass centre of the Earth) is computed using a theoretical approach (Love numbers formalism) at short timescales (where the Earth has an elastic behaviour) as well as at long timescales (where the Earth has a viscoelastic or quasi-fluid behaviour). For a Maxwell model of rheology, the degree-one relaxation modes associated with the viscoelastic Love numbers have been investigated: the M_o mode does not exist and there is only one transition mode (instead of two) generated by a viscosity discontinuity.
The translations at each interface of the incompressible layers of the earth model [surface, 670 km depth discontinuity, core-mantle boundary (CMB) and inner-core boundary (ICB)] are computed. They are elastic with an order of magnitude of about 1 mm when the excitation source is the atmospheric continental loading or a magnetic pressure acting at the CMB. They are viscoelastic when the earth is submitted to Pleistocene deglaciation, with an order of magnitude of about 1 m. In a quasi-fluid approximation (Newtonian fluid) because of the mantle density heterogeneity their order of magnitude is about 100 m (except for the ICB, which is in quasi-hydrostatic equilibrium at this timescale).     

Magnitude and Q determinations in southern Africa using Lg wave amplitudes

b
The amplitude of vertical, short period (1 s) Lg -waves from 575 shallow earthquakes recorded within the distance range 0|Mo-30|Mo by the Rhodesian seismograph network during the period 1968–77 are analysed to separate the effects of earthquake size, epicentral distance and station structure.
When corrected for geometrical spreading and Airy phase dispersion the decay of amplitude with distance yields an estimate of anelastic attenuation of 0.160 deg^-1 which gives an average value of Q (the specific quality factor) of 603 |Mp 50 for propagation paths that lie along and across the East African Rift System. Inversion of the amplitude—distance curve gives the calibration or distance normalizing function. Thereby the amplitude of Lg can be used to provide an estimate of the size of small, local earthquakes in terms of the teleseismic body wave magnitude m_b (after Henderson). The station effects of the six seismograph stations making up the network all lie within |Mp0.1 magnitude units. Since three of the stations lie on the Rhodesian craton while the remaining three lie on Precambrian mobile belts adjacent to the craton, the Precambrian basement geology does not significantly affect the amplitude of Lg     

Preliminary Results on the 3-dimensional Seismic Structure of the Lithosphere under the USGS Central California Seismic Array

About 1500 readings of teleseismic P -time residuals obtained from the US Geological Survey seismograph network in central California have been used to obtain a three-dimensional image of seismic velocity anomalies for this area by the method of Aki, Christoffersson & Husebye We found that the California network is less suitable than the LASA and NORSAR arrays for this kind of studies because of its greater proportion of peripheral blocks in which the resolution is very poor for the stochastic inverse solution and the random error effect is severe for the generalized inverse solution. Nevertheless, the resultant velocity anomalies show a remarkable correlation with the San Andreas fault zone to a depth of 75 km. The anomaly pattern changes drastically as the depth exceeds 75 km, suggesting that the asthenosphere has been reached.     

Double seismic zones and stresses of intermediate depth earthquakes

Summary. Data from Japanese local seismograph networks suggest that the stresses in double seismic zones are in-plate compression for the upper zone and in-plate tension for the lower zone; the stresses do not necessarily appear to be down-dip. It may therefore be possible to identify other double seismic zones on the basis of data which indicate that events with differing orientations of in-plate stresses occur in a given segment of slab.
A global survey of published focal mechanisms for intermediate depth earthquakes suggests that the stress in the slab is controlled, at least in part, by the age of the slab and the rate of convergence. Old and slow slabs are under in-plate tensile stresses and the amount of in-plate compression in the slab increases with increasing convergence rate or decreasing slab age. Young and fast slabs are an exception to this trend; all such slabs are down-dip tensile. Since these slabs all subduct under continents, they may be bent by continental loading. Double seismic zones are not a feature common to all subduction zones and are only observed in slabs which are not dominated by tensile or compressive stresses.
Unbending of the lithosphere and upper mantle phase changes are unlikely to be the causes of the major features of double zones, although they may contribute to producing some of their characteristics. Sagging or thermal effects, possibly aided by asthenospheric relative motion, may produce the local deviatoric stresses that cause double zones.     

An earthquake sequence studied with ocean-bottom seismographs

Summary. A tripartite ocean-bottom seismograph array at the junction of the East Pacific Rise and Rivera Fracture Zone recorded an eathquake sequence, consisting of three main shocks ( m _B= 4.3, 4.3 and 4.8) and numerous aftershocks from the fracture zone, in the distance range 35–50 km. Delineation of the rupture zones by aftershocks indicates that the first two main shocks took place on overlapping fault areas, while the third occurred over a fault area separated from the first by several kilometres. Both rupture zones were about 4 km long. Surface wave spectra indicate a shallow (about 3 km below the sea floor) source, as does OBS array phase velocity data. The seismic moments, obtained from teleseismic surface wave data, of 1.3, 2.1 and 2.8 × 10²³ dyn cm, with the fault areas as delineated by aftershocks, imply a stress drop of about 8 bars for the main shocks. Aftershock sequences of each of the main shocks are similar, with a b -value of about 0.65. Teleseismic P travel times are similar to those from near-surface sources in Nevada.     

Inferring the relative measure of principal stress components

The phase velocity and the attenuation coefficient of compressional seismic waves, propagating in poroelastic, fluid-saturated, laminated sediments, are computed analytically from first principles. The wavefield is found to be strongly affected by the medium heterogeneity. Impedance fluctuations lead to poroelastic scattering; variations of the layer compressibilities cause inter-layer flow (a 1-D macroscopic local flow). These effects result in significant attenuation and dispersion of the seismic wavefield, even in the surface seismic frequency range, 10–100 Hz. The various attenuation mechanisms are found to be approximately additive, dominated by inter-layer flow at very low frequencies. Elastic scattering is important over a broad frequency range from seismic to sonic frequencies. Biot's global flow (the relative displacement of solid frame and fluid) contributes mainly in the range of ultrasonic frequencies. From the seismic frequency range up to ultrasonic frequencies, attenuation due to heterogeneity is strongly enhanced compared to homogeneous Biot models. Simple analytical expressions for the P -wave phase velocity and attenuation coefficient are presented as functions of frequency and of statistical medium parameters (correlation lengths, variances). These results automatically include different asymptotic approximations, such as poroelastic Backus averaging in the quasi-static and the no-flow limits, geometrical optics, and intermediate frequency ranges.