Crustal Architecture of Puerto Rico Using Body-Wave Seismic Tomography and High-Resolution Earthquake Relocation

Puerto Rico is a highly seismically active island, where several damaging historical earthquakes have occurred and frequent small events persist. It situates at the boundary between the Caribbean and North American plates, featuring a complex fault system. Here, we investigate the seismotectonic crustal structure of the island by interpreting the 3D compressional-wave velocity VP and compressional- to shear-wave velocity ratio VP/VS models and by analyzing the distribution of the relocated earthquakes. The 3D velocity models are obtained by applying the simul2000 tomographic inversion algorithm based on the phase arrivals recorded by the Puerto Rico seismic network. We find high-VP and low-VP/VS anomalies in the eastern and central province between the Great Northern Puerto Rico fault zone and the Great Southern Puerto Rico fault zone, correlating with the Utuado pluton. Further, there are low-VP anomalies beneath both the Great Southern Puerto Rico fault zone and the South Lajas fault, indicating northerly dipping structures from the southwest to the northwest of the island. We relocate 19,095 earthquakes from May 2017 to April 2021 using the new 3D velocity model and waveform cross-correlation data. The relocated seismicity shows trends along the Investigator fault, the Ponce faults, the Guayanilla rift, and the Punta Montalva fault. The majority of the 2019–2021 Southwestern Puerto Rico earthquakes are associated with the Punta Montalva fault. Earthquakes forming 17° northward-dipping structures at various depths possibly manifest continuation of the Muertos trough, along which the Caribbean plate is being subducted beneath the Puerto Rico microplate. Our results show complex fault geometries of a diffuse fault network, suggesting possible subduction process accommodated by faults within a low-velocity zone.

Lessons Learned from Applying Varying Coefficient Model to Controlled Simulation Datasets

The development of site- and path-specific (i.e., non-ergodic) ground motion models (GMMs) can drastically improve the accuracy of probabilistic seismic hazard analyses (PSHA). The Varying Coefficient Model (VCM) is a novel technique for developing non-ergodic GMMs, which puts epistemic uncertainty into spatially varying coefficients. The coefficients at nearby locations are correlated by placing a Gaussian process prior on them. The correlation structure is determined by the data, and later used to predict coefficients and their epistemic uncertainties at new locations. It is important to carefully verify the technique before its results can be accepted by the engineering community. In this study, we used a series of simulation-based controlled ground motion datasets from CyberShake to test a modified VCM technique, which partitions the epistemic uncertainty into spatially varying source, site and path terms. Because the simulation parameters (inputs) are known, it is straightforward to verify what is recovered by the VCM from CyberShake simulation. We find that the site effects in CyberShake datasets can be reliably recovered by the VCM. However, the densely-located self-similar events in CyberShake datasets lead to large correlation lengths, which violates the isotropic assumption underlying the method and prevents the VCM from capturing the genuine source effects. For path effects, cell-specific attenuation approaches fail to recover the anelastic attenuation pattern of the 3D velocity model, most likely due to inappropriate assumption of point sources and straight-line wave propagation. Instead, a midpoint approach that only considers the aggregated path effects can better recover the strong attenuation within basins by fixing the correlation length of path effects. Lessons learned in this study not only provide important guidance for future applications of VCM to both simulation and empirical datasets, but also help further development of the technique, notably for the recovery of path effects.

An Investigation of Seismic Velocity Variation through a Tectonic Boundaries-Case Study in Central Iraq

A 3D velocity model was created by using stacking velocity of 9 seismic lines and average velocity of 6 wells drilled in Iraq. The model was achieved by creating a time model to 25 surfaces with an interval time between each two successive surfaces of about 100 msec. The summation time of all surfaces reached about 2400 msec, that was adopted according to West Kifl-1 well, which penetrated to a depth of 6000 m, representing the deepest well in the study area. The seismic lines and well data were converted to build a 3D cube time model and the velocity was spread on the model. The seismic inversion modeling of the elastic properties of the horizon and well data was applied to achieve a corrected velocity cube. Then, the velocity cube was converted to a time model and, finally, a corrected 3D depth model was obtained. This model shows that the western side of the study area, which is a part of the stable shelf, is characterized by relatively low thickness and high velocity layers. While the eastern side of the study area, which is a part of the Mesopotamian, is characterized by high thickness and low velocity of the Cretaceous succession. The Abu Jir fault is considered as a boundary between the stable and unstable shelves in Iraq, situated at the extreme west part of the study area. The area of relatively high velocity gradient is considered as the limit of the western side of the Mesopotamian basin. This area extends from Najaf-Karbala axis in the west to the Euphrates River in the east. It is found that the 3D stacking velocity model can be used to obtain good results concerning the tectonic boundary.

Interplay of seismic and a-seismic deformation during the 2020 sequence of Atacama, Chile.

An earthquake sequence occurred in the Atacama region of Chile throughout September 2020. The sequence initiated by a mainshock of magnitude Mw6.9, followed 17 hours later by a Mw6.4 aftershock. The sequence lasted several weeks, during which more than a thousand events larger than Ml 1 occurred, including several larger earthquakes of magnitudes between 5.5 and 6.4. Using a dense network that includes broad-band, strong motion and GPS sites, we study in details the seismic sources of the mainshock and its largest aftershock, the afterslip they generate and their aftershock, shedding light of the spatial temporal evolution of seismic and aseismic slip during the sequence. Dynamic inversions show that the two largest earthquakes are located on the subduction interface and have a stress drop and rupture times which are characteristics of subduction earthquakes. The mainshock and the aftershocks, localised in a 3D velocity model, occur in a narrow region of interseismic coupling (ranging 40%-80%), i.e. between two large highly coupled areas, North and South of the sequence, both ruptured by the great Mw~8.5 1922 megathrust earthquake. High rate GPS data (1 Hz) allow to determine instantaneous coseismic displacements and to infer coseismic slip models, not contaminated by early afterslip. We find that the total slip over 24 hours inferred from precise daily solutions is larger than the sum of the two instantaneous coseismic slip models. Differencing the two models indicates that rapid aseismic slip developed up-dip the mainshock rupture area and down-dip of the largest aftershock. During the 17 hours separating the two earthquakes, micro-seismicity migrated from the mainshock rupture area up-dip towards the epicenter of the Mw6.4 aftershocks and continued to propagate upwards at ~0.7 km/day. The bulk of the afterslip is located up-dip the mainshock and down-dip the largest aftershock, and is accompanied with the migration of seismicity, from the mainshock rupture to the aftershock area, suggesting that this aseismic slip triggered the Mw6.4 aftershock. Unusually large post-seismic slip, equivalent to Mw6.8 developed during three weeks to the North, in low coupling areas located both up-dip and downdip the narrow strip of higher coupling, and possibly connecting to the area of the deep Slow Sleep Event detected in the Copiapo area in 2014. The sequence highlights how seismic and aseismic slip interacted and witness short scale lateral variations of friction properties at the megathrust.

A new 3D velocity model of Central Apennines: insights from 2D/3D geological modelling and earthquake relocation.

<p>The Central Apennines fold-and-thrust belt (Central Italy) is characterized by the presence of several active faults, potentially capable of generating damaging earthquakes. To support seismic hazard studies over the area, a new 3D velocity model was built, integrating a wide range of surface and subsurface data.</p><p>The tectonic framework of the area (from Sulmona plain to Maiella Mt), is still debated in literature, also due to the lack of both an adequate geophysical data set and a reliable velocity model at the crustal scale.</p><p>In addition, the low number of seismic stations available for the acquisition of Vp/Vs arrival times, and the very low seismicity detected in the study area (the Sulmona and Caramanico Apennine valleys are considered as “seismic gaps”), lead to a difficult interpretation of the subsurface tectonic structures.</p><p>3D velocity modelling could well represent an important tool to support these deep crustal reconstructions as well earthquake relocation studies and could enhance the definition of seismogenic faults deep geometries, hence supporting a better risk assessment over the area of these potential locked faults.</p><p>Using the knowledge developed within the oil&gas industry as well in gas/CO<sub>2</sub> storage projects for the construction of 3D velocity models, extensively used to obtain subsurface imaging and define the geometry of the reservoirs and traps in the depth domain, a similar methodological approach was implemented over the study area.</p><p>The subsurface dataset was partially inherited by the past hydrocarbon exploration activities (e.g. seismic lines, exploration wells and sonic logs) and by the literature (e.g. time/depth regional models). Tomographic sections and relocated earthquake hypocentres were also integrated form geophysical studies. Geological maps (1:50.000 & 1:100.000 scale) represent the surface dataset that we used to create the surface interpretation of the regional geology.</p><p>As a first step, 18 2D balanced regional geological cross-sections, dip-oriented (W-E) across the Central Apennine, were built define the structural picture at regional scale. The cross-sections were built using MOVE (Petroleum Experts) and Petrel (Schlumberger) software. The following modelling step was the 3D model construction, in which the surface/subsurface data as well as all the geological sections were integrated in the final 3D structural and geological model.</p><p>The main geological layers reconstructed in the 3D model were than populated using the appropriated interval velocity values, building the final 3D velocity model in which the lateral velocity variation due to the presence of different facies/geological domains were considered.</p><p>As one of the results, we defined several 1D-velocity models coherent with the regional 3D velocity model, in which the key seismic stations and the earthquakes hypocentres dataset for the most potential seismogenic faults were included. 1D models were characterized by different degree of simplification, in order to test diverse approaches for the earthquake relocation. For this exercise, we used public dataset extracted by the analysis of microseismicity of the Sulmona basin.</p><p>We believe that the proposed approach can represents an effective method for combining geological and geophysical data to improve the subsurface and seismogenic faults interpretation, contributing to the seismic hazard assessment.</p>

Diffraction imaging and depth-velocity inversion with 3D P-Cable seismic data

<p>Most established methods for the estimation of subsurface velocity models rely on the measurements of reflected or diving waves and therefore require data with sufficiently large source-receiver offsets. For seismic data that lacks these offsets, such as vintage data, low-fold academic data or near zero-offset P-Cable data, these methods fail. Building on recent studies, we apply a workflow that exploits the diffracted wavefield for depth-velocity-model building. This workflow consists of three principal steps: (1) revealing the diffracted wavefield by modeling and adaptively subtracting reflections from the raw data, (2) characterizing the diffractions with physically meaningful wavefront attributes, (3) estimating depth-velocity models with wavefront tomography. We propose a hybrid 2D/3D approach, in which we apply the well-established and automated 2D workflow to numerous inlines of a high-resolution 3D P-Cable dataset acquired near Ritter Island, a small volcanic island located north-east of New Guinea known for a catastrophic flank collapse in 1888. We use the obtained set of parallel 2D velocity models to interpolate a 3D velocity model for the whole data cube, thus overcoming possible issues such as varying data quality in inline and crossline direction and the high computational cost of 3D data analysis. Even though the 2D workflow may suffer from out-of-plane effects, we obtain a smooth 3D velocity model that is consistent with the data.</p>

Ground motion simulations for finite-fault earthquake scenarios on the Húsavík-Flatey Fault, North Iceland

<p>The Tjörnes Fracture Zone (TFZ) in North Iceland is the largest and most complex zone of transform faulting in Iceland, formed due to a ridge-jump between two spreading centers of the Mid-Atlantic Ridge, the Northern Volcanic Zone and Kolbeinsey Ridge in North Iceland. Strong earthquakes (Ms>6) have repeatedly occurred in the TFZ and affected the North Icelandic population. In particular the large historical earthquakes of 1755 (Ms 7.0) and 1872 (doublet, Ms 6.5), have been associated with the Húsavı́k-Flatey Fault (HFF), which is the largest linear strike-slip transform fault in the TFZ, and in Iceland. We simulate fault rupture on the HFF and the corresponding near-fault ground motion for several potential earthquake scenarios, including scenario events that replicate the large 1755 and 1872 events. Such simulations are relevant for the town of Húsavı́k in particular, as it is located on top of the HFF and is therefore subject to the highest seismic hazard in the country. Due to the mostly offshore location of the HFF, its precise geometry has only recently been studied in more detail. We compile updated seismological and geophysical information in the area, such as a recently derived three-dimensional velocity model for P and S waves. Seismicity relocations using this velocity model, together with bathymetric and geodetic data, provide detailed information to constrain the fault geometry. In addition, we use this 3D velocity model to simulate seismic wave propagation. For this purpose, we generate a variety of kinematic earthquake-rupture scenarios, and apply a 3D finite-difference method (SORD) to propagate the radiated seismic waves through Earth structure. Slip distributions for the different scenarios are computed using a von Karman autocorrelation function whose parameters are calibrated with slip distributions available for a few recent Icelandic earthquakes. Simulated scenarios provide synthetic ground motion and time histories and estimates of peak ground motion parameters (PGA and PGV) at low frequencies (<2 Hz) for Húsavík and other main towns in North Iceland along with maps of ground shaking for the entire region [130 km x 110 km]. Ground motion estimates are compared with those provided by empirical ground motion models calibrated to Icelandic earthquakes and dynamic fault-rupture simulations for the HFF. Directivity effects towards or away from the coastal areas are analyzed to estimate the expected range of shaking. Thick sedimentary deposits (up to ∼4 km thick) located offshore on top of the HFF (reported by seismic, gravity anomaly and tomographic studies) may affect the effective depth of the fault's top boundary and the surface rupture potential. The results of this study showcase the extent of expected ground motions from significant and likely earthquake scenarios on the HFF. Finite fault earthquake simulations complement the currently available information on seismic hazard for North Iceland, and are a first step towards a systematic and large-scale earthquake scenario database on the HFF, and for the entire fault system of the TFZ, that will enable comprehensive and physics-based hazard assessment in the region.</p>