Frédérique Rolandone

Coseismic and postseismic slip of the 2004 Parkfield earthquake from space-geodetic data

We invert interferometric synthetic aperture radar (insar) data jointly with campaign and continuous global positioning system (gps) data for slip in the coseismic and postseismic periods of the 2004 Parkfield earthquake. The insar dataset consists of eight interferograms from data collected by the Envisat and Radarsat satellites spanning the time of the earthquake and variable amounts of the postseismic period. The two datasets complement each other, with the insar providing dense sampling of motion in the range direction of the satellite and the gps providing more sparse, but three-dimensional measurements of ground motion. The model assumes exponential decay of the postseismic slip with a decay time constant of 0.087 years, determined from time series modeling of continuous gps and creepmeter data. We find a geodetic moment magnitude of M 6.2 for a 1-day coseismic model and M w 6.1 for the entire postseismic period. The coseismic rupture occurred mainly in two slip asperities; one near the hypocenter and the other 15–20 km north. Postseismic slip occurred on the shallow portions of the fault and near the rupture areas of two M 5.0 aftershocks. A comparison of the geodetic slip models with seismic moment estimates suggests that the coseismic moment release of the Parkfield earthquake is as little as 25% of the total. This underlines the importance of aseismic slip in the slip budget for the Parkfield segment. Online material: Complete data tables and supplemental tables.

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.

Geologic versus geodetic deformation adjacent to the San Andreas fault, central California

We combine geologic and global positioning system (GPS) data to characterize the style and magnitude of off-fault deformation across the San Andreas fault system in central California. Geologic structures record ∼12 km of both fault-parallel and fault-perpendicular displacements across creeping and locked portions of the San Andreas fault. Analysis of 150 GPS site velocities suggests that the borderlands record 4–6 mm/yr of fault-parallel and 3–5 mm/yr of fault-perpendicular motion alongside the creeping segment, where elastic strain is minimized. The distribution of both long-term geologic and short-term geodetic deformation is affected by basement type, where more deformation is concentrated northeast of the San Andreas fault on Franciscan basement. We suggest that at least half the fault-parallel GPS deformation measured by GPS bordering the creeping segment must be accommodated by geologic structures; this permanent deformation needs to be incorporated into dynamic models of the fault system. Elastic modeling of the San Andreas fault in central California, which incorporates its well-known transition from locked to creeping behavior near Parkfield, predicts first-order variations in the GPS velocity field along the fault and corresponding variations in dilatational strain rates. The strain rate pattern is dominated by a large contractional region northeast of the transition from locked to creeping behavior and a large extensional region southwest of the transition. The former coincides with the Coalinga and Kettleman Hills anticlines, the growth and development of which seem to have occurred under at least two kinematic conditions. We suggest that the onset of fault creep in central California promoted the growth of these folds. By implication, fault creep has been active over geologic time scales.

Contrasted styles of rifting in the eastern Gulf of Aden: A combined wide-angle, multichannel seismic, and heat flow survey

Continental rifts and passive continental margins show fundamental along-axis segmentation patterns that have been attributed to one or a number of different processes: extensional fault geometry, variable stretching along strike, preexisting lithospheric compositional and structural heterogeneities, oblique rifting, and the presence or absence of eruptive volcanic centers. The length and width scales of the rift stage fault-bounded basin systems change during the late evolution of the new plate boundary, and the role of magmatism may increase as rifting progresses to continental rupture. Along obliquely spreading ridges, first-order mid-ocean ridge geometries originate during the synrift stage, indicating an intimate relationship between magma production and transform fault spacing and location. The Gulf of Aden rift is a young ocean basin in which the earliest synrift to breakup structures are well exposed onshore and covered by thin sediment layers offshore. This obliquely spreading rift is considered magma-poor and has several large-offset transforms that originated during late stage rifting and control the first-order axial segmentation of the spreading ridge. Widely spaced geophysical transects of passive margins that produce only isolated 2-D images of crust and uppermost mantle structure are inadequate for evaluation of competing rift evolution models. Using closely spaced new geophysical and geological observations from the Gulf of Aden we show that rift sectors between transforms have a large internal variability over short distances (∼10 km): the ocean-continent transition (OCT) evolves from a narrow magmatic transition to wider zones where continental mantle is probably exhumed. We suggest that this small-scale variability may be explained (1) by the distribution of volcanism and (2) by the along-strike differences in time-averaged extension rate of the oblique rift system. The volcanism may be associated with (1) the long-offset Alula-Fartak Fracture Zone, which may enhance magma production on its younger side, or (2) channeled flow from the Afar plume material along the newly formed OCT and the spreading ridge. Oblique extension and/or hot spot interactions may thereby have a significant control on the styles of rifting and continental breakup and on the evolution of many magma-poor margins.

Crustal and upper mantle structure beneath south-western margin of the Arabian Peninsula from teleseismic tomography

We image the lithospheric and upper asthenospheric structure of western continental Yemen with 24 broadband stations to evaluate the role of the Afar plume on the evolution of the continental margin and its extent eastward along the Gulf of Aden. We use teleseismic tomography to compute relative P wave velocity variations in south-western Yemen down to 300 km depth. Published receiver function analysis suggest a dramatic and localized thinning of the crust in the vicinity of the Red Sea and the Gulf of Aden, consistent with the velocity structure that we retrieve in our model. The mantle part of the model is dominated by the presence of a low-velocity anomaly in which we infer partial melting just below thick Oligocene flood basalts and recent off-axis volcanic events (from 15 Ma to present). This low-velocity anomaly could correspond to an abnormally hot mantle and could be responsible for dynamic topography and recent magmatism in western Yemen. Our new P wave velocity model beneath western Yemen suggests the young rift flank volcanoes beneath margins and on the flanks of the Red Sea rift are caused by focused small-scale diapiric upwelling from a broad region of hot mantle beneath the area. Our work shows that relatively hot mantle, along with partial melting of the mantle, can persist beneath rifted margins after breakup has occurred.

Partitioning of oblique convergence in the Northern Andes subduction zone: Migration history and the present-day boundary of the North Andean Sliver in Ecuador

Along the Ecuadorian margin, oblique subduction induces deformation of the overriding continental plate. For the last 15 Ma, both exhumation and tectonic history of Ecuador suggest that the northeastward motion of the North Andean Sliver (NAS) was accompanied by an eastward migration of its eastern boundary and successive progressively narrowing restraining bends. Here we present geologic data, earthquake epicenters, focal mechanisms, GPS results, and a revised active fault map consistent with this new kinematic model. All data sets concur to demonstrate that active continental deformation is presently localized along a single major fault system, connecting fault segments from the Gulf of Guayaquil to the eastern Andean Cordillera. Although secondary faults are recognized within the Cordillera, they accommodate a negligible fraction of relative motion compared to the main fault system. The eastern limit is then concentrated rather than distributed as first proposed, marking a sharp boundary between the NAS, the Inca sliver, and the Subandean domain overthrusting the South American craton. The NAS limit follows a northeast striking right-lateral transpressional strike-slip system from the Gulf of Guayaquil (Isla Puna) to the Andean Cordillera and with the north-south striking transpressive faults along the eastern Andes. Eastward migration of the restraining belt since the Pliocene, abandonment of the sutures and reactivation of north-south striking ancient fault zones lead to the final development of a major tectonic boundary south and east of the NAS, favoring its extrusion as a continental sliver, accommodating the oblique convergence of the Nazca oceanic plate toward South America.

Active tectonics in Quito, Ecuador, assessed by geomorphological studies, GPS data, and crustal seismicity

The Quito Fault System (QFS) extends over 60 km along the Interandean Depression in northern Ecuador. Multidisciplinary studies support an interpretation in which two major contemporaneous fault systems affect Quaternary volcanoclastic deposits. Hanging paleovalleys and disruption of drainage networks attest to ongoing crustal deformation and uplift in this region, further confirmed by 15 years of GPS measurements and seismicity. The resulting new kinematic model emphasizes the role of the N-S segmented, en echelon eastward migrating Quito Fault System (QFS). Northeast of this major tectonic feature, the strike-slip Guayllabamba Fault System (GFS) aids the eastward transfer of the regional strain toward Colombia. These two tectonic fault systems are active, and the local focal mechanisms are consistent with the direction of relative GPS velocities and the regional stress tensor. Among active features, inherited N-S direction sutures appear to play a role in confining the active deformation in the Interandean Depression. The most frontal of the Quito faults formed at the tip of a blind thrust, dipping 40°W, is most probably connected at depth to inactive suture to the west. A new GPS data set indicates active shortening rates for Quito blind thrust of up to 4 mm/yr, which decreases northward along the fold system as it connects to the strike-slip Guayllabamba Fault System. The proximity of these structures to the densely populated Quito region highlights the need for additional tectonic studies in these regions of Ecuador to generate further hazard assessments.

Uppermost mantle velocity from Pn tomography in the Gulf of Aden

We determine the lateral variations in seismic velocity of the lithospheric mantle beneath the Gulf of Aden and its margins by inversion of Pn (upper mantle high-frequency compressional P wave) traveltimes. Data for this study were collected by several temporary seismic networks and from the global catalogue. A least-squares tomographic algorithm is used to solve for velocity variations in the mantle lithosphere. In order to separate shallow and deeper structures, we use separate inversions for shorter and longer ray path data. High Pn velocities (8.2–8.4 km/s) are observed in the uppermost mantle beneath Yemen that may be related to the presence of magmatic underplating of the volcanic margins of Aden and the Red Sea. Zones of low velocity (7.7 km/s) are present in the shallow upper mantle beneath Sana’a, Aden, Afar, and along the Gulf of Aden that are likely related to melt transport through the lithosphere feeding active volcanism. Deeper within the upper mantle, beneath the Oman margin, a low-velocity zone (7.8 km/s) suggests a deep zone of melt accumulation. Our results provide evidence that the asthenosphere undergoes channelized flow from the Afar hotspot toward the east along the Aden and Sheba Ridges.

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.

Magmatism at continental passive margins inferred from Ambient-Noise Phase-velocity in the Gulf of Aden

Non-volcanic continental passive margins have traditionally been considered to be tectonically and magmatically inactive once continental breakup has occurred and seafloor spreading has commenced. We use ambient-noise tomography to constrain Rayleigh-wave phase-velocity maps beneath the eastern Gulf of Aden (eastern Yemen and southern Oman). In the crust, we image low velocities beneath the Jiza-Qamar (Yemen) and Ashawq-Salalah (Oman) basins, likely caused by the presence of partial melt associated with magmatic plumbing systems beneath the rifted margin. Our results provide strong evidence that magma intrusion persists after breakup, modifying the composition and thermal structure of the continental margin. The coincidence between zones of crustal intrusion and steep gradients in lithospheric thinning, as well as with transform faults, suggests that magmatism post-breakup may be driven by small-scale convection and enhanced by edge-driven flow at the juxtaposition of lithosphere of varying thickness and thermal age.