Search

The stress regime patterns of high-seismically active regions within the western part of the India–Eurasia collision, spanning from 67° E to 83° E and 27° N to 39° N, are elucidated through analysis of 684 Focal Mechanism Solutions from 1962 to 2021. Eighteen seismically active zones used for the stress tensor inversion, are defined based on the spatial extent of the seismicity, the depth distribution of seismic events, focal mechanism studies, and seismotectonics of the region. The defined regimes are: (1) Sulaiman Ranges and Lobe Region, (2) Hindukush, (3) Pamir, (4) Nanga Parbat Syntaxis, (5) Hazara Syntaxis, (6) Kashmir–Zanskar region, (7) Kangra–Chamba, (8) Kinnaur and Kaurik–Chango fault zone (KCFZ), (9) Garhwal, (10) Kumaon, (11) Karakoram fault zone, and (12) Gozha–Ashikule fault zone. Seismicity is reported only in the crust or up to mid-crust in most of the regions, except for the Pamir and Hindukush, where the seismicity can be observed down to 160 and 280 km, respectively. We report a clockwise rotation of the maximum horizontal stress (SHmax) of about 42° and 21° in the Hindukush and Pamir regions, respectively. with increasing focal depths from north west to north. The region where major and strong earthquakes occur indicates pure compressive regimes. Most of the zones support transpressive and transtensional tectonics with a few zones by normal and strike-slip fault regimes. Regions like Nanga Parbat syntaxis, Kinnaur, KCFZ, and Zanskar are exceptions, where extensional and transformational tectonic features dominate. Plate convergence force has less effect on defining the stress regime in the Karakoram fault zone and Gozha–Ashikule regions, which display transtensional and pure extensional regimes, respectively. Underthrusting of the Indian plate through complex tectonics is indicated by dominant compression stresses with evidences of normal, strike-slip, and oblique fault mechanisms. [ABSTRACT FROM AUTHOR]

A high-precision and high-resolution vertical velocity for the Chinese mainland is obtained by integrating precise levelling and Global Navigation Satellite System data, using a Helmert joint adjustment method. The results show that the surface vertical rates range between −3.0 and 3.9 mm yr−1 with continuous deformation in most areas, except the obvious subsidence at the rates of −15.0 to −94.2 mm yr−1 induced by groundwater exploitation in the North China Plain. Particularly, the central and southern Tibet, Tien Shan, Alashan, Ordos, eastern Cathaysia and Northeast China uplift at the rates of 0.5–3.9 mm yr−1; the southeastern Tibetan Plateau, Sichuan basin and Yangtze block are dominated by surface subsidence at the rates of −3.3 to −0.5 mm yr−1. Furthermore, the vertical rates vary little between the eastern and western regions of the Chinese mainland despite their pronounced differences in horizontal deformations. The effects of gravity isostasy and non-tectonic factors, including the environmental mass loads, Glacier Isostatic Adjustment (GIA), poroelastic expansion/compression, and mining operations have partially contributed to the vertical deformation of the Chinese mainland. Overall, this velocity reflects the complicated deformation features induced by the multiple geodynamic processes of the Chinese mainland. These geodynamic processes include isostasy, orogenic processes and geothermal anomalies associated with slab subduction/plate collision. [ABSTRACT FROM AUTHOR]

The 2008 M w 7.9 Wenchuan earthquake ruptured the middle and northeastern segments of the Longmenshan Fault Zone (LMSFZ), and the 2013 M w 6.6 Lushan earthquake ruptured a 50-km-long fault in the southwestern segment. Subsequently, an M w 5.8 earthquake occurred approximately 10 km distant from the M w 6.6 Lushan earthquake. Therefore, the potential risk for larger earthquakes (> M w 6.6) on the southwestern section must be considered. This study collects the latest seismological and GPS data to construct an integrated seismotectonic model for the two neighbouring earthquake sequences. The model integrates the fault planes involved, the main shock rupture processes, the main shock-caused Coulomb stress perturbation, the aftershock distribution and the 3-D velocity structure of the source region, providing information for seismic risk evaluation. We find that three fault planes were involved, two of which were related to the main shocks, and the third was generated by the aftershocks following the first main shock. The main shocks were caused by nearly pure thrust faulting on the two planes with dip angles of approximately 45° and almost opposite dipping directions, thereby forming a conjugate angle of around 90°. The third plane was located between the two main shocks, approximately parallel to the second main shock's fault plane. Each of the main shocks primarily ruptured a single asperity, displaying simple time history. The Coulomb stress change of the first main shock facilitated the generation of the second main shock and the third fault plane, and the second main shock increased the stress on the first main shock's fault plane. The aftershocks were distributed within stratified materials by spatially varying interfaces and characterized by high Vp and Vs velocity, and a low Vp / Vs ratio. The atypical dip angles of approximately 45° for thrust faults and the conjugate angle of approximately 90° are indicative of high stress state. The single asperity rupture implies simple stress accumulation. The main shock-caused Coulomb stress change did not reduce the seismic risk in the source region. The varying interfaces are interpreted as a consequence of long-term horizontal compression. All of these characteristics suggest that the two earthquake sequences were generated by the breakage of three immature faults under strong compression by background stress, and the high stress state remains within the southwestern LMSFZ. [ABSTRACT FROM AUTHOR]

Rotation of rigid blocks within continental interiors far from the plate convergence boundary is an unusual process, the dynamics of which is not clear. The Tarim block, as a rigid Precambrian block in central Asia, is surrounded by the Tibetan–Pamir plateau to the south and Tian Shan mountains to the north. Numerous geophysical data suggested that the Tarim block experienced significant clockwise rotation in the Cenozoic. Meanwhile, contrasting deformation patterns and associated topographic responses were observed between the western–central and eastern Tian Shan. The relationship among the India–Asia collision, Tarim rotation, and Tian Shan responses are poorly constrained. Here, a series of large-scale, high-resolution 3-D numerical models were constructed. The model results reveal that the collision of the indenting Indian lithosphere with the southwestern rim of the Tarim block triggers clockwise rotation of the Tarim block. Further on, the Tarim rotation produces differential deformation responses along the strike of Tian Shan, that is convergence-induced higher compression and strong uplifting in central Tian Shan but divergence-induced less compression and moderate uplifting in eastern Tian Shan. Thus, the Tarim rotation serves as an indispensable linkage between the Tibetan plateau evolution and the far-field Tian Shan activation. [ABSTRACT FROM AUTHOR]

The Gonghe Basin in the northeast Tibetan Plateau presents significant potential for hot dry rock (HDR) geothermal resources. A 1990 M w 6.4 earthquake in the basin furthers the need for an improved understanding of its sedimentary structure. In this study, we utilize data from a dense seismic array of 88 short-period seismometers deployed at an interstation spacing of approximately 3 km to scrutinize the sedimentary structure of the Gonghe Basin. By analysing teleseismic P waveforms, we identify P -to- S converted waves (Ps wave) originating from the sedimentary basement. We then determine the delay time between the Ps waves and the direct P waves (P wave) through waveform cross-correlation. By integrating this delay time with empirical velocity structure models, HDR borehole data and results from teleseismic receiver function analysis, we derive a sediment thickness model of the Gonghe Basin for the Qabqa geothermal area. Our findings reveal a gradual increase in sediment thickness from around 500 m in the east to approximately 3000 m in the west, which is consistent with other geophysical surveys and borehole data. The thick sediments in the basin could potentially serve as an excellent thermal storage cover for HDR. The strong ground motion simulation using our sediment thickness model shows that thick sediments can amplify seismic waves, increasing the risk of seismic hazards. Moreover, our study indicates that the clear Ps waves can be effectively extracted to construct a dependable sediment thickness model using teleseismic P waves recorded by a short-period dense seismic array. [ABSTRACT FROM AUTHOR]

In this study, the focal mechanism solutions and source parameters of recent earthquakes that occurred in the northeastern region of India have been determined. The region has very complex tectonics as it is subjected to the compressional forces from all sides, due to the collision of the Indian Plate with the Eurasian, Burma and Tibetan plates. Waveform data from deployed broad-band seismographs (BBS) and strong motion accelerographs (SMA) in the northeastern India are used to determine the focal mechanism solutions and source parameters of moderate earthquakes, respectively. The estimated focal mechanisms are used to understand the existing stress field in the region. It is found that the Shillong-Plateau as well as the Indo-Burma subduction zone is dominated by the compressional tectonic regime, Mikir Hills and Bengal basin are dominated by the trans-tension tectonic regime, and the easternmost Himalayan region is dominated by the strike-slip tectonic regime. The maximum horizontal stress direction Sh max is also determined for above subregions. The direction of Sh max is southeast in the Bengal basin, northeast in Mikir Hills and Indo-Burma subduction zone whereas it is NNE in Shillong Plateau and SSW in the eastern Himalayas. The estimated stress drop value of the earthquakes in the region ranges from 2.11 to 23.89 MPa. The relationship between the source parameters and focal mechanisms is also explored. It is found that the earthquakes with a strike-slip mechanism have the highest average stress drop (7.05 MPa) followed by reverse (6.82 MPa) and normal (5.12 MPa) in the northeastern region of India. According to the examined data set, the stress drop is found to be dependent on the type of focal mechanism, seismic moment and hypocentral depths. The comparison of the results with the Kachchh intraplate region in western India shows earthquakes in Kachchh have larger mean stress drop for all types of mechanisms. In both intraplate and interplate regions of India, the stress drop of earthquakes depends on the type of focal mechanism solution. [ABSTRACT FROM AUTHOR]

Due to the far-field effect of the India–Eurasia collision, the Tien Shan orogenic belt has been undergoing reactivation and modification. Two end-member models of the geodynamic mechanisms are (1) surface uplift due to crustal shortening caused by lithospheric compression and (2) mountain formation resulting from thermal upwelling of asthenospheric mantle materials generated by lithospheric subduction. However, the topography along the Tien Shan orogenic belt changes significantly, and the deep structure and dynamic process are quite different beneath the Tien Shan orogenic belt from both geological and geophysical observations. Therefore, the reactivation and modification of the Tien Shan orogenic belt are likely influenced by both geodynamic mechanisms, which also generate various thermal anomalies in the crust. Seismic Lg-wave attenuation is very sensitive to crustal material composition and status and can point to the presence of partial melting within the crust resulting from mantle upwelling. In this study, we develop a high-resolution Lg-wave attenuation model between 0.05 and 10.0 Hz in Northwest China and use lateral attenuation variations to infer thermal structures in the crust. The central Tien Shan is characterized by prominent low- Q Lg anomalies, whereas relatively high- Q Lg distributions are imaged beneath the eastern and western Tien Shan. The surface uplift and crustal deformation are mostly related to the upwelling of hot mantle materials in the central Tien Shan and are likely induced by lithospheric compression in the eastern and western Tien Shan. However, low- Q anomalies are observed in the junction between the Pamir Plateau and western Tien Shan, indicating that the uplift in the south of the western Tien Shan is related to thermal subduction-induced upwelling and intrusion into the crust due to the collision between the Indian and Eurasian plates. The Kazakh Shield, characterized by pronounced high- Q Lg values, is likely a cold and hard terrane, and hence blocks the far-field effect of the India–Eurasia collision. [ABSTRACT FROM AUTHOR]

Earthquake moment tensors and centroid locations in the catalogue of the Global CMT (gCMT) project, formerly the Harvard CMT project, have become an essential resource for studying active global tectonics, used by many solid-Earth researchers. The catalogue's quality, long duration (1976–present), ease of access and global coverage of earthquakes larger than about Mw 5.5 have transformed our ability to study regional patterns of earthquake locations and focal mechanisms. It also allows researchers to easily identify earthquakes with anomalous mechanisms and depths that stand out from the global or regional patterns, some of which require us to look more closely at accepted interpretations of geodynamics, tectonics or rheology. But, as in all catalogues that are, to some extent and necessarily, produced in a semi-routine fashion, the catalogue may contain anomalies that are in fact errors. Thus, before re-assessing geodynamic, tectonic or rheological understanding on the basis of anomalous earthquake locations or mechanisms in the gCMT catalogue, it is first prudent to check those anomalies are real. The purpose of this paper is to illustrate that necessity in the eastern Himalayas and SE Tibet, where two earthquakes that would otherwise require a radical revision of current geodynamic understanding are shown, in fact, to have gCMT depths (and, in one case, also focal mechanism) that are incorrect. [ABSTRACT FROM AUTHOR]

The Daliang Shan is a tectonic unit that connects the active southeastern Tibetan Plateau with the stable South China block. As a newly generated (formed later than the Xianshuihe–Xiaojiang active fault system) seismotectonic zone induced by the Tibetan tectonics, the detailed constraints of the crustal deformation are central to the understanding of the kinematics and dynamics of the Tibetan expansion. This paper establishes and analyses a high-spatial resolution global positioning system (GPS) velocity field from a dense GPS network in this region. Our modelling results indicate that, in contrast to the equivalent sinistral strike-slip rate of approximately 5 mm yr–1 on the Anninghe–Zemuhe and Daliangshan faults, their inferred interseismic locking depth varies within a large range. The southern segment of the Anninghe Fault and the middle segment of the Daliangshan Fault have deep locking depths of ∼13 km, indicating that the seismic risk is high in these areas. In addition, the detectable counter-clockwise rotation rate of 0.35 ± 0.12° Myr–1 of the Mabian block makes a significant contribution of ∼50 per cent to the strike-slip motion on its boundary faults. This counter-clockwise rotation may be induced by a left-lateral shear gradient with southeastward motion relative to the South China-fixed reference frame, indicating the significance of a simple-shear pattern in exploring the kinematics of the encroachment of the Tibetan tectonics upon a stable block (craton). [ABSTRACT FROM AUTHOR]

The western Balkans occupy a region influenced by two major active tectonic processes: the collision between the Adriatic Region and the Dinarides in the west, and the extension of the Aegean Region and its surroundings as they move towards the Hellenic Trench. An understanding of the kinematics and dynamics of the western Balkans has significance for our understanding of continental tectonics in general, and is the object of this paper. The region is rich in observational data, with many well-studied earthquakes, good geodetic coverage by GNSS (Global Navigation Satellite System) and abundant exposure of active faulting and its associated geomorphology, especially within the Mesozoic carbonates that cover large sectors of the extensional areas. We first use such observations to establish the regional kinematic patterns, by which we mean a clarification of how active faulting achieves the motions observed in the deforming velocity field obtained from GNSS measurements. We then use geomorphological observations on the evolution of drainage systems to establish how kinematic and faulting patterns have changed and migrated during the Late Neogene-Quaternary. The kinematics, and its evolution, can then be used to infer characteristics of the dynamics, by which we mean the origin and effect of the forces that control the overall deformation. The principal influences are: (i) the distribution and evolution of gravitational potential energy (GPE) contrasts arising from crustal thickness variations and elevation, in particular the growth of topography by shortening in the Albanides–Hellenides mountain ranges and the high elevation of mainland Greece relative to the Mediterranean seafloor and (ii) the ability of the boundaries of the region, along the Adriatic coast and in the Hellenic Trench, to support the forces arising from those GPE contrasts. The evolution in space and time indicates an interaction between the anisotropic strength fabric of the upper crust associated with faulting, and the more distributed and smoother patterns of flow that are likely to characterize the ductile deformation of the lower, aseismic part of the lithosphere—both of which influence the deformation on the scale of 100–200 km. The persistent argument about whether continental deformation is best described by a continuum or by rigid-block motions is largely a matter of scale and particular location: both are influential in establishing the patterns we see. [ABSTRACT FROM AUTHOR]

We have studied the active and recent tectonics of New Guinea, using earthquake source modelling, analysis of gravity anomalies, seismic reflection profiles, and thermal and mechanical models. Our aim is to investigate the behaviour and evolution of a young continental deformation belt, and to explore the effects of lateral variations in foreland rheology on the deformation. We find that along-strike gradients in the lithosphere thickness of the southern foreland have resulted in correlated changes in seismogenic thickness, likely due to the effects on the temperature structure of the crust. The resulting variation in the strength of the foreland means that in the east, the foreland is broken through on thrust faults, whereas in the west it is relatively intact. The lack of correlation between the elevation of the mountain belt and the seismogenic thickness of the foreland is likely to be due to the time taken to thicken the crust in the mountains following changes in the rheology of the underthrusting foreland, as the thinned passive margin of northern Australia is consumed. The along-strike variation in whether the force exerted between the mountains and the lowlands is able to break the foreland crust enables us to estimate the effective coefficient of friction on foreland faults to be in the range of 0.01–0.28. We use force-balance calculations to show that the recent tectonic re-organization in western New Guinea is likely to be due to the development of increasing curvature in the Banda Arc, and that the impingement of continental material on the subduction zone may explain the unusually low force it exerts on western New Guinea. [ABSTRACT FROM AUTHOR]

As early as the Palaeozoic era, the Tianshan Mountains constituted a vast tectonic convergent zone within the Eurasian Continent; by this time, they had already experienced multiple intense crustal movements. The Indian Plate is currently advancing northwards by more than 50 mm yr–1; it affects the Tianshan Mountains even though they are located more than 1000 km away. However, many of the geological issues of the Tianshan Mountains, especially in eastern Tianshan, which remains an active tectonic zone today, are not yet fully understood. There have been few studies into the seismogenic characteristics of this intense seismic zone. Here, the full gravity gradient tensor (GGT) within the eastern Tianshan Mountains was calculated using a stable spatial frequency domain algorithm, based on high-spatial resolution (7.2 arcsec or approximately 200 m) Global Gravity Model plus (GGMplus) surface gravity data. In addition, the features of the geological structures in the eastern Tianshan Mountains were interpreted using different combinations of GGT components. Moreover, the possibilities of using different combinations of GGT components to identify the distributions and strikes of faults was discussed regarding the potential earthquake rupture risks in this study area. The results show that (1) the distributions and strikes of the main fault zones in the study area are highly consistent with the linear geological features within the mountains, as extracted from a combination of different GGT components, including terrain effects. (2) Removing the terrain effects revealed that the different components of GGT (derived from the complete Bouguer gravity anomaly) vary gently, showing few linear geological features that could be extracted. This implies that satellite gravity data are not sufficient to analyse the characteristics of small underground geological structures (spatial scales less than ∼10 km). (3) Comparing the spatial distributions of fault traces depicted at different ages with the linear geological features of the study area revealed that discovering more unknown faults increased the coincidence ratio of linear geological features to fault traces in the area (from 27 to 40 per cent). In conclusion, the findings of this study represent a valuable reference for further revealing the seismogenic characteristics of this area, where there is a lack of detailed surface fault structure measurements due to its inaccessibility. The proposed method could also recognize faults in other similar areas with harsh conditions that are not feasible for ground-based surveys, such as forests and glaciers. [ABSTRACT FROM AUTHOR]

The northeastern Indian region is characterized by complex lithospheric structure that developed due to collision between the Indian and Eurasian tectonic plates, in the north, and to subduction beneath the Burmese arc, in the east. We report results from joint modelling of Ps and Sp receiver functions and Rayleigh wave group velocity dispersion curves in which a broad search for acceptable models is performed via simulated annealing. We identify three tectonic domains, the Shillong plateau, Brahmaputra valley and Indo-Burma convergence zone (IBCZ), sampled by teleseismic earthquake data recorded by nine broad-band seismic stations. Our results reveal that the region's thinnest crust lies beneath the Shillong plateau, where it increases slightly from the plateau's eastern edge to its centre and reaches a maximum at the western edge of the plateau. Crustal Vp / Vs ratios range between 1.69 and 1.75 for the Shillong plateau, which is consistent with a felsic composition. Deeper Moho depths beneath the Brahmaputra valley, adjacent to the northern front of the Shillong plateau, may be due to the flexure of Indian lithosphere subducting beneath Asia. Low velocity zones are indicated at ∼5–10 km depth beneath the Brahmaputra valley, which may have been developed by NE–SW trending compressional stresses from the collision at the Himalayan arc and subduction at the Burmese arc. The crust is thickest in Kohima, beneath the Naga thrust in the IBCZ, where a high velocity zone is observed for both Vp and Vs at a depth of 25–40 km. This anomaly may be associated with a high velocity slab, trending N–NE to S–SW, that comprises the subducting Indian lithosphere in the IBCZ. [ABSTRACT FROM AUTHOR]

As the largest and most active intracontinental orogenic belt on Earth, the Tien Shan (TS) is a natural laboratory for understanding the Cenozoic orogenic processes driven by the India–Asia collision. On 19 January 2020, a M w 6.1 event stuck the Kalpin region, where the southern frontal TS interacts with the Tarim basin. To probe the local ongoing orogenic processes and potential seismic hazard in the Kalpin region, both interseismic and instantaneous deformation derived from geodetic observations are used in this study. With the constraint of interseismic global navigation satellite system (GNSS) velocities, we estimate the décollement plane parameters of the western Kalpin nappe based on a 2-D dislocation model, and the results suggest that the décollement plane is nearly subhorizontal with a dip of ∼3° at a depth of 24 km. Then, we collect both Sentinel-1 and ALOS-2 satellite images to capture the coseismic displacements caused by the 2020 Kalpin event, and the interferometric synthetic aperture radar (InSAR) images show a maximum displacement of 7 cm in the line of sight near the epicentral region. With these coseismic displacement measurements, we invert the source parameters of this event using a finite-fault model. We determine the optimal source mechanism in which the fault geometry is dominated by thrust faulting with an E–W strike of 275° and a northward dip of 11.2°, and the main rupture slip is concentrated within an area 28.0 km in length and 10.3 km in width, with a maximum slip of 0.3 m at a depth of 6–8 km. The total released moment of our preferred distributed slip model yields a geodetic moment of 1.59 × 1018 N |$ \cdot $| m, equivalent to M w 6.1. The contrast of the décollement plane depth from interseismic GNSS and the rupture depth from coseismic InSAR suggests that a compression still exists in the Kalpin nappe forefront, which is prone to frequent moderate events and may be at risk of a much more dangerous earthquake. [ABSTRACT FROM AUTHOR]

The increasing density of geodetic measurements makes it possible to map surface strain rate in many zones of active tectonics with unprecedented spatial resolution. Here we show that the strain tensor rate calculated from GPS in the India–Asia collision zone represents well the strain released in earthquakes. This means that geodetic data in the India–Asia collision zone region can be extrapolated back in time to estimate strain buildup on active faults, or the kinematics of continental deformation. We infer that the geodetic strain rates can be assumed stationary through time on the timescale needed to build up the elastic strain released by larger earthquakes, and that they can be used to estimate the probability of triggering earthquakes. We show that the background seismicity rate correlates with the geodetic strain rate. A good fit is obtained assuming a linear relationship (⁠|$\dot{N} = \lambda \ \cdot \dot{\epsilon }$|⁠ , where |$\dot{N}$| is the density of the rate of M w ≥ 4 earthquakes, |$\dot{\epsilon }$| is strain rate and λ = 2.5 ± 0.1 × 10−3 m−2), as would be expected from a standard Coulomb failure model. However, the fit is significantly better for a non-linear relationship (⁠|$\dot{N} = \gamma _1 \cdot \dot{\epsilon }^{\gamma _2}$| with γ1 = 2.5 ± 0.6 m−2 and γ 2 = 1.42 ± 0.15). The b -value of the Gutenberg–Richter law, which characterize the magnitude–frequency distribution, is found to be insensitive to the strain rate. In the case of a linear correlation between seismicity and strain rate, the maximum magnitude earthquake, derived from the moment conservation principle, is expected to be independent of the strain rate. By contrast, the non-linear case implies that the maximum magnitude earthquake would be larger in zones of low strain rate. We show that within areas of constant strain rate, earthquakes above M w 4 follow a Poisson distribution in time and and are uniformly distributed in space. These findings provide a framework to estimate the probability of occurrence and magnitude of earthquakes as a function of the geodetic strain rate. We describe how the seismicity models derived from this approach can be used as an input for probabilistic seismic hazard analysis. This method is easy to automatically update, and can be applied in a consistent manner to any continental zone of active tectonics with sufficient geodetic coverage. [ABSTRACT FROM AUTHOR]

Rapid urban growth has led to large population densities in foreland basin regions, and therefore a rapid increase in the number of people exposed to hazard from earthquakes in the adjacent mountain ranges. It is well known that earthquake-induced ground shaking is amplified in sedimentary basins. However, questions remain regarding the main controls on this effect. It is, therefore, crucial to identify the main controls on earthquake shaking in foreland basins as a step towards mitigating the earthquake risk posed to these regions. We model seismic-wave propagation from range-front thrust-faulting earthquakes in a foreland-basin setting. The basin geometry (depth and width) and source characteristics (fault dip and source-to-basin distance) were varied, and the resultant ground motion was calculated. We find that the source depth determines the amount of near-source ground shaking and the basin structure controls the propagation of this energy into the foreland basin. Of particular importance is the relative length scales of the basin depth and dominant seismic wavelength (controlled by the source characteristics), as this controls the amount of dispersion of surface-wave energy, and so the amplitude and duration of ground motion. The maximum ground motions occur when the basin depth matches the dominant wavelength set by the source. Basins that are shallow compared with the dominant wavelength result in low-amplitude and long-duration dispersed waveforms. However, the basin structure has a smaller effect on the ground shaking than the source depth and geometry, highlighting the need for understanding the depth distribution and dip angles of earthquakes when assessing earthquake hazard in foreland-basin settings. [ABSTRACT FROM AUTHOR]

We present a multiscale decomposition model for analysing the widespread deformation and complex dynamics mechanism of the southeastern margin of the Tibetan Plateau. We first carry out simulations to test the performance of the model by synthetic Global Navigation Satellite System (GNSS) data for a network, whose geometry distribution is same with actual Crust Movement Observation Network of China in the southeast edge of the Tibetan Plateau. We then estimate the multiscale velocity field and strain rate field using observed data with 240 stations between 1999 and 2007. Results show that the deformation field at the different spatial scale corresponds to the different tectonic driving force and strain accumulation coverage. Tectonic background of the Tibetan Plateau, extrusion of the Sichuan–Yunnan rhombic block, rigid block motion, interior deformation of the block and segmented motion of fault zones are separately displayed in the corresponding scale, which provides a strong support keys to the understanding of the controversial deformation models for the tectonic evolution and uplift mechanisms of the Tibetan Plateau. Our results also show that scale 7 with the coverage of about 50 km can better reflect the seismogenic background of mid-strong earthquakes (> M s 5.0) of the region. The convergence area of high compression and inflation strain or the high-value area of maximum shear strain at scale 7 is obviously correlated with seismic activities, including 2008 Wenchuan M s 8.0 earthquake, no significant pre-seismic deformation across fault in other scales. Our research provides a new insight into the nature of tectonic deformation and the seismic risk assessment of the southeastern margin of the Tibetan Plateau. [ABSTRACT FROM AUTHOR]

We investigate how different crustal models can affect the stress field, velocities and associated deformation in the India–Eurasia collision zone. We calculate deviatoric stresses, which act as deformation indicators, from topographic load distribution and crustal heterogeneities coupled with density driven mantle convection constrained by tomography models. We use three different crustal models, CRUST2.0, CRUST1.0 and LITHO1.0 and observe that these models have different crustal thickness and densities. As a result, gravitational potential energy (GPE) calculated based on these densities and crustal thicknesses differ between these models and so do the associated deviatoric stresses. For GPE only models, LITHO1.0 provides a better constraint on deformation as it yields the least misfit (both orientation and relative magnitude) with the surface observations of strain rates, lithospheric stress, plate motions and earthquake moment tensors. However, when the stresses from GPE are added to those associated with mantle tractions arising from density-driven mantle convection, the coupled models in all cases provide a better fit to surface observations. The N–S tensional stresses predicted by CRUST2.0 in this area get reduced significantly due to addition of large N–S compressional stresses predicted by the tomography models S40RTS and SAW642AN leading to an overall strike-slip regime. On the other hand, the hybrid models, SINGH_S40RTS and SINGH_SAW that are obtained by embedding a regional P -wave model, Singh et al. in global models of S40RTS and SAW642AN, predict much lower compression within this area. These hybrid models provide a better constraint on surface observations when coupled with CRUST1.0 in central Tibet, whereas the combined LITHO1.0 plus mantle traction model provides a better fit in some other areas, but with a degradation of fit in central Tibet. [ABSTRACT FROM AUTHOR]

Orogenic plateaus can exist in a delicate balance in which the buoyancy forces due to gravity acting on the high topography and thick crust of the plateau interior are balanced by the compressional forces acting across their forelands. Any shortening or extension within a plateau can indicate a perturbation to this force balance. In this study, we present new observations of the kinematics, morphology and slip rates of active normal faults in the South Peruvian Altiplano obtained from field studies, high-resolution DEMs, Quaternary dating and remote sensing. We then investigate the implications of this faulting for the forces acting on the Andes. We find that the mountains are extending ∼NNE–SSW to ∼NE–SW along a normal fault system that cuts obliquely across the Altiplano plateau, which in many places reactivates Miocene-age reverse faults. Radiocarbon dating of offset late Quaternary moraines and alluvial fan surfaces indicates horizontal extension rates across the fault system of between 1 and 4 mm yr–1—equivalent to an extensional strain rate in the range of 0.5–2 × 10−8 1 yr–1 averaged across the plateau. We suggest the rate and pattern of extension implies there has been a change in the forces exerted between the foreland and the Andes mountains. A reduction in the average shear stresses on the sub-Andean foreland detachment of ≲4 MPa (20–25 per cent of the total force) can account for the rate of extension. These results show that, within a mountain belt, the pattern of faulting is sensitive to small spatial and temporal variations in the strength of faults along their margins. [ABSTRACT FROM AUTHOR]

Numerous V-shaped conjugate strike-slip fault systems distributed between the Lhasa block and the Qiangtang block serve as some of the main structures accommodating the eastward motion of the Tibetan Plateau. The Beng Co-Dongqiao conjugate fault system is a representative section, and determining its tectonic environment is a fundamental issue for understanding the dynamic mechanism of the V-shaped conjugate strike-slip fault systems throughout central Tibet. In this paper, we investigate the deformation rates of the Beng Co-Dongqiao conjugate faults using 3 yr of SAR data from both ascending and descending tracks of Sentinel-1 satellites. Only interferograms with a long temporal baseline were used to increase the proportion of the deformation signals. The external atmospheric delay product and the InSAR stacking strategy were employed to reduce various errors in the large-spatial-coverage Sentinel-1 data. The InSAR results revealed that the fault-parallel deformation velocities along the eastern and western segments of the Beng Co fault are 5 ± 1 mm/yr and 2.5 ± 1 mm/yr, respectively. The second invariant of the horizontal strain rates shows that the accumulated strain is centered on the eastern segment of the Beng Co Fault and Gulu rift. The velocity and strain rate fields show that the Anduo-Peng Co faults may be paired with the Beng Co fault to form a new conjugate system and the tectonic transformation between the Beng Co fault and Gulu rift. These results can better explain the tectonic deformation environment of the Beng Co-Dongqiao conjugate fault system and provide insights on the crustal dynamics throughout the entire plateau interior. [ABSTRACT FROM AUTHOR]

Focal mechanism solutions and their predicted stress pattern can be used to investigate tectonic deformation in seismically active zones and contribute to understanding and constraining the kinematic patterns of the outward growth and uplift of the Tibetan Plateau. Herein, we determined the focal mechanisms of 398 earthquakes in Northeast Tibet recorded by the China National Seismic Network (CNSN) by using the cut-and-paste method. The results show that the earthquakes predominately exhibited thrust and strike-slip faulting mechanisms with very few normal events. We then combined the derived focal mechanisms with global centroid moment tensor (GCMT) catalogue solutions and previously published solutions to predict the regional distribution of the stress field through a damped linear inversion. The inversion results show that most of region is dominated by a thrust faulting regime. From the southern East Kunlun fault in the west to the northern Qilian Mountains along the Altyn Tagh fault (ATF), the maximum compression axis rotates slightly clockwise; farther to the south of the Haiyuan fault in the east, there is an evident clockwise rotation of the maximum compression axis, especially at the eastern end of the Haiyuan fault. In the Qilian Mountains, the axis of the compressive stress orientation approximately trends NE–SW, which does not markedly differ from the direction of India–Eurasia convergence, emphasizing the importance of the compressive stress in reflecting the remote effects of this continental collision. The overall spatial pattern of the principal stress axes is closely consistent with the GPS-derived horizontal surface velocity. A comparison of the stress and strain rate fields demonstrated that the orientations of the crustal stress axes and the surface strain axes were almost identical, which indicates that a diffuse model is more suitable for describing the tectonic characteristics of Northeast Tibet. Additionally, the compressive stress orientation rotated to ENE–WSW in the northern Qilian Mountains along the ATF and to ENE–WSW or E–W along the eastern part of the Haiyuan fault and its adjacent area to the south, highlighting the occurrence of strain partitioning along large left-lateral strike-slip faults or the lateral variation of crustal strength across these faults. Combining geodetic, geological and seismological results, we suggest that a hybrid model incorporating both the diffuse model associated with shortening and thickening of the upper crust and the asthenospheric flow model accounting for the low-velocity zone in the middle-lower crust may reflect the primary mode of crustal deformation in Northeast Tibet. [ABSTRACT FROM AUTHOR]

The activities of frontal thrusts in the northern Qilian Shan are critical for understanding the deformation of the Qilian Shan and the northeastern Tibetan Plateau. In this study, we estimate the slip rate of the active Fodongmiao–Hongyazi thrust along the northern margin of the Qilian Shan. High-resolution satellite imagery interpretations and detailed field investigations suggest that the fault displaced late Pleistocene terraces and formed fresh prominent north-facing fault scarps. To quantify the slip rate of the fault, we measured the displacements along the fault scarps using an unmanned aerial vehicle system and dated the displaced geomorphic surfaces using optically stimulated luminescence (OSL) and 14C methods. The vertical slip rate of the fault is estimated at 1.0 ± 0.3 mm yr−1 for the western segment. The slip rates for two branches in the eastern segment are 0.3 ± 0.1 and 0.6 ± 0.1 mm yr−1. Using a fault dip of 40 ± 10°, we constrain the corresponding shortening rates to 1.4 ± 0.5 and 1.2 ± 0.4 mm yr−1, respectively. The rates are consistent with values over different timescales, which suggests steady rock uplift and northeastward growth of the western Qilian Shan. Crustal shortening occurs mainly on the range-bounding frontal thrust. [ABSTRACT FROM AUTHOR]

The formation of fold-thrust belts at convergent margins is a dynamic process. Accretion of weak sediments to the front of the overriding plate results in crustal thickening and continued flexural subsidence of the underthrusting plate. Fold-thrust belts are often treated as a Coulomb wedge having self-similar geometries with a critical taper, and either a rigid or isostatically compensated base. In this paper, we build upon this work by developing a new dynamic model to investigate both the role of the thickness and material properties of the incoming sediment, and the flexure in the underthrusting plate in controlling the behaviour and evolution of fold-thrust belts. Our analysis shows that the evolution of fold-thrust belts can be dominated by either gravitational spreading or vertical thickening, depending on the relative importance of sediment flux, material properties and flexure. We apply our model to the Makran accretionary prism and the Indo–Burman Ranges, and show that for the Makran flexure must be considered in order to explain the dip of the sediment–basement interface from seismic reflection profiles. In the Indo–Burman Ranges, we show that incoming sediment thickness has a first-order control on the variations in the characteristics of the topography from north to south of the Shillong Plateau. [ABSTRACT FROM AUTHOR]

Free air gravity anomalies are now available from Tibetan Plateau and surrounding regions. They show that the Plateau itself is isostatically compensated, and that its elevation is supported by its margins acting as barriers. Its interior has an elastic thickness of less than 4 km. The main features of the Plateau can be understood by using a simple lithostatic flexural model, which can account for the paired gravity anomalies and the thrust faulting on its margins. The resulting estimates of the magnitudes of flexural forces agree with those from other observations. The thick crust that underlies the Plateau generates radiogenic heat that greatly reduces its viscosity and accounts for its weakness. As Tibet flows south over cold Indian lithosphere, heat is conducted downwards, heating up the upper mantle beneath the Moho. Combined with isostatic compensation, the resulting thermal expansion then produces a gradient in crustal thickness, from about 75 km in the south to 65 km in the north. Though this lithostatic model provides a framework for understanding many of the features of Tibet, it does not account for the difference in the dynamical behaviour between the north and south of the Plateau. [ABSTRACT FROM AUTHOR]

This article presents a list of forthcoming papers related to geophysics to be published in Geophysical Journal International. A test of laboratory based rheological parameters of olivine from an analysis of late Cenozoic convective removal of mantle lithosphere beneath the Sierra Nevada, California, is written by P. Molnar and C.H. Jones. A one-way wave equation for modelling variations in seismic waveforms due to elastic anisotropy by D.A. Angus, C.J. Thomson and R.G. Pratt is also to be published. The list also includes a research paper on sea-level change and remagnetization of continental shelf sediments off New Jersey by H. Oda and M. Torii.

We assess the stress field that drives lithospheric deformation, focusing on the active New Zealand plate-boundary zone between the Australian and Pacific plates. Here there is a rich database for the horizontal velocity field, and also crustal structure used to derive gravitational potential energy (GPE). We solve the stress balance equations, in the context of a thin sheet model of a viscously deforming lithosphere, for characteristic deviatoric stresses and viscosities, defined as the integral of these with depth (divided by the layer thickness), using the stress method of Flesch et al. where the input parameters are the fields of strain rate and GPE that apply to the sheet. Synthetic tests show that the stress method is able to resolve the stress field even with high (20 per cent) levels of noise in the input strain rates, and the mean stress is a very robust feature of the inversions, regardless of noise levels. We invert for the stress and viscosity fields in New Zealand, calculating the field of GPE from topography/bathymetry and crustal data, and the field of strain rates from either a long-term (multimillennial) velocity field inferred from the rate and pattern of Quaternary faulting or a short-term (decadal) velocity field directly observed with Global Positioning System (GPS) measurements. In addition, we consider the effect of shear stresses on the subducted plate interface along the Hikurangi Margin (5–15 MPa), or regionally. We explore the effect of GPE on the inversion results by calculating these for a range of deforming layers (i.e. thin sheet with 35–150 km thickness), in effect sampling the lithospheric strength in the crust and mantle. The results show that the derived stress magnitudes (square root second invariant of the stress deviator) are in the range 0–35 MPa, with mean values of 13 ± 1 MPa for all models, comparable to typical earthquake stress drops. Gravitationally induced stresses account for approximately half of the full deviatoric stress. Effective characteristic viscosities are 0.5–5 × 1021 Pa s in the deforming zone, with an approximate inverse relation between strain rate and viscosity, most likely controlled by thermal structure and/or lithology. [ABSTRACT FROM AUTHOR]

Magnitudes of differential stress in the lithosphere, especially in the crust, are still disputed. Earthquake-based stress drop estimates indicate median values <10 MPa whereas the lateral variation of gravitational potential energy per unit area (GPE) across significant relief indicates stress magnitudes of ca. 100 MPa in average across a 100 km thick lithosphere between the Indian lowland and the Tibetan plateau. These standard GPE -based stress estimates correspond to membrane stresses because they are associated with a deformation that is uniform with depth. We show here with new analytical results that lateral variations in GPE can also cause bending moments and related bending stresses of several hundreds of MPa. Furthermore, we perform 2-D thermomechanical numerical simulations (1) to evaluate estimates for membrane and bending stresses based on GPE variations, (2) to quantify minimum crustal stress magnitudes that are required to maintain the topographic relief between Indian lowland and Tibetan plateau for ca. 10 Ma and (3) to quantify the corresponding relative contribution of crustal strength to the total lithospheric strength. The numerical model includes viscoelastoplastic deformation, gravity and heat transfer. The model configuration is based on density fields from the CRUST1.0 data set and from a geophysically and petrologically constrained density model based on in situ field campaigns. The numerical results indicate that values of differential stress in the upper crust must be > ca. 180 MPa, corresponding to a friction angle of ca. 10° to maintain the topographic relief between lowland and plateau for >10 Ma. The relative contribution of crustal strength to total lithospheric strength varies considerably laterally. In the region between lowland and plateau and inside the plateau the depth-integrated crustal strength is approximately equal to the depth-integrated strength of the mantle lithosphere. Simple analytical formulae predicting the lateral variation of depth-integrated stresses agree with numerically calculated stress fields, which show both the accuracy of the numerical results and the applicability of simple, rheology-independent, analytical predictions to highly variable, rheology-dependent stress fields. Our results indicate that (1) crustal strength can be locally equal to mantle lithosphere strength and that (2) crustal stresses must be at least one order of magnitude larger than median stress drops in order to support the plateau relief over a duration of ca. 10 Ma. [ABSTRACT FROM AUTHOR]

Isostatic gravity anomalies provide a measure of the Earth's gravity field free from the gravitational attractions of the topography and its isostatic compensation, most commonly represented by a variation in the depth of a compensating density contrast, for example the Moho. They are used by both geodesists and geophysicists alike, though often for different purposes. Unfortunately though, the effect of subsurface loading on the lithosphere renders transfer function (admittance) methods unusable when surface and subsurface loads coexist. Where they exist, subsurface loads are often expressed in the Bouguer anomaly but not in the topography, and it is shown here that this phase disconnect cannot be faithfully represented by either real- or complex-valued analytic admittance functions. Additionally, many studies that employ the isostatic anomaly ignore the effects of the flexural rigidity of the lithosphere, most often represented as an effective elastic thickness (T e), and assume only Airy isostasy, i.e. surface loading of a plate with zero elastic thickness. The consequences of such an omission are studied here, finding that failure to account for flexural rigidity and subsurface loading can result in (1) over- or underestimates of both inverted Moho depths and dynamic topography amplitude, and (2) underestimates of the size of topographic load that can be supported by the plate without flexure. An example of the latter is shown over Europe. Finally, it is demonstrated how low values of the isostatic anomaly variance can actually be biased by these anomalies having low power at the long wavelengths while still possessing high power at middle to short wavelengths, compared to the corresponding Bouguer anomaly power spectrum. This will influence the choice of best-fitting isostatic model if the model is chosen by minimization of the isostatic anomaly standard deviation. [ABSTRACT FROM AUTHOR]

Six large magnitude earthquakes in Central Asia which occurred at the end of the 19th century were recorded on early magnetographs in Great Britain. Scalar seismic moment estimates of the 1911 Chon-Kemin, the 1902 Atushi and the 1907 Karatag earthquakes in Central Asia were recently determined by historical seismogram modelling. For those events, we find agreement between moment magnitudes estimated from seismograms and from magnetograms. This supports the assumption of linear scaling of magnetogram amplitudes as function of M 0, which we then use to estimate the moment magnitudes for earlier large-magnitude events, that is, the 1885 Belovodskoe, 1887 Verny and 1889 Chilik earthquakes. The magnetometer data imply that the Chilik earthquake had M W7.9, slightly smaller than the Chon-Kemin event with M W8.0. The Verny earthquake, however, for which we estimate M W7.7, is likely larger than listed in catalogues (M 7.3). Similarly, we find a larger magnitude M W7.6 (instead of the previous M 6.9) for the Belovodskoe earthquake, but this remains uncertain due to measurement imprecision. [ABSTRACT FROM AUTHOR]

We demonstrate that a continental interior reverse fault is deforming by aseismic creep, presently, and likely also in the long term. The Karkara Rangefront Fault, part of the larger Main Terskey Front, forms the northern boundary of the high Terskey Tien Shan in southeastern Kazakhstan and is a mature structure with evidence for high slip rates throughout the late Cenozoic. Combining field studies with a satellite stereo-image derived digital elevation model (2 m resolution), we map a series of fluvial terraces along the rangefront which are uplifted by up to ∼300 m above the present river level. Radiocarbon ages from one catchment constrain the ages of the lowermost two terraces to be ∼4–5  ka and ∼10–15  ka, consistent with prominent, regionally extensive terraces observed elsewhere in the Tien Shan. Based on conservative estimates for the fault dip under the displaced terraces, we estimate a slip rate along the fault plane of 3.5 |$^{+1.7}_{-0.4}$|  mm yr–1 on the Karkara Rangefront Fault and a further >0.8 mm yr–1 on a fold structure in the Kegen basin that we infer is driven by a detachment from the main rangefront. We therefore estimate a minimum shortening rate across the rangefront of 1.1–3.3  mm yr–1. Elastic modelling of the regional GPS velocity field suggests that the fault is presently creeping at ∼3  mm yr–1 (horizontal shortening), consistent with the upper limit of our Late Quaternary slip rate estimate. This is the fastest known slip rate in the northern Tien Shan and the only individual structure resolved in the regional velocity field. At present the fault is accumulating minimal strain, and there is evidence in the geomorphology that this creep is sustained in the long term, but whether or not it is also capable of generating earthquakes requires further study. [ABSTRACT FROM AUTHOR]