Georges Aumaître

Didier L Bourlès

Karim Keddadouche

Franck Thomas
email: thomas@cerege.fr

Vincent Godard

Olivier Bellier

Lucilla Benedetti

Keywords: cosmogenic nuclides, denudation rates, carbonate, precipitation, climate, landscape evolution

The dynamics of carbonate-dominated landscapes is still debated, and first-order questions pertaining to the respective importance of chemical and mechanical weathering processes or the influence of climatic variability on denudation rates are largely open. To address these issues, we investigate the denudation pattern across a climatic gradient in Provence, SE France. We focused on a series of major carbonate mountain ranges with elevation spanning 200 to 2000 m and associated with more than 2-fold variations in mean annual temperature and precipitation. We collected 25 samples from 10 different mountain ranges and quantified denudation rates through the measurement of their 36 Cl concentrations in order to explore tectonic and climatic controls on these carbonate landscapes evolution. We observe that most of our samples experience denudation rates in the 30-60 mm/ka range (averaged over 5 to 40 ka), with no evident dependence on variations of mean annual temperature or precipitation. One range (Lure mountain) stands out as a significant outlier, with much higher denudation rates (>100 mm/ka) than adjacent sites (Ventoux mountain) at the same elevation and experiencing similar presentday environmental conditions. We explore in details the climatic, tectonic and geomorphological contexts of these two ranges in order to understand the reasons for this apparent 3-fold difference in denudation rates. A regionally-documented major Holocene hillslope destabilization event appears to be a likely candidate to explain the observed lower 36 Cl concentrations at the top of the Lure mountain. At longer wavelengths, our regional dataset shows no clear control of modern precipitation gradients on denudation rates, consistent with previous observations in landscapes underlain by other lithologies.

-Introduction

The respective influence of tectonic, climatic or anthropogenic forcings on landscapes evolution and their impact on the distribution of denudation rates in space and time is a major matter of debate in Earth Sciences (e.g. [START_REF] Portenga | Understanding earth's eroding surface with 10 Be[END_REF][START_REF] Riebe | Strong tectonic and weak climatic control of long-term chemical weathering rates[END_REF][START_REF] Harel | Global analysis of the stream power law parameters based on worldwide 10Be denudation rates[END_REF][START_REF] Olen | Role of climate and vegetation density in modulating denudation rates in the Himalaya[END_REF][START_REF] Scherler | Testing monsoonal controls on bedrock river incision in the Himalaya and Eastern Tibet with a stochastic-threshold stream power model[END_REF]. It is particularly interesting to understand the influence of climatic parameters such as temperature and precipitation on continental surfaces evolution and sediment fluxes (e.g. [START_REF] Tucker | Drainage basin responses to climate change[END_REF][START_REF] Roe | Orographic precipitation and the relief of mountain ranges[END_REF]. Indeed, over Quaternary time scales, landscape dynamics is believed to be strongly impacted by the drastic environmental changes associated with glacial cycles (e.g. [START_REF] Schaller | A 30 000 yr record of erosion rates from cosmogenic 10 Be in Middle European river terraces[END_REF][START_REF] Marshall | Frost for the trees : Did climate increase erosion in unglaciated landscapes during the late Pleistocene ?[END_REF][START_REF] Tofelde | 100 kyr fluvial cut-and-fill terrace cycles since the Middle Pleistocene in the southern Central Andes, NW Argentina[END_REF]. Understanding the impact of past climatic variability can help to better predict how landscapes will respond to projected changes in temperatures and precipitations in the future (Pelletier et al., 2015), across a wide range of spatial scales, from mineral surfaces to rock outcrops, landforms (hillslopes, hilltops or channels), regions and continents.

The Earth surface consists of a wide variety of lithologies including igneous, metamorphic and sedimentary rocks [START_REF] Hartmann | The new global lithological map database GLiM: A representation of rock properties at the Earth surface[END_REF]. Despite this diversity, the mineralogical composition is dominated by silicates, which chemical weathering processes are characterized by a progressive degradation of the primary minerals, leading to various secondary products that constitute a solid residual available for further transformation and transport by erosion agents (e.g. [START_REF] Hilley | Competition between erosion and reaction kinetics in controlling silicate-weathering rates[END_REF][START_REF] Mudd | Influence of Chemical Weathering on Hillslope Forms[END_REF]. From a surface processes perspective, carbonate lithologies stand out, due to the importance of congruent dissolution as a weathering agent, such that total dissolution can be achieved with only solute products left for transport (e.g. Ryb et al., 2014a). The existence of emblematic landscapes, such as karsts, associated with these specific lithologies demonstrate the role of dissolution as a potential driver of long-wavelength landscape evolution [START_REF] Kaufmann | Modelling karst denudation on a synthetic landscape[END_REF][START_REF] Fleurant | Modelling cockpit karst landforms[END_REF]. Contrastingly, a large part of our current understanding of the Earth surface dynamics relies on the analysis of processes responsible for solid mass transport by hillslopes and rivers (e.g. [START_REF] Roering | Hillslope evolution by nonlinear, slopedependent transport: Steady state morphology and equilibrium adjustment timescales[END_REF][START_REF] Whipple | Implications of sediment-flux-dependent river incision models for landscape evolution[END_REF], and how they contribute to the emergence of first order landscape attributes such as ridges and valleys topography [START_REF] Perron | Formation of evenly spaced ridges and valleys[END_REF]. In the case of carbonate landscapes, the potential dominance of a geomorphic agent whose activity is largely independent from topographic gradients question the applicability of usual Geomorphic Transport Laws to this type of environments [START_REF] Dietrich | Geomorphic transport laws for predicting landscape form and dynamics[END_REF]. Furthermore, the high dependence of calcite dissolution on temperature, water availability and pH also highlights the importance of a better assessment of the climatic controls on carbonate landscape evolution [START_REF] Kirstein | Experimental determination of natural carbonate rock dissolution rates with a focus on temperature dependency[END_REF][START_REF] Levenson | Comparison of field and laboratory weathering rates in carbonate rocks from an Eastern Mediterranean drainage basin[END_REF].

The variability and importance of these weathering processes for carbonate landscapes have lead to the development of numerous methods for the quantification of short-term carbonate dissolution (Gabrovsek, 2009), such as solute fluxes measurements (e.g. [START_REF] Gaillardet | Global silicate weathering and CO2 consumption rates deduced from the chemistry of large rivers[END_REF]Ryb et al., 2014), microtopography analysis (e.g. [START_REF] Furlani | Comparison between coastal and inland Karst limestone lowering rates in the northeastern Adriatic Region (Italy and Croatia)[END_REF] as well as in situ or laboratory calcite dissolution experiments [START_REF] Trudgill | Ten year remeasurement of chemical denudation on a magnesian limestone hillslope: Earth Surface Processes and Landforms[END_REF][START_REF] Thorn | Weathering rates of buried machine-polished rock disks[END_REF][START_REF] Liu | Comparative study of dissolution rate-determining mechanisms of limestone and dolomite[END_REF][START_REF] Plan | Factors controlling carbonate dissolution rates quantified in a field test in the Austrian alps[END_REF][START_REF] Kirstein | Experimental determination of natural carbonate rock dissolution rates with a focus on temperature dependency[END_REF]. However, these proxies focus only on direct dissolution of carbonate rock mass in controlled environments not necessarily representative of more general weathering conditions. Notably, they do not allow to assess the effect of important processes contributing to long-wavelength landscape denudation such as mechanical weathering and solid products transport. In addition, these methods provide annual or decadal information and do not integrate over Quaternary climatic variability at 1 to 10 ka time scales (e.g. [START_REF] Kirchner | Mountain erosion over 10 yr, 10 k.y., and 10 m.y. time scales[END_REF], which might make them sensitive to high-frequency forcings.

The introduction of cosmogenic nuclides in quantitative geomorphology studies has revolutionized our way to assess rates of landscape evolution, in particular by averaging processes over centuries to 10s of ka (e.g. [START_REF] Lal | Cosmic ray labeling of erosion surfaces: in situ nuclide production rates and erosion models[END_REF][START_REF] Gosse | Terrestrial in situ cosmogenic nuclides: theory and application[END_REF]Granger and Riebe, 2013). For this reason, 10 Be has become an essential tool to quantify the denudation of quartz-rich landscapes from the scale of the outcrop to that of large continental river basins (e.g. von [START_REF] Blanckenburg | The control mechanisms of erosion and weathering at basin scale from cosmogenic nuclides in river sediment[END_REF], and has led to fundamental advances in our understanding of the evolution of these landscapes such as testing the relative importance of tectonic and climatic forcings on denudation (e.g. [START_REF] Riebe | Strong tectonic and weak climatic control of long-term chemical weathering rates[END_REF][START_REF] Portenga | Understanding earth's eroding surface with 10 Be[END_REF], or probing proposed Geomorphic Transport Laws (e.g. [START_REF] Ouimet | Beyond threshold hillslopes: Channel adjustment to base-level fall in tectonically active mountain ranges[END_REF][START_REF] Scherler | Testing monsoonal controls on bedrock river incision in the Himalaya and Eastern Tibet with a stochastic-threshold stream power model[END_REF]. In comparison, our ability to quantify denudation or study the evolution of carbonate landscapes has significantly lagged behind. This is a critical blind spot because these landscapes occupy approximately 13% of continental surfaces on Earth [START_REF] Williams | Global distribution of carbonate rocks[END_REF] and are major reservoirs for the carbon cycle. 36 Cl concentration measurements in carbonate rock surfaces can be used to quantify their exposure or denudation history (e.g. [START_REF] Stone | Limestone erosion measurements with cosmogenic chlorine-36 in calcite -preliminary results from Australia[END_REF][START_REF] Schlagenhauf | Using in situ Chlorine-36 cosmonuclide to recover past earthquake histories on limestone normal fault scarps: A reappraisal of methodology and interpretations[END_REF]Ryb et al., 2014a). While it offers important opportunities to constrain the dynamics of carbonate landscapes, it has been far less used than 10 Be in quantitative, process-oriented geomorphology studies (Ryb et al., 2014a&b;[START_REF] Godard | Weather-limited hillslope evolution in carbonate landscapes[END_REF]. There are several first order pending problems that could be addressed with this approach, and in particular the question of the climatic controls on the denudation of carbonate surfaces, because the postulated major role of congruent dissolution might result in high landscape sensitivity to water availability. In this study, we document 36 Cl concentrations measured in a large dataset of carbonate rock surfaces samples collected from 11 mountain ranges across a climatic gradient in Provence, South-Eastern France, in order to assess the sensitivity of denudation to environmental parameters.

-Development and applications of cosmogenic nuclides in carbonate landscapes

The Earth surface is constantly being bombarded by cosmic ray particles, which induce nuclear reactions with atoms in the atmosphere and in rocks on the Earth surface (e.g. [START_REF] Gosse | Terrestrial in situ cosmogenic nuclides: theory and application[END_REF]. These interactions, such as spallation by neutrons, result in the production of various rare nuclides, called cosmogenic nuclides, within rock formations located a few meters below the Earth surface. The measurements of these in situ-produced cosmogenic nuclide concentrations in various substrates is of great interest for surface processes quantification. In particular, it has been demonstrated [START_REF] Lal | Cosmic ray labeling of erosion surfaces: in situ nuclide production rates and erosion models[END_REF] that for bedrock surfaces exposed to the cosmic ray particles and denuding at a constant rate for long enough, the induced cosmogenic nuclide production equilibrates the losses due to radioactive decay and mass removal linked to denudation processes. The cosmogenic nuclide concentration then reached at this steady-state is inversely proportional to the denudation rate of the surface (e.g. Granger and Riebe, 2013). This method provides a quantification of denudation rates, which are temporally averaged over the time period necessary to erode to a depth equivalent to the neutron characteristic attenuation length (approximately 60 cm in a substrate whose density is 2.5 g. cm -3 ) (Von [START_REF] Blanckenburg | The control mechanisms of erosion and weathering at basin scale from cosmogenic nuclides in river sediment[END_REF].

In quartz-rich environments, the 10 Be and 26 Al cosmogenic nuclides, mainly produced through the spallation of the constitutive O and Si atoms, are used to quantify either local denudation rates at specific sampling sites from bedrock samples or basin-averaged denudation rates from river sediments (Brown et al., 1995;[START_REF] Bierman | Estimating rates of denudation using cosmogenic isotope abundances in sediment[END_REF][START_REF] Granger | Spatially averaged long-term erosion rates measured from in situ-produced cosmogenic nuclides in alluvial sediment[END_REF][START_REF] Blanckenburg | The control mechanisms of erosion and weathering at basin scale from cosmogenic nuclides in river sediment[END_REF]. These approaches are now routinely used in quantitative geomorphology studies to determine the rates of surface processes over the 1-100 ka time window. They have provided valuable insights into first order issues as, for example, the validation and calibration of fluvial or hillslope erosion formalisms (e.g. [START_REF] Ouimet | Beyond threshold hillslopes: Channel adjustment to base-level fall in tectonically active mountain ranges[END_REF][START_REF] Dibiase | The influence of erosion thresholds and runoff variability on the relationships among topography, climate, and erosion rate[END_REF], the evaluation of the respective importance of tectonics and climate in controlling the denudation of mountain ranges (e.g. [START_REF] Bookhagen | Spatiotemporal trends in erosion rates across a pronounced rainfall gradient: Examples from the southern Central Andes[END_REF][START_REF] Godard | Dominance of tectonics over climate in Himalayan denudation[END_REF] or the assessment of the evolution of surface processes over the last Ma (e.g. [START_REF] Charreau | Paleo-erosion rates in Central Asia since 9Ma: A transient increase at the onset of Quaternary glaciations?[END_REF].

While cosmogenic nuclides have been intensively used to quantify denudation in settings dominated by quartz-rich lithologies for more than two decades, their use in areas dominated by carbonate lithologies is still in its early development stages. In such contexts, it is worth stressing again that congruent carbonate dissolution can act as a major weathering process, questioning the applicability of the classical slope-dependent Geomorphic Transport Laws [START_REF] Dietrich | Geomorphic transport laws for predicting landscape form and dynamics[END_REF] that have been validated by 10 Be derived denudation rate determinations in quartz rich settings (e.g. [START_REF] Ouimet | Beyond threshold hillslopes: Channel adjustment to base-level fall in tectonically active mountain ranges[END_REF].

In carbonate rocks, the 36 Cl cosmogenic nuclide is mainly produced through the spallation of Ca atoms and offers the possibility to investigate the exposure and denudation histories of carbonate surfaces. [START_REF] Dockhorn | Determination of erosion rates with cosmic ray produced 36Cl[END_REF] conducted one of the earliest studies proposing the measurement of 36 Cl concentrations in limestones or dolomites in order to estimate denudation rates. Later, [START_REF] Stone | Limestone erosion measurements with cosmogenic chlorine-36 in calcite -preliminary results from Australia[END_REF] used a similar approach to constrain the denudation rates of carbonate surfaces in Australia. They reported denudation rates of 5 mm/ka in various arid parts of Australia and significantly higher denudation rates of nearly 200 mm/ka in the wet tropical highlands of Papua New Guinea.

While these pioneering studies demonstrated the relevance of the use of 36 Cl to investigate surface processes in carbonated settings, the systematic use of this cosmogenic nuclide was still hampered at that time by significant uncertainties associated to some of its production pathways, including up to 50% difference in reported 36 Cl production rates from spallation of Ca [START_REF] Schimmelpfennig | Sources of in-situ 36 Cl in basaltic rocks. Implications for calibration of production rates[END_REF]. A significant methodological effort was undertaken to alleviate these issues and to calibrate the rates associated with the different 36 Cl production modes [START_REF] Schimmelpfennig | Sources of in-situ 36 Cl in basaltic rocks. Implications for calibration of production rates[END_REF][START_REF] Marrero | Quaternary Geochronology Cosmogenic nuclide systematics and the CRONUScalc program[END_REF].

These studies paved the way for further developments of the use of 36 Cl in quantitative geomorphology studies. [START_REF] Matsushi | Denudation rates of carbonate pinnacles in Japanese karst areas: estimates from cosmogenic 36Cl in calcite[END_REF] showed the importance of chemical dissolution in karst areas under subtropical climates in Japan and the greater impact of physical erosion in subarctic environments. [START_REF] Xu | In-situ cosmogenic 36Cl denudation rates of carbonates in Guizhou karst area[END_REF] came to similar conclusions studying samples collected in various climatic contexts in China and whose sampling altitudes ranged from 500 to 2300 meters above sea level. They determined denudation rates varying from 17 to 47 mm/ka that are roughly inversely proportional to the Mean Annual Temperature (MAT) at the sampling sites.

The authors suggest that, in their study area, chemical weathering is the dominant contributor to denudation, and that mechanical weathering could play a more important role in areas experiencing colder MAT, where freeze-thaw cycles occur. In their study of the Chauvet cave in the limestone karsts of Ardèche (France), Sadier et al. (2012) used 36 Cl to determine when the entrance of the cave was definitively closed by dating several rockfall events responsible for the accumulation of the sealing debris. In addition, the authors analyzed two samples from the surface located above the entrance and determined denudation rates similar to those from other studies under similar climates (18.5 and 21.5 mm/ka). Another surface 36 Cl measurement performed in South-Eastern France by [START_REF] Zerathe | Dating chert (diagenetic silica) using in-situ produced 10 Be: Possible complications revealed through a comparison with 36 Cl applied to coexisting limestone[END_REF] 36 Cl-based denudation rates increase linearly with MAP in sub-humid-Mediterranean climates but were independent from MAP in more humid climates. These results were interpreted in the sense that in semi-arid climatic conditions, chemical weathering is limited by the amount of precipitations whereas in more humid conditions, chemical weathering is apparently limited by the rate of carbonate mineral dissolution. In another study from the same area, Ryb et al. (2014b) compare long-term denudation rates to contemporary chemical weathering rates over variable time scales. Overall, their results suggest that on annual, decadal or 10 4 years time scales, denudation is mostly controlled by chemical weathering whereas mechanical weathering has a much more limited influence. Furthermore, this study highlights the importance of the type of sampled landform and the distribution of soil cover. [START_REF] Ryb | Exhumation and uplift coupled with precipitation along the western Dead Sea Rift margin[END_REF] provided paleo-precipitation gradient estimates along the western Dead Sea rift (10 4 -10 7 time scale) in order to shed light on the interactions between climate, denudation and tectonic uplift over the long term. Their study build on the idea that mass lost by dissolution over time is proportional to the time-integrated amount of precipitation. By combining total denudation estimated from a Turonian datum, denudation rates determined from 36 Cl concentrations and present-day mean annual precipitations (MAP), they propose that the effective precipitation gradient across the range was higher than the current one by at least 40% on average and that total precipitation at the 10 7 time scale correlates with both denudation rates measured at the 10 4 time scale and present-day MAP. [START_REF] Godard | Weather-limited hillslope evolution in carbonate landscapes[END_REF] studied the relationship between denudation rates and local curvatures along the bedrock crests of the Petit Luberon in South-Eastern France. Using 36 Cl derived denudation rates and a high-resolution Digital Elevation Model, they demonstrated a direct correlation between curvature and denudation at sites at which the curvature is lower than 0.02 m -1 and a roughly constant denudation rate (60-70 mm/ka) at sites characterized by higher curvatures.

These results imply that mechanical weathering and clast transport contribute significantly to hillslope evolution in this setting, and that the observed limit in the rate of denudation at the hilltop might contribute to relief growth in mediterranean settings. Focusing on the same mountain range, [START_REF] Thomas | Morphological controls on the dynamics of carbonate landscapes under a mediterranean climate[END_REF] investigated the importance of slope-dependent denudation processes in carbonate lithologies. The authors compared denudation rates from samples originating from the mountain crest and basin-averaged rates from the flanks where they observed a much faster evolution (70-200 mm/ka) than on the crest (30-60 mm/ka), implying that despite the importance of chemical weathering in carbonate lithologies, mechanical transport plays a major role in the context of this relatively low elevation range (600-700 meters a.s.l) receiving moderate precipitation.

Overall, these recent studies illustrate a complex response of carbonate landscapes to environmental conditions, both in terms of the controlling parameters and the limiting processes at work. Further research is clearly needed to pinpoint their respective influences and assess the full spectrum of their variations in the range of contexts where carbonate landscapes are present.

First order questions, in line with some of the most salient problems in Earth Surface sciences, are still open, such as the importance of climatic gradients on surface denudation. Indeed, while the existing studies have already addressed this problem, the aggregated existing data remains equivocal and could be interpreted either as showing an overall increasing trend in denudation rate with MAP or as being a transition between two regimes with low sensitivity of denudation beyond a threshold in MAP. In this study we take advantage of an important climatic gradient over a carbonate landscape in SouthEast France to further test this idea.

-Geological and geomorphological context

Provence landscapes in South-Eastern France are dominated by East-West trending ranges originating from the Pyrenean orogeny during the Late Cretaceous to Early Tertiary (figure 1).

The associated structures were subsequently reactivated as normal faults during the Oligocene, and as reverse strike-slip faults during the Miocene Alpine orogeny [START_REF] Roure | Geometry and kinematics of extensional structures in the alpine foreland basin of southeastern France[END_REF][START_REF] Baroux | Analyses of the stress field in southeastern France from earthquake focal mechanisms[END_REF][START_REF] Champion | La déformation post-Miocène en Provence occidentale[END_REF][START_REF] Chardon | Geological boundary conditions of the 1909 Lambesc (Provence, France) earthquake: Structure and evolution of the Trevaresse ridge anticline[END_REF][START_REF] Molliex | Tectonic and sedimentary inheritance on the structural framework of Provence (SE France): Importance of the Salon-Cavaillon fault[END_REF]. The latest deformation phase was followed by the uplift and rejuvenation of pre-existing ranges [START_REF] Siame | Local erosion rates versus active tectonics: cosmic ray exposure modelling in Provence (south-east France)[END_REF][START_REF] Bestani | Reconstruction of the Provence Chain evolution, southeastern France[END_REF]. However, various studies using diverse methodologies indicate that post-Miocene deformation rates are low in Provence [START_REF] Chardon | First paleoseismological constraints on the strongest earthquake in France (Provence) in the twentieth century[END_REF][START_REF] Cushing | A multidisciplinary study of a slow-slipping fault for seismic hazard assessment: the example of the Middle Durance Fault (SE France)[END_REF][START_REF] Nocquet | Present-day kinematics of the Mediterranean: A comprehensive overview of GPS results[END_REF]. The ranges are mostly made of uplifted lower Cretaceous carbonate units.

During the Late Miocene, the Provence area was heavily affected by the Messinian Salinity Crisis with a 1500 m drop of its Mediterranean base level (e.g. [START_REF] Clauzon | Le canyon messinien du Rhône : une preuve décisive du "dessicated deepbasin model[END_REF][START_REF] Gautier | Age et durée de la crise de salinité messinienne[END_REF][START_REF] Krijgsman | Chronology, causes and progression of the Messinian salinity crisis[END_REF][START_REF] Mocochain | The effect of river dynamics induced by the Messinian Salinity Crisis on karst landscape and caves: Example of the Lower Ardèche river (mid Rhône valley)[END_REF]. As described by these authors, this perturbation deeply imprinted the fluvial network, triggering intense regressive river erosion and the creation of deep canyons, which were later filled by marine and continental sediments once These sites were selected and sampled to provide a wide distribution of elevations and climatic settings, especially for Mean Annual Temperatures (ranging from 5 to 16°C) and Precipitation (ranging from 500 to 1400 mm/yr). Such diversity will allow to gain insights into the denudation pattern at the scale of a geologically coherent regional context and to specify the importance of climatic gradients in controlling the respective importance of physical and chemical weathering processes in carbonate settings.

-Methods

36 Cl concentrations were measured in carbonate samples from 25 sites distributed over 10 mountain ranges in Provence, South-Eastern France. The sampled ranges were, from North to South, the Ventoux and the Lure mountains, the Grand Luberon, the Alpilles, the Costes, the Sainte Victoire, the Sainte Baume, the La Fare mountain, the Nerthe and the Calanques. We also added 39 samples from the Petit Luberon crest to this dataset [START_REF] Godard | Weather-limited hillslope evolution in carbonate landscapes[END_REF][START_REF] Thomas | Morphological controls on the dynamics of carbonate landscapes under a mediterranean climate[END_REF] (figure 2). Plateau or flat surfaces on the crest of each range were selected to collect carbonate bedrock from roughly circular areas of approximately 10 meters in diameter using a chisel and a hammer (figure 3) with a sampling depth of 3-4 cm. The sampling sites were restricted to the crests of these ranges to avoid the influence of material flux from upslope. There is no topographic shielding in such sampling configuration. All samples were collected on bedrock surfaces. We selected sites with no soil cover and no evidence of former significant regolith development. Moreover, significant snow shielding is very unlikely at our sampling sites in Lure and Ventoux because snowfall is sporadic during the winter and these sites are located at narrow convex hilltops where snow is quickly removed by wind and downslope fluxes. Using a conservative value of 1 m of snow all over the winter duration yields a negligible shielding (>0.97, [START_REF] Gosse | Terrestrial in situ cosmogenic nuclides: theory and application[END_REF] 4). The sampling sites were specifically selected to encompass the largest possible range in MAT and MAP, which led to a wide range of elevation in the choice of sampling sites, from 200 meters at the Nerthe and the Alpilles ranges up to 1900 meters a.s.l at the Ventoux Mount (figure 4).

After collection, samples were crushed and sieved, and chlorine was chemically extracted using AgCl precipitation following the approach described by [START_REF] Schlagenhauf | Using in situ Chlorine-36 cosmonuclide to recover past earthquake histories on limestone normal fault scarps: A reappraisal of methodology and interpretations[END_REF] [START_REF] Sharma | Development of 36 Cl standards for AMS[END_REF][START_REF] Fifield | Accelerator mass spectrometry at the Australian National University's 14UD accelerator: experience and developments[END_REF].

In addition to denudation rates determination using cosmogenic nuclides, local curvature and slope at each sampling site were calculated using high-resolution Digital Elevation Models (DEM) with 1 meter resolution from the RGE ALTI database of the IGN (French National Institute of Geographical and Forest Information), in order to compare local topographic parameters and assess their influence on denudation.

We used WorldClim [START_REF] Fick | Worldclim 2: New 1-km spatial resolution climate surfaces for global land areas[END_REF] set of global climate rasters with a spatial resolution of ~1 km (figure 4) to link the obtained denudation rates with climatic parameters (figure 5) such as MAP, MAT and temperature amplitude. The WorldClim data was compared to the Météo France climate records at nine weather stations located throughout Provence to test its robustness (see figure 4A for location). A slight difference in MAT and a more significant deviation in MAP was observed (figure 5). The MétéoFrance data was used to determine correction parameters using a linear regression and to recalculate the climatic data at all the sampling sites.

-Results

The measured 36 Cl concentrations range from 5.94 x 10 5 to 3.36 x 10 6 atoms/g of rock, and the number of 36 Cl atoms in the samples range from ~6 x 10 6 to ~1.2 x 10 8 , which is 2 to 3 orders of magnitude higher than those measured in the simultaneously processed chemical blanks (10 5 atoms of 36 Cl). The amount of Clnat measured in the samples is also 10 times higher than the amount measured in the processed chemical blanks. Clnat concentrations are on average 8.6 ppm for the whole dataset, with all samples below 25 ppm except one reaching 81 ppm, which is consistent with natural Cl concentrations measured in other similar studies (e.g. [START_REF] Xu | In-situ cosmogenic 36Cl denudation rates of carbonates in Guizhou karst area[END_REF].

For our dataset the different production pathways of 36 Cl are, in average, as follows: spallation of target elements (78%), slow negative muons capture (19%) and thermal and epithermal neutron capture (3%) [START_REF] Schimmelpfennig | Calibration of cosmogenic 36 Cl production rates from Ca and K spallation in lava flows from Mt. Etna (38°N, Italy) and Payun Matru (36°S, Argentina)[END_REF]. The 36 Cl production rate from spallation of calcium has been calibrated in Sicily by [START_REF] Schimmelpfennig | Calibration of cosmogenic 36 Cl production rates from Ca and K spallation in lava flows from Mt. Etna (38°N, Italy) and Payun Matru (36°S, Argentina)[END_REF] on Ca-rich feldspar at 42.2 ± 4.8 atoms g/year. We used this production rate because the Provence area is located at almost the same latitude and elevation as this calibration site, and because it has been determined on mineral fractions. Scaling with respect to latitude and elevation was performed considering constant geomagnetic field intensity and [START_REF] Stone | Air pressure and cosmogenic isotope production[END_REF] scaling. Denudation rates were calculated using the approach of [START_REF] Schimmelpfennig | Sources of in-situ 36 Cl in basaltic rocks. Implications for calibration of production rates[END_REF]. The decay constant used (2.303 ± 0.016x10 6 year -1 ) corresponds to a 36 Cl half-life (T1/2) of 3.014x10 5 years. The analytical uncertainties include counting statistics, machine stability and blank correction. We note that calculating the denudation rates with CRONUS-Earth approach, as recently described in [START_REF] Marrero | Quaternary Geochronology Cosmogenic nuclide systematics and the CRONUScalc program[END_REF], would change the denudation rates values by at most 30% mostly because the authors used a different production rate. This won't affect our interpretation which is based on relative difference between the yielded denudation rate and thus independent of the chosen production rate or calculation scheme. However, this becomes important when comparing values with other studies and we choose [START_REF] Schimmelpfennig | Sources of in-situ 36 Cl in basaltic rocks. Implications for calibration of production rates[END_REF] calculator to be able to compare our data with previously published values which used the same calculation scheme, and also because the 36 Cl production rate from spallation of Ca, which is the most important production reaction for the samples in this study, was calibrated directly from mineral fractions in [START_REF] Schimmelpfennig | Calibration of cosmogenic 36 Cl production rates from Ca and K spallation in lava flows from Mt. Etna (38°N, Italy) and Payun Matru (36°S, Argentina)[END_REF].

Denudation rates across Provence deduced from our dataset range from 18.6 to 141.4 mm/ka including outliers (see Table 1 and figure 6A). The interquartile range is 30 to 50 mm/ka, which is in agreement with denudation rates determined on mountain crests over the Mediterranean area (e.g. Ryb et al., 2014a;[START_REF] Godard | Weather-limited hillslope evolution in carbonate landscapes[END_REF][START_REF] Thomas | Morphological controls on the dynamics of carbonate landscapes under a mediterranean climate[END_REF]. The time period over which the determined denudation rates are integrated range from approximately 5000 years, considering the highest denudation rates, and up to 32000 years considering the lowest denudation rates. The lowest denudation rates are determined in the lowermost ranges, such as the Nerthe and Calanques massifs (approximately 200 m), whereas the highest denudation rates are associated to the samples located at an altitude of 1800 m in the Lure mountain. A focus on each single range shows significant local-scale variability of the calculated denudation rates. For example, the Grand Luberon and Lure mountains display rates in the 55-117 and 111-141 mm/ka ranges, respectively. Smaller but noticeable variability is also observed at other sites, such as the Sainte Victoire (22-40 mm/ka), the Sainte Baume (34-51 mm/ka) and the Ventoux Mount (27-49 mm/ka). In the next section we develop an analysis of the observed variability in our dataset.

-Discussion

Our dataset provides the opportunity to discuss the dynamics and controls on carbonate surface denudation across a wide range of environmental conditions. We first discuss the observed variability between the different sites, as well as inside a same range, and we then focus on the influence of the tectonic and climatic boundary conditions. At last we attempt to put our data in a broader perspective by comparing with global denudation rates databases for both carbonates and quartz-rich landscapes.

-Variability of denudation rates in Provence

Our dataset provides extensive information on denudation rates from eleven mountain ranges in Provence, allowing a regional scale analysis of landscape evolution in a carbonate-dominated area. The determined denudation rates appear variable at the scale of both the individual ranges and the region. The former could be explained by local variations in the morphology at the sampling site (figure 6). Indeed, the analyzed samples have all been collected at the highest elevation of each range, but some were located on large plateau areas, such as in the Alpilles mountain, whereas others were collected on convex crests or on flat surfaces adjacent to cliffs (see figure 6D). As illustrated by data from the Petit Luberon, hilltop convexity is correlated with denudation rates determined at sampling locations on the crests [START_REF] Godard | Weather-limited hillslope evolution in carbonate landscapes[END_REF] and could explain a significant amount of variability of the denudation rates ranging from 30 to 70 mm/ka (figure 6D). Samples collected on a flat area but in the direct vicinity of a cliff display denudation rates that are in the 30-50 mm/ka range but have artificially high surface convexity due to the size of the computation window (e.g. Petit Luberon, Costes). Samples collected on relatively narrow crests can exhibit high denudation rates (e.g. Grand Luberon, Lure). We also note that the distance between sampling locations and the Mediterranean Sea does not clearly impact the denudation rates reported for individual ranges (figure 7).

Most of the samples, collected at various altitudes, display denudation rates clustering in the 30-60 mm/ka range despite contrasting climates ranging from the dry and hot coastal ranges of the Nerthe and Calanques massifs up to the highest, wettest and coldest sampling points at the top of the Ventoux Mount. This diversity of sampling sites encompasses a 3-fold and 2-fold changes in MAT and MAP, respectively. This range of denudation rates is consistent with data from other carbonate domains (e.g. [START_REF] Matsushi | Denudation rates of carbonate pinnacles in Japanese karst areas: estimates from cosmogenic 36Cl in calcite[END_REF][START_REF] Xu | In-situ cosmogenic 36Cl denudation rates of carbonates in Guizhou karst area[END_REF], including Mediterranean areas (Ryb et al., 2014a,b).

Samples from the Lure mountain and, to a lesser extent, from the Grand Luberon, whose denudation rates are comprised in the 111-141 mm/ka and 55-117 mm/ka ranges respectively, display significantly higher rates than the 30-60 mm/ka values reported for the other mountain ranges in our dataset as well as for carbonate landscapes elsewhere (e.g. [START_REF] Levenson | Comparison of field and laboratory weathering rates in carbonate rocks from an Eastern Mediterranean drainage basin[END_REF].

Considering the data distribution displayed in figure 6A, the reasoning developed in the previous paragraph could be reversed and thus the low denudation rates reported at high elevation for the Ventoux considered as outliers. In this case, the dataset including the samples from the Lure mountain and from the Grand Luberon mountain range would display a significant dependence of the denudation rates relative to their respective elevation. As previously exposed, our dataset has been acquired to question whether altitude and climatic conditions influence carbonate denudation or if, on the contrary, denudation remains constant over a wide range of climatic conditions, meaning that out of range denudation rates, such as those observed at the Lure mountain, are caused by other specific local parameters. Therefore, the interpretation of this dataset specifically hinges on the question of the signification of the large denudation variability observed along the Ventoux-Lure axis, which are areas displaying similar elevations as well as similar mean annual climatic parameters. To address this fundamental conundrum, the possible tectonic, climatic or geomorphic causes for the observed denudation contrast are discussed below.

-Importance of tectonic and climatic forcings

The Ventoux and Lure mountains are morphologically separated by the NNE Sault graben (George, 1952). This Oligocene geological feature could have induced a decoupling of the tectonic uplift between the eastern (Lure) and western (Ventoux) sides of the graben that could have had repercussions during the late Quaternary. According to [START_REF] Bestani | Reconstruction of the Provence Chain evolution, southeastern France[END_REF], this graben is linked to the deformation of the hanging wall of the Ventoux-Lure thrust system.

Considering that the Sault graben could have partially disconnected the Ventoux mountain from the Alpine deformation field, this could potentially explain the faster denudation on the Lure mountain because it is located closer to the Alps. This would have therefore led to an East-West rock uplift gradient. [START_REF] Bestani | Reconstruction of the Provence Chain evolution, southeastern France[END_REF] constructed two North-South cross sections to estimate the shortening resulting from the various tectonic phases that affected Provence. The first section crossing the Ventoux mountain shows a significantly smaller late-Cretaceous early-Tertiary shortening than the second section crossing the Lure mountain (7.8% vs. 26%). However, the post-Eocene shortening is equivalent (1.5 km) along the two cross sections. The similar crest elevations between the two mountains also plead against an important differential finite uplift accumulated over the Late Cenozoic.

To assess the possibility of a tectonic control over relevant time-scales considering the integration time of our denudation rates (~10 4 years), we investigated the variability of morphological proxies for tectonic uplift across the Ventoux-Lure system. In steady-state landscapes, a simple relationship exists between slope (S) and drainage area (A) measured along rivers [START_REF] Hack | Studies of Longitudinal Stream Profiles in Virginia and Maryland[END_REF][START_REF] Flint | Stream gradient as a function of order, magnitude, and discharge[END_REF][START_REF] Howard | Channel changes in badlands[END_REF][START_REF] Kirby | Expression of active tectonics in erosional landscapes[END_REF], which can be written as the slope-drainage area relationship:

S=ks A -θ
where ks is the channel steepness index and θ is the concavity. To allow easier inter-comparison between different river profiles, a normalized steepness index (ksn) can be calculated, using a fixed reference concavity (e.g. [START_REF] Kirby | Expression of active tectonics in erosional landscapes[END_REF]. Under the assumption of homogeneity of the erodibility of the substrate for the considered catchments, the ksn parameter can be directly linked to tectonic uplift in geomorphological studies [START_REF] Kirby | Distribution of active rock uplift along the eastern margin of the Tibetan Plateau: Inferences from bedrock channel longitudinal profiles[END_REF][START_REF] Wobus | Tectonics from topography: procedures, promise, and pitfalls[END_REF], and, in our case, a lateral variation of ksn along the Ventoux-Lure axis would evidence a gradient in rock uplift. A reference concavity of 0.5 was used to estimate averaged ksn values for catchments over several ranges, and we calculated ksn values for ~100 catchments from both the northern and southern flanks of the Ventoux and Lure ranges (figure 8).

The results obtained indicate that the northern flank of the Ventoux mountain displays relatively higher ksn values than the Lure mountain, which would imply a more rapid recent uplift of the Ventoux mountain and conflicts with the observed pattern of denudation rates. Similar observations can be made from the distribution of the average long-wavelength slopes of the flanks of both mountains (figure 8). Therefore, we conclude that a difference in Quaternary tectonic forcing along the Ventoux-Lure system is not likely to be a viable explanation for the observed denudation rate contrast between the two ranges.

The East-West trending Ventoux-Lure system is approximately 70 kilometers long and is located at the Northern limit of Provence, and at the boundary between two different climatic domains.

To the West, the Ventoux is under the influence of the Rhône valley climatic regime where the cool and dry Mistral wind dominates and results in low cloud cover throughout the year as well as scattered precipitations with sporadic storm events. To the East, the Lure mountain, located at the border of the Southern Alps region, is under the influence of the Alpine climate whose precipitations are more evenly distributed throughout the year.

The 36 Cl-derived denudation rates are averaged over approximately 5000 and 32000 years for the Lure and Ventoux mountain, respectively, and major climatic and geomorphic events occurring over these time intervals, such as the last deglaciation or Holocene climatic fluctuations, have to be considered in our interpretation of the denudation rates and could be one of the keys to understand their variability. In particular, the denudation rates calculated for the Ventoux integrate over the end of the last glacial period, whereas the denudation rates calculated at the Lure mountain only encompass the latter half of the Holocene.

High resolution DEM and orthophotos show a dense network of gullies incising the southern flank of the Lure mountain (figure 9) and suggest the recent development of widespread regressive erosion that almost reaches the crest, whereas no such geomorphic response can be observed near the summit of the Ventoux where our sampling site is located. A major postglacial valley sediment filling episode due to massive denudation of mountain flanks at the beginning of the Holocene has been recognized in the Southern Alps [START_REF] Miramont | The "main postglacial" sedimentation cycle of marly subalpine basins : an inventory[END_REF] and may explain the large regressive erosion observed on the Lure mountain. One consequence of such major hillslope destabilisation would have been to increase denudation rates resulting in present-day significantly lower 36 Cl concentrations. In this case the samples from the Lure mountain summit would be out-of-equilibrium with respect to 36 Cl accumulation and could not be directly interpreted in terms of denudation rates. Due to its position to the West, the Ventoux is disconnected from the Durance river base level and associated Alpine influence and, for that reason, was not affected by the erosion event. Most paleoenvironmental studies focusing on 14 C or pollen analysis to reconstruct past environments (e.g. [START_REF] Thinon | La végétation du mont Ventoux au cours des derniers millénaires[END_REF] suggest that surface conditions and vegetation cover has remained stable throughout the first half of the Holocene on the Ventoux, even though agricultural activities may have induced a change in vegetation during the last 4000 years mainly on the summit area. Additionally, the Durance glacier was located directly East of the Lure range and had a major impact on the local morphology until at least 20 ka BP [START_REF] Rosique | The last glacial phase in the middle Durance (region from Gap to Sisteron): research appraisal, latest chronological conclusions[END_REF], whereas the Ventoux was not directly affected by these major motion rates and gradients in the external parts of the Alps [START_REF] Champagnac | Quaternary erosioninduced isostatic rebound in the western Alps[END_REF], which could influence present-day denudation patterns.

-Global perspective on carbonate landscapes denudation

The dataset acquired during this study was compared to a global compilation of denudation rates derived from 36 Cl concentrations measured in carbonates (Levenson et al., 2016) that spans a broad spectrum of climatic settings, with MAP ranging from almost 0 to 2500 mm (figure 10).

Our dataset encompasses the intermediate 500 to 1500 mm MAP domain. The reported denudation rates are mostly within the 20-60 mm/ka range which is compatible with observations from other settings experiencing comparable or wetter climatic conditions. Notable deviations are (1) the high-convexity crests of the Petit Luberon where denudation is accelerated by regolith mass transport at the surface [START_REF] Godard | Weather-limited hillslope evolution in carbonate landscapes[END_REF] and reaches 70-80 mm/ka and (2)

the Lure mountain summit whose particular geomorphic history has been discussed above. The lack of any significant correlation with MAP in our dataset is consistent with results from other non-arid regions, where water availability does not act as a limiting factor. Samples from drier regions display significantly lower denudation rates (<20 mm/ka, Ryb et al., 2014a,b), and this contrast supports the proposition by Ryb et al (2014a) that, for steep hillslopes, a fundamental MAP threshold exists around 100-300 mm, below which mechanical weathering processes become the main drivers of landscape denudation.

Many insights in the geomorphological evolution of carbonates rocks and landscapes have been derived from other methodologies. In another Mediterranean area, [START_REF] Furlani | Comparison between coastal and inland Karst limestone lowering rates in the northeastern Adriatic Region (Italy and Croatia)[END_REF] have investigated the lowering rates of limestone surfaces using micro erosion measurements at various sites in the Adriatic region. For inland karst measurement stations, which are not subjected to coastal geomorphological forcings, they report lowering rates of 0.018 mm/yr, and as low as 0.009 mm/yr for the inland Istrian Karst. These surface lowering rates are significantly less than denudation rates derived from 36 Cl concentrations measurements under comparable climatic conditions (figure 10), but consistent with low rates derived from similar methodologies or limestone tablets dissolution experiments (e.g. [START_REF] Plan | Factors controlling carbonate dissolution rates quantified in a field test in the Austrian alps[END_REF]Gabrovsek, 2009 and references therein). It is to be noted that these experiments integrate weathering processes over at most a few decades, and thus do not incorporate major climatic and environmental perturbations occurring on Holocene timescales such as the regressive erosion event we considered to discuss the apparent high denudation rates at the top of the Lure mountain. Furthermore, due to the experimental strategy of most of these studies, they typically measure the effects of surface lowering or mass loss due mainly to chemical weathering, whereas the denudation signal obtained through cosmogenic nuclides measurements is by definition the sum of chemical and mechanical processes. Such difference could be one of the reasons for the lower rates reported by these experiments and would further highlight the importance of mechanical weathering agents breaking down bedrock mass and producing sediment clasts available for transport in the evolution of carbonate landscapes [START_REF] Thomas | Morphological controls on the dynamics of carbonate landscapes under a mediterranean climate[END_REF].

We also attempt to put the question of an eventual climatic control on carbonate landscapes into a broader context by comparing the existing 36 Cl denudation rates with a global compilation of bedrock 10 Be denudation rates in quartz-rich settings [START_REF] Portenga | Understanding earth's eroding surface with 10 Be[END_REF], figure 11).

The broad range of climatic and geologic settings investigated as well as the non-constrained influence of site-specific morphological contexts (slope, curvature) are probably major contributors to the observed high variability over given MAP intervals. Such a large data dispersion is inherent to these global compilations, and we acknowledge the associated discussion should be treated cautiously and as a general investigation of the first order trends.

Nevertheless, we observe two important features that we think worthy of consideration in this comparison of the two datasets. First, even if the datasets display numerous overlaps, when considering average values over similar ranges of MAP the 36 Cl derived denudation rates in carbonate settings are systematically at least 2 times higher than the 10 Be derived denudation rates in quartz-rich settings. This difference would imply significantly shorter residence times of carbonate terrains at the surface of the Earth with respect to other types of bedrock. Second, the exponents of the general power-law trends extracted from both datasets are similar (~0.5) and suggest a weak dependence of surface denudation on the variability of precipitation over the range of investigated MAP. Furthermore, for both datasets, the existence of this shallow trend appears to be mostly due to very low denudation in arid regions, whereas results from temperate to humid regions show no apparent correlation and suggest that water availability is not a limiting factor on denudation over a large range of precipitation conditions at the Earth surface.

Assessing the extent of climatic controls, mainly through precipitation, on landscape evolution is a long-standing issue in geomorphology, and this general result is consistent with independent observations in specific climatic and tectonic settings (e.g. [START_REF] Riebe | Strong tectonic and weak climatic control of long-term chemical weathering rates[END_REF][START_REF] Godard | Dominance of tectonics over climate in Himalayan denudation[END_REF][START_REF] Gasparini | Diagnosing climatic and tectonic controls on topography: Eastern flank of the northern Bolivian Andes[END_REF]. It is noteworthy that most of these studies downplaying the role of climate on denudation focused on silicate-rich terrains. The fact that such a result appears to hold in non-arid carbonate terrains, even though these domains are often considered to be dominated by chemical weathering, provides further support to the idea of a limited impact of precipitation gradients on landscape denudation at 1-100 ka timescales.

-Conclusion

We collected a well distributed dataset documenting denudation patterns across a carbonatedominated region. The denudation rates deduced from 36 Cl concentrations in our 25 samples and the 39 samples from the Petit Luberon (Godard et al;, 2016 and[START_REF] Thomas | Morphological controls on the dynamics of carbonate landscapes under a mediterranean climate[END_REF] are consistent with the range of denudation rates observed in other carbonate domains. We note however a discrepancy in terms of denudation rates between the Ventoux (30-40 mm/ka) and the Lure (110-140 mm/ka) mountains. It is possible that the apparent decoupling in behavior between these two ranges is of tectonic origin or results from strong morphological events following the end of the last glacial period. Nevertheless, this study is one of the first to document 36 Cl-computed denudation rates for carbonate rocks across a wide array of altitudes and variable mean annual temperatures and precipitations. Calculated denudation rates appear to be strongly impacted by changes in surface processes dynamics induced by climatic variability at the scale of the Holocene. On the other hand, the regional distribution of denudation rates suggests a low sensitivity to climatic spatial gradients over the range of conditions investigated here. We focused here on MAP as the main climatic variable to evaluate the response of carbonate surfaces, however subsequent studies should also focus on multifactorial environmental controls, including temperature [START_REF] Marshall | Frost for the trees : Did climate increase erosion in unglaciated landscapes during the late Pleistocene ?[END_REF] but also the complex influence of variations in vegetation (e.g. Torres [START_REF] Torres Acosta | Effect of vegetation cover on millennial-scale landscape denudation rates in East Africa[END_REF][START_REF] Schaller | Addressing the contribution of climate and vegetation cover on hillslope denudation, Chilean Coastal Cordillera (26-38S)[END_REF].

Indeed, the influence of vegetation on CO2 partial pressure in the subsurface is potentially an important control on carbonate denudation (e.g. [START_REF] Calmels | Sensitivity of carbonate weathering to soil CO2 production by biological activity along a temperate climate transect[END_REF]. Black circles are 36 Cl denudation rates from samples collected in this area.

  the elevation in western Provence progressively increases from South to North. The crests of the main carbonate ranges (figure1) reach altitudes ranging from 200 to 1900 meters, which induces a significant diversity in the climatic conditions. The Ventoux-Lure mountain marks the geological limit of Provence to the North. It is the highest mountain range in the region, reaching 1920 and 1800 m at the Ventoux and Lure mountains, respectively. Forty kilometers South of the Ventoux-Lure range, the Luberon and Alpilles mountains form another major EW system. From East to West, the Luberon range is composed of the Manosque fold (altitude: ~800 meters), the Grand Luberon (altitude: ~1100 meters) and the Petit Luberon (altitude: ~700 meters). The Alpilles, to the West, are separated from the Luberon by the North-South Salon-Cavaillon fault, and rise to lower altitudes (470 m), with a well-developed plateau surface. Further East, the major features defining Provence landscapes are the Sainte Victoire and Sainte Baume mountains, both higher than 1000 meters. Smaller ranges also contribute to Provence landscapes such as the Costes (450 meters) and La Fare (200 meters) ranges and, bordering the Mediterranean coast, the Nerthe (250 meters) and the Calanques (450 meters) massifs, respectively located to the West and to the East of Marseille.

  perturbations occurring in the upper and middle Durance valley. Collectively, these observations support the idea that the observed contrast in denudation rates and 36 Cl concentrations along the Ventoux-Lure System results from a major climatically-triggered geomorphological change in the Durance valley. A significant change in hillslope dynamics and associated denudation could have been triggered by two nonexclusive types of perturbations : (1) an abrupt fall in the base level of the gully network incising the northern flank of the Lure range, possibly related to the retreat of the Durance glacier, and (2) a drastic modification of the regolith transport parameters, associated with changes in temperature, precipitation and land cover. The melting of the Alpine ice cap at the end of the last deglaciation can also explain the observed high geodetic vertical

Figure 1 -

 1 Figure 1 -Geological map of South-Eastern France (modified from[START_REF] Molliex | Tectonic and sedimentary inheritance on the structural framework of Provence (SE France): Importance of the Salon-Cavaillon fault[END_REF]

Figure 2 -

 2 Figure 2 -Upper-left inset: general situation of the sampled sites (red circles) in Western

Figure 3 -

 3 Figure 3 -Field photographs of selected sampling sites. A-B-C: sampling locations on flat

Figure 4 -

 4 Figure 4 -A) Topographic map of Provence displaying the location of Météo France weather

Figure 5 -

 5 Figure 5 -A) Mean Annual Temperature (MAT, in °C) at each sampling site according to

Figure 6 -

 6 Figure 6 -Denudation rates for bedrock samples compared to various topographic and climatic

Figure 7 -

 7 Figure 7 -A) Western topographic profile (see Figure 1 for location) with 36 Cl denudation rates.

Figure 8 -

 8 Figure 8 -A) Topographic map of the Ventoux-Lure range displaying the location of samples

Figure 9 -

 9 Figure 9 -A) Shaded relief map (25 m IGN DEM) centered on the Ventoux summit area. B)

Figure 10 -

 10 Figure 10 -Denudation rates obtained in Provence (this study, Thomas et al. (2017), Godard et

Figure 11 -

 11 Figure11-36 Cl carbonate denudation rates from this study andLevenson et al. (2016)

  . Samples lithologies were homogeneous and consisted of consolidated and massive Cretaceous and Jurassic limestones. The sampled sites have Mean

Annual Temperature (MAT) ranging from 5°C at the Ventoux and the Lure mountains to 16°C at the locations closest to the Mediterranean Sea (table

1

). Mean Annual Precipitation (MAP) are roughly inversely correlated to MAT, with approximately 500-600 mm at seaside locations and up to 1400 mm over the Ventoux and Lure mountains (figure

  . Natural chlorine (hereafter Clnat) and 36 Cl concentrations were then determined via isotope dilution accelerator mass spectrometry based on the measured 35 Cl/ 37 Cl and 36 Cl/ 35 Cl ratio measurements

	performed at the French national accelerator mass spectrometry (AMS) facility ASTER
	(Accélérateur pour les Sciences de la Terre, Environnement, Risques) at CEREGE, Aix-en-
	Provence, France (Arnold et al., 2013). The 36 Cl concentrations were normalised to the 36 Cl
	standard KNSTD1600 with a nominal 36 Cl/ 35 Cl value of 2.11 ± 0.06 10 12

Table Sample

 Sample

	ID	Mountain range	Type	Longitude Latitude Altitude	[ 36 Cl]	[Cl]	[Ca] Denudation rate MAT MAP Slope Curvature
				°E	°N	m	10 5 atom/g	ppm	%	mm/ka	°C	mm/yr	°	m-1
							13.22 ±
	NE2	Nerthe	Bedrock	5.2403	43.3410	211	0.49	24.25 40.5	29 ± 3.1	14.4	558	2.39	0.011
							16.40 ±
	NE3	Nerthe	Bedrock	5.1596	43.3562	163	0.58	9.61	38.7	18.6 ± 2	14.7	541	3.27	0.003
							11.15 ±
	CA1	Calanques	Bedrock	5.4363	43.2200	230	0.37	13.96 40.8	34.4 ± 3.6	14.8	573	6.86	0.025
							11.64 ±
	CA3	Calanques	Bedrock	5.4638	43.2321	382	0.45	10.17 39.7	35.4 ± 3.8	13.9	652	15.44	0.027
	AL1	Alpilles	Bedrock	5.0245	43.7561	260	16.37 ± 0.5	4.59	47.1	21.6 ± 2.3	13.6	668	2.79	0.002
							10.33 ±
	AL4	Alpilles	Bedrock	5.0318	43.7620	240	0.54	7.14	38.3	35.8 ± 4	13.8	657	3.39	0.011
							12.12 ±
	SB2	Sainte Baume Bedrock	5.7332	43.3187	984	0.49	3.81	40.8	50.6 ± 5.5	12.2	827	13.10	0.057
	SB5	Sainte Baume Bedrock	5.6979	43.3143	997	17.19 ± 0.6	7.31	40	34.2 ± 3.6	10.9	927	0.79	0.011
							19.76 ±
	SB7	Sainte Baume Bedrock	5.6823	43.3189	865	0.84	80.83 40.2	39.4 ± 4.3	11.6	860	0.84	0.006
	GL1	Grand Luberon Bedrock	5.5158	43.8192	1051	8.83 ± 0.41	25.6	39.7	85.5 ± 9.4	9.7	1004	3.03	0.011
	GL3	Grand Luberon Bedrock	5.5277	43.8196	1049	6.58 ± 0.31 22.81 40.3	116.5 ± 12.9	9.7	1008	7.35	0.030
							12.17 ±
	GL4	Grand Luberon Bedrock	5.5361	43.8182	1032	0.48	12.06 39.6	54.6 ± 5.9	10.3	944	1.96	0.019
							11.54 ±
	LF1	La Fare	Bedrock	5.2110	43.5625	219	0.44	3.78	40	30.5 ± 3.3	14.1	580	2.41	0.007
	LF2	La Fare	Bedrock	5.1963	43.5613	192	8.93 ± 0.36	6.58	40.2	41.7 ± 4.5	14.2	562	10.53	0.009
							16.49 ±
	CO1	Costes	Bedrock	5.2725	43.6937	477	0.63	3.64	40.9	23.8 ± 2.5	12.9	709	10.22	0.063
							21.27 ±
	VE1	Ventoux	Bedrock	5.2678	44.1743	1835	0.76	12.1	39.9	48.9 ± 5.2	5.1	1484	3.44	0.064
							33.61 ±
	VE2	Ventoux	Bedrock	5.2633	44.1740	1838	1.16	4.66	39.8	27 ± 2.9	5.5	1434 13.09	0.063
							27.27 ±
	VE4	Ventoux	Bedrock	5.2749	44.1730	1851	0.95	7.38	39.6	35.7 ± 3.8	5.1	1484 11.41	0.052
							22.60 ±
	VE5	Ventoux	Bedrock	5.2604	44.1851	1440	0.91	16.21 38.4	35.5 ± 3.8	7.8	1206	8.35	0.036
	LU1	Lure	Bedrock	5.7853	44.1225	1725	9.22 ± 0.46	16.8	40.4	119.1 ± 13.3	6.2	1377 14.77	0.025
	LU2	Lure	Bedrock	5.7887	44.1231	1751	9.79 ± 0.48 14.14 39.4	111.1 ± 12.4	6.2	1377 12.51	0.031
	LU3	Lure	Bedrock	5.8000	44.1232	1795	8.30 ± 0.43 20.06 40.3	141.4 ± 15.9	5.9	1410 11.70	0.031
							21.78 ±
	SV1	Sainte Victoire Bedrock	5.6182	43.5321	1012	0.55	13.32 39.2	26.5 ± 2.7	11.8	822	2.74	0.004
							24.31 ±
	SV3	Sainte Victoire Bedrock	5.6081	43.5310	1004	0.65	6.76	41	22.3 ± 2.3	10.4	947	2.13	0.004
							16.31 ±
	SV4	Sainte Victoire Bedrock	5.6463	43.5391	1000	0.51	22.45 39.8	40.3 ± 4.2	8.3	927	5.07	0.064

Table 1 -

 1 Samples for in-situ cosmogenic 36 Cl measurement from Provence, SE France, with

Acknowledgments

This study was supported by the ECCOREV Research Federation, INSU-CNRS and the Cashima Research Project (funded by the CEA Cadarache and ITER Organisation). Franck Thomas benefitted from a PhD scholarship funded by the CEA and the Conseil Regional PACA. We thank the Parc Naturel Régional du Luberon and the Office National des Forêts and the Plan d'Aups city for facilitating access to our field sites. The 36Cl measurements were performed at the ASTER AMS national facility (CEREGE, Aix en Provence) which is supported by the INSU/CNRS, the ANR through the "Projets thematiques d'excellence" programme for the "Equipements d'excellence" ASTER-CEREGE action, IRD and CEA. This work is also a contribution to the Labex OT-Med (n ANR-11-LABX-0061) funded by the French Government "Investissements d'Avenir" programme of the French National Research Agency (ANR) through the A*MIDEX project (n ANR-11-IDEX-0001-02).