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 -