Present-day crustal deformation in China constrained by global positioning system measurements

Paleogeography of the southwestern Ordos Basin and exhumation history of the Liupan Shan area

The Liupan Shan, situated on the southwestern Ordos Basin, underwent Cenozoic uplift caused by the India-Asia collision and subsequent northeastward tectonic movements. The strata in this region record both the paleogeography of the southwestern Ordos Basin and the uplift history of the Tibetan Plateau. However, past studies have rarely analyzed the strata comprehensively, resulting in overlooked information within them. We present the first detrital zircon (U-Th)/He (ZHe) data from the Lower Cretaceous deposits of the Liupan Shan and the Upper Miocene red clay of Chaona, constraining a maximum burial depth of less than 6-7 km at 60 Ma for the Liupan area. By integrating zircon U-Pb ages, paleocurrent data, and sedimentary facies, we found a primary ZHe age peak (∼210 Ma) indicating that the Lower Cretaceous deposits were sourced from the Qinling Orogenic Belt, while the red clay was likely originated from the Liupan Shan. Although the Chaona strata are not preserved in the Liupan Shan, our results suggest that the Miocene Chaona section is crucial for understanding the Late Cenozoic uplift of the Liupan Shan. The distinct stratigraphic ages and source areas reflect the complex depositional and tectonic history of the region. Thermal modeling results revealed three stages of tectonic events (Late Cretaceous, Late Paleocene-Early Eocene, and Late Cenozoic), which are linked to plate interaction. Our findings offer new insights into the long-term tectonic evolution of the Tibetan Plateau.

Observed hydrological changes associated with active tectonic blocks before three consecutive earthquakes in Qinghai, China

In the past 2 years, three earthquakes of magnitude 6.0 and above occurred consecutively in Qinghai province, China, i.e., the 22 May 2021 M_s7.4 Maduo, 8 January 2022 M_s6.9 Menyuan, and 26 March 2022 M_s6.0 Delingha earthquakes. The hydrological observation instruments set up by the China Earthquake Administration allow us to study the dynamic processes in the well-aquifer systems during the establishment of criticality. Particularly, the observations played an important role in the prediction of the 8 January 2022 M_s6.9 Menyuan earthquake that was approved by the People's Government of Qinghai province. This work presents the hydrological data recorded by 7 stations to show the short-term anomalies before these earthquakes. To explore the performance of the hydrological observations in detecting earthquakes that occurred on different active tectonic blocks, we calculate the relative amplitudes of the pre-seismic changes. Results indicate that markedly pre-seismic change is found if the observation station and the detection earthquake are on the same block, and moderate change is found if they are on the adjacent blocks, while the precursor is hard to be identified if they are on the separated blocks. The variations in the hydrological responses may be caused by the strength weakening (or dilatancy) of source media. And the increased volumes in the crust can be evidenced by the changes in the geodetic time series in the same neighborhoods and during the same period, augmenting stress loading between the blocks.

Two-Decade GNSS Observation Processing and Analysis with the New IGS Repro3 Criteria: Implications for the Refinement of Velocity Field and Deformation Field in Continental China

Extensive observation collection, unified and rigorous data processing, and accurate construction of the station motion model are the three essential elements for the accuracy and reliability of the Global Navigation Satellite System (GNSS) velocity field. GNSS data reprocessing not only can weaken the influence of untrue nonlinear site signals caused by imperfect models but also can eliminate the displacement offset caused by frame transformation, solution strategy, and model change. Based on the new repro3 criteria of the International GNSS Service (IGS), we process rigorously GNSS observations of continental China from the period 2000 to 2020 to refine GNSS station secular velocities and analyze the present-day crustal deformation in continental China. The main contributions of this work included the followings. Firstly, the repro3 algorithm and model are used to uniformly and rigorously process the two-decade GNSS historical observations to obtain more reliable GNSS coordinate time series with mm-level precision. Combined with the historical records of major earthquakes in continental China, we build a GNSS time series model considering nonlinear factors (velocity, offset, period, co-seismic/post-seismic deformation) to extract GNSS horizontal velocity field whose root mean square (RMS) mean is 0.1 mm/a. Secondly, the GNSS horizontal grid velocity field in continental China is interpolated using the gpsgridder method (the minimum radius is set to 16, and the Poisson’s ratio is set to 0.5). Estimation and analysis of the crustal strain rate solution lead to the conclusion that the strain degree in West China (the high strain region is mainly located in the Qinghai Tibet Plateau and Tianshan Mountains) is much more intense than that in the east (the main strain rate is less than 5 nstrain/year). In addition, most strong earthquakes in the Chinese mainland occurred on active blocks and their boundary faults with large changes in the GNSS velocity field and strain field. Then, an improved K-means++ clustering analysis method is proposed to divide active blocks using GNSS horizontal velocity field. Furthermore, different relative motion models of different blocks are constructed using the block division results. Among them, the Eurasian block has the lowest accuracy (the RMS of residual velocity in the east and north directions are 5.60 and 9.65 mm/a, respectively), and the China block 7 has the highest accuracy (the RMS mean of relative velocity in the east and north directions are 2.60 and 2.65 mm/a, respectively). More observations (2260+ sites), longer time (20 years), and updated criteria (Repro3) are to finely obtain the GNSS velocity field in continental China, and depict crustal deformation and active block with the gpsgridder and improved K-means++ methods.

Strain Field Features and Three-Dimensional Crustal Deformations Constrained by Dense GRACE and GPS Measurements in NE Tibet

The continuing impact between the Eurasia Plate and India results in the thickening and shortening of the N-S Tibetan Plateau. There has been strong tectonic movement along the boundary of the zones of deformation of the NE corner of the Tibetan plateau (NET) since the new tectonic period, with its dynamic mechanisms remaining controversial. Here, we use observations of 39 Continuous Global Positioning System (CGPS) gauges and 451 Crustal Movement Observation Network of China (CMONOC) campaign-mode stations to detect the three-dimensional deformation of the crust in the NET. Improved processing procedures are implemented to strengthen the patterns of strain throughout the NET. The principal component analysis (PCA) technique is introduced to decompose the time series into spatial eigenvectors and principal components (PCs), and the first three PCs are used to estimate and rectify common mode errors (CMEs). In addition, GRACE observations are used to detect deformation changes that account for non-tidal oceanic mass loading, hydrological loading, and surface pressure. The rectified deformation of the crust indicates the anisotropic nature of both the subsidence and uplift, and that the highest uplift rate of the Longmen Shan fault uplift reaches 7.13 ± 0.53 mm/yr. Finally, the horizontal velocity is further used to enumerate the strain rates throughout the NET. The results show that the shear band retained property in line with the strike-slip fault along the Altyn Tagh fault, the Qilian Shan faults, the Haiyuan fault, the West Qinling fault, the East Kunlun fault, and the Longmen Shan fault. In addition, the results further indicate that the whole NET shows a strong relationship with the mean principal rates of horizontal shortening strain. Extension and compression of the crust reasonably describe its sinking and uplifting.

Depth-based sparse bundle adjustment

It is demonstrated in this paper that due to error model inconsistency, a certain degree of accuracy loss would be incurred to the estimated parameters when the traditional bundle adjustment method is directly applied to the scenario where a fraction of observations is implicitly error free (e.g., the reference image points in commonly used least squares matching refinement). To this end, a depth-based object point model and corresponding depth-based sparse bundle adjustment method are proposed in this paper, in which the position of an object point is represented by its 1D depth relative to its reference image. A corresponding projection model is derived, the sparse block structures of normal equations are studied depending on whether there are shared image parameters to be optimized or not, and corresponding sparse solutions of the normal equations and parameter covariance matrices are derived. Compared with the traditional sparse bundle adjustment method, simulated experiments demonstrate that our method matches the error model of the target scenario, and thus can avoid further accuracy loss. Moreover, both simulated and real data experiments demonstrate that our method can effectively improve computational efficiency.

The Current Crustal Vertical Deformation Features of the Sichuan–Yunnan Region Constrained by Fusing the Leveling Data with the GNSS Data

This study uses the least squares collocation method to fuse the leveling vertical deformation velocity in the Sichuan–Yunnan region with the GNSS observations of this region from 320 stations in the China Crustal Movement Observation Network (CMONOC) and the China Continental Tectonic Environment Monitoring Network (CMTEMN) from 1999 to 2017. Such fusion is to improve the accuracy of the vertical deformation rates in large spatial scales. The fused vertical deformation results show that: (1) the fused deformation field has a uniform spatial distribution, and shows detailed change characteristics of key regions; (2) the current vertical crustal motion in this region is featured by the contemporaneous occurrence of crustal compression, shortening and uplift and basin extensional subsidence; (3) most areas in this region experience uplifts, as the lateral push of the Qinghai–Tibet Plateau was blocked by the Sichuan Basin. The areas on the northwest side of the Longmenshan fault and the Lijiang-Xiaojinhe fault are dominated by uplifts, with the velocity of 1.5 mm/a–5.5 mm/a, and the region on the southeast side has slight uplifts, with the velocity of 1.0 mm/a–1.5 mm/a; (4) many areas have high gradient vertical deformation, especially the region close to the Wenshan fault and on the two sides of the Yarlung Zangbo fault that has the value of 3.0–4.0 × 10−8/a, deserving further attention to be paid to the long-term earthquake hazards.

Strain-Rates from GPS Measurements in the Ordos Block, China: Implications for Geodynamics and Seismic Hazards

A number of devastating earthquakes have occurred around the Ordos Block in recent history. For the purpose of studying where the next major event will occur surrounding the Ordos Block, much work has been done, particularly in the investigation of the Earth’s surface strain rates based on GPS measurements. However, there exist striking differences between the results from different authors although they used almost the same GPS data. Therefore, we validated the method for the calculation of GPS strain rates developed by Zhu et al. (2005, 2006) and found that the method is feasible and has high precision. With this approach and the updated GPS data, we calculated the strain rates in the region around the Ordos Block. The computed results show that the total strain rates in the interior of the Block are very small, and the high values are mainly concentrated on the peripheral zones of the Ordos Block and along the large-scale active faults, such as the Haiyuan fault, which are closely aligned to the results by geological and geophysical observations. Additionally, the strain rate results demonstrated that all rifted grabens on the margin of the Ordos Block exhibit extensional deformation. Finally, based on the strain rate, seismicity, and tectonic structures, we present some areas of high earthquake risk surrounding the Ordos Block in the future, which are located on the westernmost of the Weihe Graben, both the east and westernmost of the Hetao Graben, and in the middle of the Shanxi Graben. Hence, this work is significant in contributing to a better understanding of the geodynamics and seismic hazard assessment.

Analysis of Crustal Movement and Deformation in Mainland China Based on CMONOC Baseline Time Series

In this paper, we propose a method for the analysis of tectonic movement and crustal deformation by using GNSS baseline length change rates or baseline linear strain rates. The method is applied to daily coordinate solutions of continuous GNSS stations of the Crustal Movement Observation Network of China (CMONOC). The results show that: (a) The baseline linear strain rates are uneven in space, which is prominent in the Tianshan, Sichuan-Yunnan, Qinghai-Tibet Plateau, and Yanjing areas, with a maximum value of 1 × 10−7 a−1, and about two orders smaller in the South China block, the Northeast block, and the inner area of the Tarim basin, where the average baseline linear strain rates are 1.471 × 10−9 a−1, 2.242 × 10−9 a−1, and 3.056 × 10−9 a−1, respectively; (b) Active crustal deformation and strong earthquakes in the Xinjiang area are mainly located in the north and south sides of the Tianshan block; the compression deformations both inside the Tarim block and in the southern Tianshan fault zone are all increasing from east to west, and the Tarim block is not a completely “rigid block”, with the shrinkage rate in the west part at about 1~2 mm/a; (c) The principal directions of crustal deformation in the Xinjiang, Tibet, and Sichuan-Yunnan regions are generally in the north—south compression and east—west extension, indicating that the collision and wedging between the Indian and Eurasian plates are still the main source of tectonic movements in mainland China.

Adjoint traveltime tomography unravels a scenario of horizontal mantle flow beneath the North China craton

The North China craton (NCC) was dominated by tectonic extension from late Cretaceous to Cenozoic, yet seismic studies on the relationship between crust extension and lithospheric mantle deformation are scarce. Here we present a three dimensional radially anisotropic model of NCC derived from adjoint traveltime tomography to address this issue. We find a prominent low S-wave velocity anomaly at lithospheric mantle depths beneath the Taihang Mountains, which extends eastward with a gradually decreasing amplitude. The horizontally elongated low-velocity anomaly is also featured by a distinctive positive radial anisotropy (V_SH > V_SV). Combining geodetic and other seismic measurements, we speculate the presence of a horizontal mantle flow beneath central and eastern NCC, which led to the extension of the overlying crust. We suggest that the rollback of Western Pacific slab likely played a pivotal role in generating the horizontal mantle flow at lithospheric depth beneath the central and eastern NCC.

A possible precursor prior to the Lushan earthquake from GPS observations in the southern Longmenshan

Global Positioning System (GPS) stations installed in and around the epicenter of the Lushan earthquake (Mw 6.7), which occurred almost 5 years after the 2008 Wenchuan earthquake, recorded preseismic deformation corresponding to the Lushan earthquake within the southern Longmenshan thrust belt. A half-space dislocation model is used to simulate the theoretical values of the postseismic displacements caused by the 2008 Wenchuan earthquake, and after transforming the reference frame and filtering the GPS displacement time series, the theoretical and observed GPS values are compared to identify the geodetic anomaly preceding the Lushan earthquake. The abnormal extent of this geodetic anomaly decreases with increasing epicentral distance for each GPS site. This geodetic signal reflects preslip along a locked section of the 2013 seismogenic fault, which caused the accumulation of elastic strain energy until the faulting strength was overcome, thereby generating the Lushan earthquake. Hence, this anomaly might be used as an observable and identifiable precursor to forecast an impending earthquake within a period of less than two and half years before its occurrence.

Normal Faulting in the 2020 Mw 6.2 Yutian Event: Implications for Ongoing E–W Thinning in Northern Tibet

Extensional earthquakes in the Tibetan Plateau play an important role in the plateau’s orogenic evolution and cause heavy seismic hazard, yet their mechanisms remain poorly known, in particular in harsh northern Tibet. On 25 June 2020, a Mw 6.2 earthquake struck Yutian, Xinjiang, offering us a rare chance to gain insights into its mechanism and implications in the Tibetan extension. We used both descending and ascending Sentinel-1 images to generate coseismic deformation associated with this event, which indicates a typical extensional mechanism with a maximum subsidence displacement of 25 cm and minor uplift. The causative fault constrained with interferometric synthetic aperture radar (InSAR) data based on a finite fault model suggests that the fault plane has a strike of 186.4° and westward dip of 64.8°, and the main rupture is concentrated at a depth of 3.6–10.8 km with a peak slip of 0.85 m. Our source model indicates that the 2020 Yutian event ruptured an unknown high-angle blind normal fault with N–S striking. The total released geodetic moment yields 2.69 × 1018 N·m, equivalent to Mw 6.23. We used dense interseismic global positioning system (GPS) measurements to reveal an approximate 7 mm/yr extensional motion in the Yutian region, but it still does not seem large enough to support high local seismicity for normal events within 12 years, i.e., Mw 7.1 in 2008, Mw 6.2 in 2012, and this event in 2020. Combined with Coulomb stress change modeling, we speculate that the seismicity in Yutian is related to the lower lithospheric dynamics.

Analysis of the Potential Contributors to Common Mode Error in Chuandian Region of China

Common mode error (CME) in Chuandian region of China is derived from 6-year continuous GPS time series and is identified by principal component analysis (PCA) method. It is revealed that the temporal behavior of the CME is not purely random, and contains unmodeled signals such as nonseasonal mass loadings. Its spatial distribution is quite uniform for all GPS sites in the region, and the first principal component, uniformly distributed in the region, has a spatial response of more than 70%. To further explore the potential contributors of CME, daily atmospheric mass loading and soil moisture mass loading effects are evaluated. Our results show that ~15% of CME can be explained by these daily surface mass loadings. The power spectral analysis is used to assess the CME. After removing atmospheric and soil moisture loadings from the CME, the power of the CME reduces in a wide range of frequencies. We also investigate the contribution of CME in GPS filtered residuals time series and it shows the Root Mean Squares (RMSs) of GPS time series are reduced by applying of the mass loading corrections in CME. These comparison results demonstrate that daily atmosphere pressure and the soil moisture mass loadings are a part of contributors to the CME in Chuandian region of China.

Geodetic Model of the 2017 Mw 6.5 Mainling Earthquake Inferred from GPS and InSAR Data

On 17 November 2017, a Mw 6.5 earthquake occurred in Mainling County, Nyingchi City, China. The epicenter was located in the Namche Barwa region of the eastern Himalayan syntaxis. Here, we have derived coseismic deformation from Global Positioning System (GPS) data and ascending Sentinel-1A Synthetic Aperture Radar (SAR) data. Based on a joint inversion of the two datasets, we obtained the coseismic slip distribution along a curved, northeast trending, and high-angle (dip angle of 75°) thrust fault. Our results show that the seismic moment release was 7.49 × 1018 N∙m, corresponding to a moment magnitude of Mw 6.55. The maximum slip was 1.03 m and the main rupture zone extended to a 12 km depth. The earthquake may have been related to the release of strain accumulated during the subduction of the Indian plate beneath the Eurasian continent. We identified a high strain rate and low b-values around the epicentral area before the earthquake, indicating that the earthquake was nucleated under a high strain/stress state. The data indicate two regions, southwest and southeast to the epicenter (the eastern Main Himalaya Thrust and northern end of the Sagaing fault), which remain under high stress/strain conditions and pose a significant seismic hazard.

Common Mode Component and Its Potential Effect on GPS-Inferred Three-Dimensional Crustal Deformations in the Eastern Tibetan Plateau

Surface and deep potential geophysical signals respond to the spatial redistribution of global mass variations, which may be monitored by geodetic observations. In this study, we analyze dense Global Positioning System (GPS) time series in the Eastern Tibetan Plateau using principal component analysis (PCA) and wavelet time-frequency spectra. The oscillations of interannual and residual signals are clearly identified in the common mode component (CMC) decomposed from the dense GPS time series from 2000 to 2018. The newly developed spherical harmonic coefficients of the Gravity Recovery and Climate Experiment Release-06 (GRACE RL06) are adopted to estimate the seasonal and interannual patterns in this region, revealing hydrologic and atmospheric/nontidal ocean loads. We stack the averaged elastic GRACE-derived loading displacements to identify the potential physical significance of the CMC in the GPS time series. Interannual nonlinear signals with a period of ~3 to ~4 years in the CMC (the scaled principal components from PC1 to PC3) are found to be predominantly related to hydrologic loading displacements, which respond to signals (El Niño/La Niña) of global climate change. We find an obvious signal with a period of ~6 yr on the vertical component that could be caused by mantle-inner core gravity coupling. Moreover, we evaluate the CMC’s effect on the GPS-derived velocities and confirm that removing the CMC can improve the recognition of nontectonic crustal deformation, especially on the vertical component. Furthermore, the effects of the CMC on the three-dimensional velocity and uncertainty are presented to reveal the significant crustal deformation and dynamic processes of the Eastern Tibetan Plateau.

Geodetic Constraints on the Crustal Deformation along the Kunlun Fault and Its Tectonic Implications

This study focuses on the crustal deformation and interseismic fault coupling along the strike-slip Kunlun fault, northern Tibet, whose western segment ruptured in the 2001 Mw 7.8 Kokoxili earthquake. We first integrated published Global Positioning System (GPS) velocity solutions and calculated strain rate fields covering the Kunlun fault. Our results show abnormally high post-earthquake strain rate values across the ruptures; furthermore, these exceed those in pre-earthquake data. Together with two tracks of interferometric synthetic aperture radar (InSAR) observations (2003–2010) and position time-series data from two continuous GPS sites, we show that the postseismic deformation of the Kokoxili earthquake may continue up to 2014; and that the postseismic transients of the earthquake affect the 2001–2014 GPS velocity solutions. We then processed the GPS data observed in 2014–2017 and obtained a dense interseismic velocity field for the northern Tibet. Using a fault dislocation model in a Bayesian framework, we estimated the slip rates and fault coupling on the Kunlun fault in 1991–2001 and 2014–2017. Results show an increase of slip rates and eastward migration of high fault coupling on the Kunlun fault after 2001. We propose the temporal variations are a result of the eastward accelerating movement, as a whole, of the Bayanhar block, whose boundaries were decoupled by several large earthquakes since 1997. Moreover, our results show the accumulated elastic strains along the Alake Lake-Tuosuo Lake segments could be balanced by an Mw 7.4–7.7 earthquake.

Present-Day Deformation of the Gyaring Co Fault Zone, Central Qinghai–Tibet Plateau, Determined Using Synthetic Aperture Radar Interferometry

Because of the constant northward movement of the Indian plate and blockage of the Eurasian continent, the Qinghai–Tibet Plateau has been extruded by north–south compressive stresses since its formation. This has caused the plateau to escape eastward to form a large-scale east–west strike-slip fault and a north–south extensional tectonic system. The Karakorum–Jiali fault, a boundary fault between the Qiangtang and Lhasa terranes, plays an important role in the regional tectonic evolution of the Qinghai–Tibet Plateau. The Gyaring Co fault, in the middle of the Karakoram–Jiali fault zone, is a prominent tectonic component. There have been cases of strong earthquakes of magnitude 7 or greater in this fault, providing a strong earthquake occurrence background. However, current seismic activity is weak. Regional geodetic observation stations are sparsely distributed; thus, the slip rate of the Gyaring Co fault remains unknown. Based on interferometric synthetic aperture radar (InSAR) technology, we acquired current high-spatial resolution crustal deformation characteristics of the Gyaring Co fault zone. The InSAR-derived deformation features were highly consistent with Global Positioning System observational results, and the accuracy of the InSAR deformation fields was within 2 mm/y. According to InSAR results, the Gyaring Co fault controlled the regional crustal deformation pattern, and the difference in far-field deformation on both sides of the fault was 3–5 mm/y (parallel to the fault). The inversion results of the back-slip dislocation model indicated that the slip rate of the Gyaring Co fault was 3–6 mm/y, and the locking depth was ~20 km. A number of v-shaped conjugate strike-slip faults, formed along the Bangong–Nujiang suture zone in the central and southern parts of the -Tibet Plateau, played an important role in regional tectonic evolution. V-shaped conjugate shear fault systems include the Gyaring Co and Doma–Nima faults, and the future seismic risk cannot be ignored.

The interactions of fault patterns and stress fields during active faulting in Central North China Block: Insights from numerical simulations

The interaction of active faults as a factor affecting the mechanisms of large earthquakes has been observed in many places. Most aftershock and clustering earthquake sequences do not recur on the main seismogenic fault but are controlled by fault interactions with adjacent seismic structures. Four groups of conceptual models were generated in this study to determine how the geometry of the seismogenic faults controls the distributions of stress fields and earthquakes. The influences of the fault length ratio, center distance, overlap ratio, echelon distance and fault opening angle were considered in a 2D viscoelastic model. The results indicate that the interaction in the slipping zone is larger when collinear interacting faults are more closely positioned, with one fault lengthening. For noncollinear faults, the interaction is stronger as the inner tips pass each other, which impedes their growth after some degree of overlap. Additionally, fault interaction at the slipping zone becomes stronger as the opening angle approaches 180°. We further generated a 3D viscoelastic model of fault interactions in Central North China Block and applied the finite element method to analyze the relationship between distributions of earthquakes and fault geometry. The calculated results reveal well-matched higher stress and maximum shear strain concentrations in the southern part of the Fen-wei Graben Zone than in other zones in Central North China Block, which can be explained by the longer faults, shorter center distances, shorter overlap lengths and larger opening angles. The stress distributions and fault interactions should be considered in long-term seismic hazard assessment in these zones.

Bimodal seismicity in the Himalaya controlled by fault friction and geometry

There is increasing evidence that the Himalayan seismicity can be bimodal: blind earthquakes (up to Mw ~ 7.8) tend to cluster in the downdip part of the seismogenic zone, whereas infrequent great earthquakes (Mw 8+) propagate up to the Himalayan frontal thrust. To explore the causes of this bimodal seismicity, we developed a two-dimensional, seismic cycle model of the Nepal Himalaya. Our visco-elasto-plastic simulations reproduce important features of the earthquake cycle, including interseismic strain and a bimodal seismicity pattern. Bimodal seismicity emerges as a result of relatively higher friction and a non-planar geometry of the Main Himalayan Thrust fault. This introduces a region of large strength excess that can only be activated once enough stress is transferred upwards by blind earthquakes. This supports the view that most segments of the Himalaya might produce complete ruptures significantly larger than the 2015 Mw 7.8 Gorkha earthquake, which should be accounted for in future seismic hazard assessments.

There is increasing evidence that the seismicity of large Himalayan earthquakes can be bimodal, but the underlying mechanisms are unclear. Here, the authors present a model and show that the bimodal seismicity results from a relatively higher friction and a non-planar geometry of the Himalayan megathrust.

Spatiotemporal Evolution of Postseismic Deformation Following the 2001 Mw7.8 Kokoxili, China, Earthquake from 7 Years of Insar Observations

The 2001 Mw7.8 Kokoxili earthquake, which occurred in the north Tibetan Plateau, ruptured ~400 km of the westernmost portion of the Kunlun fault and produced significant time-dependent postseismic deformation over a large area around the rupture zone and nearby regions. To analyze the postseismic deformation features along different sections of the coseismic surface rupture, we describe the total cumulative postseismic deformation near the center of the rupture and produce velocity maps for the whole observation period and six sub-periods, using InSAR observations (ENVISAT/ASAR, 2003–2010) on five descending tracks. The results indicate that the postseismic deformation is asymmetrically distributed across the fault over a very broad area of ~300 km × 500 km. The south side of the fault exhibits larger displacements and a wider area of deformation that is steadily decaying from near-field to far-field, while the north side displays a narrow, rapidly diminishing deformation field. The maximum cumulative displacement in 2003–2010 reaches up to ~45–60 mm and the LOS peak-to-trough average velocity offset in 2003–2010 reaches ~13–16 mm/yr at ~92.5°E. The short-term postseismic velocity estimates in the six sub-periods reflect significant spatial variation and temporal differences on different sections. Motions to the south of the two ends of the rupture zone show more rapid velocity decay compared to near the main central rupture zone. The time- and distance-dependent timeseries of postseismic surface displacement reveal exponential decay in the near-field and a nearly linear trend in the far-field of the fault.

Fault Slip Rates and Seismic Moment Deficits on Major Faults in Ordos Constrained by GPS Observation

The Ordos Block, surrounded by numerous active faults, is a relatively rigid but dangerous area with many strong historical earthquakes. We derive the block rotation velocity and fault slip rates in this area by using GPS data recorded from 1999 to 2007 and implementing an elastic block model. Instead of assuming vertical faults, as did most previous studies in and around Ordos, we use an improved method to invert for the fault dip angles and construct a closed 3-D fault system in our inversion. The predicted slip rates range from <1 mm/yr to ~ 10 mm/yr. Our results are roughly consistent with geological and other geodetic observations. Using the estimated slip rates, we also calculate the cumulative seismic moment due to fault locking and the released moment from historical earthquake catalogues. A comparison of the two quantities indicates that the Hetao Rift has an unreleased seismic moment equal to a M_w 7.9 earthquake, which is also indicated by frequent earthquakes above M6 after 1900.

Effects of Spatiotemporal Filtering on the Periodic Signals and Noise in the GPS Position Time Series of the Crustal Movement Observation Network of China

Analysis of Global Positioning System (GPS) position time series and its common mode components (CMC) is very important for the investigation of GPS technique error, the evaluation of environmental loading effects, and the estimation of a realistic and unbiased GPS velocity field for geodynamic applications. In this paper, we homogeneously processed the daily observations of 231 Crustal Movement Observation Network of China (CMONOC) Continuous GPS stations to obtain their position time series. Then, we filtered out the CMC and evaluated its effects on the periodic signals and noise for the CMONOC time series. Results show that, with CMC filtering, peaks in the stacked power spectra can be reduced at draconitic harmonics up to the 14th, supporting the point that the draconitic signal is spatially correlated. With the colored noise suppressed by CMC filtering, the velocity uncertainty estimates for both of the two subnetworks, CMONOC-I (≈16.5 years) and CMONOC-II (≈4.6 years), are reduced significantly. However, the CMONOC-II stations obtain greater reduction ratios in velocity uncertainty estimates with average values of 33%, 38%, and 54% for the north, east, and up components. These results indicate that CMC filtering can suppress the colored noise amplitudes and improve the precision of velocity estimates. Therefore, a unified, realistic, and three-dimensional CMONOC GPS velocity field estimated with the consideration of colored noise is given. Furthermore, contributions of environmental loading to the vertical CMC are also investigated and discussed. We find that the vertical CMC are reduced at 224 of the 231 CMONOC stations and 170 of them are with a root mean square (RMS) reduction ratio of CMC larger than 10%, confirming that environmental loading is one of the sources of CMC for the CMONOC height time series.

Uplift, climate and biotic changes at the Eocene–Oligocene transition in south-eastern Tibet

The uplift history of south-eastern Tibet is crucial to understanding processes driving the tectonic evolution of the Tibetan Plateau and surrounding areas. Underpinning existing palaeoaltimetric studies has been regional mapping based in large part on biostratigraphy that assumes a Neogene modernization of the highly diverse, but threatened, Asian biota. Here, with new radiometric dating and newly collected plant-fossil archives, we quantify the surface height of part of the south-eastern margin of Tibet in the latest Eocene (∼34 Ma) to be ∼3 km and rising, possibly attaining its present elevation (3.9 km) in the early Oligocene. We also find that the Eocene–Oligocene transition in south-eastern Tibet witnessed leaf-size diminution and a floral composition change from sub-tropical/warm temperate to cool temperate, likely reflective of both uplift and secular climate change, and that, by the latest Eocene, floral modernization on Tibet had already taken place, implying modernization was deeply rooted in the Palaeogene.

Growth of the northeastern margin of the Tibetan Plateau by squeezing up of the crust at the boundaries

In classic orogenic models, the mountain range is underlain by a deep crustal root. Here we present the crustal and upper mantle structures along two receiver function profiles across Qilian, an orogen experiencing recent growth at the northern margin of the Tibetan plateau. Opposite to an expected crustal root beneath the orogen, the Moho beneath Qilian is arch-like, shallower beneath the center and deepens by up to 10 km beneath its southern and northern boundaries. Additional velocity interfaces sub-parallel to the Moho are observed in the lower crust of the basins south of Qilian, which we interpret as the top of a mechanically strong lower crust thrusting several tens of kilometers underneath Qilian. In the north, the small lateral offset between the surface and mantle traces of the thrust system reveals a steep boundary, indicating that the North China cratonic crust acts as a strong resistance to the northward growth of the plateau, forcing the development of the left-lateral strike-slip Haiyuan fault south of the northern Qilian suture. The young Qilian orogen thus has been rising and growing progressively from the boundaries to the center, squeezed up by more rigid tectonic blocks in the north and south.

Lithospheric foundering and underthrusting imaged beneath Tibet

Long-standing debates exist over the timing and mechanism of uplift of the Tibetan Plateau and, more specifically, over the connection between lithospheric evolution and surface expressions of plateau uplift and volcanism. Here we show a T-shaped high wave speed structure in our new tomographic model beneath South-Central Tibet, interpreted as an upper-mantle remnant from earlier lithospheric foundering. Its spatial correlation with ultrapotassic and adakitic magmatism supports the hypothesis of convective removal of thickened Tibetan lithosphere causing major uplift of Southern Tibet during the Oligocene. Lithospheric foundering induces an asthenospheric drag force, which drives continued underthrusting of the Indian continental lithosphere and shortening and thickening of the Northern Tibetan lithosphere. Surface uplift of Northern Tibet is subject to more recent asthenospheric upwelling and thermal erosion of thickened lithosphere, which is spatially consistent with recent potassic volcanism and an imaged narrow low wave speed zone in the uppermost mantle.

The timing and mechanism of uplift of the Tibetan plateau continues to be a source of debate. Here, the authors present a new tomographic model revealing a T-shaped high wave speed structure beneath South-Central Tibet and interpret this an upper-mantle remnant from lithospheric foundering.

Contemporary crustal movement of southeastern Tibet: Constraints from dense GPS measurements

The ongoing collision between the Indian plate and the Eurasian plate brings up N-S crustal shortening and thickening of the Tibet Plateau, but its dynamic mechanisms remain controversial yet. As one of the most tectonically active regions of the world, South-Eastern Tibet (SET) has been greatly paid attention to by many geoscientists. Here we present the latest three-dimensional GPS velocity field to constrain the present-day tectonic process of SET, which may highlight the complex vertical crustal deformation. Improved data processing strategies are adopted to enhance the strain patterns throughout SET. The crustal uplifting and subsidence are dominated by regional deep tectonic dynamic processes. Results show that the Gongga Shan is uplifting with 1-1.5 mm/yr. Nevertheless, an anomalous crustal uplifting of ~8.7 mm/yr and negative horizontal dilation rates of 40-50 nstrain/yr throughout the Longmenshan structure reveal that this structure is caused by the intracontinental subduction of the Yangtze Craton. The Xianshuihe-Xiaojiang fault is a major active sinistral strike-slip fault which strikes essentially and consistently with the maximum shear strain rates. These observations suggest that the upper crustal deformation is closely related with the regulation and coupling of deep material.

Seasonal Mass Changes and Crustal Vertical Deformations Constrained by GPS and GRACE in Northeastern Tibet

Surface vertical deformation includes the Earth’s elastic response to mass loading on or near the surface. Continuous Global Positioning System (CGPS) stations record such deformations to estimate seasonal and secular mass changes. We used 41 CGPS stations to construct a time series of coordinate changes, which are decomposed by empirical orthogonal functions (EOFs), in northeastern Tibet. The first common mode shows clear seasonal changes, indicating seasonal surface mass re-distribution around northeastern Tibet. The GPS-derived result is then assessed in terms of the mass changes observed in northeastern Tibet. The GPS-derived common mode vertical change and the stacked Gravity Recovery and Climate Experiment (GRACE) mass change are consistent, suggesting that the seasonal surface mass variation is caused by changes in the hydrological, atmospheric and non-tidal ocean loads. The annual peak-to-peak surface mass changes derived from GPS and GRACE results show seasonal oscillations in mass loads, and the corresponding amplitudes are between 3 and 35 mm/year. There is an apparent gradually increasing gravity between 0.1 and 0.9 μGal/year in northeast Tibet. Crustal vertical deformation is determined after eliminating the surface load effects from GRACE, without considering Glacial Isostatic Adjustment (GIA) contribution. It reveals crustal uplift around northeastern Tibet from the corrected GPS vertical velocity. The unusual uplift of the Longmen Shan fault indicates tectonically sophisticated processes in northeastern Tibet.

Seasonal Hydrological Loading in Southern Tibet Detected by Joint Analysis of GPS and GRACE

In southern Tibet, ongoing vertical and horizontal motions due to the collision between India and Eurasia are monitored by large numbers of global positioning system (GPS) continuous and campaign sites installed in the past decade. Displacements measured by GPS usually include tectonic deformation as well as non-tectonic, time-dependent signals. To estimate the regional long-term tectonic deformation using GPS more precisely, seasonal elastic deformation signals associated with surface loading must be removed from the observations. In this study, we focus on seasonal variation in vertical and horizontal motions of southern Tibet by performing a joint analysis of GRACE (Gravity Recovery and Climate Experiment) and GPS data, not only using continuous sites but also GPS campaign-mode sites. We found that the GPS-observed and GRACE-modeled seasonal oscillations are in good agreements, and a seasonal displacement model demonstrates that the main reason for seasonal variations in southern Tibet is from the summer monsoon and its precipitation. The biggest loading appears from July to August in the summer season. Vertical deformations observed by GPS and modeled by GRACE are two to three times larger than horizontal oscillations, and the north components demonstrate larger amplitudes than the east components. We corrected the GPS position time series using the GRACE-modeled seasonal variations, which gives significant reductions in the misfit and weighted root-mean-squares (WRMS). Misfit (χ2 divided by degree of freedom) reductions for campaign sites range between 20% and 56% for the vertical component, and are much smaller for the horizontal components. Moreover, time series of continuous GPS (cGPS) sites near the 2015 Nepal earthquakes must be corrected using appropriate models of seasonal loading for analyzing postseismic deformation to avoid biasing estimates of the postseismic relaxation.

The quasi-biennial vertical oscillations at global GPS stations: identification by ensemble empirical mode decomposition

Modeling nonlinear vertical components of a GPS time series is critical to separating sources contributing to mass displacements. Improved vertical precision in GPS positioning at stations for velocity fields is key to resolving the mechanism of certain geophysical phenomena. In this paper, we use ensemble empirical mode decomposition (EEMD) to analyze the daily GPS time series at 89 continuous GPS stations, spanning from 2002 to 2013. EEMD decomposes a GPS time series into different intrinsic mode functions (IMFs), which are used to identify different kinds of signals and secular terms. Our study suggests that the GPS records contain not only the well-known signals (such as semi-annual and annual signals) but also the seldom-noted quasi-biennial oscillations (QBS). The quasi-biennial signals are explained by modeled loadings of atmosphere, non-tidal and hydrology that deform the surface around the GPS stations. In addition, the loadings derived from GRACE gravity changes are also consistent with the quasi-biennial deformations derived from the GPS observations. By removing the modeled components, the weighted root-mean-square (WRMS) variation of the GPS time series is reduced by 7.1% to 42.3%, and especially, after removing the seasonal and QBO signals, the average improvement percentages for seasonal and QBO signals are 25.6% and 7.5%, respectively, suggesting that it is significant to consider the QBS signals in the GPS records to improve the observed vertical deformations.

Automatic Three-Dimensional Measurement of Large-Scale Structure Based on Vision Metrology

All relevant key techniques involved in photogrammetric vision metrology for fully automatic 3D measurement of large-scale structure are studied. A new kind of coded target consisting of circular retroreflective discs is designed, and corresponding detection and recognition algorithms based on blob detection and clustering are presented. Then a three-stage strategy starting with view clustering is proposed to achieve automatic network orientation. As for matching of noncoded targets, the concept of matching path is proposed, and matches for each noncoded target are found by determination of the optimal matching path, based on a novel voting strategy, among all possible ones. Experiments on a fixed keel of airship have been conducted to verify the effectiveness and measuring accuracy of the proposed methods.

Identification of DNA basepairing via tunnel-current decay

We propose a new approach for reading the sequence of a DNA molecule passing between electrodes on a nanopore, using hydrogen bond-mediated tunneling signals. The base-electrode interaction is modeled using a nucleo-base-functionalized STM probe that is pulled away from a nucleoside monolayer. Watson-Crick recognition results in slow decay of the tunnel current, which is uniquely characteristic of the basepair in over half the reads. Thirteen independent reads would yield the desired 99.99% accuracy.