Edmundo Norabuena

Decelerating Nazca-South America and Nazca-Pacific Plate motions

Space geodetic estimates of the rate of Nazca-South America convergence and Nazca-Pacific spreading averaging over several years show that present day rates are significantly slower than the 3 million year average NUVEL-1A model. The implied rates of deceleration are consistent with longer term trends extending back to at least 20 Ma, about the time of initiation of Andes growth, and may reflect consequences of ongoing subduction and construction of the Andes, e.g., increased friction and viscous drag on the subducted slab as the leading edge of South America thickens.

Paleoseismology and Global Positioning System: Earthquake-cycle effects and geodetic versus geologic fault slip rates in the Eastern California shear zone

Published slip rates for the Owens Valley fault zone in eastern California based on geodetic data and elastic half-space models (5–7 mm/yr) are faster than longer term geologic rates (2–3 mm/yr). We use Global Positioning System data spanning the central Owens Valley, a more realistic rheological model with an elastic upper crust over a viscoelastic lower crust and upper mantle, and paleoseismic data from adjacent faults, to show that this difference could reflect earthquake-cycle effects. We estimate a long-term rate (2.1 ± 0.7 mm/yr) and earthquake recurrence interval (2750 +350/−1000 yr) from the geodetic data, both in agreement with independent geologic estimates.

Geodetic and seismic constraints on some seismogenic zone processes in Costa Rica

New seismic and geodetic data from Costa Rica provide insight into seismogenic zone processes in Central America, where the Cocos and Caribbean plates converge. Seismic data are from combined land and ocean bottom deployments in the Nicoya peninsula in northern Costa Rica and near the Osa peninsula in southern Costa Rica. In Nicoya, inversion of GPS data suggests two locked patches centered at 14 ± 2 and 39 ± 6 km depth. Interplate microseismicity is concentrated in the more freely slipping intermediate zone, suggesting that small interseismic earthquakes may not accurately outline the updip limit of the seismogenic zone, the rupture zone for future large earthquakes, at least over the short (∼1 year) observation period. We also estimate northwest motion of a coastal “sliver block” at 8 ± 3 mm/yr, probably related to oblique convergence. In the Osa region to the south, convergence is orthogonal to the trench. Cocos-Caribbean relative motion is partitioned here, with ∼8 cm/yr on the Cocos-Panama block boundary (including a component of permanent shortening across the Fila Costena fold and thrust belt) and ∼1 cm/yr on the Panama block–Caribbean boundary. The GPS data suggest that the Cocos plate–Panama block boundary is completely locked from ∼10–50 km depth. This large locked zone, as well as associated forearc and back-arc deformation, may be related to subduction of the shallow Cocos Ridge and/or younger lithosphere compared to Nicoya, with consequent higher coupling and compressive stress in the direction of plate convergence.

Consistency of geologic and geodetic displacements during Andean orogenesis

Displacement and convergence rates inferred from GPS and marine magnetic data suggest that this trend may be continuing at present. Hence an increasing fraction of convergence is being absorbed by mountain building. This change may reflect increased plate coupling due to decreased sediment supply, younger subducting lithosphere, or increased normal stress at the interface from the effects of uplift. INDEX TERMS: 8102 Tectonophysics: Continental contract orogenic belts; 8157 Evolution of the Earth: Plate motions—past (3040); 8158 Evolution of the Earth: Plate motions—present and recent (3040); 8099 Structural Geology: General or miscellaneous

Geodetic, teleseismic, and strong motion constraints on slip from recent southern Peru subduction zone earthquakes

We use seismic and geodetic data both jointly and separately to constrain coseismic slip from the 12 November 1996 M_w 7.7 and 23 June 2001 M_w 8.5 southern Peru subduction zone earthquakes, as well as two large aftershocks following the 2001 earthquake on 26 June and 7 July 2001. We use all available data in our inversions: GPS, interferometric synthetic aperture radar (InSAR) from the ERS-1, ERS-2, JERS, and RADARSAT-1 satellites, and seismic data from teleseismic and strong motion stations. Our two-dimensional slip models derived from only teleseismic body waves from South American subduction zone earthquakes with M_w > 7.5 do not reliably predict available geodetic data. In particular, we find significant differences in the distribution of slip for the 2001 earthquake from models that use only seismic (teleseismic and two strong motion stations) or geodetic (InSAR and GPS) data. The differences might be related to postseismic deformation or, more likely, the different sensitivities of the teleseismic and geodetic data to coseismic rupture properties. The earthquakes studied here follow the pattern of earthquake directivity along the coast of western South America, north of 5°S, earthquakes rupture to the north; south of about 12°S, directivity is southerly; and in between, earthquakes are bilateral. The predicted deformation at the Arequipa GPS station from the seismic-only slip model for the 7 July 2001 aftershock is not consistent with significant preseismic motion.

An 8 month slow slip event triggers progressive nucleation of the 2014 Chile megathrust

The mechanisms leading to large earthquakes are poorly understood and documented. Here, we characterize the long-term precursory phase of the April 1st 2014 Mw8.1 North Chile megathrust. We show that a group of coastal GPS stations accelerated westward eight months before the mainshock, corresponding to a Mw6.5 slow slip event on the subduction interface, 80% of which was aseismic. Concurrent interface foreshocks underwent a diminution of their radiation at high frequency, as shown by the temporal evolution of Fourier spectra and residuals with respect to ground motions predicted by recent subduction models. Such ground-motions change suggests that, in response to the slow sliding of the subduction interface, seismic ruptures are progressively becoming smoother and/or slower. The gradual propagation of seismic ruptures beyond seismic asperities into surrounding metastable areas could explain these observations, and might be the precursory mechanism eventually leading to the mainshock.

Structure of the subducting Nazca Plate beneath Peru

Arrival times from intermediate-depth (110–150 km) earthquakes within the region of flat subduction beneath the subandean zone and foreland basins of east-central Peru provide constraints on the geometry and velocity structure of the subducting Nazca plate. Hypocentral locations and origin times for these events were determined using observations from a 15 station digitally recording locator array deployed in the epicentral region of eastern Peru. Observed P wave arrival times for coastal stations in Peru, some 3–6° from the epicenters, are up to 4 s early relative to predicted arrival times based on the best fit velocity-depth model used for hypocenter locations. These large negative time residuals appear to be the result of propagation paths which have long segments in the colder, higher-velocity subducting plate, P wave travel times were modeled for the effects of the slab using three-dimensional (3-D) ray tracing. Computed ray paths show that travel times to coastal stations for the eastern Peru events can be satisfactorily modeled with average velocities relative to the surrounding mantle 6% lower within the uppermost slab (assumed on the basis of other studies to be unconverted basaltic oceanic crust 6 km thick) and 8% higher within the cold uppermost mantle of the slab. Ray tracing for this plate model shows that P wave “shadow zones” can occur if the source-slab-receiver geometry results in seismic rays passing through regions in which the dip angle of the slab changes significantly. Such geometries exist for seismic waves propagating to some coastal stations from intermediate-depth earthquakes located east of the Andes. Observed first-arrival times for such cases do in fact have smaller negative residuals than those for geometries which allow for “direct” paths at similar distances. Modeling such arrivals as internally reflected waves propagating through the high-velocity part of the plate produces a significant improvement in the travel time residuals. For the slab velocities given above, we obtain a model thickness of approximately 36 km for the cold slab interior and a slight northwest component of dip in the region of subhorizontal subduction.

Tectonic inversion in the Caribbean‐South American plate boundary: GPS geodesy, seismology, and tectonics of the Mw 6.7 22 April 1997 Tobago earthquake

On 22 April 1997 the largest earthquake recorded in the Trinidad-Tobago segment of the Caribbean-South American plate boundary zone (Mw 6.7) ruptured a shallow (~9 km), ENE striking (~250° azimuth), shallowly dipping (~28°) dextral-normal fault ~10 km south of Tobago. In this study, we describe this earthquake and related foreshock and aftershock seismicity, derive coseismic offsets using GPS data, and model the fault plane and magnitude of slip for this earthquake. Coseismic slip estimated at our episodic GPS sites indicates movement of Tobago 135 ± 6 to 68 ± 6 mm NNE and subsidence of 7 ± 9 to 0 mm. This earthquake was anomalous and is of interest because (1) its large component of normal slip and ENE strike are unexpected given the active E-W dextral shearing across the Caribbean-South American plate boundary zone, (2) it ruptured a normal fault plane with a low (~28°) dip angle, and (3) it reactivated and inverted the preexisting Tobago terrrane-South America ocean-continent (thrust) boundary that formed during early Tertiary oblique plate convergence.