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

Performance Evaluation of Concrete Pavement Granular Subbase - Pavement Surface Condition Evaluation, TR-554, 2008

This research project covered a wide range of activities that allowed researchers to understand the relationship between stability, pavement distress, and recycled portland cement concrete (RPCC) subbase aggregate materials. Detailed laboratory and field tests, including pavement distress surveys, were conducted at 26 sites in Iowa. Findings show that specific gravities of RPCC are lower than those of crushed limestone. RPCC aggregate material varies from poorly or well-graded sand to gravel. A modified Micro-Deval test procedure showed that abrasion losses of virgin aggregate materials were within the maximum Micro-Deval abrasion loss of 30% recommended by ASTM D6028-06. Micro-Deval abrasion loss of RPCC aggregate materials, however, was much higher than that of virgin materials and exceeded 30% loss. Modulus of elasticity of RPCC subbase materials is high but variable. RPCC subbase layers normally have low permeability. The pavement surfaces for both virgin and RPCC subbase across Iowa were evaluated to fulfill the objectives of this study related to field evaluation. Visual distress surveys were conducted to gather the detailed current pavement condition information including the type, extent, and severity of the pavement distresses. The historical pavement condition information for the surveyed field sections was extracted from the Iowa DOT's Pavement Management Information System (PMIS). The current surface condition of existing field pavements with RPCC subbase was compared with the virgin aggregate subbase sections using two different approaches. The changes in pavement condition indices (PCI and IRI) with time for both types of pavements (subbases) were compared.

Hydro-Surface Preparation and Coating for Painted Structural Steel, TR-407, 1999

Cities and counties in Iowa have more than 8,890 steel bridges, most of which are painted with red lead paint. The Iowa Department of Transportation (Iowa DOT) maintains less than 35 bridges coated with red lead paint, including seven of the large border bridges over the Mississippi and Missouri Rivers. Because of the federal and state regulations for bridge painting, many governmental agencies have opted not to repaint, or otherwise maintain, lead paint coatings. Consequently, the paint condition on many of these bridges is poor, and some bridges are experiencing severe rusting of structural members. This research project was developed with two objectives: 1) to evaluate the effectiveness of preparing the structural steel surface of a bridge with high pressure water jetting instead of abrasive blasting and 2) to coat the structural steel surface with a moisture-cured polyurethane paint under different surface preparation conditions.

Identification of Laboratory Techniques to Optimize Superpave HMA Surface Friction Characteristics, TR-450, 2010

Wet pavement friction is known to be one of the most important roadway safety parameters. In this research, frictional properties of flexible (asphalt) pavements were investigated. As a part of this study, a laboratory device to polish asphalt specimens was refined and a procedure to evaluate mixture frictional properties was proposed. Following this procedure, 46 different Superpave mixtures, one stone matrix asphalt (SMA) mixture and one porous friction course (PFC) mixture were tested. In addition, 23 different asphalt and two concrete field sections were also tested for friction and noise. The results of both field and laboratory measurements were used to develop an International Friction Index (IFI)-based protocol for measurement of the frictional characteristics of asphalt pavements for laboratory friction measurements. Based on the results of the study, it appears the content of high friction aggregate should be 20% or more of the total aggregate blend when used with other, polish susceptible coarse aggregates; the frictional properties increased substantially as the friction aggregate content increased above 20%. Both steel slag and quartzite were found to improve the frictional properties of the blend, though steel slag had a lower polishing rate. In general, mixes containing soft limestone demonstrated lower friction values than comparable mixes with hard limestone or dolomite. Larger nominal maximum aggregate size mixes had better overall frictional performance than smaller sized mixes. In addition, mixes with higher fineness moduli generally had higher macrotexture and friction.

Examination of Curing Criteria for Cold In-Place Recycling, TR-553, 2008

Cold In-Place Recycling (CIR) has been used widely in rehabilitating the rural highways because it improves a long-term pavement performance. A CIR layer is normally covered by a hot mix asphalt (HMA) overlay in order to protect it from water ingress and traffic abrasion and obtain the required pavement structure and texture. Curing is the term currently used for the period of time that a CIR layer should remain exposed to drying conditions before an HMA overlay is placed. The industry standard for curing time is 10 days to 14 days or a maximum moisture content of 1.5 percent, which appear to be very conservative. When the exposed CIR layer is required to carry traffic for many weeks before the wearing surface is placed, it increases the risk of a premature failure in both CIR layer and overlay. This study was performed to explore technically sound ways to identify minimum in-place CIR properties necessary to permit placement of the HMA overlay. To represent the curing process of CIR pavement in the field construction, three different laboratory curing procedures were examined: 1) uncovered, 2) semi-covered and 3) covered specimens. The indirect tensile strength of specimens in all three curing conditions did not increase during an early stage of curing but increased during a later stage of curing usually when the moisture content falls below 1.5%. Dynamic modulus and flow number increased as curing time increased and moisture contents decreased. For the same curing time, CIR-foam specimens exhibited the higher tensile strength and less moisture content than CIR-emulsion. The laboratory test results concluded that the method of curing temperature and length of the curing period significantly affect the properties of the CIR mixtures. The moisture loss index was developed to predict the moisture condition in the field and, in the future, this index be calibrated with the measurements of temperature and moisture of a CIR layer in the field.

Covalent cell surface functionalization of human fetal osteoblasts for tissue engineering.

The chemical functionalization of cell-surface proteins of human primary fetal bone cells with hydrophilic bioorthogonal intermediates was investigated. Toward this goal, chemical pathways were developed for click reaction-mediated coupling of alkyne derivatives with cellular azido-expressing proteins. The incorporation via a tetraethylene glycol linker of a dipeptide and a reporter biotin allowed the proof of concept for the introduction of cell-specific peptide ligands and to follow the reaction in living cells. Tuning the conditions of the click reaction resulted in chemical functionalization of living human fetal osteoblasts with excellent cell survival.

Low-Cost Rural Surface Alternatives: Demonstration Project, TR-664, 2015

The goals of this project were to implement several stabilization methods for preventing or mitigating freeze-thaw damage to granular surfaced roads and identify the most effective and economical methods for the soil and climate conditions of Iowa. Several methods and technologies identified as potentially suitable for Iowa were selected from an extensive analysis of existing literature provided with Iowa Highway Research Board (IHRB) Project TR-632. Using the selected methods, demonstration sections were constructed in Hamilton County on a heavily traveled two-mile section of granular surfaced road that required frequent maintenance during previous thawing periods. Construction procedures and costs of the demonstration sections were documented, and subsequent maintenance requirements were tabulated through two seasonal freeze-thaw periods. Extensive laboratory and field tests were performed prior to construction, as well as before and after the two seasonal freeze-thaw periods, to monitor the performance of the demonstration sections. A weather station was installed at the project site and temperature sensors were embedded in the subgrade to monitor ground temperatures up to a depth of 5 ft and determine the duration and depths of ground freezing and thawing. An economic analysis was performed using the documented construction and maintenance costs, and the estimated cumulative costs per square yard were projected over a 20-year timeframe to determine break-even periods relative to the cost of continuing current maintenance practices. Overall, the sections with biaxial geogrid or macadam base courses had the best observed freeze-thaw performance in this study. These two stabilization methods have larger initial costs and longer break-even periods than aggregate columns, but counties should also weigh the benefits of improved ride quality and savings that these solutions can provide as excellent foundations for future paving or surface upgrades.

Assessment of diesel exhaust particulate exposure and surface characteristics in association with levels of oxidative stress biomarkers

Exposure to PM10 and PM2.5 (particulate matter with aerodynamic diameter smaller than 10 μm and 2.5 μm, respectively) is associated with a range of adverse health effects, including cancer, pulmonary and cardiovascular diseases. Surface characteristics (chemical reactivity, surface area) are considered of prime importance to understand the mechanisms which lead to harmful effects. A hypothetical mechanism to explain these adverse effects is the ability of components (organics, metal ions) adsorbed on these particles to generate Reactive Oxygen Species (ROS), and thereby to cause oxidative stress in biological systems (Donaldson et al., 2003). ROS can attack almost any cellular structure, like DNA or cellular membrane, leading to the formation of a wide variety of degradation products which can be used as a biomarker of oxidative stress. The aim of the present research project is to test whether there is a correlation between the exposure to Diesel Exhaust Particulate (DEP) and the oxidative stress status. For that purpose, a survey has been conducted in real occupational situations where workers were exposed to DEP (bus depots). Different exposure variables have been considered: - particulate number, size distribution and surface area (SMPS); - particulate mass - PM2.5 and PM4 (gravimetry); - elemental and organic carbon (coulometry); - total adsorbed heavy metals - iron, copper, manganese (atomic adsorption); - surface functional groups present on aerosols (Knudsen flow reactor). Several biomarkers of oxidative stress (8-hydroxy-2'-deoxyguanosine and several aldehydes) have been determined either in urine or serum of volunteers. Results obtained during the sampling campaign in several bus depots indicated that the occupational exposure to particulates in these places was rather low (40-50 μg/m3 for PM4). Bimodal size distributions were generally observed (5 μm and <1 μm). Surface characteristics of PM4 varied strongly, depending on the bus depot. They were usually characterized by high carbonyl and low acidic sites content. Among the different biomarkers which have been analyzed within the framework of this study, mean urinary levels of 8-hydroxy-2'-deoxyguanosine increased significantly (p<0.05) during two consecutive days of exposure for non-smoker workers. On the other hand, no statistically significant differences were observed for serum levels of hexanal, nonanal and 4- hydroxy-nonenal (p>0.05). Biomarkers levels will be compared to exposure variables to gain a better understanding of the relation between the particulate characteristics and the formation of ROS by-products. This project is financed by the Swiss State Secretariat for Education and Research. It is conducted within the framework of the COST Action 633 "Particulate Matter - Properties Related to Health Effects".

Surface modification of biomaterials for conjugation with human fetal osteoblasts.

Surface functionalization of hydroxyapatite (HA) and beta-tricalcium phosphate (TCP) bioceramics with chemical ligands containing a pyrrogallol moiety was developed to improve the adhesion of bone cell precursors to the biomaterials. Fast and biocompatible copper-free click reaction with azido-modified human fetal osteoblasts resulted in improved cell binding to both HA and TCP bioceramics, opening the way for using this methodology in the preparation of cell-engineered bone implants.

Development of a Subgrade Drainage Model for Unpaved Roads, TR-654, 2015

With over 68 thousand miles of gravel roads in Iowa and the importance of these roads within the farm-to-market transportation system, proper water management becomes critical for maintaining the integrity of the roadway materials. However, the build-up of water within the aggregate subbase can lead to frost boils and ultimately potholes forming at the road surface. The aggregate subbase and subgrade soils under these gravel roads are produced with material opportunistically chosen from local sources near the site and, many times, the compositions of these sublayers are far from ideal in terms of proper water drainage with the full effects of this shortcut not being fully understood. The primary objective of this project was to provide a physically-based model for evaluating the drainability of potential subbase and subgrade materials for gravel roads in Iowa. The Richards equation provided the appropriate framework to study the transient unsaturated flow that usually occurs through the subbase and subgrade of a gravel road. From which, we identified that the saturated hydraulic conductivity, Ks, was a key parameter driving the time to drain of subgrade soils found in Iowa, thus being a good proxy variable for accessing roadway drainability. Using Ks, derived from soil texture, we were able to identify potential problem areas in terms of roadway drainage . It was found that there is a threshold for Ks of 15 cm/day that determines if the roadway will drain efficiently, based on the requirement that the time to drain, Td, the surface roadway layer does not exceed a 2-hr limit. Two of the three highest abundant textures (loam and silty clay loam), which cover nearly 60% of the state of Iowa, were found to have average Td values greater than the 2-hr limit. With such a large percentage of the state at risk for the formation of boils due to the soil with relatively low saturated hydraulic conductivity values, it seems pertinent that we propose alternative design and/or maintenance practices to limit the expensive repair work in Iowa. The addition of drain tiles or French mattresses my help address drainage problems. However, before pursuing this recommendation, a comprehensive cost-benefit analysis is needed.

Grade and Cross Slope Estimation from LIDAR-based Surface Models, 2003

Many transportation agencies maintain grade as an attribute in roadway inventory databases; however, the information is often in an aggregated format. Cross slope is rarely included in large roadway inventories. Accurate methods available to collect grade and cross slope include global positioning systems, traditional surveying, and mobile mapping systems. However, most agencies do not have the resources to utilize these methods to collect grade and cross slope on a large scale. This report discusses the use of LIDAR to extract roadway grade and cross slope for large-scale inventories. Current data collection methods and their advantages and disadvantages are discussed. A pilot study to extract grade and cross slope from a LIDAR data set, including methodology, results, and conclusions, is presented. This report describes the regression methodology used to extract and evaluate the accuracy of grade and cross slope from three dimensional surfaces created from LIDAR data. The use of LIDAR data to extract grade and cross slope on tangent highway segments was evaluated and compared against grade and cross slope collected using an automatic level for 10 test segments along Iowa Highway 1. Grade and cross slope were measured from a surface model created from LIDAR data points collected for the study area. While grade could be estimated to within 1%, study results indicate that cross slope cannot practically be estimated using a LIDAR derived surface model.

Conservation of the Eurasian Lynx (Lynx lynx) in a fragmented landscape - habitat models, dispersal and potential distribution

Abstract: The expansion of a recovering population - whether re-introduced or spontaneously returning - is shaped by (i) biological (intrinsic) factors such as the land tenure system or dispersal, (ii) the distribution and availability of resources (e.g. prey), (iii) habitat and landscape features, and (iv) human attitudes and activities. In order to develop efficient conservation and recovery strategies, we need to understand all these factors and to predict the potential distribution and explore ways to reach it. An increased number of lynx in the north-western Swiss Alps in the nineties lead to a new controversy about the return of this cat. When the large carnivores were given legal protection in many European countries, most organizations and individuals promoting their protection did not foresee the consequences. Management plans describing how to handle conflicts with large predators are needed to find a balance between "overabundance" and extinction. Wildlife and conservation biologists need to evaluate the various threats confronting populations so that adequate management decisions can be taken. I developed a GIS probability model for the lynx, based on habitat information and radio-telemetry data from the Swiss Jura Mountains, in order to predict the potential distribution of the lynx in this mountain range, which is presently only partly occupied by lynx. Three of the 18 variables tested for each square kilometre describing land use, vegetation, and topography, qualified to predict the probability of lynx presence. The resulting map was evaluated with data from dispersing subadult lynx. Young lynx that were not able to establish home ranges in what was identified as good lynx habitat did not survive their first year of independence, whereas the only one that died in good lynx habitat was illegally killed. Radio-telemetry fixes are often used as input data to calibrate habitat models. Radio-telemetry is the only way to gather accurate and unbiased data on habitat use of elusive larger terrestrial mammals. However, it is time consuming and expensive, and can therefore only be applied in limited areas. Habitat models extrapolated over large areas can in turn be problematic, as habitat characteristics and availability may change from one area to the other. I analysed the predictive power of Ecological Niche Factor Analysis (ENFA) in Switzerland with the lynx as focal species. According to my results, the optimal sampling strategy to predict species distribution in an Alpine area lacking available data would be to pool presence cells from contrasted regions (Jura Mountains, Alps), whereas in regions with a low ecological variance (Jura Mountains), only local presence cells should be used for the calibration of the model. Dispersal influences the dynamics and persistence of populations, the distribution and abundance of species, and gives the communities and ecosystems their characteristic texture in space and time. Between 1988 and 2001, the spatio-temporal behaviour of subadult Eurasian lynx in two re-introduced populations in Switzerland was studied, based on 39 juvenile lynx of which 24 were radio-tagged to understand the factors influencing dispersal. Subadults become independent from their mothers at the age of 8-11 months. No sex bias neither in the dispersal rate nor in the distance moved was detected. Lynx are conservative dispersers, compared to bear and wolf, and settled within or close to known lynx occurrences. Dispersal distances reached in the high lynx density population - shorter than those reported in other Eurasian lynx studies - are limited by habitat restriction hindering connections with neighbouring metapopulations. I postulated that high lynx density would lead to an expansion of the population and validated my predictions with data from the north-western Swiss Alps where about 1995 a strong increase in lynx abundance took place. The general hypothesis that high population density will foster the expansion of the population was not confirmed. This has consequences for the re-introduction and recovery of carnivores in a fragmented landscape. To establish a strong source population in one place might not be an optimal strategy. Rather, population nuclei should be founded in several neighbouring patches. Exchange between established neighbouring subpopulations will later on take place, as adult lynx show a higher propensity to cross barriers than subadults. To estimate the potential population size of the lynx in the Jura Mountains and to assess possible corridors between this population and adjacent areas, I adapted a habitat probability model for lynx distribution in the Jura Mountains with new environmental data and extrapolated it over the entire mountain range. The model predicts a breeding population ranging from 74-101 individuals and from 51-79 individuals when continuous habitat patches < 50 km2 are disregarded. The Jura Mountains could once be part of a metapopulation, as potential corridors exist to the adjoining areas (Alps, Vosges Mountains, and Black Forest). Monitoring of the population size, spatial expansion, and the genetic surveillance in the Jura Mountains must be continued, as the status of the population is still critical. ENFA was used to predict the potential distribution of lynx in the Alps. The resulting model divided the Alps into 37 suitable habitat patches ranging from 50 to 18,711 km2, covering a total area of about 93,600 km2. When using the range of lynx densities found in field studies in Switzerland, the Alps could host a population of 961 to 1,827 residents. The results of the cost-distance analysis revealed that all patches were within the reach of dispersing lynx, as the connection costs were in the range of dispersal cost of radio-tagged subadult lynx moving through unfavorable habitat. Thus, the whole Alps could once be considered as a metapopulation. But experience suggests that only few disperser will cross unsuitable areas and barriers. This low migration rate may seldom allow the spontaneous foundation of new populations in unsettled areas. As an alternative to natural dispersal, artificial transfer of individuals across the barriers should be considered. Wildlife biologists can play a crucial role in developing adaptive management experiments to help managers learning by trial. The case of the lynx in Switzerland is a good example of a fruitful cooperation between wildlife biologists, managers, decision makers and politician in an adaptive management process. This cooperation resulted in a Lynx Management Plan which was implemented in 2000 and updated in 2004 to give the cantons directives on how to handle lynx-related problems. This plan was put into practice e.g. in regard to translocation of lynx into unsettled areas. Résumé: L'expansion d'une population en phase de recolonisation, qu'elle soit issue de réintroductions ou d'un retour naturel dépend 1) de facteurs biologiques tels que le système social et le mode de dispersion, 2) de la distribution et la disponibilité des ressources (proies), 3) de l'habitat et des éléments du paysage, 4) de l'acceptation de l'espèce par la population locale et des activités humaines. Afin de pouvoir développer des stratégies efficaces de conservation et de favoriser la recolonisation, chacun de ces facteurs doit être pris en compte. En plus, la distribution potentielle de l'espèce doit pouvoir être déterminée et enfin, toutes les possibilités pour atteindre les objectifs, examinées. La phase de haute densité que la population de lynx a connue dans les années nonante dans le nord-ouest des Alpes suisses a donné lieu à une controverse assez vive. La protection du lynx dans de nombreux pays européens, promue par différentes organisations, a entraîné des conséquences inattendues; ces dernières montrent que tout plan de gestion doit impérativement indiquer des pistes quant à la manière de gérer les conflits, tout en trouvant un équilibre entre l'extinction et la surabondance de l'espèce. Les biologistes de la conservation et de la faune sauvage doivent pour cela évaluer les différents risques encourus par les populations de lynx, afin de pouvoir rapidement prendre les meilleuresmdécisions de gestion. Un modèle d'habitat pour le lynx, basé sur des caractéristiques de l'habitat et des données radio télémétriques collectées dans la chaîne du Jura, a été élaboré afin de prédire la distribution potentielle dans cette région, qui n'est que partiellement occupée par l'espèce. Trois des 18 variables testées, décrivant pour chaque kilomètre carré l'utilisation du sol, la végétation ainsi que la topographie, ont été retenues pour déterminer la probabilité de présence du lynx. La carte qui en résulte a été comparée aux données télémétriques de lynx subadultes en phase de dispersion. Les jeunes qui n'ont pas pu établir leur domaine vital dans l'habitat favorable prédit par le modèle n'ont pas survécu leur première année d'indépendance alors que le seul individu qui est mort dans l'habitat favorable a été braconné. Les données radio-télémétriques sont souvent utilisées pour l'étalonnage de modèles d'habitat. C'est un des seuls moyens à disposition qui permette de récolter des données non biaisées et précises sur l'occupation de l'habitat par des mammifères terrestres aux moeurs discrètes. Mais ces méthodes de- mandent un important investissement en moyens financiers et en temps et peuvent, de ce fait, n'être appliquées qu'à des zones limitées. Les modèles d'habitat sont ainsi souvent extrapolés à de grandes surfaces malgré le risque d'imprécision, qui résulte des variations des caractéristiques et de la disponibilité de l'habitat d'une zone à l'autre. Le pouvoir de prédiction de l'Analyse Ecologique de la Niche (AEN) dans les zones où les données de présence n'ont pas été prises en compte dans le calibrage du modèle a été analysée dans le cas du lynx en Suisse. D'après les résultats obtenus, la meilleure mé- thode pour prédire la distribution du lynx dans une zone alpine dépourvue d'indices de présence est de combiner des données provenant de régions contrastées (Alpes, Jura). Par contre, seules les données sur la présence locale de l'espèce doivent être utilisées pour les zones présentant une faible variance écologique tel que le Jura. La dispersion influence la dynamique et la stabilité des populations, la distribution et l'abondance des espèces et détermine les caractéristiques spatiales et temporelles des communautés vivantes et des écosystèmes. Entre 1988 et 2001, le comportement spatio-temporel de lynx eurasiens subadultes de deux populations réintroduites en Suisse a été étudié, basé sur le suivi de 39 individus juvéniles dont 24 étaient munis d'un collier émetteur, afin de déterminer les facteurs qui influencent la dispersion. Les subadultes se sont séparés de leur mère à l'âge de 8 à 11 mois. Le sexe n'a pas eu d'influence sur le nombre d'individus ayant dispersés et la distance parcourue au cours de la dispersion. Comparé à l'ours et au loup, le lynx reste très modéré dans ses mouvements de dispersion. Tous les individus ayant dispersés se sont établis à proximité ou dans des zones déjà occupées par des lynx. Les distances parcourues lors de la dispersion ont été plus courtes pour la population en phase de haute densité que celles relevées par les autres études de dispersion du lynx eurasien. Les zones d'habitat peu favorables et les barrières qui interrompent la connectivité entre les populations sont les principales entraves aux déplacements, lors de la dispersion. Dans un premier temps, nous avons fait l'hypothèse que les phases de haute densité favorisaient l'expansion des populations. Mais cette hypothèse a été infirmée par les résultats issus du suivi des lynx réalisé dans le nord-ouest des Alpes, où la population connaissait une phase de haute densité depuis 1995. Ce constat est important pour la conservation d'une population de carnivores dans un habitat fragmenté. Ainsi, instaurer une forte population source à un seul endroit n'est pas forcément la stratégie la plus judicieuse. Il est préférable d'établir des noyaux de populations dans des régions voisines où l'habitat est favorable. Des échanges entre des populations avoisinantes pourront avoir lieu par la suite car les lynx adultes sont plus enclins à franchir les barrières qui entravent leurs déplacements que les individus subadultes. Afin d'estimer la taille de la population de lynx dans le Jura et de déterminer les corridors potentiels entre cette région et les zones avoisinantes, un modèle d'habitat a été utilisé, basé sur un nouveau jeu de variables environnementales et extrapolé à l'ensemble du Jura. Le modèle prédit une population reproductrice de 74 à 101 individus et de 51 à 79 individus lorsque les surfaces d'habitat d'un seul tenant de moins de 50 km2 sont soustraites. Comme des corridors potentiels existent effectivement entre le Jura et les régions avoisinantes (Alpes, Vosges, et Forêt Noire), le Jura pourrait faire partie à l'avenir d'une métapopulation, lorsque les zones avoisinantes seront colonisées par l'espèce. La surveillance de la taille de la population, de son expansion spatiale et de sa structure génétique doit être maintenue car le statut de cette population est encore critique. L'AEN a également été utilisée pour prédire l'habitat favorable du lynx dans les Alpes. Le modèle qui en résulte divise les Alpes en 37 sous-unités d'habitat favorable dont la surface varie de 50 à 18'711 km2, pour une superficie totale de 93'600 km2. En utilisant le spectre des densités observées dans les études radio-télémétriques effectuées en Suisse, les Alpes pourraient accueillir une population de lynx résidents variant de 961 à 1'827 individus. Les résultats des analyses de connectivité montrent que les sous-unités d'habitat favorable se situent à des distances telles que le coût de la dispersion pour l'espèce est admissible. L'ensemble des Alpes pourrait donc un jour former une métapopulation. Mais l'expérience montre que très peu d'individus traverseront des habitats peu favorables et des barrières au cours de leur dispersion. Ce faible taux de migration rendra difficile toute nouvelle implantation de populations dans des zones inoccupées. Une solution alternative existe cependant : transférer artificiellement des individus d'une zone à l'autre. Les biologistes spécialistes de la faune sauvage peuvent jouer un rôle important et complémentaire pour les gestionnaires de la faune, en les aidant à mener des expériences de gestion par essai. Le cas du lynx en Suisse est un bel exemple d'une collaboration fructueuse entre biologistes de la faune sauvage, gestionnaires, organes décisionnaires et politiciens. Cette coopération a permis l'élaboration du Concept Lynx Suisse qui est entré en vigueur en 2000 et remis à jour en 2004. Ce plan donne des directives aux cantons pour appréhender la problématique du lynx. Il y a déjà eu des applications concrètes sur le terrain, notamment par des translocations d'individus dans des zones encore inoccupées.

Simulation of surface waves in porous media

We present a novel numerical algorithm for the simulation of seismic wave propagation in porous media, which is particularly suitable for the accurate modelling of surface wave-type phenomena. The differential equations of motion are based on Biot's theory of poro-elasticity and solved with a pseudospectral approach using Fourier and Chebyshev methods to compute the spatial derivatives along the horizontal and vertical directions, respectively. The time solver is a splitting algorithm that accounts for the stiffness of the differential equations. Due to the Chebyshev operator the grid spacing in the vertical direction is non-uniform and characterized by a denser spatial sampling in the vicinity of interfaces, which allows for a numerically stable and accurate evaluation of higher order surface wave modes. We stretch the grid in the vertical direction to increase the minimum grid spacing and reduce the computational cost. The free-surface boundary conditions are implemented with a characteristics approach, where the characteristic variables are evaluated at zero viscosity. The same procedure is used to model seismic wave propagation at the interface between a fluid and porous medium. In this case, each medium is represented by a different grid and the two grids are combined through a domain-decomposition method. This wavefield decomposition method accounts for the discontinuity of variables and is crucial for an accurate interface treatment. We simulate seismic wave propagation with open-pore and sealed-pore boundary conditions and verify the validity and accuracy of the algorithm by comparing the numerical simulations to analytical solutions based on zero viscosity obtained with the Cagniard-de Hoop method. Finally, we illustrate the suitability of our algorithm for more complex models of porous media involving viscous pore fluids and strongly heterogeneous distributions of the elastic and hydraulic material properties.

Evaluación de los factores de riesgo y los tipos de superficie para el desarrollo de las úlceras por presión en el enfermo crítico = Assessment of risk factors and the types of surface for the development of pressure ulcers on critical ill patients

The main objective of this article is to assess the risk factors and the types of surface for the development of pressure ulcers (PU) on critical ill patients in an Intensive Care Unit (ICU)

Analysis of sensitivity with respect to a compositional parameter : A comment on 'Local sensitivity analysis for compositional data with application to soil texture in hydrologic modelling' by L. Loosvelt and co-authors

A comment about the article “Local sensitivity analysis for compositional data with application to soil texture in hydrologic modelling” writen by L. Loosvelt and co-authors. The present comment is centered in three specific points. The first one is related to the fact that the authors avoid the use of ilr-coordinates. The second one refers to some generalization of sensitivity analysis when input parameters are compositional. The third tries to show that the role of the Dirichlet distribution in the sensitivity analysis is irrelevant