Mountain building reorganizes drainage networks, influencing riverine biodiversity. Northern Italy offers a natural experiment in the impact of tectonic and geomorphic processes on aquatic species distributions. We combined geomorphic analysis with environmental DNA from rivers to assess the influence of tectonically driven drainage reorganization on genetic diversity, targeting an endemic fish species, Telestes muticellus (A. Risso, 1827). In the Northern Apennines, horizontal shortening and topographic advection in an orogenic wedge have been hypothesized as leading to river capture and drainage divide migration. In addition, slab rollback has produced a spatial transition from contraction to extension, which is more pronounced from north to south, with normal faulting producing range-parallel drainage only in the southern regions. In contrast, the adjacent Ligurian Alps are a remnant of the Alpine orogen with little modern deformation. We found distinct zones of geomorphic characteristics from north to south, including divide asymmetry and frequency of range-parallel drainage. Analysis of DNA sequences shows cross-divide assemblage characteristics that correlate with the geomorphic zonation. The Northern Apennines show higher values of the directional measures of assemblage change gain, loss, and turnover; the Ligurian Alps show higher values of overlap and nestedness. There is a positive correlation between divide asymmetry and genetic distance, and gain, loss, and turnover of DNA sequences from Adriatic to Ligurian sites; there is a negative correlation with overlap and nestedness. Since the species is confined to freshwater environments, tectonically driven drainage reorganization is one of the only mechanisms that can explain its spatial genetic differentiation.

To better understand the landscape dynamics and changes in habitat connectivity influenced by glacial and interglacial oscillations over the biodiversity-rich Hengduan Mountains (HM) region, high-resolution climate data for past periods are essential. We apply the non-hydrostatic limited-area model COSMO, with a resolution of 12 km over East Asia, to simulate the Last Glacial Maximum (LGM), a period characterized by a generally colder and drier climate compared to present-day conditions. We perform the downscaling with a novel approach for paleoclimate modelling, the Pseudo-Global Warming (PGW) method. The COSMO PGW simulation for the LGM shows that COSMO generally replicates the large-scale dynamics of the driving global climate model simulation in the colder climate. Both models suggest weaker Asian summer monsoon systems during this period. Consequently, regions such as the Bay of Bengal, and the South China Sea, which typically receive substantial monsoon rainfall, experience significantly reduced precipitation. However, despite these model similarities, the high-resolution COSMO simulation exhibits distinctive differences on a smaller scale—particularly over land. For instance, COSMO suggests a more pronounced southward shift of the jet stream during the LGM winter, with cooler conditions in southern China. Moreover, the COSMO simulation, despite the overall weaker summer monsoon circulation, features increased precipitation amounts for much of the HM. Additionally, COSMO suggests a more extensive increase in snowfall over the High Mountain Asia region. Our study suggests that the resource-saving PGW approach is a suitable method to bridge the gap between large-scale projections and regional climate impacts—also for past periods like the LGM.

The Hengduan Mountains (HM) are located on the southeastern edge of the Tibetan Plateau (TP) and feature high mountain ridges (> 6000 m a.s.l.) separated by deep valleys. The HM region also features an exceptionally high biodiversity, believed to have emerged from the topography interacting with the climate. To investigate the role of the HM topography on regional climate, we conduct simulations with the regional climate model COSMO at high horizontal resolutions (at ~12 km and a convection-permitting scale of ~4.4 km) for the present-day climate. We conduct one control simulation with modern topography and two idealised experiments with modified topography, inspired by past geological processes that shaped the mountain range. In the first experiment, we reduce the HM’s elevation by applying a spatially non-uniform scaling to the topography. The results show that, following the uplift of the HM, the local rainy season precipitation increases by ~25%. Precipitation in Indochina and the Bay of Bengal (BoB) also intensifies. Additionally, the cyclonic circulation in the BoB extends eastward, indicating an intensification of the East Asian summer monsoon. In the second experiment, we remove the deep valley by applying an envelope topography to quantify the effects of terrain undulation with high amplitude and frequency on climate. On the western flanks of the HM, precipitation slightly increases, while the remaining fraction of the mountain range experiences ~20% less precipitation. Simulations suggest an overall positive feedback between precipitation, erosion, and valley deepening for this region, which could have influenced the diversification of local organisms.

Understanding the intricate dynamics of biodiversity within and across riverine ecosystems, influenced by geological history and environmental factors, is crucial for effective conservation and management strategies. Italy, particularly the Ligurian region, harbors diverse freshwater fish communities and populations shaped by unique geological and hydrological conditions. Here, we investigated the suitability of eDNA-metabarcoding to identify inter- and intraspecific diversity patterns of riverine fish populations along the main drainage divide (MDD) between the Adriatic and Ligurian basins in Northern Italy. We collected 96 aquatic eDNA samples across 48 riverine sites, amplified them using a cytb primer pair and denoised the sequences to retrieve amplicon sequence variants (ASV). We calculated communities’ phylogenetic distance with betaMPD based on genetic distances derived from the ASVs, combined them with conductance-based landscape metrics and applied generalized dissimilarity models (GDM) to assess spatial genetic structure. Our results reveal genetic differentiation among populations of several fish species, with some displaying clustering patterns across the drainage divide and isolation by distance patterns. Overall, taxon richness was significantly higher in the Ligurian sites (25) than in the Adriatic side of the MDD (22), as was ASV richness (205 vs. 196). Our findings highlight the effectiveness of eDNA-metabarcoding in uncovering various facets of diversity, shedding light on hidden genetic diversity within ASVs, and revealing significant spatial genetic structuring in freshwater fish populations across multiple species.

Mixed siliciclastic and carbonate active orogens are common on Earth’s surface, yet most studies have focused on physical erosion and chemical weathering in silicate-rich landscapes. Relative to purely siliciclastic landscapes, the response of erosion and weathering to uplift may differ in mixed-lithology regions. However, our knowledge of weathering and erosion in mixed carbonate-silicate lithologies is limited and thus our understanding of the mechanistic coupling between uplift, chemical weathering, and the carbon cycle. Here, we partition the denudation fluxes into erosion and weathering fluxes of carbonates and silicates in the Northern Apennine Mountains of Italy—a mixed siliciclastic-carbonate active orogen—using dissolved solutes, the fraction of carbonate sand in sediments, and existing 10Be denudation rates. Erosion fluxes are generally an order of magnitude higher than weathering fluxes and dominate total denudation. The contribution of carbonate and silicate minerals to erosion varies between lithologic units, but weathering fluxes are systematically dominated by carbonates. Silicate weathering may be limited by reaction rates, whereas carbonate weathering may be limited by acidity of the rivers that drain the orogen. Precipitation of secondary calcite from super-saturated streams leads to the loss of up to 90% of dissolved Ca2+ from carbonate-rich catchments. Thus, in the weathering zone, [Ca2+] is exceptionally high, likely driven by high soil pCO2; however, re-equilibration with atmospheric pCO2 in rivers converts solutes back into solid grains that become part of the physical denudation flux. Limits on weathering in this landscape therefore differ between the subsurface weathering zone and what is exported by rivers.

The eastern margin of Madagascar has a prominent relief change from the flat coastal plain to the low-relief high plateau, characterizing a typical great escarpment topography at a passive margin. A quantification of the spatial distribution of erosion rates is necessary to understand the rate of landscape evolution. We present catchment-averaged erosion rates from detrital cosmogenic 10Be concentrations, systematically covering distinct morphological zones of the escarpment. Erosion rates are differentiated across the escarpment, where the high plateau and the coastal plain are slowly eroding with an average rate of 9.7 m/Ma, and the escarpment basins are eroding faster with an average rate of 16.6 m/Ma. The Alaotra-Ankay Graben related basins have the highest erosion rate with an average rate of 27 m/Ma. The spatial pattern of erosion rates indicates a retreating escarpment landscape. Retreat rates calculated from the 10Be concentrations are from 182 m/Ma to 1886 m/Ma. The rates of escarpment retreat on Madagascar are consistent with a model of a steady retreat from the coastline since the time of rifting, similar to the Western Ghats escarpment on its conjugate margin of the India Peninsula.

A great escarpment at a rift margin can correspond directly to a major water divide but at many margins, including in Madagascar, the escarpment often appears as a steep knickzone on rivers that have their main water divides in the interior of the high plateau. We hypothesize that this variability in morphology is a reflection of the frequency and size of drainage area captured from the high plateau over the escarpment. To test this hypothesis, we document morphological features and weathering conditions from river sediment of Madagascar. We propose that the existence of a weathered weak surface layer of crystalline bedrock encourages large river captures from the upper plateau, leading to a dominance of knickpoint type rivers. We demonstrate that this is feasible, using 2-D landscape evolution models and show that an easily-eroded surface layer is prone to fast divide migration through frequent river capture and reversal. A positive scaling relationship between captured area and escarpment retreat rate is found from models. We demonstrate that this scaling is also observed in the great escarpments of Madagascar and India.