No description available.

One of the key areas of climate change impacts is the water resources. Surface water has been compromised by climate change, and so has groundwater. Sub-Saharan Africa has been one of the major sufferers of climate change since the early 70s. The limitations of current global climate models in the prediction of the future climate over the continent have been a source of research challenges in the last decade. Over the decades, West Africa and the Sahel region have been subjected to major droughts and recurrent dry spells. While most studies in the region have tackled climate change effects on surface water and agriculture, a few have highlighted its effects on groundwater. This study investigates climate change impacts on both surface and groundwater in the Kou basin in Burkina Faso. The Kou River and its tributaries have experienced depletion over the last decades. Therefore, it has become necessary to investigate the deficit of the river flow and how the groundwater that forms the springs that feed the rivers is affected by climate change. The study resorts to integrated hydrologic modelling approach using the SWAT surface model and the MODFLOW groundwater model to assess the impact of climate change. Simulations from selected Regional Climate Models (RCMs) are used. Investigations from the RCMs show that the study area is expected to be drier with less precipitation and higher temperatures. Furthermore, the SWAT model results show that as rainfall reduces, future stream flows also significantly decrease. Results from the MODFLOW model also follow the trend of the SWAT model. Groundwater levels are declining whatever the RCM outputs used. Overall, all the models predicted a severely decreasing trend in surface and groundwater in the Kou basin. The study results will be particularly useful for water resources managers in the Kou River basin. / Thesis / Master of Science (MSc)

Correctly predicting water levels is key to the success of created wetlands. The Pierce method is a commonly used technique for modeling and designing mitigation wetlands that assumes minimal groundwater interaction with the wetland. This technique for mitigation wetland design relies primarily on surface water inputs, assuming a relatively impermeable substrate (perched system), and level pool routing. The Pierce method was applied utilizing two different evapotranspiration estimation methods: Thornthwaite (IPM) and FAO-56 Penman-Monteith (IPM-FAO). A second process-based model, utilizing MODFLOW-2005, was constructed to better predict water levels in mitigation wetlands. Modeled processes included groundwater movement and vegetative resistance to flow, which can be a significant factor in wetland water levels. The two versions of the Pierce method were compared to the process-based wetland representation developed in MODFLOW-2005 using data from an existing mitigation wetland. Output from these models were compared to observed data from an existing mitigation wetland near Manassas, VA, USA. Results indicate the use of Thornthwaite's method to estimate wetland evapotranspiration (ET) does not capture the timing or magnitude of wetland ET losses, leading to over-prediction of wetland water levels during the growing season. The Modflow-based approach resulted in more accurate hydroperiod predictions on a yearly basis than the Pierce Method. However, the Integrated Pierce method model, utilizing the FAO-56 Penman-Monteith method of estimating potential evapotranspiration instead of Thornthwaite's method most accurately predicted water levels during the growing season (March-October). / Master of Science

Ur ett byggtekniskt perspektiv kan en sänkning i grundvattenytan i en sluten akvifer i ett område med lerjordar ge sättningar som kan skada byggkonstruktioner. Kopplingen mellan hydrogeologi och geoteknik är tydlig men oftast görs grova uppskattningar av konsekvenserna av en grundvattensänkning. Detta beror på att sättning vanligen beräknas i enskilda punkter där data finns att tillgå vilket endast ger resultat för sättningsberäkningar i dessa punkter. Ytor emellan punkterna utelämnas ofta. Dessutom är det inte vanligt att grundvattensänkningen beräknas med etablerade mjukvaruverktyg som Modflow för att få en mer detaljerad bild av avsänkningen och påverkansområdet.Denna studie kopplar samman en numerisk grundvattenmodell simulerad med Modflow samt en jordmodell, framtagen och interpolerad med kriging, med sättningsberäkningar. Detta resulterar i en integrerad modell som har till syfte att generera översiktskartor med predikterad sättning som resultat av grundvattensänkningar i utvalt område. Den integrerade modellen och sättningsberäkningarna är programmerade till denna studie med beräkningsverktyget Octave. Den integrerade sättningsmodellen testas på en fallstudie med verkliga geotekniska och hydrogeologiska data från ett område i Mälardalen. I fallstudien har ett hypotetiskt fall av grundvattensänkning simulerats. Grundläggande hydrogeologisk teori för slutna akviferer används för att bedöma vilka laster och ökad effektiv spänning grundvatten-sänkningar ger upphov till i jorden för fallstudien.Resultaten, i form av meter sättning, från den integrerade sättningsmodellen har verifierats mot beräkningar i programvaran Geosuite Settlement som är ett erkänt verktyg för sättningsberäkningar. Det visar sig att den integrerade sättningsmodellen beräknar sättningar med stor noggrannhet. Resultatet från grundvattenmodellen jämförs med en kriginginterpolerad grundvattenyta baserad på mätningar i grund-vattenrör i området. Grundvattenmodellen har i sin tur jordmodellen integrerad samt en vattenbalans som båda är baserade på verkliga data.Resultaten från fallstudien visar att oväntat stora sättningar kan förekomma i områden som ligger relativt långt ifrån källan för grundvattenavsänkningen. Detta motiverar användandet av en sådan metod för att bedöma risken för sättning.Till fallstudien har tre olika jordmodeller använts både i simulering av grundvattenströmning och simulering med den integrerade sättningsmodellen. De tre jordmodellerna skiljer sig åt så till vida att de har olika mängd data som kriginginterpolationen grundar sig på, det innebär också olika datatäthet i jordmodellerna. Detta har till syfte att undersöka hur sättningsbilden påverkas av datatätheten i jordmodellen och även hur grundvattenmodellen respektive sättningsmodellen påverkas. Resultaten i respektive modell påverkas inte mer än försumbart av datatätheten i jordmodellen visar fallstudien. / From a construction engineering point of view groundwater drawdown in a confined aquifer can result in ground subsidence that can damage buildings and constructions. The connection between hydrogeology and soil mechanics is clear, however when estimating ground settlement as a result of groundwater drawdown the estimations are often rough. This is due to that settlement is traditionally calculated with methods that only allow calculations in single points where geotechnical data is estimated. Areas between these points are often left out of the calculations. Groundwater drawdown is seldom simulated with acknowledged software programs like Modflow when estimating groundwater lowering and the affected area.This study combines a groundwater model simulated in Modflow and a soil strata model, interpolated with Kriging, with settlement calculations. This ends up as a an integrated soil settlement model which has the purpose to generate overview maps over areas that are sensitive to settlement as a result of ground water lowering. The integrated model is programmed in Octave for this study. The model is then tested with a case study that uses data from a real construction project in the area of Mälardalen. A hypothetical case of ground water lowering is simulated for the case study. Fundamental hydro-geological theory is used to estimate loads and effective stresses from the lowering of the water table.The result from the integrated model has been validated against calculations of settlement in the software Geosuite Settlement which is an acknowledged method for settlement calculations. This shows that the integrated model calculates settlement with great precision. The modeled initial ground water table is compared with a kriginginterpolated groundwater table which is based on data from ground water pipes in the area. Based on the comparison the initial ground water conditions simulated in Modflow are accepted. This simulated ground water model has the soil model and also a water balance integrated.The results from the case study show that unexpectedly large ground settlements can occur even far from the source of the ground water lowering.For the case study three different soil models are used, both in the ground water model and in the integrated model. The soil models differ in a way that they are based on different amounts of data from which the kriging interpolation is done. The purpose for this is to investigate what effects this might have on the ground water model and the integrated model respectively. The results from these different simulations show insignificantly small differences.

A importância do estudo de aqüíferos livres e rasos se evidencia quando se leva em conta que tais aqüíferos são altamente vulneráveis a contaminação e que em geral se encontram em regiões densamente povoadas. Com o objetivo de melhor compreender a dinâmica do fluxo de água subterrânea na zona saturada destas formações aqüíferas, foi construído e calibrado um detalhado modelo tridimensional e transiente para representar o fluxo subterrâneo em uma porção de um aqüífero livre e raso que ocorre no Parque Ecológico do Tietê, no Município de São Paulo. Este modelo, baseado em dados de cargas hidráulicas de 81 poços de monitoramento instalados em uma área de 320m2 medidas diariamente, em valores diários de precipitação na área e em 75 determinações de condutividades hidráulicas, foi construído e executado através do software Visual MODFLOW e calibrado por meio do software PEST, através de um método de regressão não linear. Os resultados gerados por este modelo para um evento de recarga rápida indicam que a distribuição de cargas hidráulicas na área de estudo é governada pela recarga regional, pela distribuição das condutividades hidráulicas horizontais e pelas vazões específicas das camadas superiores da formação, mas é praticamente independente da intensidade e distribuição temporal da recarga local, que influencia apenas os níveis mais rasos da formação. Além disso, verifica-se o caráter extremamente dinâmico do fluxo, que responde rapidamente ao evento de recarga regional, apresentando alterações significativas de intensidade e direção no intervalo de apenas alguns dias. Tais resultados, verificados apenas por meio de um modelo tão detalhado como o produzido neste trabalho, sugerem que estratégias comuns de avaliação podem não ser capazes de avaliar suficientemente bem o comportamento deste tipo de aqüífero. / The importance of studying unconfined and shallow aquifers is verified when taking into account that these aquifers are highly vulnerable to contamination and that they are often located under heavily populated regions. A detailed three-dimensional and transient model, which represents the groundwater flow in an area of an unconfined and shallow aquifer at the Parque Ecológico do Tietê, in the Municipality of São Paulo, Brazil, was created and calibrated with the objective of enhancing the comprehension of the dynamics of the groundwater flow in the saturated zone of these formations. This model was based on hydraulic head data of 81 monitoring wells, which were installed in an area of 320 m2 and measured on a daily basis, on daily values of precipitation occurring in the area and 75 determinations of hydraulic conductivities. The model was created and run using the software Visual MODFLOW and calibrated using the software PEST, based on a non-linear regression method. The results generated by this method for an event of rapid recharge indicate that the distribution of the horizontal hydraulic conductivities and by the specific yield of the upper layers of the formation, but it is mostly disconnected from the intensity and temporal distribution of the local recharge, that only influences the shallowest levels of the formation. Furthermore, the extremely dynamic character of the flow can be verified, which responds rapidly to the regional recharge event, presenting significant changes of intensity and direction within only a few days. These results, which were only observable through a very detailed model such as the one produced and presented in this study, suggest that commonly employed evaluation strategies may not be capable of evaluating this kind of aquifer in a satisfactory way.

För att uppskatta inläckage till en bergförlagd tunnel finns ett antal analytiska modeller som tar hänsyn till parametrar som hydraulisk konduktivitet, storlek på tunnel och hur djupt tunneln är förlagd. Ett antagande i samtliga analytiska modeller är en homogen och isotrop geologi. I verkligheten är detta sällan fallet, som i exempelvis Kista, strax norr om Stockholm, där en VA-tunnel har drivits under 2017 och inläckaget till tunneln överskrider kravet som baserats på en analytisk modell. En modell utvecklas över områdets geologi och grundvattenmagasin för att simulera grundvattennivåer med och utan inläckage till en bergförlagd tunnel. Avsänkningens utbredning till följd av tunneln är ett s.k. influensområde. Modellen utvecklas i ett relativt nytt modelleringsverktyg, MODFLOW-USG, som utnyttjar ostrukturerade rasternät för att bättre representera tunnelsträckning och öka numerisk stabilitet. Geologin modelleras utifrån jorddjup, bergöveryta och geotekniska undersökningar som utförts i området. Hydraulisk konduktivitet i berget approximeras med hjälp av SGU:s brunnsarkiv. Dessutom utvecklas en modell som innehåller en anslutande tunnel och två st sprickzoner i berget. Jämfört med det influensområde som skattats i tunnelns MKB visar modellen på ett mindre influensområde i vissa riktningar och ett större i andra riktningar. Resultaten bekräftar att influensområdet ökar med ett ökat inläckage, och att en numerisk modell är bättre på att representera områdesspecifik geologi än en enkel vattenbalansmodell. Influensområdet i det undre magasinet i Kista bedöms vara cirka 24,4 ha stort. Enligt den numeriska modellen ger en förändring i inläckage på 0,01 l/min och m tunnel en ökning respektive miskning av influensområdet med cirka 19 %. Vidare visar arbetet på att val av analytisk modell för att skatta inläckage kan påverka skattning med en faktor 3. För att hantera osäkerheter i skattning av inläckage föreslås att flera olika värden på inläckage simuleras för att bedöma influensområdets utbredning och för att skapa en uppfattning om hur storleken på inläckaget påverkar influensområdets utbredning. / To determine leakage into rock tunnels one can choose between a range of analytical models. Common to all analytical models is their inability to reflect complex hydrogeological settings and only taking into account hydraulic conductivity, dimensions of the tunnel and tunnel depth. In Stockholm, Sweden, a rock tunnel built in 2017 currently exceeds maximum allowed leakage which was derived from an analytical model. A model is developed to represent the site’s hydrogeological settings. The model is used to simulate groundwater levels with and without representation of the tunnel, the area of influence being the difference between the two. A Voronoi-cell unstructured grid is implemented in MODFLOW-USG to improve representation of complex geology and increase numerical stability. Geology is modeled using soil depth, bed rock surface and geotechnical investigations. Rock hydraulic conductivity is derived from a well archive. A second model is developed to represent geological fault zones and a connecting tunnel. The model verifies the positive relationship between magnitude of leakage and size of the zone of impact. Results also show an improved approximation of the zone of impact as compared to an earlier approximation derived from a simple water balance. The numerical models predicts a larger zone of impacts in some areas and a smaller zone of impact in other areas. Furthermore, the numerical model shows that a change in leakage of 0,01 l/min and m tunnel can alter the size of the zone of impact by about 19%. By comparing some of the most well-known analytical models to predict leakage it is found that the approximation can differ by a factor 3 using identical parameter values. To reflect the uncertainty in predicted leakage the author recommends simulating a numerical model with varying magnitudes of leakage. This will help in risk assessment and identify a site specific relationship between magnitude of leakage and area of influence.

På en fastighet i Luthagen i Uppsala uppfördes en byggnad grundlagd på träpålar år 1936. Under 60-talet uppstod läckage på en spillvattenledning vilket ledde till en kontaminering av sprickvattenakviferen där träpålar är särskilt utsatta för påverkan av bakterie, svamp- och virusangrepp till följd av torrläggning. För att spola bort kontaminerat sprickvatten och för att hålla en jämn sprickvattennivå i akviferen installerades påfyllningsbrunnar på två platser i källaren under huset. Brunnarna är driva ner i de pålrännor där träpålarna är slagna. Den normala vattentillförseln låg under flera decennier på runt 10 m3/år och vattennivåerna övervakades av bostadsföreningens fastighetstjänst med regelbundna observationer av sprickvattennivån. Under hösten 2016 skedde ett trendbrott och förbrukningen av dricksvattenpåfyllning på över 10 m3/dygn uppmättes. Vattenförbrukningen fortsatte att öka och var vid vissa perioder uppe på nästan 20 m3/dygn. Bjerking AB fick då i uppdrag att undersöka orsaken till den ökade vattenförbrukningen och var vattnet tog vägen. I samband med upptäckten av den ökande vattenförbrukningen skedde även ett brott på en kommunal vattenledning i Kyrkogårdsgatan. Akviferen antogs vara i princip tät och borde inte haft någon hydraulisk kontakt med det vattenförande moränlagret under leran. Syftet med följande rapport är att undersöka de flöden som sprickvattnet har och försöka bestämma den förhärskande flödesriktningen. Den frågeställning som valts är baserad på möjligheten att modellera flödet i grundvattenmodelleringsprogrammet GMS-MODFLOW. Frågeställningen är följande För att sedan kunna modellera grundvattenflödet användes parameter estimation (PEST) som utgår från det framkalibrerade initialt grundvattenförhållandet. Där efter bestäms ett antal zoner som programmet sedan beräknar fram den hydrauliska konduktiviteten för respektive zon. Modelleringen med PEST gav resultat som visar på höga hydrauliska konduktiviteter i husets sydvästra kortsida. Den transienta modelleringen utfördes genom att de observerade vattennivåerna som mätts upp under avstängningsförsöket jämfördes med beräknade värden som fåtts genom modellen. Resultaten visade då på att be beräknade vattennivåerna till viss del stämmde över ens med de observerade, även om de var förskjuta från varandra med nästan 0,6 meter. Denna skillnad berodde antagligen på det gränsvärde som satts för randvillkoret i modellen på 7,25. Vilket gjorde att inga vattennivåer kunde bli högra än detta, vilket som tolkades som att det förekom trösklar i rännorna. Dessa resultat stämmer överens med den hypotes som antogs före projektstarten och innebar att den mest troliga flödesriktningen skulle vara åt sydväst och husets kortsida. Att resultaten bevisar antagandet säger att modellen som byggdes är mer eller mindre rättvisande. Det är dock mycket osäkerheter i modellen och de ingångsvärden som använts. Bland annat är materialtyperna som finns i modellen enbart antagna och det har inte gjorts några bestämningar av markens hydrauliska egenskaper. Andra möjliga fel som kan ha påverkat resultaten är de skalningsproblem som finns i MODFLOW. Slutsatsen är att det är möjligt att göra denna typ av ”småskalig” flödesmodellering i MODFLOW och att förhärskande flödesriktningen är åt sydväst och husets kortsida. / Most of Sweden’s older buildings constructed in clay rich areas are founded on wood poles. The poles are used to build the constructions in areas with soils without satisfying stability, such e.g. clay soils, where the poles are used as the “stable ground” where the building are founded upon. One of the problems with wood poles are that to prevent the wood from rotting the poles must be covered with water, to make an oxygen free environment. The problems start first when the water levels start to decrease due to dewatering or if the water is contaminated with bacteria from e.g. leaky sewer pipes. The contamination leads to decomposing of the wood, which affect the stability of the poles negatively. To prevent the leaky aquifers many houses, have water supply by wells where water is added to the aquifer to keep the water at stable levels. In Uppsala, Sweden an apartment building founded on wood poles have this problems with a contaminated and leaky aquifer. The reason to the contamination was an old sewer pipe that start leaking due to subsidence of the clay below the house. The dewatering of the aquifer was detected in 2016 when the water supply to the aquifer increased from 10th of cubic meters per year to 10th of cubic meters per day. In an attempt to find out what the reason to the high discharge from the aquifer a groundwater model was constructed to modelling the groundwater flows below the building. To modeling this problem the software Groundwater Modeling System (GMS) and MODFLOW was used. By construct a 3D grid of cells in the same dimensions as the building divided in to three different layers it was possible to simulating the groundwater flow through the aquifer. The way we did it was by knowing that as part of the foundation there was “channels” filled with gavel above the poles. The hypothesis was that the water was flowing through this high hydraulic conductivity “gravel channels” and there for we used the model to performed calculations of the hydraulic conductivity in the channels and the areas around the channels. The results told us that there were high conductivity zones in the south west part of the building. After the calculation of the conductivity, a test was performed, where the water supply where turned off and the decrease of the water levels was measured. By using this calculated hydraulic conductivity and the observed levels from the water supply test we let the model calculate the change of water level during the whole-time series. The results told us that there was a possible groundwater flow to the sought west and that some sort of threshold in the channel prevented the water to sink below a curtain limit of 7.65 meter.

Karst aquifers can be conceptualized as dual flow systems comprised of a low-conductive matrix with embedded high-conductive conduits / preferential flow zones. Discharge in conduits ranges from low-velocity laminar flow to high-velocity transitional and turbulent flow. Commonly employed continuum models do not account for the specific behavior of transitional and turbulent flow. In response to this limitation, enhancements have been made to MODFLOW, a commonly used groundwater flow model, by adding a discrete conduit network to the matrix continuum (hybrid model). The Conduit Flow Process (CFP) package is the latest realization of this model approach. CFP Mode 1 (CFPM1) computes laminar and turbulent flow in discrete conduits that are coupled to the laminar continuum model. CFP Mode 2 (CFPM2) accounts for turbulent flow in preferential flow layers by adapting the continuum model. Therefore, laminar hydraulic conduc-tivities are converted into turbulent hydraulic conductivities. CFPM2 was further modified to consider steady turbulent pipe flow. Karst models based on CFPM2 require potentially less input data and computational efforts than karst models based on CFPM1. Furthermore, CFPM2 integrates more easily into MODFLOW versions including e.g. transport models. Parameter studies for a synthetic catchment demonstrates that continuum models with turbulent flow representation and an additional flow barrier between conduits and matrix can represent karst systems similar to hybrid models. For simulation of highly transient flow processes in karst conduit systems, i.e. during flood events, it is crucial to consider dynamics such as free-surface flow, wave propagation, and changes between pressurized and non-pressurized conduit flow. The coupled overland- and groundwater flow model MODBRANCH was therefore enhanced to consider unsteady and non-uniform flow processes in karst conduits. Flow in discrete conduits is simulated using the Saint-Venant-equations for free-surface flow. Contrary to overland flow, the cross sectional area of karst conduits is finite. Accordingly, both pressurized and non-pressurized flow may occur within conduits. To simulate pressurized flow, a hypothetical, narrow, open-top slot (Preissmann slot) is added to the conduit crown, which allows the use of the free-surface flow equations for fully filled conduits. Beyond this, the model features a variable time step to consider wave speed variations, for example due to the transition from free-surface to pressurized flow. Parameter studies for a synthetic catchment demonstrate the significance of free-surface flow representation for variably filled conduits. / Karstgrundwasserleiter können als duale Fließsysteme konzeptionalisiert werden, bestehend aus einer geringdurchlässigen Matrix mit eingebundenen hochdurchlässigen Bereichen, z. B. Karströhren. Der Abfluss in den hochdurchlässigen Bereichen reicht von langsamer laminarer Strömung bis zu schneller turbulenter Strömung. Herkömmliche numerische Grundwasser-strömungsmodelle berücksichtigen nicht die spezifischen Eigenschaften von nicht-laminarer Strömung (Übergangsbereich laminar-turbulent bzw. turbulente Verhältnisse). Ein Ansatz um diese Einschränkung zu umgehen, ist die Erweiterung des laminaren Kontinuums um ein dis-kretes Röhrenmodell, das zustandsabhängig laminare und turbulente Strömung berücksichtigt (Hybridmodell). Eine aktuelle Umsetzung dieses Ansatzes ist Conduit Flow Process (CFP), ein Modul für das weitverbreitete Grundwasserströmungsmodell MODFLOW. CFP Mode 1 (CFPM1) berechnet laminare und turbulente Strömung in diskreten, mit dem Kontinuummodell gekoppelten Röhren. CFP Mode 2 (CFPM2) berücksichtigt nicht-laminare Strömung in hochdurchlässigen Schichten mit einer angepassten hydraulischen Leitfähigkeit des Kontinuummodells. CFPM2 wurde weiter modifiziert, so dass auch turbulente Strömung in Karströhren berechnet werden kann. Dadurch kann möglicherweise der Parameterbedarf sowie der Rechenaufwand gegenüber Hybrid¬modellen reduziert werden. CFPM2 lässt sich einfach in vorhandene MODFLOW Modelle einbinden, z. B. zur Berechnung von Transportprozessen. Parameterstudien für ein idealisiertes Karsteinzugsgebiet zeigen, dass Kontinuummodelle bei Berücksichtigung der turbulenten Strömung sowie des zusätzlichen hydraulischen Widerstand zwischen Röhren und Matrix, Karstsysteme ähnlich wie Hybridmodelle darstellen. Zur Simulation von instationären Prozessen in Karströhren, z. B. ausgeprägte Abflusssignale infolge pulsförmiger Grundwasserneubildung, ist es notwendig, dynamische Prozesse infolge Freispiegelabfluss, Wellenausbreitung sowie Wechsel zwischen Abfluss in teil- und vollgefüllten Röhren zu berücksichtigen. Aus diesem Grund wurde das numerische Modell MODBRANCH, welches ein diskretes Oberflächenwassermodell mit einem Kontinuummodell koppelt, so angepasst, dass instationäre und nichtgleichförmige Abflussprozesse in Karströhren berücksichtigt werden können. Der Abfluss in diskreten Röhren wird dabei mit den Saint-Venant-Gleichungen für Freispiegelabfluss berechnet. Im Gegensatz zu Oberflächengewässern ist der für den Abfluss zur Verfügung stehende Querschnitt in Karströhren limitiert, so dass sowohl Freispiegel- als auch Druckabfluss innerhalb der Röhren auftreten kann. Druckabfluss wird mit Hilfe eines schmalen virtuellen Schlitzes an der Röhrenoberkante simuliert (Preissmann Schlitz), der auch im Fall vollgefüllter Röhren die Anwendung der Gleichungen für Freispiegelabfluss erlaubt. Durch die Verwendung eines variablen Zeitschrittes kann die geänderte Dynamik beim Übergang von Freispiegel- zu Druckabfluss berücksichtigt werden. Parameterstudien für idealisierte, synthetische Karsteinzugsgebiete demonstrieren die Bedeutung der Berücksichtigung von Freispiegelabfluss in teilgefüllter Röhren.

Karst aquifers can be conceptualized as dual flow systems comprised of a low-conductive matrix with embedded high-conductive conduits / preferential flow zones. Discharge in conduits ranges from low-velocity laminar flow to high-velocity transitional and turbulent flow. Commonly employed continuum models do not account for the specific behavior of transitional and turbulent flow. In response to this limitation, enhancements have been made to MODFLOW, a commonly used groundwater flow model, by adding a discrete conduit network to the matrix continuum (hybrid model). The Conduit Flow Process (CFP) package is the latest realization of this model approach. CFP Mode 1 (CFPM1) computes laminar and turbulent flow in discrete conduits that are coupled to the laminar continuum model. CFP Mode 2 (CFPM2) accounts for turbulent flow in preferential flow layers by adapting the continuum model. Therefore, laminar hydraulic conduc-tivities are converted into turbulent hydraulic conductivities. CFPM2 was further modified to consider steady turbulent pipe flow. Karst models based on CFPM2 require potentially less input data and computational efforts than karst models based on CFPM1. Furthermore, CFPM2 integrates more easily into MODFLOW versions including e.g. transport models. Parameter studies for a synthetic catchment demonstrates that continuum models with turbulent flow representation and an additional flow barrier between conduits and matrix can represent karst systems similar to hybrid models. For simulation of highly transient flow processes in karst conduit systems, i.e. during flood events, it is crucial to consider dynamics such as free-surface flow, wave propagation, and changes between pressurized and non-pressurized conduit flow. The coupled overland- and groundwater flow model MODBRANCH was therefore enhanced to consider unsteady and non-uniform flow processes in karst conduits. Flow in discrete conduits is simulated using the Saint-Venant-equations for free-surface flow. Contrary to overland flow, the cross sectional area of karst conduits is finite. Accordingly, both pressurized and non-pressurized flow may occur within conduits. To simulate pressurized flow, a hypothetical, narrow, open-top slot (Preissmann slot) is added to the conduit crown, which allows the use of the free-surface flow equations for fully filled conduits. Beyond this, the model features a variable time step to consider wave speed variations, for example due to the transition from free-surface to pressurized flow. Parameter studies for a synthetic catchment demonstrate the significance of free-surface flow representation for variably filled conduits. / Karstgrundwasserleiter können als duale Fließsysteme konzeptionalisiert werden, bestehend aus einer geringdurchlässigen Matrix mit eingebundenen hochdurchlässigen Bereichen, z. B. Karströhren. Der Abfluss in den hochdurchlässigen Bereichen reicht von langsamer laminarer Strömung bis zu schneller turbulenter Strömung. Herkömmliche numerische Grundwasser-strömungsmodelle berücksichtigen nicht die spezifischen Eigenschaften von nicht-laminarer Strömung (Übergangsbereich laminar-turbulent bzw. turbulente Verhältnisse). Ein Ansatz um diese Einschränkung zu umgehen, ist die Erweiterung des laminaren Kontinuums um ein dis-kretes Röhrenmodell, das zustandsabhängig laminare und turbulente Strömung berücksichtigt (Hybridmodell). Eine aktuelle Umsetzung dieses Ansatzes ist Conduit Flow Process (CFP), ein Modul für das weitverbreitete Grundwasserströmungsmodell MODFLOW. CFP Mode 1 (CFPM1) berechnet laminare und turbulente Strömung in diskreten, mit dem Kontinuummodell gekoppelten Röhren. CFP Mode 2 (CFPM2) berücksichtigt nicht-laminare Strömung in hochdurchlässigen Schichten mit einer angepassten hydraulischen Leitfähigkeit des Kontinuummodells. CFPM2 wurde weiter modifiziert, so dass auch turbulente Strömung in Karströhren berechnet werden kann. Dadurch kann möglicherweise der Parameterbedarf sowie der Rechenaufwand gegenüber Hybrid¬modellen reduziert werden. CFPM2 lässt sich einfach in vorhandene MODFLOW Modelle einbinden, z. B. zur Berechnung von Transportprozessen. Parameterstudien für ein idealisiertes Karsteinzugsgebiet zeigen, dass Kontinuummodelle bei Berücksichtigung der turbulenten Strömung sowie des zusätzlichen hydraulischen Widerstand zwischen Röhren und Matrix, Karstsysteme ähnlich wie Hybridmodelle darstellen. Zur Simulation von instationären Prozessen in Karströhren, z. B. ausgeprägte Abflusssignale infolge pulsförmiger Grundwasserneubildung, ist es notwendig, dynamische Prozesse infolge Freispiegelabfluss, Wellenausbreitung sowie Wechsel zwischen Abfluss in teil- und vollgefüllten Röhren zu berücksichtigen. Aus diesem Grund wurde das numerische Modell MODBRANCH, welches ein diskretes Oberflächenwassermodell mit einem Kontinuummodell koppelt, so angepasst, dass instationäre und nichtgleichförmige Abflussprozesse in Karströhren berücksichtigt werden können. Der Abfluss in diskreten Röhren wird dabei mit den Saint-Venant-Gleichungen für Freispiegelabfluss berechnet. Im Gegensatz zu Oberflächengewässern ist der für den Abfluss zur Verfügung stehende Querschnitt in Karströhren limitiert, so dass sowohl Freispiegel- als auch Druckabfluss innerhalb der Röhren auftreten kann. Druckabfluss wird mit Hilfe eines schmalen virtuellen Schlitzes an der Röhrenoberkante simuliert (Preissmann Schlitz), der auch im Fall vollgefüllter Röhren die Anwendung der Gleichungen für Freispiegelabfluss erlaubt. Durch die Verwendung eines variablen Zeitschrittes kann die geänderte Dynamik beim Übergang von Freispiegel- zu Druckabfluss berücksichtigt werden. Parameterstudien für idealisierte, synthetische Karsteinzugsgebiete demonstrieren die Bedeutung der Berücksichtigung von Freispiegelabfluss in teilgefüllter Röhren.

The Mt. Tom Price Mine, located in the Pilbara region of Western Australia, has been the site of major iron ore mining since the 1960s by Rio Tinto Iron Ore/Pilbara Iron. The thesis project area covers approximately 121 km², covering the Mt. Tom Price Mining area and the surrounding catchment boundary. The climate in the Pilbara region is arid, with rainfall driven by seasonal cyclonic events, producing 300 mm/year net rainfall on average. The geology of the Mt. Tom Price area consists of a series of banded iron formations (BIF) and shales that are generally low in hydraulic conductivity values. Iron ore in the region is produced through the process of supergene enrichment whereby gangue minerals are dissolved and replaced with haematite and goethite. Mining is focused in a series of open cast pits including, North Deposit, West Pits, Centre Pits, Southern Ridge, South East Prongs, Section Six, Section Seven, and the proposed Marra Mamba Pits. Due to the impermeable nature and complex geology of the BIF sequence, groundwater flow is dominated by bedrock aquifer flow, with compartmentalization occurring in several areas of the mine. Highly faulted and folded units can also have increased hydraulic conductivity values. Pit floor lowering began to encounter the regional water table in early 1994. A series of dewatering bores and depressurization measurements have been utilized to ensure dry mining practice. This data was used to help understand regional groundwater flow and create the Mt. Tom Price Groundwater Model (MTPGM). A 3D geological model of the project area was created to aid visualisation of semi-regional hydrogeology. From this model, accurate template files were created so that geological detail loss is kept to a minimal when entering hydrogeological parameters into the MTPGM. The MTPGM was setup using PMWIN Pro, a graphical user interface for use with MODFLOW. Stresses such as recharge and pumping were entered via software packages within MODFLOW. The model was run to simulate measured 1994-2007 responses to dewatering and high rainfall events. A Parameter Estimation (PEST) software package and trial and error calibration was used to lower stress response variances that were observed in the model output files. This was achieved by the adjustment of hydrogeological parameters such as hydraulic conductivity and specific yield values. A prediction simulation of final pit lake recovery was created Using the calibrated MTPGM. Recovery curves predicted that full recovery of the water table of the pit voids varied from 96 to 120 years, recovering to levels close to the initial heads measured in 1994 before large-scale pumping commenced. The hydrochemistry of the groundwater in the mining area is highly influenced by geological hosts, with clearly defined hydrochemical signatures approximated for each screened geological unit. Due to the sulphur rich, acid- forming Mt. McRae Shale, regular monitoring of pit and groundwater is essential. Final pit lake water quality was estimated using final pit levels and recovery rates approximated from the MTPGM, combined with historical data and previous groundwater quality reports. Pit lake water quality is dominantly driven by evaporation concentration, caused by high evaporation rates and low throughflow. Pit waters are expected to be brine waters (greater than 100,000 mg/L TDS), with high levels of acidity values occurring in the South East Prongs and Section Six pits due to the exposure of the acid forming Mt. McRae Shale above the pit lakes at these localities. Future studies should focus on more detailed modelling of the compartmentalised aquifer systems. This would produce much more accurate final pit lake levels. Further study of the Mt. McRae Shale formation and its implications on acidity should also be undertaken. Seasonal fluctuations in lake levels will affect acidity due to the continual re-exposure and oxidation of the Mt. McRae Shale. This could be studied to help understand short term pit lake quality conditions and help to predict long term acidity conditions in the pit lakes.