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

宋贯一（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₂激发自转轴摆动过程中,在空间上自转轴的摆动中心相对自转轴中心（地心）移动造成的.上述的定量解析数据均得到实际观测资料的验证,为极移光压成因理论的正确性提供了具体详实的证据.     

Response Surface Analysis and Statistical Modeling of Sulfide Generation from Municipal Wastewater

In this study, the interactive effects of four independent variables (initial chemical oxygen demand (COD) concentration, rotational velocity, temperature, and retention time) on the sulfide generation process in municipal wastewater were investigated. The process was modeled and analyzed with the variables in a series batch experiments. Experiments were conducted based on a central composite face‐centered design and analyzed using response surface methodology. The region of exploration for the process was taken as the area enclosed by COD (250 and 650 mg/L), rotational velocity (40 and 200 rpm), temperature (16 and 28°C), and retention time (2 and 40 h) boundaries. Eleven dependent parameters were either directly measured or calculated as response. The most influential variable on H₂S (g) production was found to be retention time. The results showed that simultaneous increase in temperature and rotational velocity caused an increasing trend in the amount of H₂S (g) emitted. Total sulfide was decreased by increasing rotational velocity due to more H₂S (g) secretion resulting from increased mixing rate. The present study provides valuable information about interrelations of quality and process parameters at different values of the studied variables.     

On the characteristics of ground motion rotational components using Chiba dense array data

An effort is made to examine the properties of rotational (torsional and rocking) ground motions using Chiba dense array data. The Chiba array system, located 30 km east of Tokyo, Japan, is composed of 15 boreholes with separation distances varying from 5 to 320 m. This provides a unique opportunity to examine the characteristics of rotational components. For this purpose, 17 events are considered and rotational ground motions are evaluated using spatial derivatives of translational ones. The effects of seismological parameters and separation distances between stations on properties of rotational motions are examined, showing a sudden increase in rotational motions for the earthquakes with large magnitude or PGA and decrease of these motions with increasing separation distance. While the duration of torsional motion is found to be larger than translational ones, there is no significant difference between durations of rocking and vertical motions. The effects of separation distance and earthquake magnitude on rotational response spectra are also investigated. The normalized rotational response spectra are found to be strongly affected by separation distance. The spectral ratios of rotational and translational motions are not linearly proportional to period as suggested by the previous studies. Finally, the torsional motion is predicted from translation ones for different separation distances at the site. The comparison of the predicted and the calculated torsional motions reveals a weak estimation in close separation distances (<30m) and satisfactory predictions in other cases. Copyright © 2007 John Wiley & Sons, Ltd.     

Apparent translational component for rotational ground motions

Even though the rotational ground motion may contribute significantly to the response of certain structures, their effects are generally ignored in seismic design, because of non-availability of appropriate instruments for direct recording of the rotational components. Like many others, a simplified framework was proposed by the authors elsewhere (Rodda and Basu in Int J Earthq Impact Eng 1(3):253–288, 2016) to extract the rotational motion as a temporal derivative of an apparent translational component (ATC) followed by scaling with an apparent velocity. ATC was defined such that its time derivative is closely correlated with the respective rotational motion. But the a priori knowledge of rotational motion is required in estimating the ATC for rocking component. An empirical procedure has been proposed here to bypass the requirement of rotational motion a priori. This paper also assesses the definition of ATC through examining the similitude between the time derivative of ATC and the respective rotational motion (benchmark) quantitatively. Similitude is assessed on smoothened response spectra (by Hamming window) of the time derivative of ATC and that of rotational motion. A new definition of spectral contrast angle (SCA) based on distance correlation has been proposed to assess the spectral similitude. To differentiate the similar from non-similar spectra, SCA corresponding to an acceptable degree of similarity is proposed by studying a large ensemble of ground motions from the PEER database. This similitude study is further extended using relative energy build up and energy spectra.     

On planar seismic wavefront modeling for estimating rotational ground motions: Case of 2-D SH line-source

At large hypocentral distances, it is convenient to approximate the curved transient seismic wavefronts as planar to estimate rotational ground motions from the single-station recordings of translational ground motions. In this paper, we investigate whether and when this approximation, referred to as the ‘plane-wave’ approximation, can be considered adequate close to the source. For this, we consider a simplistic source model comprising a two-dimensional, kinematic shear dislocation SH line-source buried in a homogenous, elastic half-space and assume this to be an equivalent representation of a finite-sized fault. The ‘plane-wave’ rotational motion is then synthesized from the exact translational motion solution to the assumed model and compared with the exact rotational motion solution for this model. The comparison between the two sets of rotational amplitudes in frequency domain suggests that the plane-wave approximation may be adequate, when the wavelengths of the seismic waves are much smaller than the source depth. When this is not true, the plane-wave approximation is seen to underestimate the Fourier amplitudes close to the source by several orders, particularly when the fault planes are vertically oriented. A similar comparison in the time domain indicates that a severe underestimation may also occur when the source rise time is longer than the shear-wave arrival time at the epicenter. Significant discrepancies are also observed between the waveforms of the exact and plane-wave rotational motions.     

Experimental and theoretical simulations of seismic poundings between two adjacent structures

Shaking table tests have been carried out to investigate the pounding phenomenon between two steel towers of different natural frequencies and damping ratios, subject to different combinations of stand‐off distance and seismic excitations. Both harmonic waves and ground motions of the 1940 El Centro earthquake are used as input. Subjected to sinusoidal excitations, poundings between the two towers could appear as either periodic or chaotic. For periodic poundings, impact normally occurs once within each excitation cycle or within every other excitation cycle. A type of periodic group poundings was also observed for the first time (i.e. a group of non‐periodic poundings repeating themselves periodically). Chaotic motions develop when the difference of the natural frequency of the two towers become larger. Under sinusoidal excitations, the maximum relative impact velocity always develops at an excitation frequency between the natural frequencies of the two towers. Both analytical and numerical predictions of the relative impact velocity, the maximum stand‐off distance, and the excitation frequency range for pounding occurrences were made and found to be comparable with the experimental observations in most of the cases. The stand‐off distance attains a maximum when the excitation frequency is close to that of the more flexible tower. Pounding appears to amplify the response of the stiffer structure but suppress that of the more flexible structure; and this agrees qualitatively with previous shaking table tests and theoretical studies. Copyright © 2003 John Wiley & Sons, Ltd.     

阿克苏断层氢气浓度动态特征及其映震效能

首先对阿克苏断层氢气浓度观测资料的周期成分进行了时频分析,然后分别利用线性回归和互相关方法分析了气温、气压与氢气浓度年周期成分的关系,并利用Molchan图表法对阿克苏断层氢气的映震效能进行了定量检验。结果显示:① 阿克苏断层氢气浓度具有清晰的年周期和半日周期成分;② 阿克苏断层氢气浓度的年周期动态变化与气温具有显著的正相关性,与气压的相关性一般,表明其年周期变化主要受气温影响;③ 阿克苏断层氢气浓度异常对观测点周边中强以上地震的反应较灵敏,映震效果较好。      

川滇地区地震活动与地球自转速率变化之间的关系

利用川滇地区1962—2016年间发生的MS≥6.0地震资料,在去除余震的情况下,分析了其发生与地球自转速率变化极值点之间的关系。得到以下结果:川滇地区MS≥6.0地震的发生与地球自转速率变化极值点具有显著关系,对于不同震级范围的地震,其与不同周期的地球自转速率变化有关。对于MS≥7.0的地震,有90%发生在地球自转速率季节性变化极值点前14天至后37天时间段内;对于6.0≤MS≤6.4的地震,有80%发生在地球自转速率短周期变化极值点前68小时至后30小时的时间段内;对于6.5≤MS≤6.9的地震,有87.5%发生在地球自转速率短周期变化极值点前36小时至后64小时的时间段内。地震发生在特定时间段的显著性检验结果表明,川滇地区MS≥6.0地震与地球自转速率变化极值点的关系都可以在α=0.05的显著性水平下通过显著性检验。这个结果表明,在地球自转速率发生转折的期间容易触发地震,对川滇地区地震发生时间预测具有参考意义。     

Early to Middle Miocene rotational tectonics of the Ou Backbone Range,northeast Japan

It is well known that a counterclockwise rotation occurred in the Miocene in northeast Japan. However, the detailed timing and mechanism of the rotation has been debated. Moreover, there has been no research about the relationship between rotational tectonics and the evolution of sedimentary basins. We carried out paleomagnetic and rock magnetic analyses in Nishiwaga Town, Iwate Prefecture, northeast Japan, where the stratigraphy and sedimentary basin formation have previously been clarified. We found that there was a counterclockwise rotational movement of about 45° at about 15 Ma. From our results and previous studies on the tectonics and sedimentary basin development, we are able to ascertain the following tectonic history and sedimentary basin evolution in this area: (i) before the rotational movement, sandstone and mudstone were deposited in a tranquil environment with no volcanic activity coupled with slow tectonic subsidence; (ii) between 16.4–15.1 Ma and 14 Ma, a counterclockwise rotation occurred with rapid tectonic subsidence and continuous explosive volcanism; (iii) at about 14 Ma, the counterclockwise rotation ended and there was a reduction in both subsidence and volcanism. This result shows the impact that rotational tectonics can have on sedimentary basin formation.     

太阳绕太阳系质心运动22年周期及其与太阳活动的联系

对于太阳活动22年周期的成因机制长期存在着争论.本文借助于行星会合指数以及开普勒第三定律,对太阳绕太阳系质心运动周期进行了分析计算.结果发现,太阳绕太阳系质心运动存在22.1826年显著周期,这与太阳磁场变化的22.20年周期相吻合.并从太阳系角动量守恒的角度解释了两者之间的成因联系:在太阳绕太阳系质心运动的准22年周期中,太阳系质心与太阳质心逐步接近而后逐步分离.当两个质心之间的距离接近零的时候,太阳轨道角动量与自转角动量叠加,会导致太阳自转角速度的加快;当两个质心之间的距离逐渐远离的时候,则导致太阳自转角速度的减慢.这可能是引发太阳活动和太阳磁场变化的原因.这一新认识为太阳活动准22年周期成因机制的解释提供了新的线索和依据.     

Seismic rotational displacement of gravity walls by pseudo-dynamic method: Passive case

Prediction of the seismic rotational displacements of retaining wall under passive condition is an important aspect of design in earthquake prone region. In this paper, the pseudo-dynamic method is used to compute the rotational displacements of rigid retaining wall supporting cohesionless backfill under seismic loading for the passive earth pressure condition. The proposed method considers time, phase difference and effect of amplification in shear and primary waves propagating through both the backfill and the retaining wall. The influence of ground motion characteristics on rotational displacement of the wall is evaluated. Also the effects of variation of parameters like wall friction angle, soil friction angle, amplification factor, shear wave velocity, primary wave velocity, period of lateral shaking, horizontal and vertical seismic accelerations on the rotational displacements are studied. The rotational displacement of the wall increases substantially with increase in amplification of both shear and primary waves, time of input motion, period of lateral shaking and decreases with increase in soil friction angle, wall friction angle. The rotational displacements of the wall also increase when the effect of wall inertia is taken into account. Results are provided in graphical form.     

THEORETICAL AND EXPERIMENTAL APPROACHES TO THE GEOPHONE SPURIOUS FREQUENCY*

The geophone spurious frequency is modeled as the resonance of the planar motion of a spider spring carrying a moving mass. An analytic solution is found using the Castigliano method by assuming that the spring arm is a single-mode vibrating cantilever beam. The spring shape is found from this analysis. When the typical spring has a circumferentially varying cross-section, the Castigliano method is no longer applicable. A dimensional analysis is used as an approximate method for general design. Based on the theoretical result, a rotational fixture and a translational fixture were designed for experimental purposes. A low-noise fixture and a phase-averaging technique provide the amplitude and spurious resonance in the frequency domain. Finally, a test is run by using a rotational fixture to compare with the approximate method of spurious frequency prediction. Very good agreement between prediction and experiment is found.     

Solar Irradiance Models and Measurements: A Comparison in the 220–240?nm wavelength band

Solar irradiance models that assume solar irradiance variations to be due to changes in the solar surface magnetic flux have been successfully used to reconstruct total solar irradiance on rotational as well as cyclical and secular time scales. Modelling spectral solar irradiance is not yet as advanced, and also suffers from a lack of comparison data, in particular on solar cycle time scales. Here, we compare solar irradiance in the 220–240?nm band as modelled with SATIRE-S and measured by different instruments on the UARS and SORCE satellites. We find good agreement between the model and measurements on rotational time scales. The long-term trends, however, show significant differences. Both SORCE instruments, in particular, show a much steeper gradient over the decaying part of cycle 23 than the modelled irradiance or that measured by UARS/SUSIM.     

Equivalent time-average and effective medium for periodic layers

The propagation of acoustic waves through a periodic layered medium is analyzed by an eigenvalue decomposition of the propagator matrix. This reveals how the velocity and attenuation of the layered medium vary as function of the periodic structure, material parameters and frequency. There are two important parameters which control the wave propagation in the periodic medium: the reflection coefficient and the ratio between one‐way traveltimes of the two parts of the cyclic layered medium. For low frequencies (large values of wavelength to layer thickness), the layered structure behaves as an effective medium, then there is a transition zone, and for higher frequencies (small values of wavelength to layer thickness) the medium is described by the time‐average velocity. In this paper we mostly concentrate on the transition zone between an effective medium and time‐average medium regimes. The width of the transition zone increases with larger values of the reflection coefficient. The transition zone corresponds to a blocking regime for which the transmission response of the layered structure is close to zero. For even higher frequencies, the time‐average medium is replaced by a new transition zone, and then again a time‐average medium. This pattern is periodically repeated with higher frequencies. For small values of the reflection coefficient, the transition between effective medium and time‐average medium occurs around a value of wavelength to layer thickness equal to 4.     

A model study of the temporal behaviour of the emission intensity and rotational temperature of the OH Meinel bands for high-latitude summer conditions

The temporal variation of OH* emission and weighted rotational temperature has been studied for high-latitude summer conditions. Observations for 60°N latitude show OH weighted temperatures that always exceed 145 K even during periods of noctilucent clouds. Using a one-dimensional model the effects in excited OH concentration produced by changes in temperature, eddy diffusion, and water concentration have been analysed. We are forced to conclude that there remains a discrepancy between the OH temperatures predicted by the model and that obtained from OH^* measurements. An increase in OH^* concentration from June to the beginning of August, followed by a slow decrease during August has been obtained in agreement with the measurements. The 16-day modulation present in the measurements was simulated in a simple manner by varying the temperature in the mesopause region. This variation produces periodic modulations in both OH^* concentration and weighted temperature of 16 days. The results show the temperature leading the OH^* column concentration by three days. This phase shift is also present in the observations.     

Long‐term variations and temporal scaling of hydroclimatic time series with focus on the German part of the Elbe River Basin

Long‐term variations and temporal scaling of mean monthly time series of river flow, precipitation, temperature, relative humidity, air pressure, duration of bright sunshine, degree of cloud cover, short wave radiation, wind speed and potential evaporation within or in vicinity of the German part of the Elbe River Basin are analyzed. Statistically significant correlations between the 2–15 year scale‐averaged wavelet spectra of the hydroclimatic variables and the North Atlantic Oscillation‐ and Arctic Oscillation index are found which suggests that such long‐term patterns in hydroclimatic time series are externally forced. The Hurst parameter estimates (H) based on the Detrended Fluctuation Analysis (DFA) indicate persistence for discharge, precipitation, wind speed, air pressure and the degree of cloud cover, all having an annual cycle and a broad low‐frequency distribution. Also, DFA H parameter estimates are higher for discharge than for precipitation. The major long‐term quasi‐periodic variability modes of precipitation detected using Singular Spectrum Analysis coincide with those detected in the discharge time series. Upon subtraction of these low‐frequency quasi‐periodic modes, the DFA H parameter estimates suggest absence of the persistence for both precipitation and discharge. Copyright © 2013 John Wiley & Sons, Ltd.     

Weakly nonlinear approximation of periodic flow in phreatic aquifers

Published frequency-domain solutions of periodic flow in aquifers apply strictly to mathematically linear systems that arise when aquifer diffusivity is assumed constant in space and time. This assumption can be invalid in phreatic aquifers that experience spatiotemporal variation in the free surface position and consequent variation in saturated thickness. A weakly nonlinear approach to formulating and solving periodic flow problems in the frequency domain can be applied in situations where conventional linearized approximations break down. The weakly nonlinear equations provide robust approximations of the true nonlinear response and require much less computational effort and time to solve than the full nonlinear problem. Nondimensional rules of thumb are presented for choosing between linear, weakly nonlinear and nonlinear solution strategies.     

相邻偏心高层建筑物间的振动台碰撞试验

地震引起的相邻建筑物间的碰撞过去多次强烈地震中出现,并造成了较为严重的结构破坏。为研究此问题,设计建造了两座缩尺钢质模型来模拟香港两座相邻的21层建筑,对其进行了一系列振动台碰撞试验,试验结果表明复杂结构间的碰撞可能是周期性的、群周期性的或混沌的,这主要取决于间距、结构参数以及地震动参数。碰撞引起的能量传递可能对较柔较轻的邻近建筑造成较严重的地震反应和破坏。     

Composite periodic foundation and its application for seismic isolation

This paper proposes a novel multidimensional composite periodic foundation for seismic isolation. The composite periodic foundation achieves multidimensional attenuation by innovative arrangement of periodic structures and taking advantage of the directional attenuation zone of periodic structures. Directional attenuation zones of periodic structures are derived for the in‐plane wave, and the impact of geometrical parameters of the periodic structure on the characteristics of the directional attenuation zones is studied. The effectiveness of the proposed composite periodic foundation is demonstrated through application in seismic isolation for nuclear power plant structures. Harmonic analysis and time history analysis results show that the proposed composite periodic foundation with low‐frequency directional attenuation zones can effectively reduce vibrations of the upper structure in both horizontal and vertical directions.     

Large wavenumber convection in the rotating annulus

Abstract

The annulus model considers convection between concentric cylinders with sloping endwalls. It is used as a simplified model of convection in a rapidly rotating sphere. Large azimuthal wavenumbers are preferred in this problem, and this has been exploited to develop an asymptotic approach to nonlinear convection in the annulus. The problem is further reduced because the Taylor-Proudman constraint simplifies the dependence in the direction of the rotation vector, so that a nonlinear system dependent only on the radial variable and time results. As Rayleigh number is increased a sequence of bifurcations is found, from steady solutions to periodic solutions and 2-tori, typically ending in chaotic behaviour. Both the magnetic (MHD convection) and non-magnetic problem has been considered, and in the non-magnetic case our bifurcation sequence can be compared with those found by previous two-dimensional numerical simulations.