COMPUTATIONAL SIMULATION OF RIVER SEDIMENTATION AND MORPHOLOGY - A REVIEW OF THE STATE OF THE ART

In this paper, the state-of-the-art computational simulation technologies for river sedimentation and morphology are comprehensively reviewed. Emphasis is on computational modelling, which includes mathematical modelling equations and closures, numerical solution methods, verification and validation, model integration, and applications. Several important areas of future research are recommended.     

地壳深部岩浆岩岩基体与大型、超大型金属矿床的形成及找矿效应

基于部分大型、超大型金属矿床和多金属矿集区的地球物理研究发现:地壳内部物质在热动作用下,深部物质和能量强烈的在交换,当深部热液沿通道（破碎带、软弱岩层或断裂系）上涌时,在深部适宜的岩层部位聚积,并形成大型或较大型的酸性、中酸性或基性的岩基体（镁铁质、超镁铁质）。这些类型的岩浆岩以其携带的热能驱动着深部物质重新分异和调整,并在这一动力作用过程中与围岩进行交代、蚀变,并促使其不断变质。在这样的深部介质和构造环境下矿物元素得以在大型或较大型岩基体附近与周边介质耦合,逐使矿物元素不断聚积,并形成大型、超大型金属矿床和多金属矿集区。     

北苏鲁仰口地区变辉长岩中锆石U-Pb定年、微量元素和Hf同位素特征及其地质意义

北苏鲁仰口地区出露超高压的变辉长岩.锆石阴极发光图像和其内部矿物包体激光拉曼测试的联合研究结果表明,变辉长岩锆石具有弱发光效应的岩浆韵律环带的核和被改造的强发光效应的边.岩浆韵律环带的核部保存大量而复杂的矿物包体,包括普通辉石(Cpx)+斜方辉石(Opx)+斜长石(P1)+石英(Qtz)+黑云母(Bt)+钛铁矿(Ilm)+磷灰石(Ap);边部保存的矿物包体则相对较少,包括普通辉石(Cpx)+斜方辉石(Opx)+斜长石(Pl)+磷灰石(Ap).尽管岩浆韵律环带核部的稀土元素总合量比被改造的锆石边部明显偏高,但二者稀土元素配分模式具有明显的相似性,主要表现为轻稀土相对亏损,而重稀土明显富集,相应的(La/Yb)N=0.00015～0.00039,并具有明显的负Eu异常(Eu/Eu*=0.20～0.26)、相对明显的正Ce异常(Ce/Ce*=71.5～147.4)和较高的Th/U比值(1.97～2.69).上述特征表明,仰口地区变辉长岩中的锆石均为继承性的岩浆锆石,而没有新生的变质锆石.LA-(MC)-ICP-MS锆石原位U-Pb定年和Lu-Hf同位素分析结果表明,两件锆石样品Y1和Y2的年龄数据所构成的不一致线显示了十分接近的上交点和下交点年龄.其上交点年龄分别为785±15Ma(2σ)和784±12Ma(2σ),应代表原岩的形成时代,表明变辉长岩的原岩与Rodinia超大陆裂解的岩浆事件存在密切的成因关系;而下交点年龄分别为226±24Ma(2σ)和228±26Ma(2σ),与苏鲁其它类型超高压岩石中含柯石英锆石微区记录的变质年龄十分吻合,应代表变辉长岩的超高压变质时代.岩浆结晶锆石的核部具有明显偏高的176Lu/177Hf(0.00044～0.00291)和176Yh/177Hf(0.0165～0.1168)比值,而176Hf/177Hf比值变化于0.281956～0.282048之间,相应的εHf(t)=-8.5～-14.0,tDM2=2.03～2.32Ga,表明仰口地区变辉长岩的原岩起源于古元古代时期的富集地幔或发生部分熔融的下地壳残留体.被改造的岩浆结晶锆石的边部则具有明显偏低的176Lu/177Hf(0.00029～0.00060)和176Yh/177Hf(0.0112～0.0200)比值,而176Hf/177Hf(t)比值变化于0.281953～0.282002之间,相应的εHf(t)=-10.2～-11.9,tDM2=2.12～2.21Ga.与岩浆结晶锆石核部相比,被改造的岩浆锆石边部的176Lu/177Hf、176Yb/177Hf、176Hf/177Hf(t)比值和εHf(t)和tDM2值的变化范围更小,表明中-新三叠纪的超高压变质作用使岩浆结晶锆石边部的Lu-Hf同位素体系发生调整,更趋向于均一化.     

Topographic expression of active faults in the foothills of the Northern Apennines

Active faults that rupture the earth's surface leave an imprint on the topography that is recognized using a combination of geomorphic and geologic metrics including triangular facets, the shape of mountain fronts, the drainage network, and incised river valleys with inset terraces. We document the presence of a network of active, high-angle extensional faults, collectively embedded in the actively shortening mountain front of the Northern Apennines, that possess unique geomorphic expressions. We measure the strain rate for these structures and find that they have a constant throw-to-length ratio. We demonstrate the necessary and sufficient conditions for triangular facet development in the footwalls of these faults and argue that rock-type exerts the strongest control. The slip rates of these faults range from 0.1 to 0.3 mm/yr, which is similar to the average rate of river incision and mountain front unroofing determined by corollary studies. The faults are a near-surface manifestation of deeper crustal processes that are actively uplifting rocks and growing topography at a rate commensurate with surface processes that are eroding the mountain front to base level.     

Consequences of crystal shape and fabric on anisotropic permeability in magmatic mush

Crystals that form an interconnected porous network can become preferentially oriented both prior to and during compaction of magmatic mush. This introduces anisotropy in the melt pore-space that can reduce permeability in the direction of compaction and in turn decrease melt flux and compaction rate. Using a number of grain-scale numerical models, the consequences of end-member magmatic fabrics on the directional dependence of permeability are tested over a range in melt fraction from 22 to 77%. As the crystal aspect ratio (i.e. ratio of long to short axis length) increases from 2 to 10, isotropic permeability decreases by a factor of 2 and 5 for randomly oriented prolate and oblate-shaped crystals, respectively, at a melt fraction of 22%. With a flattening fabric, permeability is reduced in the compaction direction no more than approximately a factor of 2 relative to the isotropic permeability at the same melt fraction and crystal shape for both oblate and triaxial prisms. However, permeability is enhanced in directions orthogonal to the compaction direction. For example, permeability is enhanced up to a factor of 11 relative to the isotropic permeability at a melt fraction of 22% for oblate prisms with a ratio of the long to short axis length of 10. Anisotropy in permeability increases as the melt fraction decreases and the crystal aspect ratio increases. Ratios of the principal permeabilities are sufficiently large based on the realistic crystal shapes tested here to warrant including anisotropic permeability into macroscale melt segregation models including those for compaction.     

Thermo-hydro-mechanical simulation of ATLAS in situ large scale test in Boom Clay

Following the need for understanding and quantifying the effect of temperature on the response of a candidate host formation for radioactive waste disposal, finite element modelling of an in situ thermal experiment has been carried out. Based on a thermo-hydro-mechanical (THM) finite element approach including a consistent thermo-plastic constitutive model, it has been possible to reproduce the THM response of a clay formation submitted to in situ thermal loading. The simulated large-scale experiment, called ATLAS was designed in the underground research facility (HADES-URF) in Mol, Belgium. After an extensive literature analysis on the thermal, hydraulic and mechanical characteristics of Boom Clay, laboratory tests were simulated to calibrate model parameters. The results of the finite element modelling of the ATLAS experiment were compared with in situ measurements and revealed the necessity to account for flow diffusion in all three directions through a 2D axisymmetric analysis. Finally, those results were interpreted in the light of elasto-thermoplasticity, which emphasizes the significant role of thermo-plastic processes in the global THM response of the clay formation.     

雷暴电过程对动力发展的影响研究

利用三维时变积云动力—电过程耦合模式对雷暴中的电过程对动力发展的影响进行了模拟计算。结果表明,随着雷暴云体发展,云内形成三极性电荷结构。电场力对雷暴发展的直接作用与动力场相比很小,电涡度10-6·s-1比动力涡度小两个量级。但电场力影响了雷暴中的微物理过程,特别是水成物粒子的降落速度,并通过对三相水成物之间转换过程的调制,改变了三相水成物的时空结构,使云中水汽总量增加41%,释放潜热增加19.4%,从而改变了流场结构,电涡度达到10-4·s-1量级,和动力涡度相当。模拟结果指出,雷暴内电过程在云体成熟阶段对雷暴发展有较大影响。     

Turbidites facies response to the morphological confinement of a foredeep (Cervarola Sandstones Formation,Miocene, northern Apennines,Italy)

The Cervarola Sandstones Formation, Aquitanian–Burdigalian in age, was deposited in an elongate, north‐west stretched foredeep basin formed in front of the growing northern Apennines orogenic wedge. As other Apennine foredeep deposits, such as the Marnoso‐arenacea Formation, the stratigraphic succession of the Cervarola Sandstones Formation records the progressive closure of the basin due to the propagation of thrust fronts towards the north‐east, i.e. towards the outer and shallower foreland ramp. This process produces a complex foredeep that is characterized by syn‐sedimentary structural highs and depocentres that strongly influence lateral and vertical turbidite facies distribution. This work describes and discusses this influence, providing a high‐resolution physical stratigraphy with ‘bed by bed’ correlations of an interval ca 1000 m thick, parallel and perpendicular to the palaeocurrents and to the main structural alignments, on an area of ca 30 km that covers the proximal portion of the Cervarola basin in the northern Apennines. The main aim is to show, for the first time ever, a detailed facies analysis of the Cervarola Sandstones Formation, based on a series of bed types that have proven fundamental to understand the morphology of the basin. The knowledge of the vertical and lateral distribution of these bed types, such as contained‐reflected and slurry (i.e. hybrid) beds, together with other important sedimentary structures, i.e. cross‐bedded bypass facies and delamination structures, is the basis for better understanding of facies processes, as well as for proposing an evolutionary model of the foredeep in relation to the syn‐sedimentary growth of the main tectonic structures. This makes the Cervarola Sandstones, like the Marnoso‐arenacea Formation, a typical example of foredeep evolution.     

Mopping up leaking carbon: A natural analog at Wadi Namaleh,Jordan

Carbon capture and sequestration (CCS) is one of the important options available for partially stemming greenhouse gas emissions from large point sources. The possibility of leaking from deep storage needs to be addressed. The Wadi Namaleh area in southern Jordan provides an interesting case study of how excess CO₂ can be trapped in the form of carbonates in the near surface, even when the local geology is not obviously conducive for such a process.Carbonate veins are formed in surface alteration zones of rhyolite host rock in this arid region. The alteration zones are limited to areas where surface soil or colluvium are present. Oxygen, deuterium and carbon isotopes of the carbonates and near-surface ground water in the area suggest that the source of carbon is deep seated CO₂, and that the carbonate precipitated in local meteoric water under ambient temperature conditions. Analysis of strontium in the carbonate, fresh rhyolite and altered host shows that the source for calcium is aeolian. Trace elements show that metal and REE mobility are constrained to the alteration zone.Thus, interaction of H₂O, CO₂ and atmospheric wet and dry deposition lead to the formation of the clayey (montmorillonite) alteration zone. This zone acts to trap seeping CO₂ and water, and thus produces conditions of progressively more efficient trapping of carbon dioxide by means of a positive feedback mechanism. Replication of these conditions in other areas will minimize CO₂ leakage from man-made CCS sites.     

Compositional variations of several Early Permian magmatic sulfide deposits in the Kalatongke district,southern Altai,western China: With genetic and exploration implications

A Permian magmatic Ni-Cu sulfide deposit cluster occurs in the Kalatongke district in the Southern Chinese Altai Orogenic Belt, western China. These deposits are associated with the mafic units of the Y1, Y2, Y3, Y9 and G21 mafic-intermediate complexes. In this paper we report the first zircon U-Pb ages for the Y3 and G21 intrusions, which are 283.3 ± 1.3 Ma and 281.1 ± 1.5 Ma, respectively. Our new age data confirm that the sulfide-bearing mafic units of the Y1, Y2 (connected with Y1 at depth), Y3, Y9 and G21 intrusions all formed in Early Permian between ∼281 and ∼287 Ma. New and existing petrological-geochemical data show some important regular variations between these deposits. The host lithologies change from olivine-bearing rocks for the Y1-Y2-Y9 deposits to olivine-free rocks such as norite for the Y3 deposit and leucogabbro for the G21 deposit. The olivine Fo contents of the Y1 deposit are up to 82 mol%, which are slightly higher than those of the Y2 deposit (up to 81 mol%) and the Y9 deposit (up to 79 mol%). The average plagioclase An contents of the olivine-bearing Y1-Y2-Y9 deposits are higher than those of the olivine-free Y3-G21 deposits. Among the three deposits (Y1, Y2 and Y3) that occur closely along the same structural lineament, the Ni/Cu ratios of bulk sulfides decrease from the olivine-bearing deposits (Y1 and Y2) to the olivine-free deposit (Y3). The PGE tenors of these deposits (Y1, Y2 and Y3) and the nearby coeval deposits (Y9 and G21) are extremely low, indicating that their parental magmas are severely depleted in PGEs. The variations of PGE tenors within a single deposit as well as among the different deposits are mainly due to variable R factors. The host rocks of these deposits are all characterized by elevated initial ⁸⁷Sr/⁸⁶Sr ratios from 0.7045 to 0.7047, positive εNd values from 4.95 to 6.86, positive εHf values of zircon from 9 to 16, and elevated δ¹⁸O values of zircon from 6.15 to 6.7‰. The isotope data indicate that the parental magmas for these deposits experienced up to ∼15 wt% crustal contamination. The δ³⁴S values of the sulfide minerals from these deposits are from −3.1‰ to 0.4‰, with a peak at −2.2‰, indicating the involvement of crustal sulfur. The isotope data and mineral chemistry together indicate that both olivine fractional crystallization and addition of crustal sulfur played a role in triggering sulfide saturation in the parental magmas for these deposits. Based on higher Ni/Cu ratios of sulfide mineralization in the olivine-bearing intrusions (Y1, Y2, Y9) than in the coeval olivine-free intrusions (Y3, G21), we recommend that Ni exploration in the region focus on the olivine-bearing intrusions that were emplaced in the Early Permian.