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.

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.

Two-dimensional probabilistic inversion of plane-wave electromagnetic data: Methodology, model constraints and joint inversion with electrical resistivity data

Probabilistic inversion methods based on Markov chain Monte Carlo (MCMC) simulation are well suited to quantify parameter and model uncertainty of nonlinear inverse problems. Yet, application of such methods to CPU-intensive forward models can be a daunting task, particularly if the parameter space is high dimensional. Here, we present a 2-D pixel-based MCMC inversion of plane-wave electromagnetic (EM) data. Using synthetic data, we investigate how model parameter uncertainty depends on model structure constraints using different norms of the likelihood function and the model constraints, and study the added benefits of joint inversion of EM and electrical resistivity tomography (ERT) data. Our results demonstrate that model structure constraints are necessary to stabilize the MCMC inversion results of a highly discretized model. These constraints decrease model parameter uncertainty and facilitate model interpretation. A drawback is that these constraints may lead to posterior distributions that do not fully include the true underlying model, because some of its features exhibit a low sensitivity to the EM data, and hence are difficult to resolve. This problem can be partly mitigated if the plane-wave EM data is augmented with ERT observations. The hierarchical Bayesian inverse formulation introduced and used herein is able to successfully recover the probabilistic properties of the measurement data errors and a model regularization weight. Application of the proposed inversion methodology to field data from an aquifer demonstrates that the posterior mean model realization is very similar to that derived from a deterministic inversion with similar model constraints.

Error models in hydrogeology applications

Notre consommation en eau souterraine, en particulier comme eau potable ou pour l'irrigation, a considérablement augmenté au cours des années. De nombreux problèmes font alors leur apparition, allant de la prospection de nouvelles ressources à la remédiation des aquifères pollués. Indépendamment du problème hydrogéologique considéré, le principal défi reste la caractérisation des propriétés du sous-sol. Une approche stochastique est alors nécessaire afin de représenter cette incertitude en considérant de multiples scénarios géologiques et en générant un grand nombre de réalisations géostatistiques. Nous rencontrons alors la principale limitation de ces approches qui est le coût de calcul dû à la simulation des processus d'écoulements complexes pour chacune de ces réalisations. Dans la première partie de la thèse, ce problème est investigué dans le contexte de propagation de l'incertitude, oú un ensemble de réalisations est identifié comme représentant les propriétés du sous-sol. Afin de propager cette incertitude à la quantité d'intérêt tout en limitant le coût de calcul, les méthodes actuelles font appel à des modèles d'écoulement approximés. Cela permet l'identification d'un sous-ensemble de réalisations représentant la variabilité de l'ensemble initial. Le modèle complexe d'écoulement est alors évalué uniquement pour ce sousensemble, et, sur la base de ces réponses complexes, l'inférence est faite. Notre objectif est d'améliorer la performance de cette approche en utilisant toute l'information à disposition. Pour cela, le sous-ensemble de réponses approximées et exactes est utilisé afin de construire un modèle d'erreur, qui sert ensuite à corriger le reste des réponses approximées et prédire la réponse du modèle complexe. Cette méthode permet de maximiser l'utilisation de l'information à disposition sans augmentation perceptible du temps de calcul. La propagation de l'incertitude est alors plus précise et plus robuste. La stratégie explorée dans le premier chapitre consiste à apprendre d'un sous-ensemble de réalisations la relation entre les modèles d'écoulement approximé et complexe. Dans la seconde partie de la thèse, cette méthodologie est formalisée mathématiquement en introduisant un modèle de régression entre les réponses fonctionnelles. Comme ce problème est mal posé, il est nécessaire d'en réduire la dimensionnalité. Dans cette optique, l'innovation du travail présenté provient de l'utilisation de l'analyse en composantes principales fonctionnelles (ACPF), qui non seulement effectue la réduction de dimensionnalités tout en maximisant l'information retenue, mais permet aussi de diagnostiquer la qualité du modèle d'erreur dans cet espace fonctionnel. La méthodologie proposée est appliquée à un problème de pollution par une phase liquide nonaqueuse et les résultats obtenus montrent que le modèle d'erreur permet une forte réduction du temps de calcul tout en estimant correctement l'incertitude. De plus, pour chaque réponse approximée, une prédiction de la réponse complexe est fournie par le modèle d'erreur. Le concept de modèle d'erreur fonctionnel est donc pertinent pour la propagation de l'incertitude, mais aussi pour les problèmes d'inférence bayésienne. Les méthodes de Monte Carlo par chaîne de Markov (MCMC) sont les algorithmes les plus communément utilisés afin de générer des réalisations géostatistiques en accord avec les observations. Cependant, ces méthodes souffrent d'un taux d'acceptation très bas pour les problèmes de grande dimensionnalité, résultant en un grand nombre de simulations d'écoulement gaspillées. Une approche en deux temps, le "MCMC en deux étapes", a été introduite afin d'éviter les simulations du modèle complexe inutiles par une évaluation préliminaire de la réalisation. Dans la troisième partie de la thèse, le modèle d'écoulement approximé couplé à un modèle d'erreur sert d'évaluation préliminaire pour le "MCMC en deux étapes". Nous démontrons une augmentation du taux d'acceptation par un facteur de 1.5 à 3 en comparaison avec une implémentation classique de MCMC. Une question reste sans réponse : comment choisir la taille de l'ensemble d'entrainement et comment identifier les réalisations permettant d'optimiser la construction du modèle d'erreur. Cela requiert une stratégie itérative afin que, à chaque nouvelle simulation d'écoulement, le modèle d'erreur soit amélioré en incorporant les nouvelles informations. Ceci est développé dans la quatrième partie de la thèse, oú cette méthodologie est appliquée à un problème d'intrusion saline dans un aquifère côtier. -- Our consumption of groundwater, in particular as drinking water and for irrigation, has considerably increased over the years and groundwater is becoming an increasingly scarce and endangered resource. Nofadays, we are facing many problems ranging from water prospection to sustainable management and remediation of polluted aquifers. Independently of the hydrogeological problem, the main challenge remains dealing with the incomplete knofledge of the underground properties. Stochastic approaches have been developed to represent this uncertainty by considering multiple geological scenarios and generating a large number of realizations. The main limitation of this approach is the computational cost associated with performing complex of simulations in each realization. In the first part of the thesis, we explore this issue in the context of uncertainty propagation, where an ensemble of geostatistical realizations is identified as representative of the subsurface uncertainty. To propagate this lack of knofledge to the quantity of interest (e.g., the concentration of pollutant in extracted water), it is necessary to evaluate the of response of each realization. Due to computational constraints, state-of-the-art methods make use of approximate of simulation, to identify a subset of realizations that represents the variability of the ensemble. The complex and computationally heavy of model is then run for this subset based on which inference is made. Our objective is to increase the performance of this approach by using all of the available information and not solely the subset of exact responses. Two error models are proposed to correct the approximate responses follofing a machine learning approach. For the subset identified by a classical approach (here the distance kernel method) both the approximate and the exact responses are knofn. This information is used to construct an error model and correct the ensemble of approximate responses to predict the "expected" responses of the exact model. The proposed methodology makes use of all the available information without perceptible additional computational costs and leads to an increase in accuracy and robustness of the uncertainty propagation. The strategy explored in the first chapter consists in learning from a subset of realizations the relationship between proxy and exact curves. In the second part of this thesis, the strategy is formalized in a rigorous mathematical framework by defining a regression model between functions. As this problem is ill-posed, it is necessary to reduce its dimensionality. The novelty of the work comes from the use of functional principal component analysis (FPCA), which not only performs the dimensionality reduction while maximizing the retained information, but also allofs a diagnostic of the quality of the error model in the functional space. The proposed methodology is applied to a pollution problem by a non-aqueous phase-liquid. The error model allofs a strong reduction of the computational cost while providing a good estimate of the uncertainty. The individual correction of the proxy response by the error model leads to an excellent prediction of the exact response, opening the door to many applications. The concept of functional error model is useful not only in the context of uncertainty propagation, but also, and maybe even more so, to perform Bayesian inference. Monte Carlo Markov Chain (MCMC) algorithms are the most common choice to ensure that the generated realizations are sampled in accordance with the observations. Hofever, this approach suffers from lof acceptance rate in high dimensional problems, resulting in a large number of wasted of simulations. This led to the introduction of two-stage MCMC, where the computational cost is decreased by avoiding unnecessary simulation of the exact of thanks to a preliminary evaluation of the proposal. In the third part of the thesis, a proxy is coupled to an error model to provide an approximate response for the two-stage MCMC set-up. We demonstrate an increase in acceptance rate by a factor three with respect to one-stage MCMC results. An open question remains: hof do we choose the size of the learning set and identify the realizations to optimize the construction of the error model. This requires devising an iterative strategy to construct the error model, such that, as new of simulations are performed, the error model is iteratively improved by incorporating the new information. This is discussed in the fourth part of the thesis, in which we apply this methodology to a problem of saline intrusion in a coastal aquifer.

O impacto industrial na composição química das águas subterrâneas com enfoque de consumo humano (Rio Grande, RS)

The alteration in the quality of groundwater in the industrial zone of Brazil's southernmost state was assessed by a hydrogeochemical study. In 70-85% of the studied area, fluoride, nitrate and phosphate in groundwater surpass the maximum limits for human consumption according to Brazilian environmental legislation. The chemical spectrum of contaminants and their spatial distribution show that fertilizer production processes are responsible for groundwater pollution. The natural conditions of the region are not favorable for minimal protection against infiltration of pollutants into the aquifer.

Active sulphide mine tailings impoundments as sources of contaminated drainage: controlling factors, methods of characterisation and geochemical constraints for mitigation

Low quality mine drainage from tailings facilities persists as one of the most significant global environmental concerns related to sulphide mining. Due to the large variation in geological and environmental conditions at mine sites, universal approaches to the management of mine drainage are not always applicable. Instead, site-specific knowledge of the geochemical behaviour of waste materials is required for the design and closure of the facilities. In this thesis, tailings-derived water contamination and factors causing the pollution were investigated in two coeval active sulphide mine sites in Finland: the Hitura Ni mine and the Luikonlahti Cu-Zn-Co-Ni mine and talc processing plant. A hydrogeochemical study was performed to characterise the tailingsderived water pollution at Hitura. Geochemical changes in the Hitura tailings were evaluated with a detailed mineralogical and geochemical investigation (solid-phase speciation, acid mine drainage potential, pore water chemistry) and using a spatial assessment to identify the mechanisms of water contamination. A similar spatial investigation, applying selective extractions, was carried out in the Luikonlahti tailings area for comparative purposes (Hitura low-sulphide tailings vs. Luikonlahti sulphide-rich tailings). At both sites, hydrogeochemistry of tailings seepage waters was further characterised to examine the net results of the processes observed within the impoundments and to identify constraints for water treatment. At Luikonlahti, annual and seasonal variation in effluent quality was evaluated based on a four-year monitoring period. Observations pertinent to future assessment and mine drainage prevention from existing and future tailings facilities were presented based on the results. A combination of hydrogeochemical approaches provided a means to delineate the tailings-derived neutral mine drainage at Hitura. Tailings effluents with elevated Ni, SO₄ ^2- and Fe content had dispersed to the surrounding aquifer through a levelled-out esker and underneath the seepage collection ditches. In future mines, this could be avoided with additional basal liners in tailings impoundments where the permeability of the underlying Quaternary deposits is inadequate, and with sufficiently deep ditches. Based on the studies, extensive sulphide oxidation with subsequent metal release may already initiate during active tailings disposal. The intensity and onset of oxidation depended on e.g. the Fe sulphide content of the tailings, water saturation level, and time of exposure of fresh sulphide grains. Continuous disposal decreased sulphide weathering in the surface of low-sulphide tailings, but oxidation initiated if they were left uncovered after disposal ceased. In the sulphide-rich tailings, delayed burial of the unsaturated tailings had resulted in thick oxidized layers, despite the continuous operation. Sulphide weathering and contaminant release occurred also in the border zones. Based on the results, the prevention of sulphide oxidation should already be considered in the planning of tailings disposal, taking into account the border zones. Moreover, even lowsulphide tailings should be covered without delay after active disposal ceases. The quality of tailings effluents showed wide variation within a single impoundment and between the two different types of tailings facilities assessed. The affecting factors included source materials, the intensity of weathering of tailings and embankment materials along the seepage flow path, inputs from the process waters, the water retention time in tailings, and climatic seasonality. In addition, modifications to the tailings impoundment may markedly change the effluent quality. The wide variation in the tailings effluent quality poses challenges for treatment design. The final decision on water management requires quantification of the spatial and seasonal fluctuation at the site, taking into account changes resulting from the eventual closure of the impoundment. Overall, comprehensive hydrogeochemical mapping was deemed essential in the identification of critical contaminants and their sources at mine sites. Mineralogical analysis, selective extractions, and pore water analysis were a good combination of methods for studying the weathering of tailings and in evaluating metal mobility from the facilities. Selective extractions with visual observations and pH measurements of tailings solids were, nevertheless, adequate in describing the spatial distribution of sulphide oxidation in tailings impoundments. Seepage water chemistry provided additional data on geochemical processes in tailings and was necessary for defining constraints for water treatment.

Diagnóstico espaço-temporal da ocorrência de herbicidas nas águas superficiais e sedimentos do Rio Corumbataí e principais afluentes

Residues of herbicides from sugarcane were monitored in waters and sediments of Corumbataí River and tributaries. Ametryne, atrazine, simazine, hexazinone, glyphosate, and clomazone were detected in water samples, with negligible levels of ametryne and glyphosate in sediment samples. The area of recharge of the Guarani aquifer presented the highest triazine and clomazone levels. The triazines were detected at higher levels, with atrazine above Brazil's potability and quality standards. Total herbicide levels at some sampling points were 13 times higher than the European Community potability limit. There is no Brazilian standard for ametryne, although the risk is larger due to ametryne's higher toxicity for the aquatic biota.

Teores de nitrato (NO3-) e amônio (NH4+) nas águas do aqüífero Barreiras nos bairros do Reduto, Nazaré e Umarizal - Belém/PA

This paper evaluates the occurrence of nitrate and ammonium in the Barreiras aquifer in the metropolitan area of Belém, Pará State, Brazil. The results show that some wells display ammonium and nitrate concentrations above or close to the limits of water potability regulations. The main cause of the pollution of these waters is related the local disposal of domestic effluents and/or leakage from pipes of old sewage systems.

Analytical procedures for determining Pb and Sr isotopic compositions in water samples by ID-TIMS

Few articles deal with lead and strontium isotopic analysis of water samples. The aim of this study was to define the chemical procedures for Pb and Sr isotopic analyses of groundwater samples from an urban sedimentary aquifer. Thirty lead and fourteen strontium isotopic analyses were performed to test different analytical procedures. Pb and Sr isotopic ratios as well as Sr concentration did not vary using different chemical procedures. However, the Pb concentrations were very dependent on the different procedures. Therefore, the choice of the best analytical procedure was based on the Pb results, which indicated a higher reproducibility from samples that had been filtered and acidified before the evaporation, had their residues totally dissolved, and were purified by ion chromatography using the Biorad® column. Our results showed no changes in Pb ratios with the storage time.

Comportamento do herbicida hexazinone em área de recarga do aqüífero Guarani cultivada com cana-de-açúcar

São Paulo state (Brazil) has an important area of sugarcane production, mainly for obtaining alcohol and sugar, where there is an intensive use of pesticides. An important recharge zone of Guarani aquifer, with supplies water for the local population, is located at Ribeirão Preto city, so the local behavior of pesticides must be investigated. The GUS index was obtained by using the paramenters Koc and half-life for hezazinone herbicide, determinated in representative soil of this region. This study has demonstrated that there is potential risks of hexazinone leaching to ground water, indicating that this herbicide must be monitored in ground water.

Mikkelin Pursialan pohjavesialueen haitta-aineet

Mikkelin talousvedestä kahden kolmasosan tullessa Pursialan pohjavesialueelta on alueen suojeleminen tärkeää. Pohjaveden laatua uhkaavat etenkin alueella sattuneet pohjavedenpilaantumistapaukset. Merkittävimmät pohjaveden pilaantumistapaukset ovat VAPO Oy:n sahan aiheuttama pohjaveden pilaantuminen kloorifenoleilla (CP) ja VR:n ratapölkkykyllästämön aiheuttama pohjaveden pilaantuminen kreosoottiöljyllä sekä Rinnekadun Nesteen aiheuttama pohjaveden pilaantuminen MTBE:llä. Alueella on tehty tutkimuksia ja kunnostuksia pilaantumiin liittyen, mutta näiden tuloksia ei ole aikaisemmin koottu yhteen. Tämän työn tavoitteena oli koota tulokset samaan aineistoon. Työssä keskityttiin kloorifenolien leviämisen tarkasteluun sen Pursialan pohjavedenottamolle muodostaman suurimman uhan vuoksi. Kallioperätietojen, maanpintatietojen ja näytetietojen pohjalta laadittiin myös pienoismalli CP-pilaantuman leviämisen kokonaiskuvan hahmottamiseksi. Työn tavoitteena oli lisäksi tehdä riskitarkastelua CP-pilaantumaan liittyen ja etsiä keinoja hallita havaittuja riskejä. Riskinhallintaan liittyen työssä tutkittiin kloorifenoleilla pilaantuneen alueen maaperä- ja kalliotietoja sekä pohjaveden laatutietoja. Pursialan pohjavedessä on runsaasti rautaa ja mangaania sekä aggressiivista hiilihappoa. Pohjaveden pH on alueella noin 6,5, lämpötila noin 7,5 ºC ja happipitoisuus noin 0,7 mg/l. Pursialan kaupunginalueen kallioperässä on havaittavissa VAPO Oy:n sahalta vedenottamolle etenevä kalliopainanne, jota pitkin CP etenee. Alueen kallioperä on kiillegneissiä, jossa on pohjois–etelä-suuntaista rakoilua. Maaperätuloksien perusteella on havaittavissa vettä hyvin johtavien maakerrosten jatkuminen koko vedenottamon ja sahan välisen matkan, mikä tarkoittaa, että CP-pitoisella pohjavedellä voi olla aiemmin oletettua nopeampikin yhteys sahalta vedenottamolle. Suurin CP-pitoisuus noin 100 000 µg/l on mitattu KY-5-altaan kohdalle asennetun M14-pohjavesiputken pohjasta. Talousvesiasetuksen raja-arvo CP:lle on 10 µg/l. Sahan ja vedenottamon puolivälissä on havaittu yli 10 000 µg/l meneviä CP-pitoisuuksia. Suurin vedenottamon kaivoista (kaivo 10) mitattu pitoisuus on 149 µg/l. Jakotukilta raakavedestä otetuissa näytteissä tai talousvedessä ei ole kuitenkaan havaittu talousvesiasetuksen ylittäviä CP-pitoisuuksia. Pienoismallin perusteella CP sijaitsee sahan alueella lähellä kallionpintaa ja hajaantuu koko pohjavesipatjaan vedenottamolle päin mentäessä. CP-mittaustuloksissa on havaittavissa pulssimaisuutta. Tämä johtuu todennäköisesti Saimaan pinnan vaihtelun seurauksena muuttuvasta rantaimeytyneen pohjaveden määrästä. Saimaan pinnan nousu näyttäisi tuloksien perusteella nostavan CP-pitoisuuksia saha-alueella ja laskevan lähellä vedenottamoa. Pohjaveden pintatietojen perusteella tehdyn tarkastelun mukaan pohjavesi voi kulkeutua sahalta vedenottamolle parhaimmillaan noin vuodessa. Työssä arvioitiin KY-5–liuoksen vuosittaiseksi käyttömääräksi noin 648–970 m3. Allassakkaa arvioitiin syntyneen yhteensä noin 10–31 m3. Pohjaveteen arvioitiin joutuneen toiminnan aikana yhteensä noin 3 000–4 000 kg CP:tä. Kloorifenolit esiintyvät pohjavedessä lähes täysin kloorifenolaatteina. Kloorifenolien hajoaminen ja muuntuminen pohjavedessä on epätodennäköistä. Käsitteellisen mallin mukaan kloorifenolipilaantuman suurimmat riskit aiheutuvat kloorifenolien mahdollisuudesta pilata Pursialan vedenottamon talousvesi. Tällä hetkellä riskejä hallitaan kloorifenolien leviämisen tarkkailulla, sahan ja vedenottamon puolivälissä sijaitsevalla koepumppauksella sekä varautumalla aktiivihiilijauheen syöttöön talousvesiprosessiin. Koepumppauksen avulla on saatu ylös tällä hetkellä noin 69 kg kloorifenoleita. Tutkimuksen perusteella suositeltavimmat riskinhallintatoimet tulevaisuudessa ovat sahalla sijaitseva kunnostuspumppaus, sahan ja vedenottamon väliin sijoittuva suojapumppaus- ja vesiverhoyhdistelmä sekä sahan rannan kautta tapahtuvan rantaimeytymisen estäminen.

Hydrographic and hydrogeological basin of Entre-Ribeiros: probable recharge zone delimitation and environmental impact assessment

The probable recharge zone delimitation of Entre-Ribeiros Basin (Northwest of the state of Minas Gerais / Brazil) is proposed in this study. The delimitation is based upon stratigraphy, geomorphology, geo-environmental domains and hydrogeology studies. Combining the recharge zone map with the land use variation between 1975 and 2008, the occupation trends of possible recharge zones are identified. Concluding, the environmental impacts for this basin are discussed.

Groundwater recharge estimate at Alto Rio Grande - MG watershed

Springs are outcrops of aquifers surface, and the water cycle in this environment pass through the recharge, generally defined as the amount of water added to the aquifer, which may occur locally from rainwater infiltration. This study uses the Water Table Fluctuation (WTF) method to estimate the direct recharge and a groundwater balance to estimate the deep recharge on unconfined aquifers. The WTF method employs data of the aquifer water levels and its specific yield to estimate the direct recharge. The groundwater balance considers the direct recharge estimated by the WTF method, as the water input in the system and outputs as the base flow and deep recharge. The recharge was estimated at four areas at the watershed of Alto Rio Grande city, Minas Gerais (MG) state, in Brazil. The direct recharge estimate was 121.11; 64.62; 83.99; 152.46 (mm/year) for the L1, L2, M1 and M2 areas. The effect of the presence of forest in the recharge area can prevail over slope of relief, allowing more direct recharge, even in sources with steeper relief. The runoff from the springs in the study period exceeded the direct recharge, indicating a situation in which the saturated zone feeds the vadose zone. The annual flow was above the direct recharge pointing a situation of over exploitation of the aquifer, a non sustainable situation. The specific yield of the aquifers could also have been underestimated.

Estimation of alluvial recharge in the semiarid

In areas where there is irrigated agriculture, the recuperation of water reserves in alluvial aquifers may occur preferentially due to precipitation. Recharging can be evaluated from variation information of water depth measured in piezometers or observation wells. Thus, the aim of this research is to study the recharge in the alluvial aquifer formed by the Mimoso temporary stream in the semiarid region of Pernambuco (PE), Brazil, using the method of the fluctuation of the water level. This system is typical on the Brazilian Northeast semiarid region, using groundwater for domestic supply and for irrigation on small scale agriculture. Monthly potentiometric levels and rainfall data were used. The selected period for the study, from January 2002 to October 2009, involved extreme events of flooding and droughts as well as regular years, providing a better understanding of the behavior of the alluvial recharge. It was found that the system responds significantly to precipitation events. It was also observed that even with different soil textures in the study area, recharge factors were not significantly different. The study provided a better understanding of the behavior of aquifer recharge and its relationship with the soil and the rainfall events in the region.