Stéphanie Gautier

Segmentation of the Himalayas as revealed by arc-parallel gravity anomalies.

Lateral variations along the Himalayan arc are suggested by an increasing number of studies and carry important information about the orogen’s segmentation. Here we compile the hitherto most complete land gravity dataset in the region which enables the currently highest resolution plausible analysis. To study lateral variations in collisional structure we compute arc-parallel gravity anomalies (APaGA) by subtracting the average arc-perpendicular profile from our dataset; we compute likewise for topography (APaTA). We find no direct correlation between APaGA, APaTA and background seismicity, as suggested in oceanic subduction context. In the Himalayas APaTA mainly reflect relief and erosional effects, whereas APaGA reflect the deep structure of the orogen with clear lateral boundaries. Four segments are outlined and have disparate flexural geometry: NE India, Bhutan, Nepal & India until Dehradun, and NW India. The segment boundaries in the India plate are related to inherited structures, and the boundaries of the Shillong block are highlighted by seismic activity. We find that large earthquakes of the past millennium do not propagate across the segment boundaries defined by APaGA, therefore these seem to set limits for potential rupture of megathrust earthquakes.

Improvements of the Maximum Intersection Method for 3D Absolute Earthquake Locations

Wepresentanimprovementofthemaximumintersection(MAXI)method in which absolute earthquake location is defined by the maximum intersection number of hyperbolic equal differential time (EDT) volumes (one EDT volume being described as all grid nodes satisfying the arrival-time differences between two stations, a tolerance value known as TERR). This 3D technique is well adapted to a strongly heterogeneousenvironment,avoidsthedepthversusorigin-timetrade-off(evenusingP arrivals alone), and objectively filters possible erroneous arrival times. Improvements consist of introducing an iterative multiscale approach on the TERR parameter and conducting a hierarchy of grid discretization to refine the preliminary solution (multigrid algorithm). The outlier cleaning process is strengthened by conducting statistical analyses based on the number of EDT volume intersections per station, rather than on travel-time residues. Processing two series of synthetic arrival times computed in the 3D seismological configurationofthesouthernmostRyukyusubductionzone,offshoreTaiwan,allowsus to (1) confirm MAXI robustness, even when arrival times are highly disturbed or azimuthal coverage largely reduced, and (2) show that MAXI, by limiting trade-offs between origin time and depth and between epicentral position and depth, is less af- fected by a wrong velocity model than a technique that minimizes travel-time residues. To conclude, we apply the MAXI method to locate earthquakes recorded by the combined Central Weather Bureau of Taiwan (CWB) and the Japanese Meteorological Agency (JMA) networks (1991-2008) within a georealistic a priori 3D P-velocity model. Hypocenter determinations are selected based on MAXI confidence factors, and earthquake clustering along well-known tectonic features witnesses hypocenter determination accuracy.

Evidence of interseismic coupling variations along the Bhutan Himalayan arc from new GPS data

Although the first-order pattern of present-day deformation is relatively well resolved across the Himalayas, irregular data coverage limits detailed analyses of spatial variations of interseismic coupling. We provide the first GPS velocity field for the Bhutan Himalaya. Combined with published data, these observations show strong east-west variations in coupling between central and eastern Bhutan. In contrast with previous estimations of first-order uniform interseismic coupling along the Himalayan arc, we identify significant lateral variations: In western and central Bhutan, the fully coupled segment is 135–155 km wide with an abrupt downdip transition, whereas in eastern Bhutan the fully coupled segment is 100–120 km wide and is limited updip and downdip by partially creeping segments. This is the first observation of decoupling on the upper ramp along the Himalayan arc, with important implications for large earthquake surface rupture and seismic hazard.

Fault-magma interactions during early continental rifting: seismicity of the Magadi-Natron-Manyara basins, Africa

Although magmatism may occur during the earliest stages of continental rifting, its role in strain accommodation remains weakly constrained by largely 2D studies. We analyze seismicity data from a 13-month, 39-station broadband seismic array to determine the role of magma intrusion on state-of-stress and strain localization, and their along-strike variations. Precise earthquake locations using cluster analyses and a new 3D velocity model reveal lower crustal earthquakes beneath the central basins and along projections of steep border faults that degas CO2. Seismicity forms several disks interpreted as sills at 6-10 km below a monogenetic cone field. The sills overlie a lower crustal magma chamber that may feed eruptions at Oldoinyo Lengai volcano. After determining a new ML scaling relation, we determine a b-value of 0.87 ± 0.03. Focal mechanisms for 65 earthquakes, and 13 from a catalogue prior to our array reveal an along-axis stress rotation of ∼60° in the magmatically active zone. New and prior mechanisms show predominantly normal slip along steep nodal planes, with extension directions ∼ N90°E north and south of an active volcanic chain consistent with geodetic data, and ∼ N150°E in the volcanic chain. The stress rotation facilitates strain transfer from border fault systems, the locus of early stage deformation, to the zone of magma intrusion in the central rift. Our seismic, structural, and geochemistry results indicate that frequent lower crustal earthquakes are promoted by elevated pore pressures from volatile degassing along border faults, and hydraulic fracture around the margins of magma bodies. Results indicate that earthquakes are largely driven by stress state around inflating magma bodies.

The Campi Flegrei Blind Test: Evaluating the Imaging Capability of Local Earthquake Tomography in a Volcanic Area

During the 1982–1984 bradyseismic crises in the Campi Flegrei area (Italy), the University of Wisconsin deployed a network of seismological stations to record local earthquakes. In order to analyse the potential of the recorded data in terms of tomographic imaging, a blind test was recently set up and carried out in the framework of a research project. A model representing a hypothetical 3D structure of the area containing the Campi Flegrei caldera was also set up, and a synthetic dataset of time arrivals was in turn computed. The synthetic dataset consists of several thousand P- and S-time arrivals, computed at about fourteen stations. The tomographic inversion was performed by four independent teams using different methods. The teams had no knowledge of either the input velocity model or the earthquake hypocenters used to create the synthetic dataset. The results obtained by the different groups were compared and analysed in light of the true model. This work provides a thorough analysis of the earthquake tomography potential of the dataset recording the seismic activity at Campi Flegrei in the 1982–1984 period. It shows that all the tested earthquake tomography methods provide reliable low-resolution images of the background velocity field of the Campi Flegrei area, but with some differences. However, none of them succeeds in detecting the hypothetical structure details (i.e. with a size smaller than about 1.5–2 km), such as a magmatic chamber 4 km deep and especially the smaller, isolated bodies, which represent possible magmatic chimneys and intrusions.