Nonlinearities and Linearities in Internal Gravity Waves of the Atmosphere and Oceans

The spectrum of internal gravity waves in the atmosphere and oceans is sufficiently intense that nonlinear interactions must occur, if these waves are analyzed in Eulerian coordinates as is usually done. As it happens, however, if these waves are analyzed in Lagrangian coordinates the most important nonlinearity can be entirely avoided: it is an Eulerian mathematical construct only, not a physical process. The mathematical basis for this assertion is developed here, and some of its consequences are discussed. Among the latter is a questioning of the validity of standard Eulerian eikonal methods of calculating ray paths and related functions in a multiwave environment, discussed in an appendix.     

重磁三维反演伊犁盆地中部密度和磁性结构

伊犁盆地中部区域地处多个构造单元的交汇区,构造复杂,缺乏高精度地质资料,制约了对盆地构造演化的进一步认识.本文联合使用多种先验信息作为约束条件,利用高精度重磁数据反演获得了研究区地下海拔-10 km以浅的三维密度和磁性结构,增强了反演结果的可靠性.结果揭示以白石墩次凸为中心发育了一个"北断南超"的南西-北东向不对称型凹陷,北部沉积厚度大于南部;区内主要断裂具有高角度特征,生烃中心受断裂控制,阿吾拉勒山前凹陷浅部可能发育逆掩断裂;火成岩以侵入的基性火成岩为主,伴有部分喷出火成岩和中酸性火成岩.     

石油重磁资料处理解释系统在微机上的实现

微机技术的发展和石油重磁勘探的需要促使了石油重磁资料微机处理和解释系统的建立。系统由输入、输出辅助子系统、重磁场分离和转换子系统及五个反演解释子系统八部分组成。五个反演解释子系统针对重磁力方法在油气勘探中所能解决的地质问题的类型不同而设计。计算模块的选择注意到方法的有效和适应，地球物理综合解释的需要及平面与剖面、定量正演与反演相结合的重要性。实际应用例子说明系统建立的成功和有效     

Efficient GOCE satellite gravity field recovery based on least-squares using QR decomposition

We develop and apply an efficient strategy for Earth gravity field recovery from satellite gravity gradiometry data. Our approach is based upon the Paige-Saunders iterative least-squares method using QR decomposition (LSQR). We modify the original algorithm for space-geodetic applications: firstly, we investigate how convergence can be accelerated by means of both subspace and block-diagonal preconditioning. The efficiency of the latter dominates if the design matrix exhibits block-dominant structure. Secondly, we address Tikhonov-Phillips regularization in general. Thirdly, we demonstrate an effective implementation of the algorithm in a high-performance computing environment. In this context, an important issue is to avoid the twofold computation of the design matrix in each iteration. The computational platform is a 64-processor shared-memory supercomputer. The runtime results prove the successful parallelization of the LSQR solver. The numerical examples are chosen in view of the forthcoming satellite mission GOCE (Gravity field and steady-state Ocean Circulation Explorer). The closed-loop scenario covers 1 month of simulated data with 5 s sampling. We focus exclusively on the analysis of radial components of satellite accelerations and gravity gradients. Our extensions to the basic algorithm enable the method to be competitive with well-established inversion strategies in satellite geodesy, such as conjugate gradient methods or the brute-force approach. In its current development stage, the LSQR method appears ready to deal with real-data applications.     

Gravimetric monitoring of the first field‐wide steam injection in a fractured carbonate field in Oman – a feasibility study

Gas‐Oil Gravity Drainage is to be enhanced by steam injection in a highly fractured, low permeability carbonate field in Oman. Following a successful pilot, field‐wide steam injection is being implemented, first of this type in the world. A dedicated monitoring program has been designed to track changes in the reservoir. Various observations are to be acquired, including, surface deformation, temperature measurements, microseismic, well logs, pressure and saturation measurements to monitor the reservoir. Because significant changes in the reservoir density are expected, time‐lapse gravimetry is also being considered. In this paper we investigate the feasibility of gravimetric monitoring of the thermally enhanced gravity drainage process at the carbonate field in Oman. For this purpose, forward gravity modelling is performed. Based on field groundwater measurements, the estimates of the hydrological signal are considered and it is investigated under what conditions the groundwater influences can be minimized. Using regularized inversion of synthetic gravity data, we analyse the achievable accuracy of heat‐front position estimates. In case of large groundwater variations at the field, the gravity observations can be significantly affected and, consequently, the accuracy of heat‐front monitoring can be deteriorated. We show that, by applying gravity corrections based on local observations of groundwater, the hydrological influences can to a large extent be reduced and the accuracy of estimates can be improved. We conclude that gravimetric monitoring of the heat‐front evolution has a great potential.     

Non-stationary covariance function modelling in 2D least-squares collocation

Standard least-squares collocation (LSC) assumes 2D stationarity and 3D isotropy, and relies on a covariance function to account for spatial dependence in the observed data. However, the assumption that the spatial dependence is constant throughout the region of interest may sometimes be violated. Assuming a stationary covariance structure can result in over-smoothing of, e.g., the gravity field in mountains and under-smoothing in great plains. We introduce the kernel convolution method from spatial statistics for non-stationary covariance structures, and demonstrate its advantage for dealing with non-stationarity in geodetic data. We then compared stationary and non- stationary covariance functions in 2D LSC to the empirical example of gravity anomaly interpolation near the Darling Fault, Western Australia, where the field is anisotropic and non-stationary. The results with non-stationary covariance functions are better than standard LSC in terms of formal errors and cross-validation against data not used in the interpolation, demonstrating that the use of non-stationary covariance functions can improve upon standard (stationary) LSC.     

2013年芦山Ms7.0地震前重力场源变化分析

陆地重力观测获取的时变重力信号是研究地壳内部孕震区介质变化的重要依据。利用贝叶斯平差方法处理了川西地区2010—2012年的流动重力观测资料,发现相对重力仪的格值系数在两年内变化可达2×10^-4,而且仪器存在较小的非线性漂移率;重力场累积动态变化图像显示,2013年芦山MS7.0地震发生在重力梯度带零值线与龙门山断裂带交会部位附近。进一步采用等效源方法划定了测网内以雅安、康定和石棉为中心的场源分辨能力较高的区域,并反演了地震前的地壳视密度变化。场源视密度动态变化结果表明:2013年芦山MS7.0地震前2年,在龙门山断裂带的南侧和鲜水河断裂南段的东侧区域形成了范围较大的等效视密度增加区,而在小金至康定一带和都江堰周边地区也存在小规模的视密度增加。本文基于陆地时变重力测量结果反演得到芦山地震前的地壳内部等效视密度变化,可为研究壳内深部孕震环境变化特征提供参考。     

Analysis of Heat Budget, Dissipation and Frontogenesis in a Shallow Density Current

High frequency wind and temperature measurements, obtained in March 1995 from a 10-m tower array situated in south-east Kansas, are used to analyze the structure of a shallow density current. This current is approximately 7 m deep and exhibits a current head that is estimated to be about twice the current depth. The event lasted approximately 900 s and its origin appears to be a shallow slope 2–:5 km to the west of the site, where cold air drainage occurs. The onset of the event is marked by a 5 °C temperature decrease at the 3-m level, increased variance of temperature and of wind velocity, and increased dissipation of kinetic energy, measured by a hot-wire anemometer at the 3-m level. The primary contributors to temperature changes following the frontal passage are both horizontal and vertical advections; contributions from flux divergences of temperature and of radiation, and from dew formation, do not appear to be significant. Postulated frontogenesis, prior to the arrival of the apparent equilibrated front of approximately 176-m width at the site, is examined by means of a theoretical model. The time required to equilibrate the front, by means of kinetic energy dissipation within the frontal zone, is determined to be less than 300 s, or less than the estimated travel time from the orographic slope to the observational site. The absence of upstream data is determined, however, to be a limitation of the analysis performed.     

区域重力地改“测点高程平移法”精度探讨

在讨论测点高程平移法区域重力地改的精度问题基础上，提出了系统差概念，并推荐一种新的计算方法，即把测点高程平移法和节点高程插值法有机地结合起来实现地改，通过理论模型和实例计算，表明该方法计算精度优于节点高程插值法和测点高程平移法。     

太行山东麓焚风天气的统计特征和机理分析Ⅱ:背风波对焚风产生和传播影响的个例分析

利用中尺度模式WRF3.3对太行山东麓焚风典型个例进行了数值模拟。结果表明,太行山东麓焚风的发生和移动与山脉背风波密切相关。由此建立了太行山东麓焚风的概念模型:西北或偏西气流途经山西盆地、山西境内的山脉或高原,再越过太行山,在其东麓形成背风波。背风波的下沉气流气温按干绝热方式上升,同时下沉气流也会对低层大气产生压缩增温效应,使得太行山东麓产生焚风。背风波即为重力波,可以伴随着下沉气流向下游移动,正变温区同时也向东移动。变温区移动的速度和重力波的传播速度相同。背风波的产生,需要Scorer数向上足够的减小,而且不连续,即要求大气是稳定的且存在明显的风速切变。