Abstract 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 1906 M w 8.8 earthquake and in the segment of the subduction zone that ruptured in 1958 in a M w 7.7 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 Mw 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.
Abstract We present a crust and mantle 3‐D shear velocity model extending well offshore of New Zealand's South Island, imaging the lithosphere beneath the South Island as well as the Campbell and Challenger Plateaus. Our model is constructed via linearized inversion of both teleseismic (18–70 s period) and ambient noise‐based (8–25 s period) Rayleigh wave dispersion measurements. We augment an array of 4 land‐based and 29 ocean bottom instruments deployed off the South Island's east and west coasts in 2009–2010 by the Marine Observations of Anisotropy Near Aotearoa experiment with 28 land‐based seismometers from New Zealand's permanent GeoNet array. Major features of our shear wave velocity ( V s ) model include a low‐velocity ( V s < 4.4 km/s) body extending from near surface to greater than 75 km depth beneath the Banks and Otago Peninsulas and high‐velocity ( V s ~4.7 km/s) mantle anomalies underlying the Southern Alps and off the northwest coast of the South Island. Using the 4.5 km/s contour as a proxy for the lithosphere‐asthenosphere boundary, our model suggests that the lithospheric thickness of Challenger Plateau and central South Island is substantially greater than that of the inner Campbell Plateau. The high‐velocity anomaly we resolve at subcrustal depths (>50 km) beneath the central South Island exhibits strong spatial correlation with upper mantle earthquake hypocenters beneath the Alpine Fault. The ~400 km long low‐velocity zone we image beneath eastern South Island and the inner Bounty Trough underlies Cenozoic volcanics and the locations of mantle‐derived helium measurements, consistent with asthenospheric upwelling in the region.
ABSTRACT The west part of Java sits at the transition from oblique subduction of the Australian plate under the Sunda block of the Eurasian plate along Sumatra to orthogonal convergence along central and eastern Java. This region has experienced several destructive earthquakes, the 17 July 2006 Mw 7.7 earthquake and tsunami off the coast of Pangandaran and the 2 September 2009 Mw 7 earthquake, located off the coast of Tasikmalaya. More recently, on 15 December 2017, an Mw 6.5 earthquake occurred off the coast near Pangandaran, and, on 23 January 2018, an Mw 5.9 earthquake occurred offshore Lebak, between Pelabuhan Ratu and Ujung Kulon. Ground shaking and damage occurred locally and in Jakarta on the northern coast of Java. In this study, we use the double-difference technique to relocate both mainshocks and 10 months of seismicity (228 events) following the earthquakes. The relocation result improved the mainshock locations and depth distribution of earthquakes. Moment tensor of the December 2017 event located the hypocenter at ∼108 km depth within the subducting slab. The best-fit relocation places the depth at 61 km, close to the slab interface. Aftershocks occur between 68 and 86 km depth and align along a steeper plane than slab geometry models. The January 2018 event is located at ∼46 km depth. Aftershocks form a near-vertical, pipe-like structure from the plate interface to ∼10 km depth. A burst of aftershocks immediately following the mainshock shows a shallowing upward trend at a rate of ∼2 km/hr, suggesting that a fluid pressure wave released from the oceanic crust is causing brittle failure in the overriding plate, followed by upward migration of fluids. Five months later, shallow (<25 km) seismicity collocates with background seismicity, suggesting the January 2018 event activated the Pelabuhan Ratu fault system close to the coast.
Crust and mantle 3-D shear velocity model extending well offshore of New Zealand’s South Island, imaging the lithosphere beneath the South Island as well as the Campbell and Challenger Plateaus. Our model is constructed via linearized inversion of both teleseismic (18–70 s period) and ambient noise-based (8–25 s period) Rayleigh wave dispersion measurements. We augment an array of 4 land-based and 29 ocean bottom instruments deployed off the South Island’s east and west coasts in 2009–2010 by the Marine Observations of Anisotropy Near Aotearoa experiment with 28 land-based seismometers from New Zealand’s permanent GeoNet array. This model is completely unconstrained (all layers are allowed to vary equally). The "Moho" model using Moho depth constraints obtained from the model of Salmon et al. (2011) can be found under 'Shear Wave Velocity of New Zealand's South Island (Moho version)'.