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.     

Geomorphological site classification for the Fujian province

Geomorphological classification of the Fujian province was done based on remote-sensing imaging, digital elevation maps, and slope gradient data acquired by ArcGIS 9.2. The engineering geological units in the Fujian province were divided into five types by the geomorphologic shape and genesis. The relationship among geomorphological type, engineering geological unit, and site category was determined using engineering geological data and site category data. Then, after making a preliminary site classification adjusted by drill hole data, the site was classified into four types: I₀, I₁, II, and III according to site classification standards of equivalent shear-wave velocity and overburden thickness. The results showed that the site categories in the Fujian province mainly consist of type II, which accounts for 85.26 % of the land area. The percentage of I₀ was the smallest, which accounts for only 2.44 % of the total area.     

Application of data assimilation for improving forecast of water levels and residual currents in Singapore regional waters

Hydrodynamic models are commonly used for predicting water levels and currents in the deep ocean, ocean margins and shelf seas. Their accuracy is typically limited by factors, such as the complexity of the coastal geometry and bathymetry, plus the uncertainty in the flow forcing (deep ocean tide, winds and pressure). In Southeast Asian waters with its strongly hydrodynamic characteristics, the lack of detailed marine observations (bathymetry and tides) for model validation is an additional factor limiting flow representation. This paper deals with the application of ensemble Kalman filter (EnKF)-based data assimilation with the purpose of improving the deterministic model forecast. The efficacy of the EnKF is analysed via a twin experiment conducted with the 2D barotropic Singapore regional model. The results show that the applied data assimilation can improve the forecasts significantly in this complex flow regime.     

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.     

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应用中所涉及的地震动选取、桥梁结构建模、土体本构选择以及基于构件的桩—土—桥梁结构体系损伤评估方法,为研究震后桥梁结构经济损失分析方法提供可行性方案。     

Distribution of suspended particulate matter in the northern Bohai Bay in summer and its relation with thermocline

Field observations of suspended particulate matter (SPM) in the Bohai Bay, China have not been widely reported. The aim of this paper is to describe the horizontal and vertical distribution of mass and volume concentrations of SPM, respectively, based on observed data at 312 stations in the northern Bohai Bay during summer of 2006. A numerical model ECOMSED coupled with a sediment transport module was also established to further discuss the mechanism of the thermocline effect on the vertical distribution of SPM. The mass concentrations of SPM exhibited high inshore values and low offshore values in the horizontal distribution; while in the vertical direction, characteristics of the volume concentration of SPM can be divided into two types: one with a sharp peak at depth of 10–15 m and another without. The peak value at the depth of the thermocline was resulted from concentrated phytoplankton. A numerical experiment further displayed that the thermocline can also prevent particles from being resuspended upward.     

Multiple magnetic clouds in interplanetary space

An interplanetary magnetic cloud (MC) is usually considered the byproduct of a coronal mass ejection (CME). Due to the frequent occurrence of CMEs, multiple magnetic clouds (multi-MCs), in which one MC catches up with another, should be a relatively common phenomenon. A simple flux rope model is used to get the primary magnetic field features of multi-MCs. Results indicate that the magnetic field configuration of multi-MCs mainly depends on the magnetic field characteristics of each member of multi-MCs. It may be entirely different in another situation. Moreover, we fit the data from the Wind spacecraft by using this model. Comparing the model with the observations, we verify the existence of multi-MCs, and propose some suggestions for further work.     

About Hypothesis of the Superconducting Origin of the Saturn’s Rings

Hypothesis of possible superconductivity of the iced matter of the rings of Saturn (based on the data of Voyager and Pioneer space missions) allow us to explain many phenomena which have not been adequately understood earlier. Introducing into planetary physics force of magnetic levitation of the superconducting iced particle of the rings, which interact with magnetosphere of the planet, becomes to be possible to explain origin, evolution, and dynamics of the rings; to show how the consequent precipitation of the rings’ matter upon the planet was concluded; how the rings began their rotation; how they were compressed by the magnetic field into the thin disc, and how this disc was fractured into hundreds of thousands of separated rings; why in the ring B do exist “spokes”; why magnetic field lines have distortion near by ring F; why there is a variable azimuth brightness of the ring A; why the rings reflected radio waves so efficiently; why exists strong electromagnetic radiation of the rings in the 20.4 kHz–40.2 MHz range and Saturnian kilometric radiation; why there is anomalous reflection of circularly polarized microwaves; why there are spectral anomalies of the thermal radiation of the rings; why the matter of the various rings does not mix but preserves its small-scale color differences; why there is an atmosphere of unknown origin nearby the rings of Saturn; why there are waves of density and bending waves within Saturn’s rings; why planetary rings in the solar system appear only after the Belt of Asteroids (and may be the Belt of Asteroids itself is a ring for the Sun); why our planet Earth has no rings of its own.