Modeling of in-situ crystallization processes in the Permian mafic layered intrusion of Mont Collon (Dent Blanche nappe, western Alps)

The Mont Collon mafic complex is one of the best preserved examples of the Early Permian magmatism in the Central Alps, related to the intra-continental collapse of the Variscan belt. It mostly consists (> 95 vol.%) of ol+hy-nonnative plagioclase-wehrlites, olivine- and cpx-gabbros with cumulitic structures, crosscut by acid dikes. Pegmatitic gabbros, troctolites and anorthosites outcrop locally. A well-preserved cumulative, sequence is exposed in the Dents de Bertol area (center of intrusion). PT-calculations indicate that this layered magma chamber emplaced at mid-crustal levels at about 0.5 GPa and 1100 degrees C. The Mont Collon cumulitic rocks record little magmatic differentiation, as illustrated by the restricted range of clinopyroxene mg-number (Mg#(cpx)=83-89). Whole-rock incompatible trace-element contents (e.g. Nb, Zr, Ba) vary largely and without correlation with major-element composition. These features are characteristic of an in-situ crystallization process with variable amounts of interstitial liquid L trapped between the cumulus mineral phases. LA-ICPMS measurements show that trace-element distribution in the latter is homogeneous, pointing to subsolidus re-equilibration between crystals and interstitial melts. A quantitative modeling based on Langmuir's in-situ crystallization equation successfully duplicated the REE concentrations in cumulitic minerals of all rock facies of the intrusion. The calculated amounts of interstitial liquid L vary between 0 and 35% for degrees of differentiation F of 0 to 20%, relative to the least evolved facies of the intrusion. L values are well correlated with the modal proportions of interstitial amphibole and whole-rock incompatible trace-element concentrations (e.g. Zr, Nb) of the tested samples. However, the in-situ crystallization model reaches its limitations with rock containing high modal content of REE-bearing minerals (i.e. zircon), such as pegmatitic gabbros. Dikes of anorthositic composition, locally crosscutting the layered lithologies, evidence that the Mont Collon rocks evolved in open system with mixing of intercumulus liquids of different origins and possibly contrasting compositions. The proposed model is not able to resolve these complex open systems, but migrating liquids could be partly responsible for the observed dispersion of points in some correlation diagrams. Absence of significant differentiation with recurrent lithologies in the cumulitic pile of Dents de Bertol points to an efficiently convective magma chamber, with possible periodic replenishment, (c) 2005 Elsevier B.V. All rights reserved.

Permian-polysulphide-siderite-barite-haematite deposit Rude in Samoborska Gora Mts., Zagorje-Mid-Transdanubian zone of the Inner Dinarides

Samoborska Gora Mts. is situated within the westernmost part of the Zagorje-Mid-Transdanubian zone of the Internal Dinarides. The Samoborska Gora Mts. predominantly consists of Permian unmetamorphosed siliciclastic sediments and evaporites, overlain by Lower Triassic sediments. Rude mineralisation is hosted by Permian siliciclastic sediments, below gypsum and anhydrite strata. The central part of the deposit consists of a 1.5 km long stratabound mineralisation, grading laterally into ferruginous sandstone and protruding vertically into a gypsum-anhydrite layer. Siderite-polysulphide-barite-quartz veins are located below the stratabound mineralisation. The stratiform part of the deposit is situated above the stratabound and consists of haematite layer with barite concretions and veinlets. Late stage galena-barite veins overprint earlier types of mineralisation. The Rude ore deposit was generated by predominantly NaCl +/- CaCl(2)-H(2)O solutions. Detrital quartz from stratiform mineralisation contains fluid inclusions with salinities between 7 and 11 wt. % NaCl equ., homogenizing between 150 degrees C to 230 degrees C. Stratabound/siderite-polysulphide-barite-quartz vein type mineralisation was derived from solutions with salinities between 5 and 19 wt. % NaCl equ., homogenizing between 60 degrees C and 160 degrees C, while late stage galenabarite veins were precipitated from solutions with salinities between 11 and 16 wt. % NaCl equ., homogenizing between 100 degrees C to 140 degrees C. Fluid inclusion bulk leachate chemistry recorded Na(+)> Mg(2+)>K(+)>Ca(2+)>Li(+) and Cl-> SO(4)(2-) ions. Sulphur isotope composition of barites and overlying gypsum stems from Permian seawater sulphate, supported by increased Br(-) content, which follows successively the seawater evaporation line. The sulphur isotopic composition of sulphides varies between -0.2 and + 12.5 parts per thousand , as a result of thermal reduction of Permian marine sulphate. Ore-forming fluids were produced by hydrothermal convective cells (reflux brine model), and were derived primarily from Permian seawater, modified by evaporation and interaction with Permian sedimentary rocks. Rude deposits in Samoborska Gora Mts. may be declared as a prototype of the Permian siderite-polysulphide-barite deposits (products of rifting along the passive Gondwana margin), in the Inner Dinarides, and their equivalents extending northeastward into the Zagorje-Mid-Transdanubian Zone and the Gemerides, and southeastward to the Hellenide-Albanides.

Carbon dioxide in scree slope deposits: A pathway from atmosphere to pedogenic carbonate

A continuum of carbon, from atmospheric CO2 to secondary calcium carbonate, has been studied in a soil associ- ated with scree slope deposits in the Jura Mountains of Switzerland. This approach is based on former studies conducted in other environments. This C continuum includes atmospheric CO2, soil organic matter (SOM), soil CO2, dissolved inorganic carbon (DIC) in soil solutions, and secondary pedogenic carbonate. Soil parameters (pCO2, temperature, pH, Cmin and Corg contents), soil solution chemistry, and isotopic compositions of soil CO2, DIC, carbonate and soil organic matter (δ13CCO2, δ13CDIC, δ13Ccar and δ13CSOM values) have been monitored at different depths (from 20 to 140 cm) over one year. Results demonstrated that the carbon source in secondary carbonate (mainly needle fiber calcite) is related to the dissolved inorganic carbon, which is strongly dependent on soil respiration. The heterotrophic respiration, rather than the limestone parent material, seems to control the pedogenic carbon cycle. The correlation of δ13Corg values with Rock-Eval HI and OI indices demonstrates that, in a soil associated to scree slope deposits, the main process responsible for 13C-enrichment in SOM is related to bac- terial oxidative decarboxylation. Finally, precipitation of secondary calcium carbonate is enhanced by changes in soil pCO2 associated to the convective movement of air masses induced by temperature gradients (heat pump effect) in the highly porous scree slope deposits. The exportation of soil C-leachates from systems such as the one studied in this paper could partially explain the "gap in the European carbon budget" reported by recent studies.

Rigorous derivation of a nonlinear diffusion equation as fast-reaction limit of a continuous coagulation-fragmentation model with diffusion

Weak solutions of the spatially inhomogeneous (diffusive) Aizenmann-Bak model of coagulation-breakup within a bounded domain with homogeneous Neumann boundary conditions are shown to converge, in the fast reaction limit, towards local equilibria determined by their mass. Moreover, this mass is the solution of a nonlinear diffusion equation whose nonlinearity depends on the (size-dependent) diffusion coefficient. Initial data are assumed to have integrable zero order moment and square integrable first order moment in size, and finite entropy. In contrast to our previous result [CDF2], we are able to show the convergence without assuming uniform bounds from above and below on the number density of clusters.

A DGT Technique for Plutonium Bioavailability Measurements.

The toxicity of heavy metals in natural waters is strongly dependent on the local chemical environment. Assessing the bioavailability of radionuclides predicts the toxic effects to aquatic biota. The technique of diffusive gradients in thin films (DGT) is largely exploited for bioavailability measurements of trace metals in waters. However, it has not been applied for plutonium speciation measurements yet. This study investigates the use of DGT technique for plutonium bioavailability measurements in chemically different environments. We used a diffusion cell to determine the diffusion coefficients (D) of plutonium in polyacrylamide (PAM) gel and found D in the range of 2.06-2.29 × 10(-6) cm(2) s(-1). It ranged between 1.10 and 2.03 × 10(-6) cm(2) s(-1) in the presence of fulvic acid and in natural waters with low DOM. In the presence of 20 ppm of humic acid of an organic-rich soil, plutonium diffusion was hindered by a factor of 5, with a diffusion coefficient of 0.50 × 10(-6) cm(2) s(-1). We also tested commercially available DGT devices with Chelex resin for plutonium bioavailability measurements in laboratory conditions and the diffusion coefficients agreed with those from the diffusion cell experiments. These findings show that the DGT methodology can be used to investigate the bioaccumulation of the labile plutonium fraction in aquatic biota.

Stabilized Petrov-Galerkin methods for the convection-diffusion-reaction and the Helmholtz equations

We present two new stabilized high-resolution numerical methods for the convection–diffusion–reaction (CDR) and the Helmholtz equations respectively. The work embarks upon a priori analysis of some consistency recovery procedures for some stabilization methods belonging to the Petrov–Galerkin framework. It was found that the use of some standard practices (e.g. M-Matrices theory) for the design of essentially non-oscillatory numerical methods is not feasible when consistency recovery methods are employed. Hence, with respect to convective stabilization, such recovery methods are not preferred. Next, we present the design of a high-resolution Petrov–Galerkin (HRPG) method for the 1D CDR problem. The problem is studied from a fresh point of view, including practical implications on the formulation of the maximum principle, M-Matrices theory, monotonicity and total variation diminishing (TVD) finite volume schemes. The current method is next in line to earlier methods that may be viewed as an upwinding plus a discontinuity-capturing operator. Finally, some remarks are made on the extension of the HRPG method to multidimensions. Next, we present a new numerical scheme for the Helmholtz equation resulting in quasi-exact solutions. The focus is on the approximation of the solution to the Helmholtz equation in the interior of the domain using compact stencils. Piecewise linear/bilinear polynomial interpolation are considered on a structured mesh/grid. The only a priori requirement is to provide a mesh/grid resolution of at least eight elements per wavelength. No stabilization parameters are involved in the definition of the scheme. The scheme consists of taking the average of the equation stencils obtained by the standard Galerkin finite element method and the classical finite difference method. Dispersion analysis in 1D and 2D illustrate the quasi-exact properties of this scheme. Finally, some remarks are made on the extension of the scheme to unstructured meshes by designing a method within the Petrov–Galerkin framework.

Producing on-line ultrapure dialysis fluid.

The process of on-line generation of ultrapure dialysis fluid is a core prerequisite for the safe execution of modern renal replacement therapies such as on-line hemodiafiltration and high-flux hemodialysis. In these extracorporeal treatments with variable degrees of convection, significant volumes of plasma water are removed and replaced with dialysis fluid, which must occur without causing harm to the patient. Historically, on-line generation of sterile and pyrogen-free physiological substitution fluid by the process of membrane ultrafiltration of fresh dialysis fluid has its origin in hemofiltration, a purely convective therapy. Development of this and later therapies is described in the historical context of a successful effort over decades to overcome the above formidable challenge, which was provided jointly by pioneering clinical investigators and a resourceful dialysis industry.

Registre isotòpic d'alquenones i lípids terrestres en sediment marí: una reconstrucció paleoclimàtica a l'oceà Índic

Projecte de recerca elaborat a partir d’una estada a la Institute of mineralogy and geochemistry de la University of Lausanne, Suïssa, entre 2007 i 2009. Durant l’última dècada, la comunitat científica ha reconegut que les zones tropicals juguen un paper clau en els processos dinàmics que controlen el canvi climàtic global, probablement com a desencadenant dels canvis succeïts en altes latituds. A més a més, els sediments dels oceans tropicals, en trobar-se fora de l’impacte directe de les plaques de gel continentals creades durant les glaciacions, proporcionen un registre continu de les variacions climàtiques del planeta. Malgrat tot, encara hi ha moltes incògnites sobre el paper específic de les zones tropicals, especialment pel que fa a les variacions brusques suborbitals, degut als pocs registres d’alta resolució estudiats en aquestes àrees que abastin varis cicles glacial/interglacial. Per tal d’ajudar a clarificar el paper de les zones tropicals de l’hemisferi sud en el control del clima a escala mil•lenària s’ha estudiat la distribució i la composició isotòpica de biomarcadors moleculars marins i terrestres, a baixa resolució, en el testimoni MD98-2165 (9º39’S, 118º20’E, 2100 m de profunditat d’aigua, 42.3 m de llarg) està situat al sud-oest d’Indonèsia, on s’enregistren les temperatures superficials del mar més elevades del planeta i una elevada activitat convectiva, que té una influència en la distribució de la humitat atmosfèrica en una extensa superfície de la Terra. Les distribucions observades de biomarcadors terrígens (C23-C33 n-alcans i C20-C32 n-alcan-1-ols) són típiques del lipids de plantes superiors que arriben a l’oceà principalment per via eòlica. L’alcà de 31 àtoms de carboni i els alcohols de 28 o 32 àtoms de carboni són els homòlegs més abundants en ambdós testimonis. Cal destacar l’alcohol C32 com a homòleg principal durant les èpoques glacials, tot suggerint una expansió de les plantes tropicals C4 associada a unes condicions més àrides. La procedència d’aquests lipids queda corroborada mitjançant la seva composició isotòpica de carboni, que ens permet diferenciar la ruta fotosintètica emprada i per tant, entre el tipus de plantes.

Desenvolupament de codis numèrics per a la resolució de fluxos turbulents en ordinadors paral·lels emprant discretitzacions tipus symmetry-preserving

Projecte de recerca elaborat a partir d’una estada a la University of Groningen, Holanda, entre 2007 i 2009. La simulació directa de la turbulència (DNS) és una eina clau dins de la mecànica de fluids computacional. Per una banda permet conèixer millor la física de la turbulència i per l'altra els resultats obtinguts són claus per el desenvolupament dels models de turbulència. No obstant, el DNS no és una tècnica vàlida per a la gran majoria d'aplicacions industrials degut al elevats costos computacionals. Per tant, és necessari cert grau de modelització de la turbulència. En aquest context, s'han introduïts importants millores basades en la modelització del terme convectiu (no lineal) emprant symmetry-preserving regularizations. En tracta de modificar adequadament el terme convectiu a fi de reduir la producció d'escales més i més petites (vortex-stretching) tot mantenint tots els invariants de les equacions originals. Fins ara, aquest models s'han emprat amb èxit per nombres de Rayleigh (Ra) relativament elevats. En aquest punt, disposar de resultats DNS per a configuracions més complexes i nombres de Ra més elevats és clau. En aquest contexte, s'han dut a terme simulacions DNS en el supercomputador MareNostrum d'una Differentially Heated Cavity amb Ra=1e11 i Pr=0.71 durant el primer any dels dos que consta el projecte. A més a més, s'ha adaptat el codi a fi de poder simular el fluxe al voltant d'un cub sobre una pared amb Re=10000. Aquestes simulacions DNS són les més grans fetes fins ara per aquestes configuracions i la seva correcta modelització és un gran repte degut la complexitat dels fluxes. Aquestes noves simulacions DNS estan aportant nous coneixements a la física de la turbulència i aportant resultats indispensables per al progrés de les modelitzacións tipus symmetry-preserving regularization.

Spatio-temporal dynamics and transition from asymptotic equilibrium to bounded oscillations in Chrysomya albiceps (Diptera, Calliphoridae)

The sensitivity of parameters that govern the stability of population size in Chrysomya albiceps and describe its spatial dynamics was evaluated in this study. The dynamics was modeled using a density-dependent model of population growth. Our simulations show that variation in fecundity and mainly in survival has marked effect on the dynamics and indicates the possibility of transitions from one-point equilibrium to bounded oscillations. C. albiceps exhibits a two-point limit cycle, but the introduction of diffusive dispersal induces an evident qualitative shift from two-point limit cycle to a one fixed-point dynamics. Population dynamics of C. albiceps is here compared to dynamics of Cochliomyia macellaria, C. megacephala and C. putoria.

Arrested kinetic Li isotope fractionation at the margin of the llimaussaq complex, South Greenland: Evidence for open-system processes during final cooling of peralkaline igneous rocks

Li contents [Li] and isotopic composition (delta Li-7) of mafic minerals (mainly amphibole and clinopyroxene) from the alkaline to peralkaline Ilimaussaq plutonic complex, South Greenland, track the behavior of Li and its isotopes during magmatic differentiation and final cooling of an alkaline igneous system. [Li] in amphibole increase from < 10 ppm in Caamphiboles of the least differentiated unit to >3000 ppm in Na-amphiboles of the highly evolved units. In contrast, [Li] in clinopyroxene are comparatively low (<85 ppm) and do not vary systematically with differentiation. The distribution of Li between amphibole and pyroxene is controlled by the major element composition of the minerals (Ca-rich and Na-rich, respectively) and changes in oxygen fugacity (due to Li incorporation via coupled substitution with ferric iron) during magmatic differentiation. delta(7) Li values of all minerals span a wide range from + 17 to - 8 parts per thousand, with the different intrusive units of the complex having distinct Li isotopic systematics. Amphiboles, which dominate the Li budget of whole-rocks from the inner part of the complex, have constant delta Li-7 of + 1.8 +/- 2.2 parts per thousand (2 sigma, n = 15). This value reflects a homogeneous melt reservoir and is consistent with their mantle derivation, in agreement with published O and Nd isotopic data. Clinopyroxenes of these samples are consistently lighter, with Delta Li-7(amph-cpx). as large as 8 parts per thousand and are thus not in Li isotope equilibrium. These low values probably reflect late-stage diffusion of Li into clinopyroxene during final cooling of the rocks, thus enriching the clinopyroxene in 6 Li. At the margin of the complex delta(7) Li in the syenites increases systematically, from +2 to high values of + 14 parts per thousand. This, coupled with the observed Li isotope systematics of the granitic country rocks, reflects post-magmatic open-system processes occurring during final cooling of the intrusion. Although the shape and magnitude of the Li isotope and elemental profiles through syenite and country rock are suggestive of diffusion-driven isotope fractionation, they cannot be modeled by one-dimensional diffusive transport and point to circulation of a fluid having a high 67 Li value (possibly seawater) along the chilled contact. In all, this study demonstrates that Li isotopes can be used to identify complex fluid- and diffusion-governed processes taking place during the final cooling of such rocks. (c) 2007 Elsevier B.V All rights reserved.

Preservation of an extreme transient geotherm in the Raft River detachment shear zone

Extensional detachment systems separate hot footwalls from cool hanging walls, but the degree to which this thermal gradient is the product of ductile or brittle deformation or a preserved original transient geotherm is unclear. Oxygen isotope thermometry using recrystallized quartz-muscovite pairs indicates a smooth thermal gradient (140 degrees C/100 m) across the gently dipping, quartzite-dominated detachment zone that bounds the Raft River core complex in northwest Utah (United States). Hydrogen isotope values of muscovite (delta D-Ms similar to-100 parts per thousand) and fluid inclusions in quartz (delta D-Fluid similar to-85 parts per thousand) indicate the presence of meteoric fluids during detachment dynamics. Recrystallized grain-shape fabrics and quartz c-axis fabric patterns reveal a large component of coaxial strain (pure shear), consistent with thinning of the detachment section. Therefore, the high thermal gradient preserved in the Raft River detachment reflects the transient geotherm that developed owing to shearing, thinning, and the potentially prominent role of convective flow of surface fluids.

Greenhouse gases and climatic states of minimum entropy production

The hypothesis of minimum entropy production is applied to a simple one-dimensional energy balance model and is analysed for different values of the radiative forcing due to greenhouse gases. The extremum principle is used to determine the planetary “conductivity” and to avoid the “diffusive” approximation, which is commonly assumed in this type of model. For present conditions the result at minimum radiative entropy production is similar to that obtained by applying the classical model. Other climatic scenarios show visible differences, with better behaviour for the extremal case

Caution on the storage of waters and aqueous solutions in plastic containers for hydrogen and oxygen stable isotope analysis

RATIONALE The choice of containers for storage of aqueous samples between their collection, transport and water hydrogen (2H) and oxygen (18O) stable isotope analysis is a topic of concern for a wide range of fields in environmental, geological, biomedical, food, and forensic sciences. The transport and separation of water molecules during water vapor or liquid uptake by sorption or solution and the diffusive transport of water molecules through organic polymer material by permeation or pervaporation may entail an isotopic fractionation. An experiment was conducted to evaluate the extent of such fractionation. METHODS Sixteen bottle-like containers of eleven different organic polymers, including low and high density polyethylene (LDPE and HDPE), polypropylene (PP), polycarbonate (PC), polyethylene terephthalate (PET), and perfluoroalkoxy-Teflon (PFA), of different wall thickness and size were completely filled with the same mineral water and stored for 659?days under the same conditions of temperature and humidity. Particular care was exercised to keep the bottles tightly closed and prevent loss of water vapor through the seals. RESULTS Changes of up to +5 parts per thousand for d2H values and +2.0 parts per thousand for d18O values were measured for water after more than 1?year of storage within a plastic container, with the magnitude of change depending mainly on the type of organic polymer, wall thickness, and container size. The most important variations were measured for the PET and PC bottles. Waters stored in glass bottles with Polyseal (TM) cone-lined PP screw caps and thick-walled HDPE or PFA containers with linerless screw caps having an integrally molded inner sealing ring preserved their original d2H and d18O values. The carbon, hydrogen, and oxygen stable isotope compositions of the organic polymeric materials were also determined. CONCLUSIONS The results of this study clearly show that for precise and accurate measurements of the water stable isotope composition in aqueous solutions, rigorous sampling and storage procedures are needed both for laboratory standards and for unknown samples. Copyright (c) 2012 John Wiley & Sons, Ltd.

Effects of convection on the removal of the multiplicity of stable states in one-dimensional vertical models

Here I develop a model of a radiative-convective atmosphere with both radiative and convective schemes highly simplified. The atmospheric absorption of radiation at selective wavelengths makes use of constant mass absorption coefficients in finite width spectral bands. The convective regime is introduced by using a prescribed lapse rate in the troposphere. The main novelty of the radiative-convective model developed here is that it is solved without using any angular approximation for the radiation field. The solution obtained in the purely radiation mode (i. e. with convection ignored) leads to multiple equilibria of stable states, being very similar to some results recently found in simple models of planetary atmospheres. However, the introduction of convective processes removes the multiple equilibria of stable states. This shows the importance of taking convective processes into account even for qualitative analyses of planetary atmosphere