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.

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.

The deforming and rotating Earth——A review of the 18th International Symposium on Geodynamics and Earth Tide,Trieste 2016

The 18th International Symposium on Geodynamics and Earth Tides 2016 covered phenomena that generate temporal variations in geodetic and geophysical observations. In calculating the stress field for Earth tides, the observed geodetic response is used for defining the Earth＇s theology, the Earth internal structure, ＇Earth rotation parameters, and the functioning of the sophisticated instrumentation mounted on Earth and satellites. The instrumentation capable of observing Earth tides, measures changes generated by lithospheric plate movements, as the earthquake cycle and volcanism. Hydrology, tem- perature, and pressure, either of natural or anthropogenic origin, affect the high precision observations, and therefore must be included in this study-realm.

Earth rotation deceleration/acceleration due to semidiurnal oceanic/atmospheric tides:Revisited with new calculation

The global oceanic/atmospheric tides exert decelerating/accelerating secular torques on the Earth rotation. We developed new formulations to accurately calculate amounts of two kinds of secular tidal torques. After Melchior, we found that an additional factor 1+k-l = 1.216, which has been formerly neglected, must be multiplied unto the tidal torque integral. By using our refined formulations and the recent oceanic/atmospheric global tide models, we found that:(i) semidiurnal oceanic lunar/solar tides exert decelerating torques of about-4.462 × 10^(16)/-0.676 × 10^(16) Nm respectively and(ii) atmospheric S_2 tide exerts accelerating torque of 1.55 × 10^(15) Nm. Former estimates of the atmospheric S_2 tidal torque were twice as large as our estimate due to improper consideration of loading effect. We took the load Love number for atmospheric loading effect from Guo et al.(2004). For atmospheric loading of spherical harmonic degree two, the value of k′=-0.6031 is different from that for ocean loading as k′ =-0.3052,while the latter is currently used for both cases-ocean/atmospheric loading-without distinction. We discuss(i) the amount of solid Earth tidal dissipation(which has been left most uncertain) and(ii) secular changes of the dynamical state of the Earth-Moon-Sun system. Our estimate of the solid Earth tidal torque is-4.94×10^(15) Nm.

Zero drift and solid Earth tide extracted from relative gravimetric data with principal component analysis

Zero drift and solid Earth tide corrections to static relative gravirnetric data cannot be ignored. In this paper, a new principal component analysis （PCA） algorithm is presented to extract the zero drift and the solid Earth tide, as signals, from static relative gravimetric data assuming that the components contained in the relative gravimetric data are uncorrelated. Static relative gravity observations from Aug. 15 to Aug. 23, 2014 are used as statistical variables to separate the signal and noise with PCA to obtain desired signals. The results of the linear drift extracted by PCA are consistent with those calculated by the least squares linear fitting, and the differences only reach to 10-2μGal/day order of magnitude. Furthermore, PCA is used to estimate the solid Earth tide from the relative gravimetric data corrected by the zero drift. The statistical results are consistent with the results derived from the solid Earth tide correction provided by the internal software of the CG-5 gravimeter （SCINTREX Limited Ontario Canada）. The statistical results of the differences between the two methods are both less than 8 ,Gal, and the RMSs for 9 days are all less than 5 μGal.

Effect of solid-earth-tide on GPS time series

The solid-earth-tide models IERS1992 and IERS2003 are used to analyze some GPS-baseline,vertical-component and zenith-tropospheric-delay data from the Crustal Movement GPS Continuous Observation Net-work of Shandong and IGS stations. The results show that the differences between the baselines computed with the different models are at sub-millimeter level, and the differences in vertical component is direct proportional to station latitude. Also the amplitude of ZTD differences is about 0. 6-1.0 mm, which is 6% -8% of the amplitude of solid-earth-tide differences. Although these effects are quite small, to analyze non-tidal deformation correctly, we should still use a single standard for processing GPS data.

On the Imminent Regional Seismic Activity Forecasting Using INTERMAGNET and Sun-Moon Tide Code Data

In this paper we present an approach for forecasting the imminent regional seismic activity by using geomagnetic data and Earth tide data. The time periods of seismic activity are the time periods around the Sun-Moon extreme of the diurnal average value of the tide vector module. For analyzing the geomagnetic data behaviour we use diurnal standard deviation of geomagnetic vector components F (North, East, Down) for calculating the time variance GeomagSignal. The Sun storm influence is avoided by using data for daily A-indexes (published by NOAA). The precursor signal for forecasting the incoming regional seismic activity is a simple function of the present and previous day GeomagSignal and A-indexes values. The reliability of the geomagnetic “when, regional” precursor is demonstrated by using statistical analysis of day difference between the times of “predicted” and occurred earthquakes. The base of the analysis is a natural hypothesis that the “predicted” earthquake is the one whose surface energy density in the monitoring point is bigger than the energy densities of all occurred earthquakes in the same period and region. The reliability of the approach was tested using the INTERMAGNET stations data located in Bulgaria, Panagurishte, PAG (Jan 1, 2008-Jan 29, 2014), Romania, Surlari, SUA (Jan 1, 2008-Jan 27, 2014), Italy, L’Aquila, AQU (Jan 1, 2008-May 30, 2013) in the time of EU IRSES BlackSeaHazNet (2011-2014) project. The steps of program for solving the “when, where and how” earthquake prediction problem are shortly described.

Inversion of the Strain Accumulation State in Crustal Interior from Observations of Earth Tide Tilts

The strain accumulation state in crustal interior and its migration trend can be inverted from observational data of earth tide tilts.In China,large earthquakes frequently occur in the border area between Yunnan and Sichuan Provinces.Using the observational data for years from 7 stations in this area,γ(the amplitude factor of earth tide tilts)values that characterize the strain accumulation state in crustal interior have been inverted in this paper.It has been found out that,Ms】5 events in the area often occur when γreaches an anomaly value.Using the time-space collocational stereomodelling the migration trend of strain accumulation can be determined to make a prediction for the probable location of a forthcoming earthquake in the area.

Short-term and imminent anomalies of earthquake of load and unload response ratio of the well level to earth tides

In this paper, through the nonlinear response of rock strain and stress, we have analized the physical mechanism of loading and unloading response ratio of the well level to the earth tides,the respouse of an aquifer of confined well to bulk strain tide and showed two methods of the calculation of loading and unloading response ratio of the well level to the earth tides. We took the example of the Yu 01 well, which is near the epicenter of Heze M S 5.9 earthquake, calculated the response rate and loading and unloading response ratio of two kinds of the earth tides of it. The response rate and response ratio before the earthquake had the variation of increase.

Application Study on Correction Method for Lag of Water Level Response to Earth Tide and Atmospheric Pressure

The water level in a deep well instantly responds to the earth’s tide and atmospheric pressure, and varies accordingly, not only in terms of amplitude but also in the phase lag. Therefore, phase lag correction is used in analyzing digital groundwater observation data in eastern China. Calculation results presented by the authors in this paper show that the correction method is effective in the identification of anomalous changes for short-term seismic precursors. The correction method can also be applied to the processing of observed deformation and tilt data.

潮汐效应对沿海多条带时序InSAR地表形变监测的影响

潮汐效应对沿海大范围、高精度的合成孔径雷达干涉测量(Interferometric Synthetic Aperture Radar,InSAR)技术地表形变监测的影响不可忽略.潮汐效应中的固体潮(Solid Earth Tide,SET)位移和海洋潮汐负荷(Ocean Tide Loading,OTL)位移会在时间和空间域上对沿海大范围InSAR地表形变监测产生分米级的误差,而且不同条带中地表潮汐形变时空差异性会对多条带InSAR影像拼接产生较大误差.本文针对沿海地区多条带InSAR形变时序,重点分析了潮汐效应在不同条带InSAR地表沉降监测中的影响,并采用多种潮汐位移估计方法对多条带地表潮汐形变进行改正.结果表明,研究区域不同条带时序InSAR中地表潮汐形变具有较大的时空差异性,其造成沿海地区相邻条带时序InSAR地表沉降速率差异能达到1~2 cm·a^(-1);通过固体潮模型、海潮负荷模型或GPS参考站网海潮位移改正,能够消除不同条带中地表潮汐形变时空差异性对时序InSAR形变结果拼接产生的空间高阶非线性误差,弥补了传统拼接方法在拟合潮汐位移偏差的不足.在大气延迟误差改正的基础上,时序InSAR形变残差的标准差(Standard Deviation,STD)由潮汐改正前的2.3 cm减少至0.75 cm;本文研究揭示了沿海地表潮汐形变改正对多条带时序InSAR形变结果拼接的重要性,可大幅提高沿海地区广域时序InSAR地表形变监测的精度.

利用GPS观测中国大陆固体地球潮汐形变

利用中国陆态网242个GPS连续观测站超过10年的观测数据,采用精密单点定位静态方法获得了八个主要潮波的三维潮汐位移.与建模值对比结果表明,GPS观测潮汐位移的精度可达亚毫米级.进一步将GPS观测值与OTL(Ocean Tide Loading)建模值对比,除DTU10外其他七种海潮模型与GPS观测值之间一致性较好.综合考虑验潮站观测结果,最适合中国大陆M2、O1潮波OTL建模的海潮模型是TPXO9-Atlas,地球模型为PREM.本研究还发现,M2、O1潮波GPS观测值与建模值之间的残差存在明显的大尺度空间相关性,表明中国地区存在固体潮建模误差等系统误差.本文对198个离海岸线距离超过200 km的内陆站点的系统误差进行分析,改正系统误差后内陆站点东、北和垂直方向加权均方根误差(Weighted Root Mean Square,简称WRMS)分别为:0.04 mm、0.07 mm和0.11 mm(M2),0.05 mm、0.07 mm和0.11 mm(O1);44个沿海站点WRMS分别为:0.18 mm、0.27 mm和0.59 mm(M2),0.10 mm、0.16 mm和0.30 mm(O1).部分沿海站点仍存在较大残差,主要源于海潮负荷位移建模误差,除海潮模型外还包括海水密度、地球模型、上地幔和地壳模型等因素影响.本文首先使用最新WOA(World Ocean Atlas)18海洋数据集计算海水密度和可压缩性的影响,结果表明中国内地OTL建模对海水密度变化和可压缩性并不敏感.其次,软流圈滞弹性散耗效应可使半日潮软流圈的剪切模量降低约8.5%,周日潮软流圈的剪切模量降低约9.1%,考虑滞弹性散耗效应,M2、O1潮波垂直方向上残差最大可分别减小0.50 mm、0.25 mm.最后,本文利用GPS实测数据对四种三维地球模型进行了适用性评价,结果表明除USTC地球模型受沉积层影响外,其他三种地球模型都能有效提高GPS观测结果与理论建模结果的一致性;另外全球地壳模型CRUST1.0与GPS观测值之间的残差最小,说明中国大陆海潮负荷建模精度受地壳结构精度影响较大.本文研究结果证明现今中国大陆GPS观测潮汐位移的结果已经具有足够的精度用于约束该区域地球内部结构.

固体潮模型在GPS/BDS PPP中的影响分析

针对固体潮改正模型对GPS/BDS精密单点定位(PPP)的影响问题,借助IERS2003固体潮改正模型分别对GPS、BDS-2、BDS-3、BDS-2+BDS-3、GPS+BDS-2+BDS-3定位结果进行定量分析。结果表明,固体潮改正模型使GPS、BDS-2、BDS-3、BDS-2+BDS-3、GPS+BDS-2+BDS-3的PPP结果在N、E和U方向的平均改正量分别为1.06 cm、1.01 cm、3.77 cm,0.87 cm、1.09 cm、3.48 cm,1.30 cm、1.32 cm、5.29 cm,1.45 cm、1.85 cm、6.37 cm以及1.71 cm、1.83 cm、6.41 cm。

Apparent permeability variation of underground water aquifer induced by an earthquake:A case of the Zhouzhi well and the 2008 Wenchuan earthquake

t Taking the M2 wave as calibration signals, we extract the phase shifts of the water level relative to the Earth tide in the Zhouzhi well by utilizing the cross-correlation function. And we further obtain the apparent permeability variation in the aquifer of the Zhouzhi well in 2008. Comparison with the commonly used tidal analysis software Baytap-G shows that phase shifts obtained by cross-correlation function are more stable. The resulting apparent permeability of the Zhouzhi well aquifer fluctuates with time, indicating it is a dynamically controlled parameter. The 2008 Wenchuan earthquake caused the apparent permeability increasing drastically, which is interpreted as the combination effects of effective stress changes and the barriers removing in the flow channel due to seismic wave pressure pulse. After the Wenchuan earthquake, the effective stress began to recover and the impurities deposited gradually, causing the apparent permeability to decrease a month later and almost recover to the pre-earthquake level in six months.

Atmospheric acceleration and Earth-expansion deceleration of the Earth rotation

Previous studies suggest that tidal friction gives rise to the secular deceleration of the Earth rotation by a quantity of about 2.25 ms/cy. Here we just consider additional contributions to the secular Earth rotation deceleration. Atmospheric solar semi-diurnal tide has a small amplitude and certain amount of phase lead. This periodic global air-mass excess distribution exerts a quasi-constant torque to accelerate the Earth＇s spin rotation. Using an updated atmospheric tide model, we re-estimate the amounts of this atmospheric acceleration torque and corresponding energy input, of which the associated change rate in LOD（length of day） is-0.1 ms/cy. In another aspect, evidences from space-geodesy and sea level rise observations suggest that Earth expands at a rate of 0.35 mm/yr in recent decades, which gives rise to the increase of LOD at rate of 1.0 ms/cy. Hence, if the previous estimate due to the tidal friction is correct, the secular Earth rotation deceleration due to tidal friction and Earth expansion should be 3.15 ms/cy.

The Akp-Btk value method and the resultsfor the retrieval of the parameters of the Earth's free core nutation

The rotation of the Earth's liquid core creates the Nearly Diurnal Free Wobble (NDFW). It is one of the problems of researching the Earth's liquid core for us to retrieve the parameters of the Earth's Free Core Nutation (FCN),from the Earth's liquid core resonance of gravity tide waves on the diurnal frequency band. Since 1987, some scientists in many different countries have come to calculate the parameters of FCN by using the observational data of gravity tide waves on the diumal frequency band. They basically followed the Stacking method, which needed five diurnal waves for the resolution. In this paper, authors introduced the Aky-Btk value method with clearly geometrical meaning as a new method, which only requires three very high signal-noise-ratio waves o1, K1 and P1 for the resolution. Authors chose the observational data of the three superconducting gravimeter stations respectively located in Cantley of Canada, Wuhan of China and Brussels of Belgium, to compute the parameters of FCN. It was the first time that the observational results of the parameters of FCN obtained from gravitational tide were in accord with the parameters of FCN gained from VLBI since 1987.

Abnormal Phenomena of Volume Strain before Large Earthquakes

This paper studies the imminent anomalies observed by the Sacks volume strainmeter in Erzhangying station and Tiantanghe station before 80 earthquakes with M S≥7. 0 which took place from January 2011 to April 2014 all over the world. Then,preconditions for anomaly identification are put forward for complex earthquake cases. Statistical results show that volume strain observation has a better earthquake reflecting ability for earthquakes with magnitudes larger than 7. 0 and epicentral distance within 8000 km. In addition,these results also reflect that the volume strain observation can better reflect precursory anomalies of such earthquakes. Based on categorization and description of those anomalies,we divide the anomalies into three types,that is,earth tide distortion type,abrupt change type and slow earthquake type. Furthermore,the paper makes a statistical analysis of these types and preliminarily discusses their mechanical properties as well. According to research,volume strain anomaly has an indicative significance to future strong earthquakes in the world.

Time and space analysis of two earthquakes in the Appennines (Italy)

In this paper, we study two earthquakes: the April 6th 2009 earthquake of L’Aquila in the region of Abruzzo (Italy) and the 1997 Colfiorito earthquake in the regions of Umbria and Marche (Italy). The data sets of these two earthquakes were analysed in both time and space domains. For time domain we used statistical methods and models both parametric and non-parametric. Concerning the space domain, we used Mathematical Morphology filters. The time domain analysis provides evidence of a possible correlation between seismic activities and the tides of the crust of the Earth. The results obtained show evidence that the daily number of earthquakes of the sequences proceeding and following the April 6th 2009 earthquake of L’Aquila and that of the sequence following the 1997 Colfiorito earthquake have a periodic component of occurrence with period of about 7 days. It seems that the maxima of this component occur at a position of the Moon with respect to the Earth and the Sun corresponding to approximately 3 days before the four main Moon phases. The space domain analysis indicates that the foreshock activity in both earthquakes is clustered and concentrated. Furthermore, in each of the two earthquakes the clusters are located at about 3 kilometers from the epicentre of the main shock.

An evidence for earthquake ocurence time related with geologic structure──Relations between local mean lunar times of earthquake ocurence and their fault trends

The relation between the local mean lunar time τ of earthquake occurrence and their fault trends is studied in this paper. The local mean lunar times τ of 53 earthquakes in 24 groups are calculated. Because the tidal generation force arisen by the moon is a cyclic function of about 12 hours 25 minutes in the main, the two tidal generation forces anywhere in the earth arising by the moon are equal in general when the moon lies to the two sites of 180° interval of local mean lunar time. Based on this phenomenon the values Δ τ of τ 1- τ 2 or τ 1-τ 2±180° of two earthquakes occurring repetitiously in the same place are also calculated. The calculated results show that if the fault trends of the two earthquakes in the same place is near, the Δ τ is usually smaller and if the fault trends of the two ones is not near, the Δ τ is usually larger and the distribution of the local mean lunar time τ of earthquakes in different places is dispersive even if fault trends of these earthquakes are near, and the τ does not concentrate on the lower and upper transit of the moon. The above phenomena clear up that the triggering earthquake of earth solid tide arisen by the moon is relative with the fault trends of earthquakes and we ought to think over the difference of environmental conditions of earthquake preparation of each seismogenic zone and can not make statistics to earthquakes in different places when we study the relation between the solid earth tide arisen by the moon and earthquakes.

2023年土耳其7.8级强震前烟台台伸缩仪EW分量固体潮畸变

烟台台SS-Y伸缩仪EW分量在2023-02-06土耳其7.8级双震前2 d开始出现急速上升趋势变化,至6日凌晨开始下降,变化幅度达到9.58×10^(-8),震后维持快速上升趋势。对烟台台伸缩仪稳定性指标进行判定,认为仪器稳定;经现场核实,未发现环境干扰。最后对烟台台伸缩仪固体潮畸变特征进行分析,并基于畸变与远震前兆的关系进行初步探讨。