Refine
Year of publication
- 2020 (5) (remove)
Document Type
- Article (3)
- Doctoral Thesis (1)
- Postprint (1)
Language
- English (5)
Is part of the Bibliography
- yes (5)
Keywords
- numerical modeling (5) (remove)
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.
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.
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.
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.
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.