The styles of deformation of the fore-arc wedges along the Chilean convergent margin are observed to vary significantly, despite similar plate kinematic conditions. Here, I focus on the analysis of fore-arc deformation on two regions along the Chilean convergent margin at 20°-24°S and 37°-42°S. Although both regions are subjected to the oblique subduction of the oceanic Nazca plate and backstopped by the Andes mountain chain; they display different patterns of deformation. The northern Chilean study area (20° - 24°S) is characterized by an exceptionally thick crust of about 60 km beneath the Altiplano – Puna plateau, lack of an accretionary wedge in the fore-arc due to hyperarid climate, and consequently a sediment starved trench. Two major margin parallel strike slip faults are observed in this area, the Atacama Fault Zone (AFZ) and the Precordilleran Fault System (PFS). Both strike-slip faults do not exhibit significant recent displacement. The southern study area (37° - 42°S), compared to the northern study area, is characterized by lower topography, high precipitation rates (~2000 mm/yr), and a younger subducted oceanic plate. An active strike-slip fault, the Liquiñe-Ofqui-Fault-Zone (LOFZ), shows ~1 cm/yr recent dextral movement and shapes the surface of this area. Thus, the southern Chilean study area exhibits localized strike-slip motion. Within this area the largest earthquake ever recorded, the 1960 Valdivia earthquake, occurred with a moment magnitude of MW=9.5. I have constructed 2D thermal models and 3D mechanical models for both Chilean study areas to study processes related to active subduction. The applied numerical method is the finite element technique by means of the commercial software package ABAQUS. The thermal models are focused on the thermal conditions along the plate interface. The thermal structure along the plate interface reveals the limits of coupling but also the type of transition from coupled to uncoupled and vice versa. The model results show that shear heating at the plate interface is an important mechanism that should be taken into account. The models also show that the thermal condition at the downdip limit of the coupling zone leads to a sharp decrease of friction along the interface. Due to the different geometries of the two Chilean study areas, such as the slab dip and the thickness of the continental crust, the downdip limit of the southern study area is slightly shallower than that of the northern study area. The results of the 2D thermal models are used to constrain the spatial extent of the coupling zone in the 3D mechanical models. 3D numerical simulations are used to investigate how geometry, rheology and mechanical parameters influence strain partitioning and styles of deformation in the Chilean fore-arc. The general outline of the models is based on the fore-arc geometry and boundary conditions as derived from geophysical and geological field data. I examined the influence of different rheological approaches and varying physical properties of the fore-arc to identify and constrain the parameters controlling the difference in surface deformation between the northern and southern study area. The results of numerical studies demonstrate that a small slab dip, a high coefficient of basal friction, a high obliquity of convergence, and a high Young’s modulus favour localisation of deformation in the fore-arc wedge. This parameter study helped me to constrain preferred models for the two Chilean study areas that fit to first order observations. These preferred models explain the difference in styles of deformation as controlled by the angle of obliquity, the dip of subducting slab, and the strength of wedge material. The difference in styles can be even larger if I apply stronger coupling between plates within the southern area; however, several independent observations indicate opposite tendency showing southward decrease of intensity of coupling. The weaker wedge material of the preferred model for the northern study area is associated with advanced development of the adjacent orogen, the Central Andes. Analysis of world-wide examples of oblique subduction zones supports the conclusion that more mature subduction zones demonstrate less pronounced localization of strike-slip motion. / Die Deformationsmuster der Fore-Arc Keile entlang des chilenischen konvergenten Plattenrandes variieren beachtlich, trotz ähnlicher plattenkinematischer Randbedingungen. In dieser Arbeit konzentriere ich mich auf die Analyse der Deformation des Fore-Arcs in zwei Gebieten entlang des chilenischen konvergenten Plattenrandes zwischen 20°-24°S und 37°-42°S. Obwohl beide Gebiete durch schiefe Subduktion der ozeanischen Nazca Platte und der östlichen Begrenzung durch die Andine Gebirgskette gekennzeichnet sind, zeigen sie unterschiedliche Deformationsmuster an der Oberfläche. Das nördliche chilenische Gebiet (20° - 24°S) ist gekennzeichnet durch eine außergewöhnliche Krustendicke von ~ 60 km unterhalb des Altiplano - Puna Plateaus, dem Fehlen eines akkretionären Prismas im Fore-Arc aufgrund des trockenen Klimas und somit einer nahezu sedimentfreien Tiefseerinne. Zwei große Plattenrand-parallele Strike-Slip Störungen werden in diesem Gebiet beobachtet, die Atacama Fault Zone (AFZ) und das Precordilleran Fault System (PFS). Beide Strike-Slip Störungen zeigen keine signifikanten aktuellen Bewegungsraten. Das südliche Gebiet (37° - 42°S) ist im Vergleich zum nördlichen Gebiet durch eine niedrigere Topographie, hohe Niederschlagsraten (~2000 mm/a) und eine jüngere abtauchende ozeanische Platte gekennzeichnet. Die aktive Strike-Slip Störung, Liquiñe-Ofqui-Fault-Zone (LOFZ), ist gekennzeichnet durch aktuelle dextrale Bewegungsraten von 1 cm/a und prägt die Oberflächenstruktur in dieser Region. Folglich ist der südliche Arbeitsbereich durch lokalisierte Strike-Slip Bewegung charakterisiert. Innerhalb dieses Gebietes ereignete sich das größte instrumentell aufgezeichnete Erdbeben, das 1960 Valdivia Erdbeben, mit einer Stärke von MW=9.5. 2D thermische Modelle und 3D mechanische Modelle wurden für die beiden chilenischen Gebiete konstruiert, um Prozesse im Zusammenhang mit aktiver Subduktion zu untersuchen. Als numerisches Verfahren wurde die Finite Elemente Methode mit Hilfe des kommerziellen Softwarepakets ABAQUS angewandt. Die thermischen Modelle sind auf die thermischen Konditionen entlang der Plattengrenzfläche fokussiert. Die thermische Struktur entlang der Plattengrenzfläche zeigt die Grenzen der Kopplung an aber auch die Art des Überganges von gekoppelt zu nicht gekoppelt und umgekehrt. Die Modellergebnisse zeigen, dass Heizen infolge der Scherung an der Plattengrenzfläche ein wichtiger Faktor ist, der in Betracht gezogen werden sollte. Die Modelle zeigen auch, dass die thermische Struktur an der unteren Begrenzung der Koppelzone zu einer deutlichen Abnahme der Reibung entlang der Grenzfläche führt. Aufgrund der unterschiedlichen Geometrien der zwei chilenischen Untersuchungsgebiete, z.B. Abtauchwinkel der ozeanischen Platte und Krustendicke, ist die untere Begrenzung der Koppelzone des südlichen Untersuchungsgebietes in geringerer Tiefe als die des nördlichen Gebietes. Die Ergebnisse der thermischen 2D Modelle werden genutzt, um die räumliche Ausdehnung der Koppelzone in den mechanischen 3D Modellen festzulegen. Numerische 3D Simulationen werden genutzt, um zu verstehen, wie Geometrien, Rheologien und mechanische Parameter die Verformungspartitionierung und das Deformationsmuster im chilenischen Fore-Arc beeinflussen. Ich habe den Einfluss unterschiedlicher rheologischer Ansätze und unterschiedlicher physikalischer Eigenschaften auf den Fore-Arc untersucht, um Parameter zu identifizieren und zu bestimmen, die den Unterschied des Deformationsmusters zwischen dem nördlichen und südlichen Gebiet steuern. Die Ergebnisse der numerischen Studien stellen heraus, dass ein kleinerer Abtauchwinkel der ozeanischen Platte, ein hoher basaler Reibungskoeffizient, eine hohe Konvergenzschiefe und ein großer Elastizitätsmodul die Lokalisierung der Deformation im Fore-Arc Keil begünstigen. Basierend auf dieser Parameterstudie habe ich Modelle für die beiden chilenischen Gebiete ausgewählt, die in Beobachtungen erster Ordnung übereinstimmen. Diese ausgewählten Modelle erklären die unterschiedlichen Deformationsmuster durch eine größere Konvergenzschiefe, einen kleineren Abtauchwinkel der ozeanischen Platte und ein härteres Keilmaterial für das südliche Untersuchungsgebiet. Der Unterschied bezüglich der Deformationsmuster kann sogar größer sein, wenn ich eine größere Reibung zwischen den Platten im südlichen Gebiet anwende; jedoch zeigen einige unabhängige Beobachtungen eine umgekehrte Tendenz: eine Abnahme der Intensität der Koppelung von Nord nach Süd. Das schwächere Keilmaterial des ausgewählten Modells für das nördliche Untersuchungsgebiet steht im Zusammenhang mit der fortgeschrittenen Entwicklung des angrenzenden Orogens, der zentralen Anden. Die Analyse weltweiter Beispiele von schiefen Subduktionzonen unterstützt die Schlussfolgerung, dass ältere Subduktionzonen weniger ausgeprägte Lokalisierung von Strike-Slip Bewegung aufzeigen.

Da geologische Störungen können als Grundwasserleiter, -Barrieren oder als gemischte leitende /stauende Fluidsysteme wirken. Aufgrund dessen können Störungen maßgeblich den Grundwasserfluss im Untergrund beeinflussen, welcher deutliche Veränderungen des tiefen thermischen Feldes bewirken kann. Grundwasserdynamik und Temperaturveränderungen sind wiederum entscheidende Faktoren für die Exploration geothermischer Energie. Diese Studie untersuchte den Einfluss von Störungen auf das Fluidsystem und das thermische Feld im Untergrund. Sie erforschte die physikalischen Prozesse, welche das Fluidverhalten und die Temperaturverteilung in Störungen und in den umgebenden Gesteinen. Dazu wurden 3D Finite Elemente Simulationen des gekoppelten Fluid und Wärmetransports für synthetische sowie reale Modelszenarien auf unterschiedlichen Skalen durchgeführt. Um den Einfluss einer schräg einfallenden Störung systematisch durch die schrittweise Veränderung der hydraulischen Öffnungsweite und der Permeabilität, zu untersuchen, wurde ein klein-skaliges synthetisches Modell entwickelt. Ein inverser linearer Zusammenhang wurde festgestellt, welcher zeigt, dass sich die Fluidgeschwindigkeit in der Störung jeweils um ~1e-01 m/s verringert, wenn die Öffnungsweite der Störung um jeweils eine Magnitude vergrößert wird. Ein hoher Permeabilitätskontrast zwischen Störung und umgebender Matrix begünstigt die Fluidadvektion hin zur Störung und führt zu ausgeprägten Druck- und Temperaturveränderungen innerhalb und um die Störung herum. Bei geringem Permeabilitätskontrast zwischen Störung und umgebendem Gestein findet hingegen kein Fluidfluss in der Störung statt, wobei das hydrostatische Druck- sowie das Temperaturfeld unverändert bleiben. Auf Grundlage der synthetischen Modellierungsergebnisse wurde der Einfluss von Störungen auf einer größeren Skala anhand eines komplexeren (realen) geologischen Systems analysiert. Dabei handelt es sich um ein 3D Modell des Geothermiestandortes Groß Schönebeck, der ca. 40 km nördlich von Berlin liegt. Die Integration von einer permeablen und drei impermeablen Hauptstörungen, zeigte unterschiedlich starke Einflüsse auf Fluidzirkulation, Temperatur – und Druckfeld. Die modellierte konvektive Zirkulation in der permeablen Störung verändert das thermische Feld stark (bis zu 15 K). In den gering durchlässigen Störungen wird die Wärme ausschließlich durch Diffusion geleitet. Der konduktive Wärmetransport beeinflusst das thermische Feld nicht, bewirkt jedoch lokale Veränderungen des hydrostatischen Druckfeldes. Um den Einfluss großer Störungszonen mit kilometerweitem vertikalen Versatz auf das geothermische Feld der Beckenskala zu untersuchen, wurden gekoppelte Fluid- und Wärmetransportsimulationen für ein 3D Strukturmodell des Gebietes Brandenburg durchgeführt (Noack et al. 2010; 2013). Bezüglich der Störungspermeabilität wurden verschiedene geologische Szenarien modelliert, von denen zwei Endgliedermodelle ausgewertet wurden. Die Ergebnisse zeigten, dass die undurchlässigen Störungen den Fluidfluss nur lokal beeinflussen. Da sie als hydraulische Barrieren wirken, wird der Fluidfluss mir sehr geringen Geschwindigkeiten entlang der Störungen innerhalb eines Bereichs von ~ 1 km auf jeder Seite umgelenkt. Die modellierten lokalen Veränderungen des Grundwasserzirkulationssystems haben keinen beobachtbaren Effekt auf das Temperaturfeld. Hingegen erzeugen permeable Störungszonen eine ausgeprägte thermische Signatur innerhalb eines Einflussbereichs von ~ 2.4-8.8 km in -1000 m Tiefe und ~6-12 km in -3000 m Tiefe. Diese thermische Signatur, in der sich kältere und wärmere Temperaturbereiche abwechseln, wird durch auf- und abwärts gerichteten Fluidfluss innerhalb der Störung verursacht, der grundsätzlich durch existierende Gradienten in der hydraulischen Druckhöhe angetrieben wird. Alle Studien haben gezeigt, dass Störungen einen beachtlichen Einfluss auf den Fluid-, und Wärmefluss haben. Es stellte sich heraus, dass die Permeabilität in der Störung und in den umgebenden geologischen Schichten so wie der spezifische geologische Rahmen entscheidende Faktoren in der Ausbildung verschiedener Wärmetransportmechanismen sind, die sich in Störungen entwickeln können. Die von permeablen Störungen verursachten Temperaturveränderungen können lokal, jedoch groß sein, genauso wie die durch hydraulisch leitende und nichtleitende Störungen hervorgerufenen Veränderungen des Fluidystems. Letztlich haben die Simulationen für die unterschiedlich skalierten Modelle gezeigt, dass die Ergebnisse sich nicht aufeinander übertragen lassen und dass es notwendig ist, jeden geologischen Rahmen hinsichtlich Konfiguration und Größenskala gesondert zu betrachten. Abschließend hat diese Studie demonstriert, dass die Betrachtung von Störungen in 3D Finiten Elementen Modellen für die Simulation von gekoppeltem Fluid- und Wärmetransport auf unterschiedlichen Skalen möglich ist. Da diese Art von numerischen Simulationen sowohl die geologische Struktur des Untergrunds sowie die im Erdinnern ablaufenden physikalischen Prozesse integriert, können sie einen wertvollen Beitrag leisten, indem sie Feld- und Laborgestützte Untersuchungen vervollständigen. / Faults can act as conduits, barriers or mixed conduit/barrier systems to fluid flow. Therefore, faults may significantly influence fluid flow regimes operating in the subsurface, possibly resulting in distinct variations of the deep thermal field. Both, flow dynamics and temperature changes are in turn crucial factors that need to be taken into account for geothermal energy exploration. This study investigated the influence of faults on the subsurface fluid system and thermal field and explored the processes controlling fluid behavior and thermal distribution both within host rocks and faults. For this purpose, 3D finite element simulations of coupled fluid and heat transport have been carried out, both for synthetic and real-case model scenarios on different scales. A small-scale synthetic model was developed to systematically assess the impact of an inclined fault by changing gradually its hydraulic width and its permeability within the simulations. An observed linear inverse relationship revealed that changing the fault width by one order of magnitude results in a fluid velocity decrease (~1e-01 m/s) within the fault. A high permeability contrast between fault and matrix favors fluid advection into the fault and leads to pronounced pressure and temperature changes in and around the same domain. When the permeability contrast between fault domain and host rock is low, however, no fluid flow is observed in the fault, thus resulting in undisturbed hydrostatic pressure and temperature fields. On the basis of these synthetic fault modelling results, the influence of faults on a larger scale have been analyzed within a more complex (real-case) geological setting,- a 3D model of the geothermal site Groß Schönebeck / located ~40 km north of Berlin. The integration of one permeable and three impermeable major faults, resulted in distinct changes observed in the local fluid circulation, thermal and pressure field. Modelled convective circulation within the permeable fault decisively modifies the thermal field (up to 15 K). Within the low permeable faults, heat is transferred only by conduction, inducing no thermal imprint but local deviations of the hydrostatic pressure field. To investigate the impact of major fault zones on the basin-scale geothermal field, coupled fluid and heat transport simulations have been conducted for a 3D structural model for Brandenburg region (Noack et al. 2010; 2013). Different geological scenarios in terms of modelled fault permeability have been carried out of which two end member models are analyzed. The results showed that tight fault zones affect the flow field locally. Acting as hydraulic barriers, fluid flow is deviated with very low velocities along them within a range of ~ 1 km on either sides. The modelled local changes in the groundwater circulation system have no considerable effect on the temperature field. By contrast, permeable fault zones induce a pronounced signature on the thermal field extending over a distance of ~ 2.4-8.8 km at -1000 m depth and ~6-12 km at -3000 m depth. This thermal signature, characterized by alternating cooler and hotter temperature domains, is controlled by up- and downward directed flow within the fault domain, principally driven by existing hydraulic head gradients. All studies demonstrated that faults have a considerable impact on the fluid and heat flow. The permeability in faults and surrounding geological layers as well as the specific geological setting turned out to be crucial factors in controlling the different kinds of heat transfer mechanisms that may evolve in faults. Temperature variations caused by permeable faults may be local but significant as well as changes in fluid dynamics by both conduits and barriers. Thus, the results demonstrated the importance to consider faults in geothermal energy exploration. In the final analysis, the simulations for the small-, regional- and basin-scale models showed that the outcomes cannot be transferred by upscaling and that it is necessary to consider each geological setting separately with respect to its configuration and scale dimension. In summary, this study demonstrated that the consideration of faults in 3D finite element models for coupled fluid and heat transport simulations on different scales is feasible. As these type of numerical simulations integrate both, the structural setting of the subsurface and the physical processes controlling subsurface transport, the outcomes of this thesis may provide positive contributions in that they valuably complement field- and laboratory-based investigations.

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.

Forecasting time series data is crucial in a number of domains such as supply chain management and display advertisement. In these areas, the time series data to forecast is typically organized along multiple dimensions leading to a high number of time series that need to be forecasted. Most current approaches focus only on selection and optimizing a forecast model for a single time series. In this paper, we explore how we can utilize time series at different dimensions to increase forecast accuracy and, optionally, reduce model maintenance overhead. Solving this problem is challenging due to the large space of possibilities and possible high model creation costs. We propose a model configuration advisor that automatically determines the best set of models, a model configuration, for a given multi-dimensional data set. Our approach is based on a general process that iteratively examines more and more models and simultaneously controls the search space depending on the data set, model type and available hardware. The final model configuration is integrated into F2DB, an extension of PostgreSQL, that processes forecast queries and maintains the configuration as new data arrives. We comprehensively evaluated our approach on real and synthetic data sets. The evaluation shows that our approach significantly increases forecast query accuracy while ensuring low model costs.