Numerical simulation of a stratigraphic model

A multi-lithology diffusive stratigraphic model is considered, which simulates at large scales in space and time the infill of sedimentary basins governed by the interaction between tectonics displacements, eustatic variations, sediment supply, and sediment transport laws. The model accounts for the mass conservation of each sediment lithology resulting in a mixed parabolic, hyperbolic system of partial differential equations (PDEs) for the lithology concentrations and the sediment thickness. It also takes into account a limit on the rock alteration velocity modeled as a unilaterality constraint. To obtain a robust, fast, and accurate simulation, fully and semi-implicit finite volume discre tization schemes are derived for which the existence of stable solutions is proved. Then, the set of nonlinear equations is solved using a Newton algorithm adapted to the unilaterality constraint, and preconditioning strategies are defined for the solution of the linear system at each Newton iteration. They are based on an algebraic approximate decoupling of the sediment thickness and the concentration variables as well as on a proper preconditioning of each variable. These algorithms are studied and compared in terms of robustness, scalability, and efficiency on two real basin test cases.     

New age constraints on emplacement of the Cévenol granitoids, South French Massif Central

During the development of the Variscan orogeny, large amounts of granitic melt were produced, giving rise to the intrusion of granitoids at different structural levels. Despite numerous studies, ages available from previous work on the Cévennes granites remain largely imprecise. In order to better constrain the age and emplacement mode of these granites, we have combined U–Pb dating on monazites and zircons and ⁴⁰Ar/³⁹Ar dating on biotites with petrological observations, major element chemical analysis and SEM zircon imaging on five samples from the Aigoual–St Guiral–Liron and Mont Lozère granitic massifs. The results revealed that granitic intrusions and cooling in Southern Cévennes occurred in a short time span at ∼306 Ma after the main episode of regional metamorphism. Petrological and chemical data suggest that they result from a mixing between mantle-derived basic magmas (lamprophyres) and lower crust acid magmas. At a regional scale the production of these melts occurred at the end of crustal thickening induced by nappe stacking, at the same time as the late anatectic events recorded further north in the Velay dome and the granulite facies metamorphism recorded in metasedimentary granulite enclaves brought up by Tertiary volcanoes of the Velay area (Bournac).     

Continental rift systems and anorogenic magmatism

Precambrian Laurentia and Mesozoic Gondwana both rifted along geometric patterns that closely approximate truncated-icosahedral tessellations of the lithosphere. These large-scale, quasi-hexagonal rift patterns manifest a least-work configuration. For both Laurentia and Gondwana, continental rifting coincided with drift stagnation, and may have been driven by lithospheric extension above an insulated and thermally expanded mantle. Anorogenic magmatism, including flood basalts, dike swarms, anorthosite massifs and granite-rhyolite provinces, originated along the Laurentian and Gondwanan rift tessellations. Long-lived volcanic regions of the Atlantic and Indian Oceans, sometimes called hotspots, originated near triple junctions of the Gondwanan tessellation as the supercontinent broke apart. We suggest that some anorogenic magmatism results from decompression melting of asthenosphere beneath opening fractures, rather than from random impingement of hypothetical deep-mantle plumes.     

Transport in a 2-D saturated porous medium: A new method for particle tracking

A new method for solving the transport equation based on the management of a large numbe of particles in a discretized 2-D domain is presented. The method uses numerical variables to represent the number of particles in a given mesh and is more complex than the 1-D problem. The first part of the paper focuses on the specific management of particles in a 2-D problem. The method also would be valid for three dimensions as long as the medium can be modeled similar to a layered system. As the particles are no longer tracked individually, the algorithm is fast and does not depend on the number of particles present. The numerical tests show that the method is nearly numerical dispersion free and permits accurate calculations even for simulations of low-concentration transport. Because each mesh is considered as a closed system between two successive time steps, it is easy to add adsorption phenomenon without any problem of numerical stability. The model is tested under conditions that are extremely demanding for its operating mode and gives a good fit to analytical solutions. The conditions in which it can be used to best advantage are discussed.     

Modeling coastal tsunami hazard from submarine mass failures: effect of slide rheology,experimental validation,and case studies off the US East Coast

Natural Hazards - We perform numerical simulations to assess how coastal tsunami hazard from submarine mass failures (SMFs) is affected by slide kinematics and rheology. Two types of two-layer SMF...     

Cooling history of nested plutons from the Variscan Tichka plutonic complex (Morocco)

International Journal of Earth Sciences - Four imbricated mafic to felsic plutons of Variscan age from Morocco have been investigated for their cooling history and geochemical interactions with...     

The Banhad?o Alkaline Complex, Southeastern Brazil: source and evolution of potassic SiO2-undersaturated high-Ca and low-Ca magmatic series

The Cretaceous Banhad?o alkaline complex in southeastern Brazil presents two potassic SiO₂-undersaturated series. The high-Ca magmatic series consist of initially fractionated olivine (Fo_92-91) + diopside (Wo_48-43En_49-35Ae_0-7), as evidenced by the presence of xenocrysts and xenoliths. In that sequence, diopside (Wo_47-38En_46-37Ae_0-8) + phlogopite + apatite + perovskite (Prv_>92) crystallized to form the phlogopite melteigite and led to the Ca enrichment of the magma. Diopside (Wo_47-41En_32-24 Ae_3-14) continued to crystallize as an early mafic mineral, followed by nepheline (Ne_74.8-70.1Ks_26.3-21.2Qz_7.6-0.9) and leucite (Lc_65-56) and subsequently by melanite and potassic feldspar (Or_85-99Ab_1-7) to form melanite ijolites, wollastonite-melanite urtites and melanite-nepheline syenites. Melanite-pseudoleucite-nepheline syenites are interpreted to be a leucite accumulation. Melanite nephelinite dykes are believed to represent some of the magmatic differentiation steps. The low-Ca magmatic series is representative of a typical fractionation of aegirine-augite (Wo_36-29En_25-4Ae_39-18) + alkali feldspar (Or_57-96Ab_3-43) + nepheline (Ne_76.5-69.0Ks_19.9-14.4Qz_15.1-7.7) + titanite from phonolite magma. The evolution of this series from potassic nepheline syenites to sodic sodalite syenites and sodalitolites is attributed to an extensive fractionation of potassic feldspar, which led to an increase of the NaCl activity in the melt during the final stages forming sodalite-rich rocks. Phonolite dykes followed a similar evolutionary process and also registered some crustal assimilation. The mesocratic nepheline syenites showed interactions with phlogopite melteigites, such as compatible trace element enrichments and the presence of diopside xenocrysts, which were interpreted to be due to a mixing/mingling process of phonolite and nephelinite magmas. The geochemical data show higher TiO₂ and P₂O₅ contents and lower SiO₂ contents for the high-Ca series and different LILE evolution trends and REE chondrite-normalized patterns as compared to the low-Ca series. The ⁸⁷Sr/⁸⁶Sr, ¹⁴³Nd/¹⁴⁴Nd, ²⁰⁶Pb/²⁰⁴Pb and ²⁰⁸Pb/²⁰⁴Pb initial ratios for the high-Ca series (0.70407–0.70526, 0.51242–0.51251, 17.782–19.266 and 38.051–39.521, respectively) were slightly different from those of the low-Ca series (0.70542–0.70583, 0.51232–0.51240, 17.758–17.772 and 38.021–38.061, respectively). For both series, a CO₂-rich potassic metasomatized lithospheric mantle enriched the source with rutile-bearing phlogopite clinopyroxenite veins. Kamafugite-like parental magma is attributed to the high-Ca series with major contributions from the melting of the veins. Potassic nephelinite-like parental magma is assigned to the low-Ca series, where the metasomatized wall-rock played a more significant role in the melting process.     

Magnetic Characterization of Amphibolites from the Fomopéa Pluton (West Cameroon): Their Implication in the Pan-African Deformation of the Central African Fold Belt

The Fomopea granitic pluton is emplaced in gnessic and amphibolitic basement.These gneissic and amphibolitic basement rocks are represented in the pluton's body as sub-rounded,elongated or stretched xe...     

Numerical analysis of frost effects in porous media. Benefits and limits of the finite element poroelasticity formulation

The aim of this paper is to analyse the performance of a finite element formulation usable for predicting the mechanical consequences of frost effects on porous media. It considers the characteristics of porous media and how the frost action can be assessed. The problem is then separated into two parts: thermal and poromechanical calculations. The constitutive equations developed in the framework of poromechanics are presented and the implementation in a usual finite element poroelasticity formulation based on Zuber's method is adopted. An analysis of the time‐step influence on the convergence rate is given and leads us to propose a simple method in order to obtain objectivity of the finite element response and avoid over‐long calculations. Frost effect simulations are carried out on real porous media (two fired clays) as a case study. Although the experimental behaviour of the porous media subjected to frost action is in accordance with some observations, the calculated strains appear to be overestimated compared with measurements. The problem could be largely attributable to the difficulty of assessing permeability evolution during frost development. Copyright © 2011 John Wiley & Sons, Ltd.