Constraints on the Geometry and Frictional Properties of the Main Himalayan Thrust Using Coseismic, Postseismic, and Interseismic Deformation in Nepal

The geometry and frictional properties of a fault system are key parameters required to understand its seismic behavior. The Main Himalayan Thrust in Nepal is the type example of a continental megathrust and forms part of a fault system which accommodates a significant fraction of India‐Eurasia convergence. Despite extensive study of this zone of shortening, the geometry of the fault system remains controversial. Here, we use interseismic, coseismic, and postseismic geodetic data in Nepal to investigate the proposed downdip geometries. We use interseismic and coseismic data from previous studies, acquired before and during the 2015 urn:x-wiley:jgrb:media:jgrb53987:jgrb53987-math-0001 7.8 Gorkha earthquake. We then supplement these by processing our own postseismic deformation data, acquired following the Gorkha earthquake. We find that kinematic modeling of geodetic data alone cannot easily distinguish between the previously proposed geometries. We therefore develop a mechanical joint coseismic‐postseismic slip inversion which simultaneously solves for the distribution of coseismic slip and rate‐strengthening friction parameters. We run this inversion using the proposed geometries and find that they are all capable of explaining the majority of geodetic data. We find values for the rate parameter, urn:x-wiley:jgrb:media:jgrb53987:jgrb53987-math-0002, from the rate‐and‐state friction law that are between 0.8 and urn:x-wiley:jgrb:media:jgrb53987:jgrb53987-math-0003, depending on the geometry used. These values are in agreement with results from laboratory studies and those inferred from other earthquakes. We suggest that the limitations of earthquake cycle geodesy partly explain the continued controversy over the geometry and role of various faults in the Nepal Himalaya

Earthquake distribution patterns in Africa: their relationship to variations in lithospheric and geological structure, and their rheological implications

We use teleseismic waveform inversion, along with depth phase analysis, to constrain the centroid depths and source parameters of large African earthquakes. The majority of seismic activity is concentrated along the East African Rift System, with additional active regions along stretches of the continental margins in north and east Africa, and in the Congo Basin. We examine variations in the seismogenic thickness across Africa, based on a total of 227 well-determined earthquake depths, 112 of which are new to this study. Seismogenic thickness varies in correspondence with lithospheric thickness, as determined from surface wave tomography, with regions of thick lithosphere being associated with seismogenic thicknesses of up to 40 km. In regions of thin lithosphere, the seismogenic thickness is typically limited to ≤20 km. Larger seismogenic thicknesses also correlate with regions that have dominant tectonothermal ages of ≥1500 Ma, where the East African Rift passes around the Archean cratons of Africa, through the older Proterozoic mobile belts. These correlations are likely to be related to the production, affected by method and age of basement formation, and preservation, affected by lithospheric thickness, of a strong, anhydrous lower crust. The Congo Basin contains the only compressional earthquakes in the continental interior. Simple modelling of the forces induced by convective support of the African plate, based on long-wavelength free-air gravity anomalies, indicates that epeirogenic effects are sufficient to account for the localization and occurrence of both extensional and compressional deformation in Africa. Seismicity along the margins of Africa reflects a mixture between oceanic and continental seismogenic characteristics, with earthquakes in places extending to 40 km depth

Reconciling Geophysical and Petrological Estimates of the Thermal Structure of Southern Tibet

The thermal structure of the Tibetan plateau—the largest orogenic system on Earth—remains largely unknown. Numerous avenues provide fragmentary pressure/temperature information, both at the present (predominantly informed though geophysical observation) and on the evolution of the thermal structure over the recent past (combining petrological, geochemical, and geophysical observables). However, these individual constraints have proven hard to reconcile with each other. Here, we show that models for the simple underthrusting of India beneath southern Tibet are capable of matching all available constraints on its thermal structure, both at the present day and since the Miocene. Many parameters in such models remain poorly constrained, and we explore the various trade-offs among the competing influences these parameters may have. However, three consistent features to such models emerge: (i) that present-day geophysical observations require the presence of relatively cold underthrust Indian lithosphere beneath southern Tibet; (ii) that geochemical constraints require the removal of Indian mantle from beneath southern Tibet at some point during the early Miocene, although the mechanism of this removal, and whether it includes the removal of any crustal material, is not constrained by our models; and (iii) that the combination of the southern extent of Miocene mantle-derived magmatism and the present-day geophysical structure and earthquake distribution of southern Tibet require that the time-averaged rate of underthrusting of India relative to central Tibet since the middle Miocene has been faster than it is at present

An explanation for the age independence of oceanic elastic thickness estimates from flexural profiles at subduction zones, and implications for continental rheology

Most properties of oceanic lithosphere are widely observed to be dependent on the age of the plate, such as water depth, heat flow, and seismogenic thickness. However, estimates of the ‘effective elastic thickness' of oceanic lithosphere based on the deflection of the plate as it enters a subduction zone show little correlation with the age of the incoming lithosphere. This paradox requires reconciliation if we are to gain a full understanding of the structure, rheology, and behaviour of oceanic lithosphere. Here, we show that the permanent deformation of the plate due to outer-rise faulting, combined with uncertainties in the yield stress of the lithosphere, the in-plane forces transmitted through subduction zones, and the levels of noise in bathymetric and gravity data, prevents simple elastic plate modelling from accurately capturing the underlying rheological structure of the incoming plate. The age-independent estimates of effective elastic thickness obtained by purely elastic plate modelling are therefore not likely to represent the true rheology of the plate, and hence are not expected to correspond to the plate age. Similar effects may apply to estimates of elastic thickness from continental forelands, with implications for our understanding of continental rheology

Forearc collapse, plate flexure, and seismicity within the downgoing plate along the Sunda Arc west of Sumatra

Deformation within the downgoing oceanic lithosphere seawards 3 of subduction zones is typically characterised by regimes of shallow 4 extension and deeper compression, due to the bending of the oceanic 5 plate as it dips into the subduction zone. However, o shore Suma- 6 tra there are shallow compressional earthquakes within the down- 7 going oceanic plate outboard of the region of high slip in the 2004 8 Aceh-Andaman earthquake, occurring at the same depth as exten- 9 sional faulting further seaward from the trench. A clear separation is 10 seen in the location of intraplate earthquakes, with extensional earth- 11 quakes occurring further seawards than compressional earthquakes at 12 the same depth within the plate. The adjacent section of the fore- 13 arc prism west of Aceh is also anomalous in its morphology, charac- 14 terised by a wide prism with a steep bathymetric front and broad, 15 gradually-sloping top. This shape is in contrast to the narrower and 16 more smoothly-sloping prism to the south, and along other subduction 1 17 zones. The anomalous near-trench intraplate earthquakes and prism 18 morphology are likely to be the result of the geologically-rapid gravi- 19 tational collapse of the forearc, which leads to induced bending within 20 the subducting plate, and the distinctive plateau-like morphology of 21 the forearc. Such collapse of the forearc could be caused by changes 22 through time of the material properties of the forearc rocks, or of the 23 thickness of the sediments entering the subduction zone

Resolving the location of small intracontinental earthquakes using Open Access seismic and geodetic data: lessons from the 2017 January 18 mb 4.3, Ténéré, Niger, earthquake

A low-magnitude earthquake was recorded on 2017 January 18, in the Ténéré desert in northern Niger. This intraplate region is exceptionally sparsely covered with seismic stations and the closest open seismic station, G.TAM in Algeria at a distance of approximately 600 km, was unusually and unfortunately not operational at the time of the event. Body-wave magnitude estimates range from mb 4.2 to mb 4.7 and both seismic location and magnitude constraints are dominated by stations at teleseismic distances. The seismic constraints are strengthened considerably by array stations of the International Monitoring System for verifying compliance with the Comprehensive Nuclear Test-Ban-Treaty. This event, with magnitude relevant to low-yield nuclear tests, provides a valuable validation of the detection and location procedure for small land-based seismic disturbances at significant distances. For seismologists not in the CTBT system, the event is problematic as data from many of the key stations are not openly available. We examine the uncertainty in published routinely determined epicentres by performing multiple Bayesloc location estimates with published arrival times considering both all published arrival times and those from open stations only. This location exercise confirms lateral uncertainties in seismologically derived location no smaller than 10 km. Coherence for interferometric synthetic aperture radar in this region is exceptionally high, and allows us to confidently detect a displacement of the order 6 mm in the time frame containing the earthquake, consistent with the seismic location estimates, and with a lateral length scale consistent with an earthquake of this size, allowing location constraint to within one rupture length (≤5 km)—significantly reducing the lateral uncertainty compared with relying on seismological data only. Combining Open Access-only seismological and geodetic data, we precisely constrain the source location, and conclude that this earthquake likely had a shallow source. We then discuss potential ways to continue the integration of geodetic data in the calibration of seismological earthquake location

An enigmatic earthquake in the continental mantle lithosphere of stable North America

The existence of earthquakes within continental lithospheric mantle remains a highly controversial topic. Here, we present a detailed set of seismological analyses confirming the occurrence of a mantle earthquake beneath the Wind River Range of central Wyoming. Combining regional waveform inversion with the analysis of the delay and relative amplitudes of teleseismically-observed depth phases, we demonstrate that the 2013 Wind River earthquake – a MW 4.7 highly-oblique thrust-faulting event – occurred at 75±8 km, well beneath the base of the crust. The magnitude, mechanism, and location of this earthquake suggest that it represents simple brittle failure at relatively high temperatures within the mantle lithosphere, as a result of tectonic, rather than magmatic, processes

Improving the Resolving Power of InSAR for Earthquakes Using Time Series: A Case Study in Iran

Interferometric Synthetic Aperture Radar (InSAR) is an established method to measure earthquake surface displacements. However, due to decorrelation and atmospheric noise, only a certain fraction of earthquakes is readily observable with single interferograms. To enhance the potential of retrieving InSAR earthquake observations, we apply InSAR time series analysis and use several recent earthquakes (Mw 5.6–6.3, 2018–2019) in Iran as case studies. We find that the coseismic displacement signals of these earthquakes, which might not be discernible within single interferograms, are better resolved using our approach. We reconstruct the coseismic deformation fields by fitting surface displacements using a time series approach. We find that the reconstructed coseismic deformation fields yield more robust and seismologically consistent earthquake modeling results when compared to single coseismic interferograms. Our work suggests that a time series approach is an effective way to improve the resolving power of InSAR for earthquake studies

Imbalanced Moment Release Within Subducting Plates During Initial Bending and Unbending

Abstract: Internal deformation within the downgoing plate in subduction zones to accommodate the bending of the plate as it starts to subduct is reflected in widespread intraplate seismicity. This seismicity, extending from the outer rise and outer trench slope down to intermediate depths within the slab, exhibits a wide variety of focal mechanism types, often indicative of the accumulation and recovery of down‐dip curvature through the combination of normal‐faulting and thrust‐faulting earthquakes. In the idealized case, where the bulk internal deformation is recovered and slabs descend as a straight plate into the deeper mantle, we might expect the seismic moment released in both extension and compression to balance. However, a number of factors may complicate this: the thermal, compositional, and rheological evolution of the slab as it subducts, changes in the proportion of deformation accommodated seismically, and whether the slab undergoes any permanent deformation. These factors may result in intraslab deformation and seismicity unrelated to the slab curvature. Here, we assess earthquake moment release in intraslab settings around the world, focusing on those subduction systems with relatively simple slab geometries. Consideration of several regions around the western Pacific and eastern Indian Ocean indicates that substantially more deformation is accommodated seismically during bending than during unbending, and that in both settings, significantly more moment release reflects down‐dip extension than down‐dip compression. This suggests that, although the location of seismicity is clearly related to changes in slab curvature, there is a component of permanent, unrecovered down‐dip extension in many subducting slabs