Moho depth inversion in the Tibetan Plateau from high-precision gravity data

The Tibetan Plateau(TP)is the youngest orogenic belt resulting from a continental collision on the Earth.It is a natural laboratory for studying continental dynamics,such as continental convergence,plate subduction,and plateau uplift.Investigating the deep structure of the TP has always been a popular issue in geological research.The Moho is the boundary between the crust and the mantle and therefore plays a crucial role in the Earth’s structure.Parameters such as depth and lateral variation,as well as the fine structure of the crust-mantle interface,reveal the lithospheric dynamics in the TP.Two methods are generally employed to study the Moho surface:seismic detection and gravity inversion.Seismic detection has the characteristic of high precision,but it is limited to a few cross-sectional lines and is quite costly.It is not suitable for and cannot be carried out over a large area of the TP.The Moho depth over a large area can be obtained through gravity inversion,but this method is affected by the nature of gravity data,and the accuracy of the inversion method is lower than that of seismic detection.In this work,a high-precision gravity field model was selected.The Parker-Oldenburg interface inversion method was used,within the constraints of seismic observations,and the Bott iteration method was introduced to enhance the inversion efficiency.The Moho depth in the TP was obtained with high precision,consistent with the seismic detection results.The research results showed that the shape of the Moho in the TP is complex and the variation range is large,reaching 60−80 km.In contrast with the adjacent area,a clear zone of sharp variation appears at the edge of the plateau.In the interior of the TP,the buried depth of the Moho is characterized by two depressions and two uplifts.To the south of the Yarlung Zangbo River,the Moho inclines to the north,and to the north,the Moho depresses downward,which was interpreted as the Indian plate subducting to the north below Tibet.The Moho depression on the north side of the Qiangtang block,reaching 72 km deep,may be a result of the southward subduction of the lithosphere.The Moho uplift of the Qiangtang block has the same strike as the Bangong−Nujiang suture zone,which may indicate that the area is compensated by a low-density and low-velocity mantle.

Noise reduction and periodic signal extraction for GNSS height data in the study of vertical deformation

Global navigation satellite system(GNSS)technique has irreplaceable advantages in the continuous monitoring of surface deformation.Reducing noise to improve the signal-to-noise ratio(SNR)and extract the concerned signals is of great significance.As an improved algorithm of empirical mode decomposition(EMD),complete ensemble empirical mode decomposition with adaptive noise(CEEMDAN)algorithm has better signal processing ability.Using the CEEMDAN algorithm,the height time series of 29GNSS stations in Chinese mainland were analyzed,and good denoising effects and extraction from periodic signals were achieved.The numerical results showed that the annual signal obtained with the CEEMDAN algorithm was significantly based on Lomb_Scargle spectrum analysis,and large differences in the long-term signals were found between the stations at different locations in Chinese mainland.With respect to data denoising,compared with the EMD and wavelet denoising algorithms,the CEEMDAN algorithm respectively improved the SNR by 29.35% and 36.54%,increased the correlation coefficient by 8.67% and 11.96%,and reduced root mean square error(RMSE)by 44.68% and 43.48%,indicating that the CEEMDAN algorithm had better denoising behavior than the other two algorithms.In addition,the results demonstrated that different denoising methods had little influence on estimating the annual vertical deformation velocity.The extraction of periodic signals showed that more components were retained by using the CEEMDAN algorithm than the EMD algorithm,which indicated that the CEEMDAN algorithm had advantages over frequency aliasing.In conclusion,the CEEMDAN algorithm was recommended for processing the GNSS height time series to analyze the vertical deformation due to its excellent features of denoising and the extraction of periodic signals.

A review of the 19th International Symposium on geodynamics and earth tide, Wuhan 2021

Modern geodetic technologies such as high-precision ground gravity measurements,satellite gravity measurements,the global navigation satellite system,remote sensing methods,etc.provide rich observation data for monitoring various geodynamic processes of the global Earth and its surface.The 19th International Symposium on Geodynamics and Earth Tides brought together scientific researchers from 26 countries around the world,shared the application of various measurements in different geoscience issues,covering Earth tidal deformation,oceanic and atmospheric loading effects,earthquake cycle,hydrology,Earth rotation changes,etc.,and provided a precious exchange platform for global peers.

A review of tidal triggering of global earthquakes

Earthquake prediction remains a challenging and difficult task for scientists all over the world.The tidal triggering of earthquakes is being proven by an increasing number of investigations,most of which have shown that earthquakes are positively correlated with tides,and thus,tides provide a potential tool for earthquake prediction,especially for imminent earthquakes.In this study,publications concerning the tidal triggering of earthquakes were compiled and analyzed with regard to global earthquakes,which were classified into three main types:tectonic,volcanic,and slow earthquakes.The results reveal a high correlation between tectonic earthquakes and tides(mainly for semidiurnal and diurnal tides;14-day tides) before and after the occurrence of significant earthquakes.For volcanic earthquakes,observations of volcanoes on the seafloor and land indicate that volcanic earthquakes in near-shore volcanic areas and mid-ocean ridges have a strong correlation with tidal forces,mostly those with semidiurnal and diurnal periods.For slow earthquakes,the periodicity of the tremor duration is highly correlated with semidiurnal and diurnal tides.In conclusion,the tidal triggering of these three types of earthquakes makes a positive contribution to earthquake preparation and understanding the triggering mechanism,and thus,the prediction of these types of earthquakes should be investigated.However,there are still several inadequacies on this topic that need to be resolved to gain a definitiveanswer regarding the tidal triggering of all earthquakes.The main inadequacies are discussed in this paper from our point of view.

Editorial note for the geodesy and geodynamics journal special issue contemporary research in geodynamics and earth tides-Selection from the 19th international symposium on geodynamics and earth tides,2021,Wuhan,China

This volume aims at providing a platform for sharing valuable topics discussed at the 19th International Symposium on Geodynamics and Earth Tides,23-26 June 2021,Wuhan,China.The complete overview of all nineteen Symposia is found in Table 1,where the times and venues are listed.

Research Progress in Surface and Marine Gravimetry

Gravity field is the most basic physical field generated by the material properties of the Earth system.It reflects the spatial distribution,movement and change of materials determined by the interaction and dynamic process inside the Earth.Over the years,a variety of technical means have been used to detect the Earth’s gravity field and supported numerous studies on the global change,resource detection,geological structure movement,water resources change and other related fields of research.Here is part of the progress in surface and marine gravimetry obtained by Chinese geodesy scientists from 2019 to 2023 from the following aspects,including:①Continuous gravity network in Chinese mainland;②Application of superconducting gravity measurement;③Network adjustment for continental-scale gravity survey campaign and data quality control;④Regional time-variable gravity field and its application;⑤Research progress on novel technologies for gravity inversion;⑥Research progress on marine gravity field determination;⑦Application research on marine gravity field.

Geodesy Discipline: Progress and Perspective

The geodesy discipline has been evolving and constantly intersecting and merging with other disciplines in the last 50 years,due to the continuous progress of geodetic observation techniques and expansion of application fields.This paper first introduces the development and roles of geodesy and its formation.Secondly,the development status of geodesy discipline is analyzed from the progress of observation techniques and cross-discipline formation is analyzed from the expansion of application fields.Furthermore,the development trend of geodesy is stated from the perspective of national requirements and scientific developments.Finally,the sub-disciplines for geodesy are suggested at the present stage,based on the requirements of the National Natural Science Foundation of China and development status of geodesy itself,which can provide references for topic selection and fund application of geodetic scientific research.

2004年12月26日苏门答腊——安达曼地震激发的地球自由振荡

根据300-500s的面波估计得到的哈佛大学矩心矩张量断层机制，当周期大于1000s时，地球的地震自由振荡的振幅将异常地大。我们对一个比较陡的断层应用根据地震的体波和面波得到的更真实的破裂模型，估计了自由振荡的振幅近似等于地震矩（6．5×10^22 N·m），这相当于矩震级为9．15。由于破裂的持续时间达600s，故断层破裂模型足以描述地震的观测结果，但却低估了大地测量位移，这些位移表明在尼科巴群岛和安达曼群岛下面发生了缓慢的断层运动。

Preface to the special issue on Lushan earthquake

On April 20th,2013,a strong earthquake(Ms7,China Earthquake Network Center)struck Lushan county of Sichuan province and the quake(hereafter referred to as Lushan earthquake)caused substantial loss of life and damage to infrastructure.Just as the 2008 Ms8 Wenchuan earthquake,the Lushan earthquake also occurred on the Longmenshan fault system.After the Lushan earthquake,preliminary studies of the Lushan earthquake and its pos-

Co-seismic change in ocean bottom topography:Implication to absolute global mean sea level change

Earthquakes perturb both the ocean bottom topography due to displacements of sea floor and the geoid due to mass redistribution, which induces the relative sea level(RSL) change. However, the relative global mean sea level(GMSL) change is zero in that sea water mass is conserved. But the absolute GMSL change is not zero because earthquakes displace total ocean mass with respect to the Earth’s center of mass(CM) which remains unchanged after an earthquake. This displacement, i.e. the absolute GMSL change, may be detectable by altimetry since the satellites are orbiting around CM. In this paper, we proposed a method to estimate co-seismic absolute GMSL change caused by earthquakes based on the point dislocation theory for a spherically symmetric, non-rotating, elastic and isotropic(SNREI) Earth.This change can be directly connected to the perturbation of ocean bottom topography. We first computed co-seismic displacements as well as the change in geo-potential and solved the sea level equation to validate the insignificance of the oceans’ feedback, i.e. the loading effect due to RSL change, to co-seismic displacements. The results imply that the loading effect due to RSL change is negligible on displacements while is considerable on geoid. We then computed the absolute GMSL change caused by co-seismic vertical and horizontal displacements by making use of the integrated Green’s function method. The numerical results show that a large earthquake may raise the absolute GMSL by magnitude of sub-millimeter and the recent three large events cause GMSL to rise about one millimeter, in which the contribution from horizontal displacement is non-negligible.

Theoretical calculation of tidal Love numbers of the Moon with a new spectral element method

The tidal Love numbers of the Moon are a set of nondimensional parameters that describe the deformation responses of the Moon to the tidal forces of external celestial bodies.They play an important role in the theoretical calculation of the Moon’s tidal deformation and the inversion of its internal structure.In this study,we introduce the basic theory for the theoretical calculation of the tidal Love numbers and propose a new method of solving the tidal Love numbers:the spectral element method.Moreover,we explain the mathematical theory and advantages of this method.On the basis of this new method,using 10 published lunar internal structure reference models,the lunar surface and lunar internal tidal Love numbers were calculated,and the influence of different lunar models on the calculated Love numbers was analyzed.Results of the calculation showed that the difference in the second-degree lunar surface Love numbers among different lunar models was within 8.5%,the influence on the maximum vertical displacement on the lunar surface could reach±8.5 mm,and the influence on the maximum gravity change could reach±6μGal.Regarding the influence on the Love numbers inside the Moon,different lunar models had a greater impact on the Love numbers h_(2) and l_(2) than on k_(2) in the lower lunar mantle and core.

精密(量子)测量时代下时空基准研究中的关键科学问题和核心技术

基于第316期双清论坛“精密(量子)测量时代下时空基准研究”会议内容,本文总结了我国在精密(量子)测量时代下时空基准研究所面临的国家重大需求,回顾了时空基准在国际国内的发展历程以及我们面临的挑战,凝练了该领域未来5~10年的重大关键科学问题和核心技术,建议在时空基准理论与方法、关键技术和设备研发、时空基准应用等相关领域重点攻关,并就科学基金重点资助项目给出了建议。

Observation and research of deep underground multi-physical fields—Huainan–848 m deep experiment

Compared with the surface,the deep environment has the advantages of allowing“super-quiet and ultra-clean”-geophysical field observation with low vibration noise and little electromagnetic interference,which are conducive to the realization of long-term and high-precision observation of multi-physical fields,thus enabling the solution of a series of geoscience problems.In the Panyidong Coal Mine,where there are extensive underground tunnels at the depth of 848 m below sea level,we carried out the first deep-underground geophysical observations,including radioactivity,gravity,magnetic,magnetotelluric,background vibration and six-component seismic observations.We concluded from these measurements that(1)the background of deep subsurface gravity noise in the long-period frequency band less than 2 Hz is nearly two orders of magnitude weaker than that in the surface observation environment;(2)the underground electric field is obviously weaker than the surface electric field,and the relatively high frequency of the underground field,greater than 1 Hz,is more than two orders of magnitude weaker than that of the surface electric field;the east-west magnetic field underground is approximately the same as that at the surface;the relatively high-frequency north-south magnetic field underground,below 10 Hz,is at least one order of magnitude lower than that at the surface,showing that the underground has a clean electromagnetic environment;(3)in addition to the highfrequency and single-frequency noises introduced by underground human activities,the deep underground space has a significantly lower background vibration noise than the surface,which is very beneficial to the detection of weak earthquake and gravity signals;and(4)the underground roadway support system built with ferromagnetic material interferes the geomagnetic field.We also found that for deep observation in the“ultra-quiet and ultra-clean”environment,the existing geophysical equipment and observation technology have problems of poor adaptability and insufficient precision as well as data cleaning problems,such as the effective separation of the signal and noise of deep observation data.It is also urgent to interpret and comprehensively utilize these high-precision multi-physics observation data.

Complexities of the Turkey-Syria doublet earthquake sequence

WHAT,WHERE,AND WHEN In the early morning of February 6th,2023,an M7.8 earthquake occurred in southeastern T€urkiye near the northern border of Syria.The event initiated a complex sequence of aftershocks,including an M7.6 earthquake about 9 h later and 90 km to the north(Figures 1A and 1B).The earthquake sequence is also referred to as a strong doublet earthquake sequence.Aftershocks of the two strong earthquakes occurred along two separate branches of the East Anatolia Fault,with lengths of up to 300 km,and some aftershocks occurred in Syria(NEIC/USGS,2023).

Probing signals of atmospheric gravity waves excited by the July 29,2021 M_(W8.2) Alaska earthquake

It is commonly believed that the atmosphere is decoupled from the solid Earth.Thus,it is difficult for the seismic wave energy inside the Earth to propagate into the atmosphere,and atmospheric pressure wave signals excited by earthquakes are unlikely to exist in atmospheric observations.An increasing number of studies have shown that earthquakes,volcanoes,and tsunamis can perturb the Earth's atmosphere due to various coupling effects.However,the observations mainly focus on acoustic waves with periods of less than 10 min and inertial gravity waves with periods of greater than 1 h.There are almost no clear observations of gravity waves that coincide with observations of low-frequency signals of the Earth's free oscillation frequency band within 1 h.This paper investigates atmospheric gravity wave signals within1 h of surface-atmosphere observations using the periodogram method based on seismometer and microbarometer observations from the global seismic network before and after the July 29,2021 Mw8.2Alaska earthquake in the United States.The numerical results show that the atmospheric gravity wave signals with frequencies similar to those of the Earth's free oscillations_0S_(2) and_0T_(2) can be detected in the microbaro meter observations.The results con firm the existence of atmospheric gravity waves,indicating that the atmosphere and the solid Earth are not decoupled within this frequency band and that seismic wave energy excited by earthquakes can propagate from the interior of the Earth to the atmosphere and enhance the atmospheric gravity wave signals within 1 h.

Check of free oscillation signal with SG data

The observation of the Earth’s free oscillation is an important way to study the Earth’s inner structure and geophysical parameters. Many geophysicists have investigated the Earth’s normal modes with strainmeters and spring gravimeters. Since the superconducting gravimeters appeared, people have gained a new kind of more stable instruments to detect the normal modes. The measurement for the Earth’s free oscillation with the superconducting gravimeters is one of the main goals of the Global Geodynamics Project. In this note, we have investigated the Earth’s normal modes excited by Peru Ms 7.9 Earthquake on June 23, 2001 with the superconducting gravimeter C032 at the Wuhan station. After having completed the subtraction of gravity tides with the 20th order polynomial fitting, and the correction for pressure and the analysis of the noise spectrum of the superconducting gravimeter, we accurately observed all base normal modes from 0S0 to 0S32 and the splittings of 0S2 and 0S3. Those results show a good