Shallow Deformation Research Articles

Impact of multiple weak layers on deformation of the interior of Qaidam Basin, northern Tibetan Plateau

The Qaidam Basin is the largest sedimentary basin within the Tibetan Plateau, with up to ∼15-km-thick deposits accumulated in the Cenozoic. Understanding how it deformed in response to the far-field effects of India-Eurasia collision is critical to improving our knowledge of the mechanism underlying northward plateau growth. Unlike typical compressional basins, where upper-crustal deformation concentrates at their margins, the Qaidam Basin features the development of many NW- to WNW-striking folds across the entire basin. Why crustal shortening occurred in the interior of Qaidam Basin, ∼100 km away from the margins, together with the underground geometries beneath these folds, remains unknown. Herein, based on newly acquired three- and two-dimensional (3-D and 2-D) seismic reflection data, borehole logging, and scaled physical analog modeling, we investigated the geometries, kinematics, and formation mechanisms of the folds within the interior of Qaidam Basin. For the first time, we reveal three local weak layers in the Lulehe, Upper Xiaganchaigou, and Shangyoushashan Formations, respectively. They consist mainly of mudstone intercalated with evaporites and limestones, and they have different spatial distributions that are likely confined by major faults and folds. These mechanically weak layers became locally thickened or thinned in response to tectonic loading and/or facilitated detachment slip to form many décollement folds that were observed at the surface. The shallow deformation above the weak layers is largely decoupled from underlying basement-involved faulting and folding, which mostly terminate upward in these weak layers. Analog modeling results suggest that the lowermost and widely distributed décollement layer in the Lulehe Formation likely facilitated long-distance rapid propagation of deformation into the basin interior. In sum, our study highlights the significance of multiple weak layers during Cenozoic deformation in the Qaidam Basin interior.

Geodetic imaging of magma ascent through a bent and twisted dike during the Tajogaite eruption of 2021 (La Palma, Canary Islands)

On Sept. 19th, 2021, the largest historical eruption on the island of La Palma began, which had a significant scientific, social, and economic impact. The 2021 Tajogaite eruption was characterised by short precursors, lasting only 8 days. The seismicity started on Sept. 11th with a westward and upward migration of hypocenters. Permanent GNSS stations started recording deformation on Sept. 12th on the island's western side, which reached more than 15 cm just before the eruption. After the eruption onset, the ground deformation increased, reaching a maximum on Sept. 22nd and showing a nearly steady deflation trend in the following months. To better understand the dynamics of the eruption, we exploited a joint dataset of GNSS and Sentinel-1 SBAS time series along both ascending and descending orbits. To obtain the geometry of the causative source of the ground deformation, we combined the result of a preliminary non-linear inversion and the precise location of the seismicity. The resulting geometry of the source is that of a twisted dike bending eastward. We performed inverse modelling to obtain the spatiotemporal kinematics of the opening function of the dike. The forward modelling has been realised using a 3D finite-element approach considering the island's topography. Our findings reveal a close correspondence between the magmatic intrusion and pre-eruptive seismicity. The ascent of the magma occurred along two branches, and the rheology of a previously identified ductile layer strongly affected the magma propagation process. Finally, we found evidence of an early shallow deformation, which we interpret as the effect of ascending hydrothermal fluids. Our findings highlight the need for advanced modelling to understand pre-eruptive processes in basaltic volcanoes.

Major California faults are smooth across multiple scales at seismogenic depth

Surface traces of earthquake faults are complex and segmented on multiple scales. At seismogenic depth the detailed geometry of faults and earthquake rupture is mainly constrained by earthquake locations. Standard earthquake locations are usually too diffuse to constrain multi-scale fault geometry, while differential-timing relocation mainly improves finest scale precision. NLL-SSST-coherence, an enhanced, absolute-timing earthquake location procedure, iteratively generates traveltime corrections to improve multi-scale precision and uses waveform similarity to improve fine-scale precision. Here we apply NLL-SSST-coherence to large-earthquake sequences and background seismicity along strike-slip faults in California. Our relocated seismicity at seismogenic depth along major fault segments and around large-earthquake ruptures often defines smooth, planar or arcuate, near-vertical surfaces across the sub-km to 10’s of km scales. These results show that multi-scale smooth fault segments are characteristic of major, strike-slip fault zones and may be essential to large earthquake rupture. Our results suggest that smoothness and curvature of faults influences earthquake initiation, rupture, rupture direction and arrest, and can define earthquake hazard. The results corroborate that surface traces of strike-slip fault zones reflect complex, shallow deformation and not directly simpler, main slip surfaces at depth, and support use of planar or smoothly curved faults for modeling primary earthquake rupture.

Both Longitudinal and Transverse Extension Controlling Gas Migration Through Submarine Anticlinal Ridges, New Zealand's Southern Hikurangi Margin

AbstractSub‐seabed fluid flow, gas hydrate accumulation and seafloor methane seepage are tightly interwoven processes with implications for marine biodiversity, ocean chemistry and seafloor stability. We combine long‐offset seismic reflection data with high‐resolution seismic data to investigate shallow structural deformation and its relationship to focused gas migration and hydrate accumulation in the southern Hikurangi subduction wedge. Anticlines, effective traps for focusing free gas, are characterized by both normal faults and vertical zones of hydraulic fracturing within the hydrate stability zone. The normal faults form as a result of sediment layer folding and gravitational collapse of ridges during uplift. We document both longitudinal (ridge‐parallel) and transverse (ridge‐perpendicular) extensional structures (normal faults and elongated hydraulic fracture zones) in the sub‐seafloor of anticlinal ridges. Intriguingly, gas flow through ridges close to the deformation front of the wedge exploits longitudinal structures, while ridges further inboard are characterized by gas flow along transverse structures. This highlights pronounced changes in the shallow deformation of ridges in different parts of the wedge, associated with a switching of the least and intermediate principal stress directions. It is critical to understand these shallow stress fields because they control fluid flow patterns and methane seepage out of the seafloor.

Dueling dynamics of low-angle normal fault rupture with splay faulting and off-fault damage

Despite a lack of modern large earthquakes on shallowly dipping normal faults, Holocene Mw > 7 low-angle normal fault (LANF; dip<30°) ruptures are preserved paleoseismically and inferred from historical earthquake and tsunami accounts. Even in well-recorded megathrust earthquakes, the effects of non-linear off-fault plasticity and dynamically reactivated splay faults on shallow deformation and surface displacements, and thus hazard, remain elusive. We develop data-constrained 3D dynamic rupture models of the active Mai’iu LANF that highlight how multiple dynamic shallow deformation mechanisms compete during large LANF earthquakes. We show that shallowly-dipping synthetic splays host more coseismic slip and limit shallow LANF rupture more than steeper antithetic splays. Inelastic hanging-wall yielding localizes into subplanar shear bands indicative of newly initiated splay faults, most prominently above LANFs with thick sedimentary basins. Dynamic splay faulting and sediment failure limit shallow LANF rupture, modulating coseismic subsidence patterns, near-shore slip velocities, and the seismic and tsunami hazards posed by LANF earthquakes.

Shallow deformation on the Kirby Hills fault, Sacramento–San Joaquin Delta, California (USA), revealed from high-resolution seismic reflection data and coring in a fluvial system

Abstract The Sacramento–San Joaquin Delta (Delta) in California (USA) is an important part of the state’s freshwater system and is also a major source of agricultural and natural resources. However, the Delta is traversed by a series of faults that make up the easternmost part of the San Andreas fault system at this latitude and pose seismic hazard to this region. In this study, we use new high-resolution chirp subbottom data1 to map and characterize the shallow expression of the Kirby Hills fault, where it has been mapped to cross the Sacramento River at the western extent of the Delta. The fault is buried here, but we document a broad zone of deformation associated with the eastern strand of the fault that changes in character, along strike, across ~600 m of the river channel. Radiocarbon dates from sediment cores collected in the Sacramento River provide some minimum constraints on the age of deformation. We do not observe evidence of the western strand as previously mapped. We also discuss difficulties of conducting a paleoseismologic study in a fluvial environment.

Intracontinental deformation of the Tianshan Orogen in response to India-Asia collision

How the continental lithosphere deforms far away from plate boundaries has been long debated. The Tianshan is a type-example of ongoing lithospheric deformation in an intracontinental setting. It formed during the Paleozoic accretion of the Altaids and was rejuvenated in the Cenozoic, which might be a far-field response to the India-Asia collision. Here we present seismic images of the lithosphere across the central Tianshan, which were constructed from receiver functions and Rayleigh wave dispersions along a N–S-trending linear seismic array. We observe an extensively deformed lithosphere in the Tianshan with inherited, structurally controlled brittle deformation in the shallow crust and plastic deformation near the Moho. We find that earlier multiple accretionary structures were preserved in the crust, which was deformed by pure-shear shortening in the south and thick-skinned tectonics in the north but was limitedly underthrusted by surrounding blocks. A balanced cross-section of Moho discontinuities supports the concept that intracontinental deformation in the Tianshan intensified synchronously with the direct contact between the underthrusting Indian slab and the Tarim Craton in the Late Miocene (~10 Ma). These findings provide a robust and unified seismic model for the Tianshan Orogen, and confirm that effective delivery of the India-Asia collision stress induced the rejuvenation of this intracontinental orogen.

Unrecognized post-rifting deposition and uplift history of the Erlian Basin, Northern China: Evidence from seismic reflection data and apatite fission-track thermochronology

Early Cretaceous rift basins in Northeast Asia typically experienced varying degrees of post-rifting thermal subsidence and compressional deformation. However, previous workers believed that only a thin post-rifting sequence was developed in basins west of the Great Xing'an Range, making a proper understanding of the post-rift tectonic history in this part of Northeast Asia difficult. In this study, newly acquired shallow seismic reflection data were utilized to unravel the shallow deformation of the Erlian Basin, which to the west of the Great Xing'an Range contains a thin post-rifting sequence. Furthermore, we used apatite fission-track thermochronology to analyze the post-rifting deposition and uplift processes in the Erlian Basin, for which there are limited stratigraphic records. Our results show that the Erlian Basin experienced a typical post-rifting stage whose sequence from the Albian to the Campanian was not controlled by normal faults. The corresponding sedimentary thickness was up to 2 km, which contradicts the previously reported little or no post-rifting thermal subsidence. The current thin post-rifting sequence is the remaining stratum after 1.4–1.8 km of denudation at the end of the Erlian period. In addition, by further comparing the sedimentary thickness and time-lag in the onset of post-rifting to those of other basins in Northeast Asia, we found two driving mechanisms of rifting in Northeast Asia, which occurred simultaneously: the gravitational collapse and back-arc extension. Moreover, the basins in Northeast Asia experienced a gradual weakening in the compressional deformation from east to west during the Cenomanian to the Coniacian, and an intense compressional deformation in the east of Northeast Asia and large-scale vertical uplift in the west during the Late Cretaceous Campanian to the Eocene.

Modeling long-term volcanic deformation at Kusatsu-Shirane and Asama volcanoes, Japan, using the GNSS coordinate time series

Long-term deformation of Kusatsu-Shirane and Asama volcanoes in central Japan were investigated using Global Navigation Satellite System (GNSS) measurements. Large postseismic deformation caused by the 2011 Tohoku earthquake—which obscures the long-term volcanic deformation—was effectively removed by approximating the postseismic and other recent tectonic deformation in terms of quadrature of the geographical eastings/northings. Subsequently, deformation source parameters were estimated by the Markov Chain Monte Carlo (MCMC) method and linear inversion, employing an analytical model that calculates the deformation from an arbitrary oriented prolate/oblate spheroid. The deformation source of Kusatsu-Shirane volcano was found to be a sill-like oblate spheroid located a few kilometers northwest of the Yugama crater at a depth of approximately 4 text {km}, while that of Asama was also estimated to be a sill-like oblate spheroid beneath the western flank of the edifice at a depth of approximately 12 text {km}, along with the previously reported shallow east–west striking dike at a depth of approximately 1 text {km}. It was revealed that (1) volume changes of the Kusatsu-Shirane deformation source and the shallow deformation source of Asama were correlated with the volcanic activities of the corresponding volcanoes, and (2) the Asama deep source has been steadily losing volume, which may indicate that the volcano will experience fewer eruptions in the near future.

Differences in east-west thermo-rheological structure and lithospheric deformation in the Qaidam Basin, Northern Tibetan Plateau, China

ABSTRACT To investigate the differences in east-west lithospheric deformation in the Qaidam Basin, we present thermo-rheological models of two profiles across the western Qaidam Basin (WQB) and eastern Qaidam Basin (EQB). The differences in east-west geodynamic deformation styles are also described, involving GPS motions, focal mechanisms (P axes), seismic anisotropy (SKS-wave splitting), and low-velocity zones (LVZs). The WQB is characterized by a warm destabilized cratonic basin with a weak lower crust, and the rheological structure changes from a rather weak crème brûlée-1 regime into a strong jelly sandwich-1 regime, from the Altyn-Qaidam boundary to the northeastern corner of the WQB. Combined with the limited distribution of the LVZs and the strong crust-mantle decoupling revealed by the unmatched pattern between P axes (N20ºE) and SKS-wave splitting (N110°E), the crust-mantle mixing related to the under-thrusting of the Tarim Basin along the Altyn Tagh Fault is suggested as the primary tectonic dynamic inducement of the destabilized WQB craton, which weakens the lithospheric strength greatly and contributes to the shallow brittle deformation. However, the EQB performs as a typical cold and rigid cratonic basin characterized by a jelly sandwich-2 rheological regime, which is sufficiently strong to maintain crust-mantle coupling. The LVZs beneath the EQB revealed by recent wide-angle seismic profiles have only a small effect on the lithospheric strength drop. The EQB could be regarded as a solid basin, anchored in the NE Tibetan Plateau, which shows a strong resistance to NE extrusion of the weak plateau material, leading to clockwise rotation.

Active fault segments along the North Anatolian Fault system in the Sea of Marmara: implication for seismic hazard

A new analysis of high-resolution multibeam and seismic reflection data, collected during several oceanographic expeditions starting from 1999, allowed us to compile an updated morphotectonic map of the North Anatolian Fault below the Sea of Marmara. We reconstructed kinematics and geometries of individual fault segments, active at the time scale of 10 ka, an interval which includes several earthquake cycles, taking as stratigraphic marker the base of the latest marine transgression. Given the high deformation rates relative to sediment supply, most active tectonic structures have a morphological expression at the seafloor, even in presence of composite fault geometries and/or overprinting due to mass-wasting or turbidite deposits. In the frame of the right-lateral strike-slip domain characterizing the North Anatolian fault system, three types of deformation are observed: almost pure strike-slip faults, oriented mainly E–W; NE/SW-aligned axes of transpressive structures; NW/SE-oriented trans-tensional depressions. Fault segmentation occurs at different scales, but main segments develop along three major right-lateral oversteps, which delimit main fault branches, from east to west: (i) the transtensive Cinarcik segment; (ii) the Central (East and West) segments; and (iii) the westernmost Tekirdag segment. A quantitative morphometric analysis of the shallow deformation patterns observed by seafloor morphology maps and high-resolution seismic reflection profiles along the entire basin allowed to determine nature and cumulative lengths of individual fault segments. These data were used as inputs for empirical relationships, to estimate maximum expected Moment Magnitudes, obtaining values in the range of 6.8–7.4 for the Central, and 6.9–7.1 for the Cinarcik and Tekirdag segments, respectively. We discuss these findings considering analyses of historical catalogues and available paleoseismological studies for the Sea of Marmara region to formulate reliable seismic hazard scenarios.

Mantle upwelling beneath the Apennines identified by receiver function imaging

Magmatism, uplift and extension diffusely take place along collisional belts. Even though links between mantle dynamics and shallow deformation are becoming more evident, there is still poor understanding of how deep and surface processes are connected. In this work, we present new observations on the structure of the uppermost mantle beneath the Apennines belt. Receiver functions and seismic tomography consistently define a broad zone in the shallow mantle beneath the mountain belt where the shear wave velocities are lower than about 5% and the Vp/Vs ratio is higher than 3% than the reference values for these depths. We interpret these anomalies as a pronounced mantle upwelling with accumulation of melts at the crust-mantle interface, on top of which extensional seismicity responds to the crustal bending. The melted region extends from the Tyrrhenian side to the central part of the belt, with upraise of fluids within the crust favored by the current extension concentrated in the Apennines mountain range. More in general, mantle upwelling, following detachment of continental lithosphere, is a likely cause for elevated topography, magmatism and extension in post-collisional belts.

Stress and load-bearing structure analysis of the surrounding rock in a soft broken roadway

After roadway excavation, the surrounding rock is the main bearing structure for the original stress and forms a supporting ring with the supporting structure. The automatic load-bearing ring is the main interactional structure between the surrounding rock and the supporting system. The rock mass strength, stress, and deformation characteristics were considered to establish the theoretical calculation model for the softened broken surrounding rocks. In the model, the strength, stress, and deformations were analyzed for the elastic-plastic softening and broken zones. The analytical solution was obtained for the stress field, displacement field, and plastic damage range of the surrounding rocks. Considering the postpeak strength softening, shearing dilatation effects, original rock stress, and supporting resistance, the roadway bearing structure was divided into a shallow bearing layer and a deep bearing layer. Moreover, the crucial bearing layer was vital to the bearing capacity and stability of the layer. The influence of the shallow bearing layer, deep bearing layer, and crucial bearing layer on controlling the deformation of the surrounding rock was also analyzed. Based on the evolution characteristics of the shallow and deep buried surrounding rock, the weaker the residual strength of the surrounding rock was, the thicker the shallow and deep buried surrounding rock; the thicker the shallow and deep buried layers were, the greater the surface displacement of the surrounding rock; the crucial bearing layer was farther away from the roadway, and the roadway surface displacement was great; the greater the supporting resistance was, the thinner the shallow and deep bearing layers; and the crucial bearing layer was closer to the roadway, and the roadway deformation control was effective. Furthermore, based on the specific calculation example, the influence of the residual stress and supporting measures on the evolution characteristics of shallow and deep layers and roadway deformation was analyzed in detail.

Creep Along the Central San Andreas Fault From Surface Fractures, Topographic Differencing, and InSAR

AbstractImaging tectonic creep along active faults is critical for measuring strain accumulation and ultimately understanding the physical processes that guide creep and the potential for seismicity. We image tectonic deformation along the central creeping section of the San Andreas Fault at the Dry Lake Valley paleoseismic site (36.468°N, 121.055°W) using three data sets with varying spatial and temporal scales: (1) an Interferometric Synthetic Aperture Radar (InSAR) velocity field with an ~100‐km footprint produced from Sentinel‐1 satellite imagery, (2) light detection and ranging (lidar) and structure‐from‐motion 3‐D topographic differencing that resolves a decade of deformation over a 1‐km aperture, and (3) surface fractures that formed over the 3‐ to 4‐m wide fault zone during a drought from late 2012 to 2014. The InSAR velocity map shows that shallow deformation is localized to the San Andreas Fault. We demonstrate a novel approach for differencing airborne lidar and structure‐from‐motion topography that facilitates resolving deformation along and adjacent to the San Andreas Fault. The 40‐m resolution topographic differencing resolves a 2.5 ± 0.2 cm/yr slip rate localized to the fault. The opening‐mode fractures accommodate cm/yr of fault slip. A 90% ± 30% of the 1‐km aperture deformation is accommodated over the several meter‐wide surface trace of the San Andreas Fault. The extension direction inferred from the opening‐mode fractures and topographic differencing is 40°–48° from the local trend of the San Andreas Fault. The localization of deformation likely reflects the well‐oriented and mature fault.

Evaluation of Rainfall-Induced Landslides triggering using a multidisciplinary approach

In a large part of South America, slow landslides are triggered by extreme hydrometeorological conditions leading to, for instance, Rainfall-Induced Landslides – RILs. These RILs are common in urban areas and have a negative impact on the population and infrastructure development. Despite their importance, these events are little understood. We aimed at understanding the spatial distribution of RILs in the urban zone of Temuco, Chile (38.8°S, 72.6°W). The area has the typical hydrometeorological conditions of southern Chile. We conducted our assessment with a temporal analysis of shallow deformations, obtained by synthetic aperture radar interferometry (Sentinel 1 A/B). These shallow deformation rates were compared with satellite precipitation data (CHIRPS product) and electrical resistivity tomography (ERT). We identified active RIL-prone zones with deformation rates greater than 60 mm during the period 2014 to 2017, supporting theories of hydrometeorological control. Slow movements were observed in volcanic soils, suggesting the influence of their geotechnical characteristics. Our results can be extrapolated to the southern Andes (35°S-43°S), where a large number of volcanic-sedimentary units are susceptible to RILs. Finally, integration of our multidisciplinary approach will facilitate understanding of the local RIL dynamics, allowing a better risk management to decision-makers in South American and other developing countries.

Fluid seepage associated with slope destabilization along the Zambezi margin (Mozambique)

Evidences for active fluid seepages have been discovered along the Zambezi continental slope (offshore Southern Mozambique). These seepages are mostly associated with pockmarks which are aligned along a trend parallel to the slope and running closely upstream of the headwall scarp of a wide zone of slope destabilization. Fluid seepages are interpreted as a potential trigger for the slope destabilization. Acoustic anomalies within the water column have been interpreted as related to moderate bubble seepages mostly located outside and only punctually inside the destabilization zone. Exploration with the SCAMPI towed camera system in the widest pockmark (diameter 200 m wide) has shown fluid seepages associated to authigenic carbonate crusts and possibly bacterial mats. These fluid seepages are also associated to the presence of chemiosynthetic organisms (Vesicomyidae and Thyasiridae bivalves, Siboglinidae tubeworms). The sampled gas in the sediment corresponds mainly to CH4 of microbial origin, generated by hydrogenotrophic methanogenesis from a substrate of organic origin, i.e. a conventional process of genesis of microbial gas in the marine domain. No evidence for thermogenic gas was detected. Another type of pockmarks has been observed within the core of the slope destabilization zone. Most of these pockmarks are inactive in terms of fluid seepage at present time and are associated to carbonate buildups forming chimney geometries. They probably correspond to diagenetic chimneys of former fluid migration pathways that have been exhumed during the mass sliding and the surrounding depression are related to recurrent activity of strong lateral slope currents which have scoured the sediments around. The spatial organization of the slope destabilization features is considered as representative of the temporal evolution of the landslide giving information about the dynamics of slope instability processes. This proposed evolution started by scattered seepages of formation water with dissolved gas. Then free gas seepages appeared notably in the upper part of the slope. This was followed by progressive shallow deformation in the sediments downslope of the main gas seepages. Finally, the whole slope was destabilized forming imbricated landslides exhuming locally former diagenetic chimneys.

Extensional tectonics during the Tyrrhenian back‐arc basin formation and a new morpho‐tectonic map

AbstractWe present a new tectonic map focused upon the extensional style accompanying the formation of the Tyrrhenian back‐arc basin. Our basin‐wide analysis synthetizes the interpretation of vintage multichannel and single‐channel seismic profiles, integrated with modern seismic images, P‐wave velocity models, and high‐resolution morpho‐bathymetric data. Four distinct evolutionary phases of the Tyrrhenian back‐arc basin opening are further constrained, redefining the initial opening to Langhian/Serravallian time. Listric and planar normal faults and their conjugates bound a series of horst and graben, half‐graben and triangular basins. Distribution of extensional faults, active throughout the basin since Middle Miocene, allows us to define an arrangement of faults in the northern/central Tyrrhenian mainly related to a pure shear which evolved to a simple shear opening. At depth, faults accommodate over a Ductile‐Brittle Transitional zone cut by a low‐angle detachment fault. In the southern Tyrrhenian, normal, inverse and transcurrent faults appear to be related to a large shear zone located along the continental margin of the northern Sicily. Extensional style variation throughout the back‐arc basin combined with wide‐angle seismic velocity models allows to explore the relationships between shallow deformation, faults distribution throughout the basin, and crustal‐scale processes as thinning and exhumation.

Tectonics and seismicity in the Northern Apennines driven by slab retreat and lithospheric delamination

Understanding how long-term subduction dynamics relates to the short-term seismicity and crustal tec tonics is a challenging but crucial topic in seismotectonics. We attempt to address this issue by linking long-term geodynamic evolution with short-term seismogenic deformation in the Northern Apennines. This retreating subduction orogen displays tectonic and seismogenic behaviors on various spatiotemporal scales that also characterize other subduction zones in the Mediterranean area. We use visco-elasto-plastic seismo-thermo-mechanical (STM) modeling with a realistic 2D setup based on available geological and geophysical data. The subduction dynamics and seismicity are coupled in the numerical modeling, and driven only by buoyancy forces, i.e., slab pull. Our results suggest that lower crustal rheology and lithospheric mantle temperature modulate the crustal tectonics of the Northern Apennines, as inferred by previous studies. The observed spatial distribution of upper crustal tectonic regimes and surface displacements requires buoyant, highly ductile material in the subduction channel beneath the internal part of the orogen. This allows protrusion of the asthenosphere in the lower crust and lithospheric delamination associated with slab retreat. The resulting surface velocities and principal stress axes generally agree with present-day observations, suggesting that slab delamination and retreat can explain the dynamics of the orogen. Our simulations successfully reproduce the type and overall distribution of seismicity with thrust faulting events in the external part of the orogen and normal faulting in its internal part. Slab temperatures and lithospheric mantle stiffness affect the cumulative seismic moment release and spatial distribution of upper crustal earthquakes. The properties of deep, sub-crustal material are thus shown to influence upper crustal seismicity in an orogen driven by slab retreat, even though the upper crust is largely decoupled from the lithospheric mantle. Our simulations therefore highlight the effect of deep lower crustal rheologies, self-driven subduction dynamics and mantle properties in controlling shallow deformation and seismicity.

Tiltmeter data inversion to characterize a strain tensor source at depth: application to reservoir monitoring

Surface deformation measured by geodetic data is the sum of single-strain sources deforming at depth. A combination of volume changes from several analytical models (e.g. a point source or dislocation along a plane) can be used to model the different sources. However, solving for the best fit of volume variations, dislocations, position and orientation parameters of all sources is a nonlinear problem, and its solution is generally non-unique. This problem can be converted into a linear one by assimilating the sum of sources to a simplified model formed by three orthogonal planes of dislocations at fixed position and orientation. This strain source model is equivalent to having all neighbouring deformation sources contained in a small size volume. The determination of the strain tensor components can be performed by inverting geodetic data. Because of their high resolution, tiltmeters are well adapted to survey shallow deformation of volcanoes and geological reservoirs. However, they are known to display unknown long-term drift. We propose an approach to jointly estimate the temporal evolution of the strain source and time-dependent instrumental parameters. We verify the approach using synthetic data, giving confidence intervals for each component of the strain tensor. Finally, we link geological information to the internal deformation by interpreting the strain tensor as principal directions of deformation. This approach seems promising for the identification of fracture onset and fault reactivation in geothermal, hydrocarbon exploitations or volcanic systems.

Constraints on the shallow deformation around the Main Frontal Thrust in central Nepal from refraction velocities

The youngest fault system in the Himalayan orogeny is the Main Frontal Thrust (MFT), the frontal ramp of the Main Himalayan Thrust, which is expected to host the largest and most damaging earthquakes in Nepal. We characterize the upper few hundred meters below the surface across two MFT fault strands using ten high-resolution seismic profiles that we acquired in 2014 and 2015 with a 6-tonne Vibroseis source. We use first arrival picks from 625,416 seismic traces to derive P-wave seismic velocity models using a wavepath eikonal traveltime inversion method, and derive estimates of alluvium thickness and water table depth across these faults (the Patu and Bardibas thrusts), allowing us to constrain the subsurface geometry of the MFT.Our results show that 1) seismic velocities range from 255 to 3660 m/s, consistent with dry and saturated alluvium, and Siwalik bedrock; 2) low-velocity alluvium varies between ~20–50 and ~80–120 m thick in the hanging wall and footwall of the Bardibas thrust, respectively, corresponding to ~60–70 m of uplift of the hanging wall since deposition; 3) the two thrusts are soft-linked, and the western tip of the Bardibas thrust lies ~6 km west of its surface expression; 4) during the dry season, the water table is ~25–100 m higher in the hanging walls of the faults than in their footwalls, due to the larger thickness of permeable alluvium in the footwalls, and the water table shallows towards the east in the hanging wall of the Bardibas thrust; and 5) consistent with previous studies, the Patu thrust breaches the surface, while the Bardibas thrust is blind at Ratu Khola. These results also demonstrate that it should be possible to constrain the rate of uplift above the Bardibas thrust by drilling and dating sediments on both sides, which would complement existing measurements from terrace uplift.