Directions and possibilities of mathematical geology

This paper considers the present state of mathematical geology. Three directions are recognized: applied, theoretical, and mathematical. Applied mathematical geology includes formal use of mathematics to solve problems and computer processing of data. Success is achieved by a correspondence of mathematical methods used to the nature of geological data. This correspondence can be demonstrated by purely mathematical means. Theoretical mathematical geology uses mathematics as a language of geology; however, a number of methodological problems must be solved: formalization of initial geological concepts and creation of a strict conceptual basis, substantiation of initial principles of mathematical simulation, creation of theoretical geological models, problems of elementary and coincidence in geology, and methodological substantiations of possibilities of any mathematical model to approximate geological models. The essense and significance of these problems are considered. The main task of mathematical geology is to prove its correspondence to the nature of the geological objects studied, geological data obtained, and geological problems solvable. Finally, the main problems of mathematical geology are not so much mathematical as geological and methodological.     

岩体三维节点偶对分析及其应用

本文讨论了节点偶对分析的三维扩展。节点偶对分析是针对摩擦滑动节理单元的约束边界进行的序列矢量判定。在三维节理单元的分析中,几何约束和相应的力学判据极为复杂,采用节点偶对等效结合力方法才能迎刃而解。文中对三维转向节理模型的分析得到合理的结果,它表明本文所提出的方法可成功地应用于节理岩体的有限元分析。     

地壳演化与成矿作用——以川滇地洼系“四层楼”铜矿床序列为例

提出川滇地洼系“四层楼”铜矿床序列的形成与陆壳演化的成生联系,是与本区陆壳由前地槽—地槽—地台—地洼演化各阶段与之相匹配的成矿作用的产物.与此同时,并总结了本区“四层楼”铜矿床序列的成矿作用具有明显的继承性、新生性、旋回性及层控性四大特点和多因复成矿床的成矿模式.     

南极威尔克斯地SGR钻孔处雪冰转化过程中的密度变化特征

本文着重描述了SGR钻孔处冰盖上的积雪在密实化过程中的特征变化,并对该过程进行了分段的和全面的回归分析.结果表明,冰盖密度随深度增大,但增长幅度随深度减小.作者提出密度变化减小度的概念.计算得出的所研究冰芯钻取点的密度变化减小度为-0.15kg/m~3·m~2,粒雪成冰前的密实速率平均值为4.08kg/m~3·a.本文得到的冰盖密度变化“临界点”与以往报道的有所不同.分析这一现象时,作者强调当积雪还在活动层时冰盖温度的影响,并以此解释密度剖面的异常变化以及离差的回升.特别指出,积雪的密度变化具有气候学意义,它在一定程度上能够反映出积雪形成及变化过程中气候变化的某些信息.本文由密度变化确定的钻孔点雪冰转化深度为50米.     

Alluvial architecture of the Holocene Rhine–Meuse delta (the Netherlands)

Ancient fluvial successions often act as hydrocarbon reservoirs. Sub‐surface data on the alluvial architecture of fluvial successions are often incomplete and modelling is performed to reconstruct the stratigraphy. However, all alluvial architecture models suffer from the scarcity of field data to test and calibrate them. The purposes of this study were to quantify the alluvial architecture of the Holocene Rhine–Meuse delta (the Netherlands) and to determine spatio‐temporal trends in the architecture. Five north–south orientated cross‐sections, perpendicular to the general flow direction, were compiled for the fluvial‐dominated part of the delta. These sections were used to calculate the width/thickness ratios of fluvial sandbodies (SBW/SBT) and the proportions of channel‐belt deposits (CDP), clastic overbank deposits (ODP) and organic material (OP) in the succession. Furthermore, the connectedness ratio (CR) between channel belts was calculated for each cross‐section. Distinct spatial and temporal trends in the alluvial architecture were found. SBW/SBT ratios decrease by a factor of ca 4 in a downstream direction. CDP decreases from ca 0·7 (upstream) to ca 0·3 (downstream). OP increases from less than 0·05 in the upstream part of the delta to more than 0·25 in the downstream delta. ODP is approximately constant (0·4). CR is ca 0·25 upstream, which is approximately two times larger than in the downstream part of the delta. Furthermore, CDP in the downstream Rhine–Meuse delta increases after 3000 cal yr BP. These trends are attributed to variations in available accommodation space, floodplain geometry and channel‐belt size. For instance, channel belts tend to narrow in a downstream direction, which reduces SBW/SBT, CDP and CR. Tectonics cause local deviations in the general architectural trends. In addition, the positive correlation between avulsion frequency and the ratio of local to regional aggradation rate probably influenced alluvial architecture in the Rhine–Meuse delta. The Rhine–Meuse data set can be a great resource when developing more sophisticated models for alluvial architecture simulation, which eventually could lead to better characterizations of hydrocarbon reservoirs. To aid such usage of the Rhine–Meuse data set, constraints for relevant parameters are provided at the end of the paper.     

Detailed <Superscript>40</Superscript>Ar/<Superscript>39</Superscript>Ar dating of geologic events associated with the Mantos Blancos copper deposit,northern Chile

The ⁴⁰Ar/³⁹Ar geochronological method was applied to date magmatic and hydrothermal alteration events in the Mantos Blancos mining district in the Coastal Cordillera of northern Chile, allowing the distinction of two separate mineralization events. The Late Jurassic Mantos Blancos orebody, hosted in Jurassic volcanic rocks, is a magmatic-hydrothermal breccia-style Cu deposit. Two superimposed mineralization events have been recently proposed. The first event is accompanied by a phyllic hydrothermal alteration affecting a rhyolitic dome. The second mineralization event is related to the intrusion of bimodal stocks and sills inside the deposit. Because of the superposition of several magmatic and hydrothermal events, the obtained ⁴⁰Ar/³⁹Ar age data are complex; however, with a careful interpretation of the age spectra, it is possible to detect complex histories of successive emplacement, alteration, mineralization, and thermal resetting. The extrusion of Jurassic basic to intermediate volcanic rocks of the La Negra Formation is dated at 156.3 ± 1.4 Ma (2σ) using plagioclase from an andesitic lava flow. The first mineralization event and associated phyllic alteration affecting the rhyolitic dome occurred around 155–156 Ma. A younger bimodal intrusive event, supposed to be equivalent to the bimodal stock and sill system inside the deposit, is probably responsible for the second mineralization event dated at ca. 142 Ma. Other low-temperature alteration events have been dated on sericitized plagioclase at ca. 145–146, 125, and 101 Ma. This is the first time that two distinct mineralization events have been documented from radiometric data for a copper deposit in the metallogenic belt of the Coastal Cordillera of northern Chile. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Statistical analysis of sand grain/bed collision process recorded by high-speed digital camera

A high‐speed digital camera was employed to record the sand grain/bed collision process. With image processing and a statistical method, a series of parameters of the collision process were obtained. The results show that the collision process of a grain with rebounding can be represented by two parameters: the kinetic energy restitution coefficient and the collision angle. Both parameters satisfy a normal distribution, and they are dependent on one another. With an increase of the collision angle, the distribution of the kinetic energy restitution gradually reduces from a broad to a narrow range with low values. The percentage of vertical velocity restitution coefficients greater than 1 can reach 70% or more, which ensures that the settling time of the sand grains in the air increases and that they receive more energy from the air to progress the saltation movement.     

A seismic source zone model for the seismic hazard assessment of the Italian territory

We designed a new seismic source model for Italy to be used as an input for country-wide probabilistic seismic hazard assessment (PSHA) in the frame of the compilation of a new national reference map.

We started off by reviewing existing models available for Italy and for other European countries, then discussed the main open issues in the current practice of seismogenic zoning.

The new model, termed ZS9, is largely based on data collected in the past 10 years, including historical earthquakes and instrumental seismicity, active faults and their seismogenic potential, and seismotectonic evidence from recent earthquakes. This information allowed us to propose new interpretations for poorly understood areas where the new data are in conflict with assumptions made in designing the previous and widely used model ZS4.

ZS9 is made out of 36 zones where earthquakes with M_w > = 5 are expected. It also assumes that earthquakes with M_w up to 5 may occur anywhere outside the seismogenic zones, although the associated probability is rather low. Special care was taken to ensure that each zone sampled a large enough number of earthquakes so that we could compute reliable earthquake production rates.

Although it was drawn following criteria that are standard practice in PSHA, ZS9 is also innovative in that every zone is characterised also by its mean seismogenic depth (the depth of the crustal volume that will presumably release future earthquakes) and predominant focal mechanism (their most likely rupture mechanism). These properties were determined using instrumental data, and only in a limited number of cases we resorted to geologic constraints and expert judgment to cope with lack of data or conflicting indications. These attributes allow ZS9 to be used with more accurate regionalized depth-dependent attenuation relations, and are ultimately expected to increase significantly the reliability of seismic hazard estimates.