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.     

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

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

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

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

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.     

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.     

Simultaneous transport of water and solutes under transient unsaturated flow conditions — A case study

The imbalance between incoming and outgoing salt causes salinization of soils and sub-soils that result in increasing the salinity of stream-flows and agriculture land. This salinization is a serious environmental hazard particularly in semi-arid and arid lands. In order to estimate the magnitude of the hazard posed by salinity, it is important to understand and identify the processes that control salt movement from the soil surface through the root zone to the ground water and stream flows. In the present study, Malaprabha sub-basin (up to dam site) has been selected which has two distinct climatic zones, sub-humid (upstream of Khanapur) and semi-arid region (downstream of Khanapur). In the upstream, both surface and ground waters are used for irrigation, whereas in the downstream mostly groundwater is used. Both soils and ground waters are more saline in downstream parts of the study area. In this study we characterized the soil salinity and groundwater quality in both areas. An attempt is also made to model the distribution of potassium concentration in the soil profile in response to varying irrigation conditions using the SWIM (Soil-Water Infiltration and Movement) model. Fair agreement was obtained between predicted and measured results indicating the applicability of the model.     

西藏驱龙超大型斑岩铜矿床:地质、蚀变与成矿

驱龙超大型矿床是一个产于后碰撞伸展环境下、与大洋俯冲无关的新型斑岩铜矿。文章通过对驱龙铜矿床地质、蚀变与矿化的详细研究,建立了驱龙中新世岩浆演化序列,初步查明了岩浆浅成侵位的构造控制要素,厘定了主要的围岩蚀变类型及空间展布规律,查明了引起各期蚀变事件的地质记录及矿化的空间分布规律,并探讨了成矿物质沉淀的机制,初步建立了该矿床的成矿模型。研究表明,驱龙铜矿中新世斑岩是闪长质深部岩浆房不断演化的产物,花岗闪长岩中新发现的、结晶时间为22.2Ma左右的闪长质包体可近似代表深部岩浆房组分,依次产出的花岗闪长岩、呈岩株或岩枝产出的P斑岩、X斑岩及最晚期的闪长玢岩(15.7±0.2)Ma,均为深部岩浆房连续演化的产物,岩浆持续6Ma左右。岩浆演化过程中角闪石、斜长石不断的结晶分异,导致了岩石常量元素、稀土元素及微量元素组成的规律性变化,斑岩埃达克质的特征也因岩浆演化过程中角闪石等矿物的不断结晶分异而引起。X斑岩中锆石的Hf同位素特征表明,岩石可能形成于新生下地壳的部分熔融。大面积产出的花岗闪长岩为驱龙铜矿最主要的含矿围岩,容纳了驱龙矿床70%以上的矿体,主要由斜长石、钾长石和石英组成,具花岗结构-似斑状结构,近EW向产出,其浅成就位可能受背斜控制,其后的各期斑岩均沿该侵位中心上侵,而冈底斯地壳中新世的快速抬升与剥蚀是导致含矿斑岩浅成侵位的根本原因;矿区内的SN向裂隙带既不控岩,也不控矿。浅成侵位的斑岩及深部岩浆房均发生了流体出溶。发生了大量流体出溶的深部岩浆房,是矿区早期蚀变流体的主要来源,显微晶洞构造及单向固结结构(UST)是流体出溶的地质记录。蚀变主要有3种类型,分别为早期的钾硅酸盐化、青磐岩化以及晚期的长石分解。钾硅酸盐化可分为2个阶段,即蚀变矿物以次生钾长石为主的早期钾硅酸盐化和以次生黑云母为主的晚期钾硅酸盐化。青磐岩化因产出的岩石类型不同,蚀变矿物组合具有明显差异性:产于叶巴组地层中的青磐岩化相对较强,蚀变矿物以绿帘石为主;产于花岗闪长岩中的青磐岩化相对较弱,蚀变矿物以绿泥石为主。晚期长石分解蚀变以破坏长石类矿物为特征,蚀变矿物主要为绢云母-绿泥石-粘土等。石英和硬石膏贯穿于上述各种蚀变中。空间上,钾硅酸盐化位于斑岩体及其周围地区,青磐岩化位于钾硅酸岩化外侧。后期形成的长石分解蚀变强烈叠加了早期钾硅酸盐化,介于钾硅酸盐化带与青磐岩化带之间。与早期钾长石化有关的脉体主要为不规则石英-钾长石脉,与晚期黑云母化有关的脉体主要为不规则至板状的石英-硬石膏脉、黑云母脉,与青磐岩化有关的脉体主要为板状的绿帘石-石英脉,与晚期长石分解蚀变有关的脉体主要为板状黄铜矿-黄铁矿脉及黄铁矿脉;在早期钾硅酸盐蚀变与晚期长石分解蚀变转换阶段,发育一组板状的石英-硫化物脉。早期不规则的脉体形成于斑岩结晶早期、矿区裂隙小规模发育阶段;晚期的板状脉体形成于斑岩弱固结或固结之后、矿区大规模连通裂隙发育阶段。驱龙矿区的铜矿化分布较为均一,主体产于花岗闪长岩中,其中,铜矿化主体形成于黑云母化蚀变阶段,转变阶段及长石分解阶段也有大量铜的形成;钼主要形成于转换阶段,长石分解蚀变阶段也有产出。黑云母化阶段,铜的沉淀与角闪石黑云母化、斜长石钾长石化过程中Ca2 的大量释放有关;转换阶段,铜钼矿化可能与压力和(或)温度骤降有关;晚期铜矿化与长石矿化蚀变阶段,斜长石绿泥石化、黑云母绿帘石化过程中Ca2 及Fe2 的释放有关。     

内蒙古敖包吐萤石矿床的Sr、Nd、Pb同位素地球化学特征

敖包吐萤石矿床是内蒙古北部苏莫查干地区单一萤石矿集区中的一个代表性矿床,产于早二叠世大石寨组火山-沉积岩与早白垩世敖包吐花岗岩的接触带上。文章通过分析该矿床岩、矿石的微量元素和稀土元素,揭示出萤石的成矿作用可分为2个阶段,即交代作用和充填作用。交代作用过程中大石寨组的结晶灰岩可能为萤石的形成提供了部分Ca来源,萤石矿石的稀土元素配分模式与海水基本类似,具有Ce负异常;成矿作用后期主要表现为充填作用,形成颗粒粗大的萤石,表现为重稀土元素富集的特征,并随着萤石的沉淀析出,稀土元素总量逐渐下降,反映出成矿流体经历了较长期的演化过程。各地层单元、花岗岩体和萤石矿石的Sr、Nd、Pb同位素研究表明,萤石的放射性同位素组成具有壳、幔源混合的特点,成矿物质来源具有多源性。早白垩世敖包吐花岗岩可能是萤石中F的主要来源,而大石寨组的结晶灰岩则可能提供了Ca。另外,Pb、Nd同位素的极大不均一性,有可能是成矿流体在运移过程中对艾力格庙群放射性组分的选择性吸收的结果。萤石成矿作用与钾玄岩的时空关系暗示了萤石的成矿过程可能是中国东部岩石圈减薄和下地壳的置换地质事件的结果。在构造转型的过程中,燕山中期富碱的酸性花岗岩浆的活动分异出富含F的成矿流体,与幔源流体混合,沿区域重新活化的深大断裂和大石寨组的层间破碎带上升,交代其间的灰岩透镜体,从而形成敖包吐中型萤石矿床。     

秦岭造山带的印支运动及印支期成矿作用

秦岭碰撞造山经历了长期的板块构造的俯冲-碰撞的构造演化,于印支期最终完成对接拼合,形成了统一的中国大陆,并由此转入陆内变形。众多Au、Mo多金属矿床的同位素年龄资料表明,印支期是秦岭的重要成矿期,其成矿作用明显受到构造演化的控制,反映特定的地球动力学背景和作用过程。秦岭印支期成矿作用不仅是中国东部中生代成矿作用的先导和开始,奠定了中国东部中生代成矿大爆发的基础,而且为碰撞期和碰撞期后构造体制快速转换的研究提供了依据。重视秦岭以及中国印支期成矿作用的研究,对正确认识秦岭成矿带的区域成矿规律、造山带演化的深部动力学过程,建立符合中国和东亚实际的印支期成矿理论体系具有重要科学意义。