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 - JOUR 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) JF - Frontiers in Earth Science 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. KW - marine terrace KW - submerged sequence KW - digital bathymetric model KW - TerraceM KW - numerical modeling KW - Bay of Biscay Y1 - 2020 U6 - https://doi.org/10.3389/feart.2020.00047 SN - 2296-6463 VL - 8 IS - 47 SP - 1 EP - 20 PB - Frontiers Media CY - Lausanne 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 - Schreiber, Lisa A1 - Munz, Matthias A1 - Salzmann, Thomas A1 - Oswald, Sascha E. T1 - Coupled simulation of groundwater and drainage dynamics in a coastal fen T1 - Modellierung der Strömungsdynamik in einem revitalisierten Küstenmoorgebiet an der Ostsee JF - Grundwasser : Zeitschrift der Fachsektion Hydrogeologie in der Deutschen Gesellschaft für Geowissenschaften (FH-DGG) N2 - Coastal wetlands are characterized by continued human influence, e.g. with drainage ditches, coastal dikes or landscape restoration. In addition, it is important to understand the complex interactions with the sea to predict impacts of further development. In the present study the aim was to analyze surface and subsurface flow in a coastal wetland located at the Baltic Sea coastline near Warnemunde (Germany) to quantify water exchange with the Baltic Sea and analyze the effect of a storm flood event on saline intrusion. A 3-D transient groundwater model and a one-dimensional surface water model were set up and calibrated by using hydraulic head measurements. The results indicate that in addition to ditch flow, groundwater discharge to the Baltic Sea often has a significant influence on the overall water budget of the fen. From the transient modelling it became evident that water exchange between groundwater in the fen and the Baltic Sea depends on sea level and very often fluctuates between seaward and landward flow directions on daily to weekly time scales. N2 - Küstennahe Niedermoore wurden durch den Menschen verändert, bspw. durch das Anlegen von Entwässerungsgräben, dem Bau von Küstenschutzdeichen oder aktuell einer Renaturierung. Außerdem ist es wichtig die komplexe Interaktion mit der See zu verstehen, um Aussagen über die zukünftige Entwicklung treffen zu können. In der vorliegenden Studie wurde die ober- und unterirdische Strömung in einem Feuchtgebiet an der mecklenburgischen Ostseeküste nahe Warnemünde (Deutschland) untersucht, um dessen wechselseitigen Austausch mit der Ostsee zu quantifizieren und zu analysieren, wie sich ein Sturmhochwasserereignis auf den Salzeintrag ins Gebiet auswirkt. Hierfür wurde ein dreidimensionales instationäres Grundwassermodell erstellt, mit einem eindimensionalen Modell des Grabensystems gekoppelt und mit Messungen im Gebiet kalibriert und verglichen. Die Ergebnisse zeigen, dass neben der oberirdischen Entwässerung auch der Grundwasserabstrom in Richtung Ostsee eine wesentliche Komponente der Wasserbilanz darstellt. Das Verhalten entlang der Küste wird deutlich durch die Dynamik der Ostseewasserstände geprägt, wobei ein Grundwasserabstrom mit einem Zustrom von Ostseewasser bei hohen Küstenwasserständen innerhalb täglicher bis wöchentlicher Zeitskalen wechselt. KW - numerical modeling KW - Baltic Sea coast KW - groundwater-surface water KW - interaction KW - storm flood KW - salinization KW - Numerische Modellierung KW - Ostseeküste KW - Grundwasser-Oberflächenwasser-Interaktion KW - Sturmhochwasser KW - Versalzung Y1 - 2021 U6 - https://doi.org/10.1007/s00767-021-00486-y SN - 1430-483X SN - 1432-1165 VL - 26 IS - 3 SP - 289 EP - 304 PB - Springer CY - Heidelberg ER - TY - JOUR A1 - Heckenbach, Esther Lina A1 - Brune, Sascha A1 - Glerum, Anne C. A1 - Bott, Judith T1 - Is there a speed limit for the thermal steady-state assumption in continental rifts? JF - Geochemistry, geophysics, geosystems : G 3 ; an electronic journal of the earth sciences N2 - The lithosphere is often assumed to reside in a thermal steady-state when quantitatively describing the temperature distribution in continental interiors and sedimentary basins, but also at active plate boundaries. Here, we investigate the applicability limit of this assumption at slowly deforming continental rifts. To this aim, we assess the tectonic thermal imprint in numerical experiments that cover a range of realistic rift configurations. For each model scenario, the deviation from thermal equilibrium is evaluated. This is done by comparing the transient temperature field of every model to a corresponding steady-state model with an identical structural configuration. We find that the validity of the thermal steady-state assumption strongly depends on rift type, divergence velocity, sampling location, and depth within the rift. Maximum differences between transient and steady-state models occur in narrow rifts, at the rift sides, and if the extension rate exceeds 0.5-2 mm/a. Wide rifts, however, reside close to thermal steady-state even for high extension velocities. The transient imprint of rifting appears to be overall negligible for shallow isotherms with a temperature less than 100 degrees C. Contrarily, a steady-state treatment of deep crustal isotherms leads to an underestimation of crustal temperatures, especially for narrow rift settings. Thus, not only relatively fast rifts like the Gulf of Corinth, Red Sea, and Main Ethiopian Rift, but even slow rifts like the Kenya Rift, Rhine Graben, and Rio Grande Rift must be expected to feature a pronounced transient component in the temperature field and to therefore violate the thermal steady-state assumption for deeper crustal isotherms. KW - basin analysis KW - geodynamics KW - numerical modeling KW - rifting KW - thermal KW - modeling Y1 - 2021 U6 - https://doi.org/10.1029/2020GC009577 SN - 1525-2027 VL - 22 IS - 3 PB - Wiley CY - Hoboken, NJ ER - 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 - 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 -