Earth Rotation Parameters 1899.7–:1992.0 After Reanalysis Within The Hipparcos Frame

Earth rotation parameters (ERP) in the interval 1899.7–1992.0 are obtained from re-analysis of the observed latitude/universal time variations by optical astrometry. Hipparcos Catalogue is used to define the celestial reference frame, within which the ERP are described, with special care devoted to 'problematic' double and/or multiple stars. The terrestrial reference frame is defined by the adopted latitudes/longitudes of participating instruments and their secular motions as given by the NUVEL-1 NNR model of plate motions, and it is chosen to be very close to the International Terrestrial Reference Frame (ITRF). More than four million observations made with 48 different instruments at 31 observatories, located all over the world, are utilized to determine polar motion, celestial pole offsets and (after 1956) universal time UT1, all at 5-day intervals. Along with these parameters, the combinations of Love and Shida numbers, governing the tidal variations of the local verticals at individual observatories, are also determined. The analysis of the results covering almost a century, namely the long-periodic polar motion and length-of-day changes, is presented.     

天文观测极移运动周期变化的原因解析

宋贯一（1991,2006）相继发现日-地之间存在一种奇特的能量（动量）相耦合的自然现象:当太阳辐射光压作用于地表之后,地球表面物质的特殊物理性质会自然地把太阳辐射光压立刻分解为两部分,即P₁和P₂.其中P₁为全球各纬度带内单位海洋和陆地表面接收到的等值光压（P₁为随时间和纬度而变化的变量）.在一个回归年的时段内,由P₁在北、南两半球上的不平衡分布激发自转轴摆动,引起转动惯量的变化（该变化是目前所了解到的地球转动惯量的最大变化）并产生极移,该项极移运动的周期为12个月左右.由于P₁在北、南半球上的分布相对赤道是基本对称的和规则的,P₁对自转轴摆动的激发可称谓“规则性激发”;P₂为全球各纬度带内单位陆地和海洋表面接收到的光压值之差（P₂为随时间和纬度而变化的变量）.在一个回归年的时段内,由 P₂在北、南两半球上的不平衡分布激发自转轴摆动而产生的极移运动周期为14个月附近.由于P₂在北、南两半球上的分布相对赤道是不对称和不规则的, P₂对自转轴摆动的激发可称谓“不规则性激发”.这种天然存在的力源,恰恰是近百年来世界各国研究地球自转的地球物理学家渴望寻找的那种既能激发自转轴产生自由章动、又使其摆动含有两个不同周期（12个月和14个月）的激发力源.正是这一奇特自然现象的发现,才使长期以来困扰地球物理学领域内的极移、极移所包含两个周期的涨落变化及由此引起的地球自转速度变化等自然之谜得以破解成为可能. 本文作者仅对天文观测的极移运动周期及其极移运动所包含的两个周期分量在一定的范围内变化的成因作出了详细的解析,并得出如下结论： （1）极移运动主要是由太阳光压P₁ 和P₂共同激发引起的.天文观测的极移摆动周期的涨落变化是太阳光压（P₁＋P₂）激发自转轴摆动过程中,在空间上自转轴的摆动中心相对自转轴中心（地心）移动造成的,涨落范围在395～403±2天之间,即天文观测的极移运动的实际计算周期应在13.0～13.3个月之间变化.（2）极移所含的周年期摆动是由太阳光压P₁激发的.天文观测的周年期摆动周期涨落很小,变化于365.24～365.53天之间.在一个回归年内,由于日-地间距离的变化,使地球表面接受到的太阳辐射光能产生微小差异则是造成观测的周年摆动周期稍有拖长的原因.（3）极移所含的钱德勒摆动周期是由太阳光压P₂激发的.天文观测的钱德勒周期涨落较大,变化于426～437±2天之间,即实际计算周期应在14.0～14.4个月之间变化.观测的钱德勒摆动周期的变化是太阳光压P₂激发自转轴摆动过程中,在空间上自转轴的摆动中心相对自转轴中心（地心）移动造成的.上述的定量解析数据均得到实际观测资料的验证,为极移光压成因理论的正确性提供了具体详实的证据.     

Effect of the irregularities of shape of the Earth's core-mantle boundary on pole wandering

Summary Using a simple model, the motion of the Earth's axis of rotation due to local, axially asymmetric irregularities in the distribution of masses at the Earth's core-mantle boundary, is investigated. Under the assumption of linear displacement of the pole of inertia with time, expressions have been derived for the motion of the pole of the axis of rotation. A numerical estimate of the extent of this motion is made for the characteristic dimensions of the irregularities.     

Motion of the north magnetic pole within the scope of the dynamic model of the main geomagnetic field sources

The dipole model of the main geomagnetic field sources has been developed by the authors for several years. At present, the model includes 13 sources that existed and continuously developed during the 20th century. It has been assumed that the main dipole motion can be related to the motion of the Earth’s axis of inertia. At the same time, the known sharp changes in the direction of this motion, the so-called “wanderings” of the axis of inertia coincide in time with a change in the coordinates of the exit point of the main dipole magnetic moment vector on the Earth’s surface, dependent mostly on changes in the vector inclination. The motion of the north magnetic pole has been studied based on the model. It has been obtained that the dynamics of the main dipole parameters and, mainly, a stable variation in the inclination of the magnetic moment vector are responsible for the westward pole motion. At the same time, the observed rapid northward motion of the pole is related to the time variations in the parameters of 12 sources approximating the so-called nondipole part of the main field.     

空间局部加权回归自适应TFPF

时频峰值滤波(TFPF)算法是一种非常有效的去噪方法.但是传统的TFPF采用的单一窗长,并且仅沿时间方向进行滤波,忽略了信号的空间信息,并且TFPF近似等效成一个时不变的低通滤波器,不能追踪快速变化的信号.针对这些问题,引入空间局部加权回归自适应TFPF(SLWR-ATFPF).鉴于随机噪声在各个位置的方向随机性,以及有效信号在各个位置的方向确定性,首先利用空间局部加权回归(SLWR),对含噪信号进行空间加权,从而使加权之后的信号包含空间信息.然后,再引入凸集和Viterbi的思想,对空间加权之后的信号进行自适应滤波.从而,完成时空域二维自适应滤波.将SLWR-ATFPF应用于合成记录和实际的共炮点记录,实验结果表明,改进的方法与原算法相比,能够在压制低信噪比(SNR)随机噪声的同时更好地保留有效信号.     

On the period and quality factor of the Chandler wobble

The motion of the Earth’s geophysical pole is modeled by numerical integration of the Liouville equations. The dependence of variable parameters of these equations (the period of free nutation T and the quality factor of the shell Q) is examined as a function of time and initial data sets used for constructing geophysical perturbing functions. In particular, we used oceanic perturbing functions constructed from TOPEX/POSEIDON altimetry data. The annual and Chandler components of the pole motion were also modeled. Apart from the subtraction of quadratic trends, initial data were not filtered or averaged. The performed analysis provided model values of T = 425–440 days and Q = 20–60, and Q was found to be somewhat unstable with time.     

Amplitude-compensated Laplacian filtering of reverse time migration and its application

Reverse time migration is an advanced seismic migration imaging method. When the source wavefield and the receiver wavefield are cross-correlated, the cross-correlations of direct arrivals, backscattered waves and overturned waves will produce a lot of low-frequency noise, which will mask the final imaging results. Laplacian filtering, as a common method to suppress low-frequency noise, can adapt to any complex media, just adding a little computational cost. However, simple direct Laplacian filtering will destroy the characteristics of the useful signals. Therefore, the amplitude needs to be compensated before filtering when using the Laplacian filtering method. Zhang and Sun proposed an improved Laplacian filtering method and gave a simple calculation formula and explanation. This method can effectively suppress the low-frequency noise in reverse time migration while retaining the useful signal characteristics, but lacks detailed and strict mathematical derivation. Therefore, this paper gives a detailed and rigorous mathematical derivation of the amplitude-compensated Laplace filtering method from the point of view of amplitude-preserved filtering. The source wavelet is used instead of the source wavefield to compensate amplitude, just adding a little calculation cost. Finally, the amplitude-compensated Laplace filtering method is verified by two theoretical models and compared with the direct Laplacian filtering method.     

Effects of sunspot cycle length and CO2 on air temperature along Qinghai-Xizang railway and air temperature''s trend prediction

There are well coherences between annual averaged air temperatures at every meteorological station along the Qinghai-Xizang railway, and its 10-year moving average correlation coefficient is 0.92. Thus, the regional averaged annual mean temperature series along the Qinghai-Xizang railway (Trw) from 1935 to 2000 are constructed. The investigation is suggested that: Trw had significant responses to the 5-year lagged sunspot cycle length (SCL) and 15-year lagged concentration of atmospheric carbon dioxide (CO2), and the correlation coefficients between them are -0.76 (SCL) and 0.88 (CO2), respectively. The future SCL is predicted by the model of average generated function constructed with its main cycles of 76a, 93a, 108a, 205a and 275a. The result shows that the SCL would be becoming longer in the first half of the 21st century, and then it could be becoming shorter in the second half of the 21st century. Based on the natural change of SCL and the effect of double CO2 concentration, Trw in the 21st century is forecasted. It could warm up about 0.50℃ in the first half of the 21st century compared with the last decade of last century. The mean maximum air temperature could be likely about 0.20℃ in July and from 0.40℃ to 1.10℃ in January. The annual air temperature difference would likely reduce 0.3-1.00℃. The probability of above predictions ranges from 0.64 to 0.73.     

Which way do you lean? Using slope aspect variations to understand Critical Zone processes and feedbacks

Soil‐mantled pole‐facing hillslopes on Earth tend to be steeper, wetter, and have more vegetation cover compared with adjacent equator‐facing hillslopes. These and other slope aspect controls are often the consequence of feedbacks among hydrologic, ecologic, pedogenic, and geomorphic processes triggered by spatial variations in mean annual insolation. In this paper we review the state of knowledge on slope aspect controls of Critical Zone (CZ) processes using the latitudinal and elevational dependence of topographic asymmetry as a motivating observation. At relatively low latitudes and elevations, pole‐facing hillslopes tend to be steeper. At higher latitudes and elevations this pattern reverses. We reproduce this pattern using an empirical model based on parsimonious functions of latitude, an aridity index, mean‐annual temperature, and slope gradient. Using this empirical model and the literature as guides, we present a conceptual model for the slope‐aspect‐driven CZ feedbacks that generate asymmetry in water‐limited and temperature‐limited end‐member cases. In this conceptual model the dominant factor driving slope aspect differences at relatively low latitudes and elevations is the difference in mean‐annual soil moisture. The dominant factor at higher latitudes and elevations is temperature limitation on vegetation growth. In water‐limited cases, we propose that higher mean‐annual soil moisture on pole‐facing hillslopes drives higher soil production rates, higher water storage potential, more vegetation cover, faster dust deposition, and lower erosional efficiency in a positive feedback. At higher latitudes and elevations, pole‐facing hillslopes tend to have less vegetation cover, greater erosional efficiency, and gentler slopes, thus reversing the pattern of asymmetry found at lower latitudes and elevations. Our conceptual model emphasizes the linkages among short‐ and long‐timescale processes and across CZ sub‐disciplines; it also points to opportunities to further understand how CZ processes interact. We also demonstrate the importance of paleoclimatic conditions and non‐climatic factors in influencing slope aspect variations. Copyright © 2017 John Wiley & Sons, Ltd.     

Control Volume Method for the Dynamo Problem in the Sphere with the Free Rotating Inner Core

A model of thermally driven dynamo in the Boussinesq approximation in the spherical shell with the free rotating inner core is considered. To solve equations we use a new in dynamo modeling control volume technique (for details of this method for hydrodynamics see Patankar, 1980). The main advantage of this method over previous attempts to solve magnetohydrodynamics equations in the spherical grids is that no filtering of high harmonics in the pole regions is needed. We present the results of simulations for the self-consistent dynamo system evolution over the diffusion time and longer periods. Different ways of stabilizations of magnetohydrodynamics equations, when convective terms are of the same order (or larger) as conductive ones, are considered.     

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.     

Trends in the characteristics of the annual and semiannual variations observed in the radio wave absorption in the lower ionosphere

The continuous increase in concentration of greenhouse gases in the atmosphere is expected to cool higher levels of the atmosphere. There is some direct and indirect experimental evidence of long-term trends in temperature and other parameters in the mesosphere and lower thermosphere (MLT). Here we look for long-term trends in the annual and semiannual variations of the radio wave absorption in the lower ionosphere, which corresponds to the MLT region heights. Data from central and southeastern Europe are used. A consistent tendency to a positive trend in the amplitude of the semiannual wave appears to be observed. The reality of a similar tendency in the amplitude of the annual wave is questionable in the sense that the trend in the amplitude of the annual wave is probably induced by the trend in the yearly average values of absorption. The phases of both the annual and semiannual waves display a forward tendency, i.e. shift to an earlier time in the year. A tentative interpretation of these results in terms of changes of the seasonal variation of temperature and wind at MLT heights does not contradict the trends observed in those parameters.     

Reconstruction of precipitation variation from tree rings in recent 1000 years in Delingha, Qinghai

A major attempt of the studies on past global changes (PAGES) and climate variability and predict- ability (CLIVAR) is to reconstruct climate change us- ing high-resolution proxies[1―3]. Tree-ring data have played an important role in such studies. To date, tree rings have been extensively used to reconstruct tem- perature variations in the recent 1000 and even 2000 years of the Northern Hemisphere[4―7], to assess the effect of volcanic eruptions on temperature varia-tions[8], and to in…     

地震波双模态时变修正Kanai-Tajimi非平稳随机模型的改进及参数识别

在Vlachos等提出的双模态时变修正Kanai-Tajimi功率谱模型及其参数识别方法的基础上,利用杜修力等提出的Kanai-Tajimi功率谱滤波方法并引进遗传算法及二次优化识别技术进行改进,建立地震动时变功率谱的参数模型化方法。通过集集地震波的时变功率谱模型参数识别及模拟地震动算例,验证改进后的双模态时变修正Kanai-Tajimi功率谱模型的可行性和有效性,其方法可运用到重大工程结构抗震分析的设计地震动输入中。     

A Late Eocene paleomagnetic pole for the Pacific plate

Seamount magnetic anomaly inversions as well as DSDP paleomagnetic and equatorial sediment facies data constrain a paleomagnetic pole for the Pacific plate of Late Eocene age. The location of the pole at 77.5°N, 21.2°E implies 12.5 ± 1.6° of apparent polar wander for the Pacific plate during the last 41 ± 5 m.y. The Late Eocene pole is significantly different from the Pacific Maastrichtian pole at the 95% confidence level and indicates 7.2° of apparent polar motion of the Pacific between 69 and 41 m.y. B.P. The data source locations for the Late Eocene pole are scattered over a large area of the North Pacific and thus the consistency of the data supports the hypothesis that the north central Pacific plate has been rigid since the Eocene. The agreement of the Late Eocene pole with the motion predicted for the Pacific from hotspot models suggests that relative motion between the spin axis and hotspots has been small since that time. Additionally, this finding dictates that the significant amounts of hotspot versus spin axis motion inferred by other authors to have occurred since the Cretaceous must have instead occurred at a faster rate and concluded before the Eocene.     

Wavelet analysis of the Earth center of mass time series determined by satellite techniques

The Earth center of mass (geocenter) time series with the sampling interval of one week are determined from Satellite Laser Ranging (SLR), Global Navigation Satellite System (GNSS) and Doppler Orbitography and Radiopositioning Integrated by Satellite (DORIS) observations. The 3D geocenter time series were projected onto XY, YZ and ZX planes of the International Terrestrial Reference Frame (ITRF), thus, three complex-valued time series can be analyzed for each observation technique. The signal to noise ratio in these geocenter time series is very small and detectable oscillations are rather broadband, however, the annual oscillation can be noticed in each one of them. The wavelet transform technique with the Morlet wavelet function was applied to compute the mean and spectra-temporal polarization functions from the prograde (positive periods) and retrograde (negative periods) spectra of the examined complex-valued time series. The sign of the polarization function determines turning direction in the elliptical motion. If this function is positive or negative for oscillation with a chosen period, then this oscillation is prograde or retrograde, respectively. To estimate the significance level of polarization functions, corresponding to time series data length, the Monte Carlo experiment was performed using complex-valued white noise data. In order to detect similarity between elliptic oscillations in two different time series the spectra-temporal wavelet semblance function was computed. This function reveals that in the XY equatorial plane there is phase agreement between retrograde annual oscillation for SLR and GNSS techniques, and between prograde annual oscillation for DORIS and two other techniques. To construct a model of geocenter motion from GNSS and SLR center of mass time series the wavelet based semblance filtering method was applied. Common oscillations in the analyzed time series are dominated by the annual oscillation with amplitude less than 5 mm.     

Chinese tombs oriented by a compass: Evidence from paleomagnetic changes versus the age of tombs

Extant written records indicate that knowledge of an ancient type of compass in China is very old — dating back to before the Han dynasty (206 BC-220 AD) to at least the 4th century BC. Geomancy (feng shui) was practised for a long time (for millenia) and had a profound influence on the face of China’s landscape and city plans. The tombs (pyramids) near the former Chinese capital cities of Xi’an and Luoyang (together with their suburban fields and roads) show strong spatial orientations, sometimes along a basic South-North axis (relative to the geographic pole), but usually with deviations of several degrees to the East or West. The use of the compass means that the needle was directed towards the actual magnetic pole at the time of construction, or last reconstruction, of the respective tomb. However the magnetic pole, relative to the nearly ‘fixed’ geographic pole, shifts significantly over time. By matching paleomagnetic observations with modeled paleomagnetic history we have identified the date of pyramid construction in central China with the orientation relative to the magnetic pole positions at the respective time of construction. As in Mesoamerica, where according to the Fuson hypothesis the Olmécs and Maya oriented their ceremonial buildings and pyramids using a compass even before the Chinese, here in central China the same technique may have been used. We found a good agreement of trends between the paleodeclinations observed from tomb aligments and the available global geomagnetic field model CALS7K.2.     

Variable Chandler and Annual Wobbles in Earth’s Polar Motion During 1900–2015

The Chandler wobble (CW) and annual wobble (AW) are the two main components of polar motion, which are difficult to separate because of their very close periods. In the light of Fourier dictionary and basis pursuit method, a Fourier basis pursuit (FBP) spectrum is developed, which can reduce spectral smearing and leakage caused by the finite length of the time series. Further, a band-pass filtering method based on FBP spectrum (FBPBPF), which can effectively suppress the edge effect, is proposed in this paper. The simulation test results show that the FBPBPF method can effectively suppress the edge effect caused by spectral smearing and leakage and that its reconstruction accuracy at the boundary is approximately three times higher than the Fourier transform band-pass filtering method, which is based on Hamming windowed FFT spectrum, in extracting quasi-harmonic signals. The FBPBPF method is then applied to Earth’s polar motion data during 1900–2015. Through analyzing the amplitude and period variations of CW and AW, and calculating the eccentricity variation of the AW, we found that: (1) the amplitude of the CW is currently at a historic minimum level, and it is even possible to diminish further until a complete stop; and (2) the eccentricity of the AW has a gradually decreased fluctuation during the last 116 years.     

全球五大板块的运动和形变──用卫星激光测距数据导出的站速度分析

根据德国卫星激光测距（ＳＬＲ）数据分析中心ＧＦＺ对１９８０年１月－１９９１年７月获取的ＳＬＲ观测数据处理后得到的４１个ＳＬＲ站的站速度，解算了北美、欧亚、太平洋、南美和澳大利亚板块之间的相对运动欧拉矢量，得到了第１个ＳＬＲ实测的板块运动模型ＳＰＭＭ１．与地学板块模型ＲＭ２和ＮＵＶＥＬ－１的比较指出，ＳＰＭＭ１大体上与地学模型一致，与ＮＵＶＥＬ－１更为接近；ＳＰＭＭ１的欧亚与北美板块相对运动欧拉极与ＮＵＶＥＬ－１的相应极很接近，但旋转速率明显偏小．还分析了各板块上ＳＬＲ站的局部形变特征．欧亚板块东部和西部存在１０－２０ｍｍ／ａ的相对运动；板块边界附近ＳＬＲ站的残差站速度基本上反映了该边界的构造形变特征．