Velocity structure of upper-mantle transition zones beneath central eurasia from seismic inversion using genetic algorithms

We present velocity constraints for the upper-mantle transition zones beneath Central Siberia based on observations of the 1982 RIFT Deep Seismic Sounding (DSS) profile. The data consist of seismic recordings of a nuclear explosion in north-western Siberia along a 2600 km long seismic profile extending from the Yamal Peninsula to Lake Baikal. We invert seismic data from the mantle transition zones using a non-linear inversion scheme using a genetic algorithm for optimization and the WKBJ method to compute the synthetic seismograms. A statistical error analysis using a graph-binning technique was performed to provide uncertainty values in the velocity models.
Our best model for the upper-mantle velocity discontinuity near 410 km depth has a two-stage velocity-gradient structure, with velocities increasing from 8.70–9.25 km s⁻¹ over a depth range of 400–415 km, a gradient of 0.0433 s⁻¹, and from 9.25–9.60 km s⁻¹ over a depth range of 415–435 km, a gradient of 0.0175 s⁻¹. This derived model is consistent with other seismological observations and mineral-physics models. The model for the velocity discontinuity near 660 km depth is simple, sharp and includes velocities increasing from 10.15 km s⁻¹ at 655 km depth to 10.70 km s⁻¹ at 660 km depth, a gradient of 0.055 s⁻¹.     

Dynamic stress field of a kinematic earthquake source model with k-squared slip distribution

Source models such as the k -squared stochastic source model with k -dependent rise time are able to reproduce source complexity commonly observed in earthquake slip inversions. An analysis of the dynamic stress field associated with the slip history prescribed in these kinematic models can indicate possible inconsistencies with physics of faulting. The static stress drop, the strength excess, the breakdown stress drop and critical slip weakening distance D _c distributions are determined in this study for the kinematic k -squared source model with k -dependent rise time. Several studied k -squared models are found to be consistent with the slip weakening friction law along a substantial part of the fault. A new quantity, the stress delay, is introduced to map areas where the yielding criterion of the slip weakening friction is violated. Hisada's slip velocity function is found to be more consistent with the source dynamics than Boxcar, Brune's and Dirac's slip velocity functions. Constant rupture velocities close to the Rayleigh velocity are inconsistent with the k -squared model, because they break the yielding criterion of the slip weakening friction law. The bimodal character of D _c / D _tot frequency–magnitude distribution was found. D _c approaches the final slip D _tot near the edge of both the fault and asperity. We emphasize that both filtering and smoothing routinely applied in slip inversions may have a strong effect on the space–time pattern of the inferred stress field, leading potentially to an oversimplified view of earthquake source dynamics.     

A seismic study of deep geological structure in the Bristol Channel area, SW Britain

Summary. Results from eight seismic refraction lines, 35–90 km long, in the Bristol Channel area are presented. The data, mostly land recordings of marine shots, have been interpreted mainly by ray-tracing and time-term modelling. Upper layer velocities through Palaeozoic rocks usually fall within the range 4.8–5.2 km s⁻¹. Below the Carboniferous Limestone with a normal velocity of 5.1–5.2 kms⁻¹, the Old Red Sandstone with a velocity of 4.7–4.8 kms⁻¹ acts as a low velocity layer, as do parts of the underlying Lower Palaeozoic succession. In the central South Wales/Bristol Channel area and the Mendips, a 5.4–5.5 km s⁻¹ refractor is correlated with a horizon at or near the top of the Lower Palaeozoic succession. Under the whole area, except for north Devon, a 6.0–6.2 km s⁻¹ basal refractor has been located and is correlated with Precambrian crystalline basement rocks. In general, this refractor deepens southwards from a series of basement highs, which existed before the major movements of the Variscan orogeny in South Wales, resulting in a southerly thickening of the pre Upper Carboniferous supra-basement sequence. In north Devon, a 6.2 km s⁻¹ refractor at shallow depth, interpreted as a horizon in the Devonian or Lower Palaeozoic succession, overlies a deep reflector that may represent the Precambrian crystalline basement.     

Experimental study of the anisotropy of longitudinal and transverse waves from local earthquake records

Summary. The paper gives the results of a study of the anisotropy of seismic wave velocities within the Ashkhabad test field in Central Asia. The anisotropy was studied by analysing variations in the values of apparent velocities of first arrivals for epicentral distances ranging from 30 to 130 km and by analysing the delays (Δ t_s1-_s2 ) between the arrival times of shear waves with different polarizations.
The velocities of P -waves vary with azimuth from 5.3 to 6.27 km s^-1 and the velocities of S -waves vary from 3.15 to 3.5 km s^-1.
The delay times Δ t_S1 - _S2 depend on the direction of the propagation. The character of the variation of the propagation velocity of the longitudinal wave, the presence of two differently polarized shear waves S 1 and S 2 propagating at different velocities, and the character of the distribution of Δ t_S1 - _S2 on the stereogram suggest that the symmetry of the anisotropic medium is close to hexagonal with a nearly horizontal symmetry axis coinciding with the direction of maximal velocity. The azimuth of the symmetry axis of the medium is 140° and coincides with the direction of geological faults.     

Analysis of the Spatial Coherence of Short-Period Acoustic-Gravity Waves in the Atmosphere

Summary The coherence of atmospheric acoustic-gravity waves has been measured in the period range 10–100 s at the Large Aperture Microbarograph Array in south-eastern Montana. The acoustic-gravity waves observed were signals generated by presumed nuclear explosions. The decrease of coherence with increasing distance between pairs of microbarographs is less rapid in the direction of wave propagation than transverse to it. Variation of direction of arrival over a small range of azimuth (±5°) explains the spatial behaviour of coherence in the direction normal to the wave propagation; variation of phase velocity of ±10 ms^-1 explains the behaviour along the direction of wave propagation. Both effects may be due to inhomogeneities in the atmosphere; the velocity variation may be due to the presence in the signal of several normal modes of acoustic- gravity waves, each travelling at a slightly different phase velocity in the range 300–330 ms^-1.     

Imaging the Hikurangi subduction zone, New Zealand, using teleseismic receiver functions: crustal fluids above the forearc mantle wedge

We image the Hikurangi subduction zone using receiver functions derived from teleseismic earthquakes. Migrated receiver functions show a northwest dipping low shear wave feature down to 60 km depth, which we associate with the crust of the subducted Pacific Plate. Receiver functions (RF) at several stations also show a pair of negative and positive polarity phases with associated conversion depths of ∼20–26 km, where the subducted Pacific Plate is at a depth of ∼40–50 km beneath the overlying Australian Plate. RF inversion solutions model these phases with a thin low S -wave velocity zone less than 4 km thick, and an S -wave velocity contrast of more than ∼0.5 km s⁻¹ with the overlying crust. We interpret this phase pair as representing fluids near the base of the lower crust of the Australian Plate, directly overlying the forearc mantle wedge.     

Waveform modelling of teleseismic S, Sp, SsPmP, and shear-coupled PL waves for crust- and upper-mantle velocity structure beneath Africa

We describe a waveform modelling technique and demonstrate its application to determine the crust- and upper-mantle velocity structure beneath Africa. Our technique uses a parallelized reflectivity method to compute synthetic seismograms and fits the observed waveforms by a global optimization technique based on a Very Fast Simulated Annealing (VFSA). We match the S , Sp, SsPmP and shear-coupled PL phases in seismograms of deep (200–800 km), moderate-to-large magnitude (5.5–7.0) earthquakes recorded teleseismically at permanent broad-band seismic stations in Africa. Using our technique we produce P - and S -wave velocity models of crust and upper mantle beneath Africa. Additionally, our use of the shear-coupled PL phase, wherever observed, improves the constraints for lower crust- and upper-mantle velocity structure beneath the corresponding seismic stations. Our technique retains the advantages of receiver function methods, uses a different part of the seismogram, is sensitive to both P - and S -wave velocities directly, and obtains helpful constraints in model parameters in the vicinity of the Moho. The resulting range of crustal thicknesses beneath Africa (21–46 km) indicates that the crust is thicker in south Africa, thinner in east Africa and intermediate in north and west Africa. Crustal P - (4.7–8 km s⁻¹) and S -wave velocities (2.5–4.7  km s⁻¹) obtained in this study show that in some parts of the models, these are slower in east Africa and faster in north, west and south Africa. Anomalous crustal low-velocity zones are also observed in the models for seismic stations in the cratonic regions of north, west and south Africa. Overall, the results of our study are consistent with earlier models and regional tectonics of Africa.     

Crustal and upper mantle structure of the northwestern North Island, New Zealand, from seismic refraction data

The crustal and upper mantle structure of the northwestern North Island of New Zealand is derived from the results of a seismic refraction experiment; shots were fired at the ends and middle of a 575 km-long line extending from Lake Taupo to Cape Reinga. The principal finding from the experiment is that the crust is 25 ± 2 km thick, and is underlain by what is interpreted to be an upper mantle of seismic velocity 7.6 ± 0.1 km s⁻¹, that increases to 7.9 km s⁻¹ at a depth of about 45 km. Crustal seismic velocities vary between 5.3 and 6.36 km s⁻¹ with an average value of 6.04 km s⁻¹. There are close geophysical and geological similarities between the north-western North Island of New Zealand and the Basin and Range province of the western United States. In particular, the conditions of low upper-mantle seismic velocities, thin crust with respect to surface elevation, and high heat-flow (70–100 mW m⁻²) observed in these two areas can be ascribed to their respective positions behind an active convergent margin for about the past 20 Myr.     

Upper mantle shear structure of North America

Summary. The waveforms and travel times of S and SS phases in the range 10°–60° have been used to derive upper mantle shear velocity structures for two distinct tectonic provinces in North America. Data from earthquakes on the East Pacific Rise recorded at stations in western North America were used to derive a tectonic upper mantle model. Events on the north-west coast of North America and earthquakes off the coast of Greenland provided the data to investigate the upper mantle under the Canadian shield. All branches from the triplications due to velocity jumps near 400 and 660 km were observed in both areas. Using synthetic seismograms to model these observations placed tight constraints on heterogeneity in the upper mantle and on the details of its structure. SS–S travel-time differences of 30 s along with consistent differences in waveforms between the two data sets require substantial heterogeneity to at least 350 km depth. Velocities in the upper 170 km of the shield are about 10 per cent higher than in the tectonic area. At 250 km depth the shield velocities are still greater by about 4.5 per cent and they gradually merge near 400 km. Below 400 km no evidence for heterogeneity was found. The two models both have first-order discontinuities of 4.5 per cent at 405 km and 7.5 per cent at 695 km. Both models also have lids with lower velocities beneath. In the western model the lid is very thin and of relatively low velocity. In the shield the lid is 170 km thick with very high elocity (4.78 km s^-1); below it the velocity decreases to about 4.65 km s^-1. Aside from these features the models are relatively smooth, the major difference between them being a larger gradient in the tectonic region from 200 to 400 km.     

Coseismic slip model of the 2007 August Pisco earthquake (Peru) as constrained by Wide Swath radar observations

The Pisco earthquake ( M _w 8.0; 2007 August 15) occurred offshore of Peru's southern coast at the subduction interface between the Nazca and South American plates. It ruptured a previously identified seismic gap along the Peruvian margin. We use Wide Swath InSAR observations acquired by the Envisat satellite in descending and ascending orbits to constrain coseismic slip distribution of this subduction earthquake. The data show movement of the coastal regions by as much as 85 cm in the line-of-sight of the satellite. Distributed-slip model indicates that the coseismic slip reaches values of about 5.5 m at a depth of ∼18–20 km. The slip is confined to less than 40 km depth, with most of the moment release located on the shallow parts of the interface above 30 km depth. The region with maximum coseismic slip in the InSAR model is located offshore, close to the seismic moment centroid location. The geodetic estimate of seismic moment is 1.23 × 10²¹ Nm ( M _w 8.06), consistent with seismic estimates. The slip model inferred from the InSAR observations suggests that the Pisco earthquake ruptured only a portion of the seismic gap zone in Peru between 13.5° S and 14.5° S, hence there is still a significant seismic gap to the south of the 2007 event that has not experienced a large earthquake since at least 1687.     

Kinematic multistation discriminator between local quarry blasts and earthquakes

It has long been known that S waves on seismograms of local explosions are often accompanied by strong low-frequency, low-velocity, Rg surface wave trains, often significantly diminished for earthquakes. We utilize this fact to construct a new formal discriminator between earthquakes and explosions by measuring the S -surface-wave group velocity. The method is based on analysing the velogram; that is, the display of the envelope of ground motion versus group velocity V = R/T , where R is the epicentral distance and T  the traveltime. We examine the distribution of seismic energy in time and space using envelopes of records from the Israel Seismic Network (ISN), from which we compute the velograms and observe differences in the velograms of quarry blasts and earthquakes. The data include 143 seismic events occurring in three areas (Galilee, Dead Sea, and Gilad) monitored by the ISN; the magnitude range is M _L = 1.0–2.8 at distances of 15–310 km. From the velograms we measure the group velocity, V _{m s} , within the 1–4 km s⁻¹ range at which the velogram reaches its maximum for each available station. The resulting V _{m s} (R) function is closely fitted by the empirical relationship a + b  ln  R , with a and b coefficients varying from event to event. A simple linear function c = b + 0.33a at a threshold C = 0.69 completely separates ( a,b ) pairs for the 67 Galilee events, and, for the 76 remaining events, one earthquake and four explosions are wrongly classified. After data validation and application of the Fisher linear discriminator, adapted to the events from Galilee, only two misidentified events remain for the whole data set.     

Transitional continental–oceanic structure beneath the Norwegian Sea from inversion of surface wave group velocity data

We have analysed the fundamental mode of Love and Rayleigh waves generated by 12 earthquakes located in the mid-Atlantic ridge and Jan Mayen fracture zone. Using the multiple filter analysis technique, we isolated the Rayleigh and Love wave group velocities for periods between 10 and 50  s. The surface wave propagation paths were divided into five groups, and average group velocities calculated for each group. The average group velocities were inverted and produced shear wave velocity models that correspond to a quasi-continental oceanic structure in the Greenland–Norwegian Sea region. Although resolution is poor at shallow depth, we obtained crustal thickness values of about 18  km in the Norwegian Sea area and 9  km in the region between Svalbard and Iceland. The abnormally thick crust in the Norwegian Sea area is ascribed to magmatic underplating and the thermal blanketing effect of sedimentary layers. Maximum crustal shear velocities vary between 3.5 and 3.9  km  s⁻¹ for most paths. An average lithospheric thickness of 60  km was observed, which is lower than expected for oceanic-type structure of similar age. We also observed low shear wave velocities in the lower crust and upper mantle. We suggest that high heat flow extending to depths of about 30  km beneath the surface can account for the thin lithosphere and observed low velocities. Anisotropy coefficients of 1–5 per cent in the shallow layers and >7 per cent in the upper mantle point to the existence of polarization anisotropy in the region.     

A surge of Perseibreen, Svalbard, examined using aerial photography and ASTER high resolution satellite imagery

The identification of surge activity is important in assessing the duration of the active and quiescent phases of the surge cycle of Svalbard glaciers. Satellite and aerial photographic images are used to identify and describe the form and flow of Perseibreen, a valley glacier of 59 km² on the east coast of Spitsbergen. Heavy surface crevassing and a steep ice front, indicative of surge activity, were first observed on Perseibreen in April 2002. Examination of high resolution (15 m) Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) satellite imagery confirmed this surge activity. Perseibreen retreated by almost 750 m between 1961 and 1990. Between 1990 and the summer of 2000, Perseibreen switched from retreat and its front began to advance. Rapid advance was underway during the period June 2000 to May 2001, with terminus advance at over 400 m yr⁻¹. Between May and August 2001 the rate increased to over 750 m yr⁻¹. The observed crevasse orientation indicates that ice was in longitudinal tension, suggesting the down-glacier transfer of mass. Ice surface velocities, derived from image correlation between ASTER images, were 2-2.5 m d⁻¹ between May and August 2001. The glacier was flowing at a relatively uniform speed with sharp velocity gradients located close to its lateral margins, a velocity structure typical of ice masses in the active phase of the surge cycle. The stress regime is extensional throughout and the surge appears to be initiated low on the glacier. This is similar to the active-phase dynamics of other Svalbard tidewater glaciers. Perseibreen has probably been inactive since at least 1870, a period of about 130 years to the present surge which defines a minimum length for the quiescent phase.     

Upper crustal shear velocity models from higher mode Rayleigh wave dispersion in Scotland

Summary. Group velocities for first and second higher mode Rayleigh waves, in the frequency range 0.8–4.8 Hz, generated from a local earthquake of magnitude 3.7 M _L in western Scotland, are measured at stations along the 1974 LISPB line. These provide detailed information about the crustal structure west of the line. The data divide the region into seven apparently homogeneous provinces. Averaged higher mode velocity dispersion curves for each province are analysed simultaneously using a linearized inversion technique, yielding regionalized shear velocity profiles down to a depth of 17 km into the upper crust. Shear wave velocity is between 3.0 and 3.4 km s⁻¹ in the upper 2 km, with a slow increase to around 3.8 km s⁻¹. P -wave models computed using these results agree with profiles from the LISPB and LUST refraction experiments.     

Crustal structure in the Western Somali Basin

Summary. As part of integrated marine geophysical studies in the Western Somali Basin, we performed 118 sonobuoy experiments to define better the crustal structure of the margins and basin created by the separation of Madagascar and Africa. After using T ²/ X ², conventional slope-intercept methods, and slant-stacked t-p techniques to analyse the data, we combined our solutions with all previous velocity information for the area. Velocity functions were derived for the sediment coiumn, and we detected a high-velocity (4.58 ± 0.29 km s^–1) sediment layer overlying acoustic basement. We confirmed that the crust is indeed seismically oceanic, and that it may be considered either in terms of a layered model – layers 2B (5.42 ± 0.19 km s^–1), 2C (6.23 ± 0.22 km s^–1), 3 (7.03 ± 0.25 km s^–1), and mantle (7.85 ± 0.32 km s^–1) were identified – or a more complex gradient model in which layer 2 is marked by a steeper velocity gradient than underlying layer 3. Integrated igneous crustal thicknesses (1.62 ± 0.22 s, 5.22 ± 0.64 km) are significantly less than what is considered normal. We present a revised seismic transect across the East African margin, as well as total sediment thickness, depth to basement and crustal thickness maps.     

Comparison of the Upper Mantle Shear Wave Velocity Structure of the Indian Shield and the Tibetan Plateau and Tectonic Implications

Summary. Travel times and waveforms of long-period SH -waves recorded at distances of 10–30° and some SS waveforms are used to constrain the upper mantle velocities down to a depth of 400km beneath both the Indian Shield and the Tibetan Plateau. the shear velocity in the uppermost mantle beneath both the Indian Shield and the Tibetan Plateau is high and close to 4.7 km s⁻¹. the Indian Shield has a fairly thick high velocity lid, and the mean velocity between 40 and 250 km is between 4.58 and 4.68 km s⁻¹. In contrast, S -wave travel times and waveforms of S -waves, as well as a few for SS , show that the mean velocity between 70 and 250km beneath the central and northern part of the Tibetan Plateau is slower by 4 per cent or more than that beneath the Indian Shield and probably is between 4.4 and 4.5km s⁻¹. No large differences in the structure of the two areas below 250 km are required to explain both the arrival times and the waveforms of SH phases crossing Tibet or the Indian Shield. These results show that the structure of Tibet is not that of a shield and imply that the Indian plate is not underthrusting the whole of the Tibetan Plateau at the present time.     

Crustal structure of the Mid-Atlantic ridge crest at 37° N

Summary. A structural model of the Mid-Atlantic Ridge at 37° N is proposed on the basis of travel-time data and synthetic seismograms. At the ridge axis the crust is only 3 km thick and overlies material with an anomalously low'upper mantle'velocity of 7.2 km s⁻¹. Crustal thickening and the formation of layer 3 and a layer with velocity 7.2–7.3 km s⁻¹ takes place within a few kilometres of the axis, producing a 6–7 km thick crust by less than 10 km from the axis. A normal upper mantle velocity of 8.1 km s⁻¹ exists within 10 km of the axis. Shear waves propagate across the axis, thus precluding the existence of any sizeable magma chamber at shallow depth.     

Strain accumulation along an oblique plate boundary: the Reykjanes Peninsula, southwest Iceland

We use annual GPS observations on the Reykjanes Peninsula (RP) from 2000 to 2006 to generate maps of surface velocities and strain rates across the active plate boundary. We find that the surface deformation on the RP is consistent with oblique plate boundary motion on a regional scale, although considerable temporal and spatial strain rate variations are observed within the plate boundary zone. A small, but consistent increase in eastward velocity is observed at several stations on the southern part of the peninsula, compared to the 1993–1998 time period. The 2000–2006 velocities can be modelled by approximating the plate boundary as a series of vertical dislocations with left-lateral motion and opening. For the RP plate boundary we estimate left-lateral motion  18⁺⁴₋₃ mm yr⁻¹  and opening of  7⁺³₋₂ mm yr⁻¹  below a locking depth of  7⁺¹₋₂ km  . The resulting deep motion of  20⁺⁴₋₃ mm yr⁻¹  in the direction of  N(100⁺⁸₋₆)°E  agrees well with the predicted relative North America–Eurasia rate. We calculate the areal and shear strain rates using velocities from two periods: 1993–1998 and 2000–2006. The deep motion along the plate boundary results in left-lateral shear strain rates, which are perturbed by shallow deformation due to the 1994–1998 inflation and elevated seismicity in the Hengill–Hrómundartindur volcanic system, geothermal fluid extraction at the Svartsengi power plant, and possibly earthquake activity on the central part of the peninsula.     

Stress accumulation between volcanoes: an explanation for intra-arc earthquakes in Nicaragua?

Destructive upper crustal earthquakes in Central America are often located between active volcanic centres—a geometric relationship that we study using finite element Coulomb failure stress (CFS) models that incorporate the rheologically heterogeneous nature of the volcanic arc. Volcanoes are simulated as mechanically weak zones within a stronger crust. We find that deformation of the volcanic centres within a regional stress field dominated by dextral shear causes stress increases in surrounding crust, with a maximum CFS change between neighbouring volcanoes. This increase in CFS enhances the probability of fault slip on arc-normal faults that are located between volcanic centres; for example, the Tiscapa fault, which ruptured during the 1972 December 13,   M _s   6.2 Managua earthquake. The amount of stress increase due to long-term (100 yr) volcano shearing is on the order of 0.1–0.6 bars, similar to values estimated for subduction zone earthquakes.