Mixing-cell boundary conditions and apparent mass balance errors for advective-dispersive solute transport

In many field or laboratory situations, well-mixed reservoirs like, for instance, injection or detection wells and gas distribution or sampling chambers define boundaries of transport domains. Exchange of solutes or gases across such boundaries can occur through advective or diffusive processes. First we analyzed situations, where the inlet region consists of a well-mixed reservoir, in a systematic way by interpreting them in terms of injection type. Second, we discussed the mass balance errors that seem to appear in case of resident injections. Mixing cells (MC) can be coupled mathematically in different ways to a domain where advective-dispersive transport occurs: by assuming a continuous solute flux at the interface (flux injection, MC-FI), or by assuming a continuous resident concentration (resident injection). In the latter case, the flux leaving the mixing cell can be defined in two ways: either as the value when the interface is approached from the mixing-cell side (MC-RT -), or as the value when it is approached from the column side (MC-RT +). Solutions of these injection types with constant or-in one case-distance-dependent transport parameters were compared to each other as well as to a solution of a two-layer system, where the first layer was characterized by a large dispersion coefficient. These solutions differ mainly at small Peclet numbers. For most real situations, the model for resident injection MC-RI + is considered to be relevant. This type of injection was modeled with a constant or with an exponentially varying dispersion coefficient within the porous medium. A constant dispersion coefficient will be appropriate for gases because of the Eulerian nature of the usually dominating gaseous diffusion coefficient, whereas the asymptotically growing dispersion coefficient will be more appropriate for solutes due to the Lagrangian nature of mechanical dispersion, which evolves only with the fluid flow. Assuming a continuous resident concentration at the interface between a mixing cell and a column, as in case of the MC-RI + model, entails a flux discontinuity. This flux discontinuity arises inherently from the definition of a mixing cell: the mixing process is included in the balance equation, but does not appear in the description of the flux through the mixing cell. There, only convection appears because of the homogeneous concentration within the mixing cell. Thus, the solute flux through a mixing cell in close contact with a transport domain is generally underestimated. This leads to (apparent) mass balance errors, which are often reported for similar situations and erroneously used to judge the validity of such models. Finally, the mixing cell model MC-RI + defines a universal basis regarding the type of solute injection at a boundary. Depending on the mixing cell parameters, it represents, in its limits, flux as well as resident injections. (C) 1998 Elsevier Science B.V. All rights reserved.

Glacier mass balance reconstruction by sublimation induced enrichment of chemical species on Cerro Tapado (Chilean Andes)

A 36 m long ice core down to bedrock from the Cerro Tapado glacier (5536 m a.s.l, 30°08' S, 69°55' W) was analyzed to reconstruct past climatic conditions for Northern Chile. Because of the marked seasonality in the precipitation (short wet winter and extended dry summer periods) in this region, major snow ablation and related post-depositional processes occur on the glacier surface during summer periods. They include predominantly sublimation and dry deposition. Assuming that, like measured during the field campaign, the enrichment of chloride was always related to sublimation, the chemical record along the ice core may be applied to reconstruct the history of such secondary processes linked to the past climatic conditions over northern Chile. For the time period 1962–1999, a mean annual net accumulation of 316 mm water equivalent (weq) and 327 mm weq loss by sublimation was deduced by this method. This corresponds to an initial total annual accumulation of 539 mm weq. The annual variability of the accumulation and sublimation is related with the Southern Oscillation Index (SOI): higher net-accumulation during El-Niño years and more sublimation during La Niña years. The deepest part of the ice record shows a time discontinuity; with an ice body deposited under different climatic conditions: 290 mm higher precipitation but with reduced seasonal distribution (+470 mm in winter and –180 mm in summer) and –3°C lower mean annual temperature. Unfortunately, its age is unknown. The comparison with regional proxy data however let us conclude that the glacier buildup did most likely occur after the dry mid-Holocene.

Protein and lipid binding parameters in rainbow trout (Oncorhynchus mykiss) blood and liver fractions to extrapolate from an in vitro metabolic degradation assay to in vivo bioaccumulation potential of hydrophobic organic chemicals

Binding of hydrophobic chemicals to colloids such as proteins or lipids is difficult to measure using classical microdialysis methods due to low aqueous concentrations, adsorption to dialysis membranes and test vessels, and slow kinetics of equilibration. Here, we employed a three-phase partitioning system where silicone (polydimethylsiloxane, PDMS) serves as a third phase to determine partitioning between water and colloids and acts at the same time as a dosing device for hydrophobic chemicals. The applicability of this method was demonstrated with bovine serum albumin (BSA). Measured binding constants (K(BSAw)) for chlorpyrifos, methoxychlor, nonylphenol, and pyrene were in good agreement with an established quantitative structure-activity relationship (QSAR). A fifth compound, fluoxypyr-methyl-heptyl ester, was excluded from the analysis because of apparent abiotic degradation. The PDMS depletion method was then used to determine partition coefficients for test chemicals in rainbow trout (Oncorhynchus mykiss) liver S9 fractions (K(S9w)) and blood plasma (K(bloodw)). Measured K(S9w) and K(bloodw) values were consistent with predictions obtained using a mass-balance model that employs the octanol-water partition coefficient (K(ow)) as a surrogate for lipid partitioning and K(BSAw) to represent protein binding. For each compound, K(bloodw) was substantially greater than K(S9w), primarily because blood contains more lipid than liver S9 fractions (1.84% of wet weight vs 0.051%). Measured liver S9 and blood plasma binding parameters were subsequently implemented in an in vitro to in vivo extrapolation model to link the in vitro liver S9 metabolic degradation assay to in vivo metabolism in fish. Apparent volumes of distribution (V(d)) calculated from the experimental data were similar to literature estimates. However, the calculated binding ratios (f(u)) used to relate in vitro metabolic clearance to clearance by the intact liver were 10 to 100 times lower than values used in previous modeling efforts. Bioconcentration factors (BCF) predicted using the experimental binding data were substantially higher than the predicted values obtained in earlier studies and correlated poorly with measured BCF values in fish. One possible explanation for this finding is that chemicals bound to proteins can desorb rapidly and thus contribute to metabolic turnover of the chemicals. This hypothesis remains to be investigated in future studies, ideally with chemicals of higher hydrophobicity.

Anatomically-Driven Soft-Tissue Simulation Strategy for Cranio-Maxillofacial Surgery Using Facial Muscle Template Model

We propose a computationally efficient and biomechanically relevant soft-tissue simulation method for cranio-maxillofacial (CMF) surgery. A template-based facial muscle reconstruction was introduced to minimize the efforts on preparing a patient-specific model. A transversely isotropic mass-tensor model (MTM) was adopted to realize the effect of directional property of facial muscles in reasonable computation time. Additionally, sliding contact around teeth and mucosa was considered for more realistic simulation. Retrospective validation study with postoperative scan of a real patient showed that there were considerable improvements in simulation accuracy by incorporating template-based facial muscle anatomy and sliding contact.

Plant diversity effects on the water balance of an experimental grassland

In the literature, contrasting effects of plant species richness on the soil water balance are reported. Our objective was to assess the effects of plant species and functional richness and functional identity on soil water contents and water fluxes in the experimental grassland of the Jena Experiment. The Jena Experiment comprises 86 plots on which plant species richness (0, 1, 2, 4, 8, 16, and 60) and functional group composition (zero to four functional groups: legumes, grasses, tall herbs, and small herbs) were manipulated in a factorial design. We recorded meteorological data and soil water contents of the 0·0–0·3 and 0·3–0·7 m soil layers and calculated actual evapotranspiration (ETa), downward flux (DF), and capillary rise with a soil water balance model for the period 2003–2007. Missing water contents were estimated with a Bayesian hierarchical model. Species richness decreased water contents in subsoil during wet soil conditions. Presence of tall herbs increased soil water contents in topsoil during dry conditions and decreased soil water contents in subsoil during wet conditions. Presence of grasses generally decreased water contents in topsoil, particularly during dry phases; increased ETa and decreased DF from topsoil; and decreased ETa from subsoil. Presence of legumes, in contrast, decreased ETa and increased DF from topsoil and increased ETa from subsoil. Species richness probably resulted in complementary water use. Specific functional groups likely affected the water balance via specific root traits (e.g. shallow dense roots of grasses and deep taproots of tall herbs) or specific shading intensity caused by functional group effects on vegetation cover. Copyright © 2013 John Wiley & Sons, Ltd.

Chromium uptake and adsorption in marine phytoplankton - Implications for the marine chromium cycle

Using the radioisotope 51Cr, we investigated the controls of cellular Cr accumulation in an array of marine phytoplankton grown in environmentally relevant Cr concentrations (1–10 nM). Given the affinity of Cr(III) for amorphous Fe-hydroxide mineral surfaces, and the formation of these mineral phases on the outside of phytoplankton cells, extracellular Cr was monitored in a model diatom species (Thalassiosira weissflogii) as extracellular Fe concentrations varied. Extracellular Cr in T. weissflogii increased with increasing extracellular Fe, demonstrating that Cr may be removed from seawater via extracellular adsorption to phytoplankton. Short-term Cr(VI) and Cr(III) uptake experiments performed with T. weissflogii demonstrated that Cr(III) was the primary oxidation state adsorbing to cells and being internalized by them. Cellular Cr:C ratios (<0.5 μmol Cr mol C−1) of the eight phytoplankton species surveyed were significantly lower than previously reported Cr:C ratios in marine particles with a high biogenic component (10–300 μmol Cr mol C−1). This indicates that Cr(III) likely accumulates in marine particles due to uptake and/or adsorption. Mass balance calculations demonstrate that surface water Cr deficits can be explained via loss of Cr(III) to exported particles, thereby providing a mechanism to account for the nutrient depth profile for Cr in modern seawater. Given the large fractionation of stable Cr isotopes during Cr(VI) reduction, Cr(III) associated with exported organic carbon is likely enriched in lighter isotopes. Most sedimentary Cr isotope studies have thus far neglected internal fractionating processes in the marine Cr cycle, but our data indicate that loss of Cr to exported particles may be traced in the sedimentary d53Cr record.

Is Vesta an intact and pristine protoplanet?

It is difficult to find a Vesta model of iron core, pyroxene and olivine-rich mantle, and HED crust that can match the joint constraints of (a) Vesta's density and core size as reported by the Dawn spacecraft team; (b) the chemical trends of the HED meteorites, including the depletion of sodium, the FeO abundance, and the trace element enrichments; and (c) the absence of exposed mantle material on Vesta's surface, among Vestoid asteroids, or in our collection of basaltic meteorites. These conclusions are based entirely on mass-balance and density arguments, independent of any particular formation scenario for the HED meteorites themselves. We suggest that Vesta either formed from source material with non-chondritic composition or underwent after its formation a radical physical alteration, possibly caused by collisional processes, that affected its global composition and interior structure. (C) 2015 Elsevier Inc. All rights reserved.

Ventilation-perfusion ratio in perflubron during partial liquid ventilation

BACKGROUND: Functional magnetic resonance imaging (fMRI) of fluorine-19 allows for the mapping of oxygen partial pressure within perfluorocarbons in the alveolar space (Pao(2)). Theoretically, fMRI-detected Pao(2) can be combined with the Fick principle approach, i.e., a mass balance of oxygen uptake by ventilation and delivery by perfusion, to quantify the ventilation-perfusion ratio (Va/Q) of a lung region: The mixed venous blood and the inspiratory oxygen fraction, which are equal for all lung regions, are measured. In addition, the local expiratory oxygen fraction and the end capillary oxygen content, both of which may differ between the lung regions, are calculated using the fMRI-detected Pao(2). We investigated this approach by numerical simulations and applied it to quantify local Va/Q in the perfluorocarbons during partial liquid ventilation. METHODS: Numerical simulations were performed to analyze the sensitivity of the Va/Q calculation and to compare this approach with another one proposed by Rizi et al. in 2004 (Magn Reson Med 2004;52:65-72). Experimentally, the method was used during partial liquid ventilation in 7 anesthetized pigs. The Pao(2) distribution in intraalveolar perflubron was measured by fluorine-19 MRI. Respiratory gas fractions together with arterial and mixed venous blood samples were taken to quantify oxygen partial pressure and content. Using the Fick principle, the local Va/Q was estimated. The impact of gravity (nondependent versus dependent) of perflubron dose (10 vs 20 mL/kg body weight) and of inspired oxygen fraction (Fio(2)) (0.4-1.0) on Va/Q was examined. RESULTS: In numerical simulations, the Fick principle proved to be appropriate over the Va/Q range from 0.02 to 2.5. Va/Q values were in acceptable agreement with the method published by Rizi et al. In the experimental setting, low mean Va/Q values were found in perflubron (confidence interval [CI] 0.08-0.29 with 20 mL/kg perflubron). At this dose, Va/Q in the nondependent lung was higher (CI 0.18-0.39) than in the dependent lung regions (CI 0.06-0.16; P = 0.006; Student t test). Differences depending on Fio(2) or perflubron dose were, however, small. CONCLUSION: The results show that derivation of Va/Q from local Po(2) measurements using fMRI in perflubron is feasible. The low detected Va/Q suggests that oxygen transport into the perflubron-filled alveolar space is significantly restrained.

Effects of drilling and stress release on hydraulic properties and porewater chemistry of crystalline rocks

The experimental verification of matrix diffusion in crystalline rocks largely relies on indirect methods performed in the laboratory. Such methods are prone to perturbations of the rock samples by collection and preparation and therefore the laboratory-derived transport properties and fluid composition might not represent in situ conditions. We investigated the effects induced by the drilling process and natural rock stress release by mass balance considerations and sensitivity analysis of analytical out-diffusion data obtained from originally saturated, large-sized drillcore material from two locations drilled using traced drilling fluid. For in situ stress-released drillcores of quartz-monzodiorite composition from the Aspo HRL, Sweden, tracer mass balance considerations and 1D and 2D diffusion modelling consistently indicated a contamination of <1% of the original pore water. This chemically disturbed zone extends to a maximum of 0.1 mm into the drillcore (61.8 mm x 180.1 mm) corresponding to about 0.66% of the total pore volume (0.77 vol.%). In contrast, the combined effects of stress release and the drilling process, which have influenced granodioritic drillcore material from 560 m below surface at Forsmark. Sweden, resulted in a maximum contamination of the derived porewater Cl(-) concentration of about 8%. The mechanically disturbed zone with modified diffusion properties covers the outermost similar to 6 mm of the drillcore (50 mm x 189 mm), whereas the chemically disturbed zone extends to a maximum of 0.3 mm based on mass balance considerations, and to 0.15 mm to 0.2 mm into the drillcore based on fitting the observed tracer data. This corresponds to a maximum of 2.4% of the total pore volume (0.62 vol.%) being affected by the drilling-fluid contamination. The proportion of rock volume affected initially by drilling fluid or subsequently with experiment water during the laboratory diffusion and re-saturation experiments depends on the size of the drillcore material and will become larger the smaller the sample used for the experiment. The results are further in support of matrix diffusion taking place in the undisturbed matrix of crystalline rocks at least in the cm range.

Impacts of environmental change on water resources in the Mount Kenya Region

Water resources are becoming increasingly scarce in the Mt. Kenya region. Land use and climate change may pose additional challenges to water management in the future. In order to assess the impacts of environmental change, the NRM3 Streamflow Model, a simple, semi-distributed, grid-based water balance model, is evaluated as a tool for discharge prediction in six meso-scale catchments on the western slopes of Mt. Kenya, and used to analyse the impact of land use and climate change scenarios on water resources.

Studying the fate of non-volatile organic compounds in a commercial plasma air purifier

Degradation of non-volatile organic compounds-environmental toxins (methyltriclosane and phenanthrene), bovine serum albumin, as well as bioparticles (Legionella pneumophila, Bacillus subtilis, and Bacillus anthracis)-in a commercially available plasma air purifier based on a cold plasma was studied in detail, focusing on its efficiency and on the resulting degradation products. This system is capable of handling air flow velocities of up to 3.0m s(-1) (3200Lmin(-1)), much higher than other plasma-based reactors described in the literature, which generally are limited to air flow rates below 10Lmin(-1). Mass balance studies consistently indicated a reduction in concentration of the compounds/particles after passage through the plasma air purifier, 31% for phenanthrene, 17% for methyltriclosane, and 80% for bovine serum albumin. L. pneumophila did not survive passage through the plasma air purifier, and cell counts of aerosolized spores of B. subtilis and B. anthracis were reduced by 26- and 15-fold, depending on whether it was run at 10Hz or 50Hz, respectively. However rather than chemical degradation, deposition on the inner surfaces of the plasma air purifier occured. Our interpretation is that putative "degradation" efficiencies were largely due to electrostatic precipitation rather than to decomposition into smaller molecules.

Net accumulation rates derived from ice core stable isotope records of Pío XI glacier, Southern Patagonia Icefield

Pío XI, the largest glacier of the Southern Patagonia Icefield, reached its neoglacial maximum extent in 1994 and is one of the few glaciers in that area which is not retreating. In view of the recent warming it is important to understand glacier responses to climate changes. Due to its remoteness and the harsh conditions in Patagonia, no systematic mass balance studies have been performed. In this study we derived net accumulation rates for the period 2000–2006 from a 50 m (33.2 4 m weq) ice core collected in the accumulation area of Pío XI (2600 m a.s.l., 49°16'40"S, 73°21'14"W). Borehole temperatures indicate near temperate ice, but the average melt percent is only 16 ± 14%. Records of stable isotopes are well preserved and were used for identification of annual layers. Net accumulation rates range from 3.4–7.1 water equivalent (m weq) with an average of 5.8 m weq, comparable to precipitation amounts at the Chilean coast, but not as high as expected for the Icefield. Ice core stable isotope data correlate well with upper air temperatures and may be used as temperature proxy.

Modelling Lake Kivu water level variations over the last seven decades

This study aimed at analysing the hydrological changes in the Lake Kivu Basin over the last seven decades with focus on the response of the lake water level to meteorological factors and hydropower dam construction. Historical precipitation and lake water levels were acquired from literature, local agencies and from global databases in order to compile a coherent dataset. The net lake inflow was modelled using a soil water balance model and the water levels were reconstructed using a parsimonious lake water balance model. The soil water balance shows that 370 mm yr−1 (25%) of the precipitation in the catchment contributes to the runoff and baseflow whereas 1100 mm yr−1 (75%) contributes to the evapotranspiration. A review of the lake water balance resulted in the following estimates of hydrological contributions: 55%, 25%, and 20% of the overall inputs from precipitation, surface inflows, and subaquatic groundwater discharge, respectively. The overall losses were 58% and 42% for lake surface evaporation and outflow discharge, respectively. The hydrological model used indicated a remarkable sensitivity of the lake water levels to hydrometeorological variability up to 1977, when the outflow bed was artificially widened.

Alpine lacustrine varved record reveals summer temperature as main control of glacier fluctuations over the past 2250 years

Glacier fluctuations are a key indicator of changing climate. Their reconstruction beyond historical times unravels glacier variability and its forcing factors on long time scales, which can considerably improve our understanding of the climate–glacier relationship. Here, we present a 2250-year-long reconstruction of particle-mass accumulation rates recorded in the lacustrine sediments of Lake Trüebsee (Central Swiss Alps) that are directly related to glacier extent, thus reflecting a continuous record of fluctuations of the upstream-located Titlis Glacier. Mass accumulation rate values show strong centennial to multi-centennial fluctuations and reveal 12 well-pronounced periods of enhanced values corresponding to times of maximum extent of the neighboring Lower Grindelwald Glacier. This result supports previous studies of proglacial lake sediments that documented high mass accumulation rate values during glacier advances. The strong variability in the Lake Trüebsee mass accumulation rate record thus represents a highly sensitive paleoclimatic archive, which mirrors rapid and pronounced feedbacks of Titlis Glacier to climatic changes over the past 2250years. The comparison of our data with independent paleo-temperature reconstructions from tree rings suggests that variations in mean summer temperature were the primary driving factor of fluctuations of Titlis Glacier. Also, advances of Titlis Glacier occurred during the grand solar minima (Dalton, Maunder, Spörer, Wolf) of the last millennium. This relation of glacier extent with summer temperature reveals strong evidence that the mass balance of this Alpine glacier is primarily controlled by the intensity of glacier melting during summer.