Seismological features of island arc crust as inferred from recent seismic expeditions in Japan

Crustal studies within the Japanese islands have provided important constraints on the physical properties and deformation styles of the island arc crust. The upper crust in the Japanese islands has a significant heterogeneity characterized by large velocity variation (5.5–6.1 km/s) and high seismic attenuation (Qp=100–400 for 5–15 Hz). The lateral velocity change sometimes occurs at major tectonic lines. In many cases of recent refraction/wide-angle reflection profiles, a “middle crust” with a velocity of 6.2–6.5 km/s is found in a depth range of 5–15 km. Most shallow microearthquakes are concentrated in the upper/middle crust. The velocity in the lower crust is estimated to be 6.6–7.0 km/s. The lower crust often involves a highly reflective zone with less seismicity, indicating its ductile rheology. The uppermost mantle is characterized by a low Pn velocity of 7.5–7.9 km/s. Several observations on PmP phase indicate that the Moho is not a sharp boundary with a distinct velocity contrast, but forms a transition zone from the upper mantle to the lower crust. Recent seismic reflection experiments revealed ongoing crustal deformations within the Japanese islands. A clear image of crustal delamination obtained for an arc–arc collision zone in central Hokkaido provides an important key for the evolution process from island arc to more felsic continental crust. In northern Honshu, a major fault system with listric geometry, which was formed by Miocene back arc spreading, was successfully mapped down to 12–15 km.     

东海陆架盆地西湖凹陷岩石圈热流变性质

为了探讨东海陆架盆地西湖凹陷岩石圈热流变性质,本文以实测地温数据为依据,模拟西湖凹陷岩石圈热结构,在此基础上,应用流变学原理模拟确定西湖凹陷岩石圈流变性质。结果表明,西湖凹陷岩石圈为一个冷地壳-热地幔、强地壳-弱地幔的"奶油蛋糕"型岩石圈。西湖凹陷平均地表热流密度为71 m W/m~2,地幔热流密度为40~65 m W/m~2,对地表热流密度的贡献度达73%~79%,地表热流受地幔热流控制,莫霍面温度在700℃左右,热岩石圈平均厚度为66 km。西湖凹陷岩石圈流变分层明显,上、中地壳基本为脆性层,下地壳和岩石圈上地幔为韧性层,岩石圈总流变强度平均约为2.65′10~(12) N/m,其中地壳流变强度为2.12′10~(12) N/m,地幔流变强度为5.29′10~(11) N/m,有效弹性厚度为11.7~14.5 km,地壳的流变性质控制了岩石圈的流变行为。此外,西湖凹陷岩石圈总强度较低,在构造应力作用下易于变形,且存在壳幔解耦现象。西湖凹陷岩石圈热状态及流变性质决定了西湖凹陷东部地区主要以浅部地壳的断层滑动和地层破裂来调节深部的构造应力。     

Thermorheological model for the European Central Alps: brittle–ductile transition and lithospheric strength

A two‐dimensional thermorheological model of the Central Alps along a north–south transect is presented. Thermophysical and rheological parameters of the various lithological units are chosen from seismic and gravity information. The inferred temperature distribution matches surface heat flow and results in Moho temperatures between 500 and 800 °C. Both European and Adriatic lithospheres have a ‘jelly‐sandwich’ structure, with a 15–20 km thick brittle upper crust overlying a ductile lower crust and a mantle lid whose uppermost part is brittle. The total strength of the lithosphere is of the order of 0.5–1.0 × 10¹³ N m⁻¹ if the upper mantle is dry, or slightly less if the upper mantle is wet. In both cases, the higher values correspond to the Adriatic indenter.     

Modelling the extension of heterogeneous hot lithosphere

The consequences of weak heterogeneities in the extension of soft and hot lithosphere without significant previous crustal thickening has been analysed in a series of centrifuge models. The experiments examined the effects of i) the location of heterogeneities in the ductile crust and/or in the lithospheric mantle, and ii) their orientation, perpendicular or oblique to the direction of bulk extension. The observed deformation patterns are all relevant to the so-called “wide rifting” mode of extension. Weak zones located in the ductile crust exert a more pronounced influence on localisation of deformation in the brittle layer than those located in the lithospheric mantle: the former localise faulting in the brittle crust whereas the latter tend to distribute faulting over a wider area. This latter behaviour depends in turn upon the decoupling provided by the ductile crust. Localised thinning in the brittle crust is accompanied by ductile doming of both crust and mantle. Domains of maximum thinning in the brittle crust and ductile crust and mantle are in opposition. Lateral differences in brittle crust thinning are accommodated by lateral flow in the ductile crust and mantle. This contrasts with “cold and strong” lithospheres whose high strength sub-Moho mantle triggers a necking instability at the lithosphere-scale. This also differs from the extension of thickened hot and soft lithospheres whose ductile crust is thick enough to give birth to metamorphic core complexes. Thus, for the given lithospheric rheology, the models have relevance to backarc type extensional systems, such as the Aegean and the Tyrrhenian domains.     

Topography of the crust–mantle boundary beneath the Black Sea Basin

A map of Moho depth for the Black Sea and its immediate surroundings has been inferred from 3-D gravity modelling, and crustal structure has been clarified. Beneath the basin centre, the thickness of the crystalline layer is similar to that of the oceanic crust. In the Western and Eastern Black Sea basins, the Moho shallows to 19 and 22 km, respectively. Below the Tuapse Trough (northeastern margin, adjacent to the Caucasus orogen), the base of the crust is at 28 km, whereas in the Sorokin Trough, it is as deep as 34 km. The base of the crust lies at 29 and 33 km depths respectively below the southern and northern parts of the Mid-Black Sea Ridge. For the Shatsky Ridge (between the Tuapse Trough and the Eastern Black Sea Basin), the Moho plunges from the northwest (33 km) to the southeast (40 km). The Arkhangelsky Ridge (south of the Eastern Black Sea Basin) is characterised by a Moho depth of 32 km. The crust beneath these ridges is of continental type.     

南华北盆地群岩石圈热-流变结构

结合南华北盆地群现代地温场资料和深部地震测深资料及岩石热物性参数,对南华北盆地群的热结构进行了研究。结果表明：南华北盆地群平均热流值为53.7 mW/m2,地幔热流为30～34 mW/m2,莫霍面温度为500～550℃,热岩石圈厚度为110～130 km。在此基础上,进行了岩石圈流变模拟,探讨了研究区的岩石圈流变特征及其地球动力学意义。南华北盆地群岩石圈强度为（7.6~23.3）×1012 N/m,具有显著的 “三明治”结构。上地壳表现为脆性变形,中、下地壳为韧性的流动变形。这一分层变形机制决定了南华北盆地群的成盆演化动力学过程。     

Upper lithospheric seismic velocity structure across the Pripyat Trough and the Ukrainian Shield along the EUROBRIDGE'97 profile

We present new results on the structure resulting from Palaeoproterozoic terrane accretion and later formation of one of the aulacogens in the East European Platform. Seismic data has been acquired along the 530-km-long, N–S-striking EUROBRIDGE'97 traverse across Sarmatia, a major crustal segment of the East European Craton. The profile extends across the Ukrainian Shield from the Devonian Pripyat Trough, across the Palaeoproterozoic Volyn Block and the Korosten Pluton, into the Archaean Podolian Block. Seismic waves from chemical explosions at 18 shot points at approximately 30-km intervals were recorded in two deployments by 120 mobile three-component seismographs at 3–4 km nominal station spacing. The data has been interpreted by use of two-dimensional tomographic travel time inversion and ray trace modelling. The high data quality allows modelling of the P- and S-wave velocity structure along the profile. There are pronounced differences in seismic velocity structure of the crust and uppermost mantle between the three main tectonic provinces traversed by the profile: (i) the Pripyat Trough is a ca. 4-km-deep sedimentary basin, fully located in the Osnitsk–Mikashevichi Igneous Belt in the northern part of the profile. The velocity structure is typical for a Precambrian craton, but is underlain by a ca. 5-km-thick lowest crustal layer of high velocity. The development of the Pripyat Trough appears to have only affected the upper crust without noticeable thinning of the whole crust; this may be explained by a rheologically strong lithosphere at the time of formation of the trough. (ii) Very high seismic velocity and Vp/Vs ratio characterise the Volyn Block and Korosten Pluton to a depth of 15 km and probably also the lowest crust. The values are consistent with an intrusive body of mafic composition in the upper crust that formed from bimodal melts derived from the mantle and the lower crust. (iii) The Podolian Block is close to a typical cratonic velocity structure, although it is characterised by relatively low seismic velocity and Vp/Vs ratio. A pronounced SW-dipping mantle reflector from Moho to at least 70 km depth may represent the Proterozoic suture between Sarmatia and Volgo–Uralia, the structure from terrane accretion, or a later shear zone in the upper mantle. The sub-Moho P-wave seismic velocity is high everywhere along the profile, with the exception of the area above the dipping reflector. This velocity change further supports a plate tectonic origin of the dipping mantle reflector. The profile demonstrates that structure from Palaeoproterozoic plate tectonic processes are still identifiable in the lithosphere, even where younger metamorphic equilibration of the crust has taken place.     

青藏高原板内地震震源深度分布规律及其成因

青藏高原板内地震以浅源地震为主, 下地壳基本上没有地震, 地震震源多集中在15~40 km的深度范围, 主要在中地壳内, 呈似层状弥散分布.其中30~33 km深度是一个优势层, 与壳内分层有关.总体上青藏高原南、北部的震源面略呈相向倾斜特征.70~100 km深度区间出现了比较集中的震级较小的地震, 可能与壳幔过渡带的拆离作用有关.高原内部的正断层系与板内地震密切相关, 是板内浅源地震的主控构造.总之, 青藏高原地震震源沿着活动的上地壳脆性层与软弱层之间的脆-韧性过渡带分布.这些板内地震活动属于大陆动力学过程, 与板块碰撞和板块俯冲无关.初步认为青藏高原浅层到深层多震层的成因分别是韧性基底与脆性盖层、韧性下地壳与脆性上地壳、韧性下地壳与脆性上地幔的韧-脆性转换、拆离和解耦的产物.      

Lithosphere structure of the Ukrainian Shield and Pripyat Trough in the region of EUROBRIDGE-97 (Ukraine and Belarus) from gravity modelling

The southern segment of the seismic profile EUROBRIDGE—EUROBRIDGE-97 (EB'97)—located in Belarus and Ukraine, crosses the suture zone between two main segments of the East European Craton—Fennoscandia and Sarmatia—as well as Sarmatia itself. At the initial stage of our study, a 3-D density model has been constructed for the crust of the study region, including the major part of the Osnitsa–Mikashevichi Igneous Belt (OMIB) superimposed by sediments of the Pripyat Trough (PT), and three domains in the Ukrainian Shield—the Volhyn Domain (VD) with the anorthosite–rapakivi Korosten Pluton (KP), the Podolian Domain (PD), and the Ros–Tikich Domain (RTD). The model comprises three layers—sediments with maximum thickness (6 km) in the PT and two heterogeneous layers in the crystalline crust separated at a depth of 15 km. 3-D calculations show the main features of the observed gravity field are caused by density heterogeneities in the upper crust. Allocation of density domains deeper than 15 km is influenced by Moho topography. Fitting the densities here reveals an increase (up to 2960 kg m⁻³) in the modelled bodies accompanied by a Moho deepening to 50 km. In contrast, a Moho uplift to a level of 35–37 km below the KP and major part of the PT is associated with domains of reduced densities. An important role for the deep Odessa–Gomel tectonic zone, dividing the crust into two regions one of basically Archean consolidation in the west (PD and RTD) and one of Proterozoic crust in the east (Kirovograd Domain)—was confirmed.2-D density modelling on the EB'97 profile shows that in the upper crust three main domains of different Precambrian evolution—the OMIB (with the superimposed PT), the VD with the KP, and the PD—can be distinguished. Deeper, in the middle and lower crust, layered structures having no connection to the surface geology are dominant features of the models. Least thickness of the crust was obtained below the KP. Greatest crustal thickness (more than 50 km) was found below the PD, characterised also by maximum deviation of velocity/density relation in the rocks from a standard one. The velocity and density models along the EB'97 profile have been interpreted together with inferred V_p/V_s ratios to estimate crustal composition in terms of SiO₂ content. In the course of the modelling, the status of the PD as a centre of Archean granulitic consolidation has been confirmed. The crustal structure of the anorthosite–rapakivi KP is complex. For the first time, a complicated structure for the lower crust and lower crust–upper mantle transition zone beneath the KP has been determined. The peculiarities of the crustal structure of the KP are quite well explained in terms of formation of rapakivi–anorthosite massifs as originating from melt chambers in the upper mantle and lower crust. An important role for the South Pripyat Fault (SPF), repeatedly activated during Proterozoic–Palaeozoic times, has been ascertained. At the subplatform stage of crustal evolution the SPF was, probably, a magma channel facilitating the granitic intrusions of the KP. In the Palaeozoic the fault was reactivated during rifting in the PT.     

Rock mechanics observations pertinent to the rheology of the continental lithosphere and the localization of strain along shear zones

Emphasized in this paper are the deformation processes and rheologies of rocks at high temperatures and high effective pressures, conditions that are presumably appropriate to the lower crust and upper mantle in continental collision zones. Much recent progress has been made in understanding the flexure of the oceanic lithosphere using rock-mechanics-based yield criteria for the inelastic deformations at the top and base. At mid-plate depths, stresses are likely to be supported elastically because bending strains and elastic stresses are low. The collisional tectonic regime, however, is far more complex because very large permanent strains are sustained at mid-plate depths and this requires us to include the broad transition between brittle and ductile flow. Moreover, important changes in the ductile flow mechanisms occur at the intermediate temperatures found at mid-plate depths.Two specific contributions of laboratory rock rheology research are considered in this paper. First, the high-temperature steady-state flow mechanisms and rheology of mafic and ultramafic rocks are reviewed with special emphasis on olivine and crystalline rocks. Rock strength decreases very markedly with increases in temperature and it is the onset of flow by high temperature ductile mechanisms that defines the base of the lithosphere. The thickness of the continental lithosphere can therefore be defined by the depth to a particular isotherm T_c above which (at geologic strain rates) the high-temperature ductile strength falls below some arbitrary strength isobar (e.g., 100 MPa). For olivine T_c is about 700°–800°C but for other crustal silicates, T_c may be as low as 400°–600°C, suggesting that substantial decoupling may take place within thick continental crust and that strength may increase with depth at the Moho, as suggested by a number of workers on independent grounds. Put another way, the Moho is a rheological discontinuity. A second class of laboratory observations pertains to the general phenomenon of ductile faulting in which ductile strains are localized into shear zones. Ductile faults have been produced in experiments of five different rock types and is generally expressed as strain softening in constant-strain-rate tests or as an accelerating-creep-rate stage at constant differential stress. A number of physical mechanisms have been identified that may be responsible for ductile faulting, including the onset of dynamic recrystallization, phase changes, hydrothermal alteration and hydrolytic weakening. Microscopic evidence for these processes as well as larger-scale geological and geophysical observations suggest that ductile faulting in the middle to lower crust and upper mantle may greatly influence the distribution and magnitudes of differential stresses and the style of deformation in the overlying upper continental lithosphere.     

The structure of the crust and upper mantle to depths of 320 km beneath the Kaapvaal craton, from P wave arrivals generated by regional earthquakes and mining-induced tremors

Travel times from earthquakes recorded at two seismic networks were used to derive an average P wavespeed model for the crust and upper mantle to depths of 320 km below southern Africa. The simplest model (BPI1) has a Moho depth of 34 km, and an uppermost mantle wavespeed of 8.04 km/s, below which the seismic wavespeeds have low positive gradients. Wavespeed gradients decrease slightly around 150 km depth to give a ‘knee’ in the wavespeed-depth model, and the wavespeed reaches 8.72 km/s at a depth of 320 km. Between the Moho and depths of 270 km, the seismic wavespeeds lie above those of reference model IASP91 of Kennett [Research School of Earth Sciences, Australian National University, Canberra, Australia (1991)] and below the southern African model of Zhao et al. [Journal of Geophysical Research 104 (1999) 4783]. At depths near 300 km all three models have similar wavespeeds. The mantle P wavespeeds for southern Africa of Qiu et al. [Geophysical Journal International 127 (1996) 563] lie close to BPI1 at depths between 40 and 140 km, but become lower at greater depths. The seismic wavespeeds in the upper mantle of model BPI1 agree satisfactorily with those estimated from peridotite xenoliths in kimberlites from within the Kaapvaal craton.The crustal thickness of 34 km of model BPI1 is systematically lower than the average thickness of 41 km computed over the same region from receiver functions. This discrepancy can be partly explained by an alternative model (BPI2) in which there is a crust–mantle transition zone between depths of 35 and 47 km, below which seismic wavespeed increases to 8.23 km/s. A low-wavespeed layer is then required at depths between 65 and 125 km.     

中国东北地区地热资源及热结构分析

在我国东北地区，盆地的热流高而在山区和地台区（额尔古纳和佳木斯地台）热流较低，通过对该区的研究发现，在地壳和上地幔中莫霍面埋深，高导层埋深和热流值之间存在着密切的联系，通过对该区热流数据分析我们可以得到这样的结论：该区的地壳和上地幔的热结构引起了热流分布的变化。计算结果还表明，该区地台的地幔热流，地壳中10km以下和10km以上的热流对地表热流的贡献不同。该区松辽盆地中广泛分布着传导型中低温地热资源，对松辽盆地地热资源的开发利用有着广泛的前景。最后结合已有的地质和物探资料，还给出了长白山天池火山地区长白温泉地热系统的概念性模型。     

Seismic velocity model of the crust and uppermost mantle around the Mirnyi kimberlite field in Siberia

We present the first detailed seismic velocity models of the crust and uppermost mantle around the Mirnyi kimberlite field in Yakutia, Siberia. We have digitized vintage seismograms that were acquired in 1981 and 1983 by use of Taiga analogue seismographs along two perpendicular seismic profiles. The 370-km long, northwest striking profile I across the kimberlite pipe was covered by 41 seismographs, which recorded seismic signals from 21 chemical shots along the line, including one off-end shot. The perpendicular, 340-km long profile II across profile I ca. 30 km to the south of the Mirnyi kimberlite field was covered by 45 seismographs, which recorded seismic signals from 22 chemical shots, including four off-end shots. Each shot involved detonation of between 1.5 and 6.0 tons of TNT, distributed in individual charges of 100–200 kg in shallow water (< 2 m deep). The data is of high quality with high signal/noise ratio to the farthest offsets. We present the results from two-dimensional ray tracing, forward modelling.Both velocity models show normal cratonic structure of the ca. 45-km-thick crust with only slight undulation of the Moho. However, relatively small seismic velocity is detected to 25-km depth in a ca. 60-km wide zone around the kimberlite pipe, surrounded by elevated velocity (> 6.3 km/s) in the upper crust. The lower crust has a relatively constant velocity of 6.8–6.9 km/s. It appears relatively unaffected by the presence of the kimberlite field. Extremely large P-wave velocity (> 8.7 km/s) of the sub-Moho mantle is interpreted along profile I, except for a 70-km wide zone with a “normal” Pn velocity of 8.1 km/s below the kimberlite. Profile II mainly shows Pn velocities of 8.0–8.2 km/s, with unusually large velocity (> 8.5 km/s) in two, ca. 100-km wide zones, at its southwestern end, one zone being close to the kimberlite field. The nature of these exceptionally large, sub-Moho mantle velocities is not yet understood. The difference in velocity in the two profile directions indicates anisotropy, but the effect of unusual rock composition, e.g. from a high concentration of garnet, cannot be excluded.     

Crust–upper mantle seismic velocity structure across Southeastern China

One in-line wide-angle seismic profile was conducted in 1990 in the course of the Southeastern China Continental Dynamics project aimed at the study of the contact between the Cathaysia block and the Yangtze block. This 380-km-long profile extended in NW–SE direction from Tunxi, Anhui Province, to Wenzhou, Zhejiang Province. Five in-line shots were fired and recorded at seismic stations with spacing of about 3 km along the recording line. We have used two-dimensional ray tracing to model P- and S-wave arrivals and provide constraints on the velocity structure of the upper crust, middle crust, lower crust, Moho discontinuity, and the top part of the lithospheric mantle. P-wave velocity, S-wave velocity and V_P/V_S ratio are mapped. The crust is 36-km thick on average, albeit it gradually thins from the northwest end to the southeast end (offshore) of the profile. The average crustal velocity is 6.26 km/s for P-waves but 3.6 km/s for S-waves. A relatively narrow low-velocity layer of about 4 km of thickness, with P- and S-wave velocities of 6.2 km/s and 3.5 km/s, respectively, marks the bottom of the middle crust at a depth of 23-km northwest and 17-km southeast. At the crust–mantle transition, the P- and S-wave velocity change quickly from 7.4 to 7.8 km/s (northwest) and 8.0 to 8.2 km/s (southeast) and from 3.9 to 4.2 km/s (northwest) and 3.9 to 4.5 km/s (southeast), respectively. This result implies a lateral contrast in the upper mantle velocity along the 140 km sampled by the profile approximately. The average V_P/V_S ratio ranges from 1.68–1.8 for the upper crust to 1.75 for the middle and 1.75–1.85 for lower crust. With the interpretation of the wide-angle seismic data, Jiangshan–Shaoxin fault is considered as the boundary between the Yangtze and the Cathaysia block.     

Lower crust indentation or horizontal ductile flow during continental collision?

Conditions for indentation and channelised flow are investigated with two-dimensional thermomechanical models of Alpine-type continental collision. The models mimic the development of an orogen at an initial central portion of weakened lithosphere 150 km wide, coherent with several geological reconstructions. We study in particular the role of lower crustal strength in developing peculiar geometries after 20 Ma of shortening at 1 cm/year. Crustal layers produce geometries of imbricate layers, which result from two contrasted mechanisms of either channelised ductile lateral flow or horizontal rigid-like indentation:

– Channelised lateral flow develops when the lateral lower crust has a viscosity less than 10²¹ Pa s, exhibiting velocities opposite to the direction of convergence. This mechanism of deformation produces subhorizontal shear zones at the boundaries between the lower crust and the more competent upper crust and lithospheric mantle. It is also associated with a topographic plateau that equilibrates with a wide (about 200 km) but quasi-constant crustal root about 50 km deep.

– In contrast, indentation occurs with lateral lower crust layers that have a viscosity greater than about 10²³ Pa s, producing significant shortening and thickening of the central crust. In this case topography develops steep and narrow (around 100 km wide), associated with a thickened crust exceeding 60 km depth. A crustal-scale pop-up forms bounded by subvertical shear zones that root into the mantle lithosphere.

Crustal structure at the easternmost termination of the Variscan belt based on CELEBRATION 2000 and ALP 2002 data

The eastern margin of the Variscan belt in Europe comprises plate boundaries between continental blocks and terranes formed during different tectonic events. The crustal structure of that complicated area was studied using the data of the international refraction experiments CELEBRATION 2000 and ALP 2002. The seismic data were acquired along SW–NE oriented refraction and wide-angle reflection profiles CEL10 and ALP04 starting in the Eastern Alps, passing through the Moravo-Silesian zone of the Bohemian Massif and the Fore-Sudetic Monocline, and terminating in the TESZ in Poland. The data were interpreted by seismic tomographic inversion and by 2-D trial-and-error forward modelling of the P waves. Velocity models determine different types of the crust–mantle transition, reflecting variable crustal thickness and delimiting contacts of tectonic units in depth. In the Alpine area, few km thick LVZ with the Vp of 5.1 km s^− 1 dipping to the SW and outcropping at the surface represents the Molasse and Helvetic Flysch sediments overthrust by the Northern Calcareous Alps with higher velocities. In the Bohemian Massif, lower velocities in the range of 5.0–5.6 km s^− 1 down to a depth of 5 km might represent the SE termination of the Elbe Fault Zone. The Fore-Sudetic Monocline and the TESZ are covered by sediments with the velocities in the range of 3.6–5.5 km s^− 1 to the maximum depth of 15 km beneath the Mid-Polish Trough. The Moho in the Eastern Alps is dipping to the SW reaching the depth of 43–45 km. The lower crust at the eastern margin of the Bohemian Massif is characterized by elevated velocities and high Vp gradient, which seems to be a characteristic feature of the Moravo-Silesian. Slightly different properties in the Moravian and Silesian units might be attributed to varying distances of the profile from the Moldanubian Thrust front as well as a different type of contact of the Brunia with the Moldanubian and its northern root sector. The Moho beneath the Fore-Sudetic Monocline is the most pronounced and is interpreted as the first-order discontinuity at a depth of 30 km.     

Crustal density structure in the Spanish Central System derived from gravity data analysis (Central Spain)

Shallow and deep sources generate a gravity low in the central Iberian Peninsula. Long-wavelength shallow sources are two continental sedimentary basins, the Duero and the Tajo Basins, separated by a narrow mountainous chain called the Spanish Central System. To investigate the crustal density structure, a multitaper spectral analysis of gravity data was applied. To minimise biases due to misleading shallow and deep anomaly sources of similar wavelength, first an estimation of gravity anomaly due to Cenozoic sedimentary infill was made. Power spectral analysis indicates two crustal discontinuities at mean depths of 31.1 ± 3.6 and 11.6 ± 0.2 km, respectively. Comparisons with seismic data reveal that the shallow density discontinuity is related to the upper crust lower limit and the deeper source corresponds to the Moho discontinuity. A 3D-depth model for the Moho was obtained by inverse modelling of regional gravity anomalies in the Fourier domain. The Moho depth varies between a mean depth of 31 km and 34 km. Maximum depth is located in a NW–SE trough. Gravity modelling points to lateral density variations in the upper crust. The Central System structure is described as a crustal block uplifted by NE–SW reverse faults. The formation of the system involves displacement along an intracrustal detachment in the middle crust. This detachment would split into several high-angle reverse faults verging both NW and SE. The direction of transport is northwards, the detachment probably being rooted at the Moho.     

Thermal and mechanical structure of the central Iberian Peninsula lithosphere

The central Iberian Peninsula (Spain) is made up of three main tectonic units: a mountain range, the Spanish Central System and two Tertiary basins (those of the rivers Duero and Tajo). These units are the result of widespread foreland deformation of the Iberian plate interior in response to Alpine convergence of European and African plates. The present study was designed to investigate thermal structure and rheological stratification in this region of central Spain. Surface heat flow has been described to range from 80 to 60 mW m⁻². Highest surface heat flow values correspond to the Central System and northern part of the Tajo Basin. The relationship between elevation and thermal state was used to construct a one-dimensional thermal model. Mantle heat flow drops from 34 mW m⁻² (Duero Basin) to 27 mW m⁻² (Tajo Basin), and increases with diminishing surface heat flow. Strength predictions made by extrapolating experimental data indicate varying rheological stratification throughout the area. In general, in compression, ductile fields predominate in the middle and lower crusts and lithospheric mantle. Brittle behaviour is restricted to the first 8 km of the upper crust and to a thin layer at the top of the middle crust. In tension, brittle layers are slightly more extended, while the lower crust and lithospheric mantle remain ductile in the case of a wet peridotite composition. Discontinuities in brittle and ductile layer thickness determine lateral rheological anisotropy. Tectonic units roughly correspond to rheological domains. Brittle layers reach their maximum thickness beneath the Duero Basin and are of least thickness under the Tajo Basin, especially its northern area. Estimated total lithospheric strength shows a range from 2.5×10¹² to 8×10¹² N m⁻¹ in compression, and from 1.3×10¹² to 1.6×10¹² N m⁻¹ in tension. Highest values were estimated for the Duero Basin.Depth versus frequency of earthquakes correlates well with strength predictions. Earthquake foci concentrate mainly in the upper crust, showing a peak close to maximum strength depth. Most earthquakes occur in the southern margin of the Central System and southeast Tajo Basin. Seismicity is related to major faults, some bounding rheological domains. The Duero Basin is a relative quiescence zone characterised by higher total lithospheric strength than the remaining units.     

Crustal and upper mantle seismic structure of the Australian Plate, South Island, New Zealand

Seismic reflection and refraction data were collected west of New Zealand's South Island parallel to the Pacific–Australian Plate boundary. The obliquely convergent plate boundary is marked at the surface by the Alpine Fault, which juxtaposes continental crust of each plate. The data are used to study the crustal and uppermost mantle structure and provide a link between other seismic transects which cross the plate boundary. Arrival times of wide-angle reflected and refracted events from 13 recording stations are used to construct a 380-km long crustal velocity model. The model shows that, beneath a 2–4-km thick sedimentary veneer, the crust consists of two layers. The upper layer velocities increase from 5.4–5.9 km/s at the top of the layer to 6.3 km/s at the base of the layer. The base of the layer is mainly about 20 km deep but deepens to 25 km at its southern end. The lower layer velocities range from 6.3 to 7.1 km/s, and are commonly around 6.5 km/s at the top of the layer and 6.7 km/s at the base. Beneath the lower layer, the model has velocities of 8.2–8.5 km/s, typical of mantle material. The Mohorovicic discontinuity (Moho) therefore lies at the base of the second layer. It is at a depth of around 30 km but shallows over the south–central third of the profile to about 26 km, possibly associated with a southwest dipping detachment fault. The high, variable sub-Moho velocities of 8.2 km/s to 8.5 km/s are inferred to result from strong upper mantle anisotropy. Multichannel seismic reflection data cover about 220 km of the southern part of the modelled section. Beneath the well-layered Oligocene to recent sedimentary section, the crustal section is broadly divided into two zones, which correspond to the two layers of the velocity model. The upper layer (down to about 7–9 s two-way travel time) has few reflections. The lower layer (down to about 11 s two-way time) contains many strong, subparallel reflections. The base of this reflective zone is the Moho. Bi-vergent dipping reflective zones within this lower crustal layer are interpreted as interwedging structures common in areas of crustal shortening. These structures and the strong northeast dipping reflections beneath the Moho towards the north end of the (MCS) line are interpreted to be caused by Paleozoic north-dipping subduction and terrane collision at the margin of Gondwana. Deeper mantle reflections with variable dip are observed on the wide-angle gathers. Travel-time modelling of these events by ray-tracing through the established velocity model indicates depths of 50–110 km for these events. They show little coherence in dip and may be caused side-swipe from the adjacent crustal root under the Southern Alps or from the upper mantle density anomalies inferred from teleseismic data under the crustal root.     

Anomalous lithospheric structure of Northern Juan de Fuca plate — a consequence of oceanic rift propagation?

Anomalous crustal and upper mantle structure of northern Juan de Fuca plate is revealed from wide-angle seismic and gravity modelling. A 2-D velocity model is produced for refraction line II of the 1980 Vancouver Island Seismic Project (VISP80). The refraction data were recorded on three ocean bottom seismometers (OBSs) deployed at the ends and middle of a 110 km line oriented parallel to the North American continental margin. The velocity model is constructed via ray tracing and conforms to first-arrival amplitude observations and travel time picks of direct, converted and reflected phases. Between sub-sediment depths of 3 to 11 km, depths normally associated with the lower crust and upper oceanic mantle, the final model shows that compressional-wave velocities decrease significantly from southeast to northwest along the profile. At sub-sediment depths of 11 km at the northwestern end of the profile, P-wave velocities are as low as 7.2 km/s. A complementary 2-D gravity model using the geometry of the velocity model and velocity–density relationships characteristic of oceanic crust is produced. The high densities required to match the gravity field indicate the presence of peridotites containing 25–30% serpentine by volume, rather than excess gabbroic crust, within the deep low velocity zone. Anomalous travel time delays and unusual reflection characteristics observed from proximal seismic refraction and reflection experiments suggest a broader zone of partially serpentinized peridotites coincident with the trace of a pseudofault. We propose that partial serpentinization of the upper mantle is a consequence of slow spreading at the tip of a propagating rift.