Intermediate-depth Earthquakes Research Articles

Automatic relocation of intermediate-depth earthquakes using adaptive teleseismic arrays

SUMMARY Intermediate-depth earthquakes, accommodating intraslab deformation, typically occur within subduction zone settings at depths between 60–300 km. These events are in a unique position to inform us about the geodynamics of the subducting slab, specifically the geometry of the slab and the stress state of the host material. Improvements in the density and quality of recorded seismic data enhance our ability to determine precise locations of intermediate-depth earthquakes, in order to establish connections between event nucleation and the tectonic setting. Depth phases (near-source surface reflections, e.g. pP and sP) are crucial for the accurate determination of earthquake source depth using global seismic data. However, they suffer from poor signal-to-noise ratios in the P wave coda. This reduces the ability to systematically measure differential traveltimes to the direct P arrival, particularly for the frequent lower magnitude seismicity which highlights considerable seismogenic regions of the subducted slabs. To address this limitation, we have developed an automated approach to group globally distributed stations at teleseismic distances into ad-hoc arrays with apertures of 2.5$^\circ$, before optimizing and applying phase-weighted beamforming techniques to each array. Resultant vespagrams allow automated picking algorithms to determine differential arrival times between the depth phases and their corresponding direct P arrival. Using these differential times we can then determine the depths of earthquakes, which in turn can be used to create a catalogue of relocated events. This will allow new comparisons and insights into the governing controls on the distribution of earthquakes in subducted slabs. We demonstrate this method by relocating intermediate-depth events associated with northern Chile and the Peruvian flat slab regions of the subducting Nazca plate. The relocated Chilean catalogue contains comparable event depths to an established catalogue, calculated using a semi-automated global methodology, which serves to validate our fully automatic methodology. The new Peruvian catalogue we generate indicates three broad zones of seismicity approximately between latitudes 1–7$^\circ$S, 7–13$^\circ$S and 13–19$^\circ$S. These align with flat to steep slab dip transitions and the previously identified Pucallpa Nest. We also find a regionally deeper slab top than indicated by recent slab models, with intraslab events concentrated at points where the slab bends, suggesting a link between slab flexure and intermediate-depth earthquake nucleation.

Comparison between rupture parameters of intermediate and deep earthquakes at the Peru–Brazil–Bolivia border and northern Chile

We determined the main parameters of the source rupture process of intermediate- and deep-depth earthquakes occurring in the Peru–Brazil–Bolivia border region and northern Chile. The parameters of depth, fault-plane orientation, scalar seismic moment, slip distribution, and radiated seismic energy are obtained from seismograms. We selected 15 intermediate-depth earthquakes (100 km < h < 300 km) and 10 very deep earthquakes (h > 500 km) with magnitudes MW ≥ 6.0. For most events, the slip distribution over the rupture plane shows a single asperity, and the source time function presents a simple pulse. There are differences between intermediate-depth and deep earthquakes. The rupture areas, maximum slip and source time function (STF) duration are larger for intermediate-depth events than for deep events. Additionally, the STF’s show a sharper increase for deep earthquakes. The scaled radiated seismic energy shows larger values for deep depth events. The stress regime pattern derived from the obtained focal mechanism agrees with the geometry of the subduction of the Nazca plate. At intermediate depths, in the northern area up to 12°S, the stress pattern corresponds to a horizontal extension, while in the southern area, the tension axes dip at an angle of 30°. At deep depths, the stress regime corresponds to vertical compression in the north and dips of approximately 45° in the south.

Anisotropic thermal conductivity of antigorite along slab subduction impacts seismicity of intermediate-depth earthquakes

Double seismic zones (DSZs) are a feature of some subducting slabs, where intermediate-depth earthquakes (~70–300 km) align along two separate planes. The upper seismic plane is generally attributed to dehydration embrittlement, whereas mechanisms forming the lower seismic plane are still debated. Thermal conductivity of slab minerals is expected to control the temperature evolution of subducting slabs, and therefore their seismicity. However, effects of the potential anisotropic thermal conductivity of layered serpentine minerals with crystal preferred orientation on slab’s thermal evolution remain poorly understood. Here we measure the lattice thermal conductivity of antigorite, a hydrous serpentine mineral, along its crystallographic b- and c-axis at relevant high pressure-temperature conditions of subduction. We find that antigorite’s thermal conductivity along the c-axis is ~3–4 folds smaller than the b-axis. Our numerical models further reveal that when the low-thermal-conductivity c-axis is aligned normal to the slab dip, antigorite’s strongly anisotropic thermal conductivity enables heating at the top portion of the slab, facilitating dehydration embrittlement that causes the seismicity in the upper plane of DSZs. Potentially, the antigorite’s thermal insulating effect also hinders the dissipation of frictional heat inside shear zones, promoting thermal runaway along serpentinized faults that could trigger intermediate-depth earthquakes.

Stress drops of intermediate-depth intraslab earthquakes beneath Tohoku, northern Japan

We calculated stress drops for 2875 small intraslab earthquakes at intermediate depths beneath Tohoku, Japan. We applied an S-coda-wave spectral ratio method to almost 900,000 event pairs. Detailed velocity values for the oceanic crust (OC) were adopted from previous observational studies. The median stress drops in the OC are about half those in the oceanic mantle (OM). The median stress drop for earthquakes in the OC decreases from depths of 70 to 120 km and increases from 120 to 170 km. Our preferred interpretation is that the rigidity in the OC decreases and then increases with depth due to combined effects of the dehydration associated with the eclogite formation and the increasing temperature with depth. These depth variations are consistent with results of a similar study beneath Hokkaido. The median stress drops in the oceanic plate beneath Tohoku are generally smaller than those beneath Hokkaido. Previous studies imaging the seismic structure at shallow depths and b-value analyses of intraslab earthquakes indicate that the near-trench region of the oceanic plate off Tohoku is more hydrated than that off Hokkaido. Taken together, these results suggest that differences in the degree of hydration of the oceanic plate in the near-trench regions could produce the different behaviors of stress drops of intermediate-depth earthquakes observed in Tohoku and Hokkaido.Graphical

Vertical Elastic Acceleration Response Spectra for Vrancea Intermediate-Depth Earthquakes

This study is focused on the evaluation of vertical elastic acceleration response spectra for Vrancea intermediate-depth earthquakes in Romania in the context of the ongoing update of the national seismic design code. The ground motion database employed in this research consists of about 500 ground motions recorded during moderate and large Vrancea intermediate-depth earthquakes with moment magnitudes MW ≥ 5.2. The analysis of the dataset showed that a single ground motion recording had a peak ground acceleration larger than 0.20 g. The results of the analyses showed that no significant differences between the control periods of the vertical elastic response spectra as a function of the site class could be inferred. It was also observed that the mean value of the amplification factor computed for the entire ground motion database was about 2.5, irrespective of the earthquake magnitude, site class, or level of horizontal peak ground acceleration. However, larger-magnitude earthquakes generate larger spectral amplifications in the medium- and long-period ranges. The analysis of the ground motions recorded in Bucharest area revealed a magnitude dependency of the control period for the vertical ground motion TC,v. Finally, a single spectral shape for vertical acceleration response spectra characterized by a maximum dynamic amplification factor of 2.5 and control periods TBv and TCv of 0.05 s and 0.60 s is proposed for design purposes. This aspect allows for a major update from the current version of the national seismic design code which proposes control periods of the elastic vertical response spectrum dependent on the horizontal ones. The short-period (at T = 0.2 s) and long-period (at T = 1.0 s) ratios of the vertical acceleration response spectrum to the corresponding horizontal ones are 0.50 and 0.40.

A weak subducting slab at intermediate depths below northeast Japan.

Knowledge of the state of stress in subducting slabs is essential for understanding their mechanical behavior and the physical processes that generate earthquakes. Here, we develop a framework which uses a high-resolution focal mechanism catalog to determine the change in the position of the neutral plane before and after the M9 Tohoku-oki earthquake to determine that the deviatoric stress within the slab at intermediate depths must be very low (∼1 MPa). We show that by combining the static stress calculated from coseismic slip distributions with the stress orientations before and after the mainshock, we can determine the full deviatoric stress tensor within the subducting slab at intermediate depths. These results preclude earthquake source mechanisms that require large background driving stresses, favoring a mechanically weak subducting slab, thus providing quantitative constraints on the physical processes that generate intermediate-depth earthquakes.

A Comprehensive Stress Drop Map From Trench to Depth in the Northern Chilean Subduction Zone

AbstractWe compute stress drops for earthquakes in Northern Chile recorded between 2007 and 2021. By applying two analysis techniques, (a) the spectral ratio (SR) method and (b) the spectral decomposition (SDC) method, a stress drop map for the subduction zone consisting of 51,510 stress drop values is produced. We build an extended set of empirical Green’s functions (EGF) for the SR method by systematic template matching. Outputs are used to compare with results from the SDC approach, where we apply cell‐wise obtained global EGF's to compensate for the structural heterogeneity of the subduction zone. We find a good consistency of results of the two methods. The increased spatial coverage and quantity of stress drop estimates from the SDC method facilitate a consistent stress drop mapping of the different seismotectonic domains. Albeit only small differences of median stress drop, strike‐perpendicular depth sections clearly reveal systematic variations, with earthquakes at different seismotectonic locations exhibiting distinct values. In particular, interface seismicity is characterized by the lowest observed median value, whereas upper plate earthquakes show noticeably higher stress drop values. Intermediate depth earthquakes show comparatively high average stress drop and a rather strong depth‐dependent increase of median stress drop. Additionally, we observe spatio‐temporal variability of stress drops related to the occurrence of the two megathrust earthquakes in the study region. The presented study is the first coherent large scale 3D stress drop mapping of the Northern Chilean subduction zone. It provides an important component for further detailed analysis of the physics of earthquake ruptures.

Numerous Tibetan lower-crustal and upper-mantle earthquakes, detected by Sn/Lg ratios, suggest crustal delamination or drip tectonics

Whether intermediate-depth earthquakes beneath Tibet and the Himalaya are in the lower crust or upper mantle should provide insight into rheology, hence composition and temperature, of continental lower crust and upper mantle. Based on the waveguide theory that earthquakes respectively above or below the Moho will excite higher Lg or Sn energy, we develop an Sn/Lg amplitude-ratio analysis using a single permanent station to investigate earthquake depths with respect to the Moho. First, we use synthetics to show the ubiquitous sharp increase of the Sn/Lg ratio for hypocenters beneath the crust. A deep-crustal high-vs layer appropriate for eclogitized Indian lower crust causes a more gradual increase of Sn/Lg ratios for hypocenters within or below the eclogite layer. We measured Sn/Lg ratios of 595 Tibetan earthquakes from 1998 to 2022 with nominal (catalog) hypocentral depths >30 km and magnitudes >3.2 in five regions (west Tibet, southwest Tibet, south Tibet, southeast Tibet, and Qiangtang). As predicted by our normal-mode synthetics, earthquakes in west and south Tibet show Sn/Lg ratios that increase sharply as catalog depths approach the independently determined receiver-function Moho. Numerous earthquakes with high Sn/Lg ratios and nominal depths around the reported Moho are identified between the Karakoram fault and Altyn-Tagh fault (ATF) in west Tibet, which we attribute to an eclogitized Indian lower crust extending north to the ATF. Impingement of this underthrust Indian slab against the Tarim craton and consequent eclogite delamination or dripping may cause the observed below-Moho seismicity which occupies a region 100 km across and 200 km along orogenic strike and extends 30–50 km below published Moho depths. Similarly dense seismicity spans the Moho and occupies a similar volume below the Moho beneath the High Himalaya in south Tibet. The south Tibet Sn/Lg patterns may indicate an eclogitized Indian layer beginning to delaminate or drip south of the Yarlung-Zangpo suture. Our southeast Tibet and southwest Tibet Sn/Lg observations are less clear cut, perhaps due to less appropriate epicentral distances from the available observing stations. The Qiangtang region of Tibet likely has sparse mid-to-lower-crustal earthquakes, but no definitive below-Moho earthquakes. Our work expands the catalog of continental intraplate earthquakes below the Moho by >100 events with mb>3.2, as well as identifying tens of mb>3.2 lower-crustal events, so that Tibet fits neither an ideal crême-brulée nor an ideal jelly-sandwich crustal-strength paradigm. Two clusters of earthquakes that each span the Moho may be best explained by the delamination or dripping of a strong eclogitized lower crust into the upper mantle.

Impact of chlorite dehydration on intermediate-depth earthquakes in subducting slabs

Intermediate-depth earthquakes are common in the double seismic structures of many subduction zones under high pressures (~1–4 GPa). Serpentine dehydration exhibits well-established links with double seismic zone earthquakes. Additionally, dehydration of several hydrous minerals including lawsonite and chlorite underlying the upper and lower layers, respectively, may be responsible for intermediate-depth earthquakes. Here, we present experimental evidence suggesting that chlorite dehydration can trigger intermediate-depth earthquakes at the lower plane (~700 °C). We conducted deformation experiments on chlorite peridotite under high-pressure (0.5–2.5 GPa) and high-temperature (500–750 °C) using a modified Griggs apparatus. Experiments revealed the presence of faults in samples that had undergone partial chlorite dehydration with the presence of the dehydration product Ca-amphibole along these faults. Our findings confirm, together with correlation studies between seismicity and mineral stability, that a part of intermediate-depth seismicity in the lower plane of double seismic zones can be attributed to chlorite dehydration.

Tectonic setting of the northwestern andes Constrained by a high-resolution earthquake catalog: Block Kinematics

A high-precision earthquake catalog was generated using source-specific station terms and waveform cross-correlation techniques. This detailed catalog serves to interpret the crustal structure and deformation related to the tectonic setting in Northwestern South America. The Panamá-Chocó Block (PCB) is in contact with the Northwestern Andes causing crustal deformation and hence, fault interaction in the brittle regime leading to a high earthquake production concentrated in the Murindó Seismic Cluster (MSC) which is located in the northeastern corner of the PCB. Slip on the Uramita Fault Zone (UFZ) is a response to the accommodation of the PCB against the North Andean block (NAB), splitting the fault into two segments, marked by an abrupt change in its strike. The north segment of the UFZ, a dextral transform fault (strike = S52°E and near vertical plane), which produced the Mutatá earthquake (2016-09-14, Mw = 6.2). The central segment of the UFZ strikes N12°E and exhibits mainly reverse slip with a subvertical fault plane which produced the 1987-03-19 Mw = 5.4 and 1987-11-11 Mw = 5.3 events. The Murindó sinistral strike-slip fault (strike = N9°W), which produced the great Murindó earthquake (1992-10-18, Mw = 7.1, epicentral intensity XI) is characterized by its wide zone of brittle deformation as a response to internal deformation of the PCB under a compressional NW-SE stress. Intermediate depth earthquakes are widely observed in the area, where we find evidence of the Caribbean slab subducting southwards beneath Panamá. The subduction of the Malpelo microplate to the south, shows high productivity within the Cauca Cluster and the presence of at least three finger-shaped mantle zones of seismicity.

In situ study of the reaction phase A plus high-P clinoenstatite to forsterite plus water at reduced water activity

Abstract. We examined the reaction phase A plus high-P clinoenstatite to forsterite plus water (Reaction R1) by means of in situ X-ray diffraction measurements with the large volume press at the synchrotron PETRA III, Hamburg. Contrary to the study of Lathe et al. (2022), in which all experiments on Reaction (R1) were performed at a water activity of 1, the reversed experiments presented in this study were performed at reduced water activity with mole fractions of about XH2O= XCO2=0.5. The intention of this investigation was to test the observation made by Perrillat et al. (2005), which was that dehydration reactions are kinetically faster at reduced than under water-saturated conditions. The position of Reaction (R1) at the reduced conditions was determined by reversal brackets at 9.1 and 9.5 GPa (630 and 700 ∘C), at 9.7 and 10.0 GPa (725 and 700 ∘C), at 9.8 and 10.2 GPa (675 and 750 ∘C), and at 10.5 GPa (675 and 740 ∘C). Additionally, we performed two offline experiments with brackets at 10.0 and 10.6 GPa (750 and 700 ∘C, respectively) that are in agreement with the results of the in situ experiments. We do not observe any “intermediate” precursor phase in our experiments. The equilibrium of Reaction (R1) is shifted by about 100 ∘C to lower temperature compared to the results under water-saturated conditions. Thus, at a water activity (aH2O) below 1 the phase A plus clinoenstatite dehydration reaction can only occur in extremely cold subduction slabs. The kinetics of Reaction (R1) dehydration at reduced water activity is slower than that determined previously by Lathe et al. (2022) under water-saturated conditions. Thus, the above-mentioned hypothesis of Perrillat et al. (2005) could not be confirmed. However, in both of our studies on Reaction (R1), the newly formed dehydration product forsterite was of nanometer size, which supports earlier experimental observations, which is that product phases of dehydration reactions are generally very fine-grained and might promote the concept that intermediate-depth earthquakes in subduction zones are initiated by mechanical instabilities from extremely fine-grained materials formed during dehydration reactions.

Upper mantle melt caused by a subducted slab in the Indian-Eurasian continental subduction zone

A low-velocity layer atop the mantle transition zone has been extensively observed worldwide. In subduction zones, this layer is widely explained as partial melting triggered by slab subduction on a regional or global scale. However, direct observational evidence is still absent, and the response of the layer to slab subduction is not well known. Here, we image the seismic velocity around the mantle transition zone by matching synthetic and observed triplicated seismic P and sP waveforms in the Indian–Eurasian continental subduction zone. Our observations reveal a laterally varied low-velocity layer atop the mantle transition zone beneath the Hindu Kush, where a subducted slab extends to the mantle transition zone. The geometric morphology of the low-velocity layer documents that it is a partially molten layer induced by the subducted slab on a regional scale. Interestingly, our observations also support that the layer has a low viscosity. The decreased viscosity may facilitate slab motion in the deep mantle, contributing to slab stretching, tearing and break-off and its resultant rare recurring large intermediate-depth earthquakes in an intracontinental setting.

The rupture plane of the February 2022 Mw 6.2 Guatemala, intermediate depth earthquake

An intermediate depth intraplate earthquake with Mw 6.2 was generated in the Guatemalan subduction zone on 16 February 2022 with epicenter to Southwest of the department of Escuintla. More than 275 aftershocks were registered, which were relocated with the HypoDD algorithm, being able to identify a fault with an area of ~350 km2, which is considerably higher than expected for an earthquake of that magnitude. The moment tensor at the centroid of the main earthquake and the estimation of other focal mechanisms of the largest aftershocks, allowed us to identify extension earthquakes, related to the fault plane, and compression earthquakes that were associated with seismicity on the upper part of the slab. The region of the sequence has presented high seismic activity in recent years. It is proposed that the mainshock nucleated in the lower seismicity layer (LSL) of the double seismicity zone proposed for the region, triggering seismic activity on a pre-existing active fault, also triggering seismic activity in the upper seismicity layer (USL). The separation between these seismicity layers was estimated to be 12.2±5.0 km.

Seismicity around the boundary between eastern and western Makran

The Makran subduction zone is segmented into western and eastern parts. The small number of teleseismic and regionally-recorded earthquakes makes it difficult to detect possible seismic activity along the border between eastern and western Makran. Additionally, the scarce seismicity does not allow us to investigate seismic deformation of the fore-arc Jazmurian Depression and the Sistan Suture Zone, understand the relationship between the intermediate-depth seismicity and the subducting plate, and estimate the seismic moment release in Makran. To address these questions, we installed a temporary seismic network of 39 stations around the border between eastern and western Makran from June 2016 to November 2019. The observed seismicity shows a NNW-SSE trending alignment of shallow small-magnitude earthquakes within the northern part of the accretionary prism along the longitude of ∼ 60°. Earthquakes are more frequent to the east of the trend, implying the trend may be associated with the approximate boundary between western and eastern Makran. The new seismicity map shows that the Jazmurian Depression is surrounded by shallow crustal earthquakes, suggesting that it is an aseismic block. Scattered shallow seismicity is also observed within the Sistan Suture Zone, but seismicity is more concentrated along its border with the Lut block, which indicates a relative motion between them. We show that the intermediate-depth earthquakes occur within the subducting plate, and the rate of seismic moment release in the Makran megathrust zone is much less than those of the Alaska, Nankai, Chile, Kuril, and Mexico megathrust zones.

Velocity structure (1D) and earthquakes relocation in the flat-slab to normal plate transition zone in the Argentinean transarc

The Argentinean Andean foreland region between latitudes 31°S to 33.5°S presents an intense crustal and upper mantle seismic activity that was recorded by a local seismological network called CHARSME (Chile Argentina SeisMic Experiment), and was used to obtain new velocity determinations. In a first moment, from 4 P-wave velocity (VP) models proposed by other authors, the detail in the layers was increased, resulting in better models used as initial model in the investment. Finally, by carefully selecting a set of earthquakes, a simultaneous velocity inversion, hypocentre determination and station correction were performed using the VELEST code. In order to obtain better velocity models (VP and VS) and VP/VS ratio, different tests were carried out in all processes.The final VP model presented in this paper incorporates considerations made in other seismological studies. This model accurately determined velocities between 10 and 65 km depth. The upper mantle velocity (VP) was defined in 8.05 km/s. From the Wadatti diagram, which represents a line whose slope is related to the wave velocities and the Poisson's coefficient, a VP/VS value of 1.732 was obtained. Alternative VP/VS values were also considered to invert VP, VS and VP/VS simultaneously. The minimum detailed model of VP, VS and VP/VS decreased the RMS by 25% relative to the initial models, implying an improvement in the uncertainty of locations. The new and precise locations of crustal earthquakes allowed defining clusters within Cuyania, a tectonic terrain that collided to Gondwana in the Early Paleozoic. Although it is known that foreland seismicity is mainly concentrated in this terrain, this work shows that there are well-defined zones within Cuyania that are more seismically active. Earthquakes of intermediate depth (upper mantle) are strongly agglomerated in the flat-slab zone of the subducting Nazca Plate.The discontinuities present in our model suggest upper, middle and lower crustal structuring, as well as a marked contrast between the sediments and the underlying basement. It was possible to interpret an upper crust consisting of sediments and basement to a depth of 25 km, a middle crust between 25 and 45 km and a lower crust between 45 and 55 km. We also defined the Mohorovičić (Moho) discontinuity at 55 km depth. The crustal velocities are in agreement with previous seismological studies.

Continental subduction of Adria in the Apennines and relation with seismicity and hazard

The subduction of continental lithosphere is a complex process because the buoyancy of the crust is higher than the oceanic and should resist sinking into the mantle. Anyway, studies on the Alpine-Himalayan collision system indicate that a large portion of the continental crust is subducted, while some material is accreted in the orogens. The Apennine is a perfect case for studying how such processes evolve, thanks to high quality seismic images that illuminate a critical depth range not commonly resolved in many collisional settings. In this paper, we show the structure of the Apennines orogen, as jointly revealed by seismicity and deep structure from regional and teleseismic tomography and receiver function profiles. The westward subducting Adria lithosphere is well defined along the orogen showing a mid-crustal delamination. Seismicity within the underthrusting lower crust and velocity anomalies in the mantle wedge highlight how the subduction evolution is entangled with the liberation of fluids. The eclogitization of subducted material enhances the fluid release into the wedge, the delamination and retreat of the Adria plate. This delamination/subduction generates a coupled compression and extension system that migrates eastward following the retreat of the lithosphere, with broad sets of normal faults that invert or interfere with pre-existing compressional structures all over the roof plate. The sparseness and non-ubiquity of intermediate depth earthquakes along the subduction panel suggest that the brittle response of the subducting crust is governed by its different composition and fluid content. Therefore, the lower crust composition appears essential in conditioning the evolution of continental subduction.

A Cascadia Slab Model From Receiver Functions

AbstractWe map the characteristic signature of the subducting Juan de Fuca and Gorda plates along the entire Cascadia forearc from northern Vancouver Island, Canada, to Cape Mendocino in northern California, USA, using teleseismic receiver functions. The subducting oceanic crustal complex, possibly including subcreted material, is parameterized by three horizons capable of generating mode‐converted waves: a negative velocity contrast at the top of a low velocity zone underlain by two horizons representing positive contrasts. The amplitude of the conversions varies likely due to differences in composition and/or fluid content. We analyzed the slab signature for 298 long‐running land seismic stations, estimated the depth of the three interfaces through inverse modeling and fitted regularized spline surfaces through the station control points to construct a margin‐wide, double‐layered slab model. Crystalline terranes that act as the static backstop appear to form the major structural barrier that controls the slab morphology. Where the backstop recedes landward beneath the Olympic Peninsula and Cape Mendocino, the slab subducts sub‐horizontally, while the seaward‐protruding and thickened Siletz terrane beneath central Oregon causes steepening of the slab. A tight bend in slab morphology south of the Olympic Peninsula coincides with the location of recurring large intermediate depth earthquakes. The top‐to‐Moho thickness of the slab generally exceeds the thickness of the oceanic crust by 2–12 km, suggesting thickening of the slab or underplating of slab material to the overriding North American plate.

Ground Motion Duration Patterns for Vrancea (Romania) Intermediate-Depth Earthquakes

This study is focused on evaluating ground motion durations of Vrancea intermediate-depth earthquakes in Romania, in the context of future updates to the Romanian seismic design code P100-1/2013. The ground motion database compiled for this study consists of about 200 ground motions recorded during five moderate and large Vrancea intermediate-depth earthquakes that occurred in the period of 1977–2004 and had moment magnitudes of MW ≥ 6.0. Two empirical models were derived in this study for the significant ground motion duration considering two time intervals (5–75% and 5–95%) for the accumulation of the Arias Intensity IA. An analysis of the data shows that the mean ratio between D5-95 and D5-75 is about 2.8. Moreover, the regression also shows that the largest share of variability is due to the within-event component (site term). Among the regression coefficients, the hypocentral distance and the soil conditions appear to have a larger impact on the ground motion duration compared to the earthquake magnitude. It was also observed that the median ground motion durations predicted using the empirical model proposed in this study were much smaller than the ones from the proposed Eurocode 8 draft for the same magnitude range. Finally, geographic trends related to the distribution of residuals were also evaluated using the data from the three earthquakes with the largest number of available ground motion recordings.

Assessment of Liquefaction Hazard for Sites in Romania Using Empirical Models

This paper is focused on the evaluation of the liquefaction hazard for different sites in Romania. To this aim, a database of 139 ground motions recorded during Vrancea intermediate-depth earthquakes having moment magnitudes MW ≥ 6.0 is employed for the evaluation of the equivalent number of cycles for this seismic source. Several functional forms for the empirical evaluation of the equivalent number of cycles considering various seismological or engineering parameters are tested and evaluated. The regression analysis shows smaller uncertainties for the empirical models based on ground motion engineering parameters. Considering the lack of information in terms of engineering parameters, a simpler empirical model which accounts for the earthquake magnitude, source–site distance and soil conditions is selected for the liquefaction hazard analysis. Based on the proposed empirical model, specific magnitude scaling factors for Vrancea intermediate-depth earthquakes are proposed for the first time as well. The liquefaction hazard analysis is performed for sites whose seismic hazard is generated by either the Vrancea intermediate-depth seismic source or by local shallow crustal seismic sources. In the case of some of the selected sites, liquefaction phenomena were observed during past large-magnitude earthquakes. Unlike previous studies dealing with liquefaction analyses for sites in Romania, in this research, the hazard assessment is performed for various ground motion levels evaluated based on probabilistic seismic hazard assessment. Liquefaction hazard curves are constructed for each analyzed site. The results of the liquefaction hazard analysis show that this phenomenon is more likely to occur in the areas exposed to Vrancea intermediate-depth earthquakes, compared to the areas affected by local shallow earthquakes. In the case of the analyzed soil profiles from Bucharest, Craiova and Ianca, the minimum liquefaction safety factors less than one even for seismic hazard levels having mean return periods of 100 years and less.

Small-scale heterogeneities in the convective upper mantle beneath circum-Pacific subduction zones: Evidence for fragments of recycled basaltic crust

We show evidence for the presence of compositional heterogeneities at 10 km scale in the convective upper mantle below the circum-Pacific subduction zones. By array-processing we examine seismograms of intermediate-depth earthquakes recorded by large-aperture seismic networks at teleseismic distances. In the P coda until about 25 s after direct P waves, we observe anomalous and significant signals arriving from the directions close to the P waves. The amplitude and duration of the signals decrease with focal depth within the range from about 180 to 300 km. We show that the signals arise mainly from S-to-P scattering in the mantle below the foci. Using the scattering theory for randomly heterogeneous media we also show that the observed features of the signals are matched by 5 to 10 km-scale density and rigidity anomalies that diminish below about 400 km depth. The rigidity anomaly is of the same sign and nearly twice as large as the density anomaly. These properties of the scattering objects could be consistent with basalt embedded in depleted mantle. The results indicate that the upper mantle transition zone beneath subducting slabs is not well mixed even by vigorous convection.