# Human-triggered magnification of erosion rates in European Alps since the Bronze Age

William Rapuc, Charline Giguet-Covex, Julien Bouchez, Pierre Sabatier, Jérôme Gaillardet, Kévin Jacq, Kim Genuite, Jérôme Poulenard, Erwan Messager, Fabien Arnaud

► To cite this version:

William Rapuc, Charline Giguet-Covex, Julien Bouchez, Pierre Sabatier, Jérôme Gaillardet, et al. Human-triggered magnification of erosion rates in European Alps since the Bronze Age. Nature Communications, 2024, 15 (1), pp.1246. 10.1038/s41467-024-45123-3. hal-04465041

HAL Id: hal-04465041

https://hal.science/hal-04465041v1

Submitted on 19 Feb 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# Human-triggered magnification of erosion rates in European Alps since the Bronze Age

Received: 11 July 2023

Accepted: 12 January 2024

Published online: 10 February 2024

<figure><img src="image_3.png" /></figure>

Check for updates

William Rapuc Ⓔ<fnref n="1" />,²✉, Charline Giguet-Covex Ⓔ<fnref n="1" />, Julien Bouchez Ⓔ³, Pierre Sabatier Ⓔ<fnref n="1" />, Jérôme Gaillardet³,⁴, Kévin Jacq⁵, Kim Genuite⁶, Jérôme Poulenard Ⓔ<fnref n="1" />, Erwan Messager<fnref n="1" /> & Fabien Arnaud Ⓔ<fnref n="1" />

A major feature of the Anthropocene is the drastic increase in global soil erosion. Soil erosion is threatening Earth habitability not only as soils are an essential component of the Earth system but also because societies depend on soils. However, proper quantification of the impact of human activities on erosion over thousands of years is still lacking. This is particularly crucial in mountainous areas, where the highest erosion rates are recorded. Here we use the Lake Bourget catchment, one of the largest in the European Alps, to estimate quantitatively the impact of human activities on erosion. Based on a multi-proxy, source-to-sink approach relying on isotopic geochemistry, we discriminate the effects of climate fluctuations from those of human activities on erosion over the last 10,000 years. We demonstrate that until 3800 years ago, climate is the only driver of erosion. From that time on, climate alone cannot explain the measured rates of erosion. Thanks to an unprecedented regional paleoenvironmental reconstruction, we highlight that the development of pastoralism at high altitudes from the Bronze Age onwards and the extension of agriculture starting in the Middle Ages were key factors in the drastic increase in erosion observed in the Alps.

Although one of the major current environmental issues, soils degradation is still too often overlooked<fnref n="1" />. Yet, among other amenities, soils provide the substrate for 95% of human food and constitute the most efficient CO₂ sink at Earth's surface<fnref n="1" />⁻³. Soils are hence a vital resource for human societies, nowadays facing severe degradation due to the increasing extent of cultivated and urbanised lands, which leads to a loss in soil ecosystem functions⁴,⁵ and a dramatic rise in global erosion rates⁶⁻⁸. The deleterious potential of soil erosion and the corresponding lack of knowledge are such that Europe and the United Nations have recently called for an improved quantitative evaluation of soil loss over large spatial and long temporal scales²,⁹. Ideally, such assessment should in particular include a quantification of the effects of various factors driving erosion, in order to model their evolution and improve resource management policies²,⁶. Therefore, several

research groups have carried out modelling experiments at European and global scales to draw a picture of the current⁹,¹⁰ and future³ states of erosion rates. The prominent role of human activities in soil degradation and global erosion rise was thereby evidenced⁶,⁷,¹⁰. However, to date a critical knowledge gap exists as to when agriculture and grazing became the most significant driver of erosion, due to the compounding effects of climate shifts over the Holocene period; nor has it been possible to quantify accurately the effects of land use on erosion beyond decades or, at most, a few centuries.

Lake sediment sequences host a wealth of information to reconstruct the response of soils to climate and land use over several centuries to millennia¹¹⁻¹³. In lake sediment studies, sedimentation rates¹² or detrital particle fluxes¹⁴,¹⁵ from the corresponding catchments are commonly used as proxies in that aim. Although such data are

<footnote n="1">¹EDYTEM, CNRS, Université Savoie Mont Blanc, 73000 Chambéry, France. ²Department of Earth Sciences, Durham University, Durham DH1 3LE, UK.
³Université Paris Cité, Institut de Physique Du Globe de Paris, CNRS, 75005 Paris, France. ⁴Institut Universitaire de France, Paris, France. ⁵Laboratoire Commun SpecSolE, Envisol, CNRS, Université Savoie Mont Blanc, 73000 Chambéry, France. ⁶UMR PACEA 5199, CNRS, Université de Bordeaux, 33615 Pessac, France. ✉e-mail: william.l.rapuc@durham.ac.uk</footnote>
informative about changes in catchment erosion rates through time, they do not allow for assessing the relative roles of different forcing factors, such as climate and human activities¹⁶. Independent reconstructions of land cover and local climate obtained from models¹² or proxy data¹⁷ can help to identify the potential correlations between erosion rates and these forcing factors but do not allow for a quantification of these effects.

Here we apply an original source-to-sink approach¹⁶ to a 10,000-years-long lake sediment archive (Supplementary Fig. 4) representative of the western European Alps (Supplementary Fig. 1). The Alps presents currently the highest erosion rates over the entire European continent³,⁹, making it a hotspot of soil loss¹⁸. In addition, as one of the best documented places on Earth in terms of climate fluctuations¹⁹⁻²³ during the Holocene and of human activities from the Neolithic period²⁴⁻²⁷, this region represents an ideal environment to quantify the impact of both climate and human activities on erosion rates. Lake Bourget drains one of the largest alpine catchments (4976 km²) which is particularly well suited for such a study as it hosts a diversity of landscapes in terms of land use, geomorphology, and lithology²⁸ (Supplementary Fig. 1). The catchment culminates at 4808 m at the top of the Mont Blanc massif, where the presence of glaciers precludes any land use, such that erosion there is largely controlled by climate²⁹,³⁰. Moreover, most of this glacier-covered part of the catchment (158.1 km² situated between 2300 and 4808 m) is underlain by granites and metamorphic rocks³¹ with a geochemical signature distinct from those of the sedimentary formations underlying the lower elevation area³²⁻³⁴ (3803.9 km² excluding limestone formations; 232–3057 m), where agro-pastoral activities have developed over the last 7000 years to become a major human activity in the region²⁴,²⁷,³⁵. This configuration can be leveraged to quantify the role of various driving factors on erosion in the Lake Bourget catchments, given the following assumptions: (i) erosion in the glacier-covered part of the catchment responds only to climate (temperature and precipitation), whereas (ii) erosion in lowermost, non-glaciated regions depends both on climate and land use; and finally (iii) if climate shifts were the sole factor of changes in erosion, they would affect equally any part of the catchment area. Several studies support the first of these assumptions³⁶⁻³⁸ showing that in glaciated regions, climate and glacier thermal regime are the two main factors controlling erosion and the only factors explaining the variations on the time scales of our study³⁶. Furthermore, no human activity likely to have enhanced erosion has yet been identified on the Mont Blanc massif²¹. The second assumption is commonly accepted and form the basis of all models based on soil loss equation (e.g., RUSLE models)³,⁹,¹⁰. To demonstrate the validity of the third assumption in the context of our study area we have chosen a sedimentary sequence long enough to cover periods when anthropogenic effects on erosion were negligible leaving climate as the only controlling factor for the oldest part of the record. Based on these assumptions, we can quantify the evolution of climate- vs. human-driven erosion through time, such that an anthropogenic excess in erosion can be estimated. To ascertain that this excess erosion is indeed linked to agropastoralism, we employ an approach based on a synthesis of indices of agro-pastoral activities at the scale of the western European Alps. Agro-pastoral activities derived from pollen or environmental DNA data from 8 lakes covering an altitudinal gradient in or near the Lake Bourget catchment were synthesised through the calculation of a site-specific index of “agro-pastoral activity”. This index is blind to other human activities such as deforestation and mining. Because each altitudinal range presents its own environmental trajectory over the Holocene²⁴, these indexes were then binned by altitudinal ranges, that are named foothills (250–850 m), mid-altitude (850–2000 m) and high altitudes (>2000 m). Our source-to-sink approach involves the collection and measurement of 29 samples of fine river sediments from the flood plain of the Rhône and the Arve rivers, the main tributaries of Lake Bourget, and their tributaries (Methods). Eighty-three samples

were also collected in the sediment sequence retrieved in the depo-center of Lake Bourget (Methods). Erosion flux was estimated from the sediment sequence by the computation of the sediment dry bulk density, the sedimentation rate, and the proportion of detrital silicates. Isotopic ratios of neodymium, as well as major and trace element concentration, were measured on the entirety of these samples. These radiogenic isotope systems are commonly used to trace the sources of erosion products³⁹,⁴⁰ and due to the different geological nature of glaciated and non-glaciated regions, isotopic ratios of neodymium offer the possibility to track the respective contribution of each region to erosion through time, based on the use of a geochemical mixing model (Fig. 1, Methods).

## Results and discussion

### Recording the erosional effects of climate and land use in the Alps

Here, we obtain two erosion signals (Fig. 2) corresponding to the glaciated region and to the rest of the Lake Bourget catchment, respectively. Erosion signals are expressed as Sediment Yield (SY) value (t.km⁻²·yr⁻¹) that yields only relative information on erosion rates in sediment source regions as (i) sediment deposition may occur on floodplains during transport to the lake and (ii) the Rhône River, which supplies sediment to Lake Bourget, only does so while flooding, i.e. during a few days each year⁴¹. Therefore, no data comparison is possible with other SYᵢ values such as those obtained from soil erosion models⁹,⁴², such that only variations through time and relationship between erosion rates in the two parts of the catchment is discussed in the present study. Thanks to a precise age-depth model (Supplementary Fig. 4), with a median uncertainty of 60 years and a maximum uncertainty of no more than 100 years over the last 6700 years, it is possible to accurately describe the evolution of these erosion signals over time. Both erosion signals remain low and relatively constant from 9440 to 4880 yr cal BP (Fig. 2a, b). A first increase in erosion rate is observed for both regions from 4,880 to 4,065 yr cal BP. During this period, called the Neoglacial period⁴³, the Alps experienced a transition towards a wetter and colder period compared to the Early Holocene⁴⁴. This change in climate led to glacier advance at high altitudes⁴⁴⁻⁴⁶. The first phase of increase in the erosion signal being similarly observed in both regions is therefore linked to the effect of one of the main climatic transitions of the Holocene in the European Alps. From 3860 yr cal BP onwards, the two regions show different trends, with an increasing erosion rate in the non-glaciated part of the watershed, while the rate remains constant in the glaciated region with the exception of three periods listed hereafter, which corresponds to known phases of glacial advance in the area²¹,⁴³. Indeed, during the cold periods of the Migration period (1700 to 1550 yr cal BP), the Early Middle Ages (950 to 700 yr cal BP), and the end of the Little Ice Age (280–140 yr cal BP), the two regions show a similar increase in erosion or an excess of erosion from the glaciated area (Fig. 2a, b). Over the last thousand years, the intensity and variability of the erosion signal in the glaciated region has far exceeded anything previously recorded (Fig. 2b). This can be explained by the fact that never in recent millennia have the region's glaciers had such a volume and covered such a large area²¹,⁴⁷. This means that the region's climate is wetter and colder than during the rest of the Holocene²³,⁴⁸, favouring erosion as well as the massive advance of the glaciers, which will also reinforce the erosion of the area.

### First effects of human activities on erosion 3800 years ago

To quantify the effect of human activities on erosion, the effect of climate fluctuations must be isolated first. The correlation between the erosion signals from the two regions for the Early Holocene (Fig. 2a, b) seems to indicate that before 3860 yr cal BP climate alone was responsible for the erosion of the entire catchment. Indeed, a significant linear relationship (r = 0.86, p-value > 10⁻⁷) can be
<figure><img src="image_4.png" /></figure>

**Fig. 1 | Conceptual model, presenting a 3D-view of a catchment adapted to the approach of the present study.** Example here of Lake Bourget catchment (see Methods), hosting two regions, each of which is characterised by a specific geochemical signature, and in which erosion of one of the regions is influenced mostly by climate fluctuations. To obtain an erosion signal for each of the two regions, isotopic measurements on sediment samples from the source regions, as well as on the lake sedimentary sequence, are needed. These data are used to build a mixing

model to obtain the contribution of each region to the erosion signal recorded in the lake. Combined with the knowledge of the total erosion flux in the whole catchment, it is therefore possible to obtain an erosion signal per region over time. The signal from the glaciated region (in blue) corresponds to climate-driven erosion only, while the second one (brown to orange shading) is influenced by both climate and human activities (see Methods for detail).

identified between 9440 and 3860 yr cal BP (Methods). This observation leads us to hypothesize the existence of a linear relationship between the effect of climate on erosion in both glaciated and non-glaciated regions in the European Alps for the Early Holocene. This observation is valid both for a hot, dry period when local glacier extension is very limited (Early Holocene, 9440–4800 yr cal BP) and for a period when the area returns to a wetter, colder climate (Neoglacial, 4800–3800 yr cal BP)²³,²⁸,⁴⁴–⁴⁶,⁴⁸. This validates the third of our initial assumptions, according to which if climate change were the only factor in the evolution of erosion, it would affect any part of the catchment area in the same way. The differing trends between the two erosion signals from 3860 yr cal BP onwards indicate the appearance of a new driver of erosion in the non-glaciated part of the catchment, superimposed on climate. From that linear relationship an erosion rate can be computed for the non-glaciated region, that should reflect only the effect of climate on the erosion for the rest of the Holocene (Supplementary Fig. 7 and Supplementary Method 6). The ratio between this erosion rate expected from climate effects alone, and the measured erosion rate

corresponds to the excess in erosion in the non-glaciated region, which is not explained by climate fluctuations. Variations in the signal of erosion excess that do not exceed the variability observed over the correlation period (9.440–3.840 yr cal BP) will not be discussed, as such variations may be linked to the noise inherent in the archive used.

From 3800 yr cal BP onwards, climate alone cannot explain the erosion rates recorded in the Lake Bourget catchment. Because it is impossible that large-scale tectonic forcing could have affected only one of the two regions, we attribute the excess in erosion to human activities. The increase in the impact of agro-pastoral activities occurred steadily in the Bronze Age (Fig. 2c) before leading to a first occurrence of a high value of erosion excess at the end of the Bronze Age (3275 and 3025 yr cal BP). The other high excess erosion values are observed during: (i) the transition between Iron Age and Roman Period (2020 yr cal BP); (ii) at the end of the Early Middle Age (1350–1230 yr cal BP), at the end of the High Middle Age (750–710 yr cal BP), and at the end of the Early Modern (215 yr cal BP). A decrease is observed between 615 (1335 CE) and 490 yr cal BP (1460 CE), which, considering the
<figure><img src="image_5.png" /></figure>

**Fig. 2 | Erosion signals for non-glaciated and glaciated regions and the relative effect of human activities on erosion rates in Lake Bourget catchment.** Sediment yields (expressing erosion rates) for both the (a) non-glaciated (in orange) and (b) glaciated regions (in blue) of the Lake Bourget catchment, and (c) the effect factor of human activities on erosion in the non-glaciated region (blue to red shading), computed as a ratio between the measured erosion rate and the expected erosion rate in the absence of human impact (Supplementary Method 6; Supplementary Fig. 7). These two erosion signals were computed from the detrital silicate flux, the mixing model, the εNd data obtained on both the sedimentary sequence and the river sediments, and by calculating the areas of the glaciated and non-glaciated region in the catchment (Method). The error bars correspond to the standard deviation obtained using a Monte Carlo approach (Method). Colour shading of erosion signals in glaciated and non-glaciated regions

highlights increases in intensity compared with early Holocene values. On panel (a), the erosion signal for the glaciated region has been normalised to the correlation period values to highlight periods where erosion rates are dissimilar between the two regions. On panel (c), the magnification scale represents the effect of human activities on erosion compared with the theoretical value calculated for erosion impacted solely by climatic fluctuations. Periods: Pre-anthropic (until 5250 yr BP), Neolithic (5250–4200 yr BP), Early Bronze-Age (4200–3600 yr BP), Mid to Late Bronze-Age (3600–2600 yr BP), Iron Age (2600–1980 yr BP), Roman Period (1980–1750 yr BP), Post (P) Roman Period (1750–1474 yr BP), Early (E) Middle Age (1474–950 yr BP), High (H) Middle Age (950–650 yr BP), Late (L) Middle Age (650–458 yr BP), Early (E) Modern (458–161 yr BP), Late (L) Modern (161–4 yr BP), Co: Contemporary (4 yr BP – present). Source data are provided as a Source Data file.

uncertainties of the age model, could be explained by the decrease in population after the Black Death⁴⁹⁻⁵¹. Only the Migration Period (1700 to 1550 yr cal BP) and the end of the Little Ice Age (140 yr cal BP; 1810 CE) show excess values below 1. These two periods were characterised by a transition to a colder and wetter climate¹⁹,²¹ than the average climate of the last two millennia²³. This had two effects: (i) increasing the surface area covered by glaciers²¹,⁵², (ii) and enhancing erosion in the glaciated region through the combined effect of increased precipitation and the advance of the glaciers. These two periods are also known to have seen a marked decline in human activity, linked to invasions during the Migration Period¹⁹,²⁴ and to the reduction in agricultural practices³³,⁵⁴ probably induced, for both periods by the shortening of summer²⁴. These combined effects explain the strong decrease observed in the erosion excess signal for these two periods (Fig. 2c). Regardless of the climatic context, and except for the end of the Little Ice Age (1680 to 1810 CE), agro-pastoral activities have caused a two-to-two-and-a-half-fold increase in erosion rates in the western European Alps (Fig. 2c). Hence, within the climatic context of the western European Alps, agro-pastoral activities have dominated erosion for the last 3800 years, making it necessary to consider human activities as a principal driver of soil erosion over these timescales.

## From indices of agro-pastoral activities to effects on erosion

Any link between human activities and erosion excess remains to be better understood. For this purpose, we have divided the non-glacial region of the Bourget catchment, which is the region affected by agro-pastoral activities, into three parts: (i) the foothills; (ii) the mid-altitude level; and (iii) the high-altitude level (Fig. 1). We hypothesize that the two to three sites selected by altitudinal range illustrate the intensity of pastoral and/or agricultural activities in these areas²⁴,²⁷,⁵⁵. This assumption is considered as robust as the history of activities recorded at the sites in a given altitudinal range are very similar (Supplementary Fig. 8). Because the Lake Bourget integrates the sum of erosion products over all these altitudinal ranges (and is therefore not impacted by local physical and geological characteristics leading to variability in local erosion response), we can assess and quantify the extent of human impact on erosion at such large scale. Using a previously developed methodology²⁴, we aggregate several tracers of grazing, pastoral, and agricultural activities as recorded by lake sediments into an index of intensity/diversity of agro-pastoral activities. Due to their different natures, each marker of human activity are standardised beforehand. The standardization consists in the normalization to 0–100% of the
<figure><img src="image_6.png" /><figcaption>**Fig. 3** | Comparison between index of agro-pastoral activities (brown shading) and the magnification effect of agro-pastoral activities on erosion rates (blue to red shading). Indexes were computed for three altitudinal ranges of the western European Alps. To facilitate the comparison, the chronicle of the effect of agro-</figcaption></figure>

pastoral activities on erosion is presented in the same form as the indexes (Fig. 2). A Period used to compute the excess in erosion; (B) Period of dominant influence from high-altitude areas; (C) Period of dominant influence from foothills/mid-altitudes areas. Source data are provided as a Source Data file.

values of all proxies. Once obtained for each lake, these normalised values are integrated in the calculation of an average index by lake. Finally, these indexes are grouped together by altitudinal range and for each chrono-cultural period (Fig. 3). The integration of data from different lakes and different data-types allows us to offset specific taphonomic issues associated with each proxy and lake. Sedimentary DNA and coprophilous fungi signals are positively affected by climatic or human-triggered erosion, potentially leading to an overestimation of human activities, while pollen counting can be affected by a non-local origin of pollen grains and thus potentially reflect activities in other altitudinal ranges²⁴. However, for DNA and coprophilous fungi, because each lake catchment has its own sensitivity to erosion processes, potential biases are not expected to have had an identical pattern from one site to another. The binning of the data by chrono-cultural periods is an additional means of reducing noise in the signals generated by potential taphonomic problems.

The first erosion excess in Lake Bourget is recorded during the Mid-Bronze Age (3600–2600 yr cal BP), a period where agro-pastoral activities increase at high altitudes (Fig. 3). This period corresponds to the arrival of the first flocks of sheep at high altitude⁵,⁵⁶. The vegetation was sparse, facilitating opening of the vegetation and

establishment of mountain pastures as early as 3800 yr cal BP⁵¹. Therefore, these data suggest that high-altitude pastoralism has acted as the first factor that has favored soil erosion in the Alps. Then, and until the end of the Roman Empire (1750 yr cal BP), excess erosion remained high but stable, despite an increase of agro-pastoral index activities at high-altitude (Fig. 3). This might be explained by reaching a threshold of soil available to erosion. Until this period, only the high-altitude areas of Lake Bourget catchment seem to host human activities strong enough to impact soil erosion. However, from the Early Middle Ages (1750 yr cal BP) onwards, the excess in erosion increases simultaneously with the evidence of agro-pastoral activities in the mid-altitude and lowland areas. Although high-altitude areas remained characterised by significant pastoralism until the Modern period (460 yr cal BP; 1490 CE), agricultural activities⁵¹, with the advent of the plough from the Early Middle Ages (1750 yr cal BP) onwards⁵⁷, seem to have become the main driver of erosion in the western European Alps. It is interesting to notice that even if High (950 to 650 yr cal BP) and Late Middle Ages (650 to 460 yr cal BP) present the highest indices of agro-pastoral activities at mid and high altitude, there is no significant increase in the effect of these activities on erosion. This can correspond to a “system adjustment” phase during which human activities have deliberately (or by chance) evolved to lessen soil erosion, or when
virtually all erodible soil was removed, as previously proposed²⁴,⁵⁸. Except for foothills, the indices of human activity show a decreasing trend from the Early Modern (460 yr cal BP; 1490 CE); probably related to progressive agricultural abandonment of mountain environments⁵⁹ leading to a gradual return of forest cover¹⁷,⁶⁰, which may also have been favored by climate change⁶¹. This trend is directly reflected in the erosion excess, which decreases until the end of the Late Modern (4 yr cal BP; 1946 CE). Foothills, in contrast, present a differing trend of human activity index towards the Contemporary period (Fig. 3). This suggest that the activities conducted in this area do contribute less in terms of erosion, probably due to more gentle slopes characterizing the foothills reducing both erosion rates and the efficiency of sediment removal.

This study highlights the importance of the retrospective approach for understanding the effect of human activities¹¹. Only this type of approach, based on multi-proxy analyses and long timescales, can deliver detailed information on the compounded effects of climate and human activities on soil erosion. Here, we constrain the pace at which agro-pastoral activities led to an amplification of erosion rates in one of the largest Alpine catchment areas over the Holocene. We show that even non-intensive early agro-pastoral activities triggered a direct and major effect on soil erodibility for at least 3800 years. The multiplication of similar approaches could help quantify the impact of human activities on soils at a broader scale, and demonstrate the early impact of human activities on the environment, challenging the hypothesis of an Anthropocene starting with the Industrial revolution⁶².

# Methods

## Coring survey, logging, analyses, and dating

In June 2018 and 2019, a total of 56.5 m of sediment was retrieved in the deep basin of Lake Bourget (145 m below lake surface; Supplementary Fig. 2), using 90-mm diameter piston corers (Uwitec, Mondsee, Austria) operated from an UWITEC platform (EDYTEM/LSCE/C2FN). Two-meter-long sections from different holes were taken with a 1-m offset to ensure overlap and provide a continuous record. Sections were split, photographed at high-resolution (20 pixels.mm⁻¹), described and logged in detail. Identification of specific layers in the overlapping sections combined with correlations of XRF-core scanner signals allowed the creation of a 13.7-m-long sediment sequence called “LDB18&19” (IGSN n°TOAE0000000005 & n° TOAE0000000006, 45°44.717°N; 5°51.789°E). To characterise the high-resolution variations in major elements throughout the sediment sequence, X-ray fluorescence (XRF) geochemical analyses were performed on the EDYTEM laboratory’s Core Scanner (Avaatech XRF Technology). A continuous 5-mm step measurement was applied with a unique run at 10 kV and 0.12 mA for 15 s to detect lightweight elements, such as Al, Si, K, Ca, Ti, Mn, and Fe. The XRF core scanning results are expressed as peak intensity counting (cps). The description of the main sedimentological results is presented in Supplementary Note 1. Ti signal was used to estimate the proportion of detrital silicates in the sediment sequence (Supplementary Method 2). To produce a precise age-depth model for the sediment sequence, we combined identification of historical events such as earthquake and flood deposits with ¹⁴C. Nineteen samples of vegetal macro-organic remains were used to perform ¹⁴C measurements at the LMC14 laboratory (CNRS). Dates were calibrated at 2 sigma using the Intcal20 calibration curve⁶³, and the age-depth model was performed using the R code package “clam”⁶⁴ and “bacon”⁶⁵ in R software (Supplementary Method 1). Age-depth model results are shown in Supplementary Fig. 4, Supplementary Table 1.

## Lake and river sediment samples

Fifty-nine 2 cm-thick sediment samples were collected using graduated syringes for geochemical analyses. A constant sampling step of 10 cm was applied where no event deposit interrupted the continuous sedimentation (Source Data). Details of the Sediment

dry bulk density (DBD) was measured using sediment samples of an invariant 5 mL volume. To quantify the organic matter and carbonate contents throughout the sediment sequence, loss on ignition (LOI) analyses were conducted following the method developed by Heiri⁶⁶ (Supplementary Fig. 3). The age-depth model of the sediment section was performed with “rbacon” and “clam” packages on R software thanks to nineteen samples of terrestrial macro remain (Supplementary Table 1). To characterise the source of sediments accumulated in the Lake Bourget deep basin, 19 river sediment samples were collected between January and February 2018 from flood deposits in all sub-catchments directly contributing to the Arve and Rhône rivers, close to the confluence with the main river (Supplementary Fig. 1; Source Data). Flood deposits were selected as it makes it possible to obtain fine sediment samples covering a long period and representative of each sub-basin. Ten samples were also collected from flood deposits in the main river (Arve & Rhône) at several hundred meters downstream of each confluence to test for the mixing between sub-catchment and main river sediments (Supplementary Fig. 1; Source Data). Sub-catchments covered by this sampling set represent the entire catchment area of the Lake Bourget. Samples were collected manually, and only the fraction finer than 63 µm was used, to ensure that river and lake sediments were comparable.

## Geochemical analyses

Lake and river sediment samples were first dried at 60 °C for 72 h in a laboratory oven and crushed in an agate mortar. Then lake sediment samples were decarbonated by two successive 5 mL HCl 0.5 N leaching steps to remove authigenic and detrital carbonates. All samples then underwent a digestion using concentrated acid mixtures of HF, HNO₃, and HCl. Major and trace elements were measured using an Agilent 7900 quadrupole ICP-MS at the PARI analytical platform of IPGP. The reported uncertainties were calculated using the algebraic propagation of blank subtraction and sample count standard deviations (n = 3). The NIST®SRM®2709a reference material (San Joaquin soil) was processed as samples and analyzed repeatedly during the sequences to evaluate the accuracy of the measurements. The detection limit was between 0.2 and 0.5 ppt depending on the element, and the internal errors were 5% on average. Nd isotope measurements were all performed at the PARI analytical platform of IPGP on a MC-ICPMS; Neptune (Thermo-Fisher Scientific). After powder digestion, Nd was separated from the sample matrix by extraction chromatography following the same method as in ref. 16. More details on Nd isotope ratio measurements are given in Supplementary Method 3.

## Mixing model

The Nd isotopic composition of the Lake Bourget catchment allows for the identification of two main end members of detrital material (Supplementary Method 4; Supplementary Fig. 6):

(i) The Mont Blanc Massif (1.A sub-catchment, see Supplementary Fig. 1), where large glacier bodies are concentrated. This area is part of the external massif and presents the highest εNd values (−5.36). We refer to this part of the Lake Bourget catchment as the glaciated region.

(ii) The rest of the catchment is primarily composed of calcareous sedimentary rocks and Quaternary deposits (Supplementary Fig. 1). In this area εNd values are lower (Source Data). We refer to this part of the Lake Bourget catchment as the non-glaciated region.

We use mean εNd values and [Nd] concentrations of the 1.A and 3.A sub-catchments to track the contribution of the glaciated and non-glaciated region erosion in the total silicate sediment accumulated in
the Lake Bourget, respectively (Source Data). The fractional con- tribution of each rock source can be estimated using a mixing model:

$$[Nd]^M = [Nd]^A \times f_A + [Nd]^B \times f_B \quad (1)$$

$$f_A + f_B = 1 \quad (2)$$

$$\varepsilon Nd^M \times [Nd]^M = \varepsilon Nd^A \times [Nd]^A \times f_A + \varepsilon Nd^B \times [Nd]^B \times f_B \quad (3)$$

where the superscript M stands for lake sediment (“mixture”) values, and the superscripts A and B stand for the source values (glaciated and non-glaciated regions, respectively); $f_x$ represents the fractional contribution of each rock source. Solving this equation system led to estimates of the relative contributions of sources A and B:

$$f_A = \frac{[Nd]^B * (\varepsilon Nd^B - \varepsilon Nd^M)}{[Nd]^A * (\varepsilon Nd^M - \varepsilon Nd^A) + [Nd]^B * (\varepsilon Nd^B - \varepsilon Nd^M)} \quad (4)$$

$$f_B = \frac{[Nd]^A * (\varepsilon Nd^A - \varepsilon Nd^M)}{[Nd]^B * (\varepsilon Nd^M - \varepsilon Nd^B) + [Nd]^A * (\varepsilon Nd^A - \varepsilon Nd^M)} \quad (5)$$

A Monte Carlo method was used to compute uncertainties in the fractional contributions. Random sampling within the analytical uncertainties of each parameter was executed using a Box-Muller transform⁶⁷ for each of the 5000 iterations of the Monte Carlo pro- cedure on Excel software. From the simulated parameter distributions, D84, D16, median and standard deviation were obtained for each proportion $f_x$ (Source Data). The same procedure was repeated taking into account uncertainty in the definition of the end member com- position. Mixing model results are displayed in Supplementary Fig. 6.

## Computation of the sediment yields

To obtain a quantitative assessment of catchment-scale erosion rates for each region in Lake Bourget catchment, river loads, and corre- sponding sediment yields were estimated following the methodology of ref. 13. First, the total amount of sediment stored in the Lake Bourget per amount of time ($Stock_{sed}$, in g.yr⁻¹) is estimated:

$$Stock_{sed} = SCF \times Lake \ area \quad (6)$$

with SCF corresponding to the siliciclastic flux (g.cm⁻².yr⁻¹), which represents the accumulation of detrital material in the lake Bourget over the Holocene recorded. This flux is obtained from the multi- plication of the sedimentation rate (SR, cm.yr⁻¹), the dry bulk density (DBD, g.cm⁻³) and the proportion of detrital silicates (PDS, fraction of the total sediment accumulation):

$$SCF = SR \times DBD \times PDS \quad (7)$$

Details about DBD, PDS and SCF computation are available in the Supplementary Method 2. SCF results are displayed in Supplementary Fig. 5. As SR varies spatially over the surface of the Lake Bourget⁶⁸ we consider only the fraction of the total lake area where SR should be commensurate with that measured at the coring site. This area (Lake area, in km²) is estimated at 14.4 km² and corresponds to the area of the lake with depths deeper than 110 m b.l.s⁶⁸. The sediment stock is then multiplied by the contribution of each source ($f_x$, see above) obtained from the mixing model to calculate changes in the source-

specific sediment loads (SL, in t.yr⁻¹) over time:

$$SL = Stock_{sed} \times f_x \quad (8)$$

To obtain a region-normalised sediment yield value for each source ($SY_i$), SL is normalised to the drainage area of each rock type (Source area, in km²):

$$SY_i = \frac{SL}{Source \ area} \quad (9)$$

with $SY_i$ expressed in t.km⁻².yr⁻¹. Details on source area computation are given in Supplementary Method 5.

## Identification of the expected effect of climate on erosion

Identified periods of correlation between erosion signal of both gla- ciated and non-glaciated regions help to obtain a relationship between both signal for the Early Holocene (Supplementary Fig. 7). To that end, we used the “zoo” package and the “rollapply” function of the R soft- ware to examine whether or not the erosion signals from the glaciated and non-glaciated regions are correlated. We produced rolling corre- lations with sample widths ranging from 3 to 17. The sliding correlation curves between 7*7 and 15*15 roll correlations were then averaged (Cmean) and extracted using Matlab software. Cmean was then derived with a Savitsky Golay polynomial derivative with a sliding window of 5 points and a second order polynomial filter⁶⁹. From the derivative of the signal, peaks are estimated which gives inflection points corresponding to the breakpoints in the correlation relation- ship of the erosion signals of both regions. Those breakpoints repre- sent the major shifts in the relation of the correlation in erosion between glaciated and non-glaciated regions. More detail is provided about the results of the rolling correlation and the derived linear relationship in the Supplementary Method 6 and in Supplementary Table 2.

## Index of agro-pastoral activities

The proposed index of human activity intensity is the aggregation of all proxies of agricultural and pastoral activities from lake sediment DNA analyses (targeting plants and animals) and/or pollen and coprophi- lous fungi spore analyses. The list of taxa included for each lake is presented in the Supplementary Table 3. Because these selected proxies are produced by different methods and expressed in different units (i.e., in number of positive replicates for DNA analyses, in per- centage for pollen analyses and in number of spores per cm³ for coprophilous fungi), they are normalised to 0–100%. In detail, this standardization procedure consists of changing the highest value of each proxy to 100% and calculating the other values proportionally. Then, the synthetic index is obtained by calculating the mean of the normalised proxies. We integrated the data from several study sites for each altitudinal range (Supplementary Fig. 8). The robustness of this index is underpinned by several hypotheses. First, the integration of proxies from different methods allows for reducing potential biases caused by taphonomic or methodological issues (e.g., *Sporormiella* sp. spores is a more sensitive proxy to the record of low-intensity activ- ities, than animal DNA is⁷⁰). Second, the integration of plant taxa with different ecological preferences related to the intensity of activities allows us to better capture these changes in intensity (e.g. *Plantago* sp. for extensive and intensive activities and *Rumex* sp. for intensive activities and/or the presence of resting places). Third, the integration of several study sites by altitudinal range, which in addition show similar temporal evolution (Supplementary Fig. 8), allows for more reliable reconstructions. In order to synthesize these data and to better highlight the general trends for each altitudinal zone, they are repre- sented via box plots, using “ggplot2” and “grid” packages and “geom_boxplot” function on R software, per cultural period (Source
Data). Finally, this index considers that the more diversified the activities, the greater the intensity of the activities. One must keep in mind that, in regard to assessing human impact on soil stability, this index does not consider deforestation due to difficulties in discriminating human- and climatic-induced forest cover loss on sites around the tree line. Therefore, we contend that the intensity of human activity reflected in our index is underestimated for some periods, especially from the Early/High Middle Ages in foothills and mid-altitude regions. To improve data comparison, the excess erosion signal was also displayed as a box plot for each chrono-cultural period. If we consider that a sample is 2 cm thick, it therefore covers a period that depends on the sedimentation rate. Considering the error in the sedimentation rate, 2 out of the 83 samples span two chrono-cultural periods. The excess erosion values for these two samples have therefore been associated with the box plots for the periods they cover. (Source Data).

## Reporting summary

Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

## Data availability

The authors declare that the data supporting the findings of this study are available within the paper and its Supplementary Information files. Source data are provided with this paper.

## Code availability

The authors will make available, on request, any computer code or detailed algorithms used to generate the results presented in the article.

## References

1. Montgomery, D. *Dirt – The Erosion of Civilizations*. (University of California Press, 2012).

2. Shukla, P. R. et al. Technical Summary. in *Climate Change and Land: an IPCC special report on climate change, desertification, land degradation, sustainable land management, food security, and greenhouse gas fluxes in terrestrial ecosystems* 1, 37–74 (2019).

3. Borrelli, P. et al. Land use and climate change impacts on global soil erosion by water (2015–2070). *Proc. Natl Acad. Sci.* **117**, 21994–22001 (2020).

4. Wagg, C., Bender, S. F., Widmer, F. & van der Heijden, M. G. A. Soil biodiversity and soil community composition determine ecosystem multifunctionality. *Proc. Natl Acad. Sci.* **111**, 5266–5270 (2014).

5. de Graaff, M.-A., Hornslein, N., Throop, H. L., Kardol, P. & van Diepen, L. T. A. Chapter One - Effects of agricultural intensification on soil biodiversity and implications for ecosystem functioning: a meta-analysis. in *Advances in Agronomy* (ed. Sparks, D. L.) 155 1–44 (Academic Press, 2019).

6. FAO, I. Status of the world’s soil resources (SWSR)-main report. *Food Agric. Organ. U. N. Intergov. Tech. Panel Soils Rome Italy* 1, 1–607 (2015).

7. Montanarella, L. Agricultural policy: Govern our soils. *Nature* **528**, 32–33 (2015).

8. Olsson, L. et al. Land degradation: IPCC special report on climate change, desertification, land degradation, sustainable land management, food security, and 6 greenhouse gas fluxes in terrestrial ecosystems. in *IPCC Special Report on Climate Change, Desertification, Land Degradation, Sustainable Land Management, Food Security, and Greenhouse Gas Fluxes in Terrestrial Ecosystems* 345–436 (Intergovernmental Panel on Climate Change (IPCC), 2019).

9. Panagos, P. et al. The new assessment of soil loss by water erosion in Europe. *Environ. Sci. Policy* **54**, 438–447 (2015).

10. Borrelli, P. et al. An assessment of the global impact of 21st century land use change on soil erosion. *Nat. Commun.* **8**, 2013 (2017).

11. Dearing, J. A. & Jones, R. T. Coupling temporal and spatial dimensions of global sediment flux through lake and marine sediment records. *Glob. Planet. Change* **39**, 147–168 (2003).

12. Jenny, J.-P. et al. Human and climate global-scale imprint on sediment transfer during the Holocene. *Proc. Natl Acad. Sci.* **116**, 22972–22976 (2019).

13. Kemp, D. B., Sadler, P. M. & Vanacker, V. The human impact on North American erosion, sediment transfer, and storage in a geologic context. *Nat. Commun.* **11**, 6012 (2020).

14. Bajard, M. et al. Long-term changes in alpine pedogenetic processes: effect of millennial agro-pastoralism activities (French-Italian Alps). *Geoderma* **306**, 217–236 (2017).

15. Sabatier, P. et al. Long-term relationships among pesticide applications, mobility, and soil erosion in a vineyard watershed. *Proc. Natl Acad. Sci.* **111**, 15647–15652 (2014).

16. Rapuc, W. et al. Quantitative evaluation of human and climate forcing on erosion in the alpine Critical Zone over the last 2000 years. *Quat. Sci. Rev.* **268**, 107127 (2021).

17. Bajard, M. et al. Erosion record in Lake La Thuile sediments (Prealps, France): evidence of montane landscape dynamics throughout the Holocene. *Holocene* **26**, 350–364 (2016).

18. Golosov, V. & Walling, D. E. Erosion and sediment problems: Global Hotspots. **40** https://unesdoc.unesco.org/ark:/48223/pf0000369887/locale=en (2019).

19. Büntgen, U. et al. 2500 years of European climate variability and human susceptibility. *Science* **331**, 578–582 (2011).

20. Büntgen, U. et al. Recent European drought extremes beyond Common Era background variability. *Nat. Geosci.* **14**, 190–196 (2021).

21. Le Roy, M. et al. Calendar-dated glacier variations in the western European Alps during the Neoglacial: the Mer de Glace record, Mont Blanc massif. *Quat. Sci. Rev.* **108**, 1–22 (2015).

22. Luterbacher, J. et al. European summer temperatures since Roman times. *Environ. Res. Lett.* **11**, 024001 (2016).

23. Samartin, S. et al. Warm Mediterranean mid-Holocene summers inferred from fossil midge assemblages. *Nat. Geosci.* **10**, 207–212 (2017).

24. Giguet-Covex, C. et al. Long-term trajectories of mountain agroecosystems in the North-Western Alps. *Reg. Environ. Change* **23**, 58 (2023).

25. Garcés-Pastor, S. et al. High resolution ancient sedimentary DNA shows that alpine plant diversity is associated with human land use and climate change. *Nat. Commun.* **13**, 6559 (2022).

26. Pini, R. et al. Ecological changes and human interaction in Valcamonica, the rock art valley, since the last deglaciation. **29**, 19–34 (2016).

27. Messager, E. et al. Two millennia of complexity and variability in a perialpine socioecological system (Savoie, France): the contribution of palynology and sedaDNA analysis. *Front. Ecol. Evol.* **10**, 1–20 (2022).

28. Arnaud, F. et al. Erosion under climate and human pressures: an alpine lake sediment perspective. *Quat. Sci. Rev.* https://doi.org/10.1016/j.quascirev.2016.09.018 (2016).

29. Norton, K. P., Blanckenburg, Fvon & Kubik, P. W. Cosmogenic nuclide-derived rates of diffusive and episodic erosion in the glacially sculpted upper Rhone Valley, Swiss Alps. *Earth Surf. Process. Landf.* **35**, 651–662 (2010).

30. Wittmann, H., Blanckenburg, F. von, Kruesmann, T., Norton, K. P. & Kubik, P. W. Relation between rock uplift and denudation from cosmogenic nuclides in river sediment in the Central Alps of Switzerland. *J. Geophys. Res. Earth Surf.* **112**, 1–20 (2007).
31. Rolland, Y., Cox, S., Boullier, A.-M., Pennacchioni, G. & Mancktelow, N. Rare earth and trace element mobility in mid-crustal shear zones: insights from the Mont Blanc Massif (Western Alps). *Earth Planet. Sci. Lett.* **214**, 203–219 (2003).

32. Deconinck, J.-F., Beaudoin, B., Chamley, H., Joseph, P. & Raoult, J.-F. Contrôles tectonique, eustatique et climatique de la sédimentation argileuse du domaine subalpin français au Malm-Crétacé. *Contrôles. Tecton. Eustatique Clim. Sédimentation Argileuse Domaine Subalp. Fr. Au Malm-Crétac* **26**, 311–320 (1985).

33. Revel-Rolland, M. et al. Sr and Nd isotopes as tracers of clastic sources in Lake Le Bourget sediment (NW Alps, France) during the Little Ice Age: Palaeohydrology implications. *Chem. Geol.* **224**, 183–200 (2005).

34. Manalt, F., Beck, C., Disnar, J.-R., Deconinck, J.-F. & Recourt, P. Evolution of clay mineral assemblages and organic matter in the late glacial-Holocene sedimentary infill of Lake Annecy, (north-western Alps): paleoenvironmental implications. *J. Paleolimnol.* **25**, 179–192 (2001).

35. Bajard, M. et al. Pastoralism increased vulnerability of a subalpine catchment to flood hazard through changing soil properties. *Palaeogeogr. Palaeoclimatol. Palaeoecol.* **538**, 109462 (2020).

36. Koppes, M. et al. Observed latitudinal variations in erosion as a function of glacier dynamics. *Nature* **526**, 100–103 (2015).

37. Lehmann, B. et al. Postglacial erosion of bedrock surfaces and deglaciation timing: New insights from the Mont Blanc massif (western Alps). *Geology* **48**, 139–144 (2019).

38. Seguinot, J. & Delaney, I. Last-glacial-cycle glacier erosion potential in the Alps. *Earth Surf. Dyn.* **9**, 923–935 (2021).

39. Colin, C., Turpin, L., Bertaux, J., Desprairies, A. & Kissel, C. Erosional history of the Himalayan and Burman ranges during the last two glacial-interglacial cycles. *Earth Planet. Sci. Lett.* **171**, 647–660 (1999).

40. Cogez, A. et al. Constraints on the role of tectonic and climate on erosion revealed by two time series analysis of marine cores around New Zealand. *Earth Planet. Sci. Lett.* **410**, 174–185 (2015).

41. Arnaud, F. et al. Lake Bourget regional erosion patterns reconstruction reveals Holocene NW European Alps soil evolution and paleohydrology. *Quat. Sci. Rev.* **51**, 81–92 (2012).

42. Hinderer, M., Kastowski, M., Kamelger, A., Bartolini, C. & Schlu negger, F. River loads and modern denudation of the Alps — a review. *Earth-Sci. Rev.* **118**, 11–44 (2013).

43. Deline, P. & Orombelli, G. Glacier fluctuations in the western Alps during the Neoglacial, as indicated by the Miage morainic amphitheatre (Mont Blanc massif, Italy). *Boreas* **34**, 456–467 (2005).

44. Ivy-Ochs, S. et al. Latest Pleistocene and Holocene glacier varia tions in the European Alps. *Quat. Sci. Rev.* **28**, 2137–2149 (2009).

45. Simonneau, A. et al. Tracking Holocene glacial and high-altitude alpine environments fluctuations from minerogenic and organic markers in proglacial lake sediments (Lake Blanc Huez, Western French Alps). *Quat. Sci. Rev.* **89**, 27–43 (2014).

46. Luetscher, M., Hoffmann, D. L., Frisia, S. & Spötl, C. Holocene gla cier history from alpine speleothems, Milchbach cave, Switzerland. *Earth Planet. Sci. Lett.* **302**, 95–106 (2011).

47. Holzhauser, H., Magny, M. & Zumbuühl, H. J. Glacier and lake-level variations in west-central Europe over the last 3500 years. *Holo cene* **15**, 789–801 (2005).

48. Mauri, A., Davis, B. A. S., Collins, P. M. & Kaplan, J. O. The climate of Europe during the Holocene: a gridded pollen-based reconstruc tion and its multi-proxy evaluation. *Quat. Sci. Rev.* **112**, 109–127 (2015).

49. Büntgen, U. et al. Tree-ring indicators of German summer drought over the last millennium. *Quat. Sci. Rev.* **29**, 1005–1016 (2010).

50. Kaplan, J. O., Krumhardt, K. M. & Zimmermann, N. The prehistoric and preindustrial deforestation of Europe. *Quat. Sci. Rev.* **28**, 3016–3034 (2009).

51. Williams, M. Dark ages and dark areas: global deforestation in the deep past. *J. Hist. Geogr.* **26**, 28–46 (2000).

52. Nussbaumer, S. U. & Zumbühl, H. J. The Little Ice Age history of the Glacier des Bossons (Mont Blanc massif, France): a new high-resolution glacier length curve based on historical documents. *Clim. Change* **111**, 301–334 (2012).

53. Mather, A. S. The forest transition. *Area* **24**, 367–379 (1992).

54. Goldewijk, K. K. Estimating global land use change over the past 300 years: the HYDE database. *Glob. Biogeochem. Cycles* **15**, 417–433 (2001).

55. Doyen, E. et al. Land use development and environmental respon ses since the Neolithic around Lake Paladru in the French Pre-alps. *J. Archaeol. Sci. Rep.* **7**, 48–59 (2016).

56. Giguet-Covex, C. et al. Long livestock farming history and human landscape shaping revealed by lake sediment DNA. *Nat. Commun.* **5**, 1–7 (2014).

57. Sigaut, F. Les conditions d'apparition de la charrue. Contribution à l'étude des techniques de travail du sol dans les anciens systèmes de culture. *J. Agric. Tradit. Bot. Appliquée* **19**, 442–478 (1972).

58. Bajard, M. et al. Progressive and regressive soil evolution phases in the Anthropocene. *Catena* **150**, 39–52 (2017).

59. Carrer, F. The "Invisible" Shepherd and the "Visible" Dairyman: Ethnoarchaeology of Alpine Pastoral Sites in the Val di Fiemme (Eastern Italian Alps). *Summer Farms Seas. Exploit. Upl. Prehistory Present* **16**, 97–107 (2016).

60. Tasser, E., Leitinger, G. & Tappeiner, U. Climate change versus land- use change—what affects the mountain landscapes more? *Land Use Policy* **60**, 60–72 (2017).

61. Rumpf, S. B. et al. From white to green: snow cover loss and increased vegetation productivity in the European Alps. *Science* **376**, 1119–1122 (2022).

62. Stephens, L. et al. Archaeological assessment reveals Earth's early transformation through land use. *Science* **365**, 897–902 (2019).

63. Reimer, P. J. et al. The IntCal20 Northern Hemisphere radiocarbon age calibration curve (0–55 cal kBP). *Radiocarbon* **62**, 725–757 (2020).

64. Blaauw, M. Methods and code for 'classical'age-modelling of radiocarbon sequences. *Quat. Geochronol.* **5**, 512–518 (2010).

65. Blaauw, M. et al. rbacon: age-depth modelling using Bayesian Sta tistics. (2021).

66. Heiri, O., Lotter, A. F. & Lemcke, G. Loss on ignition as a method for estimating organic and carbonate content in sediments: reproducibility and comparability of results. *J. Paleolimnol.* **25**, 101–110 (2001).

67. Box, G. E. P. & Muller, M. E. A note on the generation of random normal deviates. *Ann. Math. Stat.* **29**, 610–611 (1958).

68. Jenny, J.-P. et al. A 4D sedimentological approach to reconstructing the flood frequency and intensity of the Rhône River (Lake Bourget, NW European Alps). *J. Paleolimnol.* **51**, 469–483 (2014).

69. Savitzky, A. & Golay, M. J. E. Smoothing and differentiation of data by simplified least squares procedures. *Anal. Chem.* **36**, 1627–1639 (1964).

70. Giguet-Covex, C. et al. New insights on lake sediment DNA from the catchment: importance of taphonomic and analytical issues on the record quality. *Sci. Rep.* **9**, 14676 (2019).

Acknowledgements

This work was financed by CRITLAKE project co-funded by EC2CO and AAP Université Savoie Mont Blanc. Parts of this work were also supported by IPGP multidisciplinary programme PARI and by Paris-IdF region SESAME Grant No. 12015903. We wish to thank the support of the entire staff of the C2FN-DT-INSU associated to the CLIMCORE project for the coring survey. ¹⁴C analyses were acquired thanks to the CNRS-INSU ARTEMIS national radiocarbon AMS measurement programme at Laboratoire de Mesure ¹⁴C (LMC14) in the CEA Institute at Saclay (French
Atomic Energy Commission). We are grateful to Pierre Burckel for help with Q-ICP-MS measurements; to Laëticia Faure for her valuable help and assistance with lab work; to Pascale Louvat, Thibaud Sondag, and Barthélémy Julien for assistance with MC-ICP-MS measurements; and to Raphaël Gallet, Bernard Fanget, Mathilde Banjan, Maude Biguenet, Matthieu Baril, Kim Genuite, Claire Blanchet, Kévin Jacq, Camille Girault, Manon Bajard, Ondřej Bábek, Qi Lin, Xiaqing Wang, Emmanuel Malet for the help during the different coring surveys and Jacques Mourey for help in the field for river sample collection. We wish to thank Damien Gui-noiseau for helpful conversations about the interpretation of isotope data. Finally, we would like to thank Carole Bégeot and Elise Doyen for providing us with the pollen data for Lac de Paladru.

## Author contributions

W.R., J.G., P.S. and F.A. designed the study. W.R. wrote the first draft of the paper. W.R. conducted the sedimentological and geochemical analyses of the Lake Bourget sediment sequence. W.R. and J.B. realised the mixing model and conducted statistical analyses on geochemical data. W.R. and C.G.-C. computed the agro-pastoral activities index. W.R. and K.J. developed the method of the human activities effect calcula-tion. W.R. and K.G. produced GIS data on the Lake Bourget catchment. W.R., C.G.-C., J.B., P.S., J.G., K.J., K.G., J.P., E.M. and F.A. contributed to interpreting the results. W.R., C.G.-C., J.B., P.S., J.G., K.J., K.G., J.P., E.M. and F.A. contributed to framing and revising the paper.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary information** The online version contains supplementary material available at https://doi.org/10.1038/s41467-024-45123-3.

**Correspondence** and requests for materials should be addressed to William Rapuc.

**Peer review information** *Nature Communications* thanks Daniel Schil-lereff, Francois Lapointe and the other, anonymous, reviewer for their contribution to the peer review of this work. A peer review file is available.

**Reprints and permissions information** is available at http://www.nature.com/reprints

**Publisher's note** Springer Nature remains neutral with regard to jur-isdictional claims in published maps and institutional affiliations.

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2024