Quantification of Short and Medium Range Order in Mixed Network Former Glasses of the System GeO2-NaPO3: A Combined NMR and X-ray Photoelectron Spectroscopy Study

Glasses in the system xGeO(2)-(1-x)NaPO3 (0 <= x <= 0.50) were prepared by conventional melting quenching and characterized by thermal analysis, Raman spectroscopy, X-ray photoelectron spectroscopy (XPS), and P-31 nuclear magnetic resonance (MAS NMR) techniques. The deconvolution of the latter spectra was aided by homonuclear J-resolved and refocused INADEQUATE techniques. The combined analyses of P-31 MAS NMR and O-1s XPS lineshapes, taking charge and mass balance considerations into account, yield the detailed quantitative speciations of the phosphorus, germanium, and oxygen atoms and their respective connectivities. An internally consistent description is possible without invoking the formation of higher-coordinated germanium species in these glasses, in agreement with experimental evidence in the literature. The structure can be regarded, to a first approximation, as a network consisting of P-(2) and P-(3) tetrahedra linked via four-coordinate germanium. As implied by the appearance of P-(3) units, there is a moderate extent of network modifier sharing between phosphate and germanate network formers, as expressed by the formal melt reaction P-(2) + Ge-(4) -> P-(3) + Ge-(3). The equilibrium constant of this reaction is estimated as K = 0.52 +/- 0.11, indicating a preferential attraction of network modifier by the phosphorus component. These conclusions are qualitatively supported by Raman spectroscopy as well as P-31{Na-23} and P-31{Na-23} rotational echo double resonance (REDOR) NMR results. The combined interpretation of O-1s XPS and P-31 MAS NMR spectra shows further that there are clear deviations from a random connectivity scenario: heteroatomic P-O-Ge linkages are favored over homoatomic P-O-P and Ge-O-Ge linkages.

Hydrodynamic analysis of particle collection efficiency: comparing downflow and upflow filtration

Models of the filtration phenomenon describe the mass balance in bed filtration in terms of particle removal mechanisms, and allow for the determination of global particle removal efficiencies. These models are defined in terms of the geometry and characteristic elements of granule collectors, particles and fluid, and also the composition of the balance of forces that act in the particle collector system. This work analyzes particles collection efficiency comparing downflow and upflow direct filtration, taking into account the contribution of the gravitational factor of the settling removal efficiency in future proposal of initial collection efficiency models for upflow filtration. A qualitative analysis is also made of the proposal for the collection efficiency models for particle removal in direct downflow and upflow filtration using a Computational Fluid Dynamics (CFD) tool. This analysis showed a strong influence of gravitational factor in initial collection efficiency (t = 0) of particles, as well as the reasons of their values to be smaller for upflow filtration in comparison with the downflow filtration.

Models of the Distribution of Persistent Organic Pollutants in the Marine Environment

Persistent organic pollutants (POPs) is a group of chemicals that are toxic, undergo long-range transport and accumulate in biota. Due to their persistency the distribution and recirculation in the environment often continues for a long period of time. Thereby they appear virtually everywhere within the biosphere, and poses a toxic stress to living organisms. In this thesis, attempts are made to contribute to the understanding of factors that influence the distribution of POPs with focus on processes in the marine environment. The bioavailability and the spatial distribution are central topics for the environmental risk management of POPs. In order to study these topics, various field studies were undertaken. To determine the bioavailable fraction of polychlorinated dibenzo-p-dioxins and dibenzofurans (PCDD/Fs), polychlorinated naphthalenes (PCNs), and polychlorinated biphenyls (PCBs) the aqueous dissolved phase were sampled and analysed. In the same samples, we also measured how much of these POPs were associated with suspended particles. Different models, which predicted the phase distribution of these POPs, were then evaluated. It was found that important water characteristics, which influenced the solid-water phase distribution of POPs, were particulate organic matter (POM), particulate soot (PSC), and dissolved organic matter (DOM). The bioavailable dissolved POP-phase in the water was lower when these sorbing phases were present. Furthermore, sediments were sampled and the spatial distribution of the POPs was examined. The results showed that the concentration of PCDD/Fs, and PCNs were better described using PSC- than using POM-content of the sediment. In parallel with these field studies, we synthesized knowledge of the processes affecting the distribution of POPs in a multimedia mass balance model. This model predicted concentrations of PCDD/Fs throughout our study area, the Grenlandsfjords in Norway, within factors of ten. This makes the model capable to validate the effect of suitable remedial actions in order to decrease the exposure of these POPs to biota in the Grenlandsfjords which was the aim of the project. Also, to evaluate the influence of eutrophication on the marine occurrence PCB data from the US Musselwatch and Benthic Surveillance Programs are examined in this thesis. The dry weight based concentrations of PCB in bivalves were found to correlate positively to the organic matter content of nearby sediments, and organic matter based concentrations of PCB in sediments were negatively correlated to the organic matter content of the sediment.

Glaciers and climate in northern Sweden during the 19th and 20th century

Our understanding of the climate of northern Sweden during the late Holocene is largely dependent on proxy-data series. These datasets remain spatially and temporally sparse and instrumental series are rare prior to the mid 19th century. Nevertheless, the glaciology and paleo-glaciology of the region has a strong potential significance for the exploration of climate change scenarios, past and future. The aim of this thesis is to investigate the 19th and 20th century climate in the northern Swedish mountain range. This provides a good opportunity to analyse the natural variability of the climate before the onset of the industrial epoch. Developing a temporal understanding of fluctuations in glacier front positions and glacier mass balance that is linked to a better understanding of their interaction and relative significance to climate is fundamental in the assessment of past climate. I have chosen to investigate previously unexplored temperature data from northern Sweden from between 1802 and 1860 and combined it with a temperature series from a synoptic station in Haparanda, which began operation in 1859, in order to create a reliable long temperature series for the period 1802 to 2002. I have also investigated two different glaciers, Pårteglaciären and Salajekna, which are located in different climatic environments. These glaciers have, from a Swedish perspective, long observational records. Furthermore, I have investigated a recurring jökulhlaup at the glacier Sälkaglaciären in order to analyse glacier-climate relationships with respect to the jökulhlaups. A number of datasets are presented, including: glacier frontal changes, in situ and photogrammetric mass balance data, in situ and satellite radar interferometry measurements of surface velocity, radar measurements, ice volume data and a temperature series. All these datasets are analysed in order to investigate the response of the glaciers to climatic stimuli, to attribute specific behaviour to particular climates and to analyse the 19th and 20th century glacier/climate relationships in northern Sweden. The 19th century was characterized by cold conditions in northern Sweden, particularly in winter. Significant changes in the amplitude of the annual temperature cycle are evident. Through the 19th century there is a marked decreasing trend in the amplitude of the data, suggesting a change towards a prevalence of maritime (westerly) air masses, something which has characterised the 20th century. The investigations on Salajekna support the conclusion that the major part of the 19th century was cold and dry. The 19th century advance of Salajekna was probably caused by colder climate in the late 18th and early 19th centuries, coupled with a weakening of the westerly airflow. The investigations on Pårteglaciären show that the glacier has a response time of ~200 years. It also suggests that there was a relatively high frequency of easterly winds providing the glacier with winter precipitation during the 19th century. Glaciers have very different response times and are sensitive to different climatic parameters. Glaciers in rather continental areas of the Subarctic and Arctic can have very long response times because of mass balance considerations and not primarily the glacier dynamics. This is of vital importance for analyzing Arctic and Subarctic glacier behaviour in a global change perspective. It is far from evident that the behaviour of the glacier fronts today reflects the present climate.

Mountain centered icefields in northern Scandinavia

Mountain centered glaciers have played a major role throughout the last three million years in the Scandinavian mountains. The climatic extremes, like the present warm interglacial or cold glacial maxima, are very short-lived compared to the periods of intermediate climate conditions, characterized by the persistence of mountain based glaciers and ice fields of regional size. These have persisted in the Scandinavian mountains for about 65% of the Quaternary. Mountain based glaciers thus had a profound impact on large-scale geomorphology, which is manifested in large-scale glacial landforms such as fjords, glacial lakes and U-shaped valleys in and close to the mountain range. Through a mapping of glacial landforms in the northern Scandinavian mountain range, in particular a striking set of lateral moraines, this thesis offers new insights into Weichselian stages predating the last glacial maximum. The aerial photograph mapping and field evidence yield evidence that these lateral moraines were overridden by glacier ice subsequent to their formation. The lateral moraines were dated using terrestrial cosmogenic nuclide techniques. Although the terrestrial cosmogenic nuclide signature of the moraines is inconclusive, an early Weichselian age is tentatively suggested through correlations with other landforms and stratigraphical archives in the region. The abundance and coherent spatial pattern of the lateral moraines also allow a spatial reconstruction of this ice field. The ice field was controlled by topography and had nunataks protruding also where it was thickest close to the elevation axis of the Scandinavian mountain range. Outlet glaciers discharged into the Norwegian fjords and major valleys in Sweden. The process by which mountain based glaciers grow into an ice sheet is a matter of debate. In this thesis, a feedback mechanism between debris on the ice surface and ice sheet growth is presented. In essence, the growth of glaciers and ice sheets may be accelerated by an abundance of debris in their ablation areas. This may occur when the debris cover on the glacier surface inhibits ablation, effectively increasing the glaciers mass balance. It is thus possible that a dirty ablation area may cause the glacier to advance further than a clean glacier under similar conditions. An ice free period of significant length allows soil production through weathering, frost shattering, and slope processes. As glaciers advance through this assemblage of sediments, significant amounts of debris end up on the surface due to both mass wastage and subglacial entrainment. Evidence that this chain of events may occur, is given by large expanses of hummocky moraine (local name Veiki moraine) in the northern Swedish lowlands. Because the Veiki moraine has been correlated with the first Weichselian advance following the Eemian, it implies a heavily debris charged ice sheet emanating from the mountain range and terminating in a stagnant fashion in the lowlands.

Dynamics of the cold surface layer of polythermal Storglaciären, Sweden

Polythermal glaciers, i.e. glaciers with a combination of ice at and below the freezing point, are widespread in arctic and subarctic environments. The polythermal structure has major implications for glacier hydrology, ice flow and glacial erosion. However, the interplay of factors governing its spatial and temporal variations such as net mass balance, ice advection and water content in the ice is poorly investigated and as yet not fully understood. This study deals with a thorough investigation of the polythermal regime on Storglaciären, northern Sweden, a small valley glacier with a cold surface layer in the ablation area. Extensive field work was performed including mapping of the cold surface layer using ground-penetrating radar, ice temperature measurements, mass balance and ice velocity measurements. Analyses of these data combined with numerical modelling were used specifically to investigate the spatial and temporal variability of the cold surface layer, the spatial distribution of the water content just below the cold surface layer transition, the effect of radar frequency on the detection of the surface layer, and the sensitivity of the cold surface layer to changes in forcing. A comparison between direct temperature measurements in boreholes and ground-penetrating surveys shows that the radar-inferred cold-temperate transition depth is within ±1 m from the melting point of ice at frequencies above ~300 MHz. At frequencies below ~155 MHz, the accuracy degrades because of reduced scattering efficiency that occurs when the scatterers become much smaller compared to the wavelength. The mapped spatial pattern of the englacial cold-temperate transition boundary is complex. This pattern reflects the observed spatial variation in net loss of ice at the surface by ablation and vertical advection of ice, which is suggested to provide the predominant forcing of the cold surface layer thickness pattern. This is further supported by thermomechanical modeling of the cold surface layer, which indicates high sensitivity of the cold surface layer thickness to changes in vertical advection rates. The water content is the least investigated quantity that is relevant for the thermal regime of glaciers, but also the most difficult to assess. Spatial variability of absolute water content in the temperate ice immediately below the cold surface layer on Storglaciären was determined by combining relative estimates of water content from ground-penetrating radar data with absolute determination from temperature measurements and the thermal boundary condition at the freezing front. These measurements indicate large-scale spatial variability in the water content, which seems to arise from variations in entrapment of water at the firn-ice transition. However, this variability cannot alone explain the spatial pattern in the thermal regime on Storglaciären. Repeated surveys of the cold surface layer show a 22% average thinning of the cold surface layer on Storglaciären between 1989 and 2001. Transient thermomechanical modeling results suggest that the cold surface layer adapts to new equilibrium conditions in only a few decades after a perturbation in the forcing is introduced. An increased winter air temperature since mid-1980s seems to be the cause of the observed thinning of the cold surface layer. Over the last decades, mass balance measurements indicate that the glacier has been close to a steady state. The quasi-steady state situation is also reflected in the vertical advection, which shows no significant changes during the last decades. Increased winter temperatures at the ice surface would result in a slow-down of the formation of cold ice at the base of the cold surface layer and lead to a larger imbalance between net loss of ice at the surface and freezing of temperate ice at the cold-temperate transition.

Geological evolution of the Monte Rosa

Eine detaillierte geochemische und geochronologische Studie an Gesteinen der Monte Rosa Decke (MR; Westalpen) wurde durchgeführt. Die MR wurde während der alpinen Orogenese zunächst eklogitfaziell und nachfolgend grünschieferfaziell überprägt.Eine detaillierte U-Pb geochronologische Studie an Zirkonen und Monaziten des MR Granits ergab ein Permisches Intrusionsalter (270 ± 4 Ma). Der MR Granit gehört zu den post-variszischen magmatischen Einheiten, welche die Instabilität der variszischen kontinentalen Kruste andeuten. Für die MR kann eine paläogeographische Position als Teil der 'Briançonnais-Schwelle' angenommen werden.Innerhalb des MR Granits treten Talk-Kyanit-Chloritoid-Gesteine ('Weißschiefer') auf. Diese stellen wesentliche Indikatoren für eine alpine Hochdruckmetamorphose in der MR dar. Massenbilanzberechnungen wurden durchgeführt, um den Massentransfer zu quantifizieren, welcher für die Bildung eines Weißschiefers aus einem granitischen Protolith notwendig ist. Ein Modell für die Entwicklung der Weißschiefer wurde entwickelt.Es wurde eine in-situ 40Ar/39Ar UV-Laser-Ablationsstudie an Hellglimmern der alpinen Mineralparagenese durchgeführt. Sie ergab eine heterogene Altersverteilung. Diese Alter können durch Glimmerrekristallisation unter relativ 'trockenen' hochdruckmetamorphen Bedingungen begleitet von partiellem Verlust von radiogenem Argon während der alpinen Metamorphose erklärt werden. Eine ähnlich komplexe Entwicklung mit partieller Homogenisierung des Isotopensystems kann in der Strontium-Isotopie beobachtet werden. Diese Isotopenstudien liefern Hinweise auf das Schließverhalten von Isotopensystemen unter hochdruckmetamorphen Bedingungen.

Quantification of deformation processes in the Torlesse accretionary wedge, New Zealand

Zusammenfassung:In dieser Studie werden Deformationsprozesse im mesozoischen Torlesse Akkretionskeil (Neuseeland) quantifiziert, um Aufschluß über die Dynamik in Akkretionskeilen zu erhalten. Absolute und relative Verformungsmessungen zeigen sowohl im lokalen als auch regionalen Maßstab eine stark heterogene Deformation des Torlesse Keils. Die regionale Deformation wurde mit Hilfe einer Tensordurchschnittsberechnung, unter Benutzung einzelner lokaler Verformungsdaten, als uniaxiale Verkürzung entlang einer subvertikalen, maximalen Verkürzungsachse charakterisiert. Absolute Verformungsmessungen an niedriggradigen Metasandsteinen belegen darüber hinaus durchschnittliche Volumenverluste von ca. 20% SiO2. Volumenveränderungen in tieferkrustalen Aufschlüssen wurden mittels einer geochemischen Massenbilanzanalyse abgeschätzt. Chemische Zusammensetzungen höhergradiger Zonen weichen je nach Grad der Volumenverformung von der Protolitzusammensetzung ab und zeigen somit Verluste von 15% SiO2 an. Da Speicherorte für das gelöste Material nicht bekannt sind, muss angenommen werden, dass das Material aus dem Keil abtransportiert wurde. Die Verformungsergebnisse geben weiterhin Aufschluß über den Grad der Kopplung zwischen Akkretionskeil und subduzierter Platte. Die ermittelten Scherwerte in den Gesteinen liegen deutlich unter den zu erwartenden Scherwerten, die mittels eines einfachen Modells berechnet wurden, das sowohl verschiedene Konvergenzgeschwindigkeiten als auch Exhumierungsraten berücksichtigt. Dies belegt, dass der Torlesse Keil stark von der subduzierten pazifischen Platte entkoppelt war und die Deformation hauptsächlich durch den Fluß der Sedimente in und aus dem Keil bestimmt wurde.

Sea Level Variations from Decades to Millennia: A few Case Studies

The main goals of this Ph.D. study are to investigate the regional and global geophysical components related to present polar ice melting and to provide independent cross validation checks of GIA models using both geophysical data detected by satellite mission, and geological observations from far field sites, in order to determine a lower and upper bound of uncertainty of GIA effect. The subject of this Thesis is the sea level change from decades to millennia scale. Within ice2sea collaboration, we developed a Fortran numerical code to analyze the local short-term sea level change and vertical deformation resulting from the loss of ice mass. This method is used to investigate polar regions: Greenland and Antarctica. We have used mass balance based on ICESat data for Greenland ice sheet and a plausible mass balance for Antarctic ice sheet. We have determined the regional and global fingerprint of sea level variations, vertical deformations of the solid surface of the Earth and variations of shape of the geoid for each ice source mentioned above. The coastal areas are affected by the long wavelength component of GIA process. Hence understanding the response of the Earth to loading is crucial in various contexts. Based on the hypothesis that Earth mantle materials obey to a linear rheology, and that the physical parameters of this rheology can be only characterized by their depth dependence, we investigate the Glacial Isostatic Effect upon the far field sites of Mediterranean area using an improved SELEN program. We presented new and revised observations for archaeological fish tanks located along the Tyrrhenian and Adriatic coast of Italy and new RSL for the SE Tunisia. Spatial and temporal variations of the Holocene sea levels studied in central Italy and Tunisia, provided important constraints on the melting history of the major ice sheets.

Structural and strain analysis in the Late Paleozoic coastal accretionary wedge of central Chile

Abstract In this study structural and finite strain data are used to explore the tectonic evolution and the exhumation history of the Chilean accretionary wedge. The Chilean accretionary wedge is part of a Late Paleozoic subduction complex that developed during subduction of the Pacific plate underneath South America. The wedge is commonly subdivided into a structurally lower Western Series and an upper Eastern Series. This study shows the progressive development of structures and finite strain from the least deformed rocks in the eastern part of the Eastern Series of the accretionary wedge to higher grade schist of the Western Series at the Pacific coast. Furthermore, this study reports finite-strain data to quantify the contribution of vertical ductile shortening to exhumation. Vertical ductile shortening is, together with erosion and normal faulting, a process that can aid the exhumation of high-pressure rocks. In the east, structures are characterized by upright chevron folds of sedimentary layering which are associated with a penetrative axial-plane foliation, S1. As the F1 folds became slightly overturned to the west, S1 was folded about recumbent open F2 folds and an S2 axial-plane foliation developed. Near the contact between the Western and Eastern Series S2 represents a prominent subhorizontal transposition foliation. Towards the structural deepest units in the west the transposition foliation became progressively flat lying. Finite-strain data as obtained by Rf/Phi and PDS analysis in metagreywacke and X-ray texture goniometry in phyllosilicate-rich rocks show a smooth and gradual increase in strain magnitude from east to west. There are no evidences for normal faulting or significant structural breaks across the contact of Eastern and Western Series. The progressive structural and strain evolution between both series can be interpreted to reflect a continuous change in the mode of accretion in the subduction wedge. Before ~320-290 Ma the rocks of the Eastern Series were frontally accreted to the Andean margin. Frontal accretion caused horizontal shortening and upright folds and axial-plane foliations developed. At ~320-290 Ma the mode of accretion changed and the rocks of the Western Series were underplated below the Andean margin. This basal accretion caused a major change in the flow field within the wedge and gave rise to vertical shortening and the development of the penetrative subhorizontal transposition foliation. To estimate the amount that vertical ductile shortening contributed to the exhumation of both units finite strain is measured. The tensor average of absolute finite strain yield Sx=1.24, Sy=0.82 and Sz=0.57 implying an average vertical shortening of ca. 43%, which was compensated by volume loss. The finite strain data of the PDS measurements allow to calculate an average volume loss of 41%. A mass balance approximates that most of the solved material stays in the wedge and is precipitated in quartz veins. The average of relative finite strain is Sx=1.65, Sy=0.89 and Sz=0.59 indicating greater vertical shortening in the structurally deeper units. A simple model which integrates velocity gradients along a vertical flow path with a steady-state wedge is used to estimate the contribution of deformation to ductile thinning of the overburden during exhumation. The results show that vertical ductile shortening contributed 15-20% to exhumation. As no large-scale normal faults have been mapped the remaining 80-85% of exhumation must be due to erosion.

Selective photo-oxidation of glucose

Biomass transformation into high-value chemicals has attracted attention according to the “green chemistry” principles. Low price and high availability make biomass one of the most interesting renewable resources as it provides the means to create sustainable alternatives to the oil-derived building blocks of the chemical industry In recent year, the need for alternative environmentally friendly routes to drive chemical reactions has in photocatalytic processes an interesting way to obtain valuable chemicals from various sources using the solar light as energy source. The purpose of this work was to use supported noble metal nanoparticles in the selective photo-oxidation of glucose through using visible light. Glucose was chosen as model molecule because it is the cheapest and the most common monosaccharide. Few studies about glucose photo oxidation have been conducted so far, and reaction mechanism is still not totally explained. The aim of this work was to systematically analyze and assess the impact of several parameters (eg. catalyst/substrate ratio, reaction time, effect of the solvent and light source) on the reaction pathway and to monitor the product distribution in order to draw a general reaction scheme for the photo oxidation of glucose under visible light. This study regards the reaction mechanism and the influence of several parameters, such as solvent, light power and substrate concentration. Furthermore, the work focuses on the influence of gold and silver nanoparticles and on the influence of metal loading. The glucose oxidation was monitored through the mass balance and the products selectivity. Reactions were evaluated in terms of glucose conversion, mass balance and selectivities towards arabinose and gluconic acid. In conclusion, this study is able to demonstrate that the photo oxidation of glucose under visible light is feasible; the full identification of the main products allows, for the first time, a comprehensive reaction mechanism scheme.

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.

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.

Modeling contaminant behavior in Lake Superior : a comparison of PCBs, PBDEs, and mercury

A mass‐balance model for Lake Superior was applied to polychlorinated biphenyls (PCBs), polybrominated diphenyl ethers (PBDEs), and mercury to determine the major routes of entry and the major mechanisms of loss from this ecosystem as well as the time required for each contaminant class to approach steady state. A two‐box model (water column, surface sediments) incorporating seasonally adjusted environmental parameters was used. Both numerical (forward Euler) and analytical solutions were employed and compared. For validation, the model was compared with current and historical concentrations and fluxes in the lake and sediments. Results for PCBs were similar to prior work showing that air‐water exchange is the most rapid input and loss process. The model indicates that mercury behaves similarly to a moderately‐chlorinated PCB, with air‐water exchange being a relatively rapid input and loss process. Modeled accumulation fluxes of PBDEs in sediments agreed with measured values reported in the literature. Wet deposition rates were about three times greater than dry particulate deposition rates for PBDEs. Gas deposition was an important process for tri‐ and tetra‐BDEs (BDEs 28 and 47), but not for higher‐brominated BDEs. Sediment burial was the dominant loss mechanism for most of the PBDE congeners while volatilization was still significant for tri‐ and tetra‐BDEs. Because volatilization is a relatively rapid loss process for both mercury and the most abundant PCBs (tri‐ through penta‐), the model predicts that similar times (from 2 ‐ 10 yr) are required for the compounds to approach steady state in the lake. The model predicts that if inputs of Hg(II) to the lake decrease in the future then concentrations of mercury in the lake will decrease at a rate similar to the historical decline in PCB concentrations following the ban on production and most uses in the U.S. In contrast, PBDEs are likely to respond more slowly if atmospheric concentrations are reduced in the future because loss by volatilization is a much slower process for PBDEs, leading to lesser overall loss rates for PBDEs in comparison to PCBs and mercury. Uncertainties in the chemical degradation rates and partitioning constants of PBDEs are the largest source of uncertainty in the modeled times to steady‐state for this class of chemicals. The modeled organic PBT loading rates are sensitive to uncertainties in scavenging efficiencies by rain and snow, dry deposition velocity, watershed runoff concentrations, and uncertainties in air‐water exchange such as the effect of atmospheric stability.