Spatiotemporal evolution of surface creep in the Parkfield region of the San Andreas Fault (1993-2004) from synthetic aperture radar.

International audienceThe Parkfield section of the San Andreas Fault (SAF) is defined as a transitional portion of the fault between slip-release behavior types in the creeping section of the SAF to the northwest and the apparently locked section to the southeast. The Parkfield section is characterized by complex frictional fault behavior because it represents a transition zone from aseismic creep to stick-slip regime. At least six historic earthquakes of Mw ~6 have occurred in this area in 1881, 1901, 1922, 1934, 1966, and 2004. It was observed in the 2004 Mw 6.0 Parkfield earthquake that ~70% of the total (coseismic and postseismic) moment release occurred aseismically. To understand the SAF behavior in this area, it is of particular interest to measure and analyze, not only the spatial evolution of the surface displacement in this area, but also its evolution over time. Using radar data acquired by the European Space Agency's European Remote Sensing (ERS1-2) satellites, we constructed descending interferograms and retrieved time series of surface displacements along the central SAF for the decade preceding the 2004 Parkfield earthquake. We focus on characterizing the space and time evolution of surface creep in the Parkfield and Cholame sections. The spatial pattern of the interseismic displacement rate indicates that tectonic strain was not uniformly distributed along the strike of the fault between 1993 and 2004. Our data indicate not only a decrease in the creep rate from the Parkfield section to south of Highway-46 from 1.4 ±0.3 cm/y to 0.6 ±0.3 cm/y, but also a small but significant creep-rate increase in the Cholame section to 0.2 ±0.1 cm/y. The evidence for episodic creep in the Cholame section of the SAF south-east of Parkfield is in contrast with previously published interpretations of GPS and trilateration data. The Cholame section of the SAF merits close monitoring because it was likely the nucleation site of the 1857 Fort Tejón earthquake and because it has shown recent evidence of deep slow slip as revealed by deep tremors

Triggered crustal earthquake swarm across subduction segment boundary after the 2016 Pedernales, Ecuador megathrust earthquake

Megathrust ruptures and the ensuing postseismic deformation cause stress changes that may induce seismicity on upper plate crustal faults far from the coseismic rupture area. In this study, we analyze seismic swarms that occurred in the north Ecuador area of Esmeraldas, beginning two months after the 2016 M$_{w}$ 7.8 Pedernales, Ecuador megathrust earthquake. The Esmeraldas region is 70 km from the Pedernales rupture area in a separate segment of the subduction zone. We characterize the Esmeraldas sequence, relocating the events using manual arrival time picks and a local a-priori 3D velocity model. The earthquake locations from the Esmeraldas sequence outline an upper plate fault or shear zone. The sequence contains one major swarm and several smaller swarms. Moment tensor solutions of several events include normal and strike-slip motion and non-double-couple components. During the main swarm, earthquake hypocenters increase in distance from the first event over time, at a rate of a few hundred meters per day, consistent with fluid diffusion. Events with similar waveforms occur within the sequence, and a transient is seen in time series of nearby GPS stations concurrent with the seismicity. The events with similar waveforms and the transient in GPS time series suggest that slow aseismic slip took place along a crustal normal fault during the sequence. Coulomb stress calculations show a positive Coulomb stress change in the Esmeraldas region, consistent with seismicity being triggered by the Pedernales mainshock and large aftershocks. The characteristics of the seismicity indicate that postseismic deformation involving fluid flow and slow slip activated upper plate faults in the Esmeraldas area. These findings suggest the need for further investigation into the seismic hazard potential of shallow upper plate faults and the potential for megathrust earthquakes to trigger slow-slip and shallow seismicity across separate segments of subduction zones

Structural Control on Megathrust Rupture and Slip Behavior: Insights From the 2016 Mw 7.8 Pedernales Ecuador Earthquake

The heterogeneous seafloor topography of the Nazca Plate as it enters the Ecuador subduction zone provides an opportunity to document the influence of seafloor roughness on slip behavior and megathrust rupture. The 2016 M$_{w}$ 7.8 Pedernales Ecuador earthquake was followed by a rich and active postseismic sequence. An internationally coordinated rapid response effort installed a temporary seismic network to densify coastal stations of the permanent Ecuadorian national seismic network. A combination of 82 onshore short and intermediate period and broadband seismic stations and six ocean bottom seismometers recorded the postseismic Pedernales sequence for over a year after the mainshock. A robust earthquake catalog combined with calibrated relocations for a subset of magnitude ≥4 earthquakes shows pronounced spatial and temporal clustering. A range of slip behavior accommodates postseismic deformation including earthquakes, slow slip events, and earthquake swarms. Models of plate coupling and the consistency of earthquake clustering and slip behavior through multiple seismic cycles reveal a segmented subduction zone primarily controlled by subducted seafloor topography, accreted terranes, and inherited structure. The 2016 Pedernales mainshock triggered moderate to strong earthquakes (5 ≤ M ≤ 7) and earthquake swarms north of the mainshock rupture close to the epicenter of the 1906M$_{w}$ 8.8 earthquake and in the segment of the subduction zone that ruptured in 1958 in a M$_{w}$ 7.7 earthquake

When do fractured media become seismically anisotropic? Some implications on quantifying fracture properties

Fractures are pervasive features within the Earth's crust and they have a significant influence on the multi-physical response of the subsurface. The presence of coherent fracture sets often leads to observable seismic anisotropy enabling seismic techniques to remotely locate and characterise fracture systems. In this study, we confirm the general scale-dependence of seismic anisotropy and provide new results specific to shear-wave splitting (SWS). We find that SWS develops under conditions when the ratio of wavelength to fracture size (λS/d) is greater than 3, where Rayleigh scattering from coherent fractures leads to an effective anisotropy such that effective medium model (EMM) theory is qualitatively valid. When 1<λS/d<3 there is a transition from Rayleigh to Mie scattering, where no effective anisotropy develops and hence the SWS measurements are unstable. When λS/d<1 we observe geometric scattering and begin to see behaviour similar to transverse isotropy. We find that seismic anisotropy is more sensitive to fracture density than fracture compliance ratio. More importantly, we observe that the transition from scattering to an effective anisotropic regime occurs over a propagation distance between 1 and 2 wavelengths depending on the fracture density and compliance ratio. The existence of a transition zone means that inversion of seismic anisotropy parameters based on EMM will be fundamentally biased. More importantly, we observe that linear slip EMM commonly used in inverting fracture properties is inconsistent with our results and leads to errors of approximately 400% in fracture spacing (equivalent to fracture density) and 60% in fracture compliance. Although EMM representations can yield reliable estimates of fracture orientation and spatial location, our results show that EMM representations will systematically fail in providing quantitatively accurate estimates of other physical fracture properties, such as fracture density and compliance. Thus more robust and accurate quantitative estimates of in situ fracture properties will require improvements to effective medium models as well as the incorporation of full-waveform inversion techniques

Interseismic coupling and seismic potential along the Central Andes subduction zone

We use about two decades of geodetic measurements to characterize interseismic strain build up along the Central Andes subduction zone from Lima, Peru, to Antofagasta, Chile. These measurements are modeled assuming a 3-plate model (Nazca, Andean sliver and South America Craton) and spatially varying interseismic coupling (ISC) on the Nazca megathrust interface. We also determine slip models of the 1996 M_w = 7.7 Nazca, the 2001 M_w = 8.4 Arequipa, the 2007 M_w = 8.0 Pisco and the M_w = 7.7 Tocopilla earthquakes. We find that the data require a highly heterogeneous ISC pattern and that, overall, areas with large seismic slip coincide with areas which remain locked in the interseismic period (with high ISC). Offshore Lima where the ISC is high, a M_w∼8.6–8.8 earthquake occurred in 1746. This area ruptured again in a sequence of four M_w∼8.0 earthquakes in 1940, 1966, 1974 and 2007 but these events released only a small fraction of the elastic strain which has built up since 1746 so that enough elastic strain might be available there to generate a M_w > 8.5 earthquake. The region where the Nazca ridge subducts appears to be mostly creeping aseismically in the interseismic period (low ISC) and seems to act as a permanent barrier as no large earthquake ruptured through it in the last 500 years. In southern Peru, ISC is relatively high and the deficit of moment accumulated since the M_w∼8.8 earthquake of 1868 is equivalent to a magnitude M_w∼8.4 earthquake. Two asperities separated by a subtle aseismic creeping patch are revealed there. This aseismic patch may arrest some rupture as happened during the 2001 Arequipa earthquake, but the larger earthquakes of 1604 and 1868 were able to rupture through it. In northern Chile, ISC is very high and the rupture of the 2007 Tocopilla earthquake has released only 4% of the elastic strain that has accumulated since 1877. The deficit of moment which has accumulated there is equivalent to a magnitude M_w∼8.7 earthquake. This study thus provides elements to assess the location, size and magnitude of future large megathurst earthquakes in the Central Andes subduction zone. Caveats of this study are that interseismic strain of the forearc is assumed time invariant and entirely elastic. Also a major source of uncertainty is due to fact that the available data place very little constraints on interseismic coupling at shallow depth near the trench, except offshore Lima where sea bottom geodetic measurements have been collected suggesting strong coupling

Impacts of the Last Glacial Cycle on ground surface temperature reconstructions over the last millennium

Borehole temperature profiles provide robust estimates of past ground surface temperature changes, in agreement with meteorological data. Nevertheless, past climatic changes such as the Last Glacial Cycle (LGC) generated thermal effects in the subsurface that affect estimates of recent climatic change from geothermal data. We use an ensemble of ice sheet simulations spanning the last 120 ka to assess the impact of the Laurentide Ice Sheet on recent ground surface temperature histories reconstructed from borehole temperature profiles over North America. When the thermal remnants of the LGC are removed, we find larger amounts of subsurface heat storage (2.8 times) and an increased warming of the ground surface over North America by 0.75 K, both relative to uncorrected borehole estimates

Segmentation and kinematics of the North America-Caribbean plate boundary offshore Hispaniola

We explored the submarine portions of the Enriquillo–Plantain Garden Fault zone (EPGFZ) and the Septentrional–Oriente Fault zone (SOFZ) along the Northern Caribbean plate boundary using high-resolution multibeam echo-sounding and shallow seismic reflection. The bathymetric data shed light on poorly documented or previously unknown submarine fault zones running over 200 km between Haiti and Jamaica (EPGFZ) and 300 km between the Dominican Republic and Cuba (SOFZ). The primary plate-boundary structures are a series of strike-slip fault segments associated with pressure ridges, restraining bends, step overs and dogleg offsets indicating very active tectonics. Several distinct segments 50–100 km long cut across pre-existing structures inherited from former tectonic regimes or bypass recent morphologies formed under the current strike-slip regime. Along the most recent trace of the SOFZ, we measured a strike-slip offset of 16.5 km, which indicates steady activity for the past ~1.8 Ma if its current GPS-derived motion of 9.8 ± 2 mm a−1 has remained stable during the entire Quaternary.Depto. de Geodinámica, Estratigrafía y PaleontologíaFac. de Ciencias GeológicasTRUEpu

Magmatism on rift flanks: insights from ambient noise phase velocity in Afar region

During the breakup of continents in magmatic settings, the extension of the rift valley is commonly assumed to initially occur by border faulting and progressively migrate in space and time toward the spreading axis. Magmatic processes near the rift flanks are commonly ignored. We present phase velocity maps of the crust and uppermost mantle of the conjugate margins of the southern Red Sea (Afar and Yemen) using ambient noise tomography to constrain crustal modification during breakup. Our images show that the low seismic velocities characterize not only the upper crust beneath the axial volcanic systems but also both upper and lower crust beneath the rift flanks where ongoing volcanism and hydrothermal activity occur at the surface. Magmatic modification of the crust beneath rift flanks likely occurs for a protracted period of time during the breakup process and may persist through to early seafloor spreading