由热带气旋MATTHEW引起的非潮汐负荷效应影响

热带气旋能在短期内造成海水和大气质量的重新分布,使得近海地表受力发生变化,进而产生非潮汐负荷形变,对现今高精度大地测量的影响已不容忽视.为了保证空间大地测量结果的精度和稳定性,热带气旋引起的地表形变必须进行有效的估计.因此本文联合NOS、GLOSS验潮站数据与海潮模型,通过获取非潮汐残余量分析了热带气旋“MATTHEW”引起的风暴潮.利用ECCO海洋环流模型、ERAin大气再分析模型、HUGO-m海洋动力学模型,分别估计了“MATTHEW”引起非潮汐海洋负荷、非潮汐大气负荷、动力学响应下非潮汐海洋负荷对地表位移的影响,结果表明热带气旋影响下的广大区域地表都不同程度受到非潮汐负荷的作用,最大位移分别达到-9.13 mm、3.31 mm、-6.11 mm,并且加入动力学响应的非潮汐海洋负荷要普遍大于IB(Inverted Barometer)响应下的结果.在对比不同位置站点所受负荷差异时,发现“大陆站”非潮汐海洋负荷形变普遍大于“岛屿站”,而“岛屿站”更易受非潮汐大气负荷的影响.     

利用GPS与环境负荷形变数据研究台风引起的垂向地表位移

台风造成的强降雨、低气压、海面高度变化均会引起地表的形变.本文利用中国大陆构造环境监测网(陆态网)7个GPS台站每日的垂向位移和环境负荷形变模型分析2018年9月10—26日台风"山竹"期间不同负荷引起的区域垂向地表形变.结果表明,台风期间大气负荷和非潮汐海洋负荷垂向形变最大分别达到5.1 mm和-9.2 mm.模型能较好地反映河流区域地表水文负荷变化造成的垂向形变,但不同模型之间存在系统偏差.由于缺少地下水等信息,模型反映负荷长期形变效应的效果不佳,且形变的量级明显小于GPS观测的结果.迅速增加的水文负荷使北海GPS站从开始下沉到最低点(-15.6 mm)5天的下沉量达到25.7 mm;珠海、广州GPS站均观测到河流汇水作用造成地表的二次下沉,且珠海站一周后才抬升到正常位置;湛江和北海GPS站能较好地反映河流水位变化,相关系数分别为-0.66和-0.50.研究结果表明,相比于形变模型,GPS能更有效地监测台风短期水文负荷形变,可为台风洪水等灾害监测与预报提供有用的信息.     

非构造形变对GPS连续站位置时间序列的影响和修正

GPS观测得到的地壳形变场通常包含有构造形变与非构造形变二类信息, 去除其中的非构造形变信息对于有效运用GPS数据研究构造形变场至关重要. 本文运用国际卫星对地观测资料及各类地球物理模型, 定量计算海潮、大气、积雪和土壤水、海洋非潮汐4项负荷效应造成的地壳非构造形变, 并以此研究和修正这些非构造形变对中国地壳运动观测网络GPS基准站位置时间序列的影响. 研究发现此4项负荷效应, 特别是大气、积雪和土壤水, 对于测站垂向位置的影响显著. 通过模型改正可以使测站垂向位置的RMS降低～1 mm, 占其总量的～11％. 对于垂向时间序列的周年项部分, 这一改正可降低其振幅的37％. 研究还表明经过地球物理模型改正和周年、半周年谐波拟合改正的时间序列比起仅经过周年、半周年谐波拟合改正的时间序列更为平滑, 表明地球物理模型改正对于消除非构造形变场的作用不是周年、半周年谐波拟合改正所能替代的.     

Determination of vertical displacements over the coastal area of Korea due to the ocean tide loading using GPS observations

This paper describes the GPS applicability for detecting the vertical displacements of ground stations caused by ocean tide loading effects. An experiment was carried out using 12 permanent GPS stations located in the coastal area of Korea using data in the period 1 July until 26 August 2003. The relative height differences were calculated from hourly DGPS data processing based on the carrier-phase observation. The power spectra of the M₂ and N₂ constituents of ocean tide loading were derived using the CLEAN algorithm. The differential vertical displacements generated by the ocean tide loading effect are typically 3–25 mm in coastal area of the Korea. We compared the results from GPS with those of the ocean tide models, NAO.99Jb regional model and GOT00.2, FES99 global models. The M₂ (N₂) amplitude differences of vertical displacements between GPS and GOT00.2 is 1.22 ± 3.61 mm (1.01 ± 1.48 mm), and that of the M₂ (N₂) amplitude difference between GPS and FES99 is 0.04 ± 4.64 mm (0.64 ± 1.75 mm), whereas the M₂ (N₂) amplitude difference between GPS and NAO.99Jb are 0.05 ± 1.03 mm (0.86 ± 1.18 mm). The highest vertical displacements at the PALM station are found for 24.5 ± 0.7 mm from GPS observation, and 22.9 mm from the regional model NAO.99Jb and 13.17 and 10.00 mm from the global models GOT00.2 and FES99, respectively. These values show that the vertical displacements derived from GPS are in good agreement with those of the regional model NAO.99Jb around Korea, more than with the global models. This result indicated that GPS is an effective tool to measure the vertical displacement caused by the ocean tide loading effect in the coastal area, and we need to use the NAO.99Jb ocean tide model rather than the global ocean tide models in and around the Korean peninsula for position determination with permanent GPS installations. This work demonstrates that vertical displacement caused by the M₂ and N₂ constituents of ocean tide loading can be measured by carrier-phase DGPS.     

Study of the Influence of the Indirect Effect of Ocean Tides and of the Direct Deformation Effect of the Atmosphere on the Results of GPS Observations

The geodynamic interpretation of the results of GPS observations among the permanently operating stations which establish the frame for defining and conserving the coordinate system and simultaneously for describing the fundamental tendencies of dynamics in the area represented by them, should be deprived of the effects of geodynamic phenomena, which can be either periodic or short-term (even non-periodic) so that secular changes may be found. Some of these influences are included in scientific processing software, others have to be modelled later. The purpose of this project is to determine the influences of the indirect effect of ocean tides and the direct deformational effect of the atmosphere on the results of GPS observations, performed at stations Pecný (Czech Republic), Wettzell (FRG) and Graz (Austria), and also for the GPS stations of the Central Europe Regional Geodynamic Project (CERGOP) in the region of Central and East Europe. On the basis of analyses and results it is possible to claim that at the present accuracy of GPS observations the influence of ocean loading and atmospheric pressure effects for geodynamic campaigns in the region of the Czech Republic is negligible. For campaigns in larger regions (i.e. Central Europe) atmospheric corrections will have to be applied to height differences. These corrections, however, due to their periodicity have negligible influence on secular changes.     

Velocity field and tectonic strain in Southern Spain and surrounding areas derived from GPS episodic measurements

The goal of this paper is to study the velocity field and deformation parameters in Southern Spain and surrounding areas (Ibero-Maghrebian region) using GPS episodic measurements. Results are compared to those previously published as well as deformation parameters derived from seismic data. For this purpose, a geodetic GPS network of 12 stations was observed during eight field campaigns from 1998 to 2005 by the San Fernando Naval Observatory (ROA), Spain. Relative GPS velocities in the Gulf of Cadiz with respect to the stable part of Eurasia are ～4.1 mm/yr in a NW–SE to NNW–SSE direction. In the Betics, Alboran Sea and North of Morocco, velocities are ～4.4 mm/yr in a NW–SE direction, and they are ～2.3 mm/yr in a N–S direction in the eastern part of the Iberian Peninsula. These results are in agreement with the anticlockwise rotation of the African plate. GPS strain tensors are determined from the velocity model, to obtain a more realistic crustal deformation model. The Gulf of Cadiz is subjected to uniform horizontal compression in a NNW–SSE direction, with a rotation to N–S in the Alboran Sea and Northern Morocco. An extensional regime in a NW–SE direction, which rotates to W–E, is present in the Internal Betics area. In the Betic, Alboran Sea and North of Morocco regions we compare seismic deformation rates from shallow earthquakes with the determined GPS deformation rates. The comparison indicates a seismic coupling of 27%, while the remaining 73% might be generated in aseismic processes. Deformations measured in the Ibero-Maghrebian region with GPS could be interpreted in terms of either elastic loading or ductile deformation.     

Annual variation detected by GPS,GRACE and loading models

Most GPS coordinate time series, surface displacements derived from the Gravity Recovery and Climate Experiment (GRACE), and loading models display significant annual signals at many regions. This paper compares the annual signals of the GPS position time series from the Crustal Dynamics Data Information System (CDDIS), estimates of loading from GRACE monthly gravity field models calculated by three processing centers (Center of Spatial Research, CSR; Jet Propulsion Laboratory, JPL; GeoForschungsZentrum, GFZ) and three geophysical fluids models (National Center for Environmental Prediction, NCEP; Estimating the Circulation and Climate of the Ocean, ECCO; Global Land Data Assimilation System, GLDAS) for 270 globally distributed stations for the period 2003-2011. The results show that annual variations derived from the level-2 products from the three GRACE product centers are very similar. The absolute difference in annual amplitude between any two centers is never larger than 1.25 mm in the vertical and 0.11 mm in horizontal displacement. The mean phase differences of the GRACE results are less than ten days for all three components. When we correct the GPS vertical coordinate time series using the GRACE annual amplitudes using the products from three GRACE analysis centers, we find that we are able to reduce the GPS annual signal in the vertical at about 80% stations and the average reduction is about 47%. In the north and the east, the annual amplitude is reduced on 77% and 72% of the stations with the average reduction 32% and 33%. We also compare the annual surface displacement signal derived from two environmental models; the two models use the same atmospheric and non-tidal ocean loading and differ only in the continental water storage model that we use, either NCEP or GLDAS. We find that the model containing the GLDAS continental water storage is able to better reduce the annual signal in the GPS coordinate time series.     

基于区域水准和跨断层资料分析云南地区地壳形变

为了获取云南地区现今地壳形变和断层活动特征,收集了云南地区1992—2016年的2期精密水准观测资料,采用基于GNSS连续站垂直形变资料约束的平差方法,对区域整体地壳垂直形变特征进行了分析,结果显示滇西北区域地壳整体表现为隆升,垂向速率为2—3mm/a,滇南区域显示地壳沉降速率为-1—-2mm/a,并在下关—永胜和思茅—景谷一带存在地壳形变高梯度带。收集了该地区30多年的跨断层场地观测资料,计算了所跨断层的三维活动量,得出了断层现今活动状态和异常特征,认为滇西北和滇中南地区的场地多数表现出异常活动的现象,其中楚雄、剑川和永胜的监测断层现今表现出与地质背景相反的走滑和拉张特征,峨山、建水和通海场地所处断层在张压性质上与地质背景相反。结合区域水准和跨断层观测分析结果,从“面”和“线”上得出滇西北和滇中南地区存在地壳形变和断层活动异常,构造应力积累较高的结果,加之缺震背景显著,认为应关注该区域中强地震的危险性。     

Modelling atmospheric and induced non-tidal oceanic loading contributions to surface gravity and tilt measurements

We investigate the contribution of atmospheric and its induced non-tidal oceanic loading effects on surface time-varying gravity and tilt measurements for several stations in Western Europe. The ocean response to pressure forcing can be modelled accordingly to the inverted barometer, i.e. assuming that air pressure variations are fully compensated by static sea height changes, or using ocean general circulation models. We validate two runs of the HUGO-m barotropic ocean model by comparing predicted sea surface height variations with hundred tide-gauge measurements along the European coasts. We then show that global surface pressure field, as well as a barotropic high-resolution ocean model forced by air pressure and winds allow in most cases a significant reduction of the variance of gravity residuals and, to a smaller extends tilt residuals.We finally show that precise gravity measurements with superconducting gravimeters allow the observation of large storm surges, occurring in the North Sea, even for inland stations. However, we also confirm that the continental hydrology contribution cannot be neglected. Thanks to their specific sensitivity feature, only tiltmeters closest to the coast can clearly detect the loading due to these storm surges.     

焉耆盆地北缘和静逆断裂-褶皱带第四纪变形

焉耆盆地是塔里木盆地东北缘天山山间的重要坳陷区,盆地北缘发育的和静逆断裂-褶皱带是一条现今活动强烈的逆断裂-褶皱带,对其第四纪以来缩短量和隆升量的计算有利于分析该区域的构造活动情况,对缩短速率和隆升速率的估计可以与天山造山带其他区域的活动速率进行横向对比,从而反映出焉耆盆地在天山晚新生代构造变形的作用。在深部资料不足的情况下,对背斜形态完整、构造样式简单的和静逆断裂-褶皱带,利用地表可获得的地层和断层产状,通过恢复褶皱几何形态,计算褶皱的缩短量、隆升量和断层滑动量,得到逆断裂-褶皱带早更新世晚期(1.8Ma)、中更新世(780ka)和晚更新世中期(80ka)以来的缩短量分别为1.79km、0.88km和26m,初步估计的缩短速率分别为0.99mm/a、1.13mm/a和0.33mm/a。显示和静逆断裂-褶皱带自开始形成以来构造活动强度并不一致。与地壳形变观测结果对比,作为南天山东段最主要的坳陷区,焉耆盆地吸收了这一区域(86°～88°E)的大部分地壳缩短,且主要表现为盆地北缘新生逆断裂-褶皱带的强烈变形。     

鲜水河断裂带库仑应力演化与强震间关系

以鲜水河断裂带为研究区,首先验证了该断裂带上1893年以来M6.7以上地震的相互触发作用,然后采用更符合实际的分层黏弹介质模型研究强震震后黏滞松弛引起的库仑应力变化对后续地震的影响,并基于负位错理论计算鲜水河断裂带10个断层段的震间长期构造加载作用引起的断层上的应力积累.在此基础上,讨论同震、震后、震间效应引起的累积库仑应力变化与区域强震活动的关系,给出断层上库仑应力随时间的演化.结果表明,鲜水河断裂带上1893年以来发生的7次强震均在其前面一系列强震及构造应力加载的驱使下发生,同震、震后、震间三方面效应均引起了鲜水河断裂带不可忽略的库仑应力变化.对断层上的库仑应力状态的研究,可进一步为揭示地震的发生规律、找寻危险断层段提供线索.     

Thermo-hydrogeomechanical modeling of vertical ground deformation during the operation of the Mutnovskii Geothermal Field

Vertical ground displacements at the Mutnovskii Geothermal Field were measured annually during the 2004–2013 period using a network of borehole markers. Three areas have been identified where the vertical displacements were different in character. Positive vertical deformation (2–5 mm/yr) was recorded in the middle of the Dachnyi Area in 2005–2006, giving way to a later stabilization. The North Test Area where the exhaust heat carrier was reinjected did not show any substantial deformation during the 2003–2006 period. This was followed by an uplift in 2006–2008 (6–7 mm/yr) and then by a subsidence (5–8 mm/yr) in 2009–2013. No vertical deformation worth mentioning was recorded in the Verkhne-Mutnovskii Area prior to 2008, but a rapid subsidence began at a rate of 6–18 mm/yr after 2008. We used the TOUGH-FLAC software to analyze vertical ground deformation. This program is used to perform thermo-hydrogeomechanical (THM) modeling based on the previous TOUGH2 model. THM modeling can explain relative vertical deformation during the exploitation of the Mutnovskii Geothermal Field by separating the geothermal field into two compartments by the Osnovnoi fault, as the two have different tectonic settings and petrophysical properties.     

青海门源6.9级地震同震及震前GNSS变形特征分析

北京时间2022年1月8日1点45分,在我国青海省门源县发生了6.9级地震.通过震中附近陆态网络GNSS连续观测数据得到的同震位移场显示,距离震中最近的QHME站同震位移最大,EW向达到20.31 mm,SW向达到-35.45 mm,震中附近五个站的同震位移反映出此次地震的左旋同震破裂特征;GNSS站间基线时间序列结果...     

Precision and accuracy of GPS-derived station displacements

Space geodetic techniques like Global Navigation Satellite Systems (GNSS), Satellite Laser Ranging (SLR) and Very Long Baseline Interferometry (VLBI) provide valuable input for, e.g., studies of Global Isostatic Adjustment (GIA). This paper discusses the current precision and accuracy of GPS-derived vertical and horizontal station displacements. The precision is evaluated by repeatabilities and solutions computed from different subintervals of the data available. However, due to systematic effects, the precision is often much better than the accuracy. The accuracy is evaluated by comparisons of the space geodetic techniques amongst each other and comparisons with geophysical models for atmospherical and hydrological loading. Besides the analysis of time series, co-located GNSS, SLR, and VLBI sites allow for a comparison of velocities estimated in Terrestrial Reference Frame (TRF) solutions of the different techniques.     

Upper crustal deformation characteristics in the northeastern Tibetan Plateau and its adjacent areas revealed by GNSS and anisotropy data

The northeastern part of the Tibetan Plateau is a region where different tectonic blocks collide and intersect, and large earthquakes are frequent. Global Navigation Satellite System(GNSS) observations show that tectonic deformation in this region is strong and manifests as non-uniform deformation associated with tectonic features. S-wave splitting studies of near-field seismic data show that seismic anisotropy parameters can also reveal the upper crustal medium deformation beneath the reporting...     

Global surface mass from a new combination of GRACE,modelled OBP and reprocessed GPS data

Weekly surface loading variations are estimated from a joint least squares inversion of load-induced GPS site displacements, GRACE gravimetry and simulated ocean bottom pressure (OBP) from the finite element sea-ice ocean model (FESOM).In this study, we directly use normal equations derived from reprocessed GPS observations, where station and satellite positions are estimated simultaneously. The OBP weight of the model in the inversion is based on a new error model, obtained from 2 FESOM runs forced with different atmospheric data sets.Our findings indicate that the geocenter motion derived from the inversion is smooth, with non-seasonal RMS values of 1.4, 0.9 and 1.9 mm for the X, Y and Z directions, respectively. The absolute magnitude of the seasonal geocenter motion varies annually between 2 and 4.5 mm. Important hydrological regions such as the Amazon, Australia, South-East Asia and Europe are mostly affected by the geocenter motion, with magnitudes of up to 2 cm, when expressed in equivalent water height.The chosen solar radiation pressure model, used in the GPS processing, has only a marginal effect on the joint inversion results. Using the empirical CODE model slightly increases the annual amplitude of the Z component of the geocenter by 0.8 mm. However, in case of a GPS-only inversion, notable larger differences are found for the annual amplitude and phase estimates when applying the older physical ROCK models. Regardless of the used radiation pressure model the GPS network still exhibits maximum radial expansions in the order of 3 mm (0.45 ppb in terms of scale), which are most likely caused by remaining GPS technique errors.In an additional experiment, we have used the joint inversion solution as a background loading model in the GPS normal equations. The reduced time series, compared to those without a priori loading model, show a consistent decrease in RMS. In terms of the annual height component, 151 of the 189 stations show a reduction of at least 10% in seasonal amplitude.On the ocean floor, we find a positive overall correlation (0.51) of the inversion solution with time series from globally distributed independent bottom pressure recorders.Even after removing a seasonal fit we still find a correlation of 0.45. Furthermore, the geocenter motion has a significant effect on ocean bottom pressure as neglecting it causes the correlation to drop to 0.42.     

PROGRESS IN APPLICATION OF GNSS TO DIVISION OF ACTIVE TECTONIC BLOCKS IN CONTINENTAL CHINA

Chinese scientists proposed that large earthquakes that occurred in mainland China are controlled by the movement and deformation of active tectonic blocks. This scientific hypothesis explains zoned phenomenon of seismicity in space. The active tectonic blocks are intense active terranes formed in late Cenozoic and late Quaternary, and the tectonic activity of block boundaries is the intensest. Global Navigation Satellite System(GNSS)has advantages of high spatio-temporal resolution, broad coverage, and high accuracy, and is utilized to monitor contemporary crustal deformation. High accuracy and resolution of GNSS velocity field within mainland China and vicinities provided by previous studies clearly demonstrate that different active tectonic blocks behave as different patterns of movement and deformation, and block interaction boundaries have intense tectonic deformation. The paper firstly introduces the GPS networks operated by the Crustal Movement Observation Network of China(CMONOC)since 1999, and GNSS data processing methods, including GAMIT, BERNESE and GIPSY/OASIS, and discusses the advantages of using South China block as a regional reference frame for GNSS velocity field, then proposes three strategies of block division, F-test, quasi-accurate detection(QUAD), and clustering analysis. Furthermore, we introduce rigid and non-rigid block motions. Rigid block motion can be denoted by translation and rotation, while non-rigid block motion can be described by rigid motion and internal strain deformation. Internal strain deformation can be divided into uniform and linear strains. We also review the usage of F-test to distinguish whether the block acts as rigid deformation or not. In addition, combining with recent GNSS velocity results, we elaborate the characteristics of present movement of rigid block, such as the South China, Tarim, Ordos, Alashan, and Northeast China, and that of non-rigid block, such as the Tibetan plateau, Tian Shan, and North China plain. Especially, the Tibetan plateau and Tian Shan seem to deform continuously with significant internal deformation. In order to enrich and perfect the active tectonic block hypothesis, we should carefully design dense GNSS networks in inner blocks and block boundaries, optimize utilizing other space geodesy technologies such as InSAR, and strengthen combining study of geodesy, seismogeology and geophysics. Through systematic summary, this paper is very useful to employing GNSS to investigate characteristics of block movement and dynamics of large earthquakes happening in block interaction boundaries.     

基于GNSS与InSAR约束的九寨沟MS7.0地震滑动模型及其库仑应力研究

2017年8月8日的九寨沟M_S7.0地震发生在岷江断裂、塔藏断裂及虎牙断裂交汇地区,地处青藏高原东北部的川甘交界地区,位于巴颜喀拉地块的东缘,地质构造复杂,对于九寨沟地震震中位置和发震断层的确定,存在不同意见.本文利用GNSS及升降轨InSAR观测,在获取九寨沟地震同震形变场的基础上,基于均匀弹性半无限位错模型,联合反演了发震断层的滑动分布模型,并计算了同震库仑应力变化.InSAR同震形变场显示,视线向最大沉降量和抬升量分别为0.21 m和0.16 m,形变场长轴为NW向,形变主要集中在断层西侧.距震中40 km和65 km的九寨和松潘两县,水平向的GNSS同震位移分别达14.31 mm和8.22 mm.联合GNSS和InSAR同震形变场反演得到的滑动分布主要集中在沿走向5~33 km,倾向2~20 km的范围内,平均滑动量为0.18 m,最大滑动量为0.91 m.发震断层长40 km,宽30 km,走向155°,倾角81°,滑动角-9.56°.同震位移场及滑移分布模型表明此次地震为一次左旋走滑为主的地震事件,地震破裂并未完全到达地表,与虎牙断裂北段的几何产状和运动学性质更为接近,结合精定位余震的分布,我们确定虎牙断裂北段为此次地震的发震断层,震中位于北纬33.25°,东经103.82°,震源深度10.86 km,矩震量为7.754×10¹⁸ Nm,相应的矩震级为M_W6.5,与美国地调局和哈佛大学给出的震源机制解基本一致.同震库仑应力导致了虎牙断裂北段延长线的东北和西南两端应力增强,其中塔藏断裂的罗叉段和马磨段未来强震的危险性值得关注.     

高频GNSS实时地震学与地震预警研究现状

为实现从注重灾后救助向注重灾前预防转变,如何提高地震灾害监测预警和风险防范能力成为我们关注的重点.本文给出了国际上GNSS位移记录、强震动加速度记录、测震速度记录在地震预警中的应用现状,并总结了各自的特点,归纳出围绕高频GNSS地震学在震级与破裂过程实时反演中的几个需要进一步研究的关键问题：（1）引入北斗系统,基于高频GNSS（GPS/BDS）双系统的实时位移解算方法来提高实时单站位移解算精度,使实时解算精度达到厘米级;（2）开展强震仪加速度记录基线偏移校正研究,弥补地震近场GNSS站密度不足问题;（3）强震仪加速度记录与GNSS位移记录特点不同,开展强震仪加速度数据与GNSS位移数据实时融合处理研究,快速获得包含丰富地震形变和速率的波形数据;（4）测震学方法可快速估算震级,但是在强震发生时会出现震级饱和现象,造成震级估算偏低.需要开展基于GNSS位移时间序列的多种方法相结合的实时震级估算方法研究,通过与地震学方法比较和结合,来得到精度高、计算快的震级估值算法;（5）基于高频GNSS、断层初始模型快速选取、断层尺度、参数自适应调整是快速判断断层破裂方向的基础,在断层破裂过程自适应准实时反演算法方面需要进一步加强.通过国内外研究现状调研、分析,表明基于高频GNSS地震学的震级快速确定、震源破裂过程准实时反演算法的发展将对我国地震预警系统从"二网融合"到"三网融合"提供坚实的技术支撑.     

On the Mechanism of the Secular Tidal Acceleration of the Solid Inner Core and the Viscosity of the Liquid Core

As is known, the secular deceleration of the Earth's diurnal rotation is explained mainly by the tidal friction in the ocean. Below we consider this mechanism in some detail, taking into account also elastic deformations of the mantle under the action of ocean loading and the interaction between the tide-generating body, ocean tidal wave, liquid outer core, and solid inner core. It is shown that elastic displacements of the core-mantle boundary under the action of ocean loading are of about the same amplitude and phase as the elastic loading displacements of the Earth's outer surface. As a result, side by side with the mechanism of secular deceleration of diurnal rotation of the mantle, there are also (1) the opposite mechanism of secular acceleration of diurnal rotation of the outer liquid core and of the solid inner core and (2) the mechanism of excitation of differential rotation in the liquid core. Taking these effects into account, we compare theoretical and modern observed data on the eastward drift of the solid inner core. It is shown that the best agreement may be obtained if the turbulent viscosity of the liquid core is about 2 × 10 ³ Poise