Efficient Modeling of Gravity Fields Caused by Sources with Arbitrary Geometry and Arbitrary Density Distribution

We present a brief review of gravity forward algorithms in Cartesian coordinate system, including both space-domain and Fourier-domain approaches, after which we introduce a truly general and efficient algorithm, namely the convolution-type Gauss fast Fourier transform (Conv-Gauss-FFT) algorithm, for 2D and 3D modeling of gravity potential and its derivatives due to sources with arbitrary geometry and arbitrary density distribution which are defined either by discrete or by continuous functions. The Conv-Gauss-FFT algorithm is based on the combined use of a hybrid rectangle-Gaussian grid and the fast Fourier transform (FFT) algorithm. Since the gravity forward problem in Cartesian coordinate system can be expressed as continuous convolution-type integrals, we first approximate the continuous convolution by a weighted sum of a series of shifted discrete convolutions, and then each shifted discrete convolution, which is essentially a Toeplitz system, is calculated efficiently and accurately by combining circulant embedding with the FFT algorithm. Synthetic and real model tests show that the Conv-Gauss-FFT algorithm can obtain high-precision forward results very efficiently for almost any practical model, and it works especially well for complex 3D models when gravity fields on large 3D regular grids are needed.     

Linear spin-up in a sliced cylinder

Abstract

spin-up and spin-down in a circular tank with a uniformly sloping bottom are studied experimentally and numerically for small values of the relative change in the angular velocity of the tank. Generally, the initial single-cell flow evolves into a number of smaller vortices. The evolution is compared with an analytical model based on an expansion of the flow field in linear Rossby waves (Pedlosky and Greenspan, 1967). Although it is possible to tune the experimental parameters in such a way that agreement with the theory is found, in most cases the experiments show shedding of vortices in the initial stage of the spin-up or spin-down, a phenomenon not described by the analytical model. Nonetheless, in such cases the analytical model still accounts for other observations: the alternating generation of cyclonic and anticyclonic vortices in the eastern part of the tank and their subsequent westward motion.     

Stochastic evaluation of simple pairing approaches to reconstruct incomplete rainfall time series

Two-station pairing approaches are routinely used to infill missing information in incomplete rainfall databases. We evaluated the performance of three simple methodologies to reconstruct incomplete time series in presence of variable nonlinear correlation between data pairs. Nonlinearity stems from the statistics describing the marginal peak-over-threshold (POT) values of rainfall events. A Monte Carlo analysis was developed to quantitatively assess expected errors from the use of chronological pairing (CP) with linear and nonlinear regression and frequency pairing (FP). CP is based on a priori selection of regression functions, while FP is based on matching the probability of non-exceedance of an event from one time series with the probability of non-exceedance of a similar event from another time series. We adopted a generalized Pareto (GP) model to describe POT events, and a t-copula algorithm to generate reference nonlinearly correlated pairs of random temporal distributions distributed according with the GP model. The results suggest that the optimal methodology strongly depends on GP statistics. In general, CP seems to provide the lowest errors when GP statistics were similar and correlation became linear; we found that a power-2 function performs well for the selected statistics when the number of missing points is limited. FP outperforms the other methods when POT statistics are different and variables are markedly nonlinearly correlated. Ensemble-based results seem to be supported by the analysis of observed precipitation at two real-world gauge stations.     

The spatial distribution of archaeal lipids in a mesoscale subtropical watershed,Southeast China

Archaea play an important role in global carbon and nitrogen cycles. Archaeal lipids, such as isoprenoid glycerol diakyl glycerol tetraethers (iGDGTs), are important biomarkers tracing change in archaeal community structure and biogeochemical processes in the natural environments. In this research, the spatial distributions of archaeal lipids in the surface sediments of the Jiulong River (JR) and the Jiulong River estuary (JRE) were examined. GDGT-0 (containing zero cyclopentyl ring) and crenarchaeol were the most abundant iGDGTs in the JR and JRE. From the rivers to the estuary, the total iGDGTs, GDGT-0, crenarchaeol and archaeol concentrations showed significant spatial variation; in particular, GDGT-0 and archaeol in the river may be predominantly derived in situ from methanogens, whereas crenarchaeol in the estuary mainly derived in situ from Thaumarchaeota. We inferred that archaeal community was dominated by methanogens in the Jiulong River and by Thaumarchaeota in the Jiulong River estuary, which are consistent with change in archaeal community structure observed in other estuarine environments.     

Long-term Geomechanical Observation in the Jeroným Mine

Geomechanical measurements have been carried out in the medieval Jeroným Mine since 2001. At first quarterly observation and later also continual monitoring by distributed measurement network were used. The network includes a number of different sensors that have been installed in selected places of underground spaces. There have not been any manifestations documenting dangerous changes of stability during the period of monitoring (2001–2015). The results of accomplished measurements from all convergence profiles confirm the total stability of all underground spaces. However, a number of places have been documented as potentially unstable. Measured values of selected parameters are presented in the paper, namely, movements along fractures, changes of stress tensor and fluctuations in the level of mine water table.     

5•12汶川地震中山区公路桥梁震害及启示

桥梁结构的振动特性受场地条件及土—结相互作用影响显著,为提高桥梁结构的抗震性能,减少因结构破坏而产生的经济损失,进一步开展考虑场地土体非线性以及土—结相互作用对桥梁结构地震反应影响的研究工作十分必要。本文以一典型桥梁结构为例,着重介绍了美国太平洋地震工程研究中心基于性能的桥梁有限元模拟平台BridgePBEE应用中所涉及的地震动选取、桥梁结构建模、土体本构选择以及基于构件的桩—土—桥梁结构体系损伤评估方法,为研究震后桥梁结构经济损失分析方法提供可行性方案。