Solute transport in crystalline rocks at Äspö — I: Geological basis and model calibration

Water-conducting faults and fractures were studied in the granite-hosted A¨ spo¨ Hard Rock Laboratory (SE Sweden). On a scale of decametres and larger, steeply dipping faults dominate and contain a variety of different fault rocks (mylonites, cataclasites, fault gouges). On a smaller scale, somewhat less regular fracture patterns were found. Conceptual models of the fault and fracture geometries and of the properties of rock types adjacent to fractures were derived and used as input for the modelling of in situ dipole tracer tests that were conducted in the framework of the Tracer Retention Understanding Experiment (TRUE-1) on a scale of metres. After the identification of all relevant transport and retardation processes, blind predictions of the breakthroughs of conservative to moderately sorbing tracers were calculated and then compared with the experimental data. This paper provides the geological basis and model calibration, while the predictive and inverse modelling work is the topic of the companion paper [J. Contam. Hydrol. 61 (2003) 175]. The TRUE-1 experimental volume is highly fractured and contains the same types of fault rocks and alterations as on the decametric scale. The experimental flow field was modelled on the basis of a 2D-streamtube formalism with an underlying homogeneous and isotropic transmissivity field. Tracer transport was modelled using the dual porosity medium approach, which is linked to the flow model by the flow porosity. Given the substantial pumping rates in the extraction borehole, the transport domain has a maximum width of a few centimetres only. It is concluded that both the uncertainty with regard to the length of individual fractures and the detailed geometry of the network along the flowpath between injection and extraction boreholes are not critical because flow is largely one-dimensional, whether through a single fracture or a network. Process identification and model calibration were based on a single uranine breakthrough (test PDT3), which clearly showed that matrix diffusion had to be included in the model even over the short experimental time scales, evidenced by a characteristic shape of the trailing edge of the breakthrough curve. Using the geological information and therefore considering limited matrix diffusion into a thin fault gouge horizon resulted in a good fit to the experiment. On the other hand, fresh granite was found not to interact noticeably with the tracers over the time scales of the experiments. While fracture-filling gouge materials are very efficient in retarding tracers over short periods of time (hours–days), their volume is very small and, with time progressing, retardation will be dominated by altered wall rock and, finally, by fresh granite. In such rocks, both porosity (and therefore the effective diffusion coefficient) and sorption Kds are more than one order of magnitude smaller compared to fault gouge, thus indicating that long-term retardation is expected to occur but to be less pronounced.

Transport of iodide in structured clay–loam soil under maize during irrigation experiments analyzed using HYDRUS model

Transport of radioactive iodide 131I− in a structured clay loam soil under maize in a final growing phase was monitored during five consecutive irrigation experiments under ponding. Each time, 27 mm of water were applied. The water of the second experiment was spiked with 200 MBq of 131I− tracer. Its activity was monitored as functions of depth and time with Geiger-Müller (G-M) detectors in 11 vertically installed access tubes. The aim of the study was to widen our current knowledge of water and solute transport in unsaturated soil under different agriculturally cultivated settings. It was supposed that the change in 131I− activity (or counting rate) is proportional to the change in soil water content. Rapid increase followed by a gradual decrease in 131I− activity occurred at all depths and was attributed to preferential flow. The iodide transport through structured soil profile was simulated by the HYDRUS 1D model. The model predicted relatively deep percolation of iodide within a short time, in a good agreement with the observed vertical iodide distribution in soil. We found that the top 30 cm of the soil profile is the most vulnerable layer in terms of water and solute movement, which is the same depth where the root structure of maize can extend.

A numerical model for inert gas transport in the lung based on fractal airway morphology

Patients suffering from cystic fibrosis (CF) show thick secretions, mucus plugging and bronchiectasis in bronchial and alveolar ducts. This results in substantial structural changes of the airway morphology and heterogeneous ventilation. Disease progression and treatment effects are monitored by so-called gas washout tests, where the change in concentration of an inert gas is measured over a single or multiple breaths. The result of the tests based on the profile of the measured concentration is a marker for the severity of the ventilation inhomogeneity strongly affected by the airway morphology. However, it is hard to localize underlying obstructions to specific parts of the airways, especially if occurring in the lung periphery. In order to support the analysis of lung function tests (e.g. multi-breath washout), we developed a numerical model of the entire airway tree, coupling a lumped parameter model for the lung ventilation with a 4th-order accurate finite difference model of a 1D advection-diffusion equation for the transport of an inert gas. The boundary conditions for the flow problem comprise the pressure and flow profile at the mouth, which is typically known from clinical washout tests. The natural asymmetry of the lung morphology is approximated by a generic, fractal, asymmetric branching scheme which we applied for the conducting airways. A conducting airway ends when its dimension falls below a predefined limit. A model acinus is then connected to each terminal airway. The morphology of an acinus unit comprises a network of expandable cells. A regional, linear constitutive law describes the pressure-volume relation between the pleural gap and the acinus. The cyclic expansion (breathing) of each acinus unit depends on the resistance of the feeding airway and on the flow resistance and stiffness of the cells themselves. Special care was taken in the development of a conservative numerical scheme for the gas transport across bifurcations, handling spatially and temporally varying advective and diffusive fluxes over a wide range of scales. Implicit time integration was applied to account for the numerical stiffness resulting from the discretized transport equation. Local or regional modification of the airway dimension, resistance or tissue stiffness are introduced to mimic pathological airway restrictions typical for CF. This leads to a more heterogeneous ventilation of the model lung. As a result the concentration in some distal parts of the lung model remains increased for a longer duration. The inert gas concentration at the mouth towards the end of the expirations is composed of gas from regions with very different washout efficiency. This results in a steeper slope of the corresponding part of the washout profile.

Gas transport in firn: multiple-tracer characterisation and model intercomparison for NEEM, Northern Greenland

Air was sampled from the porous firn layer at the NEEM site in Northern Greenland. We use an ensemble of ten reference tracers of known atmospheric history to characterise the transport properties of the site. By analysing uncertainties in both data and the reference gas atmospheric histories, we can objectively assign weights to each of the gases used for the depth-diffusivity reconstruction. We define an objective root mean square criterion that is minimised in the model tuning procedure. Each tracer constrains the firn profile differently through its unique atmospheric history and free air diffusivity, making our multiple-tracer characterisation method a clear improvement over the commonly used single-tracer tuning. Six firn air transport models are tuned to the NEEM site; all models successfully reproduce the data within a 1σ Gaussian distribution. A comparison between two replicate boreholes drilled 64 m apart shows differences in measured mixing ratio profiles that exceed the experimental error. We find evidence that diffusivity does not vanish completely in the lock-in zone, as is commonly assumed. The ice age- gas age difference (1 age) at the firn-ice transition is calculated to be 182+3−9 yr. We further present the first intercomparison study of firn air models, where we introduce diagnostic scenarios designed to probe specific aspects of the model physics. Our results show that there are major differences in the way the models handle advective transport. Furthermore, diffusive fractionation of isotopes in the firn is poorly constrained by the models, which has consequences for attempts to reconstruct the isotopic composition of trace gases back in time using firn air and ice core records.

Establishment of a confluent monolayer model with human primary trophoblast cells: novel insights into placental glucose transport

STUDY HYPOTHESIS Using optimized conditions, primary trophoblast cells isolated from human term placenta can develop a confluent monolayer in vitro, which morphologically and functionally resembles the microvilli structure found in vivo. STUDY FINDING We report the successful establishment of a confluent human primary trophoblast monolayer using pre-coated polycarbonate inserts, where the integrity and functionality was validated by cell morphology, biophysical features, cellular marker expression and secretion, and asymmetric glucose transport. WHAT IS KNOWN ALREADY Human trophoblast cells form the initial barrier between maternal and fetal blood to regulate materno-fetal exchange processes. Although the method for isolating pure human cytotrophoblast cells was developed almost 30 years ago, a functional in vitro model with primary trophoblasts forming a confluent monolayer is still lacking. STUDY DESIGN, SAMPLES/MATERIALS, METHODS Human term cytotrophoblasts were isolated by enzymatic digestion and density gradient separation. The purity of the primary cells was evaluated by flow cytometry using the trophoblast-specific marker cytokeratin 7, and vimentin as an indicator for potentially contaminating cells. We screened different coating matrices for high cell viability to optimize the growth conditions for primary trophoblasts on polycarbonate inserts. During culture, cell confluency and polarity were monitored daily by determining transepithelial electrical resistance (TEER) and permeability properties of florescent dyes. The time course of syncytia-related gene expression and hCG secretion during syncytialization were assessed by quantitative RT-PCR and enzyme-linked immunosorbent assay, respectively. The morphology of cultured trophoblasts after 5 days was determined by light microscopy, scanning electron microscopy (SEM) and transmission electron microscopy (TEM). Membrane makers were visualized using confocal microscopy. Additionally, glucose transport studies were performed on the polarized trophoblasts in the same system. MAIN RESULTS AND THE ROLE OF CHANCE During 5-day culture, the highly pure trophoblasts were cultured on inserts coated with reconstituted basement membrane matrix . They exhibited a confluent polarized monolayer, with a modest TEER and a size-dependent apparent permeability coefficient (Papp) to fluorescently labeled compounds (MW ∼400-70 000 Da). The syncytialization progress was characterized by gradually increasing mRNA levels of fusogen genes and elevating hCG secretion. SEM analyses confirmed a confluent trophoblast layer with numerous microvilli, and TEM revealed a monolayer with tight junctions. Immunocytochemistry on the confluent trophoblasts showed positivity for the cell-cell adhesion molecule E-cadherin, the tight junction protein 1 (ZO-1) and the membrane proteins ATP-binding cassette transporter A1 (ABCA1) and glucose transporter 1 (GLUT1). Applying this model to study the bidirectional transport of a non-metabolizable glucose derivative indicated a carrier-mediated placental glucose transport mechanism with asymmetric kinetics. LIMITATIONS, REASONS FOR CAUTION The current study is only focused on primary trophoblast cells isolated from healthy placentas delivered at term. It remains to be evaluated whether this system can be extended to pathological trophoblasts isolated from diverse gestational diseases. WIDER IMPLICATIONS OF THE FINDINGS These findings confirmed the physiological properties of the newly developed human trophoblast barrier, which can be applied to study the exchange of endobiotics and xenobiotics between the maternal and fetal compartment, as well as intracellular metabolism, paracellular contributions and regulatory mechanisms influencing the vectorial transport of molecules. LARGE-SCALE DATA Not applicable. STUDY FUNDING AND COMPETING INTERESTS This study was supported by the Swiss National Center of Competence in Research, NCCR TransCure, University of Bern, Switzerland, and the Swiss National Science Foundation (grant no. 310030_149958, C.A.). All authors declare that their participation in the study did not involve factual or potential conflicts of interests.

MODELING THE DYNAMICS OF AIRWAY CONSTRICTION: EFFECTS OF AGONIST TRANSPORT AND BINDING

Recent advances have revealed that during exogenous airway challenge, airway diameters can not be adequately predicted by their initial diameters. Furthermore, airway diameters can also vary greatly in time on scales shorter than a breath. In order to better understand these phenomena, we developed a multiscale model which allows us to simulate aerosol challenge in the airways during ventilation. The model incorporates agonist-receptor binding kinetics to govern the temporal response of airway smooth muscle (ASM) contraction on individual airway segments, which together with airway wall mechanics, determines local airway caliber. Global agonist transport and deposition is coupled with pressure-driven flow, linking local airway constrictions with global flow dynamics. During the course of challenge, airway constriction alters the flow pattern, redistributing agonist to less constricted regions. This results in a negative feedback which may be a protective property of the normal lung. As a consequence, repetitive challenge can cause spatial constriction patterns to evolve in time, resulting in a loss of predictability of airway diameters. Additionally, the model offers new insight into several phenomena including the intra- and inter-breath dynamics of airway constriction throughout the tree structure.

Expression, localization, and functional model of cholesterol transporters in lactating and nonlactating mammary tissues of murine, bovine, and human origin

Members of the ATP-binding cassette (ABC) transporters play a pivotal role in cellular lipid efflux. To identify candidate cholesterol transporters implicated in lipid homeostasis and mammary gland (MG) physiology, we compared expression and localization of ABCA1, ABCG1, and ABCA7 and their regulatory genes in mammary tissues of different species during the pregnancy-lactation cycle. Murine and bovine mammary glands (MGs) were investigated during different functional stages. The abundance of mRNAs was determined by quantitative RT-PCR. Furthermore, transporter proteins were localized in murine, bovine, and human MGs by immunohistochemistry. In the murine MG, ABCA1 mRNA abundance was elevated during nonlactating compared with lactating stages, whereas ABCA7 and ABCA1 mRNA profiles were not altered. In the bovine MG, ABCA1, ABCG1, and ABCA7 mRNAs abundances were increased during nonlactating stages compared with lactation. Furthermore, associations between mRNA levels of transporters and their regulatory genes LXRalpha, PPARgamma, and SREBPs were found. ABCA1, ABCG1, and ABCA7 proteins were localized in glandular MG epithelial cells (MEC) during lactation, whereas during nonlactating stages, depending on species, the proteins showed distinct localization patterns in MEC and adipocytes. Our results demonstrate that ABCA1, ABCG1, and ABCA7 are differentially expressed between lactation and nonlactating stages and in association with regulatory genes. Combined expression and localization data suggest that the selected cholesterol transporters are universal MG transporters involved in transport and storage of cholesterol and in lipid homeostasis of MEC. Because of the species-specific expression patterns of transporters in mammary tissue, mechanisms of cholesterol homeostasis seem to be differentially regulated between species.

Interleukin-6 is an efficacious marker of axonal transport disruption during experimental glaucoma and stimulates neuritogenesis in cultured retinal ganglion cells

It is increasingly recognised that chronically activated glia contribute to the pathology of various neurodegenerative diseases, including glaucoma. One means by which this can occur is through the release of neurotoxic, proinflammatory factors. In the current study, we therefore investigated the spatio-temporal patterns of expression of three such cytokines, IL-1β, TNFα and IL-6, in a validated rat model of experimental glaucoma. First, only weak evidence was found for increased expression of IL-1β and TNFα following induction of ocular hypertension. Second, and much more striking, was that robust evidence was uncovered showing IL-6 to be synthesised by injured retinal ganglion cells following elevation of intraocular pressure and transported in an orthograde fashion along the nerve, accumulating at sites of axonal disruption in the optic nerve head. Verification that IL-6 represents a novel marker of disrupted axonal transport in this model was obtained by performing double labelling immunofluorescence with recognised markers of fast axonal transport. The stimulus for IL-6 synthesis and axonal transport during experimental glaucoma arose from axonal injury rather than ocular hypertension, as the response was identical after optic nerve crush and bilateral occlusion of the carotid arteries, each of which is independent of elevated intraocular pressure. Moreover, the response of IL-6 was not a generalised feature of the gp130 family of cytokines, as it was not mimicked by another family member, ciliary neurotrophic factor. Finally, further study suggested that IL-6 may be an early part of the endogenous regenerative response as the cytokine colocalised with growth-associated membrane phosphoprotein-43 in some putative regenerating axons, and potently stimulated neuritogenesis in retinal ganglion cells in culture, an effect that was additive to that of ciliary neurotrophic factor. These data comprise clear evidence that IL-6 is actively involved in the attempt of injured retinal ganglion cells to regenerate their axons.

Diffusion-driven transport in clayrock formations

Clay mineral-rich sedimentary formations are currently under investigation to evaluate their potential use as host formations for installation of deep underground disposal facilities for radioactive waste (e.g. Boom Clay (BE), Opalinus Clay (CH), Callovo-Oxfordian argillite (FR)). The ultimate safety of the corresponding repository concepts depends largely on the capacity of the host formation to limit the flux towards the biosphere of radionuclides (RN) contained in the waste to acceptably low levels. Data for diffusion-driven transfer in these formations shows extreme differences in the measured or modelled behaviour for various radionuclides, e. g. between halogen RN (Cl-36, I-129) and actinides (U-238,U-235, Np-237, Th-232, etc.), which result from major differences between RN of the effects on transport of two phenomena: diffusion and sorption. This paper describes recent research aimed at improving understanding of these two phenomena, focusing on the results of studies carried out during the EC Funmig IP on clayrocks from the above three formations and from the Boda formation (HU). Project results regarding phenomena governing water, cation and anion distribution and mobility in the pore volumes influenced by the negatively-charged surfaces of clay minerals show a convergence of the modelling results for behaviour at the molecular scale and descriptions based on electrical double layer models. Transport models exist which couple ion distribution relative to the clay-solution interface and differentiated diffusive characteristics. These codes are able to reproduce the main trends in behaviour observed experimentally, e.g. D-e(anion) < D-e(HTO) < D-e(cation) and D-e(anion) variations as a function of ionic strength and material density. These trends are also well-explained by models of transport through ideal porous matrices made up of a charged surface material. Experimental validation of these models is good as regards monovalent alkaline cations, in progress for divalent electrostatically-interacting cations (e.g. Sr2+) and still relatively poor for 'strongly sorbing', high K-d cations. Funmig results have clarified understanding of how clayrock mineral composition, and the corresponding organisation of mineral grain assemblages and their associated porosity, can affect mobile solute (anions, HTO) diffusion at different scales (mm to geological formation). In particular, advances made in the capacity to map clayrock mineral grain-porosity organisation at high resolution provide additional elements for understanding diffusion anisotropy and for relating diffusion characteristics measured at different scales. On the other hand, the results of studies focusing on evaluating the potential effects of heterogeneity on mobile species diffusion at the formation scale tend to show that there is a minimal effect when compared to a homogeneous property model. Finally, the results of a natural tracer-based study carried out on the Opalinus Clay formation increase confidence in the use of diffusion parameters measured on laboratory scale samples for predicting diffusion over geological time-space scales. Much effort was placed on improving understanding of coupled sorption-diffusion phenomena for sorbing cations in clayrocks. Results regarding sorption equilibrium in dispersed and compacted materials for weakly to moderately sorbing cations (Sr2+, Cs+, Co2+) tend to show that the same sorption model probably holds in both systems. It was not possible to demonstrate this for highly sorbing elements such as Eu(III) because of the extremely long times needed to reach equilibrium conditions, but there does not seem to be any clear reason why such elements should not have similar behaviour. Diffusion experiments carried out with Sr2+, Cs+ and Eu(III) on all of the clayrocks gave mixed results and tend to show that coupled diffusion-sorption migration is much more complex than expected, leading generally to greater mobility than that predicted by coupling a batch-determined K-d and Ficks law based on the diffusion behaviour of HTO. If the K-d measured on equivalent dispersed systems holds as was shown to be the case for Sr, Cs (and probably Co) for Opalinus Clay, these results indicate that these cations have a D-e value higher than HTO (up to a factor of 10 for Cs+). Results are as yet very limited for very moderate to strongly sorbing species (e.g. Co(II), Eu(III), Cu(II)) because of their very slow transfer characteristics. (C) 2011 Elsevier Ltd. All rights reserved.

Enhanced Atlantic subpolar gyre variability through baroclinic threshold in a coarse resolution model

Direct observations, satellite measurements and paleo records reveal strong variability in the Atlantic subpolar gyre on various time scales. Here we show that variations of comparable amplitude can only be simulated in a coupled climate model in the proximity of a dynamical threshold. The threshold and the associated dynamic response is due to a positive feedback involving increased salt transport in the subpolar gyre and enhanced deep convection in its centre. A series of sensitivity experiments is performed with a coarse resolution ocean general circulation model coupled to a statistical-dynamical atmosphere model which in itself does not produce atmospheric variability. To simulate the impact of atmospheric variability, the model system is perturbed with freshwater forcing of varying, but small amplitude and multi-decadal to centennial periodicities and observational variations in wind stress. While both freshwater and wind-stress-forcing have a small direct effect on the strength of the subpolar gyre, the magnitude of the gyre's response is strongly increased in the vicinity of the threshold. Our results indicate that baroclinic self-amplification in the North Atlantic ocean can play an important role in presently observed SPG variability and thereby North Atlantic climate variability on multi-decadal scales.

Effects of thoracic epidural anesthesia on hemodynamics and global oxygen transport in ovine endotoxemia

Besides providing effective analgesia, thoracic epidural anesthesia (TEA) has been shown to decrease perioperative morbidity and mortality. Because of its vasodilatory properties in association with the sympathetic blockade, however, TEA may potentially aggravate cardiovascular dysfunctions resulting from sepsis and systemic inflammatory response syndrome. The objective of the present study was to assess the effects of TEA on hemodynamics, global oxygen transport, and renal function in ovine endotoxemia. After a baseline measurement in healthy sheep (n = 18), Salmonella typhosa endotoxin was centrally infused at incremental doses to induce and maintain a hypotensive-hypodynamic circulation using an established protocol. The animals were then randomly assigned to one of two groups. In the treatment group, continuous TEA was initiated with 0.1 mL.kg of 0.125% bupivacaine at the onset of endotoxemia and maintained with 0.1 mL.kg.h. In the control group, the same amount of isotonic sodium chloride solution was injected through the epidural catheter. In the animals surviving the entire experiment (n = 7 per group), cardiac index and mean arterial pressure decreased in a dose-dependent manner during endotoxin infusion. In the TEA group, neither systemic hemodynamics nor global oxygen transport were impaired beyond the changes caused by endotoxemia itself. Urinary output was increased in the TEA group as compared with the control group (P < 0.05). In this model of endotoxic shock, TEA improved renal perfusion without affecting cardiopulmonary hemodynamics and global oxygen transport.

A multiple source model for 6 MV photon beam dose calculations using Monte Carlo

A multiple source model (MSM) for the 6 MV beam of a Varian Clinac 2300 C/D was developed by simulating radiation transport through the accelerator head for a set of square fields using the GEANT Monte Carlo (MC) code. The corresponding phase space (PS) data enabled the characterization of 12 sources representing the main components of the beam defining system. By parametrizing the source characteristics and by evaluating the dependence of the parameters on field size, it was possible to extend the validity of the model to arbitrary rectangular fields which include the central 3 x 3 cm2 field without additional precalculated PS data. Finally, a sampling procedure was developed in order to reproduce the PS data. To validate the MSM, the fluence, energy fluence and mean energy distributions determined from the original and the reproduced PS data were compared and showed very good agreement. In addition, the MC calculated primary energy spectrum was verified by an energy spectrum derived from transmission measurements. Comparisons of MC calculated depth dose curves and profiles, using original and PS data reproduced by the MSM, agree within 1% and 1 mm. Deviations from measured dose distributions are within 1.5% and 1 mm. However, the real beam leads to some larger deviations outside the geometrical beam area for large fields. Calculated output factors in 10 cm water depth agree within 1.5% with experimentally determined data. In conclusion, the MSM produces accurate PS data for MC photon dose calculations for the rectangular fields specified.

The role of ocean transport in the uptake of anthropogenic CO2

We compare modeled oceanic carbon uptake in response to pulse CO2 emissions using a suite of global ocean models and Earth system models. In response to a CO2 pulse emission of 590 Pg C (corresponding to an instantaneous doubling of atmospheric CO2 from 278 to 556 ppm), the fraction of CO2 emitted that is absorbed by the ocean is: 37±8%, 56±10%, and 81±4% (model mean ±2σ ) in year 30, 100, and 1000 after the emission pulse, respectively. Modeled oceanic uptake of pulse CO2 on timescales from decades to about a century is strongly correlated with simulated present-day uptake of chlorofluorocarbons (CFCs) and CO2 across all models, while the amount of pulse CO2 absorbed by the ocean from a century to a millennium is strongly correlated with modeled radiocarbon in the deep Southern and Pacific Ocean. However, restricting the analysis to models that are capable of reproducing observations within uncertainty, the correlation is generally much weaker. The rates of surface-to-deep ocean transport are determined for individual models from the instantaneous doubling CO2 simulations, and they are used to calculate oceanic CO2 uptake in response to pulse CO2 emissions of different sizes pulses of 1000 and 5000 Pg C. These results are compared with simulated oceanic uptake of CO2 by a number of models simulations with the coupling of climate-ocean carbon cycle and without it. This comparison demonstrates that the impact of different ocean transport rates across models on oceanic uptake of anthropogenic CO2 is of similar magnitude as that of climate-carbon cycle feedbacks in a single model, emphasizing the important role of ocean transport in the uptake of anthropogenic CO2.