Estudi de la contaminació per nitrats a les aigües subterrànies de Catalunya: l'Alt Penedès, el Baix Camp, el Baix Empordà, el Baix Penedès, la Garrotxa, el Maresme, Osona, el Pla de l'Estany, la Segarra, el Segrià, el Tarragonès i l'Urgell

La contaminació per nitrats de les aigües subterrànies és un dels problemes mediambientals i de salut pública que més afecten els nostres aqüífers. És per això que un dels eixos del nostre estudi és fer la màxima difusió de l’estat de la qualitat de l’aigua dels aqüífers de diferents àrees de Catalunya. En el present estudi es tracta la contaminació per nitrats a les aigües de 945 punts de mostreig (sobretot d’aigua procedent de pous, d’algunes mines i d’algunes fonts) de Catalunya, situats en 12 comarques diferents, totes elles amb algun o una gran part dels seus municipis englobats en una zona vulnerable. D’aquests punts de mostreig s’ha fet un recull de dades de la concentració de nitrats des de l’any 2003 fins el 2012 per tal d’estudiar-ne l’evolució i la incidència de la pluviometria, així com conèixer en quines comarques les variables externes: superfície agrària útil (SAU), explotacions agràries i cabana porcina i bovina poden haver incidit directament en la contaminació. També s’estudia fins a quin punt, aquest deteriorament, ha afectat la mineralització de l’aigua o afecta a llocs d’aigües més mineralitzades. A partir dels resultats obtinguts s’ha arribat a les conclusions que la concentració mitjana de nitrats a les aigües subterrànies en la majoria de comarques estudiades és elevada, superior als 50 mg/l que estableix la normativa vigent; que ha tendit a l’estabilitat en els últims deu anys; que afecta llocs d’aigües molt mineralitzades i que la pluviometria incideix notablement en la concentració de nitrats. S’ha observat que les variables externes superfície agrària útil (SAU), explotacions agràries i cabana porcina i bovina incideixen directament en la contaminació, excepte al Maresme. Finalment es presenten diferents usos alternatius als de l’aigua de consum per aquestes aigües contaminades amb nitrats i que han posat en pràctica diferents municipis de les comarques estudiades.

Geohydrology of Muscatine Island, Muscatine County, Iowa, 1977

The availability of large quantities of high-quality ground water from the Muscatine Island aquifer has had a tremendous impact upon urban, industrial and agricultural development of this part of Iowa. Although the nonpumping level of water has been lowered significantly near major pumping centers through time, proper management of this water resource can assure a continued supply of water for all competing users. This report provides basic information for long range management. Fold out maps are included.

Estimating traveltimes and groundwater flow patterns using 3D time-lapse crosshole ERT imaging of electrical resistivity fluctuations induced by infiltrating river water

The infiltration of river water into aquifers is of high relevance to drinking-water production and is a key driver of biogeochemical processes in the hyporheic and riparian zone, but the distribution and quantification of the infiltrating water are difficult to determine using conventional hydrological methods (e.g., borehole logging and tracer tests). By time-lapse inverting crosshole ERT (electrical resistivity tomography) monitoring data, we imaged groundwater flow patterns driven by river water infiltrating a perialpine gravel aquifer in northeastern Switzerland. This was possible because the electrical resistivity of the infiltrating water changed during rainfall-runoff events. Our time-lapse resistivity models indicated rather complex flow patterns as a result of spatially heterogeneous bank filtration and aquifer heterogeneity. The upper part of the aquifer was most affected by the river infiltrate, and the highest groundwater velocities and possible preferential flow occurred at shallow to intermediate depths. Time series of the reconstructed resistivity models matched groundwater electrical resistivity data recorded on borehole loggers in the upper and middle parts of the aquifer, whereas the resistivity models displayed smaller variations and delayed responses with respect to the logging data. in the lower part. This study demonstrated that crosshole ERT monitoring of natural electrical resistivity variations of river infiltrate could be used to image and quantify 3D bank filtration and aquifer dynamics at a high spatial resolution.

Stochastic inversion of tracer test and electrical geophysical data to estimate hydraulic conductivities.

Quantifying the spatial configuration of hydraulic conductivity (K) in heterogeneous geological environments is essential for accurate predictions of contaminant transport, but is difficult because of the inherent limitations in resolution and coverage associated with traditional hydrological measurements. To address this issue, we consider crosshole and surface-based electrical resistivity geophysical measurements, collected in time during a saline tracer experiment. We use a Bayesian Markov-chain-Monte-Carlo (McMC) methodology to jointly invert the dynamic resistivity data, together with borehole tracer concentration data, to generate multiple posterior realizations of K that are consistent with all available information. We do this within a coupled inversion framework, whereby the geophysical and hydrological forward models are linked through an uncertain relationship between electrical resistivity and concentration. To minimize computational expense, a facies-based subsurface parameterization is developed. The Bayesian-McMC methodology allows us to explore the potential benefits of including the geophysical data into the inverse problem by examining their effect on our ability to identify fast flowpaths in the subsurface, and their impact on hydrological prediction uncertainty. Using a complex, geostatistically generated, two-dimensional numerical example representative of a fluvial environment, we demonstrate that flow model calibration is improved and prediction error is decreased when the electrical resistivity data are included. The worth of the geophysical data is found to be greatest for long spatial correlation lengths of subsurface heterogeneity with respect to wellbore separation, where flow and transport are largely controlled by highly connected flowpaths.

Temporal changes of an alfalfa succession and related soil physical properties on the Loess Plateau, China

The objective of this work was to investigate the relationship between changes in the plant community and changes in soil physical properties and water availability, during a succession from alfalfa (Medicago sativa L.) to natural vegetation on the Loess Plateau, China. Data from a succession sere spanning 32 years were collated, and vegetative indexes were compared to changes related to soil bulk density and soil water storage. The alfalfa yield increased for approximately 7 years, then it declined and the alfalfa was replaced by a natural community dominated by Stipa bungeana that began to thrive about 10 years after alfalfa seeding. Soil bulk density increased over time, but the deterioration of the alfalfa was mainly ascribed to a severe reduction in soil water storage, which was lowest around the time when degradation commenced. The results indicated that water consumption by alfalfa could be reduced by reducing plant density. The analysis of the data also suggested that soil water recharge could be facilitated by rotating the alfalfa with other crops, natural vegetation, or bare soil.

Iowa Source Water Protection: Workbook, November, 18, 2014

There are a few basic fundamentals you need before starting a source water protection project. These include information on your community’s wells (or intakes), aquifer, source water area, and potential contaminants. All of these essential items should be included in your community’s source water information, you may find this information in the workbook and guidebook.

Iowa Source Water Protection: Guidebook, November, 18, 2014

There are a few basic fundamentals you need before starting a source water protection project. These include information on your community’s wells (or intakes), aquifer, source water area, and potential contaminants. All of these essential items should be included in your community’s source water information, you may find this information in the workbook and guidebook.

Chemical and isotopic equilibrium between CO2 and CH4 in fumarolic gas discharges: Generation of CH4 in arc magmatic-hydrothermal systems

The chemical and isotopic composition of fumarolic gases emitted from Nisyros Volcano, Greece, and of a single gas sample from Vesuvio, Italy, was investigated in order to determine the origin of methane (CH,) within two subduction-related magmatic-hydrothermal environments. Apparent temperatures derived from carbon isotope partitioning between CH4 and CO2 of around 340degreesC for Nisyros and 470degreesC for Vesuvio correlate well with aquifer temperatures as measured directly and/or inferred from compositional data using the H2O-H-2-CO2-CO-CH4 geothermometer. Thermodynamic modeling reveals chemical equilibrium between CH4, CO2 and H2O implying that carbon isotope partitioning between CO2 and CH, in both systems is controlled by aquifer temperature. N-2/(3) He and CH4/(3) He ratios of Nisyros fumarolic gases are unusually low for subduction zone gases and correspond to those of midoceanic ridge environments. Accordingly, CH4 may have been primarily generated through the reduction of CO, by H, in the absence of any organic matter following a Fischer-Tropsch-type reaction. However, primary occurrence of minor amounts of thermogenic CH4 and subsequent re-equilibration with co-existing CO2 cannot be ruled out entirely- CO2/He-3 ratios and delta(13)C(CO2) values imply that the evolved CO2 either derives from a metasomatized mantle or is a mixture between two components, one outgassing from an unaltered mantle and the other released by thermal breakdown of marine carbonates. The latter may contain traces of organic matter possibly decomposing to CH4 during thermometamorphism. Copyright (C) 2004 Elsevier Ltd.

Growth mechanisms and timing of emplacement of a vertically layered mafic intrusion in the feeder zone of an ocean island volcano (Fuerteventura, Canary islands)

Résumé pour le grand public L'île de Fuerteventura (Canaries) offre l'occasion rare d'observer les racines d'un volcan océanique édifié il y a 25 à 30 millions d'années et complètement érodé. On y voit de nombreux petits plutons de forme et composition variées, témoignant d'autant d'épisodes de l'activité magmatique. L'un de ces plutons, appelé PX1, présente une structure inhabituelle formée d'une alternance de bandes verticales d'épaisseur métrique à hectométrique de roches sombres de composition pyroxénilique ou gabbroïque. Les pyroxénites résultent clairement de l'accumulation de cristaux de pyroxènes et non de la simple solidification d'un magma? Se pose dès lors la question de la nature du processus qui a conduit à l'accumulation verticale de niveaux concentrés en pyroxènes. En effet, les litages pyroxénitiques classiques sont subhorizontaux, car ils résultent de l'accumulation gravitaire des cristaux séparés du magma dont ils cristalli¬sent par sédimentation. Cette étude vise à identifier et comprendre les mécanismes qui ont engendré ce Iitage minéralogique vertical et l'im¬portant volume de ces faciès cumulatifs. Nous nous sommes également intéressés aux conditions de pression et de température régnant au moment de la mise en place du pluton, ainsi qu'à sa durée de vie et à sa vitesse de refroidis¬sement. Enfin une approche géochimique nous a permis de préciser la nature de la source mantellique des magmas liés à cette activité magmatique. PX1 est en réalité un complexe filonien formé à des conditions de pression et de température de 1-2 kbar et 1050- 1100°C; sa construction a nécessité au moins 150 km3 de magma. L'alternance d'horizons gabbroïques et pyroxéniti¬ques représente des injections successives de magma sous la forme de filons verticaux, mis en place dans un contexte régional en extension. L'étude des orientations des minéraux dans ces faciès révèle que les horizons gabbroïques enregistrent l'extension régionale, alors que les pyroxénites sont générées par une compaction au sein du pluton. Ceci suggère que le régime des contraintes, qui était extensif lors de l'initiation de la mise en place de PX1, est pério¬diquement devenu compressif au sein même du pluton. Cette compression serait liée à des cycles de mise en place où la vitesse de croissance du pluton dépassait celle de l'extension régionale. La différenciation observée au sein de chaque horizon, depuis des pyroxénites riches en olivine jusqu'à des pyroxé¬nites à plagioclase interstitiel et des gabbros, ainsi que la composition géochimique des minéraux qui les constituent suggèrent que chaque filon vertical s'est mis en place à partir d'un magma de composition identique, puis a évolué indépendamment des autres en fonction du régime thermique et du régime des contraintes local. Lorsque le magma en train de cristalliser s'est trouvé en compression, le liquide résiduel a été séparé des cristaux déjà formés et extrait du système, laissant derrière lui une accumulation de cristaux dont la nature et les proportions dépendaient du stade de cristallisation atteint par le magma au moment de l'extraction. Ainsi, les niveaux de pyroxénites à olivine (premier minéral à cristalliser) ont été formés lorsque le magma correspondant était encore peu cristallisé; à l'inverse, les py¬roxénites riches en plagioclase (minéral plus tardif dans la séquence de cristallisation) et certains gabbros à caractère cumulatif résultent d'une compression tardive dans le processus de cristallisation du filon concerné. Les liquides résiduels extraits des niveaux pyroxénitiques sont rarement observés dans PX1, certaines poches et filonets de com¬position anorthositique pourraient en être les témoins. L'essentiel de ces liquides a probablement gagné des niveaux supérieurs du pluton, voire la surface du volcan. L'origine du régime compressif périodique affectant les filons en voie de cristallisation est attribuée aux injections suivantes de magma au sein du pluton, qui se sont succédées à un rythme plus rapide que la vitesse de consolidation des filons. Des datations U/Pb de haute précision sur des cristaux de zircon et de baddeleyite ainsi que40Ar/39Ar sur des cris¬taux d'amphibole révèlent une initiation de la mise en place de PX1 il y a 22.1 ± 0,7 Ma; celle-ci a duré quelque 0,48 ± 0,22 à 0,52 ± 0,29 Ma. Ce laps de temps est compatible avec celui nécessaire à la cristallisation des filons individuels, qui va de moins d'une année lors de l'initiation du magmatisme à 5 ans lors du maximum d'activité de PX1. La présence de cristaux résorbés enregistrant une cristallisation complexe suggère l'existence d'une chambre mag¬matique convective sous-jacente à PX1 et périodiquement rechargée. Les compositions isotopiques des roches étu¬diées révèlent une source mantellique profonde de type point chaud avec une contribution du manteau lithosphéri- que métasomatisé présent sous les îles Canaries. Résumé L'intrusion mafique Miocène PX1 fait partie du soubassement superficiel (0.15-0.2 GPa, 1100 °Q d'un volcan d'île océanique. La particularité de ce pluton est l'existence d'alternances d'unités de gabbros et de pyroxénites qui met¬tent en évidence un litage magmatique vertical (NNE-SSW). Les horizons gabbroiques et pyroxénitiques sont constitués d'unités de différenciation métriques qui suggèrent tine mise en place par injections périodiques de filons verticaux de magma formant un complexe filonien. Chaque filon vertical a subi une différenciation parallèle à un front de solidification sub-vertical parallèle aux bords du filon. Les pyroxénites résultent du fractionnement et de l'accumulation d'olivine ± clinopyroxene ± plagioclase à partir d'un magma basaltique faiblement alcalin et sont interprétées comme étant des imités de différenciation tronquées dont le liquide interstitiel a été extrait par compaction. L'orientation préférentielle des clinopyroxènes dans ces pyroxe- nites (obtenues par analyse EBSD et micro-tomographique) révèle une composante de cisaillement simple dans la genèse de ces roches, ce qui confirme cette interprétation. La compaction des pyroxénites est probablement causée par a mise en place de filons de magma suivants. Le liquide interstitiel expulsé est probablement par ces derniers. Les clinopyroxènes des gabbros, montrent une composante de cisaillement pure suggérant qu'ils sont affectés par une déformation syn-magmatique parallèle aux zones de cisaillement NNE-SSW observées autour de PX1 et liées au contexte tectonique Miocène d'extension régionale. Ceci suggère que les gabbros sont liés à des taux de mise en place faibles à la fin de cycles d'activité magmatique et sont peu ou pas affectés par la compaction. L'initiation et la géométrie de PX1 sont donc contrôlées par le contexte tectonique régional d'extension alors que les taux et les volumes de magma dépendent de facteurs liés à la source. Des taux d'injection élevés résultent probable¬ment en une croissance du pluton supérieure à la place crée par cette extension. Dans ce cas de figure, la propagation des nouveaux dykes et l'inaptitude du magma à circuler à travers les anciens dykes cristallisés pourrait causer une augmentation de la pression non-lithostatique sur ces derniers, exprimée par un cisaillement simple et l'expulsion du liquide interstitiel qu'ils contiennent (documenté par les zones de collecte anorthositiques). Les compositions en éléments majeurs et traces des gabbros et pyroxenites de PX1 sont globalement homogènes et dépendent de la nature cumulative des échantillons. Cependant, de petites variations des concentrations en éléments traces ainsi que les teneurs en éléments traces des bordures de clinopyroxenes suggèrent que ces derniers ont subi un processus de rééquilibrage et de cristallisation in situ. L'homogénéité des compositions chimiques des échantillons, ainsi que la présence de grains de clinopyroxene résorbés suggère que le complexe filonien PX1 s'est mis en place au dessus d'une chambre magmatique périodiquement rechargée dans laquelle la convection est efficace. Chaque filon est donc issu d'un même magma, mais a subi une différenciation par cristallisation in situ (jusqu'à 70% de fraction¬nement) indépendamment des autres. Dans ces filons cristallisés, les minéraux cumulatifs subissent un rééquilibrage partiel avec les liquide interstitiel avant que ce dernier ne soit expulsé lors de la compaction (mettant ainsi un terme à la différenciation). Ce modèle de mise en place signifie qu'un minimum de 150Km3 de magma est nécessaire à la genèse de PX1, une partie de ce volume ayant été émis par le 'Central Volcanic Complex' de Fuerteventura. Les rapports isotopiques radiogéniques mesurés révèlent la contribution de trois pôles mantelliques dans la genèse du magma formant PX1. Le mélange de ces pôles HIMU, DMM et EM1 refléterai l'interaction du point chaud Cana¬rien avec un manteau lithosphérique hétérogène métasomatisé. Les petites variations de ces rapports et des teneurs en éléments traces au sein des faciès pourrait refléter des taux de fusion partielle variable de la source, résultant en un échantillonnage variable du manteau lithosphérique métasomatisé lors de son interaction avec le point chaud. Des datations U/Pb de haute précision (TIMS) sur des cristaux de zircon et de baddeleyite extraits de gabbros de PX1 révèlent que l'initiation de la cristallisation du magma a eu lieu il y a 22.10±0.07 Ma et que l'activité magmatique a duré un minimum de 0.48 à 0.52 Ma. Des âges 40Ar/39Ar obtenus sur amphibole sont de 21.9 ± 0.6 à 21.8 ± 0.3 Ma, identiques aux âges U/Pb. La combinaison de ces méthodes de datations, suggère que le temps maximum nécessaire à PX1 pour se refroidir en dessous de la température de fermeture de l'amphibole est de 0.8Ma. Ceci signifie que la durée de vie de PX1 est de 520 000 à 800 000 ans. La coexistence de cristaux de baddeleyite et de zircon dans un gabbro est attribuée à son interaction avec un fluide riche en C02 relâché par les carbonatites encaissantes lors du métamorphisme de contact généré par la mise en place de PX1 environ 160 000 ans après le début de sa mise en place. Les durées de vie obtenue sont en accord avec le modèle de mise en place suggérant une durée de cristallisation poux chaque filon allant de 1 an à 5 ans. Abstract The Miocene PX1 gabbro-pyroxenite intrusion (Fuerteventura, Canary Islands), is interpreted as the shallow-level feeder-zone (0.15-0.2 GPa and 1100-1120°C), to an ocean island volcano. The particularity of PX1 is that it displays a NNE-SSW trending vertical magmatic banding expressed by alternating gabbro and pyroxeriite sequences. The gabbro and pyroxenite sequences consist of metre-thick differentiation units, which suggest emplacement by pe¬riodic injection of magma pulses as vertical dykes that amalgamated, similarly to a sub-volcanic sheeted dyke com¬plex. Individual dykes underwent internal differentiation following a solidification front (favoured by a significant lateral/horizontal thermal gradient) parallel to the dyke edges. Pyroxenitic layers result from the fractionation and accumulation of clinopyroxene ± olivine ± plagioclase crystals from a mildly alkaline basaltic liquid and are interpre¬ted as truncated differentiation sequences, from which residual melts were extracted by compaction. Clinopyroxene mineral orientation in pyroxenites (evidenced by EBSD and micro X-ray tomography analysis) display a marked pure shear component, supporting this interpretation. Compaction and squeezing of the crystal mush is ascribed to the incoming and inflating magma pulses. The resulting expelled interstitial liquid was likely collected and erupted along with the magma flowing through the newly injected dykes. Gabbro sequences represent crystallised coalesced magma batches, emplaced at lower rates at the end of eruptive cycles, and underwent minor melt extraction as evi¬denced by clinopyroxene orientations that record a simple shear component suggesting syn-magmatic deformation parallel to observed NNF.-SSW trending shear-zones induced by the regional tensional Miocene stress-field. The initiation and geometry of PX1 is controlled by the regional extensional tectonic regime whereas rates and vo¬lumes of magma depend on source-related factors. High injection rates are likely to induce intrusion growth rates larger than could be accommodated by the regional extension. In this case, dyke tip geometry and the inability of magma to circulate through previously emplaced and crystallised dykes could result in an increase of non-lithostatic pressure on previously emplaced mushy dyke walls; generating strong pure-shear compaction and interstitial melt expulsion within the feeder-zone as recorded by the cumulitic pyroxenite bands and anorthositic collection zones. The whole-rock major and trace-element chemistry of PX1 gabbros and pyroxenites is globally homogeneous and controlled by the cumulate nature of the samples (i.e. on the modal proportions of olivine, pyroxene, plagioclase and oxides). However, small variations of whole-rock trace-element contents as well as trace-element contents of clinopyroxene rims suggest that in-situ re-equilibration and crystallisation has occurred. Additionally, the global homogeneity and presence of complex zoning of rare resorbed clinopyroxene crystals suggest that the PX1 feeder- zone overlies a periodically replenished and efficiently mixed magma chamber. Each individual dyke of magma thus originated from a compositionally constant mildly alkaline magma and differentiated independently from the others reaching up to 70% fractionation. Following dyke arrest these are affected by interaction with the trapped interstitial liquid prior to its compaction-linked expulsion (thus stopping the differentiation process). This emplacement model implies that minimum amount of approximately 150 km3 of magma is needed to generate PX1, part of it having been erupted through the overlying Central Volcanic Complex of Fuerteventura. The radiogenic isotope ratios of PX1 samples reveal the contribution on three end-members during magma genesis. This mixing of the H1MU, EMI and DMM end-members could reflect the interaction of the deep-seated Canarian mantle plume with a heterogeneous metasomatic and sepentininsed lithospheric mantle. Additionally, the observed trace-element and isotopic variations within the same fades groups could reflect varying degrees of partial melting of the source region, thus tapping more or less large areas of the metasomatised lithospheric mantle during interac¬tion with the plume. High precision ID-TIMS U/Pb zircon and baddeleyite ages from the PX1 gabbro samples, indicate initiation of magma crystallisation at 22.10 ± 0.07 Ma. The magmatic activity lasted a minimum of 0.48 to 0.52 Ma. 40Ar/39Ar amphibole ages are of 21.9 ± 0.6 to 21.8 ± 0.3, identical within errors to the U/Pb ages. The combination of the 40Ar/39Ar and U/Pb datasets imply that the maximum amount of time PX1 took to cool below amphibole Tc is 0.8 Ma, suggesting PX1 lifetime of 520 000 to 800 000 years. On top of this, the coexistence of baddeleyite and zircon in a single sample is ascribed to the interaction of PX1 with C02-rich carbonatite-derived fluids released from the host-rock carbonatites during contact metamorphism 160 000 years after PX1 initiation. These ages are in agreement with the emplacement model, implying a crystallisation time of less than 1 to 5 years for individual dykes.

Sediment budget monitoring of a debris-flow torrent (French Prealps)

The Manival near Grenoble (French Prealps) is a very active debris-flow torrent equipped with a large sediment trap (25 000 m3) protecting an urbanized alluvial fan from debris-flows. We began monitoring the sediment budget of the catchment controlled by the trap in Spring 2009. Terrestrial laser scanner is used for monitoring topographic changes in a small gully, the main channel, and the sediment trap. In the main channel, 39 cross-sections are surveyed after every event. Three periods of intense geomorphic activity are documented here. The first was induced by a convective storm in August 2009 which triggered a debris-flow that deposited ~1,800 m3 of sediment in the trap. The debris-flow originated in the upper reach of the main channel and our observations showed that sediment outputs were entirely supplied by channel scouring. Hillslope debris-flows were initiated on talus slopes, as revealed by terrestrial LiDAR resurveys; however they were disconnected to the main channel. The second and third periods of geomorphic activity were induced by long duration and low intensity rainfall events in September and October 2009 which generate small flow events with intense bedload transport. These events contribute to recharge the debris-flow channel with sediments by depositing important gravel dunes propagating from headwaters. The total recharge in the torrent subsequent to bedload transport events was estimated at 34% of the sediment erosion induced by the August debris-flow.

Sediment budget monitoring of debris-flow and bedload transport in the Manival Torrent, SE France

Steep mountain catchments typically experience large sediment pulses from hillslopes which are stored in headwater channels and remobilized by debris-flows or bedload transport. Event-based sediment budget monitoring in the active Manival debris-flow torrent in the French Alps during a two-year period gave insights into the catchment-scale sediment routing during moderate rainfall intensities which occur several times each year. The monitoring was based on intensive topographic resurveys of low- and high-order channels using different techniques (cross-section surveys with total station and high-resolution channel surveys with terrestrial and airborne laser scanning). Data on sediment output volumes from the main channel were obtained by a sediment trap. Two debris-flows were observed, as well as several bedload transport flow events. Sediment budget analysis of the two debris-flows revealed that most of the debris-flow volumes were supplied by channel scouring (more than 92%). Bedload transport during autumn contributed to the sediment recharge of high-order channels by the deposition of large gravel wedges. This process is recognized as being fundamental for debris-flow occurrence during the subsequent spring and summer. A time shift of scour-and-fill sequences was observed between low- and high-order channels, revealing the discontinuous sediment transfer in the catchment during common flow events. A conceptual model of sediment routing for different event magnitude is proposed.

Application des méthodes géophysiques à l'étude de deux sites (Kappelen et Grenchen) de l'aquifère poreux complexe du Seeland

L'aquifère du Seeland représente une richesse en ressources hydriques qu'il est impératif de préserver contre tout risque de détérioration. Cet aquifère prolifique est constitué principalement de sédiments alluviaux post-glaciaires (graviers, sables et limons). Il est soumis aux contraintes environnementales des pratiques d'agriculture intensive, du réseau routier, des villes et de leurs activités industrielles. La connaissance optimale de ces ressources est donc primordiale pour leur protection. Dans cette optique, deux sites Kappelen et Grenchen représentatifs de l'aquifère du Seeland ont été étudiés. L'objectif de ce travail est de caractériser d'un point de vue hydrogéophysique l'aquifère au niveau de ces deux sites, c'est-à-dire, comprendre la dynamique des écoulements souterrains par l'application des méthodes électriques de surface associées aux diagraphies en intégrant des méthodes hydrogéologiques. Pour le site de Kappelen, les méthodes électriques de surface ont permis d'identifier les différents faciès géoélectriques en présence et de mettre en évidence leur disposition en une structure tabulaire et horizontale. Il s'agit d'un aquifère libre constitué d'une série de graviers allant jusqu'à 15 m de profondeur reposant sur de la moraine argileuse. Les diagraphies électriques, nucléaires et du fluide ont servis à la détermination des caractéristiques pétrophysiques et hydrauliques de l'aquifère qui contrôlent son comportement hydrodynamique. Les graviers aquifères de Kappelen présentent deux minéraux dominants: quartz et calcite. Les analyses minéralogiques indiquent que ces deux éléments constituent 65 à 75% de la matrice. La porosité totale obtenue par les diagraphies nucléaires varie de 20 à 30 %, et de 22 à 29 % par diagraphies électrique. Avec les faibles valeurs de Gamma Ray ces résultats indiquent que l'aquifère des graviers de Kappelen est dépourvu d'argile minéralogique. La perméabilité obtenue par diagraphies du fluide varie de 3.10-4 à 5.10-2 m/s, et par essais de pompage de 10-4 à 10-2 m/s. Les résultats des analyses granulométriques indiquent une hétérogénéité granulométrique au niveau des graviers aquifères. La fraction de sables, sables très fins, silts et limons constitue de 10 à 40 %. Ces éléments jouent un rôle important dans le comportement hydraulique de l'aquifère. La porosité efficace de 11 à 25% estimée à partir des résultats des analyses granulométriques suppose que les zones les plus perméables correspondent aux zones les plus graveleuses du site. Etablie sur le site de Kappelen, cette méthodologie a été utilisée sur le site de Grenchen. Les méthodes électriques de surface indiquent que l'aquifère captif de Grenchen est constitué des sables silteux comprenant des passages sableux, encadrés par des silts argileux imperméables. L'aquifère de Grenchen est disposé dans une structure relativement tabulaire et horizontale. Son épaisseur totale peut atteindre les 25 m vers le sud et le sud ouest ou les passages sableux sont les plus importants. La détermination des caractéristiques pétrophysiques et hydrauliques s'est faite à l'aide des diagraphies. Les intensités Gamma Ray varient de 30 à 100 cps, les plus fortes valeurs n'indiquent qu'une présence d'éléments argileux mais pas de bancs d'argile. Les porosités totales de 15 à 25% et les densités globales de 2.25 à 2.45 g/cm3 indiquent que la phase minérale (matrice) est composée essentiellement de quartz et de calcaire. Les densités de matrice varient entre 2.65 et 2.75 g/cm3. La perméabilité varie de 2 10-6 à 5 10-4 m/s. La surestimation des porosités totales à partir des diagraphies électriques de 25 à 42% est due à la présence d'argiles. -- The vast alluvial Seeland aquifer system in northwestern Switzerland is subjected to environmental challenges due to intensive agriculture, roads, cities and industrial activities. Optimal knowledge of the hydrological resources of this aquifer system is therefore important for their protection. Two representative sites, Kappelen and Grenchen, of the Seeland aquifer were investigated using surface-based geoelectric methods and geophysical borehole logging methods. By integrating of hydrogeological and hydrogeophysical methods, a reliable characterization of the aquifer system at these two sites can be performed in order to better understand the governing flow and transport process. At the Kappelen site, surface-based geoelectric methods allowed to identify various geoelectric facies and highlighted their tabular and horizontal structure. It is an unconfined aquifer made up of 15 m thick gravels with an important sandy fraction and bounded by a shaly glacial aquitard. Electrical and nuclear logging measurements allow for constraining the petrophysical and hydrological parameters of saturated gravels. Results indicate that in agreement with mineralogical analyses, matrix of the probed formations is dominated by quartz and calcite with densities of 2.65 and 2.71 g/cc, respectively. These two minerals constitute approximately 65 to 75 % of the mineral matrix. Matrix density values vary from 2.60 to 2.75 g/cc. Total porosity values obtained from nuclear logs range from 20 to 30 % and are consistent with those obtained from electrical logs ranging from 22 to 29 %. Together with the inherently low natural gamma radiation and the matrix density values obtained from other nuclear logging measurements, this indicates that at Kappelen site the aquifer is essentially devoid of clay. Hydraulic conductivity values obtained by the Dilution Technique vary between 3.10-4 and 5.10-2 m/s, while pumping tests give values ranging from 10-4 to 10-2 m/s. Grain size analysis of gravel samples collected from boreholes cores reveal significant granulometric heterogeneity of these deposits. Calculations based on these granulometric data have shown that the sand-, silt- and shale-sized fractions constitute between 10 and 40 % of the sample mass. The presence of these fine elements in general and their spatial distribution in particular are important as they largely control the distribution of the total and effective porosity as well as the hydraulic conductivity. Effective porosity values ranging from 11 to 25% estimated from grain size analyses indicate that the zones of higher hydraulic conductivity values correspond to the zones dominated by gravels. The methodology established at the Kappelen site was then applied to the Grenchen site. Results from surface-based geoelectric measurements indicate that it is a confined aquifer made up predominantly of shaly sands with intercalated sand lenses confined impermeable shally clay. The Grenchen confined aquifer has a relatively tabular and horizontal structure with a maximum thickness of 25 m in the south and the southwest with important sand passages. Petrophysical and hydrological characteristics were performed using electrical and nuclear logging. Natural gamma radiation values ranging from 30 to 100 cps indicate presence of a clay fraction but not of pure clay layers. Total porosity values obtained from electrical logs vary form 25 to 42%, whereas those obtained from nuclear logs values vary from 15 to 25%. This over-estimation confirms presences of clays. Density values obtained from nuclear logs varying from 2.25 to 2.45 g/cc in conjunction with the total porosity values indicate that the dominating matrix minerals are quartz and calcite. Matrix density values vary between 2.65 and 2.75 g/cc. Hydraulic conductivity values obtained by the Dilution Technique vary from 2 10-6 to 5 10-4 m/s.

Presència i distribució espacial de concentracions anormals d'ió nitrat a les aigües subterrànies dels aqüífers del Baix Fluvià (Alt Empordà)

Presentem una caracterització de la distribució espacial de la contaminació per nitrats en els aqüífers de la plana del baix Fluvià. En concret, els valors de concentració més elevats s’han observat a l’aqüífer superficial, on sobrepassen, amb escreix, els límits de potabilitat aconsellats pel RD 1138/1990. L’àrea més afectada ocupa part dels termes municipals de Vilamacolum i Torroella de Fluvià, en els quals hi predominen concentracions que oscil·len entre 50 i 200 mg × L–1. Els focus de contaminació són de tipus puntual i difús. Documentada la presència anormal de nitrats en aquest aqüífer, es planteja el correcte maneig agrícola, en vistes a una doble finalitat: d’una banda, contribuir a la planificació d’unes bones pràctiques agrícoles tenint en compte la relació entre volums de reg i aportacions de nitrogen; d’altra banda, preveure una progressiva recuperació de la qualitat de l’aqüífer, a partir de l’ús de tecnologies que integrin la ubicació d’unes zones favorables on es pugui maximitzar l’extracció del recurs contaminat, i la ubicació d’uns sòls aptes per a l’aplicació d’aquest recurs, amb uns criteris agronòmics coherents.

Pohjavesihanke ja sen toteutus

Pohjavesihankkeen onnistunut toteutus vaatii tietoa useilta eri aloilta. Perustieto pohjavesistä, niiden muodostumisesta ja esiintymisestä, pohjavesiä koskevan lainsäädännön ja vedenottoon tarvittavien lupien tuntemus, pohjavesitutkimuksen eri vaiheiden hallitseminen ja tieto näiden vaiheiden kestoon, aikatauluun ja kustannuksiin vaikuttavista tekijöistä ovat kaikki onnistuneen hankkeen edellytyksiä. Tässä diplomityössä selvitetään pohjavesihankkeen eri vaiheet sekä niiden sisältö ja tarvittavat toimenpiteet vesihuoltolaitoksen kannalta.Työssä tarkastellaan yleisesti kunnallista vesihuoltoa Suomessa ja Suomen pohjavesioloja. Lisäksi tarkastellaan pohjavesiä koskevaa lainsäädäntöä sekä pohjavesialueiden kartoituksen ja luokituksen nykytilannetta. Työssä käsitellään esiintymäkohtaisen pohjavesitutkimuksen eri vaiheet menetelmineen. Pohjavesitutkimusten ja pohjavedenoton vaatimat luvat käydään läpi.Pohjavesihanke on jaettu osiin, joiden kestoa ja niihin vaikuttavia tekijöitä tarkastellaan. Lisäksi pohjavesihankkeesta on tehty esimerkkiaikataulu, jossa eri työvaiheiden vaatimat ajat esitetään suhteessa toisiinsa. Hankkeen kustannukset on jaettu henkilöstö-, materiaali- ja muihin kustannuksiin. Näiden rakenne eri työvaiheissa on selvitetty.Joensuun Veden Paavonlammen alueella Kiihtelysvaaran kunnassa tekemän esimerkkihankkeen avulla havainnollistetaan sitä, millä tavoin käyttöön otettavan pohjavesialueen ominaisuudet sekä käytettävissä oleva työvoima ja kalusto vaikuttavat kustannuksiin, aikatauluun ja eri tutkimusvaiheiden tarpeellisuuteen.Pohjavesihanketta aloitettaessa ja sen aikana on kiinnitettävä erityisesti huomiota avoimeen tiedottamiseen vaikutusalueen asukkaille ja maanomistajille. Keskeistä työn edetessä on ottaa huomioon se, että pohjavesihankkeen toteutus etenee vaiheittain aikaisempien tulosten perusteella.

Conditioning of multiple-point statistics facies simulations to tomographic images

Geophysical tomography captures the spatial distribution of the underlying geophysical property at a relatively high resolution, but the tomographic images tend to be blurred representations of reality and generally fail to reproduce sharp interfaces. Such models may cause significant bias when taken as a basis for predictive flow and transport modeling and are unsuitable for uncertainty assessment. We present a methodology in which tomograms are used to condition multiple-point statistics (MPS) simulations. A large set of geologically reasonable facies realizations and their corresponding synthetically calculated cross-hole radar tomograms are used as a training image. The training image is scanned with a direct sampling algorithm for patterns in the conditioning tomogram, while accounting for the spatially varying resolution of the tomograms. In a post-processing step, only those conditional simulations that predicted the radar traveltimes within the expected data error levels are accepted. The methodology is demonstrated on a two-facies example featuring channels and an aquifer analog of alluvial sedimentary structures with five facies. For both cases, MPS simulations exhibit the sharp interfaces and the geological patterns found in the training image. Compared to unconditioned MPS simulations, the uncertainty in transport predictions is markedly decreased for simulations conditioned to tomograms. As an improvement to other approaches relying on classical smoothness-constrained geophysical tomography, the proposed method allows for: (1) reproduction of sharp interfaces, (2) incorporation of realistic geological constraints and (3) generation of multiple realizations that enables uncertainty assessment.