An official website of the United States government
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.2Mwmegathrust 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.
more »
« less
Strong Upper‐Plate Heterogeneity at the Hikurangi Subduction Margin (North Island, New Zealand) Imaged by Adjoint Tomography
Abstract We use earthquake‐based adjoint tomography to invert for three‐dimensional structure of the North Island, New Zealand, and the adjacent Hikurangi subduction zone. The study area, having a shallow depth to the plate interface below the North Island, offers a rare opportunity for imaging material properties at an active subduction zone using land‐based measurements. Starting from an initial model derived using ray tomography, we perform iterative model updates using spectral element and adjoint simulations to fit waveforms with periods ranging from 4–30 s. We perform 28 model updates using an L‐BFGS optimization algorithm, improving data fit and introducingP‐ andS‐wave velocity changes of up to ±30%. Resolution analysis using point spread functions show that our measurements are most sensitive to heterogeneities in the upper 30 km. The most striking velocity changes coincide with areas related to the active Hikurangi subduction zone. Lateral velocity structures in the upper 5 km correlate well with New Zealand geology. The inversion reveals increased along‐strike heterogeneity on the margin. In Cook Strait we observe a low‐velocity zone interpreted as deep sedimentary basins. In the central North Island, low‐velocity anomalies are linked to surface geology, and we relate velocity structures at depth to crustal magmatic activity below the Taupō Volcanic Zone. Our velocity model provides more accurate synthetic seismograms with respect to the initial model, better constrains small (50 km), shallow (15 km) and near‐offshore velocity structures, and improves our understanding of volcanic and tectonic structures related to the active Hikurangi subduction zone.
more »
« less
- Award ID(s):
- 2052839
- PAR ID:
- 10483042
- Publisher / Repository:
- Journal of Geophysical Research
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Solid Earth
- Volume:
- 127
- Issue:
- 1
- ISSN:
- 2169-9313
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Using data from 186 stations belonging to the USArray Transportable Array, a three‐dimensional shear wave velocity model for the southeastern United States is constructed for the top 180 km by a joint inversion of receiver functions and Rayleigh wave phase velocity dispersion computed from ambient noise and teleseismic earthquake data. The resulting shear wave velocity model and the crustal thickness and Vp/Vs () measurements show a clear spatial correspondence with major surficial geological features. The distinct low velocities observed in the depth range of 0–25 km beneath the eastern Gulf Coastal Plain reflect the thick layer of unconsolidated or poorly consolidated sediments atop the crystalline crust. The low(1.70–1.74) and slow lowermost crustal velocities observed beneath the eastern Southern Appalachian Mountains (including the Carolina Terrane and Inner Piedmont) relative to the adjacent Blue Ridge Mountains and Valley and Ridge can be interpreted by lower crustal delamination followed by relamination. The Osceola intrusive complex in the central Suwannee Terrane has similar crustal characteristics as the eastern Southern Appalachian Mountains and thus can similarly be attributed to crustal delamination/relamination processes. The Grenville Province and adjacent areas possess relatively highvalues which can be attributed to mafic intrusion associated with crustal extension in a recently recognized segments of the eastern arm of the Proterozoic Midcontinent Rift.more » « less
-
Abstract Shear attenuation provides insights into the physical and chemical state of the upper mantle. Yet, observations of attenuation are infrequent in the oceans, despite recent proliferation of arrays of ocean‐bottom seismometers (OBSs). Studies of attenuation in marine environments must overcome unique challenges associated with strong oceanographic noise at the seafloor and data loss during OBS recovery in addition to untangling the competing influences of elastic focusing, local site amplification, and anelastic attenuation on surface‐wave amplitudes. We apply Helmholtz tomography to OBS data to simultaneously resolve array‐averaged Rayleigh wave attenuation and maps of site amplification at periods of 20–150 s. The approach explicitly accounts for elastic focusing and defocusing due to lateral velocity heterogeneity using wavefield curvature. We validate the approach using realistic wavefield simulations at the NoMelt Experiment and Juan de Fuca (JdF) plate, which represent endmember open‐ocean and coastline‐adjacent environments, respectively. Focusing corrections are successfully recovered at both OBS arrays, including at periods <35 s at JdF where coastline effects result in strong multipathing. When applied to real data, our observations of Rayleigh wave attenuation at NoMelt and JdF revise previous estimates. At NoMelt, we observe a low attenuation lithospheric layer (> 1,500) overlying a highly attenuating asthenospheric layer (∼ 50 to 70). At JdF, we find a broad peak in attenuation (∼ 50 to 60) centered at a depth of 100–130 km. We also report strong local site amplification at the JdF Ridge (>10% at 31 s period), which can be used to refine models of crust and shallow mantle structure.more » « less
-
Abstract The 2021MW6.0 Yangbi, Yunnan strike‐slip earthquake occurred on an unmapped crustal fault near the Weixi‐Qiaoho‐Weishan Fault along the southeast margin of the Tibetan Plateau. Using near‐source broadband seismic data from ChinArray, we investigate the spatial and temporal rupture evolution of the mainshock using apparent moment‐rate functions (AMRFs) determined by the empirical Green's function (EGF) method. Assuming a 1D line source on the fault plane, the rupture propagated unilaterally southeastward (∼144°) over a rupture length of ∼8.0 km with an estimated rupture speed of 2.1 km/s to 2.4 km/s. A 2D coseismic slip distribution for an assumed maximum rupture propagation speed of 2.2 km/s indicates that the rupture propagated to the southeast ∼8.0 km along strike and ∼5.0 km downdip with a peak slip of ∼2.1 m before stopping near the largest foreshock, where three bifurcating subfaults intersect. Using the AMRFs, the radiated energy of the mainshock is estimated as ∼. The relatively low moment scaled radiated energyof 1.5 × 10−5and intense foreshock and aftershock activity might indicate reactivation of an immature fault. The earthquake sequence is mainly distributed along a northwest‐southeast trend, and aftershocks and foreshocks are distributed near the periphery of the mainshock large‐slip area, suggesting that the stress in the mainshock slip zone is significantly reduced to below the level for more than a few overlapping aftershock to occur.more » « less
-
Abstract Silicon stable isotope ratios (30Si) of over 150 stream water samples were measured during seven storm events in six small critical zone observatory (CZO) catchments spanning a wide range in climate (sub‐humid to wet, tropical) and lithology (granite, volcanic, and mixed sedimentary). Here we report a cross‐site analysis of this dataset to gain insight into stream30Si variability across low‐order catchments and to identify potential climate (i.e., runoff), hydrologic, lithologic, and biogeochemical controls on observed stream Si chemical and isotopic signatures. Event‐based30Si exhibit variability both within and across sites (−0.22‰ to +2.27‰) on the scale of what is observed globally in both small catchments and large rivers. Notably, each site shows distinct30Si signatures that are preserved even after normalization for bedrock composition. Successful characterization of observed cross‐site behavior requires the merging of two distinct frameworks in a novel combined model describing both non‐uniform fluid transit time distributions and multiple fractionating pathways in application to low‐order catchments. The combined model reveals that site‐specific architecture (i.e., biogeochemical reaction pathways and hydrologic routing) regulates stream silicon export signatures even when subject to extreme precipitation events.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1029/2021JB022865
