Retrieving shear velocity structure from high-frequency toroidal modes of the Earth

Summary. For a smooth earth model, observations of a set of high-frequency toroidal modes at fixed slowness yield only a single piece of information, the tau value for that slowness. In this note, a procedure for obtaining the shear velocity structure from free oscillation data for an earth model with velocity discontinuities is developed, based on the method of tau inversion. The information content of the high-frequency modes is greater in this case, and the nature and depths of the discontinuities may be deduced. It is shown, for the real Earth, that the tau values obtained from free oscillation data are affected significantly by the presence of the Moho, but a simple iterative scheme may be used to remove this contamination. Brune's method of deducing mode frequencies from body wave pulses is shown to produce significant errors for a model with a pronounced Moho discontinuity, and the same iterative scheme may also be employed to correct for this effect.     

Pre-drilling reflection survey of the Black Forest, SW Germany

Summary. The unified seismic exploration program, consisting of 345 km of deep reflection profiling, a 200 km refraction profile, an expanding spread profile and near-surface high resolution reflection meaasurements, revealed a strongly differentiated crust beneath the Black Forest. The highly reflective lower crust contains numerous horizontal and dipping reflectors at depths of 13-14 km down to the crust-mantle boundary (Moho). The Moho appears as a flat horizontal first order discontinuity at a relatively shallow level of 25–27 km above a transparent upper mantle. From modelling of synthetic near-vertical and wide-angle seismograms using the reflectivity method the lower crust is supposed to be composed of laminae with an average thickness of about 100 m and velocity differences of greater than 10% increasing from top to bottom. The upper crust is characterised by mostly dipping reflectors, associated with bivergent underthrusting and accretion tectonics of Variscan age and with extensional faults of Mesozoic age. A bright spot at 9.5 km depth is characterised by low velocity material suggesting a fluid trap. It appears on all of the three profiles in the centre of the intersection region. The upper crust seems to be decoupled from the lowest crust by a relatively transparent zone which is' also identified as a low-velocity zone. This low velocity channel is situated directly above the laminated lower crust. The laminae in the Rhinegraben area are displaced vertically to greater depths indicating an origin before Tertiary rift formation and a subsidence of the whole graben wedge.     

Crustal structure of Iraq from receiver functions and surface wave dispersion: implications for understanding the deformation history of the Arabian–Eurasian collision

We report the crustal structure for two locations in Iraq estimated by joint inversion of P -wave receiver functions (RFs) and surface (Rayleigh) wave group velocity dispersion. RFs were computed from teleseismic recordings at two temporary broad-band seismic stations located in Mosul (MSL) in the Zagros Fold Belt and Baghdad (BHD) in the Mesopotamian Foredeep. Group velocity dispersion curves at the sites were derived from continental-scale tomography. The inversion results show that the crustal thicknesses are 39 km at MSL and 43 km at BHD. We observe a strong Ps _Moho at BHD consistent with a sharp Moho discontinuity. However, at MSL we observe a weak Ps _Moho suggesting a transitional Moho where crustal thickening is likely to be occurring in the deep crust. Both sites reveal low velocity surface layers consistent with sedimentary thickness of about 3 km at station MSL and 7 km at BHD and agreeing well with the previous reports. Ignoring the sediments, the crystalline crustal velocities and thicknesses are remarkably similar at both stations. The similarity of crustal structure suggests that the crust of the northeastern proto-Arabian Platform was uniform before subsidence and deposition of the sediments in the Cenozoic. If crystalline crustal structure is uniform across the northern Arabian Platform then crustal thickness variations in the Zagros Fold Belt and Thrust Zone should reveal the history of deformation and crustal shortening in the Arabian–Eurasian collision zone and not reflect pre-existing crustal thickness variations in the Arabian Plate.     

Crustal root beneath the highlands of southern Norway resolved by teleseismic receiver functions

Teleseismic data have been collected with temporary seismograph stations on two profiles in southern Norway. Including the permanent arrays NORSAR and Hagfors the profiles are 400 and 500 km long and extend from the Atlantic coast across regions of high topography and the Oslo Rift. A total of 1071 teleseismic waveforms recorded by 24 temporary and 8 permanent stations are analysed. The depth-migrated receiver functions show a well-resolved Moho for both profiles with Moho depths that are generally accurate within ±2 km.
For the northern profile across Jotunheimen we obtain Moho depths between 32 and 43 km (below sea level). On the southern profile across Hardangervidda, the Moho depths range from 29 km at the Atlantic coast to 41 km below the highland plateau. Generally the depth of Moho is close to or above 40 km beneath areas of high mean topography (>1 km), whereas in the Oslo Rift the crust locally thins down to 32 km. At the east end of the profiles we observe a deepening Moho beneath low topography. Beneath the highlands the obtained Moho depths are 4–5 km deeper than previous estimates. Our results are supported by the fact that west of the Oslo Rift a deep Moho correlates very well with low Bouguer gravity which also correlates well with high mean topography.
The presented results reveal a ca . 10–12 km thick Airy-type crustal root beneath the highlands of southern Norway, which leaves little room for additional buoyancy-effects below Moho. These observations do not seem consistent with the mechanisms of substantial buoyancy presently suggested to explain a significant Cenozoic uplift widely believed to be the cause of the high topography in present-day southern Norway.     

Lithospheric structure of the Arabian and Eurasian collision zone in eastern Turkey from S-wave receiver functions

Crustal and upper-mantle seismic discontinuities beneath eastern Turkey are imaged using teleseismic S -to- P converted phases. Three crustal phases are observed: the Moho with depth ranging between 30 and 55 km, indicating variable tectonic regimes within this continental collision zone; an upper-crustal discontinuity at approximately 10 km depth; and various crustal low-velocity zones, possibly associated with recent Quaternary volcanism. Imaging of the upper mantle is complicated by the 3-D geometry of the region, in particular due to the Bitlis–Zagros suture zone. However, several upper-mantle S -to- P converted phase are identified as being the signature of the lithosphere–asthenosphere boundary (LAB). The inferred LAB for the Eastern Anatolian Accretionary Complex indicates that eastern Turkey has an anomalously thin (between ∼60 and 80 km) lithosphere which is consistent with an oceanic slab detachment model. The observed LAB phases for the Arabian shield and Iranian plateau indicate that lithospheric thickness for these stable regions is on the order of 100 to 125 km thick, which is typical of continental margins.     

Three-dimensional image of the Moho undulations beneath the Gulf of Bothnia using wide-angle seismic data

The BABEL marine seismic experiment has been carried out to investigate the lithospheric structure and antecedent tectonic signatures of the Baltic Shield, including the Archaean-Proterozoic collisional structure in the northern part of the Gulf of Bothnia.
Lithospheric seismic-reflection streamer data and simultaneously recorded wide-angle reflection and refraction data collected in the Gulf of Bothnia as part of the BABEL project have been used for 3-D modelling. The distribution of land stations around the Gulf provides a good 3-D ray coverage of the PMP reflection data recorded at the eight stations in the area and allows an estimation of strikes and dips of the Moho boundary in the area. The traveltimes of reflected phases are calculated using a method that utilizes the finite-difference solution of the eikonal equation. The Moho wide-angle-reflection (PMP) traveltimes are modelled using an inversion method. A 2-D model from the Gulf of Bothnia extended into the third dimension is used as an initial model. During the inversion the velocity is kept constant and only the Moho boundary is allowed to vary. To estimate the strike of the Moho boundary and the stability of the inversion, two initial models with different strikes are examined.
The results indicate that the Moho depth in the Gulf of Bothnia undulates and has a maximum depth of 55 km in the south, rising to 42 km in the north. The Moho depth variations seem to be step-like. This change in the Moho depth coincides with the location of the presumed fossil subduction zone in the area. The crustal-thickness variations seem to be well approximated by a nearly 2-D structure striking parallel to a postulated subduction zone immediately to the south of the Skellefte area. The presence of the step at the crust/mantle boundary can be interpreted as a result of a plate-collision event at about 2 Ga.     

Gravity anomalies, subsidence history and the tectonic evolution of the Malay and Penyu Basins (offshore Peninsular Malaysia)

The tectonic subsidence and gravity anomalies in the Malay and Penyu Basins, offshore Peninsular Malaysia, were analysed to determine the isostatic compensation mechanism in order to investigate their origin. These continental extensional basins contain up to 14  km of sediment fill which implies that the crust had been thinned significantly during basin development. Our results suggest, however, that the tectonic subsidence in the basins cannot be explained simply by crustal thinning and Airy isostatic compensation.
The Malay and Penyu Basins are characterized by broad negative free-air gravity anomalies of between −20 and −30  mGal. To determine the cause of the anomaly, we modelled four gravity profiles across the basins using a method that combines two-dimensional flexural backstripping and gravity modelling techniques. We assumed a model of uniform lithospheric stretching and Airy isostasy in the analysis of tectonic subsidence. Our study shows that the basins are probably underlain by relatively thinned crust, indicating that some form of crustal stretching was involved. To explain the observed gravity anomalies, however, the Moho depth that we calculated based on the free-air gravity data is about 25% deeper than the Moho predicted by assuming Airy isostasy (Backstrip Moho). This suggests that the Airy model overestimates the compensation and that the basins are probably undercompensated isostatically. In other words, there is an extra amount of tectonic subsidence that is not compensated by crustal thinning, which has resulted in the discrepancy between the gravity-derived Moho and the Backstrip Moho. We attribute this uncompensated or anomalous tectonic subsidence to thin-skinned crustal extension that did not involve the mantle lithosphere. The Malay and Penyu Basins are interpreted therefore as basins that formed by a combination of whole-lithosphere stretching and thin-skinned crustal extension.     

Forward modelling receiver functions for crustal structure beneath station TBZ (Trabzon, Turkey)

We use teleseismic three-component digital data from the Trabzon, Turkey broadband seismic station TBZ to model the crustal structure by the receiver function method. The station is located at a structural transition from continental northeastern Anatolia to the oceanic Black Sea basin. Rocks in the region are of volcanic origin covered by young sediments. By forward modelling the radial receiver functions, we construct 1-D crustal shear velocity models that include a lower crustal low-velocity zone, indicating a partial melt mechanism which may be the source of surfacing magmatic rocks and regional volcanism. Within the top 5 km, velocities increase sharply from about 1.5 to 3.5 km s⁻¹. Such near-surface low velocities are caused by sedimentation, extending from the Black Sea basin. Velocities at around 20 km depth have mantle-like values (about 4.25 km s⁻¹ ), which easily correlate to magmatic rocks cropping out on the surface. At 25 km depth there is a thin low-velocity layer of about 4.0 km s⁻¹. The average Moho velocity is about 4.6 km s⁻¹, and its depth changes from 32 to 40 km. Arrivals on the tangential components indicate that the Moho discontinuity dips approximately southwards, in agreement with the crustal thickening to the south. We searched for the solution of receiver functions around the regional surface wave group velocity inversion results, which helped alleviate the multiple solution problem frequently encountered in receiver function modelling.
Station TBZ is a recently deployed broadband seismic station, and the aim of this study is to report on the analysis of new receiver function data. The analysis of new data in such a structurally complex region provides constraining starting models for future structural studies in the region.     

Deep seismic reflectors in the Campos basin, offshore Brazil

Summary. Some deep crustal features underlying the Campos basin are best recognized in a few reflection seismic sections that have been reprocessed recently to 10 s two-way traveltime. A prominent climbing-to-the-basin reflector is interpreted as the Moho, and a relatively steep fracture zone is, probably, the first example so far of an extensional fault crossing the whole crust and offsetting the Moho. Further constraints on the deep structure of the basin are provided by estimating the thinning of the crust from shallow seismic data and gravity modelling, and by cross-plotting backstripped subsidence curves against curves predicted by the lithospheric stretching model.     

Moho-offset beneath northern Scotland

Summary. Independent sets of LISPB data are presented, having as a common feature a local lateral variation in the Moho beneath northern Britain.
The evidence for this step-like feature on the Moho is taken from P and S travel times and from PS reflection times and phase velocities.
Another important observation is that a high-frequency Moho refraction is generated north of the step, whilst the refraction from the step region and south of it contains only low frequencies. This leads to the assumption that a change in the Moho structure exists in the same region as where the step has been deduced from travel-time interpretation. A numerical model is presented which generates a high-frequency refracted wave from an input signal which contains low- and high-frequency spectral energy.
The significance of this offset on the Moho, in particular its relationship to surface tectonics, is critically discussed.     

Determination of Moho dip using PS reflections

Summary. A clear PS reflection was recorded in the northern section of the LISPB line (shot-point N1, profile ALPHA). Travel times of PS and PP phases reflected at the same part of the Moho were used to model the lower crustal structure. The observed high apparent velocity of the PS phase (8.4 km/s) is interpreted as due to a localized Moho dip of 7 per cent (4°) towards the north. This is supported by observations of Pn refractions from other shots. It is shown that PS reflections could be very useful for determination of Moho dips.     

Constraining velocity and density contrasts across the crust–mantle boundary with receiver function amplitudes

This paper investigates the ability of P -wave receiver functions to constrain both the velocity and density contrasts across the Moho. Expressions for receiver function amplitudes corresponding to a layer over a half-space are extended to media with depth-dependent properties by explicitly accounting for free-surface reflection coefficients. Forward numerical computations show that receiver function amplitudes become frequency-dependent for depth-dependent structures and that, for a given frequency, wavelengths associated to Ps converted phases are larger than wavelengths associated to multiply reverberated phases. When consistent measurements of the amplitudes of Ps phase and multiples can be obtained, the extended expressions suffice to account for the observed amplitudes. Receiver function amplitudes are sensitive to both velocity and density contrasts across the Moho, and a two-step grid-search procedure is proposed to recover the contrasts from their measurement. The near-surface velocity is recovered in the first step from the amplitude of the direct P wave and then used in the second step to work out the relative density and S -wave velocity contrast from the amplitude of the secondary phases. Examples from central Spain (PAB) and the Indian Shield (HYB) are discussed and demonstrate that receiver function amplitudes can constrain contrasts across the Moho accurately enough to be utilized in geological interpretation.     

The long-term thermal consequences of rifting: implications for basin reactivation

The attenuation of the continental crust during rifting and the subsequent filling of the rift‐related accommodation alter the long‐term thermal and mechanical state of the lithosphere. This is primarily because the Moho is shallowed due to density contrasts between the sediment fill and the crust, but also reflects the attenuation of the pre‐existing crustal heat production and its burial beneath the basin, as well the thermal properties of the basin fill. Moho shallowing and attenuation of pre‐existing heat production contribute to long‐term cooling of the Moho and thus lithospheric strengthening, as has been pointed out in many previous studies. In contrast, basin filling normally contributes to significant Moho heating allowing the possibility of long‐term lithospheric weakening, the magnitude of which is dependent on the thermal properties of the basin‐fill and the distribution of heat sources in the crust. This paper focuses on the thermal property structure of the crust and basin‐fill in effecting long‐term changes in lithospheric thermal regime, with particular emphasis on the distribution of heat producing elements in the crust. The parameter space appropriate to typical continental crust is explored using a formalism for the heat production distributions that makes no priori assumptions about the specific form of the distribution. The plausible parameter space allows a wide range in potential long‐term thermal responses. However, with the proviso that the accommodation created by the isostatic response to rifting is essentially filled, the long‐term thermal response to rift basin formation will generally increase average crustal thermal gradients beneath basins but cool the Moho due to its reduction in depth. The increase in the average crustal thermal gradient induces lateral heat flow that necessarily heats the Moho along basin margins, especially in narrow rift basins. Using coupled thermo‐mechanical models with temperature sensitive creep‐parameters, we show that such heating may be sufficient to localise subsequent deformation in the vicinity of major basin bounding structures, potentially explaining the offset observed in some stacked rift basin successions.     

Analysis of gravity fields in the northwestern Himalayas and Kohistan region using deep seismic sounding data

A Bouguer gravity anomaly map of the NW Himalayas and parts of the Kohistan/Hindukush region has been prepared using all available gravity data. Analysis of the gravity field has been carried out along a profile extending from Gujranwala (located near the edge of the Indian shield) to the Haramosh massif in a NNE–SSW direction. The gravity profile is located close to the DSS profile shot under the USSR–India scientific collaborative programme. Velocity information available along different parts of the profile has been used to infer values of crustal and upper mantle density.
The observed gravity field (Bouguer) has been interpreted in terms of Moho depth and density contrast between the crust and the mantle. The Moho depth is interpreted as increasing from nearly 35 km near the edge of the Indian shield to 75 km (below sea-level) underneath the Haramosh massif. A similar model is applicable to a profile passing to the west of Nanga Parbat massif, from Gujranwala to Ghizar, through the Kohistan region. However, along this profile high-density lower-crustal rocks appear to have been emplaced in the upper part along the main mantle thrust. The nature of isostatic compensation prevailing underneath the Himalayas has been discussed, as has the theory of lithospheric flexure proposed by Karner & Watts and Lyon-Caen & Molnar. It is felt that although these ideas explain the broad features of the Moho configuration as observed in the NW Himalayas, there are significant departures. The role of tectonic forces in shaping the Moho and causing changes in the density of the crust cannot be denied.     

Regional three-dimensional gravity modelling of the NE Atlantic margin

A series of three‐dimensional models has been constructed for the structure of the crust and upper mantle over a large region spanning the NE Atlantic passive margin. These incorporate isostatic and flexural principles, together with gravity modelling and integration with seismic interpretations. An initial isostatic model was based on known bathymetric/topographic variations, an estimate of the thickness and density of the sedimentary cover, and upper mantle densities based on thermal modelling. The thickness of the crystalline crust in this model was adjusted to equalise the load at a compensation depth lying below the zone of lateral mantle density variations. Flexural backstripping was used to derive alternative models which tested the effect of varying the strength of the lithosphere during sediment loading. The models were analysed by comparing calculated and observed gravity fields and by calibrating the predicted geometries against independent (primarily seismic) evidence. Further models were generated in which the thickness of the sedimentary layer and the crystalline crust were modified in order to improve the fit to observed gravity anomalies. The potential effects of igneous underplating and variable upper mantle depletion were explored by a series of sensitivity trials. The results provide a new regional lithospheric framework for the margin and a means of setting more detailed, local investigations in their regional context. The flexural modelling suggests lateral variations in the strength of the lithosphere, with much of the margin being relatively weak but areas such as the Porcupine Basin and parts of the Rockall Basin having greater strength. Observed differences between the model Moho and seismic Moho along the continental margin can be interpreted in terms of underplating. A Moho discrepancy to the northwest of Scotland is ascribed to uplift caused by a region of upper mantle with anomalously low density, which may be associated with depletion or with a temperature anomaly.     

横穿印度洋重力剖面的基本特征及其地质意义

本文根据在印度洋实测的重力资料,用ｓｉｎｘ／ｘ反演了一些典型构造的莫氏面深度,对印度洋两种不同类型的大陆边缘、大洋中脊、无震海岭、深海盆地的重力异常特征及地质意义进行了较深入的研究,结果表明,东南非岸外以张性下沉为主,构造活动比较单一,其上的重力异常为单调的正负不对称异常,说明这儿的被动陆缘基本上处于均衡状态;而印度洋东北边缘的太平洋型活动陆缘,为一复杂的火山岛弧系,这里的重力异常面貌复杂,反映了该地区处于非均衡状态;查戈斯－拉克代夫海岭为印度洋中典型的无震海岭,它不应属于印度洋中脊的一部分;印度洋中的大部分深海盆,其上具有典型的大洋重力组合,说明印度洋存在着广大的大洋型原生地壳。     

Crustal thickness variations in the margins of the Gulf of California from receiver functions

Receiver functions (RFs) from teleseismic events recorded by the NARS-Baja array were used to map crustal thickness in the continental margins of the Gulf of California, a newly forming ocean basin. Although the upper crust is known to have split apart simultaneously along the entire length of the Gulf, little is known about the behaviour of the lower crust in this region. The RFs show clear P -to- S wave conversions from the Moho beneath the stations. The delay times between the direct P and P -to- S waves indicate thinner crust closer to the Gulf along the entire Baja California peninsula. The thinner crust is associated with the eastern Peninsular Ranges batholith (PRB). Crustal thickness is uncorrelated with topography in the PRB and the Moho is not flat, suggesting mantle compensation by a weaker than normal mantle based on seismological evidence. The approximately W–E shallowing in Moho depths is significant with extremes in crustal thickness of ∼21 and 37 km. Similar results have been obtained at the northern end of the Gulf by Lewis et al., who proposed a mechanism of lower crustal flow associated with rifting in the Gulf Extensional Province for thinning of the crust. Based on the amount of pre-Pliocene extension possible in the continental margins, if the lower crust did thin in concert with the upper crust, it is possible that the crust was thinned during the early stages of rifting before the opening of the ocean basin. In this case, we suggest that when breakup occurred, the lower crust in the margins of the Gulf was still behaving ductilely. Alternatively, the lower crust may have thinned after the Gulf opened. The implications of these mechanisms are discussed.     

Crustal thickness estimation beneath the southern central Andes at 30°S and 36°S from S wave receiver function analysis

We have used the S wave receiver function (SRF) technique to investigate the crustal thickness beneath two seismic profiles from the CHARGE project in the southern central Andes. A previous study employing the P wave receiver function method has observed the Moho interface beneath much of the profiles. They found, however, that the amplitude of the P to S conversion was diminished in the western part of the profiles and have attributed it to a reduction of the impedance contrast at the Moho due to lower crustal ecologitization. With SRF, we have successfully detected S to P converted waves from the Moho as well as possible conversions from other lithospheric boundaries. The continental South American crust reaches its maximum thickness of ∼70 km (along 30°S between 70°W and 68.5°W) beneath the Principal Cordillera and the Famatina system and becomes thinner towards the Sierras Pampeanas with a thickness of ∼40 km. Negative phases, possibly related to the base of the continental and oceanic lithosphere, can be recognized in the summation traces at different depths. By comparing our results with data obtained from previous investigations, we are able to further constrain the thickness of the crust and lithosphere beneath the central Andes.     

Tectonic evolution of the Alboran Sea basin

The Alboran Sea is an extensional basin of Neogene age that is surrounded by highly arcuate thrust belts. Multichannel seismic (MCS) reflection profile data suggest the basin has a complex tectonic fabric that includes extensional, compressional and strike-slip structures. The early Miocene history appears to be dominated by graben formation with border faults that are in large part contemporaneous with thrust movements in the external zones of the Betic and Rif mountains. Extension appears to have continued into the late Miocene although the main movements were probably completed by the time of the Messinian ‘salinity crisis’. The Pliocene and younger history of the basin is dominated by infilling of the Messinian topography, gentle subsidence, and extensional, compressional and strike-slip movements. There is evidence from the sea-floor morphology and seismicity patterns that the basin is actively deforming in response to present-day plate motions. Backstripping of well data in the basin margin suggests that the initial extensional event was accompanied by crustal and lithospheric thinning. The depth to Moho inferred from backstripping is greater than the depth expected based on seismic and gravity modelling, suggesting that backstripping underestimates the true amount of thinning. One explanation is that some of the thinning occurred while the crust was above sea level, perhaps as a result of either crustal thickening, or a period of lithospheric heating and thinning, prior to rifting. We found that a model with a ‘normal’ crustal thickness of 31.2 km, a lithospheric thickness of 50 km, and β= 1.4 predicts 0.8 km of initial uplift. These parameters fit the well subsidence data and bring the backstripped Moho into better agreement with the seismic and gravity Moho. The origin of such a thin lithosphere is not constrained by the data, but we believe that it may be a result of the detachment of a cold lithospheric ‘root’ that formed during pre-Neogene collisional orogeny in the region.     

COCORP deep seismic reflection traverses of the U.S. Cordillera

Summary. COCORP seismic reflection traverses of the U.S. Cordillera at 40°N and 48.5°N latitude reveal some fundamental similarities as well as significant differences in reflection patterns. On both traverses, autochthonous crust beneath thin-skinned thrust belts of the eastern part of the Cordillera is unreflective; immediately to the west the Cordilleran interior is very reflective above a flat, prominent reflection Moho. Mesozoic accreted terranes in the western part of the orogen are underlain on both traverses by very complex reflection patterns, in constrast to more easily deciphered patterns beneath areas of Cenozoic accretion. The prominent reflection Moho beneath the orogenic interior on both transects probably evolved through a combination of magmatic and deformational processes during Cenozoic extension. The main differences between the two traverses lie in the reflection patterns of the middle and lower crust in the Cordilleran interior; these differences are probably related to the way Cenozoic extension was accommodated at depth. Laminated middle and lower crust above the reflection Moho in the western Basin and Range (40°N) may be related to magmatism, ductile pure shear and large-scale transposition during Cenozoic extension. By contrast, beneath the eastern Basin and Range (40°N), and the orogenic interior in the NW United States (48.5°N), Cenozoic extension was probably accommodated along dipping deformation zones throughout the crust.