Azo-propofols: photochromic potentiators of GABA(A) receptors

Shine and rise! GABA(A) receptors are ligand-gated chloride ion channels that respond to γ-aminobutyric acid (GABA), which is the major inhibitory neurotransmitter of the mammalian central nervous system. Azobenzene derivatives of propofol, such as compound 1 (see scheme), increase GABA-induced currents in the dark form and lose this property upon light exposure and thus function as photochromic potentiators. Compound 1 can be employed as a light-dependent general anesthetic in translucent tadpoles.

High-resolution computer simulation of the dynamics of isoelectric focusing using carrier ampholytes: focusing with concurrent electrophoretic mobilization is an isotachophoretic process

Focusing of four hemoglobins with concurrent electrophoretic mobilization was studied by computer simulation. A dynamic electrophoresis simulator was first used to provide a detailed description of focusing in a 100-carrier component, pH 6-8 gradient using phosphoric acid as anolyte and NaOH as catholyte. These results are compared to an identical simulation except that the catholyte contained both NaOH and NaCl. A stationary, steady-state distribution of carrier components and hemoglobins is produced in the first configuration. In the second, the chloride ion migrates into and through the separation space. It is shown that even under these conditions of chloride ion flux a pH gradient forms. All amphoteric species acquire a slight positive charge upon focusing and the whole pattern is mobilized towards the cathode. The cathodic gradient end is stable whereas the anodic end is gradually degrading due to the continuous accumulation of chloride. The data illustrate that the mobilization is a cationic isotachophoretic process with the sodium ion being the leading cation. The peak height of the hemoglobin zones decreases somewhat upon mobilization, but the zones retain a relatively sharp profile, thus facilitating detection. The electropherograms that would be produced by whole column imaging and by a single detector placed at different locations along the focusing column are presented and show that focusing can be commenced with NaCl present in the catholyte at the beginning of the experiment. However, this may require detector placement on the cathodic side of the catholyte/sample mixture interface.

Quantitative magnetic resonance imaging of enzymatically induced degradation of the nucleus pulposus of intervertebral discs

STUDY DESIGN: The structural integrity of the nucleus pulposus (NP) of intervertebral discs was targeted by enzyme-specific degradations to correlate their effects to the magnetic resonance (MR) signal. OBJECTIVE: To develop quantitative MR imaging as an accurate and noninvasive diagnostic tool to better understand and treat disc degeneration. SUMMARY OF BACKGROUND DATA: Quantitative MR analysis has been previously shown to reflect not only the disc matrix composition, but also the structural integrity of the disc matrix. Further work is required to identify the contribution of the structural integrity versus the matrix composition to the MR signal. METHODS: The bovine coccygeal NPs were injected with either enzyme or buffer, incubated at 37 degrees C as static, unloaded and closed 3-disc segments, and analyzed by a 1.5-Tesla MR scanner to measure MR parameters. RESULTS: Collagenase degradation of the NP significantly decreased the relaxation times, slightly decreased the magnetization transfer ratio, and slightly increased the apparent diffusion coefficient. Targeting the proteoglycan and/or hyaluronan integrity by trypsin and hyaluronidase did not significantly affect the MR parameters, except for an increase in the apparent diffusion coefficient of the disc after trypsin treatment. CONCLUSIONS: Our results demonstrate that changes in the structural integrity of matrix proteins can be assessed by quantitative MR.

Catch-up alveolarization in ex-preterm children: evidence from (3)He magnetic resonance

RATIONALE Histologic data from fatal cases suggest that extreme prematurity results in persisting alveolar damage. However, there is new evidence that human alveolarization might continue throughout childhood and could contribute to alveolar repair. OBJECTIVES To examine whether alveolar damage in extreme-preterm survivors persists into late childhood, we compared alveolar dimensions between schoolchildren born term and preterm, using hyperpolarized helium-3 magnetic resonance. METHODS We recruited schoolchildren aged 10-14 years stratified by gestational age at birth (weeks) to four groups: (1) term-born (37-42 wk; n = 61); (2) mild preterm (32-36 wk; n = 21); (3) extreme preterm (<32 wk, not oxygen dependent at 4 wk; n = 19); and (4) extreme preterm with chronic lung disease (<32 wk and oxygen dependent beyond 4 wk; n = 18). We measured lung function using spirometry and plethysmography. Apparent diffusion coefficient, a surrogate for average alveolar dimensions, was measured by helium-3 magnetic resonance. MEASUREMENTS AND MAIN RESULTS The two extreme preterm groups had a lower FEV1 (P = 0.017) compared with term-born and mild preterm children. Apparent diffusion coefficient was 0.092 cm(2)/second (95% confidence interval, 0.089-0.095) in the term group. Corresponding values were 0.096 (0.091-0.101), 0.090 (0085-0.095), and 0.089 (0.083-0.094) in the mild preterm and two extreme preterm groups, respectively, implying comparable alveolar dimensions across all groups. Results did not change after controlling for anthropometric variables and potential confounders. CONCLUSIONS Alveolar size at school age was similar in survivors of extreme prematurity and term-born children. Because extreme preterm birth is associated with deranged alveolar structure in infancy, the most likely explanation for our finding is catch-up alveolarization.

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.

Solute transport in crystalline rocks at Äspö — II: Blind predictions, inverse modelling and lessons learnt from test STT1

Based on the results from detailed structural and petrological characterisation and on up-scaled laboratory values for sorption and diffusion, blind predictions were made for the STT1 dipole tracer test performed in the Swedish A¨ spo¨ Hard Rock Laboratory. The tracers used were nonsorbing, such as uranine and tritiated water, weakly sorbing 22Na+, 85Sr2 +, 47Ca2 +and more strongly sorbing 86Rb+, 133Ba2 +, 137Cs+. Our model consists of two parts: (1) a flow part based on a 2D-streamtube formalism accounting for the natural background flow field and with an underlying homogeneous and isotropic transmissivity field and (2) a transport part in terms of the dual porosity medium approach which is linked to the flow part by the flow porosity. The calibration of the model was done using the data from one single uranine breakthrough (PDT3). The study clearly showed that matrix diffusion into a highly porous material, fault gouge, had to be included in our model evidenced by the characteristic shape of the breakthrough curve and in line with geological observations. After the disclosure of the measurements, it turned out that, in spite of the simplicity of our model, the prediction for the nonsorbing and weakly sorbing tracers was fairly good. The blind prediction for the more strongly sorbing tracers was in general less accurate. The reason for the good predictions is deemed to be the result of the choice of a model structure strongly based on geological observation. The breakthrough curves were inversely modelled to determine in situ values for the transport parameters and to draw consequences on the model structure applied. For good fits, only one additional fracture family in contact with cataclasite had to be taken into account, but no new transport mechanisms had to be invoked. The in situ values for the effective diffusion coefficient for fault gouge are a factor of 2–15 larger than the laboratory data. For cataclasite, both data sets have values comparable to laboratory data. The extracted Kd values for the weakly sorbing tracers are larger than Swedish laboratory data by a factor of 25–60, but agree within a factor of 3–5 for the more strongly sorbing nuclides. The reason for the inconsistency concerning Kds is the use of fresh granite in the laboratory studies, whereas tracers in the field experiments interact only with fracture fault gouge and to a lesser extent with cataclasite both being mineralogically very different (e.g. clay-bearing) from the intact wall rock.

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.

Insular and caudate lesions release abnormal yawning in stroke patients.

Abnormal yawning is an underappreciated phenomenon in patients with ischemic stroke. We aimed at identifying frequently affected core regions in the supratentorial brain of stroke patients with abnormal yawning and contributing to the anatomical network concept of yawning control. Ten patients with acute anterior circulation stroke and ≥3 yawns/15 min without obvious cause were analyzed. The NIH stroke scale (NIHSS), Glasgow Coma Scale (GCS), symptom onset, period with abnormal yawning, blood oxygen saturation, glucose, body temperature, blood pressure, heart rate, and modified Rankin scale (mRS) were assessed for all patients. MRI lesion maps were segmented on diffusion-weighted images, spatially normalized, and the extent of overlap between the different stroke patterns was determined. Correlations between the period with abnormal yawning and the apparent diffusion coefficient (ADC) in the overlapping regions, total stroke volume, NIHSS and mRS were performed. Periods in which patients presented with episodes of abnormal yawning lasted on average for 58 h. Average GCS, NIHSS, and mRS scores were 12.6, 11.6, and 3.5, respectively. Clinical parameters were within normal limits. Ischemic brain lesions overlapped in nine out of ten patients: in seven patients in the insula and in seven in the caudate nucleus. The decrease of the ADC within the lesions correlated with the period with abnormal yawing (r = -0.76, Bonferroni-corrected p = 0.02). The stroke lesion intensity of the common overlapping regions in the insula and the caudate nucleus correlates with the period with abnormal yawning. The insula might be the long sought-after brain region for serotonin-mediated yawning.

Water–montmorillonite systems: Neutron scattering and tracer throughdiffusion studies

Designs for deep geological respositories of nuclear waste include bentonite as a hydraulic and chemisorption buffer material to protect the biosphere from leakage of radionuclides. Bentonite is chosen because it is a cheap, naturally occurring material with the required properties. It consists essentially of montmorillonite, a swelling clay mineral. Upon contact with groundwater such clays can seal the repository by incorporating water in the interlayers of their crystalline structure. The intercalated water exhibits significantly different properties to bulk water in the surrounding interparticle pores, such as lower diffusion coefficients (González Sánchez et. al. 2008). This doctoral thesis presents water distribution and diffusion behavior on various time and space scales in montmorillonite. Experimental results are presented for Na- and Cs-montmorillonite samples with a range of bulk dry densities (0.8 to 1.7 g/cm3). The experimental methods employed were neutron scattering (backscattering, diffraction, time-of-flight), adsorption measurements (water, nitrogen) and tracer-through diffusion. For the tracer experiments the samples were fully saturated via the liquid phase under volume-constrained conditions. In contrast, for the neutron scattering experiments, the samples were hydrated via the vapor phase and subsequently compacted, leaving a significant fraction of interparticle pores unfilled with water. Owing to these differences in saturation, the water contents of the samples for neutron scattering were characterized by gravimetry whereas those for the tracer experiments were obtained from the bulk dry density. The amount of surface water in interlayer pores could be successfully discriminated from the amount of bulk-like water in interparticle pores in Na- and Csmontmorillonite using neutron spectroscopy. For the first time in the literature, the distribution of water between these two pore environments was deciphered as a function of gravimetric water content. The amount was compared to a geometrical estimation of the amount of interlayer and interparticle water determined by neutron diffraction and adsorption measurements. The relative abundances of the 1 to 4 molecular water layers in the interlayer were determined from the area ratios of the (001)-diffraction peaks. Depending on the characterization method, different fractions of surface water and interlayer water were obtained. Only surface and interlayer water exists in amontmorillonite with water contents up to 0.18 g/g according to spectroscopic measurements and up to 0.32 g/g according to geometrical estimations, respectively. At higher water contents, bulk-like and interparticle water also exists. The amounts increase monotonically, but not linearly, from zero to 0.33 g/g for bulk-like water and to 0.43 g/g for interparticle water. It was found that water most likely redistributes between the surface and interlayer sites during the spectroscopic measurements and therefore the reported fraction is relevant only below about -10 ºC (Anderson, 1967). The redistribution effect can explain the discrepancy in fractions between the methods. In a novel approach the fractions of water in different pore environments were treated as a fixed parameter to derive local diffusion coefficients for water from quasielastic neutron scattering data, in particular for samples with high water contents. Local diffusion coefficients were obtained for the 1 to 4 molecular water layers in the interlayer of 0.5·10–9, 0.9·10–9, 1.5·10–9 and 1.4·10–9 m²/s, respectively, taking account of the different water fractions (molecular water layer, bulk-like water). The diffusive transport of 22Na and HTO through Na-montmorillonite was measured on the laboratory experimental scale (i.e. cm, days) by tracer through-diffusion experiments. We confirmed that diffusion of HTO is independent of the ionic strength of the external solution in contact with the clay sample but dependent on the bulk dry density. In contrast, the diffusion of 22Na was found to depend on both the ionic strength of the pore solution and on the bulk dry density. The ratio of the pore and surface diffusion could be experimentally determined for 22Na from the dependence of the diffusion coefficient on the ionic strength. Activation energies were derived from the temperaturedependent diffusion coefficients via the Arrhenius relation. In samples with high bulk dry density the activation energies are slightly higher than those of bulk water whereas in low density samples they are lower. The activation energies as a function of ionic strengths of the pore solutions are similar for 22Na and HTO. The facts that (i) the slope of the logarithmic effective diffusion coefficients as a function of the logarithmic ionic strength is less than unity for low bulk dry densities and (ii) two water populations can be observed for high gravimetric water contents (low bulk dry densities) support the interlayer and interparticle porosity model proposed by Glaus et al. (2007), Bourg et al. (2006, 2007) and Gimmi and Kosakowski (2011).

Dating polygenetic metamorphic assemblages along a transect across the Western Alps

Multichronometric analyses were performed on samples from a transect in the French-Italian Western Alps crossing nappes derived from the Briançonnais terrane and the Piemonte-Liguria Ocean, in an endeavour to constrain the high-pressure (HP) metamorphism and the retrogression history. 12 samples of white mica were analysed by 39Ar-40Ar stepwise heating, complemented by 2 samples from the Monte Rosa 100 km to the NE and also attributed to the Briançonnais terrane. One Sm-Nd and three Lu-Hf garnet ages from eclogites were also obtained. White mica ages decrease from ca. 300 Ma in the westernmost samples (Zone Houillère), reaching ca. 300 °C during Alpine metamorphism, to < 48 Ma in the internal units to the East, which reached ca. 500 °C during Alpine orogeny. The conventional “thermochronological” interpretation postulates Cretaceous Eo-Alpine HP metamorphism and younger “cooling ages” in the higher-temperature samples. However, Eocene Lu-Hf and Sm-Nd ages from the same samples cannot be interpreted as post-metamorphic cooling ages, which makes a Cretaceous eclogitization untenable. The age date from this transect require instead to replace conventional “thermochronology” by an approach combining age dating with detailed geochemical, petrological and microstructural investigations. Petrology reveals important mineralogical differences along the transect. Samples from the Zone Houillère mostly contain detrital mica. White mica with Si > 6.45 atoms per formula unit becomes more abundant eastward. Across the whole traverse, HP phengitic mica forms the D1 foliation. Syn-D2 mica is Si-poorer and associated with nappe stacking, exhumation, and hydrous retrogression under greenschist facies conditions. D1 phengite is very often corroded, overgrown or intergrown by syn-D2 muscovite. Most importantly, syn-D2 recrystallization is not limited to S2 schistosity domains; microchemical fingerprinting shows that it also can form pseudomorphs after crystals that could be mistaken to have formed during D1 based on microstructural arguments alone. Thereby the Cl concentration in white mica is a useful discriminator, since D2 retrogression was associated with a less saline fluid than eclogitization. Once the petrological stage is set, geochronology is straightforward. All samples contain mixtures of detrital, syn-D1 and syn-D2 mica, and retrogression phases (D3) in greatly varying proportions according to local pressure-temperature-fluid activity-deformation conditions. The correlation of age vs. Cl/K clearly identifies 47 ± 1 Ma as the age of formation of syn-D1 mica along the entire transect, including the Monte Rosa nappe samples. The inferred age of the greenschist-facies low-Si syn-D2 mica generation ranges within 39-43 Ma, with local variations. Coexistence of D1 and D2 ages, and the constancy of non-reset D1 ages along the entire transect, are strong evidence that the D1 white mica ages are very close to formation ages. Volume diffusion of Ar in white mica (activation energy E = 250 kJ/mol; pressure-adjusted diffusion coefficient D’0 < 0.03 cm2 s-1) has a subordinate effect on mineral ages compared to both prograde and retrograde recrystallization in most samples. Eocene Lu-Hf and Sm-Nd garnet ages are prograde and predate the HP peak.

A novel open channel blocker of GABA-A receptors

GABA-A receptors are chloride ion channels composed of five subunits, mediating fast synaptic and tonic inhibition in the mammalian brain. 19 different subunit isoforms have been identified, with the major receptor type in mammalian adult brain consisting of α1, β2, and γ2 subunits. GABA-A receptors are the target of numerous sedating and anxiolytic drugs such as benzodiazepines. The currently known endogenous ligands are GABA, neurosteroids and the endocannabinoid 2- arachidonoyl glycerol (2-AG). The pharmacological properties of this chloride ion channel strictly depend on receptor subunit composition and arrangement. GABA-A receptors bind and are inhibited by epileptogenic agents such as picrotoxin, and cyclodiene insecticides such as dieldrin. We screened aromatic monovalent anions with five-fold symmetry for inhibition of GABA-A receptors. One of the anions, PCCPinhibited currents elicited by GABA with comparable potency as picrotoxin. This inhibition showed all characteristics of an open channel block. The GABA-A receptor ion channel is lined by residues from the M2 membrane-spanning segment. To identify important residues of the pore involved in the interaction with the blocking molecules PCCP-, a mutation scan was performed in combination with subsequent analysis of the expressed mutant proteins using electrophysiological techniques. In a second project we characterised a light-switchable modulator of GABA-A receptors based on propofol. It was my responsibility to investigate the switching kinetics in patch clamp experiments. After its discovery in 1980, propofol has become the most widely used intravenous general anaesthetic. It is commonly accepted that the anaesthesia induced by this unusually lipophilic drug mostly results from potentiation of GABA induced currents. While GABA-A receptors respond to a variety of ligands, they are normally not sensitive towards light. This light sensitivity could be indirectly achieved by using modulators that can be optically switched between an active and an inactive form. We tested an azobenzene derivative of propofol where an aryldiazene unit is directly coupled to the pharmacophore. This molecule was termed azopropofol (AP2). The effect of AP2 on Cl- currents was investigated with electrophysiological techniques using α1β2γ2 GABA-A receptors expressed in Xenopus oocytes and HEK-cells. In the third project we wanted to investigate the functional role of GABA-A receptors in the liver, and their possible involvement in cell proliferation. GABA-A receptors are also found in a wide range of peripheral tissues, including parts of the peripheral nervous system and non-neural tissues such as smooth muscle, the female reproductive system, liver and several cancer tissues. However their precise function in non neuronal or cancerous cells is still unknown. For this purpose we investigated expression, localization and function of the hepatocytes GABA-A receptors in model cell lines and healthy and cancerous hepatocytes.

A pentasymmetric open channel blocker for Cys-loop receptor channels.

γ-Aminobutyric acid type A receptors (GABAA receptors) are chloride ion channels composed of five subunits, mediating fast synaptic and tonic inhibition in the mammalian brain. These receptors show near five-fold symmetry that is most pronounced in the second trans-membrane domain M2 lining the Cl- ion channel. To take advantage of this inherent symmetry, we screened a variety of aromatic anions with matched symmetry and found an inhibitor, pentacyanocyclopentdienyl anion (PCCP-) that exhibited all characteristics of an open channel blocker. Inhibition was strongly dependent on the membrane potential. Through mutagenesis and covalent modification, we identified the region α1V256-α1T261 in the rat recombinant GABAA receptor to be important for PCCP- action. Introduction of positive charges into M2 increased the affinity for PCCP- while PCCP- prevented the access of a positively charged molecule into M2. Interestingly, other anion selective cys-loop receptors were also inhibited by PCCP-, among them the Drosophila RDL GABAA receptor carrying an insecticide resistance mutation, suggesting that PCCP- could serve as an insecticide.

Implementation and evaluation of permeability-porosity and tortuosity-porosity relationships linked to mineral dissolution-precipitation

Changes of porosity, permeability, and tortuosity due to physical and geochemical processes are of vital importance for a variety of hydrogeological systems, including passive treatment facilities for contaminated groundwater, engineered barrier systems (EBS), and host rocks for high-level nuclear waste (HLW) repositories. Due to the nonlinear nature and chemical complexity of the problem, in most cases, it is impossible to verify reactive transport codes analytically, and code intercomparisons are the most suitable method to assess code capabilities and model performance. This paper summarizes model intercomparisons for six hypothetical scenarios with generally increasing geochemical or physical complexity using the reactive transport codes CrunchFlow, HP1, MIN3P, PFlotran, and TOUGHREACT. Benchmark problems include the enhancement of porosity and permeability through mineral dissolution, as well as near complete clogging due to localized mineral precipitation, leading to reduction of permeability and tortuosity. Processes considered in the benchmark simulations are advective-dispersive transport in saturated media, kinetically controlled mineral dissolution-precipitation, and aqueous complexation. Porosity changes are induced by mineral dissolution-precipitation reactions, and the Carman-Kozeny relationship is used to describe changes in permeability as a function of porosity. Archie’s law is used to update the tortuosity and the pore diffusion coefficient as a function of porosity. Results demonstrate that, generally, good agreement is reached amongst the computer models despite significant differences in model formulations. Some differences are observed, in particular for the more complex scenarios involving clogging; however, these differences do not affect the interpretation of system behavior and evolution.

Reconstruction of in-situ porosity and porewater compositions of low-permeability crystalline rocks: Magnitude of artefacts induced by drilling and sample recovery

Geological site characterisation programmes typically rely on drill cores for direct information on subsurface rocks. However, porosity, transport properties and porewater composition measured on drill cores can deviate from in-situ values due to two main artefacts caused by drilling and sample recovery: (1) mechanical disruption that increases porosity and (2) contamination of the porewater by drilling fluid. We investigated the effect and magnitude of these perturbations on large drill core samples (12–20 cm long, 5 cmdiameter) of high-grade, granitic gneisses obtained from 350 to 600 m depth in a borehole on Olkiluoto Island (SW Finland). The drilling fluid was traced with sodium–iodide. By combining out-diffusion experiments, gravimetry, UV-microscopy and iodide mass balance calculations, we successfully quantified the magnitudes of the artefacts: 2–6% increase in porosity relative to the bulk connected porosity and 0.9 to 8.9 vol.% contamination by drilling fluid. The spatial distribution of the drilling-induced perturbations was revealed by numerical simulations of 2D diffusion matched to the experimental data. This showed that the rims of the samples have a mechanically disrupted zone 0.04 to 0.22 cm wide, characterised by faster transport properties compared to the undisturbed centre (1.8 to 7.7 times higher pore diffusion coefficient). Chemical contamination was shown to affect an even wider zone in all samples, ranging from 0.15 to 0.60 cm, inwhich iodide enrichmentwas up to 180 mg/kgwater, compared to 0.5 mg/kgwater in the uncontaminated centre. For all samples in the present case study, it turned out that the magnitude of the artefacts caused by drilling and sample recovery is so small that no correction is required for their effects. Therefore, the standard laboratory measurements of porosity, transport properties and porewater composition can be taken as valid in-situ estimates. However, it is clear that the magnitudes strongly depend on site- and drilling-specific factors and therefore our results cannot be transferred simply to other locations. We recommend the approach presented in this study as a route to obtain reliable values in future drilling campaigns aimed at characterising in-situ bedrock properties.

Breastfeeding, lung volumes and alveolar size at school-age.

BACKGROUND Previous studies found larger lung volumes at school-age in formerly breastfed children, with some studies suggesting an effect modification by maternal asthma. We wanted to explore this further in children who had undergone extensive lung function testing. The current study aimed to assess whether breastfeeding was associated with larger lung volumes and, if so, whether all compartments were affected. We also assessed association of breastfeeding with apparent diffusion coefficient (ADC), which measures freedom of gas diffusion in alveolar-acinar compartments and is a surrogate of alveolar dimensions. Additionally, we assessed whether these effects were modified by maternal asthma. METHODS We analysed data from 111 children and young adults aged 11-21 years, who had participated in detailed lung function testing, including spirometry, plethysmography and measurement of ADC of (3)Helium ((3)He) by MR. Information on breastfeeding came from questionnaires applied in early childhood (age 1-4 years). We determined the association between breastfeeding and these measurements using linear regression, controlling for potential confounders. RESULTS We did not find significant evidence for an association between duration of breastfeeding and lung volumes or alveolar dimensions in the entire sample. In breastfed children of mothers with asthma, we observed larger lung volumes and larger average alveolar size than in non-breastfed children, but the differences did not reach significance levels. CONCLUSIONS Confirmation of effects of breastfeeding on lung volumes would have important implications for public health. Further investigations with larger sample sizes are warranted.