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.     

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

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

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

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

贵州和四川盆地云量的气候研究

贵州和四川盆地是我国云量最多的地区,素有“天无三日晴”和“蜀犬吠日”之说,盆地西缘的雅安古时还有“雅州天漏”的谚语。从全球日照百分率分布推知,川黔地区也是世界上天气最阴沉的地区之一。川黔之阴一向为国内外气象界所注意,本文根据1951—1980     

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.     

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.     

Extensional Late Paleozoic deformation on the western margin of Pangea, Patlanoaya area, Acatlán Complex, southern Mexico

New mapping in the northern part of the Paleozoic Acatlán Complex (Patlanoaya area) records several ductile shear zones and brittle faults with normal kinematics (previously thought to be thrusts). These movement zones separate a variety of units that pass structurally upwards from: (i) blueschist-eclogitic metamorphic rocks (Piaxtla Suite) and mylonitic megacrystic granites (Columpio del Diablo granite ≡ Ordovician granites elsewhere in the complex); (ii) a gently E-dipping, listric, normal shear zone with top to the east kinematic indicators that formed under upper greenschist to lower amphibolite conditions; (iii) the Middle–Late Ordovician Las Minas quartzite (upper greenschist facies psammites with minor interbedded pelites intruded by mafic dikes and a leucogranite dike from the Columpio del Diablo granite) unconformably overlain by the Otate meta-arenite (lower greenschist facies psammites and pelites): roughly temporal equivalents are the Middle–Late Ordovician Mal Paso and Ojo de Agua units (interbedded metasandstone and slate, and metapelite and mafic minor intrusions, respectively) — some of these units are intruded by the massive, 461 ± 2 Ma, Palo Liso megacrystic granite: decussate, contact metamorphic muscovite yielded a ⁴⁰Ar/³⁹Ar plateau age of 440 ± 4 Ma; (iv) a steeply-moderately, E-dipping normal fault; (v) latest Devonian–Middle Permian sedimentary rocks (Patlanoaya Group: here elevated from formation status). The upward decrease in metamorphic grade is paralleled by a decrease in the number of penetrative fabrics, which varies from (i) three in the Piaxtla Suite, through (ii) two in the Las Minas unit (E-trending sheath folds deformed by NE-trending, subhorizontal folds with top to the southeast asymmetry, both associated with a solution cleavage), (iii) one in the Otate, Mal Paso, and Ojo de Agua units (steeply SE-dipping, NE–SW plunging, open-close folds), to (iv) none in the Patlanoaya Group. ⁴⁰Ar/³⁹Ar analyses of muscovite from the earliest cleavage in the Las Minas unit yielded a plateau age of 347 ± 3 Ma and show low temperature ages of  260 Ma. Post-dating all of these structures and the Patlanoaya Group are NE-plunging, subvertical folds and kink bands. An E–W, vertical normal fault juxtaposes the low-grade rocks against the Anacahuite amphibolite that is cut by megacrystic granite sheets, both of which were deformed by two penetrative fabrics. Amphibole from this unit has yielded a ⁴⁰Ar/³⁹Ar plateau age of 299 ± 6 Ma, which records cooling through  490 °C and is probably related to a Permo-Carboniferous reheating event during exhumation. The extensional deformation is inferred to have started in the latest Devonian ( 360 Ma) during deposition of the basal Patlanoaya Group, lasting through the rapid exhumation of the Piaxtla Suite at  350–340 Ma synchronous with cleavage development in the Las Minas unit, deposition of the Patlanoaya Group with active fault-related exhumation suggested by Mississippian and Early Permian conglomerates ( 340 and 300 Ma, respectively), and continuing at least into the Middle Permian (≡ 260 Ma muscovite ages). The continuity of Mid-Continent Mississippian fauna from the USA to southern Mexico suggests that this extensional deformation occurred on the western margin of Pangea after closure of the Rheic Ocean.     

A preliminary stable isotope study on Mogok Ruby, Myanmar

The primary occurrence of ruby in the Mogok area, northern Myanmar is exclusively found in marble along with spinel–forsterite-bearing marble and phlogopite–graphite marble. These marble units are enclosed within banded biotite–garnet–sillimanite–oligoclase gneisses. Samples of these marbles collected for C–O stable isotope analysis show two trends of δ¹³C–δ¹⁸O variation resulting most likely from fluid–rock interactions. Ruby-bearing marble and phlogopite–graphite marble follow a trend with coupled C–O depletion, whereas spinel–forsterite-bearing marble follows a δ¹⁸O depletion trend with relatively constant δ¹³C values. Ruby formation might have resulted from CO₂-rich fluid–rock interaction, while spinel–forsterite-bearing marble was genetically related to CO₂-poor fluid–rock interaction. Both fluids may have arisen from external sources. Based on graphite Raman spectral thermometry, the estimated temperature for phlogopite–graphite marble, and probably ruby-bearing marble, was lower than 607 °C, and for spinel–forsterite-bearing marble, lower than 710 °C. Contrasting C/O diffusion between graphite/ruby/spinel/forsterite and calcite, local variations of isotopic compositions of newly formed minerals as a result of non-pervasive fluid infiltration, and open-system isotopic disturbance during cooling may have affected C-/O-isotopic fractionations between minerals. The estimated high formation temperatures for ruby and spinel/forsterite imply that the parental fluids may have been related to nearby igneous intrusions and/or metamorphic processes. Whether these two types of fluid were genetically related is unclear based on the present data.     

“Cretaceous black flysch” in the Pieniny Klippen Belt, West Carpathians: a case of geological misinterpretation

This is a critical assessment of the paper by Oszczypko et al. (2004: Cretaceous Research 25, 89–113), in which they tried to prove a mid-Cretaceous age for the Szlachtowa (“black flysch”) and Opaleniec Formations, in the Pieniny Klippen Belt, West Carpathians, both of which had previously been shown to be of Jurassic age. We argue that the mid-Cretaceous age assignment is a misinterpretation, primarily resulting from their field samples having been collected from some Cretaceous lithostratigraphic units, tectonically associated with the Jurassic formations, and/or from tectonic contact-breccias involving Jurassic and Cretaceous strata. In addition, we suggest that they have overlooked a number of significant palaeontological papers, published since 1962, which record the presence of in situ ammonites, aptychi, belemnites, thin-shelled bivalves (Bositra), gryphaeids, foraminifera, and ostracod assemblages, all indicating a Jurassic (mainly Aalenian), and not a Cretaceous, age for the Szlachtowa Formation, and also the in situ Jurassic (Bajocian) ammonites and thin-shelled bivalves (Bositra), Bositra-microfacies, and age-diagnostic foraminiferal assemblages of the Opaleniec Formation.Our presentation here of recently published dinocyst data from well-preserved assemblages further supports the Jurassic ages for the Szlachtowa (“black flysch”) and Opaleniec Formations.     

气相色谱法测定地下水中六六六结果的不确定度评定

依照《测量不确定度评定与表示》,对气相色谱法测定地下水中六六六(HCH)四种单体结果进行了不确定度评定。分析了测量过程中引入的不确定度来源,包括提取液体积的量取、样品提取溶液的定容体积、分析仪器的进样量、标准系列溶液的测量以及仪器重复测定等分量引入不确定度及其各参数的采集和计算方法,最后合成标准不确定度,通过乘以95%概率下的扩展因子2,获得测量结果的扩展不确定度。