Using a comparative spatial and temporal analysis on sediment cores from 8 lakes in the Mackenzie River Delta uplands region, NT, Canada, this study assessed how persistent organic pollutant (POP) deposition to lake sediments was affected by: (1) the presence of retrogressive thaw slumps on lake shores; and (2) changes occurring with increased autochthonous primary productivity. POPs examined included polychlorinated biphenyls (PCBs), penta- and hexachlorobenzenes (CBzs), and dichlorodiphenyltrichloroethane and metabolites (DDTs). Surface sediments of slump-affected lakes contained higher total organic carbon (TOC)-normalized POP concentrations than nearby reference lakes unaffected by thaw slumps. Inorganic sedimentation rates were positively related to contaminant concentrations, suggesting that the influx of siliciclastic material reducing organic carbon in slump-affected lake water indirectly results in higher concentrations of POPs on sedimentary organic matter. This explanation was corroborated by an inverse relationship between sedimentary POP concentrations and TOC content of the lake water. Deposition proxies of autochthonous carbon were not significantly correlated to POP fluxes of surface sediments, and historical profile fluctuations did not coincide with variation in POP deposition. Thus this study does not support the contention that algal-derived organic carbon increases the delivery of organic pollutants to sediments (the algal-scavenging hypothesis), as previously proposed for mercury. Higher POP concentrations observed in surface sediments of slump-affected lakes are best explained by simple solvent switching processes of hydrophobic contaminants onto a lower pool of available organic carbon when compared to neighbouring lakes unaffected by thaw slump development.

Les travaux de cette thèse portent sur des méthodes de volumes finis sur maillages quelconque pour la discrétisation de problèmes d'évolution non linéaires modélisant le transport de contaminants en milieu poreux et les écoulements diphasiques.Au Chapitre 1, nous étudions une famille de schémas numériques pour la discrétisation d'une équation parabolique dégénérée de convection-reaction-diffusion modélisant le transport de contaminants dans un milieu poreux qui peut être hétérogène et anisotrope. La discrétisation du terme de diffusion est basée sur une famille de méthodes qui regroupe les schémas de volumes finis hybrides, de différences finies mimétiques et de volumes finis mixtes. Le terme de convection est traité à l'aide d'une famille de méthodes qui s'appuient sur les inconnues hybrides associées aux interfaces du maillage. Cette famille contient à la fois les schémas centré et amont. Les schémas que nous étudions permettent une discrétisation localement conservative des termes d'ordre un et d'ordre deux sur des maillages arbitraires en dimensions d'espace deux et trois. Nous démontrons qu'il existe une solution unique du problème discret qui converge vers la solution du problème continu et nous présentons des résultats numériques en dimensions d'espace deux et trois, en nous appuyant sur des maillages adaptatifs.Au Chapitre 2, nous proposons un schéma de volumes finis hybrides pour la discrétisation d'un problème d'écoulement diphasique incompressible et immiscible en milieu poreux. On suppose que ce problème a la forme d'une équation parabolique dégénérée de convection-diffusion en saturation couplée à une équation uniformément elliptique en pression. On considère un schéma implicite en temps, où les flux diffusifs sont discrétisés par la méthode des volumes finis hybride, ce qui permet de pouvoir traiter le cas d'un tenseur de perméabilité anisotrope et hétérogène sur un maillage très général, et l'on s'appuie sur un schéma de Godunov pour la discrétisation des flux convectifs, qui peuvent être non monotones et discontinus par rapport aux variables spatiales. On démontre l'existence d'une solution discrète, dont une sous-suite converge vers une solution faible du problème continu. On présente finalement des cas test bidimensionnels.Le Chapitre 3 porte sur un problème d'écoulement diphasique, dans lequel la courbe de pression capillaire admet des discontinuité spatiales. Plus précisément on suppose que l'écoulement prend place dans deux régions du sol aux propriétés très différentes, et l'on suppose que la loi de pression capillaire est discontinue en espace à la frontière entre les deux régions, si bien que la saturation de l'huile et la pression globale sont discontinues à travers cette frontière avec des conditions de raccord non linéaires à l'interface. On discrétise le problème à l'aide d'un schéma, qui coïncide avec un schéma de volumes finis standard dans chacune des deux régions, et on démontre la convergence d'une solution approchée vers une solution faible du problème continu. Les test numériques présentés à la fin du chapitre montrent que le schéma permet de reproduire le phénomène de piégeage de la phase huile. / In Chapter 1 we study a family of finite volume schemes for the numerical solution of degenerate parabolic convection-reaction-diffusion equations modeling contaminant transport in porous media. The discretization of possibly anisotropic and heterogeneous diffusion terms is based upon a family of numerical schemes, which include the hybrid finite volume scheme, the mimetic finite difference scheme and the mixed finite volume scheme. One discretizes the convection term by means of a family of schemes which makes use of the discrete unknowns associated to the mesh interfaces, and contains as special cases an upwind scheme and a centered scheme. The numerical schemes which we study are locally conservative and allow computations on general multi-dimensional meshes. We prove that the unique discrete solution converges to the unique weak solution of the continuous problem. We also investigate the solvability of the linearized problem obtained during Newton iterations. Finally we present a number of numerical results in space dimensions two and three using nonconforming adaptive meshes and show experimental orders of convergence for upwind and centered discretizations of the convection term.In Chapter 2 we propose a finite volume method on general meshes for the numerical simulation of an incompressible and immiscible two-phase flow in porous media. We consider the case that it can be written as a coupled system involving a degenerate parabolic convection-diffusion equation for the saturation together with a uniformly elliptic equation for the global pressure. The numerical scheme, which is implicit in time, allows computations in the case of a heterogeneous and anisotropic permeability tensor. The convective fluxes, which are non monotone with respect to the unknown saturation and discontinuous with respect to the space variables, are discretized by means of a special Godunov scheme. We prove the existence of a discrete solution which converges, along a subsequence, to a solution of the continuous problem. We present a number of numerical results in space dimension two, which confirm the efficiency of the numerical method.Chapter 3 is devoted to the study of a two-phase flow problem in the case that the capillary pressure curve is discontinuous with respect to the space variable. More precisely we assume that the porous medium is composed of two different rocks, so that the capillary pressure is discontinuous across the interface between the rocks. As a consequence the oil saturation and the global pressure are discontinuous across the interface with nonlinear transmission conditions. We discretize the problem by means of a numerical scheme which reduces to a standard finite volume scheme in each sub-domain and prove the convergence of a sequence of approximate solutions towards a weak solution of the continuous problem. The numerical tests show that the scheme can reproduce the oil trapping phenomenon.

La gestion et la dépollution de sites contaminés peuvent être complexe et demandent un investissement important pour localiser les sources de contaminations, zones émettant les flux de polluants les plus importants. Les travaux réalisés proposent une stratégie pour localiser les sources de pollution à partir de mesures in situ de flux massiques et de modélisation inverse. Ainsi, dans le cadre de l’étude, un outil innovant a d’abord été développé afin de mesurer la vitesse des eaux souterraines dans un puits. L’outil appelé DVT (Direct Velocity Tool) a permis de répondre aux contraintes imposées par les outils existants et de mesurer des vitesses d’écoulement très lentes. Des essais en laboratoire et des tests en site réels ont été réalisés et comparés à d’autres outils de mesure. Le DVT permet aussi indirectement de définir la portion de source conduisant au flux de polluant maximum, en le combinant avec une mesure locale de concentration. L’étude présente ensuite l’utilisation de la modélisation inverse pour localiser une source de contaminant et d’estimer les paramètres définissant les caractéristiques du domaine. Pour cela, l'étude s'est faite sur deux cas synthétiques. Pour adapter les méthodes à une véritable gestion de sites pollués, une stratégie itérative est développée en imposant un ajout limité de nouvelles observations à chaque phase de modélisation, basée sur l’approche de type Data Worth. Les résultats de la position de la source sur les deux cas synthétiques ont permis d’évaluer la méthode mise en place et de juger son applicabilité à une problématique réelle. Cette stratégie de localisation de source est par la suite testée sur un site réel à partir (i) de mesures in situ de flux massiques avec les vitesses au DVT et les concentrations et (ii) la modélisation. Les essais ont permis de cibler les forages à mettre en place sur site aidant à localiser la source. Néanmoins, en analysant plus précisément les résultats, le champ de conductivité hydraulique estimé par l'optimisation ne correspond pas à la réalité. De plus, les flux massiques de contaminants ainsi que le ratio des polluants du site, mettent en valeur deux panaches distincts. Une phase finale de modélisation a donc été lancée afin d'estimer (i) la présence potentielle de deux sources et (ii) la chimie de la zone étudiée. Les résultats de la stratégie sont comparés aux mesures geoprobe qui a pu confirmer la présence d’une des deux sources identifiées. / Contaminated sites management and remediation can be complex and require a significant investment to locate the contaminant source, which delivers the higher pollutant mass fluxes. The study proposes a strategy for contaminant source localisation using in situ measurement and inverse modelling. First, an innovative tool was developed to measure groundwater velocity in a well. The developed tool called DVT (Direct Velocity Tool) made it possible to measure a low Darcy flux. Laboratory and field tests were performed with the DVT and compared to other velocity measurement tools. By combining the DVT with a local concentration measurement, it is possible to calculate the mass fluxes passing through wells. Then the thesis present the inverse modeling used for source localisation and parameters estimation. The study was done on two synthetics cases using the non-linear optimisation method. To adapt the method to a real management of polluted sites, an iterative strategy is developed by imposing a limited addition of new observations to each modeling phase. This strategy is base on the Data Worth approach. Source localisation results on the two synthetic cases made it possible to judge the method applicability to a real site problem. The source localisation strategy is then applied to a real site with (i) mass flux measurement with velocities (DVT) and concentrations and (ii) inverse modeling. The modeling phases made it possible to locate the new wells and helped the source localisation. Nevertheless, by analysing the results more precisely, the hydraulic conductivity field estimated by the optimisation did not correspond to reality. In addition, contaminant mass fluxes highlightes two distinct zones of flux. By analysing the pollutant ratio of the site, it appears that two plumes are potentially present. Thus, another inverse modeling phase has been tested (i) to locate the two potential sources and (ii) to estimate the chemistry of the site. Results of the strategy were compared to the geoprobe campaign which confirmed the second source location.

The disposal of waste from agricultural activities has been recognised as a source of environmental contamination by endocrine disrupting chemicals (EDCs). The New Zealand dairy industry produces a large volume of dairy farm effluent, which contains EDCs in the form of estrogens. Most of this dairy farm effluent is applied onto the land for disposal. Groundwater and soil contamination by estrogens following waste application on the land have been reported overseas, but our understanding of the processes and factors governing the fate of estrogens in the soil is poor. Therefore the main goal of the present study was to better understand the fate and transport of estrogens, in particular 17β-estradiol (E2) and estrone (E1) in soil. In order to quantify E1 and E2 in drainage water and soil samples, chemical analysis by gas-chromatography mass-spectrometry (GC-MS) was carried out. This included sample extraction, sample clean-up through silica gel and gel permeation chromatography, and sample extract derivatisation prior to analysis. In order to develop a reliable method to extract estrogens from soil, research was conducted to optimise E1 and E2 extraction conditions by adjusting the number of sonication and shaking events, as well as the volume and type of solvent. Among five solvents and solvent mixtures tested, the best recovery on spiked and aged soil was obtained using an isopropanol/water (1:1) mix. A microcosm experiment was carried out to determine the dissipation rates of E2 and E1, at 8°C and at field capacity, in the Templeton soil sampled at two different depths (5-10 cm and 30-35 cm). The dissipation rates decreased with time and half-life values of 0.6-0.8 d for E1 and 0.3-0.4 d for E2 were found for the two depths studied. A field transport experiment was also carried out in winter, over three months, by applying dairy farm effluent spiked with estrogens onto undisturbed Templeton soil lysimeters (50 cm in diameter and 70 cm deep). The hormones were applied in dairy farm effluent at 120 mg m⁻² for E2 and 137 mg m⁻² for E1. The results of the transport experiment showed that in the presence of preferential/macropore flow pathways 0.3-0.7% of E2 and 8-13% of E1 was recovered in the leachate at the bottom of the lysimeters after 3 months, and 1-7% of the recovered E2 and 3-54% of the recovered E1 was leached within 2 days of application. These results suggest that leaching of estrogens via preferential/macropore flow pathways is the greatest concern for groundwater contamination. In the absence of preferential/macropore flow pathways, a significant amount (> 99.94%) of both hormones dissipated in the top 70 cm of soil, due to sorption and rapid biodegradation. Surprisingly, in all cases, estrogen breakthrough occurred before that of an inert tracer (bromide). This could not be explained by the advection-dispersion transport of estrogens, nor by their presence as antecedent concentrations in the soil. It was therefore suggested that colloidal enhanced transport of estrogens was responsible for the earlier breakthrough of estrogens and caused the leaching of a fraction of the applied estrogens to a soil depth of 70 cm. A two-phase model, adapted from a state-space mixing cell model, was built to describe the observed estrogen transport processes under transient flow. The model takes into account 3 transport processes namely, advection-dispersion, preferential/macropore flow and colloidal enhanced transport. This model was able to successfully describe the estrogen transport observed from the lysimeters.

Le transfert et les effets de 3 éléments traces métalliques (ETMs : Cd, Pb et Zn) dans les réseaux trophiques terrestres ont été étudiés sur une zone de 40 km2 polluée par l'ancienne fonderie de Metaleurop Nord (France). Les objectifs étaient (i) d'identifier les principaux facteurs expliquant les concentrations en ETMs dans les sols et les animaux étudiés, de déterminer à quelles échelles spatiales ces facteurs interviennent et si le paysage module les contaminations; (ii) de mesurer les réponses des espèces étudiées à différents niveaux d'organisation biologique et de rechercher des liens entre réponses individuelles et/ou populationnelles. Des échantillonnages portant sur les sols d'habitats ligneux (milieux boisés favorables pour les espèces étudiées) et plusieurs espèces d'escargots et de micromammifères ont été menés durant 3 années consécutives à plusieurs saisons. L'analyse de la distribution spatiale de la contamination des sols montre que les variables expliquant les concentrations totales en ETMs dans les sols (distance à la source, fréquence du vent) agissent à des échelles spatiales relativement larges, alors que des facteurs locaux (caractéristiques des sols, notamment le pH) déterminent les concentrations extractibles au CaCl2 (considérées comme disponibles). Des méthodes géostatistiques d'interpolation spatiale, krigeage ordinaire et universel, ont permis de générer des cartes de contamination des sols sur l'ensemble de la zone d'étude et d'établir une évaluation du risque spatialisée pour les invertébrés et vertébrés des habitats ligneux (basée sur les Eco-SSLs de l'US EPA). Les relations entre les concentrations en ETMs des animaux et des sols ont été étudiées chez 2 espèces d'escargots (Cepaea sp et Oxychilus draparnaudi) et 2 espèces de micromammifères (Crocidure musette Crocidura russula et Campagnol roussâtre Myodes glareolus) sympatriques ayant des régimes alimentaires contrastés. Les concentrations internes augmentent avec la contamination des sols et sont toujours mieux expliquées par les concentrations totales qu'extractibles, sauf pour Cepaea sp. Le paysage influence les niveaux d'ETMs accumulées et leur augmentation le long du gradient de pollution. L'échelle spatiale à laquelle la corrélation entre les concentrations en ETMs dans les animaux et les sols est la meilleure varie entre espèces et entre paysages. Les facteurs physiologiques et écologiques, tels que le régime alimentaire mais aussi la préférence et l'exploitation d'habitat, modulent la bioaccumulation des ETMs. L'influence des caractéristiques du récepteur (sexe, âge, régime alimentaire...) sur l'accumulation et la séquestration des ETMs a été étudiée chez plusieurs micromammifères sympatriques (Mulot sylvestre Apodemus sylvaticus, Campagnol roussâtre, Musaraigne carrelet Sorex araneus et Musaraigne pygmée Sorex minutus). L'âge plus que le sexe, et l'espèce plus que le niveau trophique, influencent les niveaux d'ETMs et de métallothionéines (MTs) et leurs relations. Trois patrons de relations entre capacité d'accumulation des ETMs et production de MTs ont été identifiés. Des mesures de santé individuelle suggèrent que la sensibilité des récepteurs aux ETMs (musaraigne < campagnol < mulot) décroît quand leur capacité de production de MTs augmente. La mesure des MTs seule ne constitue pas un biomarqueur d'exposition pertinent, mais elle est nécessaire pour interpréter la probabilité d'occurrence d'effets délétères. La variabilité temporelle de l'accumulation des ETMs a été étudiée chez le Campagnol roussâtre, mettant en évidence des différences inter-annuelles et inter-saisonnières. Tandis que le Pb (foie et reins) montre des variations inter-saisonnières (printemps < automne), le Cd et le Zn ont des patrons différents entre organes, et des variations inter-annuelles sans effet saisonnier net. L'adaptation potentielle à la pollution métallique de 2 escargots autochtones du site d'étude, l'escargot petit-gris Cantareus aspersus et l'escargot des bois Cepaea nemoralis, a été examinée : des individus issus de 2 origines (parents exposés ou non aux ETMs) ont été exposés in situ en microcosmes dans une expérience de transplantation croisée. L'analyse de l'accumulation des ETMs (facteur de bioaccumulation, paramètres cinétiques) ne montrent pas d'adaptation. Cependant, les caractéristiques de la coquille suggèrent des processus d'acclimatation et d'adaptation à travers une modification du métabolisme du Ca. Ce travail montre que l'analyse de la réponse d'un organisme à la pollution doit s'appuyer sur la mesure de plusieurs paramètres (bioaccumulation, biomarqueur, indice de santé, population...), et que l'interprétation d'une mesure nécessite sa mise en perspective dans son contexte global en changeant de niveaux de perception (niveau d'organisation biologique, contextes temporel et spatial). Dans le cadre de l'évaluation du risque environnemental et de la biosurveillance, nos travaux soulignent la nécessité des suivis à long terme, et de la mise en place de procédures sur des étendues pertinentes par rapport aux processus écologiques qui participent aux réponses observées et tenant compte du paysage.

Two of the most important environmental challenges in the 21st century are to protect the quality of fresh water resources and to utilize renewable energy sources to lower greenhouse gas emissions. This study contributes to the solution of the first challenge by providing methodologies for optimal design of real-time water quality monitoring systems and interpretation of data supplied by the monitoring system to identify potential pollution sources in river networks. In this study, the optimal river water quality monitoring network design aspect of the overall monitoring program is addressed by a novel methodology for the analysis of this problem. In this analysis, the locations of sampling sites are determined such that the contaminant detection time is minimized for the river network while achieving maximum reliability for the monitoring system performance. The data collected from these monitoring stations can be used to identify contamination source locations. This study suggests a methodology that utilizes a classification routine which associates the observations on a contaminant spill with one or more of the candidate spill locations in the river network. This approach consists of a training step followed by a sequential elimination of the candidate spill locations which lead to the identification of potential spill locations. In order to contribute the solution of the second environmental challenge, this study suggests utilizing available excess energy in water distribution systems by providing a methodology for optimal design of energy recovery systems. The energy recovery in water distribution systems is possible by using micro hydroelectric turbines to harvest available excess energy inevitably produced to satisfy consumer demands and to maintain adequate pressures. In this study, an optimization approach for the design of energy recovery systems in water distribution networks is proposed. This methodology is based on finding the best locations for micro hydroelectric plants in the network to recover the excess energy. Due to the unsteady nature of flow in water distribution networks, the proposed methodology also determines optimum operation schedules for the micro turbines.

The research described herein is composed of four major areas of practice. It examines the overall performance of runways and provides tools designed to improve current runway operations and management with particular emphasis on contaminated surfaces. Presented in this thesis is an overview of how to design airport pavements in order to achieve optimal friction by specifically focusing on material selection and construction techniques for rigid and flexible pavements. Rubber buildup and the impact rubber accumulation has on decreasing runway friction, particularly in a range of climatic conditions, is discussed. Four commonly used rubber removal techniques are presented and evaluated. Through this research, an analytical hierarchy process (AHP) decision making protocol was developed for incorporation into airport pavement management systems (APMS). Runway surface condition reporting practices used at the Region of Waterloo International Airport are evaluated and recommendations for improving current practices are identified. Runway surface condition reporting can be improved by removing subjectivity, reporting conditions to pilots in real time, standardizing terminology and measurement techniques, and including runway pictures or sketches to identify contaminant locations where possible. Reports should be incorporated and stored in the APMS. Aircraft braking systems and their effects on landing distances under contaminated conditions are discussed. This thesis presents a proposed solution for monitoring and measuring contaminated runway surfaces and identifying the risks associated with aircraft landing through using the Braking Availability Tester (BAT). Also proposed in this thesis is a testing framework for validating the Braking Availability Tester. The proposed BAT measures interaction between aircraft antiskid braking systems and runway contaminants to determine landing distances more accurately. Finally, this thesis includes a discussion explaining how pavement design, contaminant removal, results from friction tests, and results from the BAT can be incorporated into airport pavement management systems. APMS data can be analyzed to economically optimize and prioritize scheduling of pavement maintenance, preservation and rehabilitation treatments to maintain a high level of service, thereby contributing to runway safety and optimization.

No description available.

The research described herein is composed of four major areas of practice. It examines the overall performance of runways and provides tools designed to improve current runway operations and management with particular emphasis on contaminated surfaces. Presented in this thesis is an overview of how to design airport pavements in order to achieve optimal friction by specifically focusing on material selection and construction techniques for rigid and flexible pavements. Rubber buildup and the impact rubber accumulation has on decreasing runway friction, particularly in a range of climatic conditions, is discussed. Four commonly used rubber removal techniques are presented and evaluated. Through this research, an analytical hierarchy process (AHP) decision making protocol was developed for incorporation into airport pavement management systems (APMS). Runway surface condition reporting practices used at the Region of Waterloo International Airport are evaluated and recommendations for improving current practices are identified. Runway surface condition reporting can be improved by removing subjectivity, reporting conditions to pilots in real time, standardizing terminology and measurement techniques, and including runway pictures or sketches to identify contaminant locations where possible. Reports should be incorporated and stored in the APMS. Aircraft braking systems and their effects on landing distances under contaminated conditions are discussed. This thesis presents a proposed solution for monitoring and measuring contaminated runway surfaces and identifying the risks associated with aircraft landing through using the Braking Availability Tester (BAT). Also proposed in this thesis is a testing framework for validating the Braking Availability Tester. The proposed BAT measures interaction between aircraft antiskid braking systems and runway contaminants to determine landing distances more accurately. Finally, this thesis includes a discussion explaining how pavement design, contaminant removal, results from friction tests, and results from the BAT can be incorporated into airport pavement management systems. APMS data can be analyzed to economically optimize and prioritize scheduling of pavement maintenance, preservation and rehabilitation treatments to maintain a high level of service, thereby contributing to runway safety and optimization.

In this report an equation is derived to calculate the dynamic effect of primary and secondary aircraft cabin air contamination. The equation is applied in order to understand implications and hazards. Primary contamination is from an outside source in form of normal low level contamination or high level contamination in a failure case. Secondary contamination originates from deposited material released into the cabin by a trigger event. The dynamic effect is described as an initial value problem (IVP) of a system governed by a nonhomogeneous linear first order ordinary differential equation (ODE). More complicated excitations are treated as a sequence of IVPs. The ODE is solved from first principles. Spreadsheets are provided with sample calculations that can be adapted to user needs. The method is not limited to a particular principle of the environmental control system (ECS) or contamination substance. The report considers cabin air recirculation and several locations of contamination sources, filters, and deposit points (where contaminants can accumulate and from where they can be released). This is a level of detail so far not considered in the cabin air literature. Various primary and secondary cabin contamination scenarios are calculated with plausible input parameters taken from popular passenger aircraft. A large cabin volume, high air exchange rate, large filtered air recirculation rate, and high absorption rates at deposit points lead to low contamination concentration at given source strength. Especially high contamination concentrations would result if large deposits of contaminants are released in a short time. The accuracy of the results depends on the accuracy of the input parameters. Five different approaches to reduce the contaminant concentration in the aircraft cabin are discussed and evaluated. More effective solutions involve higher implementation efforts. The method and the spreadsheets allow predicting cabin air contamination concentrations independent of confidential industrial input parameters.