Analytical investigation of the exact groundwater divide between rivers beyond the Dupuit–Forchheimer approximation

The groundwater divide is a key feature of river basins and significantly influenced by subsurface hydrological processes. For an unconfined aquifer between two parallel rivers or ditches, it has long been defined as the top of the water table based on the Dupuit–Forchheimer approximation. However, the exact groundwater divide is subject to the interface between two local flow systems transporting groundwater to rivers from the infiltration recharge. This study contributes a new analytical model for two-dimensional groundwater flow between rivers of different water levels. The flownet is delineated in the model to identify groundwater flow systems and the exact groundwater divide. Formulas with two dimensionless parameters are derived to determine the distributed hydraulic head, the top of the water table and the groundwater divide. The locations of the groundwater divide and the top of the water table are not the same. The distance between them in horizontal can reach up to 8.9% of the distance between rivers. Numerical verifications indicate that simplifications in the analytical model do not significantly cause misestimates in the location of the groundwater divide. In contrast, the Dupuit–Forchheimer approximation yields an incorrect water table shape. The new analytical model is applied to investigate groundwater divides in the Loess Plateau, China, with a Monte Carlo simulation process taking into account the uncertainties in the parameters.     

Analytical study of ground motion caused by seismic wave propagation across faulted rock masses

Studying seismic wave propagation across rock masses and the induced ground motion is an important topic, which receives considerable attention in design and construction of underground cavern/tunnel constructions and mining activities. The current study investigates wave propagation across a rock mass with one fault and the induced ground motion using a recursive approach. The rocks beside the fault are assumed as viscoelastic media with seismic quality factors, Q_p and Q_s. Two kinds of interactions between stress waves and a discontinuity and between stress waves and a free surface are analyzed, respectively. As the result of the wave superposition, the mathematical expressions for induced ground vibration are deduced. The proposed approach is then compared with the existing analysis for special cases. Finally, parametric studies are carried out, which includes the influences of fault stiffness, incident angle, and frequency of incident waves on the peak particle velocities of the ground motions.     

Canonical theory of perturbations using eccentric anomaly as independent variable

It is shown in this paper how to build a canonical transformation of variables, so that the eccentric anomaly becomes the new independent variable. In the case of eccentric elliptical orbits it changes the equations of motion so, that they can be integrated analytically to any order of approximation comparatively easy.     

A note concerning the 2:1 and the 3:2 resonances in the asteroid belt

The dynamical behavior of asteroids inside the 2:1 and 3:2 commensurabilities with Jupiter presents a challenge. Indeed most of the studies, either analytical or numerical, point out that the two resonances have a very similar dynamical behavior. In spite of that, the 3:2 resonance, a little outside the main belt, hosts a family of asteroids, called the Hildas, while the 2:1, inside the main belt, is associated to a gap (the Hecuba gap) in the distribution of asteroids.In his search for a dynamical explanation for the Hecuba gap, Wisdom (1987) pointed out the existence of orbits starting with low eccentricity and inclination inside the 2:1 commensurability and going to high eccentricity, and thus to possible encounters with Mars. It has been shown later (Henrard et al.), that these orbits were following a path from the low eccentric belt of secondary resonances to the high eccentric domain of secular resonances. This path crosses a bridge, at moderate inclination and large amplitude of libration, between the two chaotic domains associated with these resonances.The 3:2 resonance being similar in many respects to the 2:1 resonance, one may wonder whether it contains also such a path. Indeed we have found that it exists and is very similar to the 2:1 one. This is the object of the present paper.     

Simulations of multimode bolometric interferometers

We describe a procedure for the numerical modelling of astronomical interferometers, with particular relevance to far-infrared and submillimetre wavelengths. The scheme is based on identifying a set of modes that carry power from the sky to the detector. The procedure is extremely general, and can be used to model scalar or vector fields, in any state of coherence and polarization, the only limitation being that the propagation of a coherent field through the system be described by an integral transform, a constraint that is in practise always met.
We present simulations of ideal, multimode two-dimensional interferometers, and show that the modal theory reproduces the correct behaviour of both Michelson and Fizeau interferometers. We calculate simulated visibility data for a multimode bolometric Michelson interferometer, with a synthesized source, and produce a dirty map, recovering the original source with the usual artefacts associated with interferometers.     

利用自回归模型检测地球极移重力效应

同时采用4个台站的国际超导重力仪长期连续重力观测资料和国际地球自转服务中心提供的同步地球自转参数，研究了极移引起的地球重力场变化特征。利用自回归模型估计了各序列的功率谱密度和积谱密度，结果表明极移导致的重力效应的主要能量集中在Chandler摆动和周年项附近，叠积后实际重力观测与极移重力信号理论值之间的差异分别为0．4％和3．9％，说明超导重力仪可有效监测极移导致的重力变化。     

行星磁场与卫星引力及其相对运动的统计关系

从知道某些天体具有磁场起，人们就对其磁场的起源提出种种解释，例如有电池说，转子说，化石说，发电机说等等。但由于这些学说都分别与某些观测事实相抵触而未被公认，因此星球磁场的起源一直是未能解决的命题。余先河先生提出，星球的磁场起源可能与所受的引力有关，行星的磁场强度正比于其所受卫星的引力；正比于卫星与行星的引力连线转动的相对角速度。本文对这两方面的命题分别进行了相关分析，得到其相关系数分别为：ｒ＝０．８４８１和ｒ＝０．８４２５，它们都在ａ＝０．０１的信度水平上相关。结果表明余新河关于行星磁场起源的设想是有基础的。本文还对其统计结果和可能机制进行了讨论。