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).

Organic matter processing and soil evolution in a braided river system

Traditionally, braided river research has considered flow, sediment transport processes and, recently, vegetation dynamics in relation to river morphodynamics. However, if considering the development of woody vegetated patches over a time scale of decades, we must consider the extent to which soil forming processes, particularly related to soil organic matter, impact the alluvial geomorphic-vegetation system. Here we quantify the soil organic matter processing (humification) that occurs on young alluvial landforms. We sampled different geomorphic units, ranging from the active river channel to established river terraces in a braided river system. For each geomorphic unit, soil pits were used to sample sediment/soil layers that were analysed in terms of grain size (<2mm) and organic matter quantity and quality (RockEval method). A principal components analysis was used to identify patterns in the dataset. Results suggest that during the succession from bare river gravels to a terrace soil, there is a transition from small amounts of external organic matter supply provided by sedimentation processes (e.g. organic matter transported in suspension and deposited on bars), to large amounts of autogenic in situ organic matter production due to plant colonisation. This appears to change the time scale and pathways of alluvial succession (bio-geomorphic succession). However, this process is complicated by: the ongoing possibility of local sedimentation, which can serve to isolate surface layers via aggradation from the exogenic supply; and erosion which tends to create fresh deposits upon which organic matter processing must re-start. The result is a complex pattern of organic matter states as well as a general lack of any clear chronosequence within the active river corridor. This state reflects the continual battle between deposition events that can isolate organic matter from the surface, erosion events that can destroy accumulating organic matter and the early ecosystem processes necessary to assist the co-evolution of soil and vegetation. A key question emerges over the extent to which the fresh organic matter deposited in the active zone is capable of significantly transforming the local geochemical environment sufficiently to accelerate soil development.

Combined effects of endurance training and dietary unsaturated fatty acids on physical performance, fat oxidation and insulin sensitivity.

Endurance training improves exercise performance and insulin sensitivity, and these effects may be in part mediated by an enhanced fat oxidation. Since n-3 and n-9 unsaturated fatty acids may also increase fat oxidation, we hypothesised that a diet enriched in these fatty acids may enhance the effects of endurance training on exercise performance, insulin sensitivity and fat oxidation. To assess this hypothesis, sixteen normal-weight sedentary male subjects were randomly assigned to an isoenergetic diet enriched with fish and olive oils (unsaturated fatty acid group (UFA): 52 % carbohydrates, 34 % fat (12 % SFA, 12 % MUFA, 5 % PUFA), 14 % protein), or a control diet (control group (CON): 62 % carbohydrates, 24 % fat (12 % SFA, 6 % MUFA, 2 % PUFA), 14 % protein) and underwent a 10 d gradual endurance training protocol. Exercise performance was evaluated by measuring VO2max and the time to exhaustion during a cycling exercise at 80 % VO2max; glucose homeostasis was assessed after ingestion of a test meal. Fat oxidation was assessed by indirect calorimetry at rest and during an exercise at 50 % VO2max. Training significantly increased time to exhaustion, but not VO2max, and lowered incremental insulin area under the curve after the test meal, indicating improved insulin sensitivity. Those effects were, however, of similar magnitude in UFA and CON. Fat oxidation tended to increase in UFA, but not in CON. This difference was, however, not significant. It is concluded that a diet enriched with fish- and olive oil does not substantially enhance the effects of a short-term endurance training protocol in healthy young subjects.

On the diurnal soil water content dynamics during evaporation using dielectric methods

The water content dynamics in the upper soil surface during evaporation is a key element in land-atmosphere exchanges. Previous experimental studies have suggested that the soil water content increases at the depth of 5 to 15 cm below the soil surface during evapo- ration, while the layer in the immediate vicinity of the soil surface is drying. In this study, the dynamics of water content profiles exposed to solar radiative forcing was monitored at a high temporal resolution using dielectric methods both in the presence and absence of evaporation. A 4-d comparison of reported moisture content in coarse sand in covered and uncovered buckets using a commercial dielectric-based probe (70 MHz ECH2O-5TE, Decagon Devices, Pullman, WA) and the standard 1-GHz time domain reflectometry method. Both sensors reported a positive correlation between temperature and water content in the 5- to 10-cm depth, most pronounced in the morning during heating and in the afternoon during cooling. Such positive correlation might have a physical origin induced by evaporation at the surface and redistribution due to liquid water fluxes resulting from the temperature- gradient dynamics within the sand profile at those depths. Our experimental data suggest that the combined effect of surface evaporation and temperature-gradient dynamics should be considered to analyze experimental soil water profiles. Additional effects related to the frequency of operation and to protocols for temperature compensation of the dielectric sensors may also affect the probes' response during large temperature changes.

Robust Nonlinear Mapping of Soil Contamination Using Support Regression

Spatial data analysis mapping and visualization is of great importance in various fields: environment, pollution, natural hazards and risks, epidemiology, spatial econometrics, etc. A basic task of spatial mapping is to make predictions based on some empirical data (measurements). A number of state-of-the-art methods can be used for the task: deterministic interpolations, methods of geostatistics: the family of kriging estimators (Deutsch and Journel, 1997), machine learning algorithms such as artificial neural networks (ANN) of different architectures, hybrid ANN-geostatistics models (Kanevski and Maignan, 2004; Kanevski et al., 1996), etc. All the methods mentioned above can be used for solving the problem of spatial data mapping. Environmental empirical data are always contaminated/corrupted by noise, and often with noise of unknown nature. That's one of the reasons why deterministic models can be inconsistent, since they treat the measurements as values of some unknown function that should be interpolated. Kriging estimators treat the measurements as the realization of some spatial randomn process. To obtain the estimation with kriging one has to model the spatial structure of the data: spatial correlation function or (semi-)variogram. This task can be complicated if there is not sufficient number of measurements and variogram is sensitive to outliers and extremes. ANN is a powerful tool, but it also suffers from the number of reasons. of a special type ? multiplayer perceptrons ? are often used as a detrending tool in hybrid (ANN+geostatistics) models (Kanevski and Maignank, 2004). Therefore, development and adaptation of the method that would be nonlinear and robust to noise in measurements, would deal with the small empirical datasets and which has solid mathematical background is of great importance. The present paper deals with such model, based on Statistical Learning Theory (SLT) - Support Vector Regression. SLT is a general mathematical framework devoted to the problem of estimation of the dependencies from empirical data (Hastie et al, 2004; Vapnik, 1998). SLT models for classification - Support Vector Machines - have shown good results on different machine learning tasks. The results of SVM classification of spatial data are also promising (Kanevski et al, 2002). The properties of SVM for regression - Support Vector Regression (SVR) are less studied. First results of the application of SVR for spatial mapping of physical quantities were obtained by the authorsin for mapping of medium porosity (Kanevski et al, 1999), and for mapping of radioactively contaminated territories (Kanevski and Canu, 2000). The present paper is devoted to further understanding of the properties of SVR model for spatial data analysis and mapping. Detailed description of the SVR theory can be found in (Cristianini and Shawe-Taylor, 2000; Smola, 1996) and basic equations for the nonlinear modeling are given in section 2. Section 3 discusses the application of SVR for spatial data mapping on the real case study - soil pollution by Cs137 radionuclide. Section 4 discusses the properties of the modelapplied to noised data or data with outliers.

Bayesian Markov-chain-Monte-Carlo inversion of time-lapse cross hole ground-penetrating radar data to characterize the vadose zone at the Arrenaes field site, Denmark

The ground-penetrating radar (GPR) geophysical method has the potential to provide valuable information on the hydraulic properties of the vadose zone because of its strong sensitivity to soil water content. In particular, recent evidence has suggested that the stochastic inversion of crosshole GPR traveltime data can allow for a significant reduction in uncertainty regarding subsurface van Genuchten-Mualem (VGM) parameters. Much of the previous work on the stochastic estimation of VGM parameters from crosshole GPR data has considered the case of steady-state infiltration conditions, which represent only a small fraction of practically relevant scenarios. We explored in detail the dynamic infiltration case, specifically examining to what extent time-lapse crosshole GPR traveltimes, measured during a forced infiltration experiment at the Arreneas field site in Denmark, could help to quantify VGM parameters and their uncertainties in a layered medium, as well as the corresponding soil hydraulic properties. We used a Bayesian Markov-chain-Monte-Carlo inversion approach. We first explored the advantages and limitations of this approach with regard to a realistic synthetic example before applying it to field measurements. In our analysis, we also considered different degrees of prior information. Our findings indicate that the stochastic inversion of the time-lapse GPR data does indeed allow for a substantial refinement in the inferred posterior VGM parameter distributions compared with the corresponding priors, which in turn significantly improves knowledge of soil hydraulic properties. Overall, the results obtained clearly demonstrate the value of the information contained in time-lapse GPR data for characterizing vadose zone dynamics.

Soil diffuse degassing and thermal energy fluxes from the southern Lakki plain, Nisyros (Greece)

Two diffuse soil CO2 flux surveys from the southern Lakki plain show that CO2 is mainly released from the hydrothermal explosion craters. The correspondence between high CO2 fluxes and elevated soil temperatures suggests that a flux of hot hydrothermal fluids ascends towards the surface. Steam mostly condenses near the surface and the heat given off is conductively transferred to the atmosphere through the soil, accompanied by a large CO2 flux. Tt was calculated, that 68 t d(-1) of hydrothermal CO2 are released through the total surveyed area of similar to1.3 km(2) Admitting that a steam flux of 2200 t d(-1) accompanies this CO2 flux, the thermal energy released through steam condensation amounts to 58 MW.

Controls on 13C and 14C variability in soil CO2

14C dating of groundwater depends on the isotopic composition of both the solid carbonate and the soil CO2 and requires the use of 14C age correction models. To better assess the variability of the 14C activity of soil CO2 (A14Csoil-CO2) and the δ13C of soil CO2 (δ13Csoil-CO2), which are two parameters used in 14C age correction models, we studied the different processes involving carbon isotopes in the soil. The approach used experimental data from two sites in France (Fontainebleau sands and Astian sands) and a steady-state transport model. In most cases, the 14C activity (A14C) of atmospheric CO2 is directly used in the 14C age correction models as the A14Csoil-CO2. However, we demonstrate that since 1950, the evolution of the A14Csoil-CO2 reflects the competition between the fluxes of root-derived CO2 and organic matter-derived CO2. Therefore, the A14Csoil-CO2 must be used to date groundwater that is younger than 60 years old. Moreover, the δ13C of soil CO2 (δ13Csoil-CO2) showed large seasonal variations that must be taken into account in selecting the δ13Csoil-CO2 for 14C age correction models.

Testing the influence of topography and material properties on catchment-scale soil moisture patterns using remotely sensed vegetation patterns in a humid temperate catchment, northern Britain

In order to evaluate the relationship between the apparent complexity of hillslope soil moisture and the emergent patterns of catchment hydrological behaviour and water quality, we need fine-resolution catchment-wide data on soil moisture characteristics. This study proposes a methodology whereby vegetation patterns obtained from high-resolution orthorectified aerial photographs are used as an indicator of soil moisture characteristics. This enables us to examine a set of hypotheses regarding what drives the spatial patterns of soil moisture at the catchment scale (material properties or topography). We find that the pattern of Juncus effusus vegetation is controlled largely by topography and mediated by the catchment's material properties. Characterizing topography using the topographic index adds value to the soil moisture predictions relative to slope or upslope contributing area (UCA). However, these predictions depart from the observed soil moisture patterns at very steep slopes or low UCAs. Copyright (c) 2012 John Wiley & Sons, Ltd.