Structure of crust and upper mantle beneath the Ordos Block and the Yinshan Mountains revealed by receiver function analysis

A temporary seismological network of broadband three-component stations has been deployed N–S to investigate the crust and upper mantle structure across the Ordos Block and the Yinshan Mountains. P wave receiver functions reveal the Moho depth to be about 41 km beneath the central Ordos Block and down to 45 km beneath the northern Ordos Block, a slight uplifting to 42–43 km beneath the Hetao Graben, increasing to 47–48 km beneath the Yinshan Mountains and then decreasing to 44 km beneath the northern Yinshan Mountains along the profile. In the Ordos Block, the crustal Vp/Vs ratio (about 1.80) south to the Hetao Graben differs from that (about 1.75) beneath the center Ordos Block. The crustal Vp/Vs ratio is significantly lower (about 1.65–1.70) beneath the Yinshan Mountains. The P wave receiver function migration imaging suggests relatively flat discontinuities at 410 and 660 km, indicating the lack of a strong thermal anomaly beneath this profile at these depths, and a low S wave velocity anomaly in the upper mantle beneath the Hetao Graben. We suggest that the low S wave velocity anomaly may be attributable to heat and that the thermal softening advances the evolution of the Hetao Graben, while the lower-crustal ductile flows transfer from the Hetao Graben to the northern Ordos Block, resulting in crustal thickening.     

The distribution of earthquake multiplets beneath the southwest Pacific

Earthquakes beneath the southwest Pacific occur from the surface down to 700 km depth. Teleseismic waveforms created by some of these earthquakes are almost identical. We investigate Tonga–Kermadec and Vanuatu subduction zone earthquake P-coda waveforms using a cross-correlation technique and hierarchical clustering algorithm in order to determine the origin of waveform similarity and the distribution of earthquakes producing similar waveforms.We show that scatterers forming the majority of power in the P-wave coda are localised around the receiver. As a result, waveform similarity provides a much weaker constraint on source separation than in local studies. Waveform similarity can provide stronger constraints on focal mechanism.Most earthquake multiplets within the Tonga–Fiji–Kermadec Wadati–Benioff zone are found at depths between 0–60 km and 520–620 km. A significant proportion of all deep-focus events in south Pacific subduction zones have waveforms similar to those of at least one other event. Relative relocation of events within the largest identified multiplet reveals a planar zone of seismicity sub-parallel to the nodal plane of a related centroid moment tensor solution.Groups of earthquakes with similar waveforms remain active on at least the 14-year recording timescale. We equate this to repeated rupture on single or closely related shear systems within the subducting slabs.     

On the Vp/Vs–Mg# correlation in mantle peridotites: Implications for the identification of thermal and compositional anomalies in the upper mantle

We use thermodynamically self-consistent and hybrid methods to analyze the correlation of important physical parameters (e.g. bulk density, elastic moduli) with bulk Mg# and modal composition in mantle peridotites at upper mantle conditions. Temperature (anharmonic and anelastic), pressure and compositional derivatives for all these parameters are evaluated. The results show that the widely used correlations between Vp/Vs and Mg# in peridotites are strictly valid only for garnet-bearing assemblages at temperatures < 900 °C. The correlation breaks down when: i) spinel is the stable Al-rich phase in the assemblage and ii) when anelastic attenuation of seismic velocities becomes important (T ? 900 °C). This implies that the range of applicability of published Vp/Vs–Mg# correlations for the upper mantle is limited to a depth interval between the spinel–garnet phase transition and the 900 °C isotherm. We use numerical simulations to show that this depth interval is virtually nonexistent in lithospheres thinner than ～ 140 km and can comprise up to ～ 50% of the lithospheric mantle in thick (> 220 km) lithospheric domains. In addition, we show that for most of the upper mantle the expected Δ(Vp/Vs) values associated with compositional variations are smaller than the resolution limit of current seismological methods. All these considerations suggest that the Vp/Vs ratio is not a reliable measure of compositional variations and that for large parts of the upper mantle compositional anomalies cannot be separated from thermal anomalies on the basis of seismological studies only. We further confirm that the only reliable indicator of compositional anomalies in a peridotitic mantle is the ratio of density to shear wave velocities (ρ/Vs). Our results demonstrate that geophysical–petrological models (forward or inverse) that model these two fields (i.e. density and Vs) self-consistently within a robust thermodynamic framework are necessary for characterizing the small-scale thermal and compositional structure of the lithosphere and sublithospheric upper mantle.     

Structural heterogeneity of an ultra-low-velocity zone beneath the Philippine Islands: Implications for core–mantle chemical interactions induced by massive partial melting at the bottom of the mantle

An ScP phase reflected and converted at the core–mantle boundary (CMB) beneath the region east of the Philippine Islands shows clear pre- and postcursors, recorded on short-period seismic networks in Japan. These waveform variations can be explained by interaction of the ScP wavefield with thin layers at the CMB. The results of forward modeling of double-array stacks reveal two different structural heterogeneities in the lowermost mantle beneath the region east of the Philippine Islands. One of the structures represents a decreased velocity, and increased density across the reflector at the lowermost ～10 km of the mantle, with P- and S-wave velocity reductions of 5–10% and ～30%, respectively, and an increase in density of 5–10%. Another structure consists of a pair of reflectors at ～10 km and ～5 km above the CMB, both of which are characterized by reduced P- and S-wave velocities. The upper reflector is the interface of a low-velocity zone in which P- and S-wave velocities decrease of 10% and 30%, respectively, accompanied by an extremely large increase in density (20–25%). The lower reflector is characterized by a 25% reduction in S-wave velocity relative to the above low-velocity layer, as well as a 5% decrease in P-wave velocity and no change in density. The nature of the low-velocity zone detected locally at the CMB is comparable with that of ultra-low-velocity zones (ULVZs) observed by various seismic probes in the South Pacific and Central America. Extensive observations of the ULVZ beneath the region east of the Philippine Islands indicate massive partial melting at the bottom of the mantle. Low-S-velocity basal layer partly detected within the ULVZ may be resulting from core–mantle chemical interactions, driven by massive partial melting.     

Crustal velocity structure in the southern edge of the Central Alborz (Iran)

Inversion of local earthquake travel times and joint inversion of receiver functions and Rayleigh wave group velocity measurements were used to derive a simple model for the velocity crustal structure beneath the southern edge of the Central Alborz (Iran), including the seismically active area around the megacity of Tehran. The P and S travel times from 115 well-located earthquakes recorded by a dense local seismic network, operated from June to November 2006, were inverted to determine a 1D velocity model of the upper crust. The limited range of earthquake depths (between 2 km and 26 km) prevents us determining any velocity interfaces deeper than 25 km. The velocity of the lower crust and the depth of the Moho were found by joint inversion of receiver functions and Rayleigh wave group velocity data. The resulting P-wave velocity model comprises an upper crust with 3 km and 4 km thick sedimentary layers with P wave velocities (V_p) of ～5.4 and ～5.8 km s^?1, respectively, above 9 km and 8 km thick layers of upper crystalline crust (V_p ～6.1 and ～6.25 km s^?1 respectively). The lower crystalline crust is ～34 km thick (V_p ～ 6.40 km s^?1). The total crustal thickness beneath this part of the Central Alborz is 58 ± 2 km.     

Crustal structure variation along 30°N in the eastern Tibetan Plateau and its tectonic implications

We determined crustal structure along the latitude 30°N through the eastern Tibetan Plateau using a teleseismic receiver function analysis. The data came mostly from seismic stations deployed in eastern Tibet and western Sichuan region from 2004 to 2006. Crustal thickness and Vp/Vs ratio at each station were estimated by the H–k stacking method. On the profile, the mean crustal thickness and Vp/Vs ratio were found to be 62.3 km and 1.74 in the Lhasa block, 71.2 km and 1.79 near the Bangong–Nujiang suture, 66.3 km and 1.80 in the Qiangtang block, 59.8 km and 1.81 in the Songpan–Garze block, and 42.9 km and 1.76 in the Yangtze block, respectively. The estimated crustal thicknesses are consistent with predictions based on the topography and the Airy isostasy, except near the Bangong–Nujiang suture and in the Qiangtang block where the crust is 5–10 km thicker than predicted, indicating that the crust may be denser, possibly due to mafic underplating. We also inverted receiver functions for crustal velocity structure along the profile, which reveals a low S-wave velocity zone in the lower crust beneath the eastern Tibetan Plateau, although the extent of the low-velocity zone varies considerably. The low-velocity zone, together with previous results, suggests limited partial melting and localized crustal flow in the lower crust of the eastern Tibetan Plateau.     

Crustal seismic structure and depth distribution of earthquakes in the Archean Kuusamo region,Fennoscandian Shield

Two-dimensional crustal velocity models are derived from passive seismic observations for the Archean Karelian bedrock of north-eastern Finland. In addition, an updated Moho depth map is constructed by integrating the results of this study with previous data sets. The structural models image a typical three-layer Archean crust, with thickness varying between 40 and 52 km. P wave velocities within the 12–20 km thick upper crust range from 6.1 to 6.4 km/s. The relatively high velocities are related to layered mafic intrusive and volcanic rocks. The middle crust is a fairly homogeneous layer associated with velocities of 6.5–6.8 km/s. The boundary between middle and lower crust is located at depths between 28 and 38 km. The thickness of the lower crust increases from 5–15 km in the Archean part to 15–22 km in the Archean–Proterozoic transition zone. In the lower crust and uppermost mantle, P wave velocities vary between 6.9–7.3 km/s and 7.9–8.2 km/s. The average V_p/V_s ratio increases from 1.71 in the upper crust to 1.76 in the lower crust.The crust attains its maximum thickness in the south-east, where the Archean crust is both over- and underthrust by the Proterozoic crust. A crustal depression bulging out from that zone to the N–NE towards Kuusamo is linked to a collision between major Archean blocks. Further north, crustal thickening under the Salla and Kittilä greenstone belts is tentatively associated with a NW–SE-oriented collision zone or major shear zone. Elevated Moho beneath the Pudasjärvi block is primarily explained with rift-related extension and crustal thinning at ∼2.4–2.1 Ga.The new crustal velocity models and synthetic waveform modelling are used to outline the thickness of the seismogenic layer beneath the temporary Kuusamo seismic network. Lack of seismic activity within the mafic high-velocity body in the uppermost 8 km of crust and relative abundance of mid-crustal, i.e., 14–30 km deep earthquakes are characteristic features of the Kuusamo seismicity. The upper limit of seismicity is attributed to the excess of strong mafic material in the uppermost crust. Comparison with the rheological profiles of the lithosphere, calculated at nearby locations, indicates that the base of the seismogenic layer correlates best with the onset of brittle to ductile transition at about 30 km depth.We found no evidence on microearthquake activity in the lower crust beneath the Archean Karelian craton. However, a data set of relatively well-constrained events extracted from the regional earthquake catalogue implies a deeper cut-off depth for earthquakes in the Norrbotten tectonic province of northern Sweden.     

Density modeling of the Escollos Alijos Seamount from inversion of its geoid undulation anomaly

Escollos Alijos is a large seamount located in the NE Pacific Ocean about 300 km off the Baja California Peninsula. Geochronology and geochemical analysis of volcanic rocks capping the seamount indicate recent magmatism that resulted from extensive differentiation of a mildly alkalic basalt parent magma.Escollos Alijos is located towards the eastern edge of a long-wavelength geoid undulation minimum, of up to -47 m with respect of the WGS84 ellipsoid, which extends over the northeastern Pacific Ocean. Subtracting from the geoid undulation its long-wavelength component and the undulation due to the seamount topography itself, a negative undulation anomaly persists that indicates a mass deficit at depth. Linear inversion of the undulation anomaly yields a region characterized by a negative density contrast, localized under the seamount at a depth between 9 and 13 km.The age and chemical composition of Escollos Alijos, and the inferred mass deficit suggest magma trapped between the oceanic crust and the uppermost mantle, which explains the magmatic activity in recent times.     

Crustal intrusion beneath the Louisville hotspot track

We report here the first detailed 2D tomographic image of the crust and upper mantle structure of a Cretaceous seamount that formed during the interaction of the Pacific plate and the Louisville hotspot. Results show that at ～ 1.5 km beneath the seamount summit, the core of the volcanic edifice appears to be dominantly intrusive, with velocities faster than 6.5 km/s. The edifice overlies both high lower crustal (> 7.2–7.6 km/s) and upper mantle (> 8.3 km/s) velocities, suggesting that ultramafic rocks have been intruded as sills rather than underplated beneath the crust. The results suggest that the ratio between the volume of intra-crustal magmatic intrusion and extrusive volcanism is as high as ～ 4.5. In addition, the inversion of Moho reflections shows that the Pacific oceanic crust has been flexed downward by up to ～ 2.5 km beneath the seamount. The flexure can be explained by an elastic plate model in which the seamount emplaced upon oceanic lithosphere that was ～ 10 Myr at the time of loading. Intra-crustal magmatic intrusion may be a feature of hotspot volcanism at young, hot, oceanic lithosphere, whereas, magmatic underplating below a pre-existing Moho may be more likely to occur where a hotspot interacts with oceanic lithosphere that is several tens of millions of years old.     

Time-Domain Moment Tensors for shallow (h ≤ 40 km) earthquakes in the broader Aegean Sea for the years 2006 and 2007: The database of the Aristotle University of Thessaloniki

We present a catalog of moment tensor (MT) solutions and moment magnitudes, M_w, for 119 shallow (h ≤ 40 km) earthquakes in Greece and its surrounding lands (34°N–42°N, 19°E–30°E) for the years 2006 and 2007, computed with the 1D Time-Domain Moment Tensor inversion method (TDMT_INV code of Dreger, 2003). Magnitudes range from 3.2 ≤ M_w ≤ 5.7. Green's functions (GF) have been pre-computed to build a library, for a number of velocity profiles applicable to the broader Aegean Sea region, to be used in the inversion of observed broad band waveforms (10–50 s). All MT solutions are the outcome of a long series of tests of different reported source locations and hypocenter depths. Quality factors have been assigned to each MT solution based on the number of stations used in the inversion and the goodness of fit between observed and synthetic waveforms. In general, the focal mechanisms are compatible with previous knowledge on the seismotectonics of the Aegean area. The new data provide evidence for strike-slip faulting along NW–SE trending structures at the lower part of Axios basin, close to the heavily industrialized, and presently subsiding, region of the city of Thessaloniki. Normal faulting along E–W trending planes is observed at the Strimon basin, and in Orfanou Gulf in northern Greece. A sequence of events in the east Aegean Sea close to the coastline with western Anatolia sheds light on an active structure bounding the north coastline of Psara–Chios Islands about 20–25 km in length exhibiting right lateral strike-slip faulting.     

Imaging Pn velocities beneath the Pannonian basin

The Pannonian depression is an extensional back-arc basin in central Europe and is an integral part of the Alpine–Carpathian orogenic mountain belts. It can be characterized by thinned lower crust, shallow Moho discontinuity, high surface heat flow and Moho temperature, implying recent active tectonic processes. Imaging the velocity structure of the upper mantle may help us to better understand the structure and formation of the Pannonian region.In this paper, P_n traveltimes from regional earthquakes are used to tomographically image the lateral velocity variations in the uppermost mantle beneath the Pannonian basin. The set of linear tomographic equations, built up of the time term equation for each source–receiver pair, is solved by a truncated singular value decomposition algorithm. The explicit computation of the generalized inverse of the tomographic equations makes it possible to deduce both the resolution matrix and the model covariance matrix, allowing us to estimate the resolution and reliability of the solution.The mean compressional wave velocity in the uppermost mantle beneath the Pannonian basin is 7.9 km/s, substantially lower than the average continental P_n velocity of 8.1 km/s. It is mostly due to the high Moho temperature having values on average 400–500 °C more than those in the surrounding areas. The velocity anomalies range from −0.3 to 0.3 km/s relative to the mean velocity of 7.9 km/s. Due to high Moho temperature, below the North Hungarian range low (7.6–7.7 km/s) velocities can be found. High-velocity anomalies of around 8.1 km/s can be detected along the W-SW boundaries of Hungary and at the junction of the Pannonian basin and the Southern Carpathians. The Great Hungarian Plain shows average (7.9 km/s) P_n velocities.     

Estimates of eruption velocity and plume height from infrasonic recordings of the 2006 eruption of Augustine Volcano,Alaska

The 2006 eruption of Augustine Volcano, Alaska, began with an explosive phase comprising 13 discrete Vulcanian blasts. These events generated ash plumes reaching heights of 3–14 km. The eruption was recorded by a dense geophysical network including a pressure sensor located 3.2 km from the vent. Infrasonic signals recorded in association with the eruptions have maximum pressures ranging from 13–111 Pa. Eruption durations are estimated to range from 55–350 s. Neither of these parameters, however, correlates with eruption plume height. The pressure record, however, can be used to estimate the velocity and flux of material erupting from the vent, assuming that the sound is generated as a dipole source. Eruptive flux, in turn, is used to estimate plume height, assuming that the plume rises as a buoyant thermal. Plume heights estimated in this way correlate well with observations. Events that exhibit strongly impulsive waveforms are underestimated by the model, suggesting that flow may have been supersonic.     

Estimates of the transition zone temperature in a mechanically mixed upper mantle

Partial melting of mantle peridotite generates a physically and chemically layered oceanic lithosphere that is cycled back into the mantle in subduction zones. Stirring times of the mantle are too long to allow for complete re-homogenization of subducted basalt and harzburgite, given the low chemical diffusivity of the solid mantle. This suggests that the Earth's mantle is a mechanical mixture of basaltic and harzburgitic components. Using a recently developed thermodynamic formulism we determine the phase equilibria and the seismic properties of a mantle comprised of a mechanical mixture of basalt and harzburgite (MM) and a homogeneous mantle (EA) with identical pyrolitic bulk chemistry. We use the theoretical shear velocity profiles as a new thermometer of the mantle below the magma-genetic zone by modeling the difference ΔT_410-660 between traveltimes of shear wave reflections off the 410-km and 660-km with the potential temperature T_P. ΔT_410-660 are measured from waveform stacks. They indicate that, over 1000+ km wave lengths, the temperature varies by about 200 K. Lowest and highest temperatures are resolved for the western Pacific subduction zones and the central Pacific, respectively. This variation is similar for the EA and MM and is in excellent agreement with estimates of transition zone thickness and shear velocity variations. The median value of T_P for the EA is 1720 K. It is about 1625 K for the MM, a value that is in better agreement with the Normal-MORB values of 1610 ± 40 K inferred from olivine-liquid equilibria given that our sampling region encompasses the Western Pacific subduction zones and the oldest parts of the Pacific Plate. We argue therefore that a mechanical mixed mantle, with generally higher velocities and steeper velocities gradients, represents a better physical reference model than a model based on a fully equilibrated assemblage.     

Temperature dependence of the elastic moduli of ringwoodite

The elastic moduli of polycrystalline ringwoodite, (Mg_0.91Fe_0.09)₂SiO₄, were measured up to 470 K by means of the resonant sphere technique. The adiabatic bulk (K_S) and shear (μ) moduli were found to be 185.1(2) and 118.22(6) GPa at room temperature, and the average slopes of dK_S/dT and dμ/dT in the temperature range of the study were determined to be −0.0193(9) and −0.0148(3) GPa/K, respectively. Using these results, we estimate seismic wave velocity jumps for a pure olivine mantle model at 520 km depth. We find that the jump for the S-wave velocity is about 1.5 times larger than that for the P-wave velocity at this depth. This suggests that velocity jumps at the 520 km discontinuity are easier to detect using S-waves than P-waves.     

Observations of narrow bipolar events in East China

Narrow bipolar events (NBEs) are a distinct class of intra-cloud lightning discharge. In this paper we present observations of 10 negative and 67 positive such events in East China. Positive NBEs occurred at 7–12 km altitude above mean sea level (MSL) with a mean altitude of 9.5 km, and negative NBEs occurred at 14–16 km altitude. Electrical/channel characteristics of these events were derived from NBE pulse waveforms based on the transmission-line model. On average, the peak current moment and the charge moment change of a NBE event is 15 kA km, and 0.12 C km, respectively. The mean time for the propagation of current front along the channel is 2.2 μs. The upper limit on channel length for NBEs in this study is 510–1060 m, the lower limit on discharge current amplitude is 12.5–43.2 kA, and the minimum charge transfer is 0.1–0.3 C.     

Depth constraints on azimuthal anisotropy in the Great Basin from Rayleigh-wave phase velocity maps

We present fundamental-mode Rayleigh-wave azimuthally anisotropic phase velocity maps obtained for the Great Basin region at periods between 16 s and 102 s. These maps offer the first depth constraints on the origin of the semi-circular shear-wave splitting pattern observed in central Nevada, around a weak azimuthal anisotropy zone. A variety of explanations have been proposed to explain this signal, including an upwelling, toroidal mantle flow around a slab, lithospheric drip, and a megadetachment, but no consensus has been reached. Our phase velocity study helps constrain the three-dimensional anisotropic structure of the upper mantle in this region and contributes to a better understanding of the deformation mechanisms taking place beneath the western United States. The dispersion measurements were made using data from the USArray Transportable Array. At periods of 16 s and 18 s, which mostly sample the crust, we find a region of low anisotropy in central Nevada coinciding with locally reduced phase velocities, and surrounded by a semi-circular pattern of fast seismic directions. Away from central Nevada the fast directions are ～ N–S in the eastern Great Basin, NW–SE in the Walker Lane region, and they transition from E–W to N–S in the northwestern Great Basin. Our short-period phase velocity maps, combined with recent crustal receiver function results, are consistent with the presence of a semi-circular anisotropy signal in the lithosphere in the vicinity of a locally thick crust. At longer periods (28–102 s), which sample the uppermost mantle, isotropic phase velocities are significantly reduced across the study region, and fast directions are more uniform with an ～ E–W fast axis. The transition in phase velocities and anisotropy can be attributed to the lithosphere–asthenosphere boundary at depths of ～ 60 km. We interpret the fast seismic directions observed at longer periods in terms of present-day asthenospheric flow-driven deformation, possibly related to a combination of Juan de Fuca slab rollback and eastward-driven mantle flow from the Pacific asthenosphere. Our results also provide context to regional SKS splitting observations. We find that our short-period phase velocity anisotropy can only explain ～ 30% of the SKS splitting times, despite similar patterns in fast directions. This implies that the origin of the regional shear-wave splitting signal is complex and must also have a significant sublithospheric component.     

High temperatures and inverted metamorphism in the schist of Sierra de Salinas,California

New field and thermobarometric work in the Californian Salinian block clarifies current and pre-Tertiary relationships between the schist of Sierra de Salinas and Cretaceous arc-related granitic rocks. The contact is variably preserved as a brittle fault and high-temperature mylonite zone, the Salinas shear zone, which represents the contact between North America and sediments accreted above the Farallon slab between ∼ 76 Ma and ∼ 70 Ma. Near granulite facies, prograde replacement of hornblende with clinopyroxene is associated with deformation of plutonic rocks at the base of the upper plate. In the lower plate, the schist of Sierra de Salinas, garnet–biotite thermometry indicates decreasing temperatures down-section from at least 714 °C to ∼ 575 °C over an exposed thickness of ∼ 2.5 km, consistent with petrologic evidence of an inverted metamorphic gradient. The measured temperatures are significantly higher than observed at shallow levels above subducting slabs or predicted by 2D computational models assuming low shear stresses. Previous workers have called upon shear heating to explain similar observations in the correlative Pelona schist, an unlikely scenario given the results of recent rock deformation experiments which predict that feldspar–quartz–mica aggregates are far too weak to withstand stresses of ∼ 70 MPa required by the shear heating hypothesis. As an alternative, we propose that high temperatures resulted from conductive heating while the leading edge of the schist traveled ∼ 150 km beneath the recently active Salinian continental arc during the initiation of shallow subduction. Weakening of the schist due to high temperatures helped facilitate the collapse of the Salinian arc as the schist was emplaced. Schist emplacement coincided with loss of lower, mafic portions of the arc, and therefore evolution of the Southern California crust towards a more felsic composition.     

Estimation of Shallow S-Wave Velocity Structure in Damascus City,Syria, Using Microtremor Exploration

Array measurements of microtremors were carried out at thirty sites in Damascus city, Syria to estimate S-wave velocity structures of shallow soil formations for site effect analysis. The microtremor data were recorded by 6 vertical-component seismometers distributed along the circumferences of two circles as well as a 3-component seismometer deployed in the center. The phase velocities were estimated at each site from the vertical components of recorded microtremor data by using the Spatial Autocorrelation method. Then, Genetic Simulated Annealing Algorithm technique was applied for inversion of the phase velocities to estimate 1-D S-wave velocity structures beneath the sites. The inverted Vs profiles are not uniform in Damascus city and the results show that a shallow soft layer (∼200 m/s) appears in the eastern part of the city as well as the central part along Barada River. This layer controls the amplification distribution in the city with a high amplification mainly observed at the locations having this layer. The inversion results also show that the depth to the engineering bedrock (∼750 m/s) is very shallow along the foothills of Mt. Qasyoun in the north-west. Then the depth increases towards the east and the south. The maximum depth to the engineering bedrock (∼80 m) was observed in the southern part of Damascus. To validate the results of the inversions, the spectral ratios between the horizontal and vertical components (H/V) of the recorded microtremor data at the central seismometer were compared with the computed ellipticities of the fundamental-mode Rayleigh-waves based on the respective Vs structure. The results show a good agreement in a period range of 0.05 s to 0.5 s. In this period range, the dominant peaks of the H/V ratios are due to the overall effect of the velocity contrasts between the shallow layers representing the subsurface S-wave velocity structure. Moreover, the average S-wave velocity for the top 10 m of soils (V_S10) shows a better correlation with the averaged site amplification in a period range of 0.05 s to 0.5 s than V_S30 which indicates that V_S10 can be a better proxy for high-frequency site amplification in the case of Damascus city.     

Ancient Os isotope signatures from the Ontong Java Plateau lithosphere: Tracing lithospheric accretion history

In order to better understand the nature and formation of oceanic lithosphere beneath the Early Cretaceous Ontong Java Plateau, Re–Os isotopes have been analysed in a suite of peridotite xenoliths from Malaita, Solomon Islands. Geological, thermobarometric and petrological evidence from previous studies reveal that the xenoliths represent virtually the entire thickness of the southern part of subplateau lithospheric mantle (< 120 km). This study demonstrates that vertical Os isotopic variations correlate with compositional variations in a stratified lithosphere. The shallowest plateau lithosphere (< 85 km) is dominated by fertile lherzolites showing a restricted range of ¹⁸⁷Os/¹⁸⁸Os (0.1222 to 0.1288), consistent with an origin from ~ 160 Ma Pacific lithosphere. In contrast, the basal section of subplateau lithospheric mantle (~ 95–120 km) is enriched in refractory harzburgites with highly unradiogenic ¹⁸⁷Os/¹⁸⁸Os ratios ranging from 0.1152 to 0.1196, which yield Proterozoic model ages of 0.9–1.7 Ga. Although the whole range of Os isotope compositions of Malaita peridotites is within the variations seen in modern abyssal peridotites, the contrasting isotopic compositions of shallow and deep plateau lithosphere suggest their derivation from different mantle reservoirs. We propose that the subplateau lithosphere forms a genetically unrelated two-layered structure, comprising shallower, typical oceanic lithosphere underpinned by deeper impinged material, which included a component of recycled Proterozoic lithosphere. The impingement of residual but chemically heterogeneous mantle, mechanically coupled to the recently formed, thin lithosphere, may have a bearing on the anomalous initial uplift and late subsidence history of the seismically anomalous plateau root.     

Crustal structure of the Reykjanes Ridge near 62°N,on the basis of seismic refraction and gravity data

Explosion deep seismic sounding data sections of high quality had been obtained with RV Meteor in the Reykjanes Iceland Seismic Project (RRISP77 [Angenheister, G., Gebrande, H., Miller, H., Goldflam, P., Weigel, W., Jacoby, W.R., Pálmason, G., Björnsson, S., Einarsson, P., Pavlenkova, N.I., Zverev, S., Litvinenko, I.V., Loncarecic, B., Solomon, S., 1980. Reykjanes Ridge Iceland Seismic Experiment (RRISP 77). J. Geophys. 47, 228–238]) which close an information gap near 62°N. Preliminary results were presented by Weigel [Weigel, W., 1980. Aufbau des Reykjanes Rückens nach refraktionsseismischen Messungen. In: Weigel, W. (Ed.), Reykjanes Rücken, Island, Norwegischer Kontinentalrand. Abschlusskolloquium, Hamburg zur Meteor-Expedition, vol. 45. DFG, Bonn, pp. 53–61], and here we report on the data and results of interpretation. Clear refracted phases to 90 km distance permit crustal and uppermost mantle structure to be modelled by ray tracing. The apparent P-wave velocities are around 4.5, 6–6.5, 7–7.6 and 8.2–8.7 km/s, but no wide-angle reflections have been clearly seen. Accompanying sparker reflection data reveal thin sediment ponds in the axial zone and up to 400 m thick sediments at 10 Ma crustal age. Ray tracing reveals the following model below the sediments: (1) a distinct, 1–2 km thick upper crust (layer 2A) with Vp increasing with age (to 10 Ma) from <3.4 to 4.9 km/s and with a vertical gradient of 0.1–0.2 km/s/km, (2) a lower crust or layer 3 beginning at depths of 2 (axis) to 4 km (10 Ma age) below sea level with 6.1–6.8 km/s and similar vertical gradients as above, (3) the lower crust bottoms at 5.2–9.5 km depth below sea level (0–10 Ma) with a marked discontinuity, underneath which (4) Vp rises from about 7.5–7.8 km/s (0–10 Ma) with a positive vertical gradient of, again, 0.1–0.2 km/s/km such that 8 km/s would be reached at 12 km and deeper near the axis. Our preferred interpretation is that the mantle begins at the distinct discontinuity (“Moho”), but a deeper “Moho” of Vp ≈ 8 km/s cannot be excluded. From Iceland southward to 60°N several experiments show a decrease of crustal thickness from 14 to 8 km. Velocity trends with age across the ridge reflect cooling and filling of cracks, and thickness trends probably suggest volcanic productivity variations as previously suggested.Gravity inversion concentrates on a profile across the ridge with the above seismic a priori information; with 0.2–0.5 km depth uncertainty it leads to a good fit (±2.5 mGal where seismic data exist). Best fitting densities are (in kg/m³) for sediments, 2180; upper crust, 2450–2570; lower crust, 2850–2940; mantle lithosphere, 3215–3240 with a deficit for an asthenospheric wedge of no more than −100 kg/m³. The morphological ridges and troughs superimposed on the SE ridge flank are partly correlated, partly anti-correlated with the Bouguer anomaly and suggest that variable crustal density variations accompany the morphology variations.