TECTONIC CONTROLS ON FLUVIAL LANDSCAPE DEVELOPMENT IN CENTRAL-EASTERN PORTUGAL: INSIGHTS FROM LONG PROFILE TRIBUTARY STREAM ANALYSES

This study examines the long profiles of tributaries of the Tejo (Tagus) and Zêzere rivers in central eastern Portugal (West Iberia) in order to provide new insights into the patterns, timing and controls on drainage development during the Pleistocene to Holocene incision stage. The long profiles were extracted from lower order tributary streams associated with the trunk drainage of the Tejo River and one main tributary, the Zêzere River (Fig. 1). These streams flow through a landscape strongly influenced by variations in bedrock lithology (mainly granites and metasediments), fault structures delimiting crustal blocks with distinct uplift rates, and a base-level lowering history (tectonic uplift / eustatic). The long profiles of the tributaries of the Tejo and Zêzere rivers record a series of transient and permanent knickpoints. The permanent knickpoints have direct correlation with the bedrock strength, corresponding to the outcropping of very hard quartzites or to the transition from softer (slates/metagreywaques) to harder (granite) basement. The analyzed streams/rivers record also an older transient knickpoint/knickzone separating: a) an upstream relict graded profile, with lower steepness and higher concavity, that reflects a long period of quasi-equilibrium conditions reached after the beginning of the incision stage; and b) a downstream reach displaying a rejuvenated long profile, with steeper gradient and lower concavity, particularly for the final segment, which is often convex (Fig. 2). The rejuvenated reaches testify the upstream propagation of several incision waves that are the response of each stream to continuous or increasing crustal uplift and dominant periods of base-level lowering by the trunk drainages, coeval of low sea level conditions. The long profiles and their morphological configurations enabled spatial and relative temporal patterns of incision to be quantified for each individual tributary stream. The incision values of streams flowing in uplifted blocks of the Portuguese Central Range (PCR) (ca.380-280 m) indicate differential uplift and are higher than the incision values of streams flowing on the adjacent South Portugal planation surface – the Meseta (ca. 200 m). The normalized steepness index, calculated using the method of Wobus et al. (2006), proved to be sensitive to active tectonics, as lower ksn values were found in relict graded profiles of streams located in less uplifted blocks, (e.g. Sertã stream in the PCR), or in those flowing through tectonic depressions. Fig. 1 – Geological map of the study area. 1 – fluvial terraces (Pleistocene); 2 – sedimentary cover (Paleogene and Neogene); 3 – slates and metasandstones (Devonian); 4 – slates and quartzites (Silurian); 5 – quartzites (Ordovician); 6 – slates and metagreywackes (Precambrian to Cambrian); 7 – slates, metagreywackes and limestones (Precambrian); 8 – granites and ortogneisses; 9 – diorites and gabros; 10 - fault. SFf – Sobreira Formosa fault; Sf – Sertã fault; Pf – Ponsul fault; Gf – Grade fault. The differential uplift indicated by the distribution of the ksn values and by the fluvial incision was likely accumulated on a few major faults, as the Sobreira Formosa fault (SFf), thus corroborating the tectonic activity of these faults. Due to the fact that the relict graded profiles can be correlated with other geomorphic references documented in the study area, namely the T1 terrace of the Tagus River (with an age of ca. 1 Myr), the following incision rates can be estimated: a) for the studied streams located in uplifted blocks of the PCR, 0.38 m/kyr to 0.28 m/kyr; b) for the streams flowing on the South Portugal planation surface, 0.20 m/kyr. The differential uplift inferred between crustal blocks in the study area corroborates the neotectonic activity of the bordering faults, which has been proposed in previous studies based upon less robust data. Fig. 2 – Longitudinal profile of the Nisa stream a tributary of the Tejo River. Note the equilibrium relict profile upstream the older transient knickpoint (hatched line) and the downstream rejuvenated profile (continuous line). Legend: tKP – transient knickpoint; rKp – resistant knickpoint; Mt – schist and phyllite; Gr – granite; Hf – hornfels; Og – orthogneisse. In the inset Distance – Slope plots, fill circles correspond to the relict graded profile, crosses correspond to the rejuvenated profile located downstream the older transient knickpoint (tKP).

Lithosphere Thickness, Heat Flow and Moho Depth in the South of Portugal

Abstract In the last years, several models have been presented trying to obtain lithosphere and Moho thickness in the Iberian Peninsula, using data related to geoid elevation and topography, gravity, seismicity and thermal analysis. The results obtained show a decrease in the thickness of the crust and the lithosphere in the SW part of the Iberian Peninsula. Density anomalies in the crust are also referred. Data obtained in the region was collected and deviations from average values used were detected. In this work, models were made taking into account the specific characteristics of the region. Heat flow, thermal conductivity, heat production, topography, gravity, seismic and geological data available for the region, were used to adjust the model. The results show that this region is different from other parts of the Iberian Peninsula and a special attention must be given to it. This work shows the importance of trying to know and understand the thermal structure of the region.

Lithosphere Thickness, Heat Flow and Moho Depth in the South of Portugal

ABSTRACT In the last years, several models were presented trying to obtain lithosphere and Moho thickness in the Iberian Peninsula, using data related to geoid elevation and topography, gravity, seismicity and thermal analysis. The results obtained show a decrease in the thickness of the crust and the lithosphere in the SW part of the Iberian Peninsula. Density anomalies in the crust are also referred. The work I intend to present is related with the south of the Ossa Morena Zone, the South Portuguese Zone and the Algarve, in the south of Portugal. Data obtained in the region was collected and deviations from average values used were detected. Models were made taking into account the specific characteristics of the region. Heat flow, thermal conductivity, heat production, topography, gravity, seismic and geological data available for the region, are used to adapt the models. A special attention will be given to the spatial variation of heat flow values and to Moho depth in the region. The results show that this region is different from other parts of the Iberian Peninsula and a special attention must be given to it. The different values obtained using seismic, gravity, and geoid height data, and the results obtained with models using thermal data shows the importance of trying to know and understand the thermal structure of the regions. Problems related with the use of average values will be focused.

Tectonic and lithological controls on fluvial landscape development in central-eastern Portugal: Insights from long profile tributary stream analyses

This study examines the long profiles of tributaries of the Tagus and Zêzere rivers in Portugal (West Iberia) in order to provide new insights into patterns, timing, and controls on drainage development during the Quaternary incision stage. The studied streams are incised into a relict culminant fluvial surface, abandoned at the beginning of the incision stage. The streams flow through a landscape with bedrock variations in lithology (mainly granites and metasediments) and faulted blocks with distinct uplift rates. The long profiles of the analyzed streams record an older transitory knickpoint/knickzone separating (1) an upstream relict graded profile, with lower steepness and higher concavity, that reflects a long period of quasi-equilibrium conditions reached after the beginning of the incision stage, and (2) a downstream rejuvenated long profile, with steeper gradient and lower concavity, particularly for the final reach, which is often convex. The rejuvenated reaches testify to the upstream propagation of several incision waves, interpreted as the response of each stream to increasing crustal uplift and prolonged periods of base-level lowering by the trunk drainages, coeval with low sea level conditions. The morphological configurations of the long profiles enabled spatial and relative temporal patterns of incisions to be quantified. The incision values of streams flowing on the Portuguese Central Range (PCR; ca. 380–150 m) are variable but generally higher than the incision values of streams flowing on the adjacent South Portugal Planation Surface (SPPS; ca. 220–110 m), corroborating differential uplift of the PCR relative to the SPPS. Owing to the fact that the relict graded profiles can be correlated with the Tagus River T1 terrace (1.1–0.9 My) present in the study area, incision rates can be estimated (1) for the streams located in the PCR, 0.38–0.15 m/ky and (2) for the streams flowing on the SPPS, 0.22–0.12 m/ky. The differential uplift inferred in the study area supports the neotectonic activity of the bordering faults, as proposed in previous studies based upon other geological evidence.

Tectonic and lithological controls on fluvial landscape development in central-eastern Portugal: Insights from long profile tributary stream analyses

This study examines the long profiles of tributaries of the Tagus and Zêzere rivers in Portugal (West Iberia) in order to provide new insights into patterns, timing, and controls on drainage development during the Quaternary incision stage. The studied streams are incised into a relict culminant fluvial surface, abandoned at the beginning of the incision stage. The streams flow through a landscape with bedrock variations in lithology (mainly granites and metasediments) and faulted blocks with distinct uplift rates. The long profiles of the analyzed streams record an older transitory knickpoint/knickzone separating (1) an upstream relict graded profile, with lower steepness and higher concavity, that reflects a long period of quasi-equilibrium conditions reached after the beginning of the incision stage, and (2) a downstream rejuvenated long profile, with steeper gradient and lower concavity, particularly for the final reach, which is often convex. The rejuvenated reaches testify to the upstream propagation of several incision waves, interpreted as the response of each stream to increasing crustal uplift and prolonged periods of base-level lowering by the trunk drainages, coeval with low sea level conditions. The morphological configurations of the long profiles enabled spatial and relative temporal patterns of incisions to be quantified. The incision values of streams flowing on the Portuguese Central Range (PCR; ca. 380–150 m) are variable but generally higher than the incision values of streams flowing on the adjacent South Portugal Planation Surface (SPPS; ca. 220–110 m), corroborating differential uplift of the PCR relative to the SPPS. Owing to the fact that the relict graded profiles can be correlated with the Tagus River T1 terrace (1.1–0.9 My) present in the study area, incision rates can be estimated (1) for the streams located in the PCR, 0.38–0.15 m/ky and (2) for the streams flowing on the SPPS, 0.22–0.12 m/ky. The differential uplift inferred in the study area supports the neotectonic activity of the bordering faults, as proposed in previous studies based upon other geological evidence.

Extensional orogenic collapse captured by strike-slip tectonics:

The late Paleozoic collision between Gondwana and Laurussia resulted in the polyphase deformation and magmatism that characterizes the Iberian Massif of the Variscan orogen. In the Central Iberian Zone, initial con- tinental thickening (D1; folding and thrusting) was followed by extensional orogenic collapse (D2) responsible for the exhumation of high-grade rocks coeval to the emplacement of granitoids. This study presents a tectonometamorphic analysis of the Trancoso-Pinhel region (Central Iberian Zone) to ex- plain the processes in place during the transition froman extension-dominated state (D2) to a compression-dom- inated one (D3).Wereveal the existence of low-dipping D2 extensional structures later affected by several pulses of subhorizontal shortening, each of them typified by upright folds and strike-slip shearing (D3, D4 and D5, as identified by superimposition of structures). The D2 Pinhel extensional shear zone separates a low-grade domain from an underlying high-grade domain, and it contributed to the thermal reequilibration of the orogen by facil- itating heat advection from lower parts of the crust, crustal thinning, decompression melting, and magma intru- sion. Progressive lessening of the gravitational disequilibrium carried out by this D2 shear zone led to a switch from subhorizontal extension to compression and the eventual cessation and capture of the Pinhel shear zone by strike-slip tectonics during renewed crustal shortening. High-grade domains of the Pinhel shear zone were folded together with low-grade domains to define the current upright folded structure of the Trancoso-Pinhel re- gion, the D3 Tamames-Marofa-Sátão synform. Newdating of syn-orogenic granitoids (SHRIMP U\\Pb zircon dat- ing) intruding the Pinhel shear zone, together with the already published ages of early extensional fabrics constrain the functioning of this shear zone to ca. 331–311 Ma, with maximum tectonomagmatic activity at ca. 321–317 Ma. The capture and apparent cessation of movement of the Pinhel shear zone occurred at ca. 317– 311 Ma.

Carbonates and sulfides in orogenic garnet-bearing peridotites from the Ulten Zone as witnesses for Carbon and Sulfur mobility in orogenic settings

Convergent plate boundaries are sites of sustained chemical exchanges between the Earth’s surface and deep geochemical reservoirs, playing a major role in the global cycle of carbon and sulfur. However, carbon and sulfur recycling processes continue to be hotly debated. A critical gap in the knowledge of the whole subduction factory is given by the limited accessibility to the upper mantle residing above the subducting plate, the so-called mantle wedge. This thesis investigates the carbonate and sulfide metasomatism taking place during the whole metamorphic evolution of a mantle wedge involved in the Variscan continental collision. We integrate different detailed geochemical and petrological techniques to orogenic carbonated spinel and garnet peridotites from the Ulten Zone of the Eastern Italian Alps. Our data show that the Ulten Zone peridotite experienced multiple stages of addition and removal of carbon and sulfur throughout its metamorphic evolution, as follows: (1) The Variscan lithospheric mantle was initially depleted and sulfide-poor. It subsequently inherited a sulfur and carbon component during an early metasomatic stage, when hot, H2S-CO2-bearing melts leaving a subduction-modified source percolated the overlying spinel-facies peridotite in the mantle wedge; (2) Under peak eclogite-facies P-T conditions, pervasive carbonation and sulfidation occurred. Heterogeneous melt and fluid sources variably enriched in carbon, isotopically heavy sulfur and radiogenic Sr were involved; (3) Shortly after the attainment of peak-P conditions, peridotite bodies were incorporated in a tectonic mélange with the neighboring gneisses. Here, the Ulten Zone peridotite was exposed to channelized infiltration of hybridized C-O-H fluids that promoted the formation of veinlets of carbonates locally associated with sulfide grains. (4) Upon late retrogression, infiltration of serpentinizing fluids promoted C and S remobilization at shallow crustal levels.

The Moho reflectivity from ambient seismic noise autocorrelations beneath the Western Alps

The lower crustal structure beneath the Western Alps -- including the Moho -- bears the signature of past and present geodynamic processes. It has been the subject of many studies until now. However, its current knowledge still leaves significant open questions. In order to derive new information, independent from previous determinations, here I wish to address this topic using a different method --- ambient seismic noise autocorrelation --- that is for the first time applied to reveal Moho depth in the Western Alps. Moho reflections are identified by picking reflectivity changes in ambient seismic noise autocorrelations. The seismic data is retrieved from more than 200 broadband seismic stations, from the China--Italy--France Alps (CIFALPS) linear seismic network, and from a subset of the AlpArray Seismic Network (AASN). The automatically-picked reflectivity changes along the CIFALPS transect in the southwestern Alps show the best results in the 0.5--1 Hz frequency band. The autocorrelation reflectivity profile of the CIFALPS transect shows a steeper subduction profile,~55 to ~70 km, of the European Plate underneath the Adriatic Plate. The dense spacing of the CIFALPS network facilitates the detection of lateral continuity of crustal structure, and of the Ivrea mantle wedge reaching shallow crustal depths in the southwestern Alps. The data of the AASN stations are filtered in the 0.4--1 and 0.5--1 Hz frequency bands. Although the majority of the stations give the same Moho depth for the different frequency bands, the few stations with different Moho depths shows the care that has to be taken when choosing the frequency band for filtering the autocorrelation stacks. The new Moho depth maps by using the AASN stations are a compilation of the first and second picked reflectivity changes. The results show the complex crust-mantle structure with clear differences between the northwestern and southwestern Alps.

Early-stage deformation localisation in thrust systems: new perspectives from carbonates in the Italian Southern Alps and Oman Mountains

This thesis investigates mechanisms and boundary conditions that steer the early localisation of deformation and strain in carbonate multilayers involved in thrust systems, under shallow and mid-crustal conditions. Much is already understood about deformation localisation, but some key points remain loosely constrained. They encompass i) the understanding of which structural domains can preserve evidence of early stages of tectonic shortening, ii) the recognition of which mechanisms assist deformation during these stages and iii) the identification of parameters that actually steer the beginning of localisation. To clarify these points, the thesis presents the results of an integrated, multiscale and multi-technique structural study that relied on field and laboratory data to analyse the structural, architectural, mineralogical and geochemical features that govern deformation during compressional tectonics. By focusing on two case studies, the Eastern Southern Alps (northern Italy), where deformation is mainly brittle, and the Oman Mountains (northeastern Oman), where ductile deformation dominates, the thesis shows that the deformation localisation is steered by several mechanisms that mutually interact at different stages during compression. At shallow crustal conditions, derived conceptual and numerical models show that both inherited (e.g., stratigraphic) and acquired (e.g., structural) features play a key role in steering deformation and differentiating the seismic behaviour of the multilayer succession. At the same time, at deeper crustal conditions, strain localises in narrow domains in which fluids, temperature, shear strain and pressure act together during the development of the internal fabric and the chemical composition of mylonitic shear zones, in which localisation took place under high-pressure (HP) and low-temperature (LT) conditions. In particular, results indicate that those shear zones acted as “sheltering structural capsules” in which peculiar processes can happen and where the results of these processes can be successively preserved even over hundreds of millions of years.

Assessment of sea-level variability for the Emilia-Romagna coastal area in the framework of the Mediterranean Sea

Sea–level change is one of the ocean characteristics closely connected to climate change. Understanding its variation is essential since a large portion of the world’s population is located in low–lying locations. Two main techniques are employed to measure sea level: satellite altimetry and tide gauges. Satellite altimetry monitors sea–level relative to a geocentric reference, is unaffected by crustal processes and covers nearly the entire surface of the oceans since 1993. Conversely, tide gauges measure sea level at specific coastal locations and relative to a local ground benchmark, therefore are impacted by vertical land movements. In this study, the linear and non–linear geocentric and relative sea–level trends along the Emilia–Romagna coast (Northern Italy) have been analyzed over different periods. In order to assess the local sea–level variability, data from satellite altimetry and tide gauges have been compared over a common time interval (1993–2019), hence disentangling the contribute of vertical land movements. Non–linearity has been also evaluated at the broader scale of the Mediterranean Sea, in order to depict the main variability in geocentric sea–level trends from regional to sub–basin scale. Furthermore, the anthropogenic and natural influence at the river basin scale has been addressed, in order to shed light on the factors inducing the drastic reduction of riverine sediment supply to the Emilia–Romagna coast over the period 1920–2020. The findings of this analysis indicate that the sediment delivery reduction to the coast by rivers has been driven by several anthropogenic processes, acting on various spatiotemporal scales. Moreover, the local absolute sea–level trend is far from linear and appear "contaminated" by the presence of natural oscillations that act at the multi–decadal, quasi–decadal and inter–annual scale, mainly driven by both large–scale climatic modes and variations in local oceanography.