Lower Jurassic sedimentary successions from the Neuquén Basin, Argentina are unique in the abundance of radiometrically datable material (ash-beds) present, which can be tied to bio- and chemostratigraphic (carbon-isotope) zonations. Here, we present new U-Pb radio isotopic dates, integrated with carbon-isotope and Hg/TOC data, from three localities in Argentina (Arroyo Lapa, Arroyo Serrucho/Las Overas and Chacay Melehue) to generate a biostratigraphically calibrated composite carbon-isotope curve and geochronological framework for the Pliensbachian–Toarcian transition in South America. Using a Bayesian framework we present an age-depth model for this composite record and estimate the age and duration of key intervals extending from the Latest Pliensbachian carbon isotope excursion (CIE) through the Early Toarcian negative CIE. Using a  statistical analysis of all available Karoo and Ferrar Large Igneous Province (LIP ) U-Pb and Ar-Ar radioisotopic ages we create a timeline of the key events and examine the timing of the carbon cycle perturbations specifically looking at potential links to peaks of extrusive emplacement of the Karoo and Ferrar LIP. The geochronological framework is further compared with other available radioisotopic dates from correlative sections, allowing for a more precise constraint and validation of the timing and duration of these Early Jurassic events.

How to cite: Al-Suwaidi, A., Ruhl, M., Kemp, D., Storm, M., Hesselbo, S., Jenkyns, H., Mather, T., Percival, L., and Condon, D.: A Southern Hemisphere Chronostratigraphic Framework for the Pliensbachian–Toarcian Carbon Cycle Perturbations, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-15374, https://doi.org/10.5194/egusphere-egu24-15374, 2024.

Late Paleozoic Ice Age (LPIA), which peaked during the mid Permian, resulting in widespread ice centers across Gondwana during its coldest periods. Assessing the climate change across this glaciation and the following deglaciation interval contribute important data not only in terms of understanding the end-Permian extinction event but also present-day global change. Tasmania, located in a high-latitude setting, forming a bridge between the continents Antarctica and Australia provides a valuable record of the environmental and climatic shifts that occurred in areas proximal to glaciation during the acme and waning stages of the LPIA.

At the time of glaciation, Tasmania was a distinct landmass located within the Paleo-Antarctic Circle at a paleo-latitude of ~78°S. Here we present new high-resolution bulk organic carbon isotope analyses (δ¹³C_TOC), mercury, bulk and trace elemental and sedimentological data combined with palynology and conodont biostratigraphy from the late stage (P3 and P4) of the LPIA Glacial-Deglacial Episode III. We base the data on samples from the Bicheno-5 core, from Eastern Tasmania, which contains approximately 83 meters of middle Permian glaciomarine sediments.

Three negative carbon isotope excursions (CIEs) have been identified in the middle Permian (mPN1, mPN2, and mPN3). The latter two are correlated with the deglaciation episodes in Eastern Australia's glacial intervals P3 and P4 phases. Diamictites and dropstones are typically present in the sediments that record the most positive carbon isotope values, which likely correspond to the peak of the glaciation period. Elemental proxies indicate two cycles of increased weathering and terrestrial sediment influx to the marine system. These cycles coincide with the most negative carbon isotope values and are associated with deglacial cycles in Tasmania and within the paleo-Antarctic circle. The first deglacial cycle (mPN1) coincides with elevated mercury (Hg/TOC) which may hint at a link between deglaciation and volcanism, possibly from the Tarim III LIP.

Comparisons with similar records from the marine Pingdingshan (PDS) section in South China confirms that our data from Tasmania reflect global carbon cycle perturbations providing new insights into the significant global climatic shifts that occurred during the middle Permian. The end stage of the LPIA offers a unique comparison to modern environmental and climatic change in Antarctica associated with anthropogenic global warming.

How to cite: Lestari, W., Al Suwaidi, A., Fox, C., Vajda, V., Ceriani, A., Sun, Y., Frieling, J., and Mather, T.: Unraveling Tasmania's Late Paleozoic Ice Age: Carbon Isotopic and Stratigraphic Signatures in Response to Glacial-Deglacial Cycles and Large Igneous Province (LIP) Events, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-12015, https://doi.org/10.5194/egusphere-egu24-12015, 2024.

Gondwana Basins of eastern India preserves sediment record from Carboniferous to Triassic, which were possibly sourced partially from Antarctica through a radial drainage system. This study attempts to test the hypothesis regarding source of sediments based on petrographical and mineral chemical analysis of siliciclastic Paleo-Mesozoic sediments of eastern India. We present an integrated provenance and paleodrainage analysis on the sediments of Bokaro and Raniganj basins, outcrops of which occur along E-W trend along eastern part of India. The sedimentation in these Gondwana basins initiates with basal Talchir Formation, consisting of alternation between conglomerate and fine- to medium-grained sandstones, and is succeeded by Barakar-Barren Measures-Raniganj and Panchet Formation with sandstone-mudstone alternation, with or without coal.  Petrographic study of sandstones reveals moderate sorting, with angular to sub-rounded quartz and feldspar and rounded to well-rounded lithic fragments; however, the abundance of lithic fragments drastically reduces from Talchir to Panchet formation. Feldspar grains shows the dominance of K-feldspar over plagioclase.  Most of the sandstones are classified as feldspatho-quartzose arenite. The Qm-F-Lt plot indicates that these sandstones were derived primarily from transitional continental sources. Heavy minerals in sandstones include garnets, tourmaline, epidote, rutile, zircon, monazite in order of decreasing abundance. Mineral chemistry of garnet in sandstones points their source to metasedimentary amphibolite facies rocks and granitoid. The tourmaline mineral chemistry suggests the derivation of sediments from various sources, including Li-poor granitoids associated with pegmatites, aplites and Ca-poor metapelites. Rutile chemistry in sandstones indicates the predominance of metapelitic source over metamafic source. Insights from heavy mineral analysis indicates that the Gondwana sediments were derived from multiple sources, and such variation in sources bears information about paleogeographic and paleotectonic evolution of the depositional basin.  The mineral composition of source rocks and paleocurrent data tracks the source of sediments to Eastern Granulite-Schist belt and the Eastern Ghat mobile belt, situated to the east and southwest parts of the Gondwana succession.  This study when integrated with geochronological data would reveal the extent to which a particular source provided sediments to these basins, evolution of sediment sources from bottom to top and ultimately will lead to a refined understanding of timing and evolution of East Gondwana assembly.

How to cite: Dutta, A. and Banerjee, S.: Facies Analysis, Petrography, heavy mineral analysis of paleo-Mesozoic sediments of Eastern India: Implications on provenance and basin evolution, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-18433, https://doi.org/10.5194/egusphere-egu24-18433, 2024.

The basal Cambrian sandstone unit in the North China craton (NCC) is an example of globally widespread siliciclastic succession that resides on the Great Unconformity and deposited on a hypothesized low-lying peneplain during the Cambrian global eustatic sea-level rise. Detrital zircon age signatures from this distinct sequence enable recognition of the ancient drainage system of the NCC in deep time and track its potential linkage with the Gondwana landmass. LA-ICP-MS U–Pb dating of the fossil-calibrated basal Cambrian (Series 2) detrital zircon samples from seven measured sections reveal marked spatial changes in their age signatures that can be divided into three distinct types. The first is generally characterized by the bimodal age populations with broad peaks at ~1.85 Ga and ~2.5 Ga that correlate with the Archean to Paleoproterozoic basement inboard of the NCC. The second is featured by multi-modal distribution with diagnostic Neoproterozoic peaks that correspond to subregional magmatic record. The third also shows multiple-zircon age populations, but yields significant crystallization ages close to the early Cambrian age. Comparing our new data with existing age spectra for the Cambrian strata across the NCC and the northern Gondwana demonstrates that separate drainage systems did exist in the peneplained basement during global Cambrian transgression and the basal Cambrian unit in the NCC was not a part of the far-travelled sand sheet across the northern margin of Gondwana. The most suitable source for Cambrian-aged grains constrained by paleogeographic restoration is the arc terrane developed along the northern margin of the NCC as a result of subduction of the Paleo-Asian oceanic plate. Our new continental-scale detrital zircon provenance signatures in the basal Cambrian unit suggest that the NCC should be considered a discrete continental block separated by the Proto-Tethys Ocean in the Cambrian, rather than an integral part of the northern Gondwana.

How to cite: Wei, R. and Duan, L.: Detrital zircon age signatures of the basal Cambrian sandstone unit in North China: implications for drainage divides during global Sauk transgression and separation from Gondwanaland, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-2721, https://doi.org/10.5194/egusphere-egu24-2721, 2024.

The Transantarctic Basin is a continental basin system developed for ca 200 Myr, from the Devonian to the Early Jurassic, along the Panthalassa margin of Gondwana and above the peneplained Ross Orogeny rocks. The Beacon Supergroup strata form the clastic sedimentary infill of the Transantarctic Basin.

Geodynamic interpretation of the Transantarctic Basin is not univocal and likely geodynamic conditions accounting for basin subsidence have been changed in space and time. Involvement of the basin into the Gondwanian orogenic deformation is a key question for defining the geodynamic setting and the tectonic environment during the Beacon Supergroup deposition. Nonetheless, involvement of  Beacon Supergroup  in orogenic shortening is  poorly assessed for the limited exposed rocks and because formation of the Cenozoic  Transantarctic rift shoulder modified the Paleozoic geometry of the Beacon strata.

Here we explored the exhumation pattern and thermal evolution of the basement rocks and the immediately overlain Beacon sandstones through low-temperature apatite fission track and (U-Th)/He zircon thermochronology along the Prince Albert Mts, where Beacon deposits are particularly thin to suppose a relevant erosional event during the Paleozoic. Thermochronological data and thermal modelling pointed out that basement rocks and Beacon sandstones of the Prince Albert Mts have preserved evidence of a Late Paleozoic erosional event that allows to infer an actively eroding topographic high that lasted from Early to Late Paleozoic times, as a consequence of the far-field stress transmitted from the active Gondwanian convergent margin.

How to cite: Olivetti, V., Cattò, S., Balsamo, F., Zurli, L., Perotti, M., Cornamusini, G., Fioraso, M., Rossetti, F., and Zattin, M.: Uplift and erosion of an intraforeland topographic high: implication for the evolution of the Gondwanian Transantarctic foreland, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-15516, https://doi.org/10.5194/egusphere-egu24-15516, 2024.

A combination of the geochemical and isotopic compositions of the Jurassic igneous rocks exposed along the Antarctic Peninsula and Patagonia suggest that they may have formed in a continental arc setting without a significant input of a mantle plume, as presented in Bastias et al. (2021). This interpretation revises previous interpretations (e.g. Pankhurst et al., 2000; Riley et al., 2001) that suggested the early deposits of the Chon Aike magmatic province in eastern and central Patagonia formed by the coupled action of an active margin and the peripheral thermal effect of the Karoo mantle plume. We present a tectono-magmatic evolution for the Chon Aike magmatic province following the argument of Bastias et al. (2021). Furthermore, we argue that a supra-subduction zone origin is supported by: (i) enriched LILE and LREE, with negative Nb, Ta, Sr and Ti anomalies in all of the Jurassic igneous rocks, which are typical of slab-dehydration reactions and thus active margins, (ii) the trace element compositions of Jurassic igneous rocks in Patagonia and the Antarctic Peninsula are indistinguishable throughout the region, and thus it is likely that they evolved via the same magmatic processes, (iii) the whole rock Nd and Sr, and zircon Hf isotopic compositions of Late Triassic-Jurassic igneous rocks show that the magmas formed from mixed sources that resided within the continental crust and (iv) numerical modelling. Therefore, unlike the igneous units of the Karoo and Ferrar LIPs, the involvement of a mantle plume is not necessary to generate any of the chronological, geochemical or isotopic characteristics of the Chon Aike magmatic province. This contribution serves as further information on the debate searching for the responsible mechanism generating the Chon Aike magmatic province: active margin vs. mantle plume.

How to cite: Bastias-Silva, J., Spikings, R., Riley, T., and Sanhueza, J.: The Chon Aike magmatic province: an active margin origin?, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-21597, https://doi.org/10.5194/egusphere-egu24-21597, 2024.

Mantle convection is driven by buoyancy forces in Earth’s interior. The resulting radial stresses generate vertical deflections of the surface, leaving traces in the geological record. Utilizing new data assimilation techniques, geodynamic inverse models of mantle flow can provide theoretical estimates of these surface processes, which can be tested against geological observations. These inverse models are emerging as powerful tools, providing the potential for tighter constraints on the relevant physical parameters governing mantle flow.

The geodynamic inversions mentioned above require an estimate of the present-day thermal state of the mantle, which can be derived from seismic observations. Using thermodynamically self-consistent models of mantle mineralogy, it is possible to convert the seismic structure imaged by global tomographic models to temperature. However, both seismic and mineralogical models are significantly affected by inherent limitations and different sources of uncertainty. In addition, owing to the complexity of the mineralogical models, the relation between temperature and seismic velocities is highly non-linear and not strictly bijective: In the presence of phase transitions, different temperatures can result in the same seismic velocity, making the conversion from seismic heterogeneity to thermal structure non-unique.

We investigate the theoretical ability to estimate the present-day thermal state of the mantle based on tomographic models in the case of isochemical convection. The temperature distribution from a 3-D mantle circulation model with earth-like convective vigour serves as the “true” temperature field. Using a closed-loop experiment, we aim to recover this initial model after: 1) mineralogical mapping from the “true” temperatures to seismic velocities, 2) application of a tomographic filter to mimic the effect of limited tomographic resolution, and 3) mapping of the “imaged” seismic velocities back to temperatures. We test and quantify the interplay of smoothed seismic structure due to tomographic filtering with different approximations for the conversion from seismic to thermal structure. Additionally, owing to imperfect knowledge of the parameters governing mineral anelasticity, we test the effects of changes to the anelastic correction applied in the mineralogical mapping. The observed mismatch between the recovered and initial temperature field is dominated by the effect of tomographic filtering, with a depth-dependent average error of up to 200 K. Additionally, we observe systematic large errors in the vicinity of phase transitions. Our results highlight that, given the current limitations of tomographic models and the incomplete knowledge of mantle mineralogy, magnitudes and spatial scales of a temperature field obtained from global seismic models will deviate significantly from the true state, even under the assumption of purely thermally driven mantle flow. Strategies to estimate the present-day thermodynamic state of the mantle must be carefully selected to minimize additional uncertainties.

How to cite: Robl, G., Schuberth, B., Papanagnou, I., and Thomas, C.: Deriving mantle temperatures from global seismic models: A quantitative analysis in the light of uncertain mineralogy and limited tomographic resolution, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-2858, https://doi.org/10.5194/egusphere-egu24-2858, 2024.

Interpretations of seismic tomography and applications of the resulting tomographic images, e.g. for estimating present-day mantle temperatures, require information on their resolution and uncertainty. Assessing these model properties is often difficult due to the large size of tomographic systems on global scales. In consequence, there have been only few attempts to consistently analyse the spatially variable quality of tomographic images of deep mantle structure. For linear problems, both resolution and uncertainty can be quantified with the tools provided by classic Backus–Gilbert (B–G) inversion. In this theory, averaging kernels define the local resolving power at each model parameter, while uncertainties represent the propagation of data errors into the model values. By using a more efficient variant of B–G inversion, the method of 'Subtractive Optimally Localized Averages' (SOLA), global tomography can be performed with complete information for model appraisal.

Based on the SOLA framework, we present a concept for the assessment of the 3-D resolution information contained in a global set of averaging kernels. It is based on the rigorous estimation of resolution lengths from a 3-D Gaussian parametrization of the averaging kernels, together with a test for the robustness of this approximation. This is a necessary step because a perfectly bell-shaped or delta-like behaviour of resolution can not always be guaranteed in global tomography due to the inhomogeneous data coverage. Therefore, we also develop a classification scheme, which enables a basic identification of those averaging kernels that are too complex to be sufficiently described by the chosen definition of resolution length. We note that this approach is more generally applicable, i.e. it can be used with any explicitly available set of averaging kernels or point-spread functions, but also with alternative parametrizations.

In the context of the SOLA method, our resolution analysis can be further used to locally calibrate the inversion parameters. This involves on the one hand the specification of a target (resolution) kernel. On the other hand, a trade-off parameter needs to be selected that regulates the fit of the averaging to the target kernel, and the conversely affected propagation of data errors.

To this end, we apply our concept for robust resolution estimation to different sets of averaging kernels from SOLA inversions with varying parameter combinations. Most notably, we systematically increase the spatial extent of the target kernels (taken as 3-D Gaussian functions here as well). The final maps of global (and classified) resolution and uncertainty can be viewed together for a complete picture of the model quality. They reveal where, and for which target size and amount of error propagation, resolution lengths are meaningful and model values can be interpreted appropriately.

How to cite: Freissler, R., Schuberth, B. S. A., and Zaroli, C.: Global assessment of tomographic resolution and uncertainty with the SOLA method, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-10566, https://doi.org/10.5194/egusphere-egu24-10566, 2024.

Fluid dynamics simulations are a powerful tool for understanding processes in the Earth's deep interior. Mantle circulation models (MCMs), for example, provide important insight into the present-day structure of the mantle and its thermodynamic state when coupled with mineralogical models, which is essential information for other fields in the geosciences. The evolution of the heat flux through the core-mantle boundary, for instance, is a prerequisite for geodynamo simulations that aim to model the reversal frequency pattern of the Earth's magnetic field on geologic time scales. However, geodynamical modelling requires extensive knowledge of deep Earth properties and plate motions over time. Uncertainties in these model inputs propagate into the MCMs, which subsequently have to be evaluated with independent data, such as the seismological or geological record. Although state-of-the-art MCMs typically explain statistical properties of seismological data, they do not consistently reproduce the location of features in the mantle.

In this contribution, we explore the effect of varying the absolute position of mantle structure on seismic data by applying first-order modifications to an initial MCM. Normal mode data are particularly well suited for assessing the resulting changes in the location of mantle structure, as they capture its long-wavelength component throughout the entire mantle. In addition, the global sensitivity of normal modes reduces the drawbacks of uneven data coverage. Specifically, we use two different seismic forward modelling approaches, an iterative direct solution method for computing full-coupling spectra and a splitting function calculation that is based on the self-coupling approximation. Our goal is to quantify the effects of a limited number of large-magnitude earthquakes, the adequacy of the self-coupling approximation, and the resolvability of relevant model differences through a comprehensive data analysis. Our synthetic forward modelling framework is moreover well suited for testing the depth sensitivity associated with specific frequency intervals in the spectrum that generally is inferred from seismic 1-D profiles within the splitting function approximation.

How to cite: Schneider, A., Schuberth, B., Koelemeijer, P., Shephard, G., and Al-Attar, D.: Exploring the potential of normal mode seismology for the assessment of geodynamic hypotheses, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-10128, https://doi.org/10.5194/egusphere-egu24-10128, 2024.

Understanding the interaction between oceanic plates and the underlying asthenosphere and its impact on plate thickness is essential for explaining plate motions and mantle convection patterns. While sub-lithospheric small-scale convection provides an explanation for why oceanic plates do not continue to thicken after a certain age, many open questions still surround this process. Here, we link dynamic models of mantle flow, grain-scale processes, seismic imaging, and surface observations to gain new insights into the mechanisms of asthenospheric small-scale convection and its surface expressions.

We have performed a series of high-resolution 3D numerical models of the evolution of oceanic plates and the development of thermal instabilities at their base using the open-source geomodeling software ASPECT. These simulations use an Earth-like rheology that includes coupled diffusion and dislocation creep as well as their interplay with an evolving olivine grain size. Our models quantify how the effective asthenospheric viscosity and the balance between diffusion and dislocation creep affect the morphology and temporal stability of small-scale sub-lithospheric convection, including the age of its onset, the average depth and wavelength of the small-scale convection rolls, and the amplitude of the temperature and grain size anomalies within the rolls.

All of these quantities predicted by the dynamic models can be directly related to both geophysical observables and to surface manifestations such as dynamic topography and heat flux. To accurately compare our model outputs to geophysical data, we convert them to seismic velocity and attenuation using laboratory-derived constitutive relations and taking into account variations in temperature, pressure, grain size, water content and calculated stable melt fraction. We then create synthetic seismic tomography models of different dynamic scenarios and analyze their fit to observations from the Pacific OBS Research into Convecting Asthenosphere (ORCA) experiment. Comparison with both seismic imaging and surface expressions allows us to determine the parameter range in which geodynamic models fit these observations, providing new constraints on the convection patterns and the rheology of the oceanic asthenosphere beneath the Pacific Plate.

How to cite: Dannberg, J., Eilon, Z., Russell, J. B., and Gassmoeller, R.: Sub-Lithospheric Small-Scale Convection as a Window into the Asthenosphere: Insights from Integrating Models Of Mantle Convection, Grain Size Evolution and Seismic Tomography, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-13895, https://doi.org/10.5194/egusphere-egu24-13895, 2024.

The Cascadia Subduction Zone is characterized by young subducting lithosphere, its isolation from other subducting systems, and its ability to produce megathrust earthquakes (M>9.0) and devastating tsunamis. Due to its high potential hazard and risk, it is also a well-studied subduction zone where modern, diverse and detailed observational datasets are available through the USGS and initiatives like GeoPrisms and EarthScope. These datasets include high quality GPS, onshore and offshore geophysical imaging, geochemical and seismic anisotropy data. Integrating these data sets with geodynamic modeling presents an opportunity to gain insight into outstanding questions regarding slab structure, tectonic evolution, seismic hazards, and the physical processes that can self-consistently explain all these observations. For example, geologic and geophysical data suggest that there may be one or two prominent slab gaps or tears, while tomographic data does not fully constrain the depth extent of the slab. Furthermore, the overriding plate is composed of different terranes and contains numerous active and slowly moving faults, complicating efforts to accurately constrain variations in present-day stress and deformation rates.

In this study we test whether comparison of observations to geodynamic model predictions can distinguish between different slab geometries for the Cascadia Subduction Zone. To this end, we have created regional 3D geodynamic models of Cascadia including the slab based on the Slab 2.0 dataset. The model setup is built with the Geodynamic World Builder, and the models are run with the geodynamics code ASPECT. We present results which compare the Juan de Fuca plate velocities against the present day Euler poles. We have found that matching the plate velocity magnitude and direction is sensitive to the rheological model overall, while at the same time being insensitive to certain aspects of the plate boundary rheologies. During the evolution of these models we track the development of the CPO (Crystal Preferred Orientation) with an implementation of the DREX algorithm, so we can compare it against observations of seismic anisotropy in the region. Our presentation will focus on the importance of the geometry of the slab and the strength of different sections of the interface. Furthermore, these models and demonstrate workflows for linking the model results to surface tectonics.

How to cite: Fraters, M., Billen, M., Naliboff, J., Staisch, L., and Watt, J.: Exploring the structure of the Cascadia Subduction Zone by coupling 3D thermomechanical modeling and CPO evolution with observations., EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-16531, https://doi.org/10.5194/egusphere-egu24-16531, 2024.

The topography of the Iberian Meseta located between the Pyrenees-Cantabrian Mountains and the Betics Chain is moderately high (660 m on average) compared to the high plateaus on Earth (> 1 km) albeit higher than the topography of the surrounding western European plate. The Iberian Meseta encompasses Cenozoic sedimentary basins, active or inactive Alpine mountain ranges, and low-relief erosional surfaces represented by plateaus with elevations between 600 and 1400 m asl. It is commonly accepted that ~600-700 m surface uplift occurred during the Cenozoic, but the underlaying processes and the precise timing of the onset of the plateau growth are strongly debated. The main objective of the present study is to find out to what extent the topography of the Iberian Meseta has a crustal, lithospheric or a sublithospheric origin. We used the results of a recent modelling based on the joint inversion of both the crustal and lithospheric mantle structure. It encompasses an integrated geophysical-lithological multi-data modelling. The inversion is framed within an integrated geophysical-petrological setting where mantle seismic velocities and densities are computed as a function of temperature and composition whereas crustal density, shear and compressional wave velocities are lithologically linked based on empirical relationships from global petrophysical databases.

We computed the relative contribution to topography of crustal and lithospheric mantle thickness variations and density structure. The topography of the Alpine mountain belts in Iberia is largely associated with thickened crust. We find that the elevated topography in the NW Iberian Meseta (elevation > 700 m) is mostly related to the lithospheric mantle thinning. This is in agreement with Inversion of topographic data, landform dates, and erosion rates suggesting a late Cenozoic mantle-related surface uplift of several hundreds of meters in NW Iberia (EGU24-11613) and in the central Iberia (EGU24-16382). Similarly, a thin and warm lithospheric mantle is responsible for the positive elevation of the onshore Mediterranean margins. A negatively buoyant lithospheric mantle causes >1 km subsidence in the Gibraltar Arc and western Pyrenees.  We solved the Stokes flow to evaluate the contribution of temperature-related buoyancy forces at asthenospheric depths. These forces cause a long wavelength topographic response located in the centre of the Iberian Peninsula reaching a maximum value of only 100-150 m, which is much lower than the values reported in previous works assuming an isostatic balance.

How to cite: Negredo, A. M., Fullea, J., Ortega-Gelabert, O., Clemente, C., and Babault, J.: On the different contributions to the peculiar topography of the Iberian Peninsula, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-18172, https://doi.org/10.5194/egusphere-egu24-18172, 2024.

Mantle convection is a fundamental process governing the evolution of our planet. Buoyancies in the mantle induce horizontal and vertical motion of the Earth’s lithosphere, which can be mapped using independent geological datasets. Positive surface deflections induced by mantle convection create erosional/non-depositional environments which lead to gaps (hiatuses) in the stratigraphic record, while negative deflections provide accommodation space for sedimentation to occur. We use continental- and country-scale digital geological maps and regional and local stratigraphic studies at the temporal resolution of geological series (ten to tens of millions of years) to map the distribution of hiatus through geological time.

Here we present global continental hiatus surfaces since the Upper Jurassic. We find that they vary inter-regionally at timescales of geological series and that they correlate with known mantle dynamic events. For example, we tend to observe the appearance of a hiatus surface before the arrival of a mantle plume. In Europe, we mapped a large-scale sedimentary hiatus during the Paleocene, prior to the arrival of the Iceland plume. In Africa and South America, we found a widespread absence of the Upper Jurassic prior to the arrival of the Tristan plume. This pattern can be seen as characteristic of plume-induced dynamic uplift. We observe a sea level signal during some geological series, such as in the Oligocene, when there is a global increase of hiatus areas, coinciding with the onset of Antarctic glaciation and associated sea level drop. At other times, we find that the hiatus areas evolve differently for different continents, precluding their interpretation as an eustatic signal. Spectral analysis shows that hiatus surfaces have shorter wavelengths than no hiatus surfaces, requiring higher spherical harmonic degrees to describe the geological series with larger amounts of hiatus. These include the Upper Jurassic, the Paleocene, the Oligocene and the Pleistocene.

Our results imply that a key property of time-dependent geodynamic Earth models must be a difference in timescale between mantle convection itself and resulting dynamic topography. Moreover, they highlight the importance of continental-scale compilations of geological data to map the temporal evolution of mantle flow beneath the lithosphere, which can provide powerful constraints for global geodynamic models.

How to cite: Vilacís, B., Brown, H., Carena, S., Bunge, H.-P., Hayek, J. N., Stotz, I. L., and Friedrich, A. M.: Global continental hiatus surfaces as a proxy for tracking dynamic topography since the Upper Jurassic, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-12118, https://doi.org/10.5194/egusphere-egu24-12118, 2024.

How to cite: Brown, H., Vilacís, B., Stotz, I., Chen, Y.-W., and Bunge, H.-P.: Synthetic Hiatus Maps as a Tool for Constraining Global Mantle Circulation Models, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-4552, https://doi.org/10.5194/egusphere-egu24-4552, 2024.

The East Africa - Arabia topographic swell is an anomalously high-elevation region of ~4000 km long (from southern Ethiopia to Jordan) and ~1500 km wide (from Egypt to Saudi Arabia) extent. The swell is dissected by the Main Ethiopian, Red Sea, and Gulf of Aden rifts, and characterized by widespread basaltic volcanic deposits emplaced from the Eocene to the present. Although most agree that mantle plumes play a role in generating the swell, several issues including the number and locations of plumes and the uplift signatures remain debated. We seek to address these questions and provide a general evolutionary model of the region. To this end, we conduct a quantitative analysis of topography to infer isostatic and dynamic contributions. When interpreted jointly with geological data including volcanic deposits, the constraints imply causation by a single process which shaped the past and present topography of the study area: the upwelling of the Afar superplume. Once hot mantle material reached the base of the lithosphere below the Horn of Africa during the Late Eocene, the plume flowed laterally toward the Levant area guided by pre-existing discontinuities in the Early Miocene. Plume material reached the Anatolian Plateau in the Late Miocene after slab break-off and the consequent formation of a slab window. During plume material advance, buoyancy forces led to the formation of the topographic swell and tilting of the Arabia Peninsula. The persistence of mantle support beneath the study area for tens of million years also affected the formation and evolution of the Nile and Euphrates-Tigris fluvial networks. Subsequently, surface processes, tectonics, and volcanism partly modified the initial topography and shaped the present-day landscape.

How to cite: Sembroni, A., Faccenna, C., Becker, T. W., and Molin, P.: The uplift of East Africa-Arabia swell: the signature of the mantle upwelling and spreading, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-18733, https://doi.org/10.5194/egusphere-egu24-18733, 2024.

Recent advances in computational capabilities make it possible to compute global geodynamic earth models at near earthlike convective vigor. This paves the way to systematically obtain a range of synthetic data from such models in an approach that is known as closed loop experiments. Here we present results from closed loop experiments in geodynamic earth models targeted at three classes of data that are sensitive to the mantle convection process, namely seismic data, global stress patterns as reflected by the world stress map, and continent scale stratigraphy processed for the distribution of conformable and unconformable successions in recently developed so called hiatus maps. Our results reveal effects from spatially variable data collection and quality (as expected), mantle flow geometries (less expected) and (still poorly known) histories of paleo mantle flow. We conclude that the derivation of process based synthetic data from geodynamic earth models provides crucial information for data interpretion, that closed loop experiments are
a powerful tool to link geodynamic earth models to data, and that closed loop experiments could be helpful to guide future data collection efforts.

How to cite: Bunge, H.-P., Stotz, I. L., Hayek, N., vilacis, B., brown, H., freissler, R., schuberth, B., carena, S., and friedrich, A.: Closed loop experiments in global geodynamic earth models, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-21895, https://doi.org/10.5194/egusphere-egu24-21895, 2024.

Global plate reconstructions that constrain the surface plate motions provide crucial boundary conditions for mantle circulation models. Earth-like plate kinematics could reduce the impact of the uncertain mantle initial conditions and generate slab structures in the mantle that are comparable with seismic tomographies. However, due to the subduction of the oceanic plates, uncertainty increases in global plate reconstruction over time. Here, we utilize a novel slab unfolding technique to retrodeform mantle slabs imaged in the MITP08 seismic topography back to the pre-subduction states at Earth’s surface. Such a technique provides additional constraints on plate reconstructions, especially in regions dominated by intra-oceanic subductions, such as Southeast Asia.

Our reconstruction shows a significant trench retreat along Southeast Asia and Northern Australia between 90 and 65 Ma that opened a gigantic, >3,000 km wide backarc basin. This basin, named the East Asian Sea plate, was later consumed by the west-moving Philippine Sea plate and North-moving Australian plate in the Cenozoic. We then embed our reconstruction in a mantle circulation model, TERRA, testing the fidelity of the reconstruction in a closed-loop experiment.

We found that fragmented, sub-horizontal East Asian Sea slabs can be reproduced in the mantle circulation model. These slabs lying underneath the current Philippine Sea plate and northern Australia are similar to the MITP08 tomography on which the reconstruction is built. Moreover, these slabs at 800-1000 km depths result in a more negative dynamic topography on the present Philippine Sea plates comparable with the observed residual topography. On the contrary, the traditional, Andean-style reconstruction can only produce positive dynamic topography. Other Mesozoic, intra-oceanic subductions in NE Asia and western North America embedded in our reconstruction also produce negative, yet smaller magnitude, dynamic topography, possibly due to the older subduction history and deeper slabs. We conclude the negative dynamic topography within the present Pacific plate is the result of ancient intra-oceanic subductions. 

How to cite: Chen, Y.-W., Bunge, H.-P., Stotz, I., and Wu, J.: Testing tomography-based plate reconstructions from a paired, inverse-forward closed-loop experiment in a mantle circulation model, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-3083, https://doi.org/10.5194/egusphere-egu24-3083, 2024.

It is widely recognized that mantle dynamics and plate flexure both contribute to Earth’s topography and gravity fields at different wavelengths, yet the actual transition wavelength between them is not well quantified, ranging from ~100 km to ~1000 km. Here we use the observed relationship between topography and the free-air gravity anomaly fields (admittance) to infer the relative contribution of plate flexure and mantle dynamics based on mantle flow models which incorporate an essentially elastic plate. Global and regional Pacific Ocean data studies show that plate flexure and mantle convection potentially contribute to the topography and gravity for wavelengths larger than ~600 km. Plate flexure mainly contributes at wavelengths shorter than ~600 km and is consistent with the support of uncompensated topography for wavelengths shorter than ~200 km. To investigate the admittance associated with mantle dynamics at long wavelengths we have constructed mantle flow models based on a number of different seismic tomography models. The finite element software CitcomS was used to calculate mantle flow and related surface dynamic topography and associated free-air gravity anomaly. Admittance analysis in the Pacific Ocean from different mantle flow models show that the dynamic admittance is generally larger than the observed admittance, while the admittance from plate flexure is smaller than the observed admittance, suggesting that both mantle dynamics and plate flexure contribute to Earth’s topography and gravity at long-wavelengths. The difference between the dynamic admittance and the observed admittance is smallest for cases with temperature-dependent viscosity and weak asthenosphere, and the combined admittance in the presence of both flexure and mantle convection for these cases is generally consistent with the observed admittance. We use a new method to separate the effects of plate flexure and mantle convection to the topography and gravity fields at long wavelengths which has been developed from the plate flexure and dynamic admittances. We assume that the topography and gravity at long wavelength are the combination of plate flexure and mantle dynamics and further assume that the topography and gravity are linearly related through the admittance. The final separated dynamic topography in the Pacific Ocean is generally consistent with previously published residual topography studies at long-wavelengths.

How to cite: Yang, A., Watts, A., and Zhong, S.: Long-wavelength gravity and topography of the Pacific Ocean: Relative contribution of mantle dynamics and plate flexure , EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-9346, https://doi.org/10.5194/egusphere-egu24-9346, 2024.

While the isostatic compensation of crustal thickness and density heterogeneity provides the dominant contribution to Earth’s observed topography, there nonetheless remains a substantial difference (the ‘residual topography’) between these two fields. This difference is a consequence of dynamic processes occurring within the mantle, most notably due to time-dependent vertical surface stresses driven by mantle convection. Mantle convection dynamics also produce differences between the observed geoid and the isostatic geoid generated by crustal heterogeneity: the ‘residual geoid’. The joint consideration of both residual geoid and topography anomalies provides unique and fundamental global constraints on the amplitude and spatial distribution of density anomalies in the convecting mantle.
Despite the crucial role of isostasy in determining residual geoid and residual topography, accurate constraints depend heavily on the quality of the isostasy calculations. The classical theory of isostasy relies on a 1^st-order treatment of hydrostatic equilibrium, which is not sufficiently accurate for the calculation of isostatic geoid anomalies on a compressible, self-gravitating mantle. Consequently, we present a geodynamically consistent approach that is based on the surface loading response (via dynamic kernels) calculated with a viscous flow model that incorporates a fully compressible mantle and core (given by the PREM reference model) with self-gravitation.
Another critical issue that remains outstanding is the accuracy inherent in global crustal heterogeneity models. Here we show that the differences between the residual geoid and topography fields predicted using CRUST1.0 (Laske et al. 2012) and the most recent ECM1 (Mooney et al. 2023) crustal heterogeneity models are substantial. We discuss the importance and implications of these differences in the context of determining the most accurate constraints on density anomalies in the convecting mantle.

How to cite: Kamali Lima, S., Forte, A. M., and Greff, M.: A Geodynamically Consistent Approach to Residual Topography and Geoid Anomalies on the Convecting Mantle: Importance of Global Crustal Models , EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-12457, https://doi.org/10.5194/egusphere-egu24-12457, 2024.

Deflection of the Earth’s surface supported by mantle flow, known as dynamic topography, is associated with a free-air gravity anomaly because such topography is not isostatically compensated. Consequently, the ratio of the gravity anomaly to the dynamic topography, known as the admittance, has been used to estimate the amplitude of dynamic topography, which can be difficult to measure directly. However, at long wavelengths (e.g., spherical harmonic degrees 2 to 6) both dynamic topography and gravity anomalies, and thus the admittance, are sensitive to the mantle’s viscosity structure. Previous studies [e.g., Colli et al., 2016] demonstrate a reversal in sign of the free-air gravity anomaly resulting from lower mantle structures as the viscosity of the lower mantle is increased. This indicates potential complexity for inferring long-wavelength dynamic topography from observations of gravity anomalies, because the upper-lower mantle viscosity contrast is poorly constrained. We further investigated the relationship between dynamic topography and gravity anomalies by introducing lateral viscosity variations into a finite element model of global mantle flow. We find that the gravity anomaly above lower mantle density heterogeneity can change dramatically as we begin to introduce different models for lateral viscosity variations into the upper and lower mantle viscosity fields. In such models we find that the sign of the admittance varies laterally, with the horizontal gradients in mantle viscosity perturbing mantle flow patterns in ways that produce large changes gravity anomalies and smaller changes in dynamic topography. A spatially-varying admittance will greatly complicate estimation of dynamic topography from observed gravity, and may help to explain mismatches between observations of dynamic topography and predictions made using global mantle flow models. On the other hand, the reconciliation of such mismatches may help to constrain viscosity heterogeneity in the lower mantle.

Colli, L., S. Ghelichkhan, and H. P. Bunge (2016), On the ratio of dynamic topography and gravity anomalies in a dynamic Earth, Geophysical Research Letters, 43(6), 2510-2516, doi:10.1002/2016GL067929.

How to cite: Conrad, C. P. and Ramirez, F.: Sensitivity of Long-Wavelength Dynamic Topography and Free-Air Gravity to Lateral Variations in Lower Mantle Viscosity, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-6721, https://doi.org/10.5194/egusphere-egu24-6721, 2024.

Mantle convection is an essential component of the Earth system. Yet its history is not well known , in part, as the strength of tectonic plates conceals the underlying flow. To date, global mantle convection models have reached an impressive level of sophistication due to significant advancement of computational infrastructures and numerical techniques. This ultimately allows geodynamicists to reconstruct past mantle states through using, for instance, inverse geodynamic models based on adjoint equations. However, key input parameters of these models --- such as thermo-chemical flow properties and rheology --- are complex and poorly known. This in turn limits their ability to effectively interpret the reconstructed mantle flow, thus motivating one to pursue an approach that aims to conceptualize paleo-mantle-flow at a simple analytical level.

To this end, the existence of thin, mechanically weak asthenosphere permits one to develop an analytical Couette/Poiseuille model of asthenospheric flow, where flow is associated with moving tectonic plates, and with lateral pressure gradients due to rising plumes and sinking slabs. Here we present estimates of the Cenozoic asthenospheric flow history from such models in the Atlantic realm. We, moreover, link them to azimuthal seismic anisotropy as well as mantle flow retrodiction simulated by inverse geodynamic models. Our analytically derived asthenospheric flow indicates that it is in broad agreement with the orientation of seismic azimuthal anisotropy, and with the large-scale flow patterns from mantle flow retrodictions. In light of these results, our study suggests exploiting a hierarchy of geodynamic models together with growing observational constraints on mantle flow induced surface expressions to gain a better understanding of paleo-mantle-flow.

How to cite: Wang, Z. R., Stotz, I. L., Bunge, H.-P., Vilacís, B., Hayek, J. N., Ghelichkhan, S., and Lebedev, S.: Cenozoic asthenospheric flow history in the Atlantic realm: Insights from Couette/Poiseuille flow models, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-2321, https://doi.org/10.5194/egusphere-egu24-2321, 2024.

Peninsula India and Scandinavia are elevated passive continental margins (EPCMs) characterized by asymmetric relief with high mountains in the west and a gentle slope towards lowlands in the east.

New AFTA data from southern India reveal major Phanerozoic episodes of cooling, reflecting exhumation. Here we focus on the early Miocene episode possibly related to the hard India-Asia collision (van Hinsbergen et al. 2012). The Miocene exhumation resulted in a low-relief landscape; e.g. the Karnataka and Mysore plateaus (Gunnel and Fleitout, 1998) with residual regions of higher ground (e.g. Palani Hills). Today, these plateaus reach an elevation of 1 km along the coast of SW India, sloping towards the east. The Miocene peneplains were graded towards the base level of the adjacent ocean (Green et al. 2013), and therefore reached their present elevation after their formation. Thick piles of Pliocene sediments off SW India (Campanile et al. 2008) suggests that this happened during the Pliocene.

Richards et al. (2016) studied river profiles in Peninsula India and concluded that the regional tilt grew since 25 Ma, maintained by sub-lithospheric processes. However, we find that the relief is the result of two episodes: 1) Miocene peneplanation related to far-field stress. 2) Late Neogene, asymmetric uplift driven by sub-lithospheric processes.

We identified a similar development in SW Scandinavia, where two Neogene episodes of uplift and erosion define main features of the relief (Japsen et al. 2018): 1) Early Miocene uplift leading to formation of the Hardangervidda peneplain (possibly related to the hard India-Asia collision). 2) Uplift beginning in the Pliocene, raising Hardangervidda to its present elevation at 1.2 km. Pliocene uplift raised margins around the NE Atlantic with maximum elevations reached close to Iceland. This suggests support from the Iceland Plume due to outward-flowing asthenosphere extending beneath the conjugate margins (Rickers et al. 2013; Japsen et al. 2024). 
Lithospheric as well as sub-lithospheric processes appear to shape main features of EPCMs.

Campanile et al. 2008. Basin Research. Green et al. 2013. GEUS Bull. Gunnell, Fleitout 1998. ESPL. Japsen et al. 2018. JGSL. Japsen et al. 2024. ESR. Richards et al. 2016. G cubed. Rickers et al. 2013. EPSL.

How to cite: Japsen, P., Green, P. F., Bonow, J. M., and Chalmers, J. A.: The large-scale landscapes in SW Scandinavia and in SW India are the result of two episodes of Neogene uplift, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-6544, https://doi.org/10.5194/egusphere-egu24-6544, 2024.

Calderas are steep morphological and collapse basins that continue to reshape long after initial structural collapse. While large landslides are associated with caldera collapse and widen the basin, little is known about the morphological and structural changes that occur long after caldera formation. Here, we investigate the shape and slope of the Tambora caldera in Indonesia, which formed in 1815 during one of the most devastating eruptions of the past  centuries. The release of over 150 cubic kilometers of volcanic ash created a caldera 6 kilometers wide and 1250 meters deep, causing climatic effects worldwide. Here we explore an ultra-high-resolution dataset we generated from Pleiades, a tri-stereo satellite, that now allows us to apply computer vision approaches to study the morphology and geometry of the Tambora caldera. We generated a 12 million pixel point cloud resampled to a 1 m resolution Digital Elevation Model and a 0.5 m orthomosaic. We explore the dimension, slope, and outline of the caldera and find localized open fissures, tension cracks, and morphological scars. We also apply an unsupervised image classification approach to the stereo multispectral data and find locations of fumarole activity and hydrothermal alteration in close proximity to these structural features. Hydrothermal alteration sites are commonly located in the caldera wall below the scars and open fissures. We explore this proximity of alteration, scarring, and faulting using newDistinct Element Method models, emphasizing that caldera morphology and structure is strongly influenced by hydrothermal weakening that causes flank instability, localized shedding of material, and large-scale morphological changes.

How to cite: Walter, T. R., Harnett, C. E., Troll, V. R., and Heap, M. J.: Morphology, structure and gravitational instability of the steep caldera walls of Tambora (Indonesia) influenced by hydrothermal alteration., EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-2023, https://doi.org/10.5194/egusphere-egu24-2023, 2024.

The Cheakamus basalts are a set of ~31 km long, 1.65 km³ valley-filling olivine basalts which erupted into the glaciated Callaghan-Cheakamus Valley system of the Garibaldi Volcanic Belt (GVB), British Columbia. Combined paleomagnetic and radiometric (⁴⁰Ar/³⁹Ar) analysis dates the lavas as a short-lived (< 2 ka duration) eruption at 15.95 ± 7.9 ka (2σ); additional field evidence, including well-glaciated lava flow surfaces overlain by till, indicate the eruption coincided with the early stages of the Fraser Glaciation (LGM) at ~20-18 ka. The lavas preserve features indicative of a landscape hosting diverse and dynamic paleoenvironments. Subaerial eruption of basalt lava filled an ice-free Callaghan Creek drainage system before inundating and damming of the paleo-Cheakamus River, creating an upstream rising body of water. Periodic overtopping of the lava dam resulted in syn-eruptive intermittent flooding and overtopping of lavas expressed by discontinuous lenses of interflow sediment and well-developed entablatures in the upper portions of lava flows. Rare instances of enigmatic cooling columns indicate localized ice contact with glaciers that partially filled the Cheakamus Valley. Emplacement features and morphologies in the Cheakamus Valley have been heavily altered by the erosional overprinting of a glacial lake outburst flood (GLOF) that preliminary ¹⁰Be analysis dates at the close of the LGM (11-10 ka). Lavas in the Callaghan valley remain untouched by the GLOF. Their aerial extent and distribution, especially at high elevations in tributary valley mouths suggest bottlenecking and backing-up of lavas due to narrowing in valley topography. Current work combines field mapping and analogue modelling and aims to provide insight into the emplacement dynamics of effusive lavas in the steep, confined terrain of the BC Coast Mountains. Despite, and in part because of, their heavily modified morphology, the Cheakamus basalts act as an excellent recorder of both effusive volcanic processes and the paleoenvironments into which they erupted. Their thorough and accurate analysis is especially pertinent temporally, as they erupted during a time of glacial flux and can provide additional evidence for the timing and location of the advancing Cordilleran ice sheet. Spatially, the Cheakamus basalts are proximal to population centers and transport infrastructure and thus have implications for potential volcanic hazards and attendant risks, as any future effusive, valley-filling basaltic eruption from the GVB will likely share similar emplacement characteristics and processes.

How to cite: Borch, A., Russell, J. K., Barendregt, R., and Friele, P.: Emplacement and erosion of valley-filling basalt lavas in shifting Quaternary environments of the Garibaldi Volcanic Belt, British Columbia, Canada, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-4176, https://doi.org/10.5194/egusphere-egu24-4176, 2024.

The influence of climate on landscape evolution in natural settings remains debated. Here, we focus on tropical hotspot volcanic islands because they exhibit relatively uniform lithology and experience significant precipitation and climate gradients. Furthermore, intermittent volcanic flows effectively “reset” the landscape that begins to evolve post-eruption. Thus, we can constrain initial conditions by reconstructing the initial geometry of radiometrically dated volcanic flows. We constrain landscape evolution through time in several volcanic islands with strong climate gradients to assess the role of climate on incision. We perform topographic reconstructions to calculate long-term basin-averaged erosion rates in two islands of the Réunion hotspot (Réunion, Mauritius) and compile published erosion rates on Réunion and Kaua’i (Hawaii hotspot). We define the time since incision started as the age of the surface incised lava flow. To calibrate incision parameters on all three islands, we use the stream power model and apply a data-driven Bayesian approach to obtain the erodibility, K, the drainage area exponent, m, and the slope exponent, n. We also calculate a normalized erodibility index, K_n, using n = 1 to directly compare results among the different study sites. Erosion rates of Réunion Island range from 9.9 ± 0.5 mm/yr to 5.2 x 10^-3 ± 2 x 10^-4 mm/yr and erosion rates in Mauritius Island range from 6.5 x 10^-2 ± 8 x 10^-3 mm/yr to 5.1 x 10^-3 ± 4 x 10^-4 mm/yr. Incision efficiency seems to decrease with time since incision started from 63 ka to ~300 ka and then does not vary significantly with time since incision started from ~300 ka to 4300 ka. This is likely due to the covariation between the age of volcanism repaving and precipitation rates on Réunion, which is related to the configuration of the island’s two volcanoes – the active Piton de la Fournaise located on the windward side and the dormant Piton des Neiges on the center and leeward side. Our empirical calibration of the stream power law shows high dispersion in n and K_n on each individual island. m ranges from 0.2 to 2.9, and K_n ranges from 2.3 x 10^-7 to 9.8 x 10^-4 m^1-2m/yr. For Réunion, we identify a positive trend between mean annual precipitation and erosion rates, and between mean annual precipitation and K_n, for low to moderate erosion rates (<1 mm/yr). For all basins of Réunion, there is also a positive trend between mean annual cyclonic precipitation rates and erosion rates, and between mean annual cyclonic precipitation rates and n. In Kaua’i, there is a positive trend between erosion rates and mean annual precipitation, consistent with previous studies. In Réunion, the proportionality coefficient between erosion and mean annual precipitation is three times greater than in Kaua’i. In addition, considering all three islands, a nonlinear relationship exists between channel slope and incision rate: best-fit n values range from 0.5 to 6, with n generally lower than one on Kaua’i. Our results highlight different sensitivities of fluvial relief to incision, and of incision to climate, between Kaua’i and Réunion islands.

How to cite: Gourbet, L., Gallen, S. F., Famin, V., Michon, L., Ramanitra, M. O., and Gayer, E.: River incision on hotspot volcanoes: insights from paleotopographic reconstructions and numerical modelling, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-2631, https://doi.org/10.5194/egusphere-egu24-2631, 2024.

 Active volcanic catchments within low-latitude tropical humid climates constitute some of the most dynamic and sediment-rich fluvial systems globally. The combination of factors such as active explosive vulcanism, frequent high-magnitude earthquakes, dense biodiverse vegetation, and intense rainfall leads to very high sediment supply and very active channel morphodynamics within such fluvial systems.

Our research focused on channel dynamics in three remote and challenging-to-access reaches of the Sucio River, situated within the Irazú stratovolcanic structure in Costa Rica's central volcanic chain. This unique setting experiences extreme conditions, including rainfall exceeding 8000 mm annually, infrequent but significant volcanic eruptions (occurring over three times per century), high-magnitude earthquakes (>Mw5), and dense pristine vegetation. The mapped river reaches within this active volcano exhibit a distinctive confined multi-thread channel morphology predominantly comprised of coarse sediments, notably boulders, displaying exceptional dynamism. These reaches showcase rapid shifts between braided, island-braided, and anabranching morphologies in relatively short periods (<20 yrs.). Additionally, the primary sediment sources located in crateric areas have undergone rapid and substantial changes, resulting in the emergence of large landslides and drastic alterations in vegetation.

Employing remote sensing techniques, geostatistical analysis, and fieldwork, we investigated the impacts of eruptions, earthquakes, and rainfall on the Sucio River's channel morphology and its primary sediment sources (craters) from 1961 to 2023. Over 65 images were processed to generate various derived raster products, including supervised classification datasets, change detection outputs, and morphometric parameters (such as channel width, braided index, anabranching index, and area of bars and islands). Moreover, we constructed a precipitation database spanning the study period to assess the frequency, magnitude, and duration of extreme rainfall events. Historical seismic data was utilized to compile a database of relevant earthquakes that might have affected the catchment, given the river's proximity to several active faults. Subsequently, exploratory statistical analysis employing linear regression models helped discern the influential factors behind channel dynamics and changes.

Our findings provide a novel understanding of how this specific fluvial volcanic environment responds to external perturbations and adapts its channel morphology over time. Key outcomes include the rarity of the multithread boulder morphologies observed in these reaches, rapid morphological transformations occurring within this multithread system in short intervals, the significant role of dense pristine vegetation in stabilizing banks and islands, and the cyclic stability-instability phases (erosion-deposition) triggered by pivotal events like eruptions, hurricanes, and high-magnitude earthquakes.

This study presents novel insights into channel morphology dynamics in one of Central America's most extensively studied active volcanoes, likely having the river transporting the most sediments in Costa Rica's volcanic regions.

How to cite: Granados, S., Surian, N., and Alvarado, G. E.: River dynamics in an active volcano with tropical humid conditions, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-11817, https://doi.org/10.5194/egusphere-egu24-11817, 2024.

A consequence of alluvial processes acting on scoria cones is the development of a drainage network composed of radially distributed rills and gullies parallel to the volcanic edifice's downslope direction. This work focuses on the quantification of the degree of development of the drainage network by applying the Average Erosion Index method to scoria cones from the arid to semi-arid Lunar Crater volcanic field and comparing with previously obtained results from two tropical volcanic fields (Sierra Chichinautzin volcanic field and Paricutin-Tancitaro region, both in central Mexico). The results show that the method helps to determine geomorphic age relations when calibrated separately for each field. Furthermore, the differences in the resultant rates at which AEI varies as a function of time obtained for the three studied fields indicate that the method provides a tool to quantify the effects of different alluvial rates at volcanic fields across various environments.

How to cite: Zarazúa-Carbajal, M. C., Valentine, G. A., and De la Cruz-Reyna, S.: Gully incision development on scoria cones: different behaviors in three volcanic fields reflect environmental conditions., EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-6783, https://doi.org/10.5194/egusphere-egu24-6783, 2024.

Analogue experiments can enhance our understanding of complex natural volcanic landscapes formed by eruptions, intrusions, remobilisation of volcanic material and erosional processes. Experimental setup in a laboratory offers a controlled setting to investigate the development of rainfall-induced radial drainage basins on scaled volcano cones. It allows to simulate surface runoff, a prominent sediment transport process in volcanic landscapes, primarily influenced by climate, lithology, and topography. By controlling the flowrate within the setup, maintaining a uniform lithology, and using initial axisymmetric cones with the same size, this study aims to record the variations in erosion patterns caused by systematic cone slope and shape changes.

Analogue volcanic cones made up of water-saturated 70 μm silica powder were built upon a drainage layer of coarse sand at the VUB volcanology analogue laboratory. The cones were scaled based on the height/basal width ratios of natural pristine composite volcanoes with a scaling factor of 4.5*10^-6 for the basal width. Initial cones had basal widths of 33 cm with two sets of cones with heights ranging from 4.2 to 6.9 cm. For the highest cones, the lower flank was 21 degrees and the upper slope 30 degrees, with a break-in-slope at 45% of the cone height. Experiments included cones with and without summit craters, the craters were 5 cm wide and 0.5 cm deep. Rainfall-induced erosion was simulated with two atomizer sprinklers, creating a mist of droplets of circa 30 μm. Experiments were run for 3 to 5 hours, simulating erosion taking place over several millions of years at natural volcanoes. We generated a minimum of ten Digital Elevation Models (DEMs) by photogrammetry with sub-millimetre spatial resolution, enabling the estimation of volume loss and erosion rates. The automated algorithms MorVolc and DrainageVolc were used to extract morphometric and drainage parameters (e.g., height/basal width ratio, drainage density, irregularity index) from the DEM of each timestep.

The analogue models' drainage networks and morphological characteristics replicate those found on natural volcanoes. Having a steeper slope for the upper flank of the cone delayed the forming of erosional features on the lower flank, while the top part of the volcano incised deeper than the cones with one slope gradient. The cones without a summit crater develop a radial drainage network from an initial set of narrow gullies to a more stable pattern with fewer valleys that gradually widen. The introduction of a summit crater substantially modifies the resulting erosional patterns: the incision of the crater rim forms two to four dominant watersheds that widen faster than the basins of cones lacking a crater. Migration of drainage divides ceases when equilibrium in the landscape is attained, with cones featuring a summit crater reaching this equilibrium later than those without. Analogue experiments are a valuable tool for studying erosional processes in a controlled manner and give insight into complex volcanic landscapes, thereby improving our understanding of long-term volcanic landscape evolution.

How to cite: van Wees, R. M. J., Paguican, E., O'Hara, D., Kereszturi, G., Grosse, P., Lahitte, P., and Kervyn, M.: Modelling volcano degradation through analogue experiments: the impact of volcano slopes and summit craters on erosion patterns, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-13012, https://doi.org/10.5194/egusphere-egu24-13012, 2024.

We present a new book aimed at graduate students, and academics, as well as all volcano enthusiasts. We chose to publish on the topic of volcano geomorphology for two main reasons. Firstly, although several geomorphology textbooks have been published, ones focussed on volcano geomorphology are scarce or outdated and out of print (Cotton, 1944, McDonald, 1972; Ollier, 1969 and 1988). Secondly, many volcanology books have been published over the past few decades, but they do not describe landforms and geomorphic processes in sufficient detail (as stated by the late P. Francis in 1993). To our knowledge, only five modern books on volcanology include a chapter on volcanoes as landforms and landscapes (Francis, 1993 and Francis & Oppenheimer, 2004; Chester, 1993; Scarth, 1993; Lockwood and Hazlett, 2010). They are less process-oriented than modern books on geomorphology.

Our book on Volcano Geomorphology is organised into five main themes, and contains 10 chapters. The first theme is an overview of the geodynamic environments in which the Earth's volcanoes are created.  The second is a detailed account of elementary “constructional” landforms, from lava forms to monogenetic volcanoes, both terrestrial and subaqueous, reflecting a variable degree of magma-water interaction. The third deals with polygenetic volcanic edifices including shield volcanoes, composite cones and volcanic clusters. This is followed by landforms and processes that form calderas, caldera complexes, and volcano-tectonic depressions. The fourth is oriented towards the degradation of volcanoes by short- and long-term erosion processes. The fifth and final theme is twofold: first, we deal with geomorphic hazards on active and dormant volcanoes, along with five specific case studies of recent events; and lastly, we conclude with a chapter presenting a wide array of methods (morphometry, simulations of processes, structural geology, age determination, etc.) that are used to unravel processes on active and dormant volcanoes.
            In summary, our textbook aims to: (1) review the most recent research in geomorphology and physical volcanology, e.g. an improved classification and understanding of volcanic landforms, with respect to geodynamic settings, lithology, and climate; (2) update our knowledge of processes and rates of growth and destruction of volcanic landforms and landscapes by integrating recent results from the expanding sector of volcanology both in the field and in the laboratory; (3) consider how volcanic landforms, landscapes, and processes can be studied by reviewing classical and modern methods.  

In this way, we hope that our compilation, which provides a richly referenced and illustrated piece of work on volcano geomorphology, will be of interest to a broad audience. It is expected to be published later this year (2024).  

How to cite: Karátson, D. and Thouret, J.-C.: ‘Volcano Geomorphology: landforms, processes and hazards’  ̶̶ A new book in ‘Advances in Volcanology’ Collection, Springer Verlag , EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-11744, https://doi.org/10.5194/egusphere-egu24-11744, 2024.