Explore the vast range of titles available.

The western segment of the suture zone between the Yangtze and Cathaysia blocks, which is the most important tectonic boundary related to the formation and rifting of south China, is enigmatic and not fully understood due to the sporadic exposure of Precambrian strata and ophiolites. Three‐dimensional electrical resistivity models derived from inversion of magnetotelluric data identified a lithospheric‐scale conductive zone extending northeastwards beneath the Youjiang basin, which was interpreted as the western segment of the suture zone. The high conductivity and coincident high magnetic anomalies closely match the location of Carlin‐type gold deposits, which can be explained by fluids and gold‐bearing sulfide minerals in a fossil suture zone. Inconsistent with the southeast‐dip resolved at the eastern segment of the suture zone (the Jiangshan‐Shaoxing fault) in the earlier study, the slightly north‐dipping geometry at the western suture zone implies the reactivation by northward subduction and closure of the Paleo‐Tethys Ocean. Plain Language Summary Suture zones are tectonic scars that mark the subduction, collision, and amalgamation of different tectonic units. Imaging suture zones is critical because they are first‐order tectonic boundaries and can play a guiding role in continental rifting and breakup. The suture zone between the Yangtze and Cathaysia blocks is the most important tectonic boundary in south China. While the eastern segment of this suture zone is agreed to be located along the Jiangshan‐Shaoxing fault, the western segment of this suture zone is under debate. This study used magnetotelluric method to image the deep earth structure beneath the western boundary of the Yangtze and Cathaysia blocks. We identified a high‐conductivity zone from the subsurface to a depth of 80 km and interpreted this zone as being caused by fluids and gold‐bearing sulfide minerals in a fossil suture zone. The conductivity pattern observed in this suture zone differs from that in the earlier study about the eastern segment of the suture zone. This implies that the western suture zone has been modified after the initial collision of the two continents, perhaps by northward subduction and closure of the Paleo‐Tethys Ocean. Key Points Magnetotelluric array data were inverted to obtain the lithospheric resistivity structure in southwestern China Three‐dimensional resistivity model images a Precambrian suture beneath the Youjiang basin The geometry resolved at the suture zone implies the reactivation by northward subduction and closure of the Paleo‐Tethys Ocean

Regional relative sea level changes are most relevant for coastal communities and remain challenging to understand. China’s adjacent seas are among the world’s most vulnerable regions to sea level rise. This paper investigates the sea level budgets in China’s adjacent seas over the past 20 years. We use multiple time-varying gravity field data and steric data to assess the uncertainties of some components in the sea level budget and the contributions of mass loss from ice sheets, glaciers, and terrestrial water storage changes to regional relative sea level changes were estimated using sea level fingerprints. The sea level budget results based on ensemble mean data show that the root mean square errors of the budget residuals in the Bohai Sea, Yellow Sea, East China Sea, South China Sea, and Northwest Pacific are 40±3 mm, 52±4 mm, 36±2 mm, 23±2 mm, and 11±1 mm, respectively. A single dataset fails to close the long-term sea level trends for all regions within a 65% confidence interval. We discussed the impacts of each component on the budget residuals and identified steric data and the ocean dynamics model as the main reasons for the excessive residuals. The de-aliasing product of the GRACE satellite, AOD1B model, is primarily responsible for the strong interannual signals in the residuals of the sea level budget in the Bohai Sea, Yellow Sea, and East China Sea.

Interplay between seismic and aseismic slip could shed light on the frictional properties and seismic potential of faults. The well-recorded 2023 Kahramanmaraş earthquake doublet provides an excellent opportunity to understand their partitioning on strike-slip faults. Here, we utilize InSAR and strong motion data to derive the coseismic rupture during the doublet, ~4-month postseismic afterslip, and slip distributions of two Mw>6.0 aftershocks. Our results show that afterslip appears to be complementary to coseismic slip and aftershocks, accounting for ~11.3% of the coseismic moment. Aftershocks mainly fall within the regions of positive Coulomb stresses caused by afterslip and follow a temporal decay similar to that of afterslip, indicating that aftershock production is the failure of small asperities loaded by the afterslip. The early postseismic afterslip is released ~93.7% aseismically and ~6.3% seismically by aftershocks. Our modeling results thus depict a complex fault system with highly variable slip patterns and stresses. The study examines interplay between seismic and aseismic slip associated with the 2023 Kahramanmaraş earthquake doublet, constrained by InSAR and strong motion data, shedding light on the frictional properties and seismic potential of faults.

The East China Sea, situated at the intersection of the Eurasian, Philippine Sea, and Pacific plates, is characterized by complex geology influenced by tectonic phenomena such as plate movements, volcanism, faults, and uplifts. Crustal density structure inversion provides a thorough understanding of the region's geological history as well as Earth's dynamical evolution, providing critical insights into seismic disaster mitigation, resource exploration, marine environmental protection, and maritime safety. The inversion process, on the other hand, presents challenges in data quality, quantity, model complexity, uncertainty, and computational resources. With the advancement of next-generation satellite gravity measurements and developing inversion techniques, the inversion of marine crustal density structures promises to be more precise and comprehensive. We explored the density distribution in the East China Sea and surrounding areas using an innovative wavenumber domain three-dimensional density imaging method along with high-precision global satellite gravity data. By overcoming data quality and computing resource constraints, wavenumber domain three-dimensional density imaging has transformed the discipline of marine geophysics, successfully delivering accurate density distributions in the study area. We were able to get a more precise and comprehensive characterization of the crustal density structure by combining high-precision satellite gravity data and cutting-edge imaging methods. Our investigation has unveiled previously unknown details about density distribution in the East China Sea and its environs. The East China Sea shelf displays smooth low-density perturbations at 18 km depth, whereas the trench–arc–basin region exhibits increasing density perturbations. Notably, the Okinawa Trough, which is surrounded by the Tokara Volcanic Ridge and the Ryukyu Trench, displays strong positive anomalies with a north–northeastern to northeastern orientation. In contrast, the Ryukyu Ridge and the Philippine Sea Basin exhibit smaller negative values and substantial northwestward positive density trends, respectively. These findings indicate diverse material distribution, which provides important insights into the area’s geological evolution and tectonic processes. This study adds new insights into density distribution in the East China Sea and adjacent regions, offering information on the geological complexity of the region. The research lays the groundwork for future research on crustal dynamics and enhances the field of marine geophysics and related disciplines. Graphical abstract

The noise level of gravity stations is an important indicator for measuring the operating status of a station and is a prerequisite for evaluating whether the station’s observations can be used to extract weak geodynamic signals. With the continuous expansion of areas of human activity, gravity stations originally located in the wild may become increasingly closer to cities. Whether their noise levels change is an important issue that is worthy of attention. Based on power spectrum analyses and probability density function methods, the noise level of the superconducting gravimeter (SG) at Jiufeng station in Wuhan in the seismic frequency band of 0.001–0.04 Hz was calculated, and its time-varying characteristics were analyzed. The noise level of Jiufeng station did not change significantly before and after the lockdown of Wuhan due to the COVID-19 epidemic in 2020. No significant changes in the noise level were found before and after the official operation of Wuhan Metro Line 19 at the end of 2023. From October 2016 to April 2017, the noise level showed an abnormal trend of suddenly rapidly rising and then slowly declining, which was found to be caused by a tilt problem in the gravity sensor. Overall, in the seismic frequency band of 0.001–0.04 Hz, the noise level at Jiufeng station showed seasonal variation characteristics, and the noise was stronger in winter than in summer, which is consistent with the characteristics of Earth’s hum. Since January 2022, the noise level has shown an increasing trend year by year. The results of this study can provide an important reference for the operation of gravity stations and the extraction of weak geodynamic signals.

As a major grain-producing region in China, the North China Plain (NCP) faces serious challenges such as water shortage and land subsidence. In late 2014, the Central Route of the South-to-North Water Diversion Project (SNWD-C) began to provide NCP with water resources. However, the effectiveness of this supply in mitigating land subsidence remains a pivotal and yet unassessed aspect. In this paper, we utilized various geodetic datasets, including the Gravity Recovery and Climate Experiment (GRACE) and GRACE Follow On (GRACE-FO), Global Navigation Satellite System (GNSS) and leveling data, to conduct a spatial-temporal analysis of the equivalent water height (EWH) and vertical ground movement in the NCP. The results reveal a noteworthy decline in EWH from 2011 to 2015, followed by a slight increase with minor fluctuations from 2015 to 2020, demonstrating a strong correlation with the water resources supplied by the SNWD-C. The GRACE-derived surface deformation rate induced by hydrological loading is estimated to be <1 mm/yr. In comparison, GNSS-derived vertical ground movements exhibit considerable regional differences during the 2011–2020 period. Substantial surface subsidence is evident in the central and eastern NCP, contrasting with a gradual uplift in the front plain of the Taihang Mountains. Three-stage leveling results indicate that the rate of subsidence in the central and eastern plains is gradually increasing with the depression area expanding from 1960 to 2010. Based on these geodetic results, it can be inferred that the SNWD-C’s operation since 2014 has effectively mitigated the reduction in terrestrial water storage in the NCP. However, land subsidence in the NCP persists, as the subsidence rate does not turn around in sync with the change in EWH following the operation of SNWD-C. Consequently, it’s necessary to maintain and enforce existing policies, including controlling groundwater exploitation and water resources supply (e.g., SNWD-C) to curtail the exacerbation of land subsidence in the NCP. Additionally, continuous monitoring of land subsidence by GRACE, GNSS, leveling and other geodetic techniques is crucial to enable timely policy adjustments based on monitoring results.

Satellite altimetry is the main tool for constructing global or regional marine gravity fields. To improve the accuracy and spatial resolution, it is necessary to fuse multi-mission altimeters. How to determine the weights of multi-mission altimeters is a crucial issue, making the conventional calculation process very complex. In addition, traditional satellite inversion methods are often independent of shipborne gravity, which is used only as validation data, thus not take full advantages of high accuracy and resolution of shipborne gravity. In this study, we introduce a convolutional neural network (CNN) to merge the vertical deflections (DOVs) obtained from multi-altimeter missions to construct a marine gravity model in the South China Sea. High-accuracy shipborne gravity and a dataset comprising DOVs and geo-locations are employed as input data for neural network training. For the validation of CNN method, the gravity model is also computed by conventional Inverse Vening Meinesz (IVM) method. Independent shipborne gravity measurements and SIO V32.1, DTU17 models are used as validation data. The evaluation results show that the CNN-derived model achieves a higher level of accuracy, yielding a standard deviation (STD) of 3.21 mGal, with an improvement of 36.56% compared to IVM-derived model. More than 92% of the differences between the CNN-derived model and shipborne gravity are less than 5 mGal. In addition, spectral analysis results further show that the CNN-derived model has stronger energy at short wavelengths (less than 25 km) compared to other models. These findings reveal that CNN method is feasible for marine gravity recovery and the CNN-derived model can achieve higher accuracy. The CNN method can improve the accuracy and spectral characteristics of the constructed gravity model by taking advantage of the high accuracy and high resolution of shipborne gravity. Graphical Abstract

Advancements in satellite altimetry have significantly enhanced high-resolution mean sea surface (MSS) models, enabling the computation of high-resolution vertical gravity anomaly gradient (VGAG) models. This study focused on the methodology for computing VGAG models using MSS models, introducing innovative improvements to established techniques. Using the SDUST2020 MSS model within the Arabian Sea research area, the DTU22 and CNES-CLS22 mean dynamic topography (MDT) models, and the XGM2019e_2159 Earth gravity field model for the remove–restore process, the short-wavelength geoid was derived. To harness the extensive marine gravity field information within the MSS model, the study considered the complex marine environment and calculated the second-order derivatives of the geoid in multiple directions. These derivatives were then used to determine their north–south and east–west components through the least squares method, resulting in the computation of the short-wavelength VGAG. By restoring the long-wavelength VGAG, a VGAG model for the study area was established. Finally, the results were analyzed using the SIO V32.1 VGAG model (named curv). Experimental results demonstrated that this approach effectively extracted marine gravity field information from the MSS model using multidirectional data, mitigating the amplification of geoid uncertainties caused by second-order derivatives. Graphical Abstract

The precise prediction of polar motion parameters is needed for the astrogeodynamics, navigation and positioning of the deep space probe. However, the current prediction methods are limited to predicting polar motion for specific periods, either short- or long-term. In this study, a sliding multilayer perceptron (MLP) method combined singular spectrum analysis (SSA) and autoregressive moving average (ARMA) for short- and long-term polar motion prediction was proposed. MLP was introduced into PM prediction due to its automatic learning characteristics and its ability to effectively process nonlinear and multi-dimensional data. The SSA was used to extract and predict the principal components of polar motion, while the remaining components were predicted using ARMA. In the meantime, SSA and ARMA were used to provide training data and target learning data for the MLP model. MLP input data were constructed by sliding processing with a window of 7 days, composed of n series of the same length (18 years). Finally, MLP was employed to predict the residuals generated during SSA and ARMA prediction. To evaluate the accuracy of the proposed method, the polar motion prediction was applied for a 364-day lead time based on the IERS EOP 14C04 product. The method outperformed the IERS Bulletin A, as demonstrated by the mean-absolute errors of the x and y components of polar motion on the 30th day, which were lower (5.14 mas and 3.37 mas, respectively) than those predicted by IERS Bulletin A (6.66 mas and 3.94 mas). Similarly, the mean-absolute errors on the 364th day were 17.79 mas and 16.29 mas, respectively, compared to the 19.24 mas and 18.81 mas predicted by IERS Bulletin A.