Clear evidence for the sign-reversal of the pressure admittance to gravity near 3 mHz

The influence of the atmosphere on gravity measurements by Newtonian attraction and vertical displacements due to surface loading is well known and studied especially at frequencies below 1 mHz. Less work has been done at the higher seismic frequencies where inertial effects come into play. The sensor mass of a vertical accelerometer responds to several forces caused by the atmosphere, on global, regional and local scales. For the seismic frequencies the atmosphere above the observation site has the largest influence. A simple “Gedankenexperiment” demonstrates that the gravitational effect on one hand and the free air and inertial effects due to deformation on the other have opposite sign. Since the inertial effect is strongly frequency-dependent there should be a crossover-frequency where the pressure admittance to gravity changes sign. Simple analytical models clearly show this property near frequencies of a few mHz and are used here to amplify the variance reduction of the gravity residuals. The crossover-frequency depends on the properties of the models of the atmospheric phenomena and the elasticity of the Earth’s crust. Therefore in reality it must be expected to vary in time and space.We show one extremely clear example of this sign-reversal detected in the gravity data from the superconducting gravimeter GWR-C025 in Vienna among other examples which demonstrate the reality of this effect.     

A poincaré model for the earth's fluid core

Abstract

The superconducting-gravimeter data of Melchior and Ducarme (1986) has been interpreted as internal motion in the Earth's core by Aldridge and Lumb (1987) using a Poincaré model. Several low-order modes with periods of 13–16 hours have been tentatively identified in the core which is taken to be an incompressible, homogeneous fluid within a rigid, rotating container. The identification is based on asymptotic values of the frequencies which change slowly with time while the modes decay with an e-folding time of about 280 days. The slow change in frequency with time implies a small temporal variation in the rotation rate of the core. This mean flow is a nonlinear effect often observed in laboratory experiments designed to excite Poincaré modes. Interaction among modes during free ringdown is also observed in those experiments and apparently in the data of Melchior et al. (1988) as well. Laboratory work thus provides the link to extend the Poincaré model to include viscous and nonlinear effects in order to interpret the gravimetric observations as core modes.     

Can zonally symmetric inertial waves drive an oscillating zonal mean flow?

In the present paper zonal mean flow excitation by inertial waves is studied in analogy to mean flow excitation by gravity waves that plays an important role for the quasi-biennial oscillation in the equatorial atmosphere. In geophysical flows that are stratified and rotating, pure gravity and inertial waves correspond to the two limiting cases: gravity waves neglect rotation, inertial waves neglect stratification. The former are more relevant for fluids like the atmosphere, where stratification is dominant, the latter for the deep oceans or planet cores, where rotation dominates. In the present study a hierarchy of simple analytical and numerical models of zonally symmetric inertial wave-mean flow interactions is considered and the results are compared with data from a laboratory experiment. The main findings can be summarised as follows: (i) when the waves are decoupled from the mean flow they just drive a retrograde (eastward) zonal mean flow, independent of the sign of the meridional phase speed; (ii) when coupling is present and the zonal mean flow is assumed to be steady, the waves can drive vertically alternating jets, but still, in contrast to the gravity wave case, the structure is independent of the sign of the meridional phase speed; (iii) when coupling is present and time-dependent zonal mean flows are considered the waves can drive vertically and temporarily oscillating mean flows. The comparison with laboratory data from a rotating annulus experiment shows a qualitative agreement. It appears that the experiment captures the basic elements of the inertial wave mean flow coupling. The results might be relevant to understand how the Equatorial Deep Jets can be maintained against dissipation, a process currently discussed controversially.     

Higher-mode force response in multi-story strongback-braced frames

Strongback-braced frames employ an essentially elastic steel truss, or strongback, that distributes demands more uniformly to delay or prevent story mechanisms. Because inertial forces are no longer limited by the formation of a story mechanism, strongback-braced frames can exhibit large elastic force demands, particularly in the higher modes. This paper characterizes the higher-mode force response of strongback-braced frames. Four-story archetypes were designed using nonlinear dynamic analyses to incorporate higher-mode force demands into the design process. The response of the archetypes was compared with that of reference buckling-restrained braced frames that were allowed to form story mechanisms. The force demands in the strongback were then described using equivalent-static forces to represent the inertial forces induced by the higher modes. Force demands in the strongback arise from a yielding first-mode ‘pivoting’ and elastic higher-mode ‘bending’ response. These higher-mode force demands are elastic, ill-constrained by the strength of the yield mechanism, and depend significantly on the choice of ground motion record used for the analysis. In remaining elastic in the higher modes, the strongback distributes demands more uniformly and mitigates the formation of story mechanisms. Consequently, design and analysis methods for strongback-braced frames need to include estimates for these near-elastic higher-mode force demands.     

On the completeness of inertial wave modes in rotating annular channels

An important unanswered mathematical question in the theory of rotating fluids has been the completeness of the inviscid eigenfunctions which are usually referred to as inertial waves or inertial modes. We provide for the first time a mathematical proof for the completeness of the inertial modes in a rotating annular channel by establishing the completeness relation, or Parseval’s equality, for any piecewise continuous, differentiable velocity of an incompressible fluid.     

重力波在中层大气温度波导中的传播模式研究

本文给出了重力波在中层大气温度波导中的导制传播模型,并在此模型的基础上详细讨论了重力波部分导制传播下的对称模式与非对称模式,导出了不同模式下相应的特征函数和色散方程,进一步用离散的方法对两类色散方程进行了求解;同时还利用二维全隐欧拉格式（FICE）对重力波在温度波导中的传播进行了模拟,模拟的结果也成功地展现了对称与非对称两种传播模式.研究表明,下边界的扰动能量在向上传播进入波导区域后被俘获,形成导制传播.不同周期的初始扰动,在波导内均会形成对称与非对称形式两种模式的导制传播,由于两者的行进速度不一致,最终会引起两种不同模式的分离.数值模拟中重力波的水平行进速度与线性模型预测值非常接近.波导中不同模式下重力波的水平波长与初始扰动的水平波长非常一致,然而波导区域内重力波的频率与初始扰动的频率无关,频率不同的初始扰动会激发出相同频率的重力波对称与非对称导制传播模式.这表明在确定的温度波导中,水平波数才是决定重力波传播特性的决定因素.进一步的分析显示,初始扰动的水平波数-频率分布越接近完全导制传播的色散关系时,温度波导中更易于生成以该种模式部分导制传播的重力波.     

Characteristics of disturbances in the atmosphere and oceans

The characteristics of the disturbances in the atmosphere and oceans and in other stably stratified and rotating fluids are analyzed according to their phase and group velocities. It is shown that both stable stratification and rotation augment the velocity of the sound waves, and that the internal gravity waves and inertial waves are mutually exclusive when the Brunt-Väisälä frequency is different from the Coriolis parameter. It is also shown that both the barotropic and the internal Rossby waves are well separated from the gravity waves and that they can be represented accurately by the quasi-geostrophic potential vorticity equation, even close to the equator, except for the one member withn=0 which is coupled with an eastward propagating gravity wave.     

Ringdown of inertial waves in a spheroidal shell of rotating fluid

The role of inertial waves in the dynamies of the Earth's fluid core has been investigated through laboratory experiments on a spheroidal shell of rotating fluid. In these experiments inertial waves of azimuthal wavenumber 1, Ekman number 0(10^?5), Rossby number 0(10^?1) were excited by precession of an inner spheroidal body. Proximity to resonance was achieved by adjusting the ratio of the frequency of precession of the inner body to the rotational speed of the container to be near the eigenfrequency of the inertial wave mode being studied. Once the system was near resonance the perturbation was stopped and ringdown records were obtained. Amplitude, eigenfrequency and decay rate were recovered simultaneously for the principal and neighbouring modes excited using an iterative linearized least squares procedure.The recovery of complex eigenfrequencies for non-axisymmetric inertial waves in this shell geometry has given experimental verification of their existence. For those waves of azimuthal wavenumber one, a significant nonlinear interaction among modes is inferred from the simultaneous recovery of neighbouring modes. Other non linear effects include a mean azimuthal flow which appears to be stable for the low spatial order modes studied. These results contrast with highly unstable mean flow found experimentally in similar experiments carried out in cylindrical geometry.     

Mode coupling analysis and differential rotation in a flow driven by a precessing cylindrical container

We present a theoretical weakly nonlinear analysis of the dynamics of an inviscid flow submitted to both rotation and precession of an unbounded cylindrical container, by considering the coupling of two Kelvin (inertial) waves. The parametric centrifugal instability known for this system is shown to saturate when one expands the Navier–Stokes equation to higher order in the assumed small precession parameter (ratio of precession to rotation frequencies) with the derivation of two coupled Landau equations suitable to describe the dynamics of the modes. It is shown that an azimuthal mean flow with differential rotation is generated by this modes coupling. The time evolution of the associated dynamical system is studied. These theoretical results can be compared with water experiments and also to some numerical simulations where viscosity and finite length effects cannot be neglected.     

On the basic structure of oceanic gravity currents

Results from numerical simulations of idealised, 2.5-dimensional Boussinesq, gravity currents on an inclined plane in a rotating frame are used to determine the qualitative and quantitative characteristics of such currents. The current is initially geostrophically adjusted. The Richardson number is varied between different experiments. The results demonstrate that the gravity current has a two-part structure consisting of: (1) the vein, the thick part that is governed by geostrophic dynamics with an Ekman layer at its bottom, and (2) a thin friction layer at the downslope side of the vein, the thin part of the gravity current. Water from the vein detrains into the friction layer via the bottom Ekman layer. A self consistent picture of the dynamics of a gravity current is obtained and some of the large-scale characteristics of a gravity current can be analytically calculated, for small Reynolds number flow, using linear Ekman layer theory. The evolution of the gravity current is shown to be governed by bottom friction. A minimal model for the vein dynamics, based on the heat equation, is derived and compares very well to the solutions of the 2.5-dimensional Boussinesq simulations. The heat equation is linear for a linear (Rayleigh) friction law and non-linear for a quadratic drag law. I demonstrate that the thickness of a gravity current cannot be modelled by a local parameterisation when bottom friction is relevant. The difference between the vein and the gravity current is of paramount importance as simplified (streamtube) models should model the dynamics of the vein rather than the dynamics of the total gravity current. In basin-wide numerical models of the ocean dynamics the friction layer has to be resolved to correctly represent gravity currents and, thus, the ocean dynamics.     

基于B spline和正则化算法的低轨卫星轨道平滑

本文提出了一个利用纯几何轨道和力模型的新算法来计算精确且相对平滑的卫星轨道. 该法将一个纯几何轨道表达为一个B spline的线性组合，线性组合的系数可以由最小二乘法估计获得. 力模型通过计算加速度来附加约束. 为了平衡几何轨道的点位误差和加速度的不精确，一个基于“广义交互确认（GCV，generalized cross validation）”的正则化算法运用其中. 由于B spline的本地控制性，该方法的计算效率相当高. 本文的数值分析表明了该法的有效性. 模拟计算的结论是：带加速度约束较不带加速度约束的平滑效果好. 力模型越精确，平滑的轨道就越精确. 三个月的CHAMP实测轨道数据处理结果表明，平滑后的轨道改进了重力场模型.     

The self-similar,non-linear evolution of rotating magnetic flux ropes

We study, in the ideal MHD approximation, the non-linear evolution of cylindrical magnetic flux tubes differentially rotating about their symmetry axis. Our force balance consists of inertial terms, which include the centrifugal force, the gradient of the axial magnetic pressure, the magnetic pinch force and the gradient of the gas pressure. We employ the “separable” class of self-similar magnetic fields, defined recently. Taking the gas to be a polytrope, we reduce the problem to a single, ordinary differential equation for the evolution function. In general, two regimes of evolution are possible; expansion and oscillation. We investigate the specific effect rotation has on these two modes of evolution. We focus on critical values of the flux rope parameters and show that rotation can suppress the oscillatory mode. We estimate the critical value of the angular velocity _crit, above which the magnetic flux rope always expands, regardless of the value of the initial energy. Studying small-amplitude oscillations of the rope, we find that torsional oscillations are superimposed on the rotation and that they have a frequency equal to that of the radial oscillations. By setting the axial component of the magnetic field to zero, we study small-amplitude oscillations of a rigidly rotating pinch. We find that the frequency of oscillation is inversely proportional to the angular velocity of rotation ; the product being proportional to the inverse square of the Alfvén time. The period of large-amplitude oscillations of a rotating flux rope of low beta increases exponentially with the energy of the equivalent 1D oscillator. With respect to large-amplitude oscillations of a non-rotating flux rope, the only change brought about by rotation is to introduce a multiplicative factor greater than unity, which further increases the period. This multiplicative factor depends on the ratio of the azimuthal speed to the Alfvén speed. Finally, considering interplanetary magnetic clouds as cylindrical flux ropes, we inquire whether they rotate. We find that at 1 AU only a minority do. We discuss data on two magnetic clouds where we interpret the presence in each of vortical plasma motion about the symmetry axis as a sign of rotation. Our estimates for the angular velocities suggest that the parameters of the two magnetic clouds are below critical values. The two clouds differ in many respects (such as age, bulk flow speed, size, handedness of the magnetic field, etc.), and we find that their rotational parameters reflect some of these differences, particularly the difference in age. In both clouds, a rough estimate of the radial electric field in the rigidly rotating core, calculated in a non-rotating frame, yields values of the order mV m⁻¹.     

Equation of a Barotropic Fluid on a Rotating Spherical Surface and its Inertial Manifold

The inertial manifold is a positive invariant set which exponentially attracts all the trajectories of a dissipative dynamic system. It was introduced for the purpose of studying the asymptotic behaviour of such systems. The initial infinitely dimensional dynamic system, generated by a partial differential evolution equation, can be projected on to it, in order to obtain the final system of ordinary differential equations (inertial equations). These equations then simulate the initial object. Although the inertial manifold is an object relatively simpler than the attractor (a very complicated set of non-integer dimension may be an attractor) it is more difficult to prove its existence than that of the attractor. The equation of a barotropic fluid on a rotating spherical surface is one of the examples of dissipative dynamic systems with an inertial manifold. This kindles the hope that also the equations of the dynamics of the real atmosphere will have an inertial manifold. The reduction of the sample system to this Lipschitz manifold of finite dimension thus justifies us in analysing the behaviour of the atmosphere on non-linear models of finite dimensions and few parameters, in a finite system of ordinary differential equations.     

Effects of inertial and kinematic interaction on seismic behavior of pile with embedded foundation

Effects of inertial and kinematic forces on pile stresses are studied based on large shaking table tests on pile-structure models with a foundation embedded in dry and liquefiable sand deposits. The test results show that, if the natural period of the superstructure, T_b, is less than that of the ground, T_g, the ground displacement tends to be in phase with the inertial force from the superstructure, increasing the shear force transmitted to the pile. In contrast, if T_b is greater than T_g, the ground displacement tends to be out of phase with the inertial force, restraining the pile stress from increasing. With the effects of earth pressures on the embedded foundation and pile incorporated in, pseudo-static analysis is conducted to estimate maximum moment distribution in pile. It is assumed that the maximum moment is equal to the sum of the two stresses caused by the inertial and kinematic effects if T_b<T_g or the square root of the sum of the squares of the two if T_b>T_g. The estimated pile stresses are in good agreement with the observed ones regardless of the occurrence of soil liquefaction.     

On the coupled torsional and sway vibrations of a class of shear buildings

This paper is concerned with the free vibrations of a restricted class of multi-storey shear buildings in which inertial coupling exists between the torsional and the two sway vibrations. The restrictions imposed are that (a) the shear centres of all storeys lie on a vertical straight line, (b) the principal axes of shear are in the same directions in all storeys, (c) the centres of mass of all floors lie on another vertical straight line, (d) the radius of gyration about the shear centre of every floor mass is the same and (e) the ratios of the two shear stiffnesses to the torsional stiffness do not vary from storey to storey. In consequence of the last restriction it is proved that the 3n natural frequencies, normal modes and generalized masses, where n is the number of storeys, are expressible very simply in terms of products of the three natural frequencies, normal modes and generalized masses of the single-storey, three-dimensional building formed by removing everything above the first floor, with the n natural frequencies, normal modes and generalized masses of a certain n-storey, two-dimensional shear frame. In the special case of a uniform building, a simple closed form solution, valid for any number of storeys, is given.     

Time-domain modeling of global ocean tides generated by the full lunisolar potential

Traditionally, ocean tides have been modeled in frequency domain with a forcing from selected tidal constituents. It is a natural approach; however, it implicitly neglects non-linearities of ocean dynamics. An alternative approach is time-domain modeling with a forcing given by the full lunisolar potential, i.e., all tidal waves are a priori included. This approach has been applied in several ocean tide models; however, some challenging tasks still remain, for example, assimilation of satellite altimetry data. In this paper, we introduce the assimilative scheme applicable in a time-domain model, which is an alternative to existing techniques used in assimilative ocean tide models. We present results from DEBOT, a global barotropic ocean tide model, which has two modes: DEBOT-h, a purely hydrodynamical mode, and DEBOT-a, an assimilative mode. The accuracy of DEBOT in both modes is assessed through a series of tests against tide gauge data which demonstrate that DEBOT is comparable to state-of-the-art global ocean tide models for major tidal constituents. Furthermore, as signals of all tidal frequencies are included in DEBOT, we also discuss modeling of minor tidal constituents and non-linear compound tides. Our modeling approach can be useful for those applications where the frequency domain approach is not suitable.     

Ambient vibration test of an anchorage of South Bisan-seto suspension bridge

Dynamical properties of an anchorage of South Bisan-seto Suspension Bridge are determined by means of in situ measurement of the ambient vibrations. This anchorage, being located in the sea, is tall and slender; its natural frequencies are expected to be low. On the other hand, the elasticity of the ground estimated from the conventional geological data is highly dispersed. It even suggests the possibility of resonance between the towers and this anchorage. It is therefore necessary to know its natural frequencies precisely. Seismic motions and microtremors are acquired and used properly. The observation is repeated every year during the period of construction of the anchorage so that its inertial properties change significantly. First modes of the longitudinal and the transversal vibrations of the anchorage are identified and their frequencies determined. The elasticity modulus estimated from these natural frequencies is about five times as large as that of the geological survey while it is consistent with that obtained from the velocity of the elastic waves of the ground.     

Slow quasigeostrophic unstable modes of a lens vortex in a continuously stratified flow

This work addresses the linear dynamics underlying the formation of density interfaces at the periphery of energetic vortices, well outside the vortex core, both in the radial and axial directions. We compute numerically the unstable modes of an anticyclonic Gaussian vortex lens in a continuously stratified rotating fluid. The most unstable mode is a slow mode, associated with a critical layer instability located at the vortex periphery. Although the most unstable disturbance has a characteristic vertical scale which is comparable to the vortex height, interestingly, the critical levels of the successively fastest growing modes are closely spaced at intervals along the axial direction that are much smaller than the vortex height.     

The search for weak harmonic signals in a spectrum with application to gravity data

When considering the search for discovery or amplitude estimation of a spectral line with a probabilistic approach, great attention must be paid to the meaning of each step. We give the probability law for the amplitude of a spectral peak in the presence of random noise appearing in a periodogram and discuss the effective probability of the existence of the corresponding wave. We find that the estimated amplitude of a spectral peak is biased and should be corrected when the signal-to-noise ratio is small. As a first application to gravity data, it results in a re-estimation of the gravimetric amplitude factors (delta factors) provided by least-squares tidal analysis. We also estimate the probability of observing a spectral line above a given level in the spectrum of a purely random noise. This allows us to compute for given spectrum the number of peaks expected to overcross the classical levels used in statistical analysis (like nσ, where σ is the standard deviation of the temporal noise distribution and n is an integer with typical values equal to 2 or 3). A specific application to real data is investigating the gravity spectrum derived from a 5 year record of the French superconducting gravimeter and we show that the predicted statistics are indeed in agreement with the observations. We also show the statistical consequence of using longer observing periods to obtain the spectral estimations. The problem of detecting translational motion of the Earth's solid inner core (Slichter modes) in a gravity spectrum is analyzed and the probabilities of having a triplet of random peaks thresholding specific levels in a given frequency window are computed. We show that, in the case of a typical gravity spectrum (1 year of hourly data and a frequency window of 0.03 cycle h⁻¹), the probability of having a random set of three peaks exceeding a level of 3 σ, is very high. This emphasizes the need for a very careful analysis of spectral lines before inferring the existence of a true physical signal.     

Free and forced shallow-water oscillations in a rotating channel of parabolic cross-section

Free and forced oscillations of shallow water in an infinitely long rotating channel of parabolic cross-section are analyzed. The pure cross-channel oscillations ofChrystal (1905) and solutions for zero rotation first discussed byProudman (1925) andHidaka (1932) are special asymptotic solutions for the free modes of this model. However, for increasingly large, along-shore wave number, our solutions donot uniformly approach those ofReid (1958) andBall (1967) for a single shore-line and semi-infinite ocean. A method of computing eigen frequencies and eigen functions for the general problem is described, and a sufficient number of these are exhibited graphically to permit visualization of the transitions between the asymptotic regions.The forced problem consists of an incoming wave-train or surge generated at the center of the channel. Amplitude and transports near the shore are computed for a wide range of dimensionless incoming-wave frequencies and rotational frequencies.