Modelling Surface Water-Groundwater Exchange
Total Page：16
File Type：pdf, Size：1020Kb
Modelling Surface Water-Groundwater Exchange Evaluating Model Uncertainty from the Catchment to Bedform-Scale Dissertation der Mathematisch-Naturwissenschaftlichen Fakultät der Eberhard Karls Universität Tübingen zur Erlangung des Grades eines Doktors der Naturwissenschaften (Dr. rer. nat.) vorgelegt von Reynold Chow, M.Sc., P.Geo aus Toronto, Kanada Tübingen 2019 Gedruckt mit Genehmigung der Mathematisch-Naturwissenschaftlichen Fakultät der Eberhard Karls Universität Tübingen. Tag der mündlichen Qualifikation: 21.05.2019 Dekan: Prof. Dr. Wolfgang Rosenstiel 1. Berichterstatter: Prof. Dr.-Ing. Wolfgang Nowak 2. Berichterstatter: Dr. rer. nat. Thomas Wöhling 3. Berichterstatter: Prof. Dr.-Ing. Olaf A. Cirpka Abstract My dissertation focuses on evaluating model uncertainties when numerically simulating surface water-groundwater (sw-gw) interactions at different scales. To do so, I mainly use HydroGeoSphere, a physically-based distributed finite element model that fully couples variably saturated subsurface flow with surface water flow. I evaluate three predominant model uncertainty types at three different scales of sw-gw interaction. For each of these investigations, I selected a corresponding study site. Firstly, I evaluate structural (conceptual) uncertainty from delineating baseflow contribution areas to gaining stream reaches, or stream capture zones, at the catchment-scale. I investigate how the delineated stream capture zone (in the Alder Creek watershed) can differ due to the chosen model code and delineation method (Chow et al. 2016). The results indicate that different models can calibrate acceptably well to the same data and produce very similar distributions of hydraulic head, but can produce different capture zones. The stream capture zone is highly sensitive to the post-processing particle tracking algorithm. Reverse transport is an alternative and more reliable approach that accounts for local-scale parameter uncertainty and provides probability intervals for the stream capture zone. The two methods can be combined to enhance the overall confidence in the delineated stream capture zone. Secondly, I evaluate parameter uncertainty when simulating meander-scale hyporheic exchange by conducting a model-based sensitivity (bathymetry) and uncertainty (subsurface K-distribution) analysis. I select the Steinlach River Test Site for demonstration. I conduct a sensitivity analysis to determine the aspects of river bathymetry that have the greatest influence on the predictive biases (Chow et al. 2018). Results indicate that simulating hyporheic exchange with a high-resolution detailed bathymetry using a 3D fully coupled sw-gw model leads to nested multi-scale hyporheic exchange systems. A poorly resolved bathymetry will underestimate the small-scale hyporheic exchange, biasing the simulated hyporheic exchange towards larger scales, thus leading to overestimates of hyporheic exchange residence times. The detailed river slope alone is not enough to accurately simulate the locations and magnitudes of losing and gaining river reaches. Thus, local bedforms in terms of bathymetric highs and lows within the river are required. Incorporating local bedforms will likely capture the nested nature of hyporheic exchange, leading to more physically meaningful simulations of hyporheic exchange fluxes and transit times. v Additionally, I conduct an uncertainty analysis to evaluate the trade-offs between intrinsic (irreducible) and epistemic (reducible) model errors when choosing between homogeneous and highly complex subsurface parameter structures (Chow et al. 2019). Results indicate that, if the parameter structure is too simple, it will be limited by intrinsic model errors. By increasing subsurface complexity through the addition of zones or heterogeneity, we can begin to exchange intrinsic for epistemic errors. Thus, choosing the appropriate level of detail to represent the subsurface parameter distributions depends on the acceptable range of intrinsic structural errors for the given modelling objectives and the available site data. Thirdly, I evaluate data uncertainty at the bedform-scale in a fractured rock setting by modelling a conservative tracer experiment at the Eramosa Bedrock River Site. I use a stochastic discrete fracture network framework to represent the subsurface fractured bedrock connectivity, and in doing so produce a probabilistic distribution of the potential hyporheic exchange extents and residence times. The results indicate that the coincidence of fractures and hydraulic gradients determine the spatial extents of bedform-scale hyporheic exchange. Furthermore, hyporheic exchange residence times in bedrock rivers at the bedform-scale are potentially orders of magnitude longer when compared to fluvial rivers (i.e., months to years vs. minutes to hours). In the age of highly-parameterized integrated hydrological models, there is an increasing need to understand whether we are getting the right forecasts for the right reasons. Modelling is only a single tool in the scientific toolbox that works best when combined with others, e.g., field and laboratory experiments. The next generation of integrated hydrological modellers will face evermore challenging objectives, which to achieve will require multi-disciplinary teams. Regardless of the origins of your knowledge-base, it is important to always approach numerical models with a healthy level of skepticism. Ultimately, scientists should approach models with a relentless fearlessness to falsify them, which is an absolute necessity if we wish to continually push the envelope of scientific knowledge. vi Zusammenfassung Meine Dissertation beschäftigt sich mit der Einschätzung von Modellunsicherheiten bei der numerischen Simulation von Oberflächenwasser-Grundwasser-Wechselwirkungen auf unterschiedlichen Skalen. Dafür benutze ich hauptsächlich HydroGeoSphere, ein physikalisch- basiertes finite Elemente Modell, das die veränderlich gesättigte Untergrundströmung mit der Oberflächenströmung vollständig koppelt. Ich erkunde drei vorherrschende Unsicherheitstypen auf drei Skalen, indem ich die Oberflächenwasser-Grundwasser- Wechselwirkungen simuliere. Für jede dieser Untersuchungen wähle ich ein geeignetes Untersuchungsgebiet aus. Zunächst untersuche ich die strukturbedingte (konzeptionelle) Unsicherheit auf der Einzugsgebietsskala bei der Bestimmung von Basisabflussbeitragsflächen exfiltrierender Flussabschnitte oder Flusseinzugsgebiete. Dafür wähle ich das Einzugsgebiet des “Alder Creek” aus. Ich erforsche, inwieweit die so bestimmten Flusseinzugsgebiete von der Wahl des Modellcodes und der Abgrenzungsmethode abhängen (Chow et al. 2016). Die Ergebnisse deuten darauf hin, dass verschiedene Modelle hinreichend gut auf dieselben Daten kalibriert werden können, wobei sehr ähnliche Verteilungen des hydraulischen Potentials, aber unterschiedliche Geometrien der Einzugsgebiete entstehen. Das Flusseinzugsgebiet ist sehr sensitiv gegenüber der Wahl des Particle-Tracking-Nachverarbeitungsalgorithmus. Umgekehrter Transport stellt einen alternativen und zuverlässigeren Ansatz dar, der lokale Parameterunsicherheit berücksichtigt und Wahrscheinlichkeitsintervalle für das Flusseinzugsgebiet liefert. Beide Methoden können kombiniert werden, um die gesamte Konfidenz des abgeleiteten Flusseinzugsgebiets zu erhöhen. Als Zweites widme ich mich der Parameterunsicherheit, indem ich eine Sensitivitäts- (Bythymetrie) und Unsicherheitsanalyse (K-Verteilung des Untergrunds) für Simulationen des hyporheischen Austauschs auf der Mäanderskala durchführe. Dies untersuche ich am Beispiel des Steinlach-Testgeländes. Ich vollziehe eine Sensitivitätsanalyse um diejenigen Aspekte der Flussbathymetrie zu identifizieren, die den größten Einfluss auf systematische Vorhersageabweichungen haben (Chow et al. 2018). Die Ergebnisse zeigen, dass die Simulation des hyporheischen Austausches mit einer hochaufgelösten und detaillierten Bathymetrie unter Verwendung eines dreidimensionalen vollständig gekoppelten Oberflächenwasser-Grundwasser-Modells zu verschachtelten Multi-Skalen-Systemen von hyporheischem Austausch führt. Eine schlecht aufgelöste Bathymetrie führt zur vii Unterschätzung von kleinskaligem hyporheischen Austausch, wodurch derselbe auf größere Skalen verschoben wird. Dies führt zur Überschätzung von Aufenthaltszeiten in der hyporheischen Zone. Eine detaillierte Beschreibung der Flusssteigung allein reicht dabei nicht aus, um die Lage und den Betrag des Wasseraustauschs von exfiltrierenden und infiltrierenden Flussabschnitten richtig zu simulieren. Folglich sind lokale Flussgeometrien (Sohlformen) in Form von bathymetrischen Hoch- und Tiefpunkten des Flussverlaufs notwendig. Bei Berücksichtigung dieser Sohlform wird die verschachtelte Natur des hyporheischen Austauschs besser erfasst, wodurch physikalisch sinnvollere Simulationen von Austauschflüssen und Verweilzeiten erreicht werden können. Zusätzlich führe ich eine Unsicherheitsanalyse durch, um den Trade-off zwischen intrinsischen (nicht reduzierbaren) und epistemischen (reduzierbaren) Modellfehlern abzuwägen, der bei der Wahl zwischen homogenen und hoch-komplexen Untergrundparameterstrukturen entsteht (Chow et al. 2019). Die Ergebnisse deuten darauf hin, dass die Aussagekraft eines Modells bei einer zu simplen Parameterstruktur durch intrinsische Fehler eingeschränkt ist. Mit zunehmender Untergrundkomplexität durch die Berücksichtigung von Zonen oder Heterogenität können intrinsische Fehler durch epistemische Fehler ersetzt werden. Daher hängt die Wahl