不同模式模拟的降水稳定同位素效应的比较和评估

利用引入水稳定同位素循环的ECHAM4、GISS E、HadCM3、MUGCM以及iAWBM的模拟数据,分析了全球降水中稳定同位素效应的空间分布特征,对不同模式的模拟结果之间以及模拟结果与GNIP（Global Network of Isotopes in Precipitation）的实际监测结果之间进行了比较,旨在对稳定同位素大气环流模式的模拟有效性进行评价,改善对水循环中水稳定同位素效应的理解和认识。结果显示,5个模式均很好地再现了全球降水中平均δ18O和平均δ18O季节差的空间分布特征,降水中稳定同位素的温度效应、降水量效应的分布特点以及全球大气水线GMWL（Global Meteoric Water Line）均被很好地模拟出。比较而言,ECHAM4模拟的降水中的平均δ18O以及δ18O平均季节差的空间分布与GNIP的实际分布最接近,拟合水平也最高;ECHAM4、GISS E、MUGCM和iAWBM再现全球温度效应空间分布的能力较强,拟合水平大致相当;由iAWBM模拟的降水量效应空间分布与实际分布之间的相关性最强,5个模式模拟的与实测的δ18O/P相关系数符号相同的站点数大致位于同一水平;GISS E和iAWBM模拟的全球大气水线与实测的GMWL最接近。     

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.     

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.     

帕米尔东北缘地区构造变形特征与盆山结构

帕米尔东北缘乌泊尔地区是正确认识帕米尔北缘盆山结构和构造变形特征非常关键的地区,本文利用连续电磁剖面(CEMP)资料和地震资料,并结合野外地质调查资料和钻井资料,对帕米尔东北缘乌泊尔地区的盆山结构和构造变形特征进行了研究。认为帕米尔东北缘及其以北地区的盆山结构表现为帕米尔造山带向北冲断和南天山向南冲断所形成的对冲结构; 帕米尔山前为基底卷入式构造,古生界—中生界沿高角度的逆冲断层推覆到新近系和第四系之上,形成山前的古生界—中生界逆冲推覆带; 北侧由受乌泊尔断裂控制的深部隐伏冲断体系和浅部的第四纪背驮盆地所构成。研究区的新生代构造变形时间开始于上新世晚期,并持续变形至今,形成了下更新统西域组（Q1x）与下伏上新统、Q2与Q1和Q3 4与Q2之间的不整合。研究区最小构造缩短量为486 km,缩短率为481％。     

An iterative scheme for use in MT 2D forward modelling — to speed up production of the Jacobian and to improve accuracy of the inverse matrix

An important stage in two-dimensional magnetotelluric modelling is the calculation of the Earth's response functions for an assumed conductivity model and the calculation of the associated Jacobian relating those response functions to the model parameters. The efficiency of the calculation of the Jacobian will affect the efficiency of the inversion modelling. Rodi (1976) produced all the Jacobian elements by inverting a single matrix and using an approximate first-order algorithm. Since only one inverse matrix required calculation the procedure speeded up the inversion. An iterative scheme to improve the approximation to the Jacobian information is presented in this paper. While this scheme takes a little longer than Rodi's algorithm, it enables a more accurate determination of the Jacobian information. It is found that the Jacobian elements can be produced in 10% of the time required to calculate an inverse matrix or to calculate a 2D starting model. A modification of the algorithm can further be used to improve the accuracy of the original inverse matrix calculated in a 2D finite difference program and hence the solution this program produces. The convergence of the iteration scheme is found to be related both to the originally calculated inverse matrix and to the change in the newly formed matrix arising from perturbation of the model parameter. A ridge regression inverse algorithm is used in conjunction with the iterative scheme for forward modelling described in this paper to produce a 2D conductivity section from field data.     

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.