Paramesh Banerjee

The 2010 Mw 7.8 Mentawai earthquake: Very shallow source of a rare tsunami earthquake determined from tsunami field survey and near-field GPS data

[1] The Mw 7.8 October 2010Mentawai, Indonesia, earthquake was a“tsunami earthquake,” a rare type of earthquake that generates a tsunami much larger than expected based on the seismicmagnitude.Itproducedalocallydevastatingtsunami,withrunupcommonlyinexcess of 6 m. We examine this event using a combination of high-rate GPS data, from instruments located on the nearby islands, and a tsunami field survey. The GPS displacement time series are deficient in high-frequency energy, and show small coseismic displacements ( 16 m. Our modeling results show that the combination of the small GPS displacements and large tsunami can only be explained by high fault slip at very shallow depths, far from the islands and close to the oceanic trench. Inelastic uplift of trench sediments likely contributed to the size of the tsunami. Recent results for the 2011 Mw 9.0 Tohoko-Oki earthquake have also shown shallow fault slip, but the results from our study, which involves a smaller earthquake, provide much stronger constraints on how shallow the rupture can be, with the majority of slip for the Mentawai earthquake occurring at depths of <6 km. This result challenges the conventional wisdom that the shallow tips of subduction megathrusts are aseismic, and therefore raises important questions both about the mechanical properties of the shallow fault zone and the potential seismic and tsunami hazard of this shallow region.

The 2012 Mw 8.6 Wharton Basin sequence: A cascade of great earthquakes generated by near-orthogonal, young, oceanic mantle faults

We improve constraints on the slip distribution and geometry of faults involved in the complex, multisegment, Mw 8.6 April 2012 Wharton Basin earthquake sequence by joint inversion of high-rate GPS data from the Sumatran GPS Array (SuGAr), teleseismic observations, source time functions from broadband surface waves, and far-field static GPS displacements. This sequence occurred under the Indian Ocean, ∼400 km offshore Sumatra. The events are extraordinary for their unprecedented rupture of multiple cross faults, deep slip, large strike-slip magnitude, and potential role in the formation of a discrete plate boundary between the Indian and Australian plates. The SuGAr recorded static displacements of up to ∼22 cm, along with time-varying arrivals from the complex faulting, which indicate that the majority of moment release was on young, WNW trending, right-lateral faults, counter to initial expectations that an old, lithospheric, NNE trending fracture zone played the primary role. The new faults are optimally oriented to accommodate the present-day stress field. Not only was the greatest moment released on the younger faults, but it was these that sustained very deep slip and high stress drop (>20 MPa). The rupture may have extended to depths of up to 60 km, suggesting that the oceanic lithosphere in the northern Wharton Basin may be cold and strong enough to sustain brittle failure at such depths. Alternatively, the rupture may have occurred with an alternative weakening mechanism, such as thermal runaway.

Stress-driven relaxation of heterogeneous upper mantle and time-dependent afterslip following the 2011 Tohoku earthquake

Stress-driven Relaxation of Heterogeneous Upper Mantle and Time-dependent Afterslip Following the 2011 Tohoku Earthquake Yan Hu 1 , Roland Burgmann 1 , Naoki Uchida 2 , Paramesh Banerjee 3 , and Jeffrey T. Freymueller 4 University of California Berkeley, Berkeley, CA, United States. Graduate School of Science, Tohoku University, Japan Earth Observatory of Singapore, Nanyang Technological University, Singapore, Singapore. University of Alaska Fairbanks, Fairbanks, AK, United States. Corresponding author: Yan Hu, Earth and Planetary Science, University of California Berkeley, 307 McCone Hall, Berkeley, CA, 94720, USA. (yhu@seismo.berkeley.edu) This article has been accepted for publication and undergone full peer review but has not been through the copyediting, typesetting, pagination and proofreading process which may lead to differences between this version and the Version of Record. Please cite this article as doi: 10.1002/2015JB012508 ©2015 American Geophysical Union. All rights reserved.

Asthenosphere rheology inferred from observations of the 2012 Indian Ocean earthquake

The concept of a weak asthenospheric layer underlying Earth’s mobile tectonic plates is fundamental to our understanding of mantle convection and plate tectonics. However, little is known about the mechanical properties of the asthenosphere (the part of the upper mantle below the lithosphere) underlying the oceanic crust, which covers about 60 per cent of Earth’s surface. Great earthquakes cause large coseismic crustal deformation in areas hundreds of kilometres away from and below the rupture area. Subsequent relaxation of the earthquake-induced stresses in the viscoelastic upper mantle leads to prolonged postseismic crustal deformation that may last several decades and can be recorded with geodetic methods. The observed postseismic deformation helps us to understand the rheological properties of the upper mantle, but so far such measurements have been limited to continental-plate boundary zones. Here we consider the postseismic deformation of the very large (moment magnitude 8.6) 2012 Indian Ocean earthquake to provide by far the most direct constraint on the structure of oceanic mantle rheology. In the first three years after the Indian Ocean earthquake, 37 continuous Global Navigation Satellite Systems stations in the region underwent horizontal northeastward displacements of up to 17 centimetres in a direction similar to that of the coseismic offsets. However, a few stations close to the rupture area that had experienced subsidence of up to about 4 centimetres during the earthquake rose by nearly 7 centimetres after the earthquake. Our three-dimensional viscoelastic finite-element models of the post-earthquake deformation show that a thin (30–200 kilometres), low-viscosity (having a steady-state Maxwell viscosity of (0.5–10) × 1018 pascal seconds) asthenospheric layer beneath the elastic oceanic lithosphere is required to produce the observed postseismic uplift.

Upper-mantle water stratification inferred from observations of the 2012 Indian Ocean earthquake

Water, the most abundant volatile in Earth’s interior, preserves the young surface of our planet by catalysing mantle convection, lubricating plate tectonics and feeding arc volcanism. Since planetary accretion, water has been exchanged between the hydrosphere and the geosphere, but its depth distribution in the mantle remains elusive. Water drastically reduces the strength of olivine and this effect can be exploited to estimate the water content of olivine from the mechanical response of the asthenosphere to stress perturbations such as the ones following large earthquakes. Here, we exploit the sensitivity to water of the strength of olivine, the weakest and most abundant mineral in the upper mantle, and observations of the exceptionally large (moment magnitude 8.6) 2012 Indian Ocean earthquake to constrain the stratification of water content in the upper mantle. Taking into account a wide range of temperature conditions and the transient creep of olivine, we explain the transient deformation in the aftermath of the earthquake that was recorded by continuous geodetic stations along Sumatra as the result of water- and stress-activated creep of olivine. This implies a minimum water content of about 0.01 per cent by weight—or 1,600 H atoms per million Si atoms—in the asthenosphere (the part of the upper mantle below the lithosphere). The earthquake ruptured conjugate faults down to great depths, compatible with dry olivine in the oceanic lithosphere. We attribute the steep rheological contrast to dehydration across the lithosphere–asthenosphere boundary, presumably by buoyant melt migration to form the oceanic crust.

Hunt for slow slip events along the Sumatran subduction zone in a decade of continuous GPS data

10 pages, 4 figures, supporting information http://dx.doi.org/10.1002/2015JB012503.-- This work is Earth Observatory of Singapore paper number 104

Postseismic relaxation in Kashmir and lateral variations in crustal architecture and materials

Thirty horizontal displacement time series from GPS sites in the area around the 2005 Kashmir earthquake show lateral spatial variations in displacement magnitude and relaxation time for the postseismic interval from 2005 to 2012. The observed spatial pattern of surface displacements can only be reproduced by finite element models of postseismic deformation in elastic over viscoelastic crust that include lateral differences in both the thickness of the elastic layer and the viscosity of the viscoelastic layer. Solutions reproducing the sign of horizontal displacements everywhere in the epicentral region also require afterslip on the portion of the fault dislocation in the viscoelastic layer but not in the elastic lid. Although there are substantial tradeoffs among contributions to postseismic displacements of the surface, the observations preclude both crustal homogeneity and shallow afterslip. In the best family of solutions, the thickness of the elastic upper crust differs by a factor of 5 and the viscosity of the middle and lower crust by an order of magnitude between domains north and south of a suture zone containing the Main Boundary Thrust and Main Mantle Thrust.

Piecemeal Rupture of the Mentawai Patch, Sumatra: The 2008 Mw 7.2 North Pagai Earthquake Sequence

The 25 February 2008 Mw 7.2 North Pagai earthquake partially ruptured the middle section of the Mentawai patch of the Sunda megathrust, offshore Sumatra. The patch has been forecast to generate a great earthquake in the next few decades. However, in the current cycle the patch has so far broken in a sequence of partial ruptures, one of which was the 2008 event, illustrating the potential of the patch to generate a spectrum of earthquake sizes. We estimate the coseismic slip distribution of the 2008 event by jointly inverting coseismic offsets from GPS and InSAR. We then estimate afterslip with 5.6 years of cumulative GPS displacements. Our results suggest that the estimated afterslip partially overlaps the coseismic rupture. The overlap of coseismic rupture and afterslip can be explained conceptually by a simple rate-and-state model where the degree of overlapping is controlled by the dynamic weakening and the critical nucleation size in the velocity-weakening area. Comparing our rate-and-state model results with our geodetic inversion results, we suggest that the part of the coseismic rupture that does not overlap with the afterslip may represent a velocity-weakening region, while the overlapping part may represent a velocity-strengthening region.