First day of an oil spill on the open sea: early mass transfers of hydrocarbons to air and water.

During the first hours after release of petroleum at sea, crude oil hydrocarbons partition rapidly into air and water. However, limited information is available about very early evaporation and dissolution processes. We report on the composition of the oil slick during the first day after a permitted, unrestrained 4.3 m(3) oil release conducted on the North Sea. Rapid mass transfers of volatile and soluble hydrocarbons were observed, with >50% of ≤C17 hydrocarbons disappearing within 25 h from this oil slick of <10 km(2) area and <10 μm thickness. For oil sheen, >50% losses of ≤C16 hydrocarbons were observed after 1 h. We developed a mass transfer model to describe the evolution of oil slick chemical composition and water column hydrocarbon concentrations. The model was parametrized based on environmental conditions and hydrocarbon partitioning properties estimated from comprehensive two-dimensional gas chromatography (GC×GC) retention data. The model correctly predicted the observed fractionation of petroleum hydrocarbons in the oil slick resulting from evaporation and dissolution. This is the first report on the broad-spectrum compositional changes in oil during the first day of a spill at the sea surface. Expected outcomes under other environmental conditions are discussed, as well as comparisons to other models.

Evaporation from a shallow water table: Diurnal dynamics of measured and simulated water and heat regime at the vicinity of a drying sand surface

Accurate estimates of water losses by evaporation from shallow water tables are important for hydrological, agricultural, and climatic purposes. An experiment was conducted in a weighing lysimeter to characterize the diurnal dynamics of evaporation under natural conditions. Sampling revealed a completely dry surface sand layer after 5 days of evaporation. Its thickness was <1 cm early in the morning, increasing to reach 4?5 cm in the evening. This evidence points out fundamental limitations of the approaches that assume hydraulic connectivity from the water table up to the surface, as well as those that suppose monotonic drying when unsteady conditions prevail. The computed vapor phase diffusion rates from the apparent drying front based on Fick's law failed to reproduce the measured cumulative evaporation during the sampling day. We propose that two processes rule natural evaporation resulting from daily fluctuations of climatic variables: (i) evaporation of water, stored during nighttime due to redistribution and vapor condensation, directly into the atmosphere from the soil surface during the early morning hours, that could be simulated using a mass transfer approach and (ii) subsurface evaporation limited by Fickian diffusion, afterward. For the conditions prevailing during the sampling day, the amount of water stored at the vicinity of the soil surface was 0.3 mm and was depleted before 11:00. Combining evaporation from the surface before 11:00 and subsurface evaporation limited by Fickian diffusion after that time, the agreement between the estimated and measured cumulative evaporation was significantly improved.

Large calcite and bulk-rock volume loss in metacarbonate xenoliths from the Qu,rigut massif (French Pyrenees)

Chemical mass transfer was quantified in a metacarbonate xenolith enclosed within the granodiorite of the Qu,rigut massif (Pyrenees, France). Mass balance calculations suggest a strong decrease of CaO, SrO and CO(2) contents (up to -90%), correlated with a decrease of modal calcite content as the contact is approached. Most other chemical elements behave immobile during metasomatism. They are therefore passively enriched. Only a small increase of SiO(2), Al(2)O(3) and Fe(2)O(3) contents occurs in the immediate vicinity of the contact. Hence, in this study, skarn formation is characterized by the lack of large chemical element influx from the granitoid protolith. A large decrease of the initial carbonate volume (up to -86%) resulted from a combination of decarbonation reactions and loss of CaO and CO(2). The resulting volume change has potentially important consequences for the interpretation of stable isotope profiles: the isotope alteration could have occured over greater distances than those observed today.

Formation of aluminosilicate-bearing quartz veins in the Simano nappe (Central Alps): structural, thermobarometric and oxygen isotope constraints

We combined structural analysis, thermobarometry and oxygen isotope geochemistry to constrain the evolution of kyanite and/or andalusite-bearing quartz veins from the amphibolite facies metapelites of the Simano nappe, in the Central Alps of Switzerland. The Simano nappe records a complex polyphase tectonic evolution associated with nappe stacking during Tertiary Alpine collision (D1). The second regional deformation phase (132) is responsible for the main penetrative schistosity and mineral lineation, and formed during top-to-the-north thrusting. During the next stage of deformation (D3) the aluminosilicate-bearing veins formed by crystallization in tension gashes, in tectonic shadows of boudins, as well as along shear bands associated with top-to-the-north shearing. D2 and D3 are coeval with the Early Miocene metamorphic peak, characterised by kyanite + staurolite + garnet + biotite assemblages in metapelites. The peak pressure (P) and temperature (T) conditions recorded are constrained by multiple-equilibrium thermobarometry at 630 +/- 20 degrees C and 8.5 +/- 1 kbar (similar to 27 km depth), which is in agreement with oxygen isotope thermometry indicating isotopic equilibration of quartz-kyanite pairs at 670 +/- 50 degrees C. Quartz-kyanite pairs from the aluminosilicate-bearing quartz veins yield equilibration temperatures of 645 +/- 20 degrees C, confirming that the veins formed under conditions near metamorphic peak. Quartz and kyanite from veins and the surrounding metapelites have comparable isotopic compositions. Local intergranular diffusion in the border of the veins controls the mass-transfer and the growth of the product assemblage, inducing local mobilization of SiO2 and Al2O3. Andalusite is absent from the host rocks, but it is common in quartz veins, where it often pseudomorphs kyanite. For andalusite to be stable at T-max, the pressure in the veins must have been substantially lower than lithostatic. An alternative explanation consistent with structural observations would be inheritance by andalusite of the kyanite isotopic signature during polymorphic transformation after the metamorphic peak.

Fluid-rock-strain feedbacks in extensional shear zones

La présence de fluide météorique synchrone à l'activité du détachement (Farmin, 2003 ; Mulch et al., 2007 ; Gébelin et al., 2011), implique que les zones de cisaillement sont des systèmes ouverts avec des cellules de convections à l'échelle crustale et un intense gradient géothermique au sein du détachement (Morrison et Anderson, 1998, Gottardi et al., 2011). De plus, les réactions métamorphiques liées à des infiltrations fluides dans les zones de cisaillement extensionnel peuvent influencer les paramètres rhéologiques du système (White and Knipe, 1978), et impliquer la localisation de la déformation dans la croûte. Dans ce manuscrit, deux zones de cisaillement infiltrées par des fluides météoriques sont étudiées, l'une étant largement quartzitique, et l'autre de nature granitique ; les relations entre déformation, fluides, et roches s'appuient sur des approches structurales, microstructurales, chimiques et isotopiques. L'étude du détachement du Columbia river (WA, USA) met en évidence que la déformation mylonitique se développe en un million d'années. La phase de cisaillement principal s'effectue à 365± 30°C d'après les compositions isotopiques en oxygène du quartz et de la muscovite. Ces minéraux atteignent l'équilibre isotopique lors de leur recristallisation dynamique contemporaine à la déformation. La zone de cisaillement enregistre une baisse de température, remplaçant le mécanisme de glissement par dislocation par celui de dissolution- précipitation dans les derniers stades de l'activité du détachement. La dynamique de circulation fluide bascule d'une circulation pervasive à chenalisée, ce qui engendre localement la rupture des équilibres d'échange isotopiques. La zone de cisaillement de Bitterroot (MT, USA) présente une zone mylonitique de 600m d'épaisseur, progressant des protomylonites aux ultramylonites. L'intensité de la localisation de la déformation se reflète directement sur l'hydratation des feldspaths, réaction métamorphique majeure dite de « rock softening ». Une étude sur roche totale indique des transferts de masse latéraux au sein des mylonites, et d'importantes pertes de volume dans les ultramylonites. La composition isotopique en hydrogène des phyllosilicates met en évidence la présence (1) d'une source magmatique/métamorphique originelle, caractérisée par les granodiorites ayant conservé leur foliation magmatique, jusqu'aux protomylonites, et (2) une source météorique qui tamponne les valeurs des phyllosilicates des fabriques mylonitiques jusqu'aux veines de quartz non-déformées. Les compositions isotopiques en oxygène des minéraux illustrent le tamponnement de la composition du fluide météorique par l'encaissant. Ce phénomène cesse lors du processus de chloritisation de la biotite, puisque les valeurs des chlorites sont extrêmement négatives (-10 per mil). La thermométrie isotopique indique une température d'équilibre isotopique de la granodiorite entre 600-500°C, entre 500-300°C dans les mylonites, et entre 300 et 200°C dans les fabriques cassantes (cataclasites et veines de quartz). Basé sur les résultats issus de ce travail, nous proposons un modèle général d'interactions fluide-roches-déformation dans les zones de détachements infiltrées par des fluides météoriques. Les zones de détachements évoluent rapidement (en quelques millions d'années) au travers de la transition fragile-ductile ; celle-ci étant partiellement contrôlée par l'effet thermique des circulations de fluide météoriques. Les systèmes de détachements sont des lieux où la déformation et les circulations fluides sont couplées ; évoluant rapidement vers une localisation de la déformation, et de ce fait, une exhumation efficace. - The presence of meteoric fluids synchronous with the activity of extensional detachment zones (Famin, 2004; Mulch et al., 2007; Gébelin et al., 2011) implies that extensional systems involve fluid convection at a crustal scale, which results in high geothermal gradients within active detachment zones (Morrison and Anderson, 1998, Gottardi et al., 2011). In addition, the metamorphic reactions related to fluid infiltration in extensional shear zones can influence the rheology of the system (White and Knipe, 1978) and ultimately how strain localizes in the crust. In this thesis, two shear zones that were permeated by meteoric fluids are studied, one quartzite-dominated, and the other of granitic composition; the relations between strain, fluid, and evolving rock composition are addressed using structural, microstructural, and chemical/isotopic measurements. The study of the Columbia River detachment that bounds the Kettle core complex (Washington, USA) demonstrates that the mylonitic fabrics in the 100 m thick quartzite- dominated detachment footwall developed within one million years. The main shearing stage occurred at 365 ± 30°C when oxygen isotopes of quartz and muscovite equilibrated owing to coeval deformation and dynamic recrystallization of these minerals. The detachment shear zone records a decrease in temperature, and dislocation creep during detachment shearing gave way to dissolution-precipitation and fracturing in the later stages of detachment activity. Fluid flow switched from pervasive to channelized, leading to isotopic disequilibrium between different minerals. The Bitterroot shear zone detachment (Montana, USA) developed a 600 m thick mylonite zone, with well-developed transitions from protomylonite to ultramylonite. The localization of deformation relates directly to the intensity of feldspar hydration, a major rock- softening metamorphic reaction. Bulk-rock analyses of the mylonitic series indicate lateral mass transfer in the mylonite (no volume change), and significant volume loss in ultramylonite. The hydrogen isotope composition of phyllosilicates shows (1) the presence of an initial magmatic/metamorphic source characterized by the granodiorite in which a magmatic, and gneissic (protomylonite) foliation developed, and (2) a meteoric source that buffers the values of phyllosilicates in mylonite, ultramylonite, cataclasite, and deformed and undeformed quartz veins. The mineral oxygen isotope compositions were buffered by the host-rock compositions until chloritization of biotite started; the chlorite oxygen isotope values are negative (-10 per mil). Isotope thermometry indicates a temperature of isotopic equilibrium of the granodiorite between 600-500°C, between 500-300°C in the mylonite, and between 300 and 200°C for brittle fabrics (cataclasite and quartz veins). Results from this work suggest a general model for fluid-rock-strain feedbacks in detachment systems that are permeated by meteoric fluids. Phyllosilicates have preserved in their hydrogen isotope values evidence for the interaction between rock and meteoric fluids during mylonite development. Fluid flow generates mass transfer along the tectonic anisotropy, and mylonites do not undergo significant volume change, except locally in ultramylonite zones. Hydration of detachment shear zones attends mechanical grain size reduction and enhances strain softening and localization. Self-exhuming detachment shear zones evolve rapidly (a few million years) through the transition from ductile to brittle, which is partly controlled by the thermal effect of circulating surface fluids. Detachment systems are zones in the crust where strain and fluid flow are coupled; these systems. evolve rapidly toward strain localization and therefore efficient exhumation.

A new green fluorescent protein-based bacterial biosensor for analysing phenanthrene fluxes.

The polycyclic aromatic hydrocarbon (PAH)-degrading strain Burkholderia sp. RP007 served as host strain for the design of a bacterial biosensor for the detection of phenanthrene. RP007 was transformed with a reporter plasmid containing a transcriptional fusion between the phnS putative promoter/operator region and the gene encoding the enhanced green fluorescent protein (GFP). The resulting bacterial biosensor--Burkholderia sp. strain RP037--produced significant amounts of GFP after batch incubation in the presence of phenanthrene crystals. Co-incubation with acetate did not disturb the phenanthrene-specific response but resulted in a homogenously responding population of cells. Active metabolism was required for induction with phenanthrene. The magnitude of GFP induction was influenced by physical parameters affecting the phenanthrene flux to the cells, such as the contact surface area between solid phenanthrene and the aqueous phase, addition of surfactant, and slow phenanthrene release from Model Polymer Release System beads or from a water-immiscible oil. These results strongly suggest that the bacterial biosensor can sense different phenanthrene fluxes while maintaining phenanthrene metabolism, thus acting as a genuine sensor for phenanthrene bioavailability. A relationship between GFP production and phenanthrene mass transfer is proposed.

Development of fluid conduits in the auriferous shear zones of the Hutti Gold Mine, India: evidence for spatially and temporally heterogeneous fluid flow

The gold mineralization of the Hutti Mine is hosted by nine parallel, N - S trending, steeply dipping, 2 - 10 m wide shear zones, that transect Archaean amphibolites. The shear zones were formed after peak metamorphism during retrograde ductile D, shearing in the lower amphibolite facies. They were reactivated in the lower to mid greenschist facies by brittle-ductile D-3 shearing and intense quartz veining. The development of a S-2-S-3 crenulation cleavage facilitates the discrimination between the two deformation events and contemporaneous alteration and gold mineralization. Ductile D, shearing is associated with a pervasively developed distal chlorite - sed cite alteration assemblage in the outer parts of the shear zones and the proximal biotite-plagioclase alteration in the center of the shear zones. D3 is characterized by development of the inner chlorite-K-feldspar alteration, which forms a centimeter-scale alteration halo surrounding the laminated quartz veins and replaces earlier biotite along S-3. The average size of the laminated vein systems is 30-50 m along strike as well as down-dip and 2-6 m in width. Mass balance calculations suggest strong metasomatic changes for the proximal biotite-plagioclase alteration yielding mass and volume increase of ca. 16% and 12%, respectively. The calculated mass and volume changes of the distal chlorite-sericite alteration (ca. 11%, ca. 8%) are lower. The decrease in 6180 values of the whole rock from around 7.5 parts per thousand for the host rocks to 6-7 parts per thousand for the distal chlorite-sericite and the proximal biotite-plagioclase alteration and around 5 parts per thousand for the inner chlorite-K-feldspar alteration suggests hydrothermal alteration during two-stage deformation and fluid flow. The ductile D-2 deformation in the lower amphibolite facies has provided grain scale porosities by microfracturing. The pervasive, steady-state fluid flow resulted in a disseminated style of gold-sulfide mineralization and a penetrative alteration of the host rocks. Alternating ductile and brittle D3 deformation during lower to mid greenschist facies conditions followed the fault-valve process. Ductile creep in the shear zones resulted in a low permeability environment leading to fluid pressure build-up. Strongly episodic fluid advection and mass transfer was controlled by repeated seismic fracturing during the formation of laminated quartz(-gold) veins. The limitation of quartz veins to the extent of earlier shear zones indicate the importance of preexisting anisotropies for fault-valve action and economic gold mineralization. (C) 2003 Elsevier B.V. All rights reserved.

Plutonium from above-ground nuclear tests in milk teeth: investigation of placental transfer in children born between 1951 and 1995 in Switzerland.

BACKGROUND: Occupational risks, the present nuclear threat, and the potential danger associated with nuclear power have raised concerns regarding the metabolism of plutonium in pregnant women. OBJECTIVE: We measured plutonium levels in the milk teeth of children born between 1951 and 1995 to assess the potential risk that plutonium incorporated by pregnant women might pose to the radiosensitive tissues of the fetus through placenta transfer. METHODS: We used milk teeth, whose enamel is formed during pregnancy, to investigate the transfer of plutonium from the mother's blood plasma to the fetus. We measured plutonium using sensitive sector field inductively coupled plasma mass spectrometry techniques. We compared our results with those of a previous study on strontium-90 ((90)Sr) released into the atmosphere after nuclear bomb tests. RESULTS: Results show that plutonium activity peaks in the milk teeth of children born about 10 years before the highest recorded levels of plutonium fallout. By contrast, (90)Sr, which is known to cross the placenta barrier, manifests differently in milk teeth, in accordance with (90)Sr fallout deposition as a function of time. CONCLUSIONS: These findings demonstrate that plutonium found in milk teeth is caused by fallout that was inhaled around the time the milk teeth were shed and not from any accumulation during pregnancy through placenta transfer. Thus, plutonium may not represent a radiologic risk for the radiosensitive tissues of the fetus.

Modelling changes in stable isotope compositions of minerals during net transfer reactions in a contact aureole: Wollastonite growth at the northern Hunter Mountain Batholith (Death Valley National Park, USA)

One of the world's largest wollastonite deposits was formed at the contact of the northern Hunter Mountain Batholith (California, USA) in Paleozoic sediments. Wollastonite occurs as zones of variable thickness surrounding layers or nodules of quartzite in limestones. A minimum formation temperature of 650 degrees C is estimated from isolated periclase-bearing lenses in that area. Contact metamorphism of siliceous carbonates has produced mineral assemblages that are consistent with heterogeneous, and partly limited infiltration of water-rich fluids, compatible with O-18/O-16 and C-13/C-12 isotopic patterns recorded in carbonates. Oxygen isotope compositions of wollastonites in the study area may also not require infiltration of large quantities of externally-derived fluids that were out of equilibrium with the rocks. 8180 values of wollastonite are high (14.8 parts per thousand to 25.0 parts per thousand; median: 19.7 parts per thousand) and close to those of the host limestone (19.7 parts per thousand to 28 parts per thousand; median: 24.9 parts per thousand) and quartz (18.0 parts per thousand. to 29.1 parts per thousand; median: 22.6 parts per thousand). Isotopic disequilibrium exists at quartz/wollastonite and wollastonite/calcite boundaries. Therefore, classical batch/Rayleigh fractionation models based on reactant and product equilibrium are not applicable to the wollastonite rims. An approach that relies on local instantaneous mass balance for the reactants, based on the wollastonite-forming reaction is suggested as an alternative way to model wollastonite reaction rims. This model reproduces many of the measured delta O-18 values of wollastonite reaction rims of the current study to within +/- 1 parts per thousand, even though the wollastonite compositions vary by almost 10 parts per thousand. (C) 2011 Elsevier B.V. All rights reserved.

Generic approach for the sensitive absolute quantification of large undigested peptides in plasma using a particular liquid chromatography-mass spectrometry setup.

A generic LC-MS approach for the absolute quantification of undigested peptides in plasma at mid-picomolar levels is described. Nine human peptides namely, brain natriuretic peptide (BNP), substance P (SubP), parathyroid hormone 1-34 (PTH), C-peptide, orexines A and B (Orex-A and -B), oxytocin (Oxy), gonadoliberin-1 (gonadothropin releasing-hormone or luteinizing hormone-releasing hormone, LHRH) and α-melanotropin (α-MSH) were targeted. Plasma samples were extracted via a 2-step procedure: protein precipitation using 1vol of acetonitrile followed by ultrafiltration of supernatants on membranes with a MW cut-off of 30 kDa. By applying a specific LC-MS setup, large volumes of filtrates (e.g., 2×750 μL) were injected and the peptides were trapped on a 1mm i.d.×10 mm length C8 column using a 10× on-line dilution. Then, the peptides were back-flushed and a second on-line dilution (2×) was applied during the transfer step. The refocalized peptides were resolved on a 0.3mm i.d. C18 analytical column. Extraction recovery, matrix effect and limits of detection were evaluated. Our comprehensive protocol demonstrates a simple and efficient sample preparation procedure followed by the analysis of peptides with limits of detection in the mid-picomolar range. This generic approach can be applied for the determination of most therapeutic peptides and possibly for endogenous peptides with latest state-of-the-art instruments.

Determination of Vasopressin and Desmopressin in urine by means of liquid chromatography coupled to quadrupole time-of-flight mass spectrometry for doping control purposes.

The anti-diuretic neurohypophysial hormone Vasopressin (Vp) and its synthetic analogue Desmopressin (Dp, 1-desamino-vasopressin) have received considerable attention from doping control authorities due to their impact on physiological blood parameters. Accordingly, the illicit use of Desmopressin in elite sport is sanctioned by the World Anti-Doping Agency (WADA) and the drug is classified as masking agent. Vp and Dp are small (8-9 amino acids) peptides administered orally as well as intranasally. Within the present study a method to determine Dp and Vp in urinary doping control samples by means of liquid chromatography coupled to quadrupole high resolution time-of-flight mass spectrometry was developed. After addition of Lys-Vasopressin as internal standard and efficient sample clean up with a mixed mode solid phase extraction (weak cation exchange), the samples were directly injected into the LC-MS system. The method was validated considering the parameters specificity, linearity, recovery (80-100%), accuracy, robustness, limit of detection/quantification (20/50 pg mL(-1)), precision (inter/intra-day<10%), ion suppression and stability. The analysis of administration study urine samples collected after a single intranasal or oral application of Dp yielded in detection windows for the unchanged target analyte for up to 20 h at concentrations between 50 and 600 pg mL(-1). Endogenous Vp was detected in concentrations of approximately 20-200 pg mL(-1) in spontaneous urine samples obtained from healthy volunteers. The general requirements of the developed method provide the characteristics for an easy transfer to other anti-doping laboratories and support closing another potential gap for cheating athletes.

Simultaneous determination of antidementia drugs in human plasma: Procedure transfer from HPLC-MS to UPLC-MS/MS.

A previously developed high performance liquid chromatography mass spectrometry (HPLC-MS) procedure for the simultaneous determination of antidementia drugs, including donepezil, galantamine, memantine, rivastigmine and its metabolite NAP 226-90, was transferred to an ultra performance liquid chromatography system coupled to a tandem mass spectrometer (UPLC-MS/MS). The drugs and their internal standards ([(2)H(7)]-donepezil, [(13)C,(2)H(3)]-galantamine, [(13)C(2),(2)H(6)]-memantine, [(2)H(6)]-rivastigmine) were extracted from 250μL human plasma by protein precipitation with acetonitrile. Chromatographic separation was achieved on a reverse phase column (BEH C18 2.1mm×50mm; 1.7μm) with a gradient elution of an ammonium acetate buffer at pH 9.3 and acetonitrile at a flow rate of 0.4mL/min and an overall run time of 4.5min. The analytes were detected on a tandem quadrupole mass spectrometer operated in positive electrospray ionization mode, and quantification was performed using multiple reaction monitoring. The method was validated according to the recommendations of international guidelines over a calibration range of 1-300ng/mL for donepezil, galantamine and memantine, and 0.2-50ng/mL for rivastimgine and NAP 226-90. The trueness (86-108%), repeatability (0.8-8.3%), intermediate precision (2.3-10.9%) and selectivity of the method were found to be satisfactory. Matrix effects variability was inferior to 15% for the analytes and inferior to 5% after correction by internal standards. A method comparison was performed with patients' samples showing similar results between the HPLC-MS and UPLC-MS/MS procedures. Thus, this validated UPLC-MS/MS method allows to reduce the required amount of plasma, to use a simplified sample preparation, and to obtain a higher sensitivity and specificity with a much shortened run-time.

Practical considerations for improving the productivity of mass spectrometry-based proteomics.

Mass spectrometry (MS) is currently the most sensitive and selective analytical technique for routine peptide and protein structure analysis. Top-down proteomics is based on tandem mass spectrometry (MS/ MS) of intact proteins, where multiply charged precursor ions are fragmented in the gas phase, typically by electron transfer or electron capture dissociation, to yield sequence-specific fragment ions. This approach is primarily used for the study of protein isoforms, including localization of post-translational modifications and identification of splice variants. Bottom-up proteomics is utilized for routine high-throughput protein identification and quantitation from complex biological samples. The proteins are first enzymatically digested into small (usually less than ca. 3 kDa) peptides, these are identified by MS or MS/MS, usually employing collisional activation techniques. To overcome the limitations of these approaches while combining their benefits, middle-down proteomics has recently emerged. Here, the proteins are digested into long (3-15 kDa) peptides via restricted proteolysis followed by the MS/MS analysis of the obtained digest. With advancements of high-resolution MS and allied techniques, routine implementation of the middle-down approach has been made possible. Herein, we present the liquid chromatography (LC)-MS/MS-based experimental design of our middle-down proteomic workflow coupled with post-LC supercharging.

Fate of the herbicide glyphosate and its metabolite AMPA in soils and their transfer to surface waters: A multi-scale approach in the Lavaux vineyards, western Switzerland

The use of herbicides in agriculture may lead to environmental problems, such as surface water pollution, with a potential risk for aquatic organisms. The herbicide glyphosate is the most used active ingredient in the world and in Switzerland. In the Lavaux vineyards it is nearly the only molecule applied. This work aimed at studying its fate in soils and its transfer to surface waters, using a multi-scale approach: from molecular (10-9 m) and microscopic scales (10-6 m), to macroscopic (m) and landscape ones (103 m). First of all, an analytical method was developed for the trace level quantification of this widely used herbicide and its main by-product, aminomethylphosphonic acid (AMPA). Due to their polar nature, their derivatization with 9-fluorenylmethyl chloroformate (FMOC-Cl) was done prior to their concentration and purification by solid phase extraction. They were then analyzed by ultra performance liquid chromatography coupled with tandem mass spectrometry (UPLC-MS/MS). The method was tested in different aqueous matrices with spiking tests and validated for the matrix effect correction in relevant environmental samples. Calibration curves established between 10 and 1000ng/l showed r2 values above 0.989, mean recoveries varied between 86 and 133% and limits of detection and quantification of the method were as low as 5 and 10ng/l respectively. At the parcel scale, two parcels of the Lavaux vineyard area, located near the Lutrive River at 6km to the east of Lausanne, were monitored to assess to which extent glyphosate and AMPA were retained in the soil or exported to surface waters. They were equipped at their bottom with porous ceramic cups and runoff collectors, which allowed retrieving water samples for the growing seasons 2010 and 2011. Results revealed that the mobility of glyphosate and AMPA in the unsaturated zone was likely driven by the precipitation regime and the soil characteristics, such as slope, porosity structure and layer permeability discrepancy. Elevated glyphosate and AMPA concentrations were measured at 60 and 80 cm depth at parcel bottoms, suggesting their infiltration in the upper parts of the parcels and the presence of preferential flow in the studied parcels. Indeed, the succession of rainy days induced the gradual saturation of the soil porosity, leading to rapid infiltration through macropores, as well as surface runoff formation. Furthermore, the presence of more impervious weathered marls at 100 cm depth induced throughflows, the importance of which for the lateral transport of the herbicide molecules was determined by the slope steepness. Important rainfall events (>10 mm/day) were clearly exporting molecules from the soil top layer, as indicated by important concentrations in runoff samples. A mass balance showed that total loss (10-20%) mainly occurred through surface runoff (96%) and, to a minor extent, by throughflows in soils (4%), with subsequent exfiltration to surface waters. Observations made in the Lutrive River revealed interesting details of glyphosate and AMPA dynamics in urbanized landscapes, such as the Lavaux vineyards. Indeed, besides their physical and chemical properties, herbicide dynamics at the catchment level strongly depend on application rates, precipitation regime, land use and also on the presence of drains or constructed channels. Elevated concentrations, up to 4970 ng/l, observed just after the application, confirmed the diffuse export of these compounds from the vineyard area by surface runoff during main rain events. From April to September 2011, a total load of 7.1 kg was calculated, with 85% coming from vineyards and minor urban sources and 15% from arable crops. Small vineyard surfaces could generate high concentrations of herbicides and contribute considerably to the total load calculated at the outlet, due to their steep slopes (~10%). The extrapolated total amount transferred yearly from the Lavaux vineyards to the Lake of Geneva was of 190kg. At the molecular scale, the possible involvement of dissolved organic matter (DOM) in glyphosate and copper transport was studied using UV/Vis fluorescence spectroscopy. Combined with parallel factor (PARAFAC) analysis, this technique allowed characterizing DOM of soil and surface water samples from the studied vineyard area. Glyphosate concentrations were linked to the fulvic-like spectroscopic signature of DOM in soil water samples, as well as to copper, suggesting the formation of ternary complexes. In surface water samples, its concentrations were also correlated to copper ones, but not in a significant way to the fulvic-like signature. Quenching experiments with standards confirmed field tendencies in the laboratory, with a stronger decrease in fluorescence intensity for fulvic-like fluorophore than for more aromatic ones. Lastly, based on maximum concentrations measured in the river, an environmental risk for these compounds was assessed, using laboratory tests and ecotoxicity data from the literature. In our case and with the methodology applied, the risk towards aquatic species was found negligible (RF<1).