TY - JOUR A1 - Landgraf, Angela A1 - Zielke, Olaf A1 - Arrowsmith, J. Ramón A1 - Ballato, Paolo A1 - Strecker, Manfred A1 - Schildgen, Taylor F. A1 - Friedrich, Anke M. A1 - Tabatabaei, Sayyed-Hassan T1 - Differentiating simple and composite tectonic landscapes using numerical fault slip modeling with an example from the south central Alborz Mountains, Iran JF - Journal of geophysical research : Earth surface N2 - The tectonically driven growth of mountains reflects the characteristics of the underlying fault systems and the applied tectonic forces. Over time, fault networks might be relatively static, but stress conditions could change and result in variations in fault slip orientation. Such a tectonic landscape would transition from a simple to a composite state: the topography of simple landscapes is correlated with a single set of tectonic boundary conditions, while composite landscapes contain inherited topography due to earlier deformation under different boundary conditions. We use fault interaction modeling to compare vertical displacement fields with topographic metrics to differentiate the two types of landscapes. By successively rotating the axis of maximum horizontal stress, we produce a suite of vertical displacement fields for comparison with real landscapes. We apply this model to a transpressional duplex in the south central Alborz Mountains of Iran, where NW oriented compression was superseded by neotectonic NE compression. The consistency between the modeled displacement field and real landforms indicates that the duplex topography is mostly compatible with the modern boundary conditions, but might include a small remnant from the earlier deformation phase. Our approach is applicable for various tectonic settings and represents an approach to identify the changing boundary conditions that produce composite landscapes. It may be particularly useful for identifying changes that occurred in regions where river profiles may no longer record a signal of the change or where the spatial pattern of uplift is complex. KW - fault interaction KW - landscape evolution KW - numerical modeling KW - Alborz Mountains KW - Iran Y1 - 2013 U6 - https://doi.org/10.1002/jgrf.20109 SN - 2169-9003 SN - 2169-9011 VL - 118 IS - 3 SP - 1792 EP - 1805 PB - American Geophysical Union CY - Washington ER - TY - THES A1 - Mulyukova, Elvira T1 - Stability of the large low shear velocity provinces T1 - Stabilität der basalen Melange im untersten Erdmantel BT - numerical modeling of thermochemical mantle convection BT - numerische Modellierung thermochemischer Mantelkonvektion N2 - We study segregation of the subducted oceanic crust (OC) at the core mantle boundary and its ability to accumulate and form large thermochemical piles (such as the seismically observed Large Low Shear Velocity Provinces - LLSVPs). Our high-resolution numerical simulations suggest that the longevity of LLSVPs for up to three billion years, and possibly longer, can be ensured by a balance in the rate of segregation of high-density OC-material to the CMB, and the rate of its entrainment away from the CMB by mantle upwellings. For a range of parameters tested in this study, a large-scale compositional anomaly forms at the CMB, similar in shape and size to the LLSVPs. Neutrally buoyant thermochemical piles formed by mechanical stirring - where thermally induced negative density anomaly is balanced by the presence of a fraction of dense anomalous material - best resemble the geometry of LLSVPs. Such neutrally buoyant piles tend to emerge and survive for at least 3Gyr in simulations with quite different parameters. We conclude that for a plausible range of values of density anomaly of OC material in the lower mantle - it is likely that it segregates to the CMB, gets mechanically mixed with the ambient material, and forms neutrally buoyant large scale compositional anomalies similar in shape to the LLSVPs. We have developed an efficient FEM code with dynamically adaptive time and space resolution, and marker-in-cell methodology. This enabled us to model thermochemical mantle convection at realistically high convective vigor, strong thermally induced viscosity variations, and long term evolution of compositional fields. N2 - Es wird allgemein akzeptiert, dass Mantelkonvektion - das langsame Fließen der Mantelgesteine, das mutmaßlich ein wichtiger Antrieb der Plattentektonik ist - von Dichteunterschieden verursacht wird, die thermischen aber auch chemischen Ursprungs sind. Es fehlen aber Kenntnisse über die thermochemischen Prozesse im Erdinneren, vor allem wegen Schwierigkeiten bei der Beobachtung. Eines der zuverlässigsten Resultate von tomographischen Beobachtungen ist die Existenz von zwei Haufen einer basalen Melange (BAM, LLSVP auf Englisch), die sich auf gegenüber liegenden Seiten in 3000 km Tiefe am Boden des Mantels unter Afrika bzw dem Pazifik befinden. Die niedrige Scherwellengeschwindigkeit in der BAM scheint eine thermischen (heiß) sowie einen chemischen (Material mit hoher Dichte) Ursprung zu haben. Aufgrund von plattentektonischen Rekonstruktionen wird angenommen dass die BAM langlebig und stabil sind, und dass sie von überwiegend von ihren Rändern hochquellenden Manteldiapiren beprobt werden. Die Hauptfrage meiner Doktorarbeit ist, wie solche großen chemischen Speicher wie die BAM sich bilden und über hunderte von Millionen Jahren überleben können, ohne dass sie von der Mantelkonvektion zerstört werden. Was sind die physikalischen Eigenschaften des BAM-Materials, z.B. Dichte, die dazu beitragen? Ich benutze numerische Modellierung um zu erforschen, wie sich eine dichte Bodenschicht bildet und wie die Mantelkonvektion Material daraus mitnimmt. Mein Ziel ist, die langfristige thermochemische Entwicklung des Erdmantels zu verstehen, insbesondere die Rolle der Dichteheterogeintäten Viskosität im untersten Mantel. KW - earth's mantle KW - thermochemical mantle convection KW - numerical modeling KW - Erdmantel KW - thermochemischer Mantelkonvektion KW - numerische Modellierung Y1 - 2015 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-82228 ER - TY - JOUR A1 - Brune, Sascha A1 - Corti, Giacomo A1 - Ranalli, Giorgio T1 - Controls of inherited lithospheric heterogeneity on rift linkage: Numerical and analog models of interaction between the Kenyan and Ethiopian rifts across the Turkana depression JF - Tectonics N2 - Inherited rheological structures in the lithosphere are expected to have large impact on the architecture of continental rifts. The Turkana depression in the East African Rift connects the Main Ethiopian Rift to the north with the Kenya rift in the south. This region is characterized by a NW-SE trending band of thinned crust inherited from a Mesozoic rifting event, which is cutting the present-day N-S rift trend at high angle. In striking contrast to the narrow rifts in Ethiopia and Kenya, extension in the Turkana region is accommodated in subparallel deformation domains that are laterally distributed over several hundred kilometers. We present both analog experiments and numerical models that reproduce the along-axis transition from narrow rifting in Ethiopia and Kenya to a distributed deformation within the Turkana depression. Similarly to natural observations, our models show that the Ethiopian and Kenyan rifts bend away from each other within the Turkana region, thus forming a right-lateral step over and avoiding a direct link to form a continuous N-S depression. The models reveal five potential types of rift linkage across the preexisting basin: three types where rifts bend away from the inherited structure connecting via a (1) wide or (2) narrow rift or by (3) forming a rotating microplate, (4) a type where rifts bend towards it, and (5) straight rift linkage. The fact that linkage type 1 is realized in the Turkana region provides new insights on the rheological configuration of the Mesozoic rift system at the onset of the recent rift episode. Plain Language Summary The Turkana depression in the Kenya/Ethiopia borderland is most famous for its several million years old human fossils, but it also holds a rich geological history of continental separation. The Turkana region is a lowland located between the East African and Ethiopian domes because its crust and mantle have been stretched in a continent-wide rift event during Cretaceous times about 140-120 Ma ago. This thin lithosphere exerted paramount control on the dynamics of East African rifting in this area, which commenced around 15 Ma ago and persists until today. Combining analog "sandbox" experiments with numerical geodynamic modeling, we find that inherited rift structures explain the dramatic change in rift style from deep, narrow rift basins north and south of the Turkana area to wide, distributed deformation within the Turkana depression. The failed Cretaceous rift is also responsible for the eastward bend of the Ethiopian rift and the westward bend of the Kenyan rift when entering the Turkana depression, which generated the characteristic hook-shaped form of present-day Lake Turkana. Combing two independent modeling techniques-analog and numerical experiments-is a very promising approach allowing to draw robust conclusions about the processes that shape the surface of our planet. KW - continental rifting KW - east African rift KW - numerical modeling KW - analog modeling KW - tectonic inheritance KW - Turkana depression Y1 - 2017 U6 - https://doi.org/10.1002/2017TC004739 SN - 0278-7407 SN - 1944-9194 VL - 36 SP - 1767 EP - 1786 PB - American Geophysical Union CY - Washington ER - TY - THES A1 - Muldashev, Iskander T1 - Modeling of the great earthquake seismic cycles T1 - Modellierung der seismischen Zyklen von Mega-Erdbeben N2 - The timing and location of the two largest earthquakes of the 21st century (Sumatra, 2004 and Tohoku 2011, events) greatly surprised the scientific community, indicating that the deformation processes that precede and follow great megathrust earthquakes remain enigmatic. During these phases before and after the earthquake a combination of multi-scale complex processes are acting simultaneously: Stresses built up by long-term tectonic motions are modified by sudden jerky deformations during earthquakes, before being restored by multiple ensuing relaxation processes. This thesis details a cross-scale thermomechanical model developed with the aim of simulating the entire subduction process from earthquake (1 minute) to million years’ time scale, excluding only rupture propagation. The model employs elasticity, non-linear transient viscous rheology, and rate-and-state friction. It generates spontaneous earthquake sequences, and, by using an adaptive time-step algorithm, recreates the deformation process as observed naturally over single and multiple seismic cycles. The model is thoroughly tested by comparing results to those from known high- resolution solutions of generic modeling setups widely used in modeling of rupture propagation. It is demonstrated, that while not modeling rupture propagation explicitly, the modeling procedure correctly recognizes the appearance of instability (earthquake) and correctly simulates the cumulative slip at a fault during great earthquake by means of a quasi-dynamic approximation. A set of 2D models is used to study the effects of non-linear transient rheology on the postseismic processes following great earthquakes. Our models predict that the viscosity in the mantle wedge drops by 3 to 4 orders of magnitude during a great earthquake with magnitude above 9. This drop in viscosity results in spatial scales and timings of the relaxation processes following the earthquakes that are significantly different to previous estimates. These models replicate centuries long seismic cycles exhibited by the greatest earthquakes (like the Great Chile 1960 Earthquake) and are consistent with the major features of postseismic surface displacements recorded after the Great Tohoku Earthquake. The 2D models are also applied to study key factors controlling maximum magnitudes of earthquakes in subduction zones. Even though methods of instrumentally observing earthquakes at subduction zones have rapidly improved in recent decades, the characteristic recurrence interval of giant earthquakes (Mw>8.5) is much larger than the currently available observational record and therefore the necessary conditions for giant earthquakes are not clear. Statistical studies have recognized the importance of the slab shape and its surface roughness, state of the strain of the upper plate and thickness of sediments filling the trenches. In this thesis we attempt to explain these observations and to identify key controlling parameters. We test a set of 2D models representing great earthquake seismic cycles at known subduction zones with various known geometries, megathrust friction coefficients, and convergence rates implemented. We found that low-angle subduction (large effect) and thick sediments in the subduction channel (smaller effect) are the fundamental necessary conditions for generating giant earthquakes, while the change of subduction velocity from 10 to 3.5 cm/yr has a lower effect. Modeling results also suggest that having thick sediments in the subduction channel causes low static friction, resulting in neutral or slightly compressive deformation in the overriding plate for low-angle subduction zones. These modeling results agree well with observations for the largest earthquakes. The model predicts the largest possible earthquakes for subduction zones of given dipping angles. The predicted maximum magnitudes exactly threshold magnitudes of all known giant earthquakes of 20th and 21st centuries. The clear limitation of most of the models developed in the thesis is their 2D nature. Development of 3D models with comparable resolution and complexity will require significant advances in numerical techniques. Nevertheless, we conducted a series of low-resolution 3D models to study the interaction between two large asperities at a subduction interface separated by an aseismic gap of varying width. The novelty of the model is that it considers behavior of the asperities during multiple seismic cycles. As expected, models show that an aseismic gap with a narrow width could not prevent rupture propagation from one asperity to another, and that rupture always crosses the entire model. When the gap becomes too wide, asperities do not interact anymore and rupture independently. However, an interesting mode of interaction was observed in the model with an intermediate width of the aseismic gap: In this model the asperities began to stably rupture in anti-phase following multiple seismic cycles. These 3D modeling results, while insightful, must be considered preliminary because of the limitations in resolution. The technique developed in this thesis for cross-scale modeling of seismic cycles can be used to study the effects of multiple seismic cycles on the long-term deformation of the upper plate. The technique can be also extended to the case of continental transform faults and for the advanced 3D modeling of specific subduction zones. This will require further development of numerical techniques and adaptation of the existing advanced highly scalable parallel codes like LAMEM and ASPECT. N2 - Zeitpunkt und Ort der zwei größten Erdbeben des 21. Jahrhunderts (Sumatra 2004 und Tohoku 2011) überraschten die wissenschaftliche Gemeinschaft, da sie darauf hindeuten, dass die einem Megathrust-Erdbeben vorangehenden und folgenden Deformationsprozesse weiterhin rästelhaft bleiben. Ein Problem ist die komplexe Art der Subduktionsprozesse auf unterschiedlichen Skalen. Spannungen werden durch langzeitliche, tektonische Bewegungen aufgebaut, von plötzlicher, ruckartiger Deformation während Erdbeben modifiziert und anschließend durch verschiedene Entspannungsprozesse wiederhergestellt. In dieser Arbeit wird ein skalen-übergreifendes thermomechanisches Modell entwickelt mit dem Ziel den vollständigen Subduktionsprozess von kleiner Skala (Minuten, z.B. Erdbeben) bis zu Millionen Jahren zu simulieren. Dabei bleibt nur das dynamische Prozess der unberücksichtigt. Das Modell nutzt Elastizität, nicht-lineare transient viskose Rheologie und „rate-and-state“ Reibungsgesetz. Es erzeugt spontane Erdbeben-Sequenzen und stellt durch das Nutzen eines adaptiven Zeitschritt Algorithmus den Deformationsprozess nach, wie er in der Natur während einzelner und mehrerer seismischer Zyklen beobachtet wird. Anhand der 2D-Modell Reihe werden die Effekte der nicht-linearen transient (viskosen) Rheologie auf postseismische Prozesse nach großen Erdbeben untersucht. Die Modelle sagen eine Verringerung der Mantelkeil-Viskosität um 3 bis 4 Größenordnungen während der großen Erdbeben (Mw > 9) vorher. Im Gegensatz zur momentanen Meinung, ergeben sich andere räumliche und zeitliche Verteilungen des Entspannungsprozesses nach großen Erdbeben. Jahrhunderte lange seismische Zyklen, wie zum Beispiel das große Erdbeben in Chile 1960, werden reproduziert und stimmen mit den Hauptmerkmalen der postseismischen Oberflächenverschiebungen überein, die nach dem großen Tohoku Erdbeben aufgezeichnet wurden. Mithilfe der 2D-Modelle werden außerdem Schlüsselfaktoren untersucht, die die maximale Magnitude von Erdbeben bestimmen. Obwohl die instrumentellen Methoden zur Beobachtung von Erdbeben in Subduktionszonen in den letzten Jahrzehnten stetig verbessert wurden, kann das Auftreten der größten Erdbeben (Mw > 8.5) und die dafür vorrauszusetzenden Bedingungen nicht genau bestimmt werden, da die charackeristische Zeit für ein Wiederauftreten deutlich größer als der Beobachtungszeitraum ist. Wir versuchen in dieser Arbeit diese Beobachtungen zu erklären und die kontrollierenden Schlüsselfaktoren zu bestimmen. Eine 2D-Modell Reihe mit großen, seismischen Erdbebenzyklen untersucht unterschiedliche Geometrien von Subduktionszonen sowie Reibungskoeffizieten an der Subduktions Platten-Grenzoberfläche und Konvergenzraten. Subduktion mit geringem Winkel (großer Effekt) und mächtige Sedimente im Subduktionskanal (kleiner Effekt) sind grundlegend notwendige Bedingungen für riesige Erdbeben, während eine Änderung des Reibungsparameters nur geringe Auswirkungen hat. Die Modellierungsergebnisse sind in viiÜbereinstimmung mit den Beobachtungen von den größten Erdbeben. Die maximale vohergesagte Magnituden liegen an der oberen Grenze alle historische Beben der 20. und 21. Jahrhunderden. Die größte Einschränkung der entwickelten Modelle ist ihre 2-dimensionale Natur. Um die Interaktion zwischen zwei großen Unebenheiten („asperities“) auf der Subduktions Platten-Grenzoberfläche, die von einer aseismischen Lücke („gap“) voneinander getrennt sind, zu untersuchen, wurde eine 3D-Modell Reihe mit geringer Auflösung durchgeführt. Neu an diesem Modell ist, dass das Verhalten der Unebenheiten während mehrerer seismischer Zyklen berücksichtigt wird. Wie erwartet zeigt das Modell, dass eine schmale, aseismische Lücke ein Propagieren von Brüchen von einer Unebenheit zur nächsten nicht verhindern kann und diese Brüche das gesamte Modell durchkreuzen. Wird die Lücke zu groß interagieren die Unebenheiten nicht mehr miteinander und brechen unabhängig voneinander. Allerdings wurde in dem Modell mit mittlerer Breite der aseismischen Lücke eine interessante Art von Interaktion beobachtet: Nach einigen seismischen Zyklen beginnen die Unebenheiten kontinuerlich und zwar in Gegen-Phasen zu brechen. Diese Ergebnisse sind andeutend, aber können aufgrund der geringen Auflösung des Models nur als vorläufig angesehen werden und erfordern weitere Bestätigung mit hoch-auflösenden Modellen. KW - earthquake modeling KW - seismic cycle modeling KW - subduction KW - numerical modeling KW - Erdbeben Modellierung KW - Modellierung des seismischen Zyklus KW - Subduktion KW - numerische Modellierung Y1 - 2017 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-398926 ER - TY - THES A1 - Munz, Matthias T1 - Water flow and heat transport modelling at the interface between river and aquifer T1 - Wasserfluss- und Wärmetransportmodellierung an der Schnittstelle zwischen Fluss und Grundwasserleiter N2 - The functioning of the surface water-groundwater interface as buffer, filter and reactive zone is important for water quality, ecological health and resilience of streams and riparian ecosystems. Solute and heat exchange across this interface is driven by the advection of water. Characterizing the flow conditions in the streambed is challenging as flow patterns are often complex and multidimensional, driven by surface hydraulic gradients and groundwater discharge. This thesis presents the results of an integrated approach of studies, ranging from the acquisition of field data, the development of analytical and numerical approaches to analyse vertical temperature profiles to the detailed, fully-integrated 3D numerical modelling of water and heat flux at the reach scale. All techniques were applied in order to characterize exchange flux between stream and groundwater, hyporheic flow paths and temperature patterns. The study was conducted at a reach-scale section of the lowland Selke River, characterized by distinctive pool riffle sequences and fluvial islands and gravel bars. Continuous time series of hydraulic heads and temperatures were measured at different depths in the river bank, the hyporheic zone and within the river. The analyses of the measured diurnal temperature variation in riverbed sediments provided detailed information about the exchange flux between river and groundwater. Beyond the one-dimensional vertical water flow in the riverbed sediment, hyporheic and parafluvial flow patterns were identified. Subsurface flow direction and magnitude around fluvial islands and gravel bars at the study site strongly depended on the position around the geomorphological structures and on the river stage. Horizontal water flux in the streambed substantially impacted temperature patterns in the streambed. At locations with substantial horizontal fluxes the penetration depths of daily temperature fluctuations was reduced in comparison to purely vertical exchange conditions. The calibrated and validated 3D fully-integrated model of reach-scale water and heat fluxes across the river-groundwater interface was able to accurately represent the real system. The magnitude and variations of the simulated temperatures matched the observed ones, with an average mean absolute error of 0.7 °C and an average Nash Sutcliffe Efficiency of 0.87. The simulation results showed that the water and heat exchange at the surface water-groundwater interface is highly variable in space and time with zones of daily temperature oscillations penetrating deep into the sediment and spots of daily constant temperature following the average groundwater temperature. The average hyporheic flow path temperature was found to strongly correlate with the flow path residence time (flow path length) and the temperature gradient between river and groundwater. Despite the complexity of these processes, the simulation results allowed the derivation of a general empirical relationship between the hyporheic residence times and temperature patterns. The presented results improve our understanding of the complex spatial and temporal dynamics of water flux and thermal processes within the shallow streambed. Understanding these links provides a general basis from which to assess hyporheic temperature conditions in river reaches. N2 - Die Interaktion zwischen Oberflächenwasser und Grundwasser hat einen entscheidenden Einfluss auf die Wasserqualität und die ökologische Beschaffenheit von Seen, Flüssen und aquatischen Ökosystemen. Der Austausch von Wärme und gelösten Substanzen zwischen diesen beiden Kompartimenten ist maßgeblich durch die Austauschraten und die Strömungsrichtung des Wassers beeinflusst. Somit ist die Charakterisierung dieser beiden Größen in dem Übergangsbereich zwischen Oberflächenwasser und Grundwasser von besonderer Bedeutung. Diese Arbeit präsentiert die Entwicklung und Anwendung von Methoden zur Untersuchung der zeitlichen und räumlichen Dynamik des Wasser- und Wärmeflusses an der Schnittstelle zwischen Oberflächenwasser und Grundwasser. Die Arbeit besteht im Wesentlichen aus zwei Schwerpunkten. Der erste Schwerpunkt beinhaltet die Entwicklung und Bewertung von analytischen und numerischen Methoden zur Bestimmung der horizontalen Strömungsrichtung und Austauschraten unter Verwendung von kontinuierlich gemessenen Temperaturzeitreihen entlang vertikaler Profile im gesättigten Sediment. Flussbetttemperaturen können relativ einfach und kostengünstig entlang eines Flussabschnittes in verschiedenen Tiefen und unterschiedlichsten Flussbettsedimenten (organisch, sandig bis grob kiesig) gemessen werden. Die Hauptverwendung solcher Temperaturprofile ist bisher auf die analytische Quantifizierung vertikaler Austauschraten limitiert. Im Rahmen dieser Arbeit wurde ein Ansatz entwickelt, der eine punktuelle Ermittlung der horizontalen Strömungs-komponente, basierend auf der Veränderung der täglichen Temperaturamplitude mit zunehmender Tiefe, ermöglicht. Weiterhin wurde ein inverser, numerischer Ansatz entwickelt, der die ein-dimensionale Wärmetransportgleichung numerisch löst und mittels inverser Optimierungsalgorithmen die simulierten Temperaturen bestmöglich an die gemessenen Flussbetttemperaturen anpasst. Diese Methode ermöglicht die automatische, zeitlich variable Quantifizierung vertikaler Austauschraten an der Schnittstelle zwischen Oberflächenwasser und Grundwasser sowie eine einfache Unsicherheitsbetrachtung aufgrund der zugrunde liegenden Parameterunsicherheiten. Der zweite Schwerpunkt der Arbeit liegt auf der Anwendung der entwickelten Methoden und der Erstellung eines dreidimensionalen Simulationsmodelles entlang eines Flussabschnittes der Selke. Numerische Strömungs- und Stofftransportmodelle ermöglichen die gekoppelte Simulation von Fließprozessen im offenen Gerinne und im darunter liegenden porösen Medium. Die Parametrisierung des Modells erfolgte anhand empirischer Daten die im Untersuchungsgebiet detailliert erhoben wurden. Die Simulationsergebnisse zeigten zum einen gebietsspezifische Gegebenheiten auf, ermöglichten darüber hinaus jedoch auch die Beschreibung allgemeiner Muster und Zusammenhänge welche die Wasserfluss- und Wärmetransportprozesse an der Schnittstelle zwischen Oberflächenwasser und Grundwasser beeinflussen. So zeigten die Ergebnisse dieser Arbeit, dass maßgeblich die natürlich vorhandenen Flussbettstrukturen die Austauschraten und die Strömungsrichtung zwischen Oberflächenwasser und Grundwasser beeinflussen und somit den hyporheischen Austausch bestimmen. Wichtige Einflussfaktoren auf die untersuchten Austauschprozesse waren die Lage im Gerinne relativ zur Flussbettstruktur und der vorherrschende Wasserstand (Abfluss). Bedingt durch den Wasser- und Wärmeaustausch prägten sich im Untersuchungsgebiet Bereiche aus in denen die täglichen Temperaturschwingungen tief in das Sediment eindringen (Anstrombereich der Flussbettstrukturen), als auch Bereiche in denen relativ konstante Temperaturen, nahe der Grundwassertemperatur, vorherrschten. Die durchschnittliche Temperatur in der hyporheischen Zone wurde durch die saisonalen Temperaturschwankungen im Oberflächenwasser dominiert, wobei die Temperaturen entlang einzelner Fließpfade stark von der Verweilzeit des Oberflächen- oder Grundwassers im gesättigten Sediment und dem Temperaturgradienten zwischen Fluss und Grundwasser abhängig waren. Trotz der Komplexität dieser Zusammenhänge, ermöglichten die Simulationsergebnisse die Ableitung einer allgemeinen empirischen Beziehung zwischen den hyporheischen Verweilzeiten und Temperaturmustern. Sowohl die Verweilzeiten als auch die Temperatur im gesättigten Sediment haben einen entscheiden Einfluss auf biogeochemische Prozesse in dem Übergangsbereich zwischen Oberflächenwasser und Grundwasser und sind somit von besonderer Bedeutung für die Wasserqualität von Seen, Flüssen und aquatischen Ökosystemen. KW - surface water-groundwater interaction KW - Oberflächenwasser-Grundwasser Interaktion KW - hyporheic zone KW - hyporheische Zone KW - numerical modeling KW - numerische Modellierung KW - heat transport KW - Wärmetransport KW - temperature KW - Temperatur Y1 - 2017 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-404319 ER - TY - THES A1 - Korges, Maximilian T1 - Constraining the hydrology of intrusion-related ore deposits with fluid inclusions and numerical modeling T1 - Bestimmung der Hydrologie von Erzlagerstätten an Intrusionen mit Flüssigkeitseinschlüssen und numerischer Modellierung N2 - Magmatic-hydrothermal fluids are responsible for numerous mineralization types, including porphyry copper and granite related tin-tungsten (Sn-W) deposits. Ore formation is dependent on various factors, including, the pressure and temperature regime of the intrusions, the chemical composition of the magma and hydrothermal fluids, and fluid rock interaction during the ascent. Fluid inclusions have potential to provide direct information on the temperature, salinity, pressure and chemical composition of fluids responsible for ore formation. Numerical modeling allows the parametrization of pluton features that cannot be analyzed directly via geological observations. Microthermometry of fluid inclusions from the Zinnwald Sn-W deposit, Erzgebirge, Germany / Czech Republic, provide evidence that the greisen mineralization is associated with a low salinity (2-10 wt.% NaCl eq.) fluid with homogenization temperatures between 350°C and 400°C. Quartzes from numerous veins are host to inclusions with the same temperatures and salinities, whereas cassiterite- and wolframite-hosted assemblages with slightly lower temperatures (around 350°C) and higher salinities (ca. 15 wt. NaCl eq.). Further, rare quartz samples contained boiling assemblages consisting of coexisting brine and vapor phases. The formation of ore minerals within the greisen is driven by invasive fluid-rock interaction, resulting in the loss of complexing agents (Cl-) leading to precipitation of cassiterite. The fluid inclusion record in the veins suggests boiling as the main reason for cassiterite and wolframite mineralization. Ore and coexisting gangue minerals hosted different types of fluid inclusions where the beginning boiling processes are solely preserved by the ore minerals emphasizing the importance of microthermometry in ore minerals. Further, the study indicates that boiling as a precipitation mechanism can only occur in mineralization related to shallow intrusions whereas deeper plutons prevent the fluid from boiling and can therefore form tungsten mineralization in the distal regions. The tin mineralization in the Hämmerlein deposit, Erzgebirge, Germany, occurs within a skarn horizon and the underlying schist. Cassiterite within the skarn contains highly saline (30-50 wt% NaCl eq.) fluid inclusions, with homogenization temperatures up to 500°C, whereas cassiterites from the schist and additional greisen samples contain inclusions of lower salinity (~5 wt% NaCl eq.) and temperature (between 350 and 400°C). Inclusions in the gangue minerals (quartz, fluorite) preserve homogenization temperatures below 350°C and sphalerite showed the lowest homogenization temperatures (ca. 200°C) whereby all minerals (cassiterite from schist and greisen, gangue minerals and sphalerite) show similar salinity ranges (2-5 wt% NaCl eq.). Similar trace element contents and linear trends in the chemistry of the inclusions suggest a common source fluid. The inclusion record in the Hämmerlein deposit documents an early exsolution of hot brines from the underlying granite which is responsible for the mineralization hosted by the skarn. Cassiterites in schist and greisen are mainly forming due to fluid-rock interaction at lower temperatures. The low temperature inclusions documented in the sphalerite mineralization as well as their generally low trace element composition in comparison to the other minerals suggests that their formation was induced by mixing with meteoric fluids. Numerical simulations of magma chambers and overlying copper distribution document the importance of incremental growth by sills. We analyzed the cooling behavior at variable injection intervals as well as sill thicknesses. The models suggest that magma accumulation requires volumetric injection rates of at least 4 x 10-4 km³/y. These injection rates are further needed to form a stable magmatic-hydrothermal fluid plume above the magma chamber to ensure a constant copper precipitation and enrichment within a confined location in order to form high-grade ore shells within a narrow geological timeframe between 50 and 100 kyrs as suggested for porphyry copper deposits. The highest copper enrichment can be found in regions with steep temperature gradients, typical of regions where the magmatic-hydrothermal fluid meets the cooler ambient fluids. N2 - Magmatisch-hydrothermale Fluide sind verantwortlich für zahlreiche Mineralisationstypen, wie porphyrische Kupferlagerstätten und granitgebundene Zinn-Wolfram (Sn-W) Lagerstätten. Die Lagerstättenbildung ist abhängig von unterschiedlichen Faktoren, z.B. dem Druck- und Temperaturregime der Intrusion, der chemischen Zusammensetzung des Magmas und der hydrothermalen Fluide sowie den Fluid-Gesteinsreaktionen während des Aufstiegs der Fluide. Flüssigkeitseinschlüsse haben das Potential, direkte Informationen zur Temperatur, zum Salzgehalt, zum Druck und der Chemie der Fluide, welche für die Lagerstättenbildung verantwortlich sind, zu liefern. Außerdem erlauben numerische Modellierungen die Parametrisierung der Plutoneigenschaften, die nicht direkt anhand von geologischen Beobachtungen analysiert werden können. Mikrothermometrie von Flüssigkeitseinschlüssen der Zinnwald Sn-W Lagerstätte zeigen, dass die Vergreisung an ein schwach salzhaltiges (2-10 wt.% NaCl eq.) Fluid gebunden ist, das zwischen 350°C und 400°C homogenisiert. Quarze der diversen Gänge beinhalten Einschlüsse mit den gleichen Temperaturen und Salzgehalten, wohingegen Kassiterit und Wolframit Einschlüsse mit niedrigeren Temperaturen (um 350°C) und höheren Salzgehalten zeigen (ca. 15 wt. NaCl eq.). Seltene Quarzproben enthalten kochende Einschluss-Ansammlungen, die aus koexistierenden salzreichen Lösungen und gasreichen Phasen bestehen. Die Bildung der Erzminerale des Greisens entsteht durch tiefgreifende Fluid-Gesteinsreaktionen, die den Verlust des Komplexbildners (Cl-) zur Folge haben, wodurch Kassiterit ausgefällt wird. Die Einschlüsse in den Gängen verdeutlichen, dass kochende Fluide der Hauptgrund für die Kassiterit– und Wolframit–Mineralisation sind. Erz- und benachbarte Gangminerale beinhalten unterschiedliche Einschlusstypen, wobei die beginnende Phasenseparation ausschließlich in den Erzmineralen erhalten ist, wodurch die Bedeutung der Mikrothermometrie in Erzmineralen hervorgehoben wird. Die Studie verdeutlicht weiterhin, dass Phasenseparation als Ausfällungsmechanismus nur in Lagerstätten gefunden werden können, die an flache Intrusionen gebunden sind, wohingegen tiefsitzende Granite die Phasenseparation verhindern. Dies hat zur Folge, dass eine Wolfram–Vererzung eher distal zur Intrusion auftritt. Die Zinn-Mineralisation der Hämmerlein Lagerstätte tritt sowohl in einem Skarn–Horizont als auch im darunterliegenden Schiefer auf. Fluideinschlüsse in Kassiteriten des Skarns enthalten Fluide mit hohem Salzgehhalt (30-50 wt% NaCl eq.) und Homogenisierungstemperaturen von bis zu 500°C, wohingegen Kassiterite des Schiefers (und von zusätzlichen Greisenproben) Einschlüsse mit geringerem Salzgehalt (~5 wt% NaCl eq.) und geringeren Temperaturen (zwischen 350 und 400°C) enthalten. Einschlüsse in Gangmineralen (Quarz, Fluorit) zeigen Homogenisierungstemperaturen von unter 350°C und Einschlüsse in Sphaleriten konservieren die niedrigsten gemessenen Temperaturen (ca. 200°C). Allerdings haben Flüssigkeitseinschlüsse in allen Mineralen (Kassiterite der Schiefer und Greisen, Gangminerale und Sphalerit) einen annähernd gleichen Salzgehalt (2-5 wt% NaCl eq.). Ähnliche Spurenelementgehalte und lineare Trends in der Chemie der Einschlüsse deuten auf ein gemeinsames Ursprungsfluid hin. Die Einschlüsse in der Hämmerlein-Lagerstätte dokumentieren eine frühe Entmischung von heißen Fluiden mit hohem Salzgehalt aus dem darunterliegenden Granit, die für die Mineralisation im Skarn verantwortlich ist. Die Kassiterite der Schiefer und der Greisen formen sich hauptsächlich durch Fluid-Gesteinsreaktionen bei niedrigeren Temperaturen. Die Niedrigtemperatureinschlüsse in Sphalerit und die im Vergleich zu den anderen Mineralen niedrigen Spurenelementgehalte deuten auf eine Bildung aufgrund von Mischungsprozessen mit meteorischen Fluiden hin. Numerische Simulationen von Magmenkammern und deren darüber gelegenen Kupferverteilungen dokumentieren die Wichtigkeit des schrittweisen Wachstums einer Intrusion durch Sills. Wir untersuchten das Abkühlungsverhalten bei unterschiedlichen Injektionsintervallen sowie bei unterschiedlicher Mächtigkeit des Sills. Die Modelle deuten darauf hin, dass für die Bildung der Magmakammer eine Injektionsrate von mindestens ca. 4 x 10-4 km³/y benötigt wird. Solche Raten sind ebenfalls nötig um eine kontinuierliche Bildung von magmatisch-hydrothermalen Fluiden zu garantieren, denn nur dann können sich hoch vererzte Bereiche in einem kurzen geologischen Zeitraum von 50.000 bis 100.000 Jahren bilden, so wie es für porphyrische Kupferlagerstätten angenommen wird. Die höchsten Kupfergehalte bilden sich in Regionen mit steilem Temperaturgradient, also vor allem in Bereichen, wo das magmatisch-hydrothermale Fluid auf kältere meteorische Fluide trifft. KW - fluid inclusions KW - numerical modeling KW - Flüssigkeitseinschlüsse KW - numerische Modellierung Y1 - 2019 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-434843 ER - TY - JOUR A1 - Baes, Marzieh A1 - Sobolev, Stephan V. A1 - Gerya, Taras V. A1 - Brune, Sascha T1 - Subduction initiation by Plume-Plateau interaction BT - insights from numerical models JF - Geochemistry, geophysics, geosystems N2 - It has recently been demonstrated that the interaction of a mantle plume with sufficiently old oceanic lithosphere can initiate subduction. However, the existence of large lithospheric heterogeneities, such as a buoyant plateau, in proximity to a rising plume head may potentially hinder the formation of a new subduction zone. Here, we investigate this scenario by means of 3-D numerical thermomechanical modeling. We explore how plume-lithosphere interaction is affected by lithospheric age, relative location of plume head and plateau border, and the strength of the oceanic crust. Our numerical experiments suggest four different geodynamic regimes: (a) oceanic trench formation, (b) circular oceanic-plateau trench formation, (c) plateau trench formation, and (d) no trench formation. We show that regardless of the age and crustal strength of the oceanic lithosphere, subduction can initiate when the plume head is either below the plateau border or at a distance less than the plume radius from the plateau edge. Crustal heterogeneity facilitates subduction initiation of old oceanic lithosphere. High crustal strength hampers the formation of a new subduction zone when the plume head is located below a young lithosphere containing a thick and strong plateau. We suggest that plume-plateau interaction in the western margin of the Caribbean could have resulted in subduction initiation when the plume head impinged onto the oceanic lithosphere close to the border between plateau and oceanic crust. KW - subduction zone KW - plume KW - plateau KW - numerical modeling KW - plume-induced KW - subduction initiation (PISI) Y1 - 2020 U6 - https://doi.org/10.1029/2020GC009119 SN - 1525-2027 VL - 21 IS - 8 PB - American Geophysical Union CY - Washington ER - TY - JOUR A1 - Korges, Maximilian A1 - Weis, Philipp A1 - Andersen, Christine T1 - The role of incremental magma chamber growth on ore formation in porphyry copper systems JF - Earth & planetary science letters N2 - Porphyry copper deposits are formed by fluids released from felsic magmatic intrusions of batholithic dimensions, which are inferred to have been incrementally built up by a series of sill injections. The growth of the magma chamber is primarily controlled by the volumetric injection rate from deeper-seated magma reservoirs, but can further be influenced by hydrothermal convection and fluid release. To quantify the interplay between magma chamber growth, volatile expulsion and hydrothermal fluid flow during ore formation, we used numerical simulations that can model episodic sill injections in concert with multi-phase fluid flow. To build up a magma chamber that constantly maintains a small region of melt within a period of about 50 kyrs, an injection rate of at least 1.3 x 10(-3) km(3)/y is required. Higher magma influxes of 1.9 to 7.6 x 10(-3) km(3)/y are able to form magma chambers with a thickness of 2 to 3 km. Such an intrusion continuously produces magmatic volatiles which can precipitate a copper ore shell in the host rock about 2 km above the fluid injection location. The steady fluid flux from such an incrementally growing magma chamber maintains a stable magmatic fluid plume, precipitating a copper ore shell in a more confined region and resulting in higher ore grades than the ones generated by an instantaneous emplacement of a voluminous magma chamber. Our simulation results suggest that magma chambers related to porphyry copper deposits form by rapid and episodic injection of magma. Slower magma chamber growth rates more likely result in barren plutonic rocks, although they are geochemically similar to porphyry-hosting plutons. However, these low-frequency sill injection events without a significant magma chamber growth can generate magmatic fluid pulses that can reach the shallow subsurface and are typical for high-sulfidation epithermal deposits. KW - pluton KW - magma chamber KW - porphyry copper deposits KW - magmatic sill KW - numerical modeling KW - ore deposit Y1 - 2020 U6 - https://doi.org/10.1016/j.epsl.2020.116584 SN - 0012-821X SN - 1385-013X VL - 552 PB - Elsevier CY - Amsterdam [u.a.] ER - TY - THES A1 - Liu, Sibiao T1 - Controls of foreland-deformation patterns in the orogen-foreland shortening system N2 - The Andean Plateau (Altiplano-Puna Plateau) of the southern Central Andes is the second-highest orogenic plateau on our planet after Tibet. The Andean Plateau and its foreland exhibit a pronounced segmentation from north to south regarding the style and magnitude of deformation. In the Altiplano (northern segment), more than 300 km of tectonic shortening has been recorded, which started during the Eocene. A well-developed thin-skinned thrust wedge located at the eastern flank of the plateau (Subandes) indicates a simple-shear shortening mode. In contrast, the Puna (southern segment) records approximately half of the shortening of the Altiplano - and the shortening started later. The tectonic style in the Puna foreland switches to a thick-skinned mode, which is related to pure-shear shortening. In this study, carried out in the framework of the StRATEGy project, high-resolution 2D thermomechanical models were developed to systematically investigate controls of deformation patterns in the orogen-foreland pair. The 2D and 3D models were subsequently applied to study the evolution of foreland deformation and surface topography in the Altiplano-Puna Plateau. The models demonstrate that three principal factors control the foreland-deformation patterns: (i) strength differences in the upper lithosphere between the orogen and its foreland, rather than a strength difference in the entire lithosphere; (ii) gravitational potential energy of the orogen (GPE) controlled by crustal and lithospheric thicknesses, and (iii) the strength and thickness of foreland-basin sediments. The high-resolution 2D models are constrained by observations and successfully reproduce deformation structures and surface topography of different segments of the Altiplano-Puna plateau and its foreland. The developed 3D models confirm these results and suggest that a relatively high shortening rate in the Altiplano foreland (Subandean foreland fold-and-thrust belt) is due to simple-shear shortening facilitated by thick and mechanically weak sediments, a process which requires a much lower driving force than the pure-shear shortening deformation mode in the adjacent broken foreland of the Puna, where these thick sedimentary basin fills are absent. Lower shortening rate in the Puna foreland is likely accommodated in the forearc by the slab retreat. N2 - Das Andenplateau (Altiplano-Puna-Plateau) in den südlichen Zentralanden ist nach Tibet das zweithöchste orogene Plateau auf unserem Planeten. Dieses Plateau und sein Vorland weisen eine ausgeprägte Segmentierung von Nord nach Süd hinsichtlich Art und Ausmaß der Verformung auf. Im Altiplano (nördliches Segment) wird seit der im Eozän stattfindenden Deformation mehr als 300 km tektonische Verkürzung dokumentiert. Ein gut entwickelter sedimentärer Schubkeil bzw. Vorland-Überschiebungsgürtel (Subandin) an der Ostflanke des Plateaus (thin-skinned foreland deformation) deutet in dieser Region des Vorlandes auf Prozesse einfacher Scherung hin (simple-shear modus). Im Gegensatz dazu weist die Puna (südliches Plateausegment) ungefähr die Hälfte der Verkürzung des Altiplano auf - und die Verkürzung begann später. Außerdem geht der tektonische Stil im Puna-Vorland zu einem zerbrochenen Vorland mit Kristallinblöcken (thick-skinned foreland) über, der mit der Verkürzung durch reine Scherung (pure-shear modus) erklärt werden kann. In dieser Studie, die im Rahmen des StRATEGy-Projekts durchgeführt wurde, wurden hochauflösende thermomechanische 2D-Modelle entwickelt, um systematisch die Kontrolle von Verformungsmustern im Orogen-Vorland-Paar zu untersuchen. Die 2D- und 3D-Modelle wurden anschließend angewendet, um die Entwicklung der Vorlanddeformation und der Oberflächentopographie im Altiplano-Puna-Plateau zu verstehen. Die Modelle zeigen, dass drei Hauptfaktoren die Deformationsmuster des Vorlandes steuern: (i) Festigkeitsunterschiede in der oberen Lithosphäre zwischen dem Orogen und seinem Vorland - und nicht Festigkeitsunterschiede in der gesamten Lithosphäre; (ii) die gravitationsbezogene potentielle Energie des Orogens (GPE), die durch die Krusten- und Lithosphärenmächtigkeit gesteuert wird und (iii) die Festigkeit sowie Mächtigkeiten der Vorlandbeckensedimente. Die hochauflösenden 2D-Modelle sind auf tatsächliche Daten aus Beobachtungen beschränkt und reproduzieren erfolgreich Deformationsstrukturen sowie die topographischen Verhältnisse der verschiedenen Segmente des Altiplano-Puna-Plateaus und seines Vorlandes. Die entwickelten 3D-Modelle bestätigen diese Ergebnisse und legen nahe, dass die relativ hohe Verkürzungsrate im Altiplano-Vorland (Subandin) bei den vorhandenen mächtigen Sedimentabfolgen geringer mechanischer Festigkeit weniger Kraftaufwand erfordert als die Deformation des Puna-Vorlandes, wo diese Sedimente weitgehend fehlen. Die geringeren Verkürzungsbeträge im Puna-Vorland werden wahrscheinlich durch das Zurückweichen der Subduktionszone im Forearc-Bereich ausgeglichen. KW - geodynamics KW - numerical modeling KW - Central Andes KW - foreland deformation KW - geophysics KW - Geodynamik KW - numerische Modellierung KW - Zentralanden KW - Vorlanddeformation KW - Geophysik Y1 - 2020 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-445730 ER - TY - GEN A1 - Bilbao-Lasa, Peru A1 - Jara-Muñoz, Julius A1 - Pedoja, Kevin A1 - Álvarez, Irantzu A1 - Aranburu, Arantza A1 - Iriarte, Eneko A1 - Galparsoro, Ibon T1 - Submerged marine terraces identification and an approach for numerical modeling the sequence formation in the Bay of Biscay (Northeastern Iberian Peninsula) T2 - Zweitveröffentlichungen der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe N2 - Submerged sequences of marine terraces potentially provide crucial information of past sea-level positions. However, the distribution and characteristics of drowned marine terrace sequences are poorly known at a global scale. Using bathymetric data and novel mapping and modeling techniques, we studied a submerged sequence of marine terraces in the Bay of Biscay with the objective to identify the distribution and morphologies of submerged marine terraces and the timing and conditions that allowed their formation and preservation. To accomplish the objectives a high-resolution bathymetry (5 m) was analyzed using Geographic Information Systems and TerraceM(R). The successive submerged terraces were identified using a Surface Classification Model, which linearly combines the slope and the roughness of the surface to extract fossil sea-cliffs and fossil rocky shore platforms. For that purpose, contour and hillshaded maps were also analyzed. Then, shoreline angles, a geomorphic marker located at the intersection between the fossil sea-cliff and platform, were mapped analyzing swath profiles perpendicular to the isobaths. Most of the submerged strandlines are irregularly preserved throughout the continental shelf. In summary, 12 submerged terraces with their shoreline angles between approximately: -13 m (T1), -30 and -32 m (T2), -34 and 41 m (T3), -44 and -47 m (T4), -49 and 53 m (T5), -55 and 58 m (T6), -59 and 62 m (T7), -65 and 67 m (T8), -68 and 70 m (T9), -74 and -77 m (T10), -83 and -86 m (T11) and -89 and 92 m (T12). Nevertheless, the ones showing the best lateral continuity and preservation in the central part of the shelf are T3, T4, T5, T7, T8, and T10. The age of the terraces has been estimated using a landscape evolution model. To simulate the formation and preservation of submerged terraces three different scenarios: (i) 20-0 ka; (ii) 128-0 ka; and (iii) 128-20 ka, were compared. The best scenario for terrace generation was between 128 and 20 Ka, where T3, T5, and T7 could have been formed. T3 - Zweitveröffentlichungen der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe - 1414 KW - marine terrace KW - submerged sequence KW - digital bathymetric model KW - TerraceM KW - numerical modeling KW - Bay of Biscay Y1 - 2020 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-517815 SN - 1866-8372 IS - 47 ER -