Many types of contaminants of emerging concern (CECs), including per- and poly-fluoroalkyl substances (PFAS), have been found in leachate of operating municipal landfills. However, information on CECs in leachate of historic landfills (≥3 decades since closure, often lacking engineered liners or leachate collection systems) and the related risk posed from groundwater plumes discharging to nearby aquatic ecosystems is limited. In this study, 48 samples of leachate-impacted groundwater were collected from 20 historic landfills in Ontario, Canada. The CECs measured included artificial sweeteners (ASs), PFAS, organophosphate esters (OPE), pharmaceuticals, bisphenols, sulfamic acid, perchlorate, and substituted phenols. Several landfills, including ones closed in the 1960s, had total PFAS concentrations similar to those previously measured at modern landfills, with a maximum observed here of 12.7 μg/L. Notably elevated concentrations of several OPE, cotinine, and bisphenols A and S were found at many 30-60 year-old landfills. There was little indication of declining concentrations with landfill age, suggesting historic landfills can be long-term sources of CECs to groundwater. A full-year field study was performed on a 0.5-km reach of an urban stream receiving contaminated groundwater from nearby historic landfills. Elevated concentrations of ammonium, the AS saccharin, an indicator of old landfill leachate, and CECs (e.g., maximum total PFAS of 31 μg/L) in the shallow discharging groundwater were relatively stable across the seasons but were spatially restricted by hyporheic exchange and discharge of other groundwater. This indicates a patchy but long-term exposure for endobenthic organisms, which are rarely monitored. Stream water concentrations were more dilute, but increased markedly across the landfill stretch, and showed signs of increases in winter and after rain/snowmelt events. These findings provide guidance on which CECs may require monitoring at historic landfill sites and suggest how landfill monitoring programs could be improved to fully capture the risk to receiving water bodies. / Thesis / Master of Science (MSc) / Historic landfills are a known source of groundwater contamination. This study investigated whether these landfills contain new groups of chemicals, called contaminants of emerging concern (CECs), which are suspected to pose serious environmental and human health risks. This study found many CECs at high concentrations in most of the 20 historic landfill sites investigated, even those closed up to 60 years. A full-year investigation at one historic landfill site showed that organisms living in the sediments of a nearby stream are exposed to high concentrations all year long. Concentrations in the stream increased as it flowed past the landfill, and may be higher in winter and after rains, times monitoring is rarely done. The elevated concentrations of harmful contaminants in this water are potentially threatening the stream ecosystem. Operators of historic sites should consider testing for CECs and ensure that monitoring strategies accurately evaluate the risk posed to the environment.

Decision makers' conflicts about the validity of a single simulation model and inefficiencies of existing response matrix methods (RMM) hinder adopting successful groundwater management plans. We speed up the process by proposing a hybrid RMM that is most efficient for situations in which optimizable stimuli can vary through consecutive periods of uniform duration interspersed with periods of different duration. We use the hybrid RMM within Simulation-Optimization (S-O) models to develop optimal water management strategies. For the tested problems, the hybrid RMM requires as much or 63-89% less computation time than other RMMs. Second, we propose Multi-Conceptual Model Optimization (MCMO) that can help stakeholders reach a compromise strategy instead of agreeing on the validity of a single model. MCMO computes optimal strategies that simultaneously satisfy analogous constraints and bounds in multiple numerical models differing by more than parameter values. Applying MCMO to Cache Valley (Utah, USA) reveals that protecting local ecosystem limits the increased groundwater pumping to satisfy only 40% of projected water demand increase using two models. To successfully and sustainably manage Cache Valley aquifer, we evaluate sustained yield strategies (SYS) and quantify the resilience of a computed SYS. We maximize the number of new residents who can have their indoor and outdoor uses satisfied, subject to constraints on aquifer-surface waters conditions, and limiting new residents to projected increases in population (PIiP). furthermore, we examine the effect of optimization approach and sequiencing, temporally-lagged spatially distributed return flow that is a function of optimal groundwater use, and the acceptability time evaluation on the optimal yield strategy. Cache Valley aquifer can sustainably satisfy the outdoor water demand of 74%-83% and the indoor water demand of 83%-100% of the PIiP. We quantify deterministic resilience Rd(A,T,SV)=P to evaluate how completely an aquifer condition (SV) recovers after the end of climatic anomaly (A), by recovery time (T). Simulation predicts that Cache Valley aquifer system resiliences to a 2-year drought are Rd(2YD, 3 yrs, Overall) = 93% and Rd (2YD,≥8,Overall) ≥ 95%. Proportionally reducing pumping rates by 25% through the time horizon of the simulation increases the overall resilience to 96% within 3 years.

Groundwater is an important source of baseflow. Baseflow supports minimum flows and living area through dry periods, and moderates surface water temperature. The reductions in baseflow after urbanization can cause degradation of the stream ecosystem. The purpose of this study was to investigate the groundwater-surface water interaction and thermal regime of Clythe Creek, Guelph to illustrate the importance of groundwater/streamflow interaction in determining the health of a stream ecosystem. Piezometers were used to quantify vertical and lateral groundwater flow direction, and surface water temperature measurements were used to characterize and quantify the thermal regime. Groundwater-surface water interaction varied both temporally and spatially between the two geomorphic units. Average summer surface water temperatures were consistently cooler in the downstream portion of the study reach. The importance of groundwater-surface water interactions in supporting stream ecosystems was concluded from this study, and restoration strategies to address threats from urbanization were made. / Ontario Research Fund (ORF)

No description available.

Abstract Water resources management calls for methods to simultaneously manage groundwater (GW) and surface water (SW) systems. These have traditionally been considered separate units of the hydrological cycle, which has led to oversimplification of exchange processes at the GW-SW interface. This thesis studied GW hydrology and the previously unrecognised connection of the Rokua esker aquifer with lakes and streams in the area, with the aim of identifying reasons for lake water level variability and eutrophication in the Rokua esker. GW-SW interactions in the aquifer were first studied with field methods. Seepage meter measurements showed substantial spatial variability in GW-lake interaction, whereas transient variability was more modest, although present and related to the surrounding aquifer. Environmental tracers suggested that water exchange occurs in all lakes in the area, but is of varying magnitude in different lakes. Finally, GW-SW interaction was studied in peatland catchments, where drainage channels in the peat soil presumably increased groundwater outflow from the aquifer. Amount and rate of GW recharge were then estimated with a simulation approach developed explicitly to account for the physical characteristics of the Rokua esker aquifer. This produced a spatially and temporally distributed recharge estimate, which was validated by independent field techniques. The results highlighted the impact of canopy characteristics, and thereby forestry management, on GW recharge. The data collected and the new understanding of site hydrology obtained were refined into a fully integrated surface-subsurface flow model of the Rokua aquifer. Simulation results compared favourably to field observations of GW, lake levels and stream discharge. A major finding was of good agreement between simulated and observed GW inflow to lakes in terms of discharge locations and total influx. This thesis demonstrates the importance of using multiple methods to gain a comprehensive understanding of esker aquifer hydrology with interconnected lakes and streams. Importantly, site-specific information on the reasons for water table variability and the trophic status of Rokua lakes, which is causing local concern, is provided. As the main outcome, various field and modelling methods were tested, refined and shown to be suitable for integrated GW and SW resource management in esker aquifers. / Tiivistelmä Vesivarojen hallinnassa tarvitaan menetelmiä pohja- ja pintaveden kokonaisvaltaiseen huomioimiseen. Pohja- ja pintavesiä tarkastellaan usein erillisinä osina hydrologista kiertoa, mikä on johtanut niiden välisten virtausprosessien yksinkertaistamiseen. Tässä työssä selvitettiin Rokuan pohjavesiesiintymän hydrologiaa ja hydraulista yhteyttä alueella oleviin järviin ja puroihin. Tutkimuksessa pyrittiin osaltaan selvittämään syitä harjualueen järvien pinnanvaihteluun ja veden laatuongelmiin. Kenttätutkimuksissa todettiin voimakasta alueellista vaihtelua järven ja pohjaveden vuorovaikutuksessa. Pohjaveden suotautumisen ajallinen vaihtelu puolestaan oli vähäisempää, mutta havaittavissa, ja kytköksissä järveä ympäröivän pohjavesipinnan vaihteluihin. Merkkiaineet vesinäytteistä viittasivat vastaavan vuorovaikutuksen olevan läsnä myös muissa alueen järvissä, mutta suotautuvan pohjaveden määrän vaihtelevan järvittäin. Turvemailla tehdyt mittaukset osoittivat pohjaveden purkautuvan ojaverkostoon ja ojituksen mahdollisesti lisäävän ulosvirtaamaa pohjavesiesiintymästä. Pohjaveden muodostumismäärää ja -nopeutta tutkittiin numeerisella mallinnuksella, joka kehitettiin huomioimaan harjualueelle ominaiset fysikaaliset tekijät. Mallinnus tuotti arvion ajallisesti ja alueellisesti vaihtelevasta pohjaveden muodostumisesta, joka varmennettiin kenttämittauksilla. Tuloksissa korostui kasvillisuuden, ja sitä kautta metsähakkuiden, vaikutus pohjaveden muodostumismääriin. Hydrologiasta kerätyn aineiston ja kehittyneen prosessiymmärryksen avulla Rokuan harjualueesta muodostettiin täysin integroitu numeerinen pohjavesi-pintavesi virtausmalli. Mallinnustulokset vastasivat mittauksia pohjaveden ja järvien pinnantasoista sekä purovirtaamista. Työn merkittävin tulos oli, että mallinnetut pohjaveden purkautumiskohdat ja purkautumismäärät alueen järviin vastasivat kenttähavaintoja. Tämä työ havainnollisti, että ymmärtääkseen pohjaveden ja siitä riippuvaisten järvien ja purojen vuorovaikutusta harjualueella on käytettävä monipuolisia tutkimusmenetelmiä. Työ toi lisätietoa Rokuan harjualueen vesiongelmien syihin selittäen järvien vedenpinnan vaihtelua ja vedenlaatua pohjavesihydrologialla. Väitöstyön tärkein anti oli erilaisten kenttä- ja mallinnus-menetelmien soveltaminen, kehittäminen ja hyödylliseksi havaitseminen harjualueiden kokonaisvaltaisessa pinta- ja pohjavesien hallinnassa.

Periodic releases from an upstream dam cause rapid stage fluctuations in the Colorado River near Austin, Texas. These daily pulses modulate fluid exchange and residence times in the hyporheic region, where biogeochemical reactions are pronounced. We installed two transects of wells perpendicular to the river to examine in detail the reactions occurring in this zone of surface-water and groundwater exchange. One well transect recorded physical water level fluctuations and allowed us to map hydraulic head gradients and fluid movement. The second transect allowed for water sample collection at three discrete depths. Samples were collected from 12 wells every 2 hours for a 24-hour period and were analyzed for nutrients, carbon, major ions, and stable isotopes. The results provide a detailed picture of biogeochemical processes in the bank environment during low flow/drought conditions in a regulated river. Findings indicate that a pulse that causes a change in river stage of approximately 16-centimeters does not cause significant mixing in the bank. Under these conditions, the two systems act independently and exhibit only slight mixing at the interface. / text

Master of Science / Department of Civil Engineering / David R. Steward / Accurate modeling of the High Plains Aquifer depends on the availability of good data that represents and quantities properties and processes occurring within the aquifer. Thanks to many previous studies there is a wealth of good data available for the High Plains Aquifer but one key component, groundwater-surface water interaction locations and rates, is generally missing. Without these values accurate modeling of the High Plains Aquifer is very difficult to achieve. This thesis presents methods for simplifying the modeling of the High Plains Aquifer using a sloping base method and then applying mathematical optimization techniques to locate and quantify points of groundwater-surface water interaction. The High Plains Aquifer has a base that slopes gently from west to east and is approximated using a one-dimensional stepping base model. The model was run under steady-state predevelopment conditions using readily available GIS data representing aquifer properties such as hydraulic conductivity, bedrock elevation, recharge, and the predevelopment water level. The Levenberg-Marquardt and particle swarm algorithms were implemented to minimize error in the model. The algorithms reduced model error by finding locations in the aquifer of potential groundwater-surface water interaction and then determining the rate of groundwater to surface water exchange at those points that allowed for the best match between the measured predevelopment water level and the simulated water level. Results from the model indicate that groundwater-surface water interaction plays an important role in the overall water balance in the High Plains Aquifer. Findings from the model show strong groundwater-surface water interaction occurring in the northern basin of the aquifer where the water table is relatively shallow and there are many surface water features. In the central and southern basins the interaction is primarily limited to river valleys. Most rivers have baseflow that is a net sink from groundwater.

The San Joaquin Valley is part of the Great Central Valley of California, a major agricultural centre and food supplier for the United States. This area has significant water management concerns given the very high water demand for an increasing state population and for intense irrigation in a hot, temperate to semi-arid climate where the overall rate of evapotranspiration (ET) is high, and the overall rate of precipitation is low. Irrigation heavily relies upon groundwater and surface water extractions. Through the historical and current concerns of regional water resources reliability, land surface subsidence, water quality issues, and the health of ecosystems, a need for regional-scale water resource management and planning has developed. The physically-based surface-subsurface HydroGeoSphere (HGS) model is used to examine the regional-scale hydrologic budget of a large portion of the San Joaquin Valley. The objective of this investigation is to develop a steady-state groundwater-surface water model of the San Joaquin Valley representative of predevelopment hydrologic conditions. The groundwater-surface water system has undergone drastic changes since the employment of groundwater and surface water extractions for irrigation and mining, and is still responding to past and present stresses. The only certain stable initial condition must therefore be that of the natural system. The model input parameters were constrained by all relevant available hydrologic data. The model was not calibrated to subsurface hydraulic heads or river flows. However, the model does provide a fair match between simulated and actual estimated water table elevations. Historic river flow estimates were not used to calibrate the model, because data consistent with that collected by Hall (1886) and representative of the natural system were not available. For this investigation, water enters through precipitation and the inflow of major rivers only. The subsurface domain is bounded by no-flow boundaries, and groundwater is therefore only able to exit the subsurface through discharge to surface water features or through ET. Surface water is only able to exit the model through discharge via the San Joaquin River and through ET. Average river inflows circa 1878 to 1884 documented by Hall (1886) were applied where the rivers enter into the valley. The spatially variable average rate of precipitation (years 1971 to 2000) from a PRISM dataset was applied to the top of the model. The spatially variable long term average potential ET rates from the California Department of Water Resources (DWR) et al. (1999) were applied to the top of the model. Averaged overland flow parameters and vegetation factors needed to calculate actual ET were specified at the top of the model based on literature values and the 1874 spatial distribution of natural vegetation provided by California State University at Chico et al. (2003). Hydrogeological data including hydraulic conductivities, porosities, specific storage, and unsaturated zone properties are based on literature values from other relevant studies. The resulting steady state model is therefore characterized by historical long term average data assumed to be representative (as close as possible) of the flow system circa 1848. Results indicate that the natural hydrologic setting of the San Joaquin Valley is a complex one. Complex hydrologic processes, including significant groundwater-surface water interaction along the major rivers and within wetland areas formed by flooded surface water, as well as ET and impacted root zone processes were identified in the model domain. Identification and simulation of the complex recharge and discharge relationships in the model domain sheds insight into the hydrologic nature of some historic natural wetlands. Evapotranspiration is a very significant sink of both surface water and groundwater (44.8 % of the water balance input), and has a major impact on hydrologic processes in the root zone. The presence and path of the major rivers in the domain are well defined in the model output and agree well with their actual locations. The model simulates gaining and losing reaches of the major rivers, replicating the historic recharge-discharge relationship documented by others. The general location, formation, and hydrologic processes of some significant wetlands simulated by the model have a fair agreement with historical records. As mentioned above, there is also a fair match between simulated and actual estimated water table elevations. Successful simulation of the complex hydrologic processes and features that characterize the predevelopment hydrologic conditions of the San Joaquin Valley and that resolve the water balance of the natural system underscores the importance and necessity of using an integrated model. This steady state model should serve as a reasonable initial condition for future transient runs that bring the model up to current hydrologic conditions capable of estimating present and future water budgets.

The San Joaquin Valley is part of the Great Central Valley of California, a major agricultural centre and food supplier for the United States. This area has significant water management concerns given the very high water demand for an increasing state population and for intense irrigation in a hot, temperate to semi-arid climate where the overall rate of evapotranspiration (ET) is high, and the overall rate of precipitation is low. Irrigation heavily relies upon groundwater and surface water extractions. Through the historical and current concerns of regional water resources reliability, land surface subsidence, water quality issues, and the health of ecosystems, a need for regional-scale water resource management and planning has developed. The physically-based surface-subsurface HydroGeoSphere (HGS) model is used to examine the regional-scale hydrologic budget of a large portion of the San Joaquin Valley. The objective of this investigation is to develop a steady-state groundwater-surface water model of the San Joaquin Valley representative of predevelopment hydrologic conditions. The groundwater-surface water system has undergone drastic changes since the employment of groundwater and surface water extractions for irrigation and mining, and is still responding to past and present stresses. The only certain stable initial condition must therefore be that of the natural system. The model input parameters were constrained by all relevant available hydrologic data. The model was not calibrated to subsurface hydraulic heads or river flows. However, the model does provide a fair match between simulated and actual estimated water table elevations. Historic river flow estimates were not used to calibrate the model, because data consistent with that collected by Hall (1886) and representative of the natural system were not available. For this investigation, water enters through precipitation and the inflow of major rivers only. The subsurface domain is bounded by no-flow boundaries, and groundwater is therefore only able to exit the subsurface through discharge to surface water features or through ET. Surface water is only able to exit the model through discharge via the San Joaquin River and through ET. Average river inflows circa 1878 to 1884 documented by Hall (1886) were applied where the rivers enter into the valley. The spatially variable average rate of precipitation (years 1971 to 2000) from a PRISM dataset was applied to the top of the model. The spatially variable long term average potential ET rates from the California Department of Water Resources (DWR) et al. (1999) were applied to the top of the model. Averaged overland flow parameters and vegetation factors needed to calculate actual ET were specified at the top of the model based on literature values and the 1874 spatial distribution of natural vegetation provided by California State University at Chico et al. (2003). Hydrogeological data including hydraulic conductivities, porosities, specific storage, and unsaturated zone properties are based on literature values from other relevant studies. The resulting steady state model is therefore characterized by historical long term average data assumed to be representative (as close as possible) of the flow system circa 1848. Results indicate that the natural hydrologic setting of the San Joaquin Valley is a complex one. Complex hydrologic processes, including significant groundwater-surface water interaction along the major rivers and within wetland areas formed by flooded surface water, as well as ET and impacted root zone processes were identified in the model domain. Identification and simulation of the complex recharge and discharge relationships in the model domain sheds insight into the hydrologic nature of some historic natural wetlands. Evapotranspiration is a very significant sink of both surface water and groundwater (44.8 % of the water balance input), and has a major impact on hydrologic processes in the root zone. The presence and path of the major rivers in the domain are well defined in the model output and agree well with their actual locations. The model simulates gaining and losing reaches of the major rivers, replicating the historic recharge-discharge relationship documented by others. The general location, formation, and hydrologic processes of some significant wetlands simulated by the model have a fair agreement with historical records. As mentioned above, there is also a fair match between simulated and actual estimated water table elevations. Successful simulation of the complex hydrologic processes and features that characterize the predevelopment hydrologic conditions of the San Joaquin Valley and that resolve the water balance of the natural system underscores the importance and necessity of using an integrated model. This steady state model should serve as a reasonable initial condition for future transient runs that bring the model up to current hydrologic conditions capable of estimating present and future water budgets.

Natürliche Niederungsgebiete wirken als effektive Phosphor Senke. Landwirtschaftliche Nutzung und Drainage führt zur Transformation zu Phosphor Quellen. Dem Spreewald, ein Feuchtgebiet und seiner Funktionsweise als Phosphor Quelle oder –Senke, kommt im Einzugsgebiet der Spree eine wichtige Rolle zu. Die vorhandenen Kenntnisse der Phosphor Umsatzprozesse und des Phosphor Austausches zwischen den Flächen und dem Fließgewässernetz sind jedoch gering. Praktikable Ansätze zur Beschreibung des Phosphor Austauschverhaltens von Nährstoffen in den ausgedehnten Polderregionen liegen nicht vor. Anhand hydrologischer, geohydraulischer und biogeochemischer Prozessuntersuchungen in einer Polderlandschaft mit typischer Stauhaltung konnten die Phosphor Freisetzungs- und Retentionsprozesse ausgewiesen werden. Unter Nutzung verschiedener Modellansätze (geohydraulische Modellierung, Stofftransportmodellierung und Statistische Modelle) und der Berechnung von Phosphor Prozessraten werden horizontale und vertikale Phosphor Fluxe in den Flächen-Wasser Übergangszonen quantifiziert. Die Ergebnisse gehen in ein Phosphor Bilanzmodell ein, das die P-Quellen und Senkenfunktion von Polder Teilflächen (überstaut, genutzt) und des gesamten Polders in Monatsschritten abbildet. Die biogeochemischen Phosphor Umsatzprozesse sowie der horizontale und vertikale Transport werden wesentlich von den hydrologischen und den klimatischen Bedingungen gesteuert. Sie stellen die primären Einflussgrößen der P Senken oder -Quellenfunktion dar. Im Polder wird die Phosphor Netto Freisetzung der genutzten Bereiche durch eine hohe Phosphor Netto Retention in den überstauten Flächen abgemindert. Szenario Untersuchungen zeigen, dass sinkende Grabenwasserstände zu einer erheblichen Erhöhung der Phosphor Emissionen aus dem Grundwasserpfad führen. Ansteigende Phosphor Fluxe bei sinkenden Wasserständen legen ein zunehmendes Eutrophierungsrisiko in den Gräben sowie für unterhalb gelegene aquatische Systeme nahe. / Natural wetlands effectively retain phosphorus. Agricultural cultivation and drainage by ditches transform them to phosphorus sources. In the Spree catchment, the Spreewald, a large scaled lowland has a strategic importance operating as a phosphorus sink or phosphorus source. Anyhow, knowledge of its phosphorus turnover processes and the phosphorus exchange behavior between the plain and the river and ditch network are marginal. Practicable approaches to reproduce the phosphorus exchange behavior in its typical polder areas are missing. Based on process investigations in a polder area with typical weir regulation, phosphorus retention and remobilization processes are characterized. Combining different model approaches (groundwater modeling, matter transport modeling and statistical models) and calculating process rates, the horizontal and vertical phosphor fluxes in the water soil/sediment transition zones are quantified. The outcomes are used as input data for a phosphorus balance model reproducing the phosphorus source and sink character of used and rewetted polder areas and for the whole polder in monthly time steps. Results from process and transport investigations point out that biogeochemical turnover processes and horizontal or vertical phosphorus transport are driven by hydrological and climatological conditions. Net phosphorus release found in the extensive used polder areas is counteracted by significant net phosphorus retention in the rewetted parts. Scenario analyses identify decreasing ditch water levels to cause a considerable increase of phosphorus emissions from the groundwater pathway. Consequently, the increasing phosphorus fluxes at decreasing water levels provoke a growing risk for eutrophication in the ditches but also in the downstream aquatic systems.