Low-low satellite-to-satellite tracking: a comparison between analytical linear orbit perturbation theory and numerical integration

Low-low satellite-to-satellite tracking (ll-SST) range-rate observations have been predicted by two methods: one based on a linear perturbation theory in combination with the Hill equations, and one based on solving the equations of motion of two low-flying satellites by numerical integration. The two methods produce almost equivalent Fourier spectra of the range-rate observations after properly taking into account a few resonant terms. For a typical GRACE-type configuration, where the two satellites trail each other at a distance of 300 km at an altitude of 460 km, and in the presence of the EGM96 gravity field model, complete to spherical harmonic degree and order 70, the agreement between the Fourier spectra is about 1 mm/s compared to a root-mean-square (RMS) value of more than 220 mm/s for the range-rate signal. The discrepancy of 1 mm/s can be reduced significantly when not taking into account perturbations caused by the J₂ term. Excluding the J₂ term, the agreement between the two methods improves to 0.4 mm/s compared to a RMS value of 6 mm/s for the range-rate signal. These values are 0.01 and 2.3 mm/s when ignoring the spectrum for frequencies below two cycles per orbital revolution, reducing the discrepancy even further to about 0.5% of the signal. The selected linear perturbation theory is thus capable of modeling gravity field induced range-rate observations with very high precision for a large part of the spectrum.     

不同辐射传输方案对中尺度降水影响的对比分析

在MM5非静力稳定中尺度气象模式中引进了建立在δ-4流近似和相关-k分布基础上的对云水、雨水、冰晶和霰的辐射特性进行详细描述的辐射传输方案。新建立的辐射传输方案和MM5中原有的辐射传输方案在华南暴雨中的模拟结果相互比较，并与天气实况的对比表明:辐射在中尺度暴雨中起着重要的作用；辐射传输方案对云辐射特性描述的准确程度对于地面降水影响是明显的；不同的辐射传输方案对地面降水的影响存在较大的差异，并且这些差异在白天比在夜间明显；辐射传输过程对地面降水影响的差异主要表现在降水中心上，而对降水的地理分布改变很小；相对而言，不同的辐射传输方案之间对短波描述的差异对地面降水的影响明显，而对长波描述差异的影响不大；新辐射传输方案能够在一定程度上改进MM5对中尺度降水的模拟能力。     

大气热力学的若干新发展

对自20世纪80年代以来的大气热力学发展进行了综合评述, 阐述了如城市热岛效应、全球增暖等非传统大气热力学问题; 基于大气非平衡态热力学的非线性理论, 提出应发展突显热力学作用的数值模式, 并阐述其必要性和可能性。     

人造热带气旋垂直结构问题初探

在人造热带气旋中引入了倾斜垂直结构，并将其用于西折后9806号热带气旋登陆路径的数值模拟。结果表明：应用该方案所得到的模拟结果优于应用正压结构的人造热带气旋方案的模拟结果要准确。因此，此类热带气旋的预报中，其初始涡旋可采用倾斜垂直结构，从而提出改进热带气旋初始涡旋的一种新思路。若对其作进一步研究、推广，则具有广阔的业务应用前景。     

东亚南亚夏季降水对海温异常响应的数值模拟

利用IAP AGCM-Ⅱ大气环流模式，模拟了赤道中太平洋地区海温异常对东亚和南亚夏季（6～8月）降水的影响。结果表明：当1～4月该海温呈正异常时，中南半岛和我国东部降水将明显偏多，孟加拉湾附近降水偏少；而该海温呈负异常时，则我国华中和孟加拉湾附近降水偏多；对以上海温异常最敏感处为中南半岛和孟加拉湾附近地区，其中前者降水与该海温异常呈正相关，后者则呈负相关；该海温异常还会在中高纬地区激发出强迫波列。     

Groundwater-flow modeling in the Yucatan karstic aquifer, Mexico

The current conceptual model of the unconfined karstic aquifer in the Yucatan Peninsula, Mexico, is that a fresh-water lens floats above denser saline water that penetrates more than 40 km inland. The transmissivity of the aquifer is very high so the hydraulic gradient is very low, ranging from 7–10 mm/km through most of the northern part of the peninsula. The computer modeling program AQUIFER was used to investigate the regional groundwater flow in the aquifer. The karstified zone was modeled using the assumption that it acts hydraulically similar to a granular, porous medium. As part of the calibration, the following hypotheses were tested: (1) karstic features play an important role in the groundwater-flow system; (2) a ring or belt of sinkholes in the area is a manifestation of a zone of high transmissivity that facilitates the channeling of groundwater toward the Gulf of Mexico; and (3) the geologic features in the southern part of Yucatan influence the groundwater-flow system. The model shows that the Sierrita de Ticul fault, in the southwestern part of the study area, acts as a flow barrier and head values decline toward the northeast. The modeling also shows that the regional flow-system dynamics have not been altered despite the large number of pumping wells because the volume of water pumped is small compared with the volume of recharge, and the well-developed karst system of the region has a very high hydraulic conductivity. Electronic Publication     

Strike-slip fault bridge fluid pumping mechanism: insights from field-based palaeostress analysis and numerical modelling

We present a finite-element study of stress perturbation in evolving compressive and extensional strike-slip fault bridges. The results are compared with a fracture study of a compressive bridge at St Donats, UK. Horizontally interbedded calcareous mudstone and bioclastic calcilutite at St Donats have a distinct vertical permeability anisotropy. This sedimentary sequence behaves as a set of horizontal aquifers. The fluid flow in these aquifers is sensitive to mean stress gradients. Paleostress analysis of field fracture data, verified by finite-element modelling, indicates a rotation of σ₁ towards parallelism with boundary faults inside the growing compressive bridge. Boundary faults and bridge faults recorded numerous fluid flow events. The modelled mean stress pattern shows a regional maximum within the bridge and local maxima/minima pairs at boundary fault tips.Finite-element modelling of an extensional bridge indicates that σ₃ rotates towards parallelism with boundary faults. The mean stress pattern is similar to the pattern in compressive bridge but with maxima and minima locations interchanged. The stress patterns are reestablished by each stress build-up preceding the rupturation of the boundary faults throughout the development stages of strike-slip fault bridges. Mean stress gradients developed pre-failure control the fluid flow in fractures of the strike-slip fault system at and after the end of each stress build-up and the fluid flow in boundary faults post-failure. Fracture reactivation and new fracture generation within an evolving bridge is a process consisting of multiple successive events that retain the storage capacity of the bridge. Rupture and sealing of the main bounding-faults is a step-wise process that opens and closes fluid conduits between areas with different pressures.     

A HORIZONTAL 2-D HYDRAULIC NUMERICAL MODEL AND IT’S APPLICATIONS TO FLOOD FORECAST

1 INTRODUCTION There were many flood disasters in China in recent years. When the water level in a river is very high, weak parts of its dike may be destroyed resulted in the submersion of the protected land and properties. It is of significance for decision-makers to exactly predict the processes of flood propagation during flood control. There are many modes of dike bursting, such as seepage destroying by overflow on top of dike caused by dike body sinking induced by piping and soil fl…     

Hybrid vibration experiments with a bridge foundation system model

In order to improve seismic design technology of bridges, it is necessary to evaluate the vibration characteristics of a bridge–soil system that consists of soil, foundation structure, pier and superstructure. However, there have been few experimental studies on seismic behavior of bridge–soil system. In this paper, we conducted the hybrid vibration experiment on seismic behavior of bridge–soil system, and examined the applicability of hybrid vibration experiment to study seismic response of bridge–soil system. Based on the experiment results, seismic response of bridge was quantitatively studied.     

Efficient methods for solving water flow in variably saturated soils under prescribed flux infiltration

The non-linear solvers in numerical solutions of water flow in variably saturated soils are prone to convergence difficulties. Many aspects can give rise to such difficulties and in this paper we address the gravity term and the prescribed-flux boundary in the Picard iteration. The problem of the gravity term in the Picard iteration is iteration-to-iteration oscillation as the gravity term is treated, by analogy with the time-step advance technique, ‘explicitly’ in the iteration. The proposed method for the gravity term is an improvement of the ‘implicit’ approach of Zhang and Ewen [Water Resour. Res. 36 (2000) 2777] by extending it to heterogeneous soil and approximating the inter-nodal hydraulic conductivity in the diffusive term and the gravity term with the same scheme. The prescribed-flux boundary in traditional methods also gives rise to iteration-to-iteration oscillation because there is no feedback to the flux in the solution at the new iteration. To reduce such oscillation, a new method is proposed to provide such a feedback to the flux. Comparison with traditional Picard and Newton iteration methods for a wide range of problems show that a combination of these two proposed methods greatly improves the stability and consequently the computational efficiency, making the use of small time step and/or under-relaxation solely for convergence unnecessary.