Longitudinal River Profiles Research Articles

Fluvial inverse modeling for inferring the timing of Quaternary uplift in the Simbruini range (Central Apennines, Italy)

AbstractThe timing of Quaternary uplift evolution of the Simbruini range (Central Apennines, Italy) is poorly known due to the lack of reliable chronological constraints of the post‐orogenic continental clastic units, deposited after the upper Messinian thrust‐top facies, and tectonic events. In this regard, we identified the main geomorphic markers, including several levels of continental units, whose plano‐altimetric projection along the river longitudinal profile has been correlated with the main non‐lithological knickpoints. Furthermore, we inferred the uplift history of the range by applying the inverse modeling of the river longitudinal profiles. Assuming a block uplift model, the drainage network cutting the Simbruini range recorded 2.4 ± 0.5 Ma of tectonic history, characterized by variable base‐level fall rates (corresponding to variable uplift rates).

Relationship Between Dams, Knickpoints and the Longitudinal Profile of the Upper Indus River

Mass movements in mountainous areas are capable of damming rivers and can have a lasting effect on the river longitudinal profile. The long profile is commonly used to retrieve regional tectonic information, but how much dams may compromise geomorphometry-based tectonic analysis has not been systematically researched. In this study, we investigate the relationship between river dams and the longitudinal profile of the upper Indus River basin, based on interpretation and analysis of remote sensing imagery and digital elevation models (DEMs) and local field work. We identified 178 landslide, glacier and debris flow dams. Using TopoToolbox, we automatically extracted the river longitudinal profile from the 30 m SRTM DEM, determined the location of convex knickpoints and calculated the channel steepness index. One hundred and two knickpoints were detected with heights above 148 m, of which 55 were related to dams. There is good spatial correspondence between dams, convexities in the river longitudinal profile and relatively high steepness index. Different dam types have different impacts on the river profile; on the upper Indus, debris flow dams have a greater impact than landslide and glacier dams and can form knickpoints of up to 900 m. Therefore, dams may have a significant influence on the river longitudinal profile, knickpoints and steepness index, and should be considered when extracting information on regional tectonics using these indices.

Research on the Longitudinal Section of River Restoration Using Probabilistic Theory.

Since the 1960s, many rivers have been destroyed as a consequence of the process of rapid urbanization. As accurate figures are important to repair rivers, there have been many research reports on methods to obtain the exact river slope and elevation. Until now, many research efforts have analyzed the river using measured river topographic factors, but when the flow velocity changes rapidly, such as during a flood, surveying is not easy; and due to cost, frequent measurements are difficult. Previous research has focused on the cross section of the river, so the information on the river longitudinal profile is insufficient. In this research, using informational entropy theory, equations are presented that can calculate the average river slope, river slope, and river longitudinal elevation for a river basin in real time. The applicability was analyzed through a comparison with the measured data of river characteristic factors obtained from the river plan. The parameters were calculated using informational entropy theory and nonlinear regression analysis using actual data, and then the longitudinal elevation entropy equation for each river and the average river slope were calculated. As a result of analyzing the applicability of the equations presented in this study by R2 and Root Mean Square Error, all R2 values were over 0.80, while RMSE values were analyzed to be between 0.54 and 2.79. Valid results can be obtained by calculating river characteristic factors.

Author Correction: Global dominance of tectonics over climate in shaping river longitudinal profiles

Author Correction: Global dominance of tectonics over climate in shaping river longitudinal profiles

Tectonically driven drainage reorganization in the Eastern Cordillera, Colombia

High-elevation plateaus that are positioned in between topographic barriers are common orogenic features in the South American continent, formed under a range of evolving environmental conditions. For example, in the central Andes (Bolivia-Argentina), the Puna-Altiplano is arid and endorheic with a poorly developed fluvial system, while in the northern Andes (Colombia) the Chiquinquirà and Tunja highlands are characterized by a humid equatorial exorheic fluvial system. In addition to a plateau-like low-relief surface at 2500 m, the landscape of the northern Eastern Cordillera and Santander Massif (northern Colombia) displays a lower elevation (~1500 m) low-relief landscape (Mesas) comprising river captures, windgaps, and a disconnected alluvial fan that collectively record a transient state. This configuration has been achieved through a combination of compressive deformation and sub-crustal processes. The compressive shortening started to occur in the Paleogene and is still active, whereas regional surface uplift related to slab flattening and mantle wedge hydration started in the Late Miocene/Pliocene. To disentangle the crustal vs sub-crustal forcing and to investigate the relative timing of drainage network evolution we combine the analysis of topography, hydrography (river longitudinal profiles, morphometric parameters, drainage divide stability), knickpoint migration (celerity model), paleo-longitudinal profile modeling, satellite images, and field observations.In particular, we show that during the development of the low-relief Mesas landscape the older Chiquinquirà highland was a closed drainage and that the lower portion of the Suárez River flowed northward into the Bucaramanga depression forced by the Los Cobardes Anticline topographic barrier. The Suárez River collected waters from the southern Santander Massif and the upper reach of the Chicamocha River, which was draining the Tunja highland. An abandoned windgap deposit on the eastern edge of the Mesa de Barichara suggests that the lower portion of the Chicamocha River was not yet formed. Subsequent to the Chiquinquirà highland drainage opening, two main tributaries of the Magdalena River, the Lebrija and Sogamoso, captured the Suárez River in a short temporal sequence. A knickpoint celerity model allows us to date the Lebrija capture of the Bucaramanga depression at ~260–270 ka and the subsequent Sogamoso capture at 190–220 ka. Only during this final stage, the lowermost Chicamocha River section formed and the drainage network developed to its present configuration.Finally, we suggest that the early Cenozoic rift inversion has controlled the drainage network pattern and the late Miocene sub-crustal-induced surface uplift has driven the main fluvial network reorganization.

A Geomorphic Examination of the Calabrian Forearc Translation

AbstractThe evolution of subduction forearcs as archived in the tectonic and geomorphic record serves to constrain geodynamic models of many convergent plate margins. This paper presents geomorphic stream channel metrics and river longitudinal profile modeling to assess geodynamic models describing the late Cenozoic horizontal translation of a Mediterranean forearc. The northern flank of Sicily is drained by roughly parallel north‐flowing streams arranged across the hypothesized west‐to‐east migration path of the Calabrian forearc high. These streams present a snapshot of current geomorphic processes driven by spatiotemporal variations in crustal deformation, including river incision and landscape steepness. A stream power rule‐based numerical modeling of the fluvial response time and knickpoint migration provides an independent prediction for the history of rock uplift associated with the translating forearc over the past ∼4 Myr. The results serve as a test for two, crustal wave versus transform‐fault, end‐member geodynamic models for the Calabrian forearc evolution, suggesting that our modeling supports the former. The river profile modeling shows that northern Sicily rivers experienced two pulses of time‐transgressive rock uplift that progressively sweeps west to east. This analysis suggests the most rapid contemporary rock uplift rates are focused in NE‐Sicily, where the forearc is currently located and where Pleistocene marine terraces indicate the highest rock uplift rates. These results are consistent with the west‐to‐east passage of a forearc high embedded in the subduction wedge, followed by a dynamic wave of rock uplift in its wake similar to what has been proposed for slab windows in other plate boundary settings.

Cenozoic Dynamic Topography of Madagascar

AbstractIt has been proposed that Oligo‐Miocene regional uplift of Madagascar was generated and is maintained by mantle dynamical processes. Expressions of regional uplift include flat‐lying Upper Cretaceous‐Paleogene marine limestones that crop out at elevations of hundreds of meters along the western seaboard and emergent Quaternary coral‐rich terraces that rim the coastline. Here, we explore the history of subcrustal topographic support through a combined analysis of four sets of observational constraints. First, we exploit published receiver function estimates of crustal thickness and spectral admittance between gravity and topography. An admittance value of ∼+40 ± 10 mGal km−1 at wavelengths >500 km implies that ∼1 km of topography is supported by subcrustal processes. Secondly, new apatite fission‐track and helium measurements from 18 basement samples are inverted, constraining temperature and denudation histories. Results suggest that 0.5–1.6 km of regional uplift occurred after ∼30 Ma. Thirdly, we calculate a history of regional uplift by minimizing the misfit between observed and calculated longitudinal river profiles. Results suggest that topography was generated during Neogene times. Finally, inverse modeling of rare earth element concentrations in Neogene mafic rocks indicates that melting of the asthenospheric source occurred at depths of ≤65 km with potential temperatures of 1300–1370 °C. Melting occurred at higher temperatures beneath Réunion Island and northern Madagascar and at lower temperatures beneath the Comores and southern Madagascar. These inferences are consistent with shear wave velocities obtained from tomographic models. We conclude that Madagascar is underlain by thinned lithospheric mantle and that a thermal anomaly lies within an asthenospheric layer beneath northern Madagascar.

Scale‐Dependent Contributors to River Profile Geometry

AbstractA range of complex hydraulic and geomorphic processes shape terrestrial landscapes. It remains unclear how these processes act to generate observed drainage networks across scales of interest. To address this issue, we transform observed and synthetic longitudinal river profiles into the spectral domain with a view to interrogating the different scales at which fluvial landscapes are generated. North American river profiles are characterized by red noise (i.e., spectral power, ϕ ∝ k−2, where k is wave number) at wavelengths >100 km and pink noise (ϕ ∝ k−1) at shorter wavelengths. This observation suggests that river profile geometries are scale‐dependent and using small‐scale observations to develop a general understanding of large‐scale landscape evolution is not straightforward. At wavelengths >100 km, river profile geometries appear to be controlled by smoothly varying patterns of regional uplift and slope‐dependent incision. Landscape simulations, based upon stream power that are externally forced by regional uplift do not exhibit a spectral transition from red to pink noise because these simulations do not incorporate heterogeneous erodibility. Spectral analysis of erodibility extracted from patterns of lithologic variation along river profiles suggests that the missing spectral transition is accounted for by heterogeneous substrates, which are characterized by white or blue noise (ϕ ∝ k0 or k1). Our results have implications for the way by which rivers record large‐scale tectonic forcing while incising through complex lithologic patterns.

Along-Strike Rapid Structural and Geomorphic Transition From Transpression to Strike-Slip to Transtension Related to Active Microplate Rotation, Papua New Guinea

The area of southeastern Papua New Guinea includes three active microplates – the Trobriand, Woodlark, and Solomon Sea plates – that are being deformed by regional convergence between the much larger Pacific and Australian Plates. The landward extent of the plate boundary between the Trobriand and Australian Plates corresponds to the Owen-Stanley Fault Zone (OSFZ), an onland and continuous 510 km-long left-lateral strike-slip fault that forms a linear, intermontane valley within the elongate Owen-Stanley Range (OSR) and continues as a 250 km-long low-angle normal fault along the margins of Goodenough and Woodlark basins. GPS geodesy reveals that the Trobriand microplate has undergone rapid counter-clockwise rotation since the Late Miocene (8.4 Ma) and that this rotation about a nearby pole of rotation predicts transpressional deformation along the 250 km-long northwestern segment of the OSFZ, strike-slip motion along a 100 km-long central segment, and transtension along the 270 km-long ESE-trending southeastern segment of OSFZ. In order to illustrate the along-strike variations in neotectonic uplift resulting from the changing structure of the OSFZ, we delineated 3903 river segments in the northeastern side of the OSR drainage divide and derived river longitudinal profiles along each river segment. Normalized steepness indices (ksn) and knickpoint clusters are the highest and most concentrated, respectively, for the northwestern transpressional segment of the OSR, moderately high and concentrated along the southeastern segment of the OSR, and the lowest and least concentrated along the central strike-slip segment. These geomorphological indices indicate that most of the plate boundary uplift occurs along the transpressional and transtensional segments that are connected by the central strike-slip zone. Within this overall pattern of structural variation, abrupt changes in the azimuth of the OSFZ create more localized anomalies in the geomorphological indices.

Global dominance of tectonics over climate in shaping river longitudinal profiles

River networks are striking features engraved into Earths surface, shaped by uplift and erosion under the joint influence of climate and tectonics. How a river descends along its course – its longitudinal profile – varies greatly from one basin to the next, reflecting the interplay between uplift and erosional processes. It has recently been argued that climatic aridity should be a first-order control on river profile concavity, but the importance of climate relative to other factors has not been tested at global scale. Here we show, using recent global datasets of river profiles and tectonic activity, that tectonics is much more strongly expressed than climate in global patterns in river profile concavity. River profiles tend to be more strongly concave in tectonically active regions along plate boundaries, reflecting tectonically induced spatial variations in uplift rates. Rank correlations between river profile concavity and global tectonic proxies (basin slope and two indices of seismic risk) are much stronger than those between river concavity and precipitation, potential evapotranspiration, or aridity. These results show that tectonics, and not climate, exerts first-order control on the shape of river longitudinal profiles globally.

Evaluating the Vertical Accuracy of DEM Generated from ZiYuan-3 Stereo Images in Understanding the Tectonic Morphology of the Qianhe Basin, China

Currently available high-resolution digital elevation model (DEM) is not particularly useful to geologists for understanding the long-term changes in fluvial landforms induced by tectonic uplift, although DEMs that are generated from satellite stereo images such as the ZiYuan-3 (ZY3) satellite include characteristics with significant coverage and rapid acquisition. Since an ongoing analysis of fluvial systems is lacking, the ZY3 DEM was generated from block adjustment to describe the mountainous area of the Qianhe Basin that have been induced by tectonic uplift. Moreover, we evaluated the overall elevation difference in ZY3 DEM, Shuttle Radar Topography Mission (1″ × 1″) (SRTM1), and Reflection Radiometer Global Digital Elevation Model (ASTER GDEM) by using the Ice Cloud and Land Elevation Satellite/Geoscience Laser Altimeter (ICESat/GLAH14) point cloud and a DEM of 1:50,000 scale. The values of the root mean square error (RMSE) of the elevation difference for ZY3 DEM were 9.31 and 9.71 m, respectively, and are in good agreement with SRTM1. The river long profiles and terrace heights were also extracted to compare the differences in channel steepness and the incision rates with SRTM1 and ASTER GDEM. Our results prove that ZY3 DEM would be a good alternative to SRTM1 in achieving the 1:50,000 scale for DEM products in China, while ASTER GDEM is unsuitable for extracting river longitudinal profiles. In addition, the northern and southern river incision rates were estimated using the ages and heights of river terraces, demonstrating a range from 0.12–0.45 to 0.10–0.33 m/kyr, respectively. Collectively, these findings suggest that ZY3 DEM is capable of estimating tectonic geomorphological features and has the potential for analyzing the continuous evolutionary response of a landscape to changes in climate and tectonics.

Evaluating the performance of multisource digital elevation models using morphometric parameters and field survey data over the mountainous landscapes of northwest Himalaya, India

In this study, seven morphometric characteristics were assessed from three freely available multisource DEMs and correlated with extensive field observations to assess the best-fit DEM over the mountainous landscapes of Kashmir Valley, northwest Himalaya, India. The morphometric parameters were assessed using 30 m ASTER GDEM v2 (ASTERDEM), SRTMDEM, and CartoDEM v3 R1 (CARTODEM) in the Ferozpora watershed of the Jhelum Basin, Kashmir Himalaya. Our findings indicated the closeness of mean bifurcation ratio (Rbm) values between manual digitization (Rbm = 4.12) and CARTODEM (Rbm = 4.07). The GPS-based values of drainage basin asymmetry, stream gradient index, basin relief, and longitudinal river profile also indicated a strong resemblance with CARTODEM-derived values. It is hence concluded that CARTODEM is the best-fit DEM representing the topographic and morphometric characteristics over the Kashmir region, however, the results from this analysis need to be tested and validated over larger geographic domains with contrasting lithological and topographic characteristics.

Lithostructural controls in Douala-Buea Region landscape (SW Cameroon margin): Insights from morphometric analysis

The Douala-Buea Region (DBR; 4°–5°N, 9°–10°30′E), a portion of the southern Cameroon passive margin, developed on the Proterozoic granito-gneissic basement and the Cretaceous to Recent volcano-sedimentary cover, was explored in order to evaluate the lithology and tectonics controls on its landscape using the qualitative and quantitative geomorphological analyses based on 30 m resolution Digital Elevation Models (DEM) and field evidences. This area displays an uneven and complex landscape differentiated into four morphological units: Douala, Nkondjock, Kumba and Buea Units. The Douala Unit corresponds to the Coastal Plain and is extended on the Douala-Kribi/Campo sedimentary basin. The Nkondjock Unit represents a middle plateau developed on a Proterozoic granito-gneissic basement. The Kumba Unit matches with the southernmost volcanic Western Cameroon High Plateau. The Buea Unit corresponds to the Mount Cameroon stratovolcano. The relief analysis through mountain-front sinuosity (Smf, 1.19–1.43) and relative relief ratio (Rr, 0.10–2.31) vary from one unit to another. The drainage network analysis including drainage patterns, stepped longitudinal river profiles, the associated morphometric indexes, notably the concavity index (IC, 0.56–0.88), deviation (D, 267–2912 m), basin asymmetry factor (A f, 6.37–42.12), hypsometric integral (HI, 0.10–0.29), valley floor width to valley height ratio (Vf, 0.36–2.75), elongation ratio (Re, 0.30–0.50), watershed slope (S, 6.88–88.88) and stream length gradient index (SL, 0–3332) vary from one basin to another too. These results also show the heterogeneous and uneven DBR landscape controlled by its lithology and active tectonics. The drainage system reveals concave to concavo-convex longitudinal river profile shapes with knick points, asymmetric basins, and fault reactivations, and highlights the tectonics activity that used to occur through the time and their control on the morphology. Together, they document the lithostructural controls on the DBR landscape evolution, extendable in the whole SW Cameroon margin.

Use of HAND terrain descriptor for estimating flood-prone areas in river basins

The flood hazard mapping in a river basin is crucial for flooding risk management, mitigation strategies, and flood forecasting and warning systems, among other benefits. One approach for this mapping is based on the HAND (Height Above Nearest Drainage) terrain descriptor, directly derived from the Digital Elevation Model (DEM), in which each pixel represents the elevation difference of this point in relation to the river drainage network to which it is connected. Considering the Mamanguape river basin (3,522.7 km²; state of Paraíba, Brazil) as the study location, the present research applied this method and verified it as for five aspects: consideration of a spatially variable minimum drainage area for denoting the river drainage initiation; the impact of considering a depressionless DEM; evaluation of hydrostatic condition; effect of incorporating an existing river vector network; and comparative analysis of basin morphology regarding longitudinal river profiles. According to the results, adopting a uniform minimum drainage area for the river network initiation is a simplification that should be avoided, using a spatially variable approach, which influences the amount and spatial distribution of flooded areas. Additionally, considering the depressionless DEM leads to higher values of HAND and to a smaller flooded area (difference ranging between 3% and 99%), when compared with the use of DEM with depression, despite 3.1% of the pixels representing depressions. The use of the depressionless DEM is recommended, whereas the DEM pre-processing by incorporating a vector network (stream burning) generates dubious results regarding the relation between HAND and the morphological pattern presented in the DEM. Moreover, the estimation of flooded areas based on HAND does not guarantee the hydrostatic condition, but this disagreement comprises a negligible area for practical purposes.

Tectonic control on drainage network evolution and evidence of neotectonic activities in of the Wan River sub-basin, Central India

The study area lies in the Northeastern part of the Gavilgrah fault vicinity, which is one of the essential tectonic elements of the Son-Narmada-Tapti lineament, trending ENE-WSE direction. In this study field, data and satellite data observation are used for drainage characteristics of Wan river sub-basins. The course of Wan River is mostly controlled by the Gavilgrah fault and associated tectonic elements. The Wan River, which crosses the North’s fault, exhibits two river terraces that rise and continue downstream, where the river flows on the down-thrown block. The drainage characteristics of the study area were investigated by using the LISS- III satellite data. The study of neotectonism is based on field evidence and drainage pattern. The nature of the drainage pattern in the study area shows dendritic, rectangular, and parallel together and shows evidence of neotectonic activity in the study area. The longitudinal river profile, the shape of the basin, gradient fluctuations, and valley incision, and their result also indicate neotectonic activities in the study area.

Active Deformation and Relief Evolution in the Western Lurestan Region of the Zagros Mountain Belt: New Insights From Tectonic Geomorphology Analysis and Finite Element Modeling

AbstractTo unravel how and where coseismic and interseismic deformation impacts the spatial and temporal patterns of rock uplift of the Lurestan sector of the Zagros Mountains, we performed an investigation of the large‐scale features of topography and river network coupled with 2‐D finite element modeling. Geomorphological analysis and constraints from parameters such as elevation, local relief, normalized channel steepness index (ksn), river longitudinal profiles, and transformed river profiles (chi plots) were used to unravel the time‐space distribution of vertical motions. Whereas the much longer timescale over which topography grows and/or rivers respond to tectonic or climatic perturbations with respect to even multiple seismic cycles, the outputs of the finite element model yield fundamental information on the source of the late part of the spatiotemporal evolution of surface uplift recorded by the geomorphology. Model outputs shed new light into the processes controlling relief evolution in an actively growing mountain belt underlain by a major blind thrust. The outputs illustrate how coseismic slip controls localized uplift of a prominent topographic feature—the Mountain Front Flexure—located above the main upper crustal ramp of the principal basement thrust fault of the region, while continuous displacement along the deeper, aseismic portion of the same basement fault controls generalized uplift of the whole crustal block located farther to the NE, in the interior of the orogen.

Transport-limited fluvial erosion – simple formulation and efficient numerical treatment

Abstract. Most of the recent studies modeling fluvial erosion in the context of tectonic geomorphology focus on the detachment-limited regime. One reason for this simplification is the simple relationship of the constitutive law used here – often called stream-power law – to empirical results on longitudinal river profiles. Another no less important reason lies in the numerical effort that is much higher for transport-limited models than for detachment-limited models. This study proposes a formulation of transport-limited erosion where the relationship to empirical results on river profiles is almost as simple as it is for the stream-power law. As a central point, a direct solver for the fully implicit scheme is presented. This solver requires no iteration for the linear version of the model, allows for arbitrarily large time increments, and is almost as efficient as the established implicit solver for detachment-limited erosion. The numerical scheme can also be applied to linear hybrid models that cover the range between the two end-members of detachment-limited and transport-limited erosion.

Active Tectonics of the Frontal Himalayas: An Example from the Manzai Ranges in the Recess Setting, Western Pakistan

The Himalayan main frontal thrust (MFT) accommodates most of the present-day Indo–Asia convergence with related periodic earthquakes. The seismicity and deformation mechanism varies considerably across the frontal Himalayas. We mapped a segment (Manzai Ranges) of the MFT at the western margin of the Himalayas and analyzed its deformation mechanism and active tectonics using geomorphic indices and the Interferometric Synthetic Aperture Radar (InSAR) Small Baseline Subset (SBAS) technique. Two frontal thrust faults (Khirgi and Jandola) were mapped using Sentinel-2B band ratios in the study area. Water gaps were present in the form of deflected streams at the tip of the growing anticlines. The C-band RADAR interferometry (Sentinel-1A) showed an average uplift of 5–9 mm/year in the satellite line of sight (LOS) from May 2018 to October 2019. The velocity profiles show an uplift variation across the anticlines and may be related to the displacement transfer from the zone of compression in the Manzai Ranges to the zone of transpression in the Pezu–Bhittani Ranges. Four types of morphometric analyses were carried out to assess the relative tectonic activity, namely mountain front sinuosity index (Smf), valley floor width to height ratio (Vf), normalized longitudinal river profile, and normalized channel steepness index (Ksn). The landscape response to active tectonics in the study area was recorded as a deep fluvial incision in V-shaped valleys, convex river profiles, topographic breaks as knickpoints, and a high Ksn index. The geomorphic parameters show a relative increase in tectonic uplift and deformation from the Kundi anticline to the Khirgi and Manzai anticline. We concluded that the frontal structures in the western Himalayas are still going through an active phase of deformation and landscape development with both seismic and aseismic creep.

Assessment of tectonic control on the development of low mountains moderate relief in the Outer Carpathians (Southern Poland)

Inherited tectonic structures, ongoing tectonic deformation, and variations in relative rock uplift rates play an important role in conditioning the processes of relief development. Their influence among other factors, such as climate and lithology, can be quantified using landscape analysis, and geomorphometric indices, in particular. The usage of landscape analysis in recent years is increasing systematically due to the constant improvement of the digital elevation models and GIS software that significantly facilitate this approach. In this study, we aim to recognize the influence of tectonic structures and processes on relief development in the low mountains with moderate relief of the Soła River catchment in the Western Outer Carpathians. To this end, we calculated geomorphometric indices (river longitudinal profile, stream-length gradient index, minimum bulk erosion, relief ratio, circulatory ratio, elongation ratio, and hypsometric integral) for the Sola River and its 47 sub-catchments using a 25-m spatial resolution Digital Terrain Elevation Data Level 2. Additionally, we identified lineaments and knickpoints and correlated the computed results with local and regional fault networks, variations in lithology, and climate fluctuations. Obtained results indicate a significant impact of inherited tectonic structures on the relief development of the Soła River catchment, i.e., directions of principal ridges and valleys follow the orientation of main folds and faults recorded in this area. Anomalously high values of minimum bulk erosion, river gradient, and stream-length gradient index allowed us to define two areas with higher relative uplift rates: 1) the Sola Gorge and 2) the Beskid Żywiecki Mts. Polish Outer Carpathians are generally considered as an area of low strain rate and low seismic activity. However, the possibility of neotectonic processes should be considered in geohazard estimations. Observed bends in the direction of river valleys that do not correspond with changes in lithology could be related to active strike-slip faults. These are probably the reactivated basement structures, copied in the thin-skinned nappe cover, as a result of the accommodation of the Mur-Žilina Fault Zone resulting from the tectonic push of the Alcapa (Alpine-Carpathian-Pannonian) microplate against the European plate. Thus, the role of recent tectonic activity in relief development of the Sola River catchment even though appears to be subsidiary at the most, should not be excluded.

Calibrating fluvial erosion laws and quantifying river response to faulting in Sardinia, Italy

It is now widely accepted that rivers modify their erosion rates in response to variable rock uplift rates, resulting in changes in channel slope that propagate upstream through time. Therefore, present-day river morphology may contain a record of tectonic history. The simple stream power incision model can, in principle, be used to quantify past uplift rates over a variety of spatial and temporal scales. Nonetheless, the erosional model's exponents of area and slope (m and n respectively) and ‘bedrock erodibility’ (k) remain poorly constrained. In this paper, we will use a geologically and geomorphically well constrained Plio-Pleistocene volcanic landscape in central Sardinia, Italy, to calibrate the stream power erosion equation and to investigate the slip rate of faults that have been seismically quiescent in the historic past. By analysing digital elevation models, geological maps and Landsat imagery, we have identified the geomorphic expression of several volcanic features (eruption centres and basaltic lava flows) and three normal faults with 6 to 8 km fault traces within the outcrop. Downstream, river longitudinal profiles show a similar transient response to relative base level fall, probably as a result of relief inversion at the edge of the volcanic outcrop. From measurements of incision, local slope and upstream catchment area across eight different rivers, we calculate n ≈ 1, m = 0.50 ± 0.02 and, using a landscape age from literature of 2.7 Ma, bedrock erodibility k = 0.10 ± 0.04 m(1−2m) Myr−1. There are also knickpoints on rivers upstream of two normal faults, and we used numerical inverse modelling of the longitudinal profiles to predict the slip rate of these faults since 2.7 Ma. The results from the inverse model show that the erosional parameter values derived in this study can produce theoretical longitudinal profiles that closely resemble observed river profiles upstream of the faults. The lowest misfit theoretical longitudinal profiles were generated by a model of temporally discontinuous footwall uplift with consistently low throw rates (<0.1 mm yr−1). The predicted footwall uplift history is similar for the two faults, both showing periods of fault slip and no fault movement since 2.7 Ma.