Characterization of magmatic sulfur in the Aegean island arc by means of the δ34S values of fumarolic H2S, elemental S, and hydrothermal gypsum from Nisyros and Milos islands

A delta(34)S value of +6.3 +/- 1.5% was estimated for the rhyodacitic degassing magma present underneath the hydrothermal system of Nisyros, based on the S isotope ratios of H2S in fumarolic vapors. This value was estimated by modeling the irreversible water-rock mass transfers occurring during the generation of the hydrothermal liquid which separates these fumarolic vapors. The S isotope ratio of the rhyodacitic degassing magma of Nisyros is consistent with fractional crystallization of a parent basaltic magma with an initial delta(34)S value of +4% (+/-at least 1.5%). This positive value could be explained by mantle contamination due to by either transference of fluids derived from subducted materials or involvement of altered oceanic crust, whereas contribution of biogenic sulfides from sediments seems to be negligible or nil. This conclusion agrees with the lack of N-2 and CO2 from thermal decomposition of organic matter contained in subducted sediments, which is a characteristic of the whole Aegean arc system. Since hydrothermal S at Milos and Santorini has isotope ratios similar to those determined at Nisyros, it seems likely that common controlling processes are active throughout the Aegean island arc. (C) 2002 Elsevier, Science B.V. All rights reserved.

Sex differences in urinary levels of several biological indicators of exposure: a simulation study using a compartmental based toxicokinetic model

Toxicokinetic modeling is a useful tool to describe or predict the behavior of a chemical agent in the human or animal organism. A general model based on four compartments was developed in a previous study in order to quantify the effect of human variability on a wide range of biological exposure indicators. The aim of this study was to adapt this existing general toxicokinetic model to three organic solvents, which were methyl ethyl ketone, 1-methoxy-2-propanol and 1,1,1,-trichloroethane, and to take into account sex differences. We assessed in a previous human volunteer study the impact of sex on different biomarkers of exposure corresponding to the three organic solvents mentioned above. Results from that study suggested that not only physiological differences between men and women but also differences due to sex hormones levels could influence the toxicokinetics of the solvents. In fact the use of hormonal contraceptive had an effect on the urinary levels of several biomarkers, suggesting that exogenous sex hormones could influence CYP2E1 enzyme activity. These experimental data were used to calibrate the toxicokinetic models developed in this study. Our results showed that it was possible to use an existing general toxicokinetic model for other compounds. In fact, most of the simulation results showed good agreement with the experimental data obtained for the studied solvents, with a percentage of model predictions that lies within the 95% confidence interval varying from 44.4 to 90%. Results pointed out that for same exposure conditions, men and women can show important differences in urinary levels of biological indicators of exposure. Moreover, when running the models by simulating industrial working conditions, these differences could even be more pronounced. In conclusion, a general and simple toxicokinetic model, adapted for three well known organic solvents, allowed us to show that metabolic parameters can have an important impact on the urinary levels of the corresponding biomarkers. These observations give evidence of an interindividual variablity, an aspect that should have its place in the approaches for setting limits of occupational exposure.

Prediction of Infant Drug Exposure Through Breastfeeding: Population PK Modeling and Simulation of Fluoxetine Exposure.

The likelihood of significant exposure to drugs in infants through breast milk is poorly defined, given the difficulties of conducting pharmacokinetics (PK) studies. Using fluoxetine (FX) as an example, we conducted a proof-of-principle study applying population PK (popPK) modeling and simulation to estimate drug exposure in infants through breast milk. We simulated data for 1,000 mother-infant pairs, assuming conservatively that the FX clearance in an infant is 20% of the allometrically adjusted value in adults. The model-generated estimate of the milk-to-plasma ratio for FX (mean: 0.59) was consistent with those reported in other studies. The median infant-to-mother ratio of FX steady-state plasma concentrations predicted by the simulation was 8.5%. Although the disposition of the active metabolite, norfluoxetine, could not be modeled, popPK-informed simulation may be valid for other drugs, particularly those without active metabolites, thereby providing a practical alternative to conventional PK studies for exposure risk assessment in this population.

Accurate and efficient simulation of multiphase flow in a heterogeneous reservoir by using error estimate and control in the multiscale finite-volume framework

The multiscale finite-volume (MSFV) method is designed to reduce the computational cost of elliptic and parabolic problems with highly heterogeneous anisotropic coefficients. The reduction is achieved by splitting the original global problem into a set of local problems (with approximate local boundary conditions) coupled by a coarse global problem. It has been shown recently that the numerical errors in MSFV results can be reduced systematically with an iterative procedure that provides a conservative velocity field after any iteration step. The iterative MSFV (i-MSFV) method can be obtained with an improved (smoothed) multiscale solution to enhance the localization conditions, with a Krylov subspace method [e.g., the generalized-minimal-residual (GMRES) algorithm] preconditioned by the MSFV system, or with a combination of both. In a multiphase-flow system, a balance between accuracy and computational efficiency should be achieved by finding a minimum number of i-MSFV iterations (on pressure), which is necessary to achieve the desired accuracy in the saturation solution. In this work, we extend the i-MSFV method to sequential implicit simulation of time-dependent problems. To control the error of the coupled saturation/pressure system, we analyze the transport error caused by an approximate velocity field. We then propose an error-control strategy on the basis of the residual of the pressure equation. At the beginning of simulation, the pressure solution is iterated until a specified accuracy is achieved. To minimize the number of iterations in a multiphase-flow problem, the solution at the previous timestep is used to improve the localization assumption at the current timestep. Additional iterations are used only when the residual becomes larger than a specified threshold value. Numerical results show that only a few iterations on average are necessary to improve the MSFV results significantly, even for very challenging problems. Therefore, the proposed adaptive strategy yields efficient and accurate simulation of multiphase flow in heterogeneous porous media.

The roles of flux- and decompression melting and their respective fractionation lines for continental crust formation: Evidence from the Kohistan arc

Delamination and foundering of the lower continental crust (LCC) into the mantle is part of the crust-forming mechanism. However, knowledge of the composition and mineralogy of the preserved or delaminated LCC over geological timescales remains scarce. We provide a synopsis of recent research within the Kohistan arc (Pakistan) and demonstrate that hydrous and less hydrous liquid lines of descent related to flux assisted and decompression mantle melting, respectively, produce compositionally different lower crustal rocks. The argument refers to two lower crustal sections exposed in Kohistan, the older Southern Plutonic Complex (SPC) and the younger Chilas Complex. The SPC typifies a hydrous, high-pressure fractionation sequence of olivine-pyroxenes-garnet-Fe/Ti-oxide-amphibole-plagioclase. The Chilas Complex illustrates a less hydrous fractionation sequence of olivine-clinopyroxene-orthopyroxene-plagioclase-amphibole. Despite the similarity of the Chilas Complex rocks to proposed lower crust compositions, the less hydrous fractionation results in unrealistically small volumes of silica-rich rocks, precluding the Chilas Complex gabbros to represent the magmatic complement to the upper crust. The composition of the SPC lower crust differs markedly from bulk lower crust estimates, but is complementary to silica-rich rocks exposed along this section and in the Kohistan batholith. These observations inspire a composite model for the formation of continental crust (CC) where the negatively buoyant delaminated and the buoyant preserved lower continental crusts (LCC) differ in genesis, mineralogy, and composition. We propose that the upper, non-sedimentary subsequent removal of the complementary, negatively buoyant garnet-pyroxene-amphibole-plagioclase-rich cumulates. In contrast, the LCC, which is buoyant and preserved over geological timescales, is formed by less hydrous parental mantle melts. We suggest that the bulk continental crust composition is related to mixing of these petrologically not directly related end members. Published by Elsevier B.V.

Monte Carlo simulation of a whole-body counter using IGOR phantoms.

Whole-body counting is a technique of choice for assessing the intake of gamma-emitting radionuclides. An appropriate calibration is necessary, which is done either by experimental measurement or by Monte Carlo (MC) calculation. The aim of this work was to validate a MC model for calibrating whole-body counters (WBCs) by comparing the results of computations with measurements performed on an anthropomorphic phantom and to investigate the effect of a change in phantom's position on the WBC counting sensitivity. GEANT MC code was used for the calculations, and an IGOR phantom loaded with several types of radionuclides was used for the experimental measurements. The results show a reasonable agreement between measurements and MC computation. A 1-cm error in phantom positioning changes the activity estimation by >2%. Considering that a 5-cm deviation of the positioning of the phantom may occur in a realistic counting scenario, this implies that the uncertainty of the activity measured by a WBC is ∼10-20%.

Mesozoic arc magmatism along the southern Peruvian margin during Gondwana breakup and dispersal

A high-resolution U-Pb zircon geochronological study of plutonic units along the south Peruvian margin between 17 degrees and 18 degrees S allows the integration of the geochemical, geodynamic and tectonic evolution of this part of the Andean margin. This study focuses on the composite Jurassic-early Cretaceous Ilo Batholith that was emplaced along the southern Peruvian coast during two episodes of intrusive magmatism; a first period between 173 and 152 Ma (with a peak in magmatic activity between roughly 168 and 162 Ma) and a second period between 110 and 106 Ma. Emplacement of the Jurassic part of the composite Ilo Batholith shortly post-dated the accumulation of the volcanosedimentary succession it intruded (Chocolate formation), which allows to estimate a subsidence rate for this unit of similar to 3.5 km/Ma. The emplacement of the main peak of Jurassic plutonism of the Ilo Batholith was also closely coeval with widespread and repeated slumping (during deposition of the Cachios Formation) in the back-arc region, suggesting a common causal link between these phenomena, which is discussed in the context of an observed 100 km trenchward arc migration at similar to 175 Ma, and the relation with extensional tectonics that prevailed along the Central Andean margin during Pangaea break-up. (C) 2012 Elsevier B.V. All rights reserved.

Mineral inventory of continuously erupting basaltic andesites at Arenal volcano, Costa Rica: implications for interpreting monotonous, crystal-rich, mafic arc stratigraphies

Except for the first 2 years since July 29, 1968, Arenal volcano has continuously erupted compositionally monotonous and phenocryst-rich (similar to35%) basaltic andesites composed of plagioclase (plag), orthopyroxene (opx), clinopyroxene (cpx), spinel olivine. Detailed textural and compositional analyses of phenocrysts, mineral inclusions, and microlites reveal comparable complexities in any given sample and identify mineral components that require a minimum of four crystallization environments. We suggest three distinct crystallization environments crystallized low Mg# (<78) silicate phases from andesitic magma but at different physical conditions, such as variable pressure of crystallization and water conditions. The dominant environment, i.e., the one which accounts for the majority of minerals and overprinted all other assemblages near rims of phenocrysts, cocrystallized clinopyroxene (Mg# similar to71-78), orthopyroxene (Mg# similar to71-78), titanomagnetite and plagioclase (An(60) to An(85)). The second environment cocrystallized clinopyroxene (Mg# 71-78), olivine (<Fo(78)), titanomagnetite, and very high An (similar to90) plagioclase, while the third cocrystallized clinopyroxene (Mg# 71-78) with high (>7) Al/Ti and high (>4 wt.%) Al2O3, titanomagnetite with considerable Al2O3 (10-18 wt.%) and possibly olivine but appears to lack plagioclase. A fourth crystallization environment is characterized by clinopyroxene (e.g., Mg#=similar to78-85; Cr2O3=0.15-0.7 wt.%), Al-, Cr-rich spinel olivine (similar toFo(80)), and in some circumstances high-An (>80) plagioclase. This assemblage seems to record mafic inputs into the Arenal system and crystallization at high to low pressures. Single crystals cannot be completely classified as xenocrysts, antecrysts (cognate crystals), or phenocrysts, because they often contain different parts each representing a different crystallization environment and thus belong to different categories. Bulk compositions are mostly too mafic to have crystallized the bulk of ferromagnesian minerals and thus likely do not represent liquid compositions. On the other hand, they are the cumulative products of multiple mixing events assembling melts and minerals from a variety of sources. The driving force for this multistage mixing evolution to generate erupting basaltic andesites is thought to be the ascent of mafic magma from lower crustal levels to subvolcanic depths which at the same time may also go through compositional modification by fractionation and assimilation of country rocks. Thus, mafic magmas become basaltic andesite through mixing, fractionation and assimilation by the time they arrive at subvolcanic depths. We infer new increments of basaltic andesite are supplied nearly continuously to the subvolcanic reservoir concurrently to the current eruption and that these new increments are blended into the residing, subvolcanic magma. Thus, the compositional monotony is mostly the product of repetitious production of very similar basaltic andesite. Furthermore, we propose that this quasi-constant supply of small increments of magma is the fundamental cause for small-scale, decade-long continuous volcanic activity; that is, the current eruption of Arenal is flux-controlled by inputs of mantle magmas. (C) 2004 Elsevier B.V. All rights reserved.

Simulation of interface dynamics across pore discontinuities

We study the dynamics of a water-oil meniscus moving from a smaller to a larger pore. The process is characterised by an abrupt change in the configuration, yielding a sudden energy release. A theoretic study for static conditions provides analytical solutions of the surface energy content of the system. Although the configuration after the sudden energy release is energetically more convenient, an energy barrier must be overcome before the process can happen spontaneously. The energy barrier depends on the system geometry and on the flow parameters. The analytical results are compared to numerical simulations that solve the full Navier-Stokes equation in the pore space and employ the Volume Of Fluid (VOF) method to track the evolution of the interface. First, the numerical simulations of a quasi-static process are validated by comparison with the analytical solutions for a static meniscus, then numerical simulations with varying injection velocity are used to investigate dynamic effects on the configuration change. During the sudden energy jump the system exhibits an oscillatory behaviour. Extension to more complex geometries might elucidate the mechanisms leading to a dynamic capillary pressure and to bifurcations in final distributions of fluid phases in porous