Dietrich Lange

Structure of the seismogenic zone of the southcentral Chilean margin revealed by local earthquake traveltime tomography

We use traveltime data of local earthquakes and controlled sources observed by a large, temporary, amphibious seismic network to reveal the anatomy of the southcentral Chilean subduction zone (37–39°S) between the trench and the magmatic arc. At this location the giant 1960 earthquake (M = 9.5) nucleated and ruptured almost 1000 km of the subduction megathrust. For the three-dimensional tomographic inversion we used 17,148 P wave and 10,049 S wave arrival time readings from 439 local earthquakes and 94 shots. The resolution of the tomographic images was explored by analyzing the model resolution matrix and conducting extensive numerical tests. The downgoing lithosphere is delineated by high seismic P wave velocities. High vp/vs ratio in the subducting slab reflects hydrated oceanic crust and serpentinized uppermost oceanic mantle. The subducting oceanic crust can be traced down to a depth of 80 km, as indicated by a low velocity channel. The continental crust extends to approximately a 50-km depth near the intersection with the subducting plate. This suggests a wide contact zone between continental and oceanic crust of about 150 km, potentially supporting the development of large asperities. Eastward the crustal thickness decreases again to a minimum of about a 30-km depth. Relatively low vp/vs at the base of the forearc does not support a large-scale serpentinization of the mantle wedge. Offshore, low vp and high vp/vs reflect young, fluid-saturated sediments of forearc basins and the accretionary prism.

Interaction between forearc and oceanic plate at the south‐central Chilean margin as seen in local seismic data

We installed a dense, amphibious, temporary seismological network to study the seismicity and structure of the seismogenic zone in southern Chile between 37° and 39°S, the nucleation area of the great 1960 Chile earthquake. 213 local earthquakes with 14.754 onset times were used for a simultaneous inversion for the 1-D velocity model and precise earthquake locations. Relocated artificial shots suggest an accuracy of the earthquake hypocenter of about 1 km (horizontally) and 500 m (vertically). Crustal events along trench-parallel and transverse, deep-reaching faults reflect the interseismic transpressional deformation of the forearc crust due to the subduction of the Nazca plate. The transverse faults seems to accomplish differential lateral stresses between subduction zone segments. Many events situated in an internally structured, planar seismicity patch at 20 to 40 km depth near the coast indicate a stress concentration at the plates interface at 38°S which might in part be induced by the fragmented forearc structure.

Seismicity and geometry of the south Chilean subduction zone (41.5°S-43.5°S) : Implications for controlling parameters

In 2005 an amphibious seismic network was deployed on the Chilean forearc between 41.75°S and 43.25°S. 364 local events were observed in a 11-month period. A subset of the P and S arrival times were inverted for hypocentral coordinates, 1-D velocity structure and station delays. Main seismic activity occurred predominantly in a belt parallel to the coast of Chiloe Island in a depth range of 12–30 km presumably related to the plate interface. The 30° inclination of the shallow part of the Wadati-Benioff zone is similar to observations further north indicating that oceanic plate age is not controlling the subduction angle of the shallower part for the Chilean subduction zone. The down-dip termination of abundant intermediate depth seismicity at approximately 70 km depth seems to be related to the young age (and high temperature) of the oceanic plate. Crustal seismicity is associated with the Liquine-Ofqui fault zone and active volcanoes.

Splay fault activity revealed by aftershocks of the 2010 Mw 8.8 Maule earthquake, central Chile

Splay faults, large thrust faults emerging from the plate boundary to the seafloor in subduction zones, are considered to enhance tsunami generation by transferring slip from the very shallow dip of the megathrust onto steeper faults, thus increasing vertical displacement of the seafloor. These structures are predominantly found offshore, and are therefore difficult to detect in seismicity studies, as most seismometer stations are located onshore. The Mw (moment magnitude) 8.8 Maule earthquake on 27 February 2010 affected ∼500 km of the central Chilean margin. In response to this event, a network of 30 ocean-bottom seismometers was deployed for a 3 month period north of the main shock where the highest coseismic slip rates were detected, and combined with land station data providing onshore as well as offshore coverage of the northern part of the rupture area. The aftershock seismicity in the northern part of the survey area reveals, for the first time, a well-resolved seismically active splay fault in the submarine forearc. Application of critical taper theory analysis suggests that in the northernmost part of the rupture zone, coseismic slip likely propagated along the splay fault and not the subduction thrust fault, while in the southern part it propagated along the subduction thrust fault and not the splay fault. The possibility of splay faults being activated in some segments of the rupture zone but not others should be considered when modeling slip distributions.

Separating rapid relocking, afterslip, and viscoelastic relaxation: An application of the postseismic straightening method to the Maule 2010 cGPS

The postseismic deformation captured with continuous Global Positioning System (cGPS) monitoring following many recent mega-thrust events has been shown to be a signal composed of two dominant processes: afterslip on the plate interface and viscoelastic relaxation of the continental and oceanic mantles in response to the coseismic stress perturbation. Following the south-central Chile 2010 Maule Mw 8.8 earthquake, the time series from the regional cGPS network show a distinct curvature in the pathway of the horizontal motion that is not easily fit by a stationary decaying pattern of afterslip in combination with viscoelastic relaxation. Here we show that with realistic assumptions about the long-term decay of the afterslip signal, the postseismic signal can be decomposed into three first-order contributing processes: plate interface re-locking, plate interface afterslip, and mantle viscoelastic relaxation. From our analyses we conclude that the plate interface recovers its interseismic locking state rapidly (model space ranges between an instant recovery and a period of 1 year); a finding that supports laboratory experimental evidence as well as some recent studies of aftershocks and postseismic surface deformation. Furthermore, re-locking is the main cause of the curvature in the cGPS signal, and this study presents a plausible range of geodetic re-locking rates following a megathrust earthquake.

No significant steady state surface creep along the North Anatolian Fault offshore Istanbul: Results of 6 months of seafloor acoustic ranging

The submarine Istanbul-Silivri fault segment, within 15km of Istanbul, is the only portion of the North Anatolian Fault that has not ruptured in the last 250years. We report first results of a seafloor acoustic ranging experiment to quantify current horizontal deformation along this segment and assess whether the segment is creeping aseismically or accumulating stress to be released in a future event. Ten transponders were installed to monitor length variations along 15 baselines. A joint least squares inversion for across-fault baseline changes, accounting for sound speed drift at each transponder, precludes fault displacement rates larger than a few millimeters per year during the 6month observation period. Forward modeling shows that the data better fit a locked state or a very moderate surface creep-less than 6mm/yr compared to a far-field slip rate of over 20mm/yr-suggesting that the fault segment is currently accumulating stress.

The structure of the Sumatran Fault revealed by local seismicity

The combination of the Sunda megathrust and the (strike-slip) Sumatran Fault (SF) represents a type example of slip-partitioning. However, superimposed on the SF are geometrical irregularities that disrupt the local strain field. The largest such feature is in central Sumatra where the SF splits into two fault strands up to 35 km apart. A dense local network was installed along a 350 km section around this bifurcation, registering 1016 crustal events between April 2008 and February 2009. 528 of these events, with magnitudes between 1.1 and 6.0, were located using the double-difference relative location method. These relative hypocentre locations reveal several new features about the crustal structure of the SF. Northwest and southeast of the bifurcation, where the SF has only one fault strand, seismicity is strongly focused below the surface trace, indicating a vertical fault that is seismogenic to ∼15 km depth. By contrast intense seismicity is observed within the bifurcation, displaying streaks in plan and cross-section that indicate a complex system of faults bisecting the bifurcation. In combination with analysis of topography and focal mechanisms, we propose that the bifurcation is a strike-slip duplex system with complex faulting between the two main fault branches.

Structure of the oceanic lithosphere and upper mantle north of the Gloria fault in the eastern mid-Atlantic by receiver function analysis

Receiver functions (RF) have been used for several decades to study structures beneath seismic stations. Although most available stations are deployed on-shore, the number of ocean bottom station (OBS) experiments has increased in recent years. Almost all OBSs have to deal with higher noise levels and a limited deployment time (∼1 year), resulting in a small number of usable records of teleseismic earthquakes. Here, we use OBSs deployed as mid-aperture array in the deep ocean (4.5-5.5 km water depth) of the eastern mid-Atlantic. We use evaluation criteria for OBS data and beam forming to enhance the quality of the RFs. Although some stations show reverberations caused by sedimentary cover, we are able to identify the Moho signal, indicating a normal thickness (5-8 km) of oceanic crust. Observations at single stations with thin sediments (300-400 m) indicate that a probable sharp lithosphere-asthenosphere boundary (LAB) might exist at a depth of ∼70-80 km which is in line with LAB depth estimates for similar lithospheric ages in the Pacific. The mantle discontinuities at ∼410 km and ∼660 km are clearly identifiable. Their delay times are in agreement with PREM. Overall the usage of beam formed earthquake recordings for OBS RF analysis is an excellent way to increase the signal quality and the number of usable events.

Fault zone controlled seafloor methane seepage in the rupture area of the 2010 Maule Earthquake, Central Chile

Seafloor seepage of hydrocarbon-bearing fluids has been identified in a number of marine forearcs. However, temporal variations in seep activity and the structural and tectonic parameters that control the seepage often remain poorly constrained. Subduction-zone earthquakes for example, are often discussed to trigger seafloor seepage but causal links that go beyond theoretical considerations have not yet been fully established. This is mainly due to the inaccessibility of offshore epicentral areas, the infrequent occurrence of large earthquakes, and challenges associated with offshore monitoring of seepage over large areas and sufficient time periods. Here, we report visual, geochemical, geophysical, and modelling results and observations from the Concepcion Methane Seep Area (offshore Central Chile) located in the rupture area of the 2010 Mw. 8.8 Maule earthquake. High methane concentrations in the oceanic water column and a shallow sub-bottom depth of sulfate penetration indicate active methane seepage. The stable carbon isotope signature of the methane and hydrocarbon composition of the released gas indicate a mixture of shallow-sourced biogenic gas and a deeper sourced thermogenic component. Pristine fissures and fractures observed at the seafloor together with seismically imaged large faults in the marine forearc may represent effective pathways for methane migration. Upper-plate fault activity with hydraulic fracturing and dilation is in line with increased normal Coulomb stress during large plate-boundary earthquakes, as exemplarily modelled for the 2010 earthquake. On a global perspective our results point out the possible role of recurring large subduction-zone earthquakes in driving hydrocarbon seepage from marine forearcs over long timescales.

Systematic Changes of Earthquake Rupture with Depth: A Case Study from the 2010Mw 8.8 Maule, Chile, Earthquake Aftershock Sequence

The very shallow part of subduction megathrusts occasionally hosts tsunami earthquakes, with unusually slow rupture propagation. The aftershock sequence of the 2010 Mw 8.8 Maule earthquake, offshore Chile, provides us with the opportunity to study systematic changes in source properties for smaller earthquakes within a single segment of a subduction zone. We invert amplitude spectra for double-couple moment tensors and centroid depths of 71 aftershocks of the Maule earthquake down to magnitudes Mw 4.0. In addition, we also derive average source durations.We find that shallower earthquakes tend to have longer normalized source durations on average, similar to the pattern observed previously for larger magnitude events. This depth dependence is observable for thrust and normal earthquakes. The normalized source durations of normal- faulting earthquakes are at the lower end of those for thrust earthquakes, probably because of the higher stress drops of intraplate earthquakes compared to interplate earthquakes. We suggest from the similarity of the depth dependence of normal and thrust events and between smaller and larger magnitude earthquakes that the depth-dependent variation of rigidity, rather than frictional conditional stability at the plate interface, is primarily responsible for the observed pattern. Tsunami earthquakes probably require both low rigidity and conditionally stable frictional conditions; the presence of longduration moderate-magnitude events is therefore a helpful but not sufficient indicator for areas at risk of tsunami earthquakes.