Search for: All records

Abstract The Indo‐Burman subduction zone represents a global endmember for extreme sediment accretion and is a region characterized by ambiguous tectonic structure. The recent collection of broadband seismic data across the Indo‐Burman accretionary margin as part of the Bangladesh‐India‐Myanmar Array (BIMA) experiment provides an opportunity to investigate the subsurface velocity structure across the incoming plate of an endmember subduction system. We construct a three‐dimensional model for seismic shear velocity using a joint inversion of surface‐ and scattered‐wave constraints. Rayleigh‐wave phase velocities measured from ambient‐noise (12–25 s) and teleseismic earthquakes (20–80 s) constrain absolute shear velocities, while we constrain the locations of and relative contrasts across significant discontinuities in the subsurface using observations from scattered‐wave imaging. From the resulting inversion, we observe two model classes that characterize the evolution of consolidation within the markedly slow uppermost sediments and metasediments along a predominantly southwest‐to‐northeast trend. We interpret variations in deeper seismic structure under two proposed scenarios: (a) a Moho of ∼21–26 km depth underlying a package of metasediments and a thinned basement component, with a slow mantle lithosphere (4.2 km/s) that may contain retained melt from the onset of India‐Antarctica seafloor spreading; or (b) a Moho of ∼51–59 km depth underlying a package of metasediments, basement, and a thick slug of mafic material, which may correspond to significant Kerguelen‐plume‐related underplating. By combining constraints from highly resolved phase‐velocity estimates and scattered‐wave images, we successfully characterize the lateral transitions across the Indo‐Burman forearc margin. 

Abstract Earthquakes present severe hazards for people and economies and can be primary drivers of landscape change yet their impact to river-channel networks remains poorly known. Here we show evidence for an abrupt earthquake-triggered avulsion of the Ganges River at ~2.5 ka leading to relocation of the mainstem channel belt in the Bengal delta. This is recorded in freshly discovered sedimentary archives of an immense relict channel and a paleo-earthquake of sufficient magnitude to cause major liquefaction and generate large, decimeter-scale sand dikes >180 km from the nearest seismogenic source region. Precise luminescence ages of channel sand, channel fill, and breached and partially liquefied floodplain deposits support coeval timing of the avulsion and earthquake. Evidence for reorganization of the river-channel network in the world’s largest delta broadens the risk posed by seismic events in the region and their recognition as geomorphic agents in this and other tectonically active lowlands. The recurrence of comparable earthquake-triggered ground liquefaction and a channel avulsion would be catastrophic for any of the heavily populated, large river basins and deltas along the Himalayan arc (e.g., Indus, Ganges, Brahmaputra, Ayeyarwady). The compounding effects of climate change and human impacts heighten and extend the vulnerability of many lowlands worldwide to such cascading hazards. 

Abstract The shallow portion of a megathrust represents the zone of first contact between two colliding plates, and its rheological properties control the seismic and tsunami hazards generated by the fault. The high cost of underwater geodetic data collection results in sparse observations, leading to limited constraints on the interseismic behavior of megathrusts. The Rakhine‐Bangladesh megathrust offers a unique opportunity to probe the behavior of the shallow megathrust as it is the only ocean‐continent subduction zone where the near‐trench region is fully accessible on land. Here, we use observations from ALOS‐2 wide‐swath imagery spanning 2015 to 2022 to conduct an InSAR timeseries analysis of the overriding plate within Bangladesh and the Indo‐Myanmar Ranges. We identify a narrow pattern of alternating uplift and subsidence associated with mapped anticlines but show that it cannot be explained by slip on the megathrust or other fault structures. Instead, we argue that the deformation is likely caused by active aseismic folding within the wedge above a shallow decollement. We show that estimates of the decollement depth derived from a viscous folding model and the observed anticline spacing are in agreement with previous seismic observations of the decollement depth across the fold belt. We suggest that the role of ductile deformation in the overriding plate in subduction zones may be more important than previously recognized. 

Abstract The Indo‐Burma subduction zone is a highly oblique subduction system where the Indian plate is converging with the Eurasian plate. How strain is partitioned between the Indo‐Burma interface and upper plate Kabaw Fault, and whether the megathrust is a locked and active zone of convergence that can generate great earthquakes are ongoing debates. Here, we use data from a total of 68 Global Navigation Satellite System (GNSS) stations, including newly installed stations across the Kabaw Fault and compute an updated horizontal and vertical GNSS velocity field. We correct vertical rates for fluctuating seasonal signals by accounting for the elastic response of monsoon water on the crust. We model the geodetic data by inverting for 11,000 planar and non‐planar megathrust fault geometries and two geologically viable structural interpretations of the Kabaw Fault that we construct from field geological data, considering a basin‐scale wedge‐fault and a crustal‐scale reverse fault. We demonstrate that the Indo‐Burma megathrust is locked, converging at a rate ofmm/yr, and capable of hosting >8.2M_wmegathrust events. We also show that the Kabaw Fault is locked and accommodating strike‐slip motion at a rate ofmm/yr and converging at a rate ofmm/yr. Our interpretation of the geological, geophysical, and geodetic datasets indicates the Kabaw Fault is a crustal‐scale structure that actively absorbs a portion of the convergence previously ascribed to the Indo‐Burma megathrust. This reveals a previously unrecognized seismic hazard associated with the Kabaw Fault and slightly reduces the estimated hazard posed by megathrust earthquakes in the region. 

We utilized shear wave splitting analysis of teleseismic SKS, SKKS, and PKS phases to infer upper mantle deformational fabrics across a substantial area of Southeast Asia, where splitting measurements were previously limited. We used newly available permanent and temporary broadband seismic networks deployed across the Indo-Burma subduction zone and the eastern Indochina peninsula. The resulting 492 well-constrained splitting and 654 null measurements from 185 stations reveal clear large-scale patterns in the mantle deformational fabrics in response to the highly oblique active subduction and a large transform plate boundary. We identified two distinct domains of mantle deformation fabrics in the western Burma microplate and the eastern Indochina peninsula. In the former, trench parallel N-S fast polarization directions with an average lag time (δt) of 1.9 s are observed beneath the Indo-Burman Ranges. We suggest the observed splitting is partly due to anisotropy in the sub-slab region and relates to shear induced by the north moving Indian plate. The lithospheric fabric within the Indo-Burman Ranges and underlying subducting slab fabric contribute to produce the observed average δt of 1.9 s. The δt value decreases to an average of 1.0 s towards the back-arc until we reach the dextral Sagaing fault. In the second domain, starting approximately 100 km east of the Sagaing fault, we observe a consistent E-W fast direction with an average δt of 1.10 s in the eastern Shan-Thai and Indochina blocks. We interpret the E-W fabric as due to the deformation associated with the westward spreading of the Hainan mantle plume, possibly driven by overriding plate motion. Low velocities in the shallow mantle and late Cenozoic intraplate volcanism in this region support the plume-driven asthenospheric flow model in the Indochina peninsula. The sudden transition of the fast polarization direction from N-S to E-W along the eastern edge of the Burma microplate indicates the Sagaing fault acts as a mantle flow boundary between the subduction dominated trench parallel flow to the west and plume induced asthenospheric flow to the east. We also observed no net splitting beneath the Bengal basin which is most likely due to the presence of frozen vertical fabric resulting from the Kerguelen plume activity during Early Cretaceous. 

Abstract The Thai Meteorological Department (TMD) seismic network began development in 2008. There are a total of 71 seismic stations consisting of 26 borehole stations and 45 surface stations currently installed. The three-component data from the TMD seismic network have been widely used in previous seismological studies. In a recent analysis, we have found that sensor orientation as reported in the site metadata is sometimes significantly incorrect, especially for borehole stations. In this study, we analyze P-wave polarization data from regional and teleseismic earthquakes recorded in the network to estimate the true instrument orientation relative to geographic north and compare that to station metadata. Of the 45 surface stations, we found that at present, ~ 82% are well oriented (i.e., aligned within 0–15° of true north). However, 8 sites have sensors misoriented by more than 15°, and some stations had a temporal change in sensor orientation during an upgrade to the seismic system with replacement of the sensor. We also evaluated sensor orientations for 26 TMD borehole seismic stations, from 2018 to the 2022. For many of the borehole stations, the actual sensor orientation differs significantly from the TMD metadata, especially at short-period stations. Many of those stations have sensor misorientations approaching 180°, due to errors in the ambient noise analysis calibration techniques used during installation. We have also investigated how this sensor misorientation affects previous seismic studies, such as regional moment tensor inversion of earthquakes sources and receiver function stacking. We have found that the large deviations in sensor orientation can result in erroneous results and/or large measurement errors. A cause of the orientation error for borehole sites could be a combination of strong background surface ambient seismic noise coupled with an incorrect reference instrument response.