2D thermo-mechanical-chemical coupled numerical models of interactions between a cooling magma chamber and a visco-elastic host rock

Author(s):  
Dániel Kiss ◽  
Evangelos Moulas ◽  
Lisa Rummel ◽  
Boris Kaus

<div> <p>A recent focus of studies in geodynamic modeling and magmatic petrology is to understand the coupled behavior between deformation and magmatic processes. Here, we present a 2D numerical model of an upper crustal magma (or mush) chamber in a visco-elastic host rock, with coupled thermal, mechanical and chemical processes, accounting for thermodynamically consistent material parameters. The magma chamber is isolated from deeper sources of magma (at least periodically) and it is cooling, and thus shrinking. We quantify the changes of pressure and stress around a cooling magma chamber and a warming host rock, using a compressible visco-elastic formulation, considering both simplified idealized and more complex and realistic geometries of the magma chamber.</p> </div><div> <p>We present solutions based on a self-consistent system of the conservation equations for coupled thermo-mechanical-chemical processes, under the assumptions of slow (negligible inertial forces), visco-elastic deformation and constant chemical bulk composition. The thermodynamic melting/crystallization model is based on a pelitic melting model calculated with Perple_X, assuming a granitic composition and is incorporated as a look-up table. We will discuss the numerical implementation, show the results of systematic numerical simulations, and illustrate the effect of volume changes due to temperature changes (including the possibility melting and crystallization) on stress and pressure evolution in magmatic systems.</p> </div>

2020 ◽  
Author(s):  
Dániel Kiss ◽  
Evangelos Moulas ◽  
Lisa Rummel ◽  
Boris Kaus

<p>A recent focus of studies in geodynamic modeling and magmatic petrology is to understand the coupled behavior between deformation and magmatic processes. Here, we present a 2D numerical model of an upper crustal magma (or mush) chamber in a visco-elastic host rock, with coupled thermal, mechanical and chemical (TMC) processes. The magma chamber is isolated from deeper sources of magma and it is cooling, and thus shrinking. We quantify the mechanical interaction between the shrinking magma chamber and the surrounding host rock, using a compressible visco-elastic formulation, considering several geometries of the magma chamber.</p><p>We present a self-consistent system of the conservation equations for coupled TMC processes, under the assumptions of slow (negligible inertial forces), visco-elastic deformation and constant chemical bulk composition. The thermodynamic melting/crystallization model is based on a pelitic melting model calculated with Perple_X, assuming a granitic composition and is incorporated as a look-up table. We will discuss the numerical implementation, show the results of systematic numerical simulations, and illustrate the effect of volume changes due to crystallization on stresses in the host rocks.</p>


1989 ◽  
Vol 126 (6) ◽  
pp. 659-673 ◽  
Author(s):  
Håvard Gautneb ◽  
Agust Gudmundsson ◽  
Niels Oskarsson

AbstractStrike, dip, and thickness were measured for 504 sheets (inclined sheets and dykes) in the 4–6 Ma old Hafnarfjall central volcano in southwest Iceland. The average dip of sheets is 65°, 80% are less than 1.2 m thick, and the thickness tends to decrease with decreasing dip. In 0.5 km long traverses perpendicular to the average strike of sheets, the percentage of sheets ranges from about 6 to 11.Of 140 chemically analysed sheets most are quartz-tholeiites; a few are intermediate or acid. The sheets are chemically more evolved than the host rock and were generated by a shallow crustal magma chamber at a mature stage of the central volcano, whereas the host rock was generated earlier before the chamber was established. Trace element results suggest that the sheet magmas evolved by low-pressure fractional crystallization as well as by mixing of primitive magmas and crustal melts.A model is proposed where most of the sheets are generated by a growing shallow magma chamber. As the chamber grows its shape changes, and so does the local stress field associated with it. Because the sheets follow the stress trajectories of the local stress field, the potential pathways of the sheets change with the growth of the chamber, which may explain the common occurrence of cross-cutting sheets. From the evolved chemistry of the sheets, as well as from the pattern of the stress trajectories, it is concluded that the bulk of the sheets were injected from the upper part of the shallow magma chamber.


1975 ◽  
Vol 12 (5) ◽  
pp. 721-730 ◽  
Author(s):  
Giorgio Rivalenti

In the Fiskenaesset region (West Greenland), there are three generations of postorogenic doleritic dikes of tholeiitic affinity. Two types of differentiation are evident: (a) laterally from the contacts to center and vertically, with the upper centres of the youngest generation of dikes attaining an andesitic or rhyolitic composition; and (b) between the different generations of dikes.Major and trace element geochemistry and calculations of the cumulus composition indicate that the differentiation within dikes is due not to flow, but to a shallow crustal fractionation of an olivine tholeiite magma. The differentiation between the various generations is attributed to fractionation of an olivine tholeiite magma during its upward displacement from a deep crustal magma chamber.


2021 ◽  
Author(s):  
Matteo Desiderio ◽  
Anna J. P. Gülcher ◽  
Maxim D. Ballmer

<p>According to geochemical and geophysical observations, Earth's lower mantle appears to be strikingly heterogeneous in composition. An accurate interpretation of these findings is critical to constrain Earth's bulk composition and long-term evolution. To this end, two main models have gained traction, each reflecting a different style of chemical heterogeneity preservation: the 'marble cake' and 'plum pudding' mantle. In the former, heterogeneity is preserved in the form of narrow streaks of recycled oceanic lithosphere, stretched and stirred throughout the mantle by convection. In the latter, domains of intrinsically strong, primordial material (enriched in the lower-mantle mineral bridgmanite) may resist convective entrainment and survive as coherent blobs in the mid mantle. Microscopic scale processes certainly affect macroscopic properties of mantle materials and thus reverberate on large-scale mantle dynamics. A cross-disciplinary effort is therefore needed to constrain present-day Earth structure, yet countless variables remain to be explored. Among previous geodynamic studies, for instance, only few have attempted to address how the viscosity and density of recycled and primordial materials affect their mutual mixing and interaction in the mantle.</p><p>Here, we apply the finite-volume code <strong>STAGYY</strong> to model thermochemical convection of the mantle in a 2D spherical-annulus geometry. All models are initialized with a lower, primordial layer and an upper, pyrolitic layer (i.e., a mechanical mixture of basalt and harzburgite), as is motivated by magma-ocean solidification studies. We explore the effects of material properties on the style of mantle convection and heterogeneity preservation. These parameters include (i) the intrinsic strength of basalt (viscosity), (ii) the intrinsic density of basalt, and (iii) the intrinsic strength of the primordial material.</p><p>Our preliminary models predict a range of different mantle mixing styles. A 'marble cake'-like regime is observed for low-viscosity primordial material (~30 times weaker than the ambient mantle), with recycled oceanic lithosphere preserved as streaks and thermochemical piles accumulating near the core-mantle boundary. Conversely, 'plum pudding' primordial blobs are also preserved when the primordial material is relatively strong, in addition to the 'marble cake' heterogeneities mentioned above. Most notably, however, the rheology and the density anomaly of basalt affect the appearance of both recycled and primordial heterogeneities. In particular, they control the stability, size and geometry of thermochemical piles, the enhancement of basaltic streaks in the mantle transition zone, and they influence the style of primordial material preservation. These results indicate the important control that the physical properties of mantle constituents exert on the style of mantle convection and mixing over geologic time. Our numerical models offer fresh insights into these processes and may advance our understanding of the composition and structure of Earth's lower mantle.</p>


2021 ◽  
Author(s):  
Vanessa Montoya ◽  
Jaime Garibay-Rodriguez ◽  
Olaf Kolditz

<p>By 2080, Germany will have to store around 600 000 m<sup>3</sup> of low and intermediate-level nuclear waste (L-ILW) with negligible heat generation. This kind of waste is largely made up of used parts of nuclear power stations such as pumps, pipelines, filters, etc. placed in various types of waste containers made from either steel, cast iron, or reinforced concrete in different designs and sizes (i.e. cylindrical or box shaped). It is already decided that a total of 303 000 of the 600 000 m<sup>3</sup> L-ILW will be disposed in a final storage facility in the former iron ore mine Schacht Konrad which is under construction. However, it is still not clear where the L-ILW emplaced in in the old salt mine Asse (200 000 m<sup>3</sup>) will be stored in the future. The situation is particularly critical, as the waste have to be retrieved from the instable mine shafts partially flooded with groundwater, causing strong socio-political concerns as radioactive waste could contaminate the water nearby. For this reason, the new search for a nuclear waste repository for high-level waste (HLW), started in 2017, should also consider the possibility to accommodate the waste from Asse. Obviously, this is still subject to critics as this will make finding a final repository more difficult as storing HLW and L-ILW together requires different concepts and designs for each other and, above all, much more space.</p><p>In this context, in this contribution we have defined conceptual and numerical models to assess the hydro-chemical evolution of a L-ILW disposal cell in indurated clay rocks, involving the interaction of different components/materials and the expected hydraulic and/or chemical gradients over 100 000 years. The L-ILW disposal cell leverages a multi-barrier concept buried 400 m below the surface. The multi-barrier system is comprised of the waste matrix (i.e. backfilling the waste drums), the disposal container, the mortar backfill in the emplacement tunnel (where the disposal containers are located) and the clay host rock. The dimensions and design of the emplacement tunnel (e.g. 11 × 13 m) and disposal cells represent and consider some aspects taken into account in the designs of some European countries. In addition, tunnel walls reinforced with a shotcrete liner and the Excavation Damaged Zone is considered in the concept. The model is implemented in OpenGeoSys-6, an open-source version-controlled scientific software based on Finite Element Method which is capable of handling fully coupled hydro-chemical models by coupling OpenGeoSys to iPHREEQC. First calculation results, demonstrate that the most important processes affecting the near-field chemical evolution are i) the degradation of the concrete and cementitious grouts with porewater migrating inwards from the host rock and ii) the significant quantities of reactive and non-reactive gases (i.e. hydrogen, carbon dioxide and methane) that are generated as a result of: i) the anaerobic corrosion of metals present in the waste and containers and ii) the degradation of organic compounds by microbial and chemical processes. As a first approximation, some assumptions and simplifications have been considered, probably resulting in a wort case scenario.</p>


2020 ◽  
Author(s):  
Anna J. P. Gülcher ◽  
Maxim D. Ballmer ◽  
Paul J. Tackley ◽  
Paula Koelemeijer

<p>Despite stirring by vigorous convection over billions of years, the Earth’s lower mantle appears to be chemically heterogeneous on various length scales. Constraining this heterogeneity is key for assessing Earth’s bulk composition and thermochemical evolution, but remains a scientific challenge that requires cross-disciplinary efforts. On scales below ~1 km, the concept of a “marble cake” mantle has gained wide acceptance, emphasising that recycled oceanic lithosphere, deformed into streaks of depleted and enriched compositions, makes up much of the mantle. On larger scales (10s-100s of km), compositional heterogeneity may be preserved by delayed mixing of this marble cake with either intrinsically-dense or intrinsically-strong materials. Intrinsically dense materials may accumulate as piles at the core-mantle boundary, while intrinsically viscous domains (e.g., enhanced in the strong mineral bridgmanite) may survive as “blobs” in the mid-mantle for large timescales, such as plums in the mantle “plum pudding”<sup>1,2</sup>. While many studies have explored the formation and preservation of either intrinsically-dense (recycled) or intrinsically-strong (primordial) heterogeneity, only few if any have quantified mantle dynamics in the presence of different types of heterogeneity with distinct physical properties.<span> </span></p><p>To address this objective, we use state-of-the-art 2D numerical models of global-scale mantle convection in a spherical-annulus geometry. We explore the effects of the <em>(i)</em> physical properties of primordial material (density, viscosity), <em>(ii)</em> temperature/pressure dependency of viscosity, <em>(iii)</em> lithospheric yielding strength, and <em>(iv)</em> Rayleigh number on mantle dynamics and mixing. Models predict that primordial heterogeneity is preserved in the lower mantle over >4.5 Gyr as discrete blobs for high intrinsic viscosity contrast (>30x) and otherwise for a wide range of parameters. In turn, recycled oceanic crust is preserved in the lower mantle as “marble cake” streaks or piles, particularly in models with a relatively cold and stiff mantle. Importantly, these recycled crustal heterogeneities can co-exist with primordial blobs, with piles often tending to accumulate beneath the primordial domains. This suggests that the modern mantle may be in a hybrid state between the “marble cake” and “plum pudding” styles.<span> </span></p><p>Finally, we put our model predictions in context with recent discoveries from seismology. We calculate synthetic seismic velocities from predicted temperatures and compositions, and compare these synthetics to tomography models, taking into account the limited resolution of seismic tomography. Convection models including preserved bridgmanite-enriched domains along with recycled piles have the potential of reconciling recent seismic observations of lower-mantle heterogeneity<sup>3</sup> with the geochemical record from ocean-island basalts<sup>4,5</sup>, and are therefore relevant for assessing Earth’s bulk composition and long-term evolution.<span> </span></p><p><sup>1</sup> Ballmer et al. (2017), <em>Nat. Geosci</em>., 10.1038/ngeo2898<br><sup>2</sup> Gülcher et al. (in review), <em>EPSL</em>: Variable dynamic styles of primordial heterogeneity preservation in Earth’s lower mantle <br><sup>3</sup> Waszek et al. (2018), <em>Nat. Comm., </em>10.1038/s41467-017-02709-4 <br><sup>4</sup> Hofmann (1997), <em>Nature, </em>10.1038/385219a0; <br><sup>5</sup> Mundl et al. (2017), <em>Science, </em>10.1126/science.aal4179</p>


2014 ◽  
Vol 55 (9) ◽  
pp. 1685-1717 ◽  
Author(s):  
Wendy A. Bohrson ◽  
Frank J. Spera ◽  
Mark S. Ghiorso ◽  
Guy A. Brown ◽  
Jeffrey B. Creamer ◽  
...  

Sign in / Sign up

Export Citation Format

Share Document