Léonel Fouédjeu
email: leonel.fouedjeu-foumou@univ-tlse2.fr

Mélanie Saulnier

Mathieu Lejay

Martin Dušátko

Vincent Labbas

Alistair S Jump

Sylvain Burri

Sandrine Buscaino

Vanessa Py
email: vanessa.py@univ-tlse2.fr

Vanessa Py-Saragaglia

High resolution reconstruction of modern charcoal production kilns: an integrated approach combining dendrochronology, micromorphology and anthracology in the French Pyrenees

Keywords: ancient charcoal production, anthracology, micromorphology, dendrochronology, Northern Pyrenees, modern radiocarbon plateau

come   L'archive ouverte pluridisciplinaire

Introduction

Worldwide, many Charcoal Kiln (CK) terrace or platform remains still exist in present day open landscapes (fields and meadows), forests and even supra-forest areas above the current timberline, providing evidence of the wide distribution of the practice of charcoal production in the past [START_REF] Hart | Legacy of Charcoaling in a Western Highland Rim Forest in Tennessee[END_REF][START_REF] Ludemann | Past fuel wood exploitation and natural forest vegetation in the Black Forest, the Vosges and neighbouring regions in western Central Europe[END_REF][START_REF] Dupin | Anthracology of charcoal kilns in the forest of Chailluz (France) as a tool to understand Franche-Comte forestry from the mid-15th to the early 20th century AD[END_REF][START_REF] Benatti | Charcoal Kilns in the Northern Apennines (Italy): Forest Exploitation by Past Societies in Mountain Areas[END_REF][START_REF] De Oliveira | A hidden landscape, the paleo-territory of charcoal makers in Rio de Janeiro, Brazil[END_REF][START_REF] Durand | Le charbonnage en fosse : approche ethnoarchéologique d'une pratique méconnue dans le Rif marocain[END_REF][START_REF] Patzlaff | Historical charcoal kilns. A method to compare the surrounding vegetation with the anthracological data in forested slopes of Rio de Janeiro Estate, RJ, Brazil[END_REF] Indeed, charcoal was and in some locations still is a crafting resource used as fuel for ironmaking and smelting of other metals. Studying of the spatial distribution and chronology of former CK terraces [START_REF] Raab | Preindustrial charcoal production in Lower Lusatia (Brandenburg, Germany): Detection and evaluation of a large charcoal-burning field by combining archaeological studies, GIS-based analyses of shaded-relief maps and dendrochronological age determination[END_REF]Deforce et al., 2020), as well as their taxonomic composition [START_REF] Davasse | Forêts, charbonniers et paysans dans les Pyrénées de l'est du moyen âge á nos jours: une approche géographique de l'histoire de l'environnement[END_REF][START_REF] Rouaud | Les forêts de pente de la haute vallée de la Dordogne : enjeux écologiques et énergétiques d'une ancienne forêt charbonnée[END_REF][START_REF] Dupin | Anthracology of charcoal kilns in the forest of Chailluz (France) as a tool to understand Franche-Comte forestry from the mid-15th to the early 20th century AD[END_REF] or diameter of charred woods [START_REF] Ludemann | Die Wälder am Schauinsland und ihre Nutzung durch Bergbau und Köhlerei[END_REF][START_REF] Knapp | Charcoal usage in medieval and modern times in the Harz Mountains Area, Central Germany: Wood selection and fast overexploitation of the woodlands[END_REF][START_REF] Paradis-Grenouillet | Sustainable management of metallurgical forest on Mont Lozère (France) during the Early Middle Ages[END_REF][START_REF] Tolksdorf | Forest exploitation for charcoal production and timber since the 12th century in an intact medieval mining site in the Niederpöbel Valley (Erzgebirge, Eastern Germany)[END_REF], has already provided crucial insights into the history of pre-industrial charcoal production and its associated forest management systems (silvicultural treatments), as well as enabling the characterisation of their legacy in current landscapes [START_REF] Samojlik | Quantifying Historic Human Impacts on Forest Environments: A Case Study in Bialowieza Forest, Poland[END_REF][START_REF] Schmidt | Charcoal kiln sites, associated landscape attributes and historic forest conditions: DTM-based investigations in Hesse (Germany)[END_REF][START_REF] Pescini | Multi-proxy record of environmental changes and past land use practices in a Mediterranean landscape: The Punta Mesco cape (Liguria -Italy) between the 15th and 20th century[END_REF]. In addition, approaches combining historical, archaeological and ethnographic data have also enabled assessment of the art of burning charcoal, and especially its technical operational sequence (chaîne opératoire) and its seasonality [START_REF] Burri | Enquête ethnoarchéologique sur le mode de vie et le savoir-faire des derniers charbonniers de Calabre (Italie) : entre tradition et modernité[END_REF]2018;Paradis-Grenouillet et al., 2018).

However, in the northern hemisphere, the study of CK remains, which testify to a major development of charcoal manufacturing in the last centuries, i.e. between the seventeenth and nineteenth century [START_REF] Davasse | Forêts, charbonniers et paysans dans les Pyrénées de l'est du moyen âge á nos jours: une approche géographique de l'histoire de l'environnement[END_REF][START_REF] Allée | L'exploitation médiévale du plomb argentifère sur le mont Lozère: Archéologie spatiale d'un territoire proto-industriel montagnard[END_REF][START_REF] Ludemann | Past fuel wood exploitation and natural forest vegetation in the Black Forest, the Vosges and neighbouring regions in western Central Europe[END_REF][START_REF] Strachan | Application of Dendrochronology to Historical Charcoal-Production Sites in the Great Basin, United States[END_REF][START_REF] Raab | Preindustrial charcoal production in Lower Lusatia (Brandenburg, Germany): Detection and evaluation of a large charcoal-burning field by combining archaeological studies, GIS-based analyses of shaded-relief maps and dendrochronological age determination[END_REF][START_REF] Dupin | Anthracology of charcoal kilns in the forest of Chailluz (France) as a tool to understand Franche-Comte forestry from the mid-15th to the early 20th century AD[END_REF][START_REF] Karimi Moayed | Bypassing the Suess-effect: Age determination of charcoal kiln remains using OSL dating[END_REF] always faces the major issue that is the weak accuracy of dating obtained since the 1650s [START_REF] Davasse | Forêts, charbonniers et paysans dans les Pyrénées de l'est du moyen âge á nos jours: une approche géographique de l'histoire de l'environnement[END_REF][START_REF] Strachan | Application of Dendrochronology to Historical Charcoal-Production Sites in the Great Basin, United States[END_REF][START_REF] Py-Saragaglia | Late Holocene history of woodland dynamics and wood use in an ancient mining area of the Pyrenees (Ariège, France)[END_REF][START_REF] Deforce | 2500 years of charcoal production in the Low Countries: The chronology and typology of charcoal kilns and their relation with early iron production[END_REF][START_REF] Karimi Moayed | Bypassing the Suess-effect: Age determination of charcoal kiln remains using OSL dating[END_REF]. Poor accuracy leads to strong wiggles in the calibration curve creating a "radiocarbon plateau", which are aggravated first by the "De vries effect", i.e. high-frequency preindustrial secular variations in 14C activity and secondly, by the "Suesseffect", i.e. the dilution of atmospheric 14C levels resulting from the widespread burning of fossil fuels since the beginning of the 1900s (coal, gas, oil) [START_REF] De Vries | Variation in concentration of radiocarbon with time and location on Earth[END_REF][START_REF] De Vries | Measurement and use of natural radiocarbon[END_REF][START_REF] Stuiver | Variations in radiocarbon concentration and sunspot activity[END_REF][START_REF] Tans | Natural atmospheric 14C variation and the Suess effect[END_REF].

New experiments are increasingly implemented to refine the chronology of modern charcoal production. Optically Stimulated Luminescence (OSL) method, for example, has been recently successfully applied to refine the dating of modern CKs in Belgium [START_REF] Karimi Moayed | Bypassing the Suess-effect: Age determination of charcoal kiln remains using OSL dating[END_REF]. Other testing approaches are based on dendrochronological analysis applied on remaining charcoals from CKs [START_REF] Nelle | Charcoal burning remains and forest stand structure -Examples from the Black Forest (southwest Germany) and the Bavarian Forest (south-east Germany)[END_REF][START_REF] Strachan | Dendrochronological Dating of Historic Charcoal Oven Platforms in the Great Basin, North America[END_REF][START_REF] Strachan | Application of Dendrochronology to Historical Charcoal-Production Sites in the Great Basin, United States[END_REF][START_REF] Raab | Preindustrial charcoal production in Lower Lusatia (Brandenburg, Germany): Detection and evaluation of a large charcoal-burning field by combining archaeological studies, GIS-based analyses of shaded-relief maps and dendrochronological age determination[END_REF]. Dendrochronological study requires charcoal pieces that preserve a sufficient number of rings (i.e.>30) to enable robust cross-dating [START_REF] Backmeroff | Dendrochronological dating of charcoal kilns: a new method for dating historical land use at the upper timberline[END_REF][START_REF] Pichler | Dendrochronological dating of charcoal originating from an Early Iron Age fire-set pit in the mining area of Schwaz/Brixlegg (Tyrol, Austria)[END_REF][START_REF] Rouaud | Les forêts de pente de la haute vallée de la Dordogne : enjeux écologiques et énergétiques d'une ancienne forêt charbonnée[END_REF]. However, even the largest pieces found in CKs often preserve a low number of annual growth rings, because of charcoal fragmentation during carbonisation and the post-depositional processes. In order to increase the length of ring chronologies, it is therefore necessary to sample a higher number of charcoal fragments [START_REF] Blondel | Impact of carbonization on growth rings: Dating by dendrochronology experiments on oak charcoals collected from archaeological sites[END_REF]. Moreover, the progress of formation of terminal tree-ring (early wood, late wood) on charcoal fragments that kept their bark is also highly valuable to identify the felling season and thus reconstruct the seasonality of wood harvesting practices [START_REF] Py-Saragaglia | Late Holocene history of woodland dynamics and wood use in an ancient mining area of the Pyrenees (Ariège, France)[END_REF][START_REF] Fouédjeu | Charcoal analysis of kiln terraces to reconstruct practices and wood uses' temporalities for charcoal manufacturing activity : the case of the Bernadouze forest[END_REF].

Micromorphology is often used to characterise depositional history of archaeological sites [START_REF] Béguier | Étude micromorphologique du Chemin Saulnier : apport à l'archéologie des paysages de la Plaine de Caen (Calvados, Basse-Normandie)[END_REF]Nicosia and Stoops, 2017) but has been very rarely applied to the study of CK sites [START_REF] Hirsch | Soils on Historic Charcoal Hearths: Terminology and Chemical Properties[END_REF]Raab et al., 2017). Micromorphological analysis can provide valuable and complementary information about CK terrace by characterising their micro-stratigraphy. Such analyses may highlight different phases of reuse, which are often difficult to detect by standard archaeological and anthracological methods. Indeed, it is difficult to clearly identify reuse phases of a terrace stratigraphically (the stratigraphic succession of charcoal-rich levels and levels of surrounding sediment leaching following periods of long abandonment) although several historical sources [START_REF] Duhamel Du Monceau | Art du charbonnier, ou Manière de faire le charbon de bois. Description des arts et métiers, faits ou approuvés par Messieurs de l'Académie royale des Sciences[END_REF][START_REF] Marié De L'isle | Manuel complet des marchands de bois et de charbons, ou Traité de ce commerce en général[END_REF] and ethnographic surveys testify to such sequences [START_REF] Battesti | Une charbonnière expérimentale, Méthodologie de l'anthracologie sur charbonnière pour une approche de la gestion du milieu forestier en Languedoc[END_REF][START_REF] Lugli | Modelli e finalità nello studio della produzione di carbone di legna in archeologia[END_REF][START_REF] Montanari | The study of charcoal-burning sites in the Apennine mountains of Liguria (NW Italy) as a tool for forest history[END_REF][START_REF] Scipioni | Dati antracologici da due pratiche attuali : la produzione del carbone vegetale nelle alte valli Scrivia e Trebbia (Appennino Ligure)[END_REF]. CK reuse is driven by constraints that are both topographic, i.e. the low availability of flat areas in mountainous landscape that leads to the reuse of existing sites, and technical, i.e the mixture of charcoals and earth originating from previous CK burning is considered to be better to cover the kilns than newly excavated soil [START_REF] Poggi | Carbonai dell'Antola, video-cassetta 40[END_REF][START_REF] Burri | Enquête ethnoarchéologique sur le mode de vie et le savoir-faire des derniers charbonniers de Calabre (Italie) : entre tradition et modernité[END_REF]. The terrace was usually cleaned before each use and, then, the removed soil was recycled to cover the new wood pile, leading to the mixture of old and new charcoals during each carbonisation. Despite complex reuse processes, the terrace stratigraphy (also termed the charcoal platform soil profile) is often visually homogeneous (forming a continuous black layer) with thickness usually varying between 5 to 30 cm, and sometimes up to 80-100 cm [START_REF] Bonhôte | Charcoal kilns and environmental history in the eastern Pyrenees (France). A methodological approach[END_REF][START_REF] Dupin | Anthracology of charcoal kilns in the forest of Chailluz (France) as a tool to understand Franche-Comte forestry from the mid-15th to the early 20th century AD[END_REF]. In France, only [START_REF] Gebhardt | Impact of charcoal production activities on soil profiles: the micromorphological point of view[END_REF] was interested in stratigraphic analysis based on a micromorphological approach until more recent applications by [START_REF] Dupin | The taphonomic characterization of a charcoal production platform. Contribution of an innovative pair of methods: Raman analysis and micromorphology[END_REF]. [START_REF] Dupin | The taphonomic characterization of a charcoal production platform. Contribution of an innovative pair of methods: Raman analysis and micromorphology[END_REF] highlighted that stratified CKs usually result from an occasional use of the terrace while CKs characterised by a "homogeneous" profile are usually linked to periods of repeated intensive charcoal production. This methodological paper aims to evaluate the relevance of an interdisciplinary approach combining archaeological, micromorphological and dendrochronological analyses in order to better-assess the chronology of modern charcoal production. More specifically, our objectives are to detect and date, firstly, potential different phases of charcoal production, i.e. several years of CK uses, and secondly, within these phases, the different episodes, i.e. a single or several burnings within a short period of time. Consequently, we decided to apply this original approach on an identified CK exploited with certainty during the last century (in the 1900s). This CK is located in the Northern Pyrenean beech forest known as "Forêt de Bernadouze", within the Vicdessos steel valley, where charcoal production is well-documented by previous historical ecology studies [START_REF] Izard | Historique et dynamique du paysage forestier de la vallée de Suc et Sentenac (Vicdessos, Ariège)[END_REF][START_REF] Py-Saragaglia | Chapter 10: Long-term forest evolution and woodland uses in an ancient charcoal-production forest of the eastern French Pyrenees: an interdisciplinary approach at high spatio-temporal resolution[END_REF][START_REF] Saulnier | A study of late Holocene local vegetation dynamics and responses to land use changes in an ancient charcoal making woodland in the central Pyrenees (Ariège, France), using pedoanthracology[END_REF]. The importance of working with a recently exploited CK terrasse is justified by the need (i) to compare oral testimony with the results of this study and validate them, (ii) to obtain charcoal fragments whose ring chronologies can be cross-dated with a reference chronology based on living beech trees. Micromorphological analysis has been applied to detect whether the CK terrace has been reused for several years or even decades. In addition, dendrochronological analysis has been applied on charcoal fragments in order to accurately date the CK use phases. In order to cross-date the master charcoal ring width chronology, we first construct a robust chronology for living beech trees in the "Forêt de Bernadouze", considering that no beech chronology was currently available for the Vicdessos valley or its surrounding area. Subsequently, we analysed terminal growth rings of charcoal fragments that retained their bark, in order to identify the felling season and assess CK intensity of use. Beyond these strictly methodological aspects, this study raises other issues related to charcoal making practices, i.e. the ways of using a terrace and associated fuel wood logging rhythms, and their impacts on the stratigraphy of the terrace and on silvicultural practices (coppice and rotation).

Material and methods

Study area

The study area is located on the northern slope of the Pyrenees, in the high valley of Vicdessos (administrative department of Ariège), on ancient limestone rocks between 1300 to 1450 m a.s.l (Fig. 1) (see [START_REF] Saulnier | A study of late Holocene local vegetation dynamics and responses to land use changes in an ancient charcoal making woodland in the central Pyrenees (Ariège, France), using pedoanthracology[END_REF] for more details). The woody vegetation (Bernadouze catchment area) is mainly composed of Fagus sylvatica L. (beech) with a minor component of Abies alba type (fir) localized mainly above 1500 m a.s.l. Three types of stands are currently dominant: (i) a high beech forest (derived from the conversion of former coppices), (ii) an old fastgrowing coppice, and (iii) a line of stunted coppice around the Bernadouze peat bog.

In the 46ha investigated area of the site, 80 former CK terraces were found, which were geolocated and integrated in a GIS database and sampled for analysis (Fig. 1). Combination of anthracological study of 28 of them with soil charcoal analysis (pedoanthracology) and historical data provided a reconstruction of the forest composition, structure, humans uses and system of management from the late Bronze Age to contemporary times [START_REF] Py-Saragaglia | Chapter 10: Long-term forest evolution and woodland uses in an ancient charcoal-production forest of the eastern French Pyrenees: an interdisciplinary approach at high spatio-temporal resolution[END_REF][START_REF] Saulnier | A study of late Holocene local vegetation dynamics and responses to land use changes in an ancient charcoal making woodland in the central Pyrenees (Ariège, France), using pedoanthracology[END_REF].

According to oral testimony, we identified one of the last used CK terraces (57°09'.51''N and 61°90'.45''E) dated from at least the late 1930s (Fig. 1). This CK terrace was located within the current high beech forest on a relatively flat area (approx. 5% slope), at the edge of the Bernadouze peat bog (called in the past "Freychinède"; in [START_REF] Jalut | A palaeoecological approach to the last 21 000 years in the pyrenees: The peat bog of Freychinede (alt. 1350 m, Ariege, South France)[END_REF], on a topographic terrace in direct proximity to a limestone outcrop almost completely covered with mosses. Around this CK we located 4 other CKs during archaeological field survey including only one with a radiocarbon date in the mid 17th-20th c. (Lyon-14957 -SacA-53113-, 140 ± 30 BP, 1669-1945 (95.4%) Cal AD with a high probability in the 1669-1891 period (79.9%)) (Fig. 1).

Sampled living beech trees were mainly located in the northwestern part of the woods along an altitudinal transect between 1345 and 1410 m a.s.l. (Fig. 1). The two sampling areas were located in the heart of the current high forest separated from each other by a poorly accessible escarpment area and a small forest corridor mainly composed of young silver fir trees (20-30 years old) associated with a few other broadleaf trees.

Sampling strategies

Archaeological survey and charcoal sampling

CK terrace was cleaned of its litter (leaves and plant debris) to define its extent. A surface of 16 m², including most of the CK surface, was divided into 16 squares of 1 x 1m. For this experimental study, a 4 m long and 1 m wide survey trench was dug from west to east, across the centre of the terrace, in the northern part of the grid. The CK fill was then gradually stripped layer by layer (Stratigraphic Unit, SU) until the ochre mineral layer was reached. In charcoal-rich layers, at least one 10-litre bucket of sediment (= one sample) was collected in each square (4 in total). In total, 10 samples (100 litre of sediment) were collected during the investigation (Fig. 2). Finally, the northern stratigraphic profile of the trench was described and drawn (Table 1, Fig. 3).

Micromorphological sampling

Directly on the stratigraphic profile, oriented sediment columns were sampled at 2 points: (i) A a continuous column here referred to as "Bern1" was collected at the west end of the survey trench (square 4) (Fig. 3 a,b);

(ii) A second column was taken in the east half of the profile, close to the middle of the CK fill (square 2), (Fig. 3 a,c). This column, referred to as "Bern 2-3", was divided into 2 segments because it was not possible to remove another continuous column from its soft sediments.

The sampling of both columns required digging into the ochre mineral layer on which the archaeological survey had stopped, to a depth of 10 cm.

Dendrochronological sampling

We took advantage of a forest logging managed by the National Forests Office (ONF) to carry out sampling of living beech. In the autumn 2016, 64 discs were extracted by lumberjacks directly from the stumps of freshly felled trees. Every sampled stump was geolocated (Fig. 1). Sampled stumps were selected according to their estimated age (i.e selection of oldest/largest trees), their type (high beech trees) and their location in the forest.

Laboratory Analysis

Anthracological analysis

The sampled sediments were wet-sieved with a 4 mm mesh sieve and then dried in the shade. Charcoal fragments were identified using a reflected light optical microscope (Leica DM4 M) with magnifications of 100 to 1000 according to the three anatomical planes of wood (transverse, tangential and radial). All charcoal fragments were identified at the level species, genus, family or sub-family [START_REF] Schweingruber | Anatomie microscopique du bois: identification de matériel récent et subfossile d'essences de l'Europe centrale en tenant compte de la variabilité de structure du bois de tronc et de branche[END_REF][START_REF] Vernet | Guide d'identification des charbons de bois préhistoriques et récents: sud-ouest de l'Europe : France, péninsule ibérique, et îles Canaries[END_REF]. The transverse sections were carefully cut (see below) and set aside for measuring ring width in a second step. In addition, when bark was preserved, its position in the last ring was recorded.

Micromorphological analysis

After drying and resin impregnation of the columns, a total of 5 large thin sections were carried out in PACEA laboratory (CNRS-Univ. Bordeaux). The thin sections were visually observed to distinguish their general characteristics and scanned to obtain digital images for meso-scale observation [START_REF] Haaland | Documenting archaeological thin sections in high-resolution: A comparison of methods and discussion of applications[END_REF].

At the microscopic scale, we used a Zeiss Axio Imager A2m microscope allowing observations in transmitted light (TL), Plane Polarised Light (PPL), Cross Polarised Light (XPL), and Oblique Incident Light (OIL) at magnification x25 -x500. The OIL observations were made using additional Zeiss gooseneck fibre optic illumination. The microphotographs were taken using a Zeiss Axiocam ICc5 digital camera controlled by Zeiss Axiovision software. The descriptive terminology follows [START_REF] Stoops | Guidelines for analysis and description of soil and regolith thin sections[END_REF] and the interpretative framework is based on [START_REF] Stoops | Interpretation of micromorphological features of soils and regoliths[END_REF].

Dendrochronological analysis

Wood and charcoal sample preparation

The transverse surface of the living tree discs was prepared by fine sanding to grit size of 600 using an electric sander to magnify tree ring features. Then discs were divided into several pieces to facilitate measurement.

While some authors have suggested that charcoal fragments with at least 50-60 rings are more suitable for conducting a dendrochronological analysis [START_REF] Kuniholm | How to Collect Archaeological Wood and Charcoal for Dendrochronological Analysis[END_REF][START_REF] Blondel | Impact of carbonization on growth rings: Dating by dendrochronology experiments on oak charcoals collected from archaeological sites[END_REF], other studies have already successfully constructed robust chronologies using charcoal fragments with 30 rings [START_REF] Backmeroff | Dendrochronological dating of charcoal kilns: a new method for dating historical land use at the upper timberline[END_REF][START_REF] Pichler | Dendrochronological dating of charcoal originating from an Early Iron Age fire-set pit in the mining area of Schwaz/Brixlegg (Tyrol, Austria)[END_REF][START_REF] Rouaud | Les forêts de pente de la haute vallée de la Dordogne : enjeux écologiques et énergétiques d'une ancienne forêt charbonnée[END_REF]. Consequently, we selected fragments preserving at least 30 rings to maximise sample size while also enabling obtaining robust measurements.

To facilitate charcoal-growth-ring measurement, the sampled pieces were simply broken by hand or with a cutting tool, i.e. a scalpel, to obtain good cross-sections. This method is particularly suitable for hard charcoal pieces [START_REF] Pichler | Charcoal from a prehistoric copper mine in the Austrian Alps: dendrochronological and dendrological data, demand for wood and forest utilisation[END_REF].

Tree-ring measurement

For both living tree and charcoal samples, we applied a visual cross-dating method to assign to each ring a year of formation [START_REF] Stokes | An introduction to tree-ring dating[END_REF]. Tree rings were then measured with an accuracy of 0.01 mm using the LINTAB™ measuring table combined with TSAPWIN™ Scientific software [START_REF] Rinn | TSAP-Win User Reference (Version 0.53[END_REF]. Two to three radii (elementary series) were measured per sample. This technique provides a more representative individual average of tree growth [START_REF] Bourquin-Mignot | Construction d'une longue chronologie de hêtres au Pays-basque. La forêt d'Iraty et le Petit Âge Glaciaire[END_REF].

Cross-dating of series (based on the identification of characteristic years such as narrow rings) was conducted in order to detect the presence of false rings or missing rings during measurement [START_REF] Latte | La dendroécologie: un outil pour affiner nos connaissances sur l'autoécologie des essences forestières[END_REF]) that could reduce the common signal of all the samples. Visual checks and crossdating were conducted using CDendro 7.8.1 to support the validation or rejection of potential adjustments to measurements. The accuracy of the cross-dating was evaluated by calculating the cross-correlation values between series. We considered a cross-correlation value of 0.5 as the threshold to decide to keep or exclude a series from the dataset.

The best correlated series were then averaged in order to produce a master chronology for (1) living tree and (2) charcoal. Finally, master chronologies from both the living tree and charcoal were detrended using a cubic fixed spline in order to remove the effects of age, competition and local specific conditions. This method is well-adapted to trees growing in closed-canopy forest and/or in frequently disturbed forests (competition, fire, etc.) [START_REF] Cook | Methods of Dendrochronology: Applications in the Environmental Sciences[END_REF][START_REF] Lambert | Dendrochronologie et dendroclimatologie du chêne en France, Questions posées par le transfert de données de bois historiques vers la dendroclimatologie[END_REF]. The detrended chronology from charcoal fragments was then cross-dated with the detrended chronology of the Forêt de Bernadouze (felled trees) in order to propose a date for the operation of the CK.

In order to increase the number of charcoal fragments that kept the bark, which are mandatory to effectively assess charcoal making practices (seasonality of felling), we subsequently calculated new cross-correlation values between the master chronology of charcoals and part of the previous, excluded charcoal fragments. For this additional analysis, we set the new acceptable correlation value at 0.4. The new selected charcoal dataset was only used to study charcoal making practices and these additional fragments were not added to the master chronology.

Results

Archeological evidence and charcoal features

The black circular area corresponding to the location of the last CK(s) built on the terrace measured 5 m in diameter (Fig. 2). The area formed a slight basin-shaped depression which was filled with organic matter and charcoal-rich black earth up to a maximum depth of 35 cm within the trench survey. We distinguished six stratigraphic units (SUs) (Table 1; Fig. 3).

In total, we selected 303 charcoal fragments with 104 fragments extracted from the SU 2 and 199 from the SU 3. All these charcoal fragments were identified as beech (Fagus sylvatica L.). Among the 94 fragments that had preserved their bark, the state of progress in the terminal ring formation varied with 42 fragments exhibiting a terminal ring with final wood, 9 fragments exhibiting a terminal ring with only the early wood, and 3 fragments exhibiting a terminal ring at the end of the early wood. The minimum diameter measurement with the calliper shows diameters raging between 0.8 and 3.6 cm after correction (systematic addition of +25% to compensate for shrinkage losses during charring (Paradis-Grenouillet, 2012).

Microstratigraphic features

Sample Bern 1 was processed into 3 thin sections covering about 35 cm in sequence. The succession of 7 facies (Bern1_a-g) (Fig. 4) were distinguished. Sample Bern2-3 allowed the realization of two thin sections, covering a sequence of about 28 cm, in which 4 facies (Bern2-3_a-d) were distinguished (Fig. 5).

The description and the observation of facies both macroscopically (color, coarse constituents, etc.) and microscopically (fine constituents, pedofeatures, microstructure, etc.) highlighted important similarities between Bern 1 and Bern 2-3. We were able to synthesize these observations as 5 micro-facies types (MFT) (Fig. 6).

The detailed micromorphological description of Bern1 and Bern2-3 is presented in Table 1 of the Supplementary material (Table S1)

MFT_A (e,g. This MFT corresponds to a forest soil deep mineral horizon formed by the alteration of underlying dolomitic marlstone formation and, to a lesser extent, input from the erosion of the surroundings metamorphic and crystalline formations [START_REF] Trenet | Alus-les-Bains. Notice[END_REF]. These sediments were excavated to create the depression in which subsequent CK burning episodes took place.

MFT_B (e.g. Bern1_a,Bern1_c, This MFT_B is formed by rounded aggregates of material similar to the local sediment that may correspond to colluvial input into the depression. These materials are regularly affected by biopedological agents (e.g. soil mesofauna, roots) [START_REF] Kooistra | Features Related to Faunal Activity[END_REF][START_REF] Bertran | Géoarchéologie et taphonomie des vestiges archéologiques : impacts des processus naturels sur les assemblages et méthodes d'analyse[END_REF], leading to blending with combustion residues from over-and underlying deposits.

MFT_C ( e.g. Bern1_b,Bern1_d,Bern1_f, Most of the CK terrace depression fill is related to MFT_C, which was formed by a mix of mineral and organic residues related to charcoal making, including charcoal and redeposited oxidized soil aggregates. Soil mesofauna and root action disturbance were important and sometimes complemented by climatic processes such as freeze/thaw cycles (Vliet-Lanoë, 2010), all of which led to mixing and homogenization of what were likely originally differentiated deposits.

MFT_D (e.g. Bern1_e, Bern2-3_c)

This MFT is mainly characterised by in-situ thermally altered deposits, corresponding to original kiln surfaces. These deposits may correspond to mainly mineral material such as those of MFT_B (Bern2-3_c) or those formed of previous combustion residues such as MFT_C (Bern1_e).

MFT_E (e.g. Bern1_g)

The summit of the terrace depression fills a humic horizon, with mostly organic constituents and virtually no mineral content. Organic constituents are regularly altered and fragmented.

Felled trees master chronology

The mean age of felled trees was 145 years, with a maximum of 162 years and a minimum of 89 (Supplementary material, Table S1). The ring widths varied from 0 (missing rings) to 31 mm, with a mean of 10 mm. We used the 24 best cross-dated individual chronologies to build the master chronology of the Forêt de Bernadouze, resulting in a master chronology 137 years long, covering the 1879-2016 period. The average correlation between all individual series was 0.6 and the average Student's t-test value is 11.9. The square correlation matrix highlights the very good synchronization of the individuals (Fig. 8a).

Our chronology was compared with two other beech chronologies from the Pyrenees: the Iraty forest chronology (Bourquin-Mignot and Girardclos, 2001) and the Baish Aran [START_REF] Cavin | Highest drought sensitivity and lowest resistance to growth suppression are found in the range core of the tree Fagus sylvatica L. not the equatorial range edge[END_REF]. Correlation with the Iraty forest chronology (located about 200 km as the crow flies on the same slope of the French Pyrenees), was very strong with a Student t-test around 6.85, and a mean correlation of 0.48. It was weaker with the Baish Aran chronology located on a southern slope in the Spanish Pyrenees (53 km as the crow flies) with values of 3.6 for the Student's t-test and of 0.33 for the mean correlation (Fig. 7).

Charcoal kiln master chronology

In total, 104 charcoal fragments from SU 2 and 199 from SU 3 were measured. The mean length of individual charcoal chronologies was 31, with a maximum of 66 years and a minimum of 6 (Supplementary material, Table S2). The charcoal ring widths varied between 0 (missing rings) and 25 mm, with a mean of about 4 mm.

Only 167 fragments presented at least 30 rings i.e 55.1% of total. The analysis and comparison of these series allowed the selection of 49 individuals (20 from SU 2 and 29 from SU 3) exhibiting very high cross-correlations for constructing a master chronology of 65 years. The average correlation between individual series was 0.6 and the average Student's t-test was 7.1. The square correlation matrix showed the good quality of cross-dating of these series (Fig. 8b). Using the master chronology (felled trees) of the "Forêt de Bernadouze", we were able to crossdate the CK master chronology for the 1878-1942 period, but we preferred to keep only the 1881-1941 period, for which the accuracy of the cross-dating is higher, with a Student t-test value of 10.4 and a mean correlation of 0.7 (Fig. 9).

Using the charcoal master chronology, we were able to cross-date 33 additional charcoal fragments that were first excluded. These 33 new samples enabled us to obtain more charcoal fragments with bark. On the 82 total accepted fragments, 16 had retained the bark. Based on this dataset, 8 different felling years between 1924 and 1942 were obtained with certainty (Fig. 10a,b).

Discussion

Characterisation of the CK record

The combination of archaeological and micromorphological observations allowed us to better characterise (i) the CK fill, i.e. the archaeological record used to reconstruct past charcoal production and related woodland harvesting, (ii) the relative chronology of charcoal deposits. The two approaches also highlight some of the taphonomic processes that have affected the CK infilling and evidence the multiperiod nature of the CK use.

The CK terrace was established on the local regolith observed in MFT_A (SU 6). This MFT corresponds to a forest soil deep mineral horizon formed by the alteration of the underlying dolomitic marlstone and, to a lesser extent, by input from the erosion of the surroundings metamorphic and crystalline formations [START_REF] Trenet | Alus-les-Bains. Notice[END_REF]. Subsequent deposits (MFT_ B to E, SU 5, 4, 3, 2) recorded a succession of phases related to the use, maintenance and abandonment of the CK.

Phases of actual use of the CK are best exemplified by in-situ thermally altered sedimentary deposits (MFT_D) and reworked charcoal fragments, corresponding to surfaces on which the kiln was set and operated. Similar, but reworked, oxidized materials are also observed in other MFTs and suggest that in several instances these surfaces were disturbed. Indeed most of the CK depression fill is formed by a mix of redeposited mineral and organic residues related to charcoal production (MFT_C, SU 3 and 2). This result was expected considering that most of the produced charcoal was usually recovered and the terrace cleaned after the long process of charring [START_REF] Battesti | Une charbonnière expérimentale, Méthodologie de l'anthracologie sur charbonnière pour une approche de la gestion du milieu forestier en Languedoc[END_REF][START_REF] Izard | L'art du charbonnier: contributions ethno-botanique et géographico-historique à l'étude des paysages métallurgiques d'après l'anthracoanalyse des charbonnières[END_REF][START_REF] Lugli | Modelli e finalità nello studio della produzione di carbone di legna in archeologia[END_REF][START_REF] Montanari | The study of charcoal-burning sites in the Apennine mountains of Liguria (NW Italy) as a tool for forest history[END_REF][START_REF] Raab | Architecture of relict charcoal hearths in northwestern Connecticut, USA[END_REF]. Therefore, it is difficult to accurately evaluate the number of clearance and reuse phases, partly because of their own destructive nature, but also because their remains include decaying organic material in a biopedologicaly active environment and subject to bioturbation by soil mesofauna and root action [START_REF] Kooistra | Features Related to Faunal Activity[END_REF][START_REF] Bertran | Géoarchéologie et taphonomie des vestiges archéologiques : impacts des processus naturels sur les assemblages et méthodes d'analyse[END_REF] sometimes complemented by climatic processes such as freeze/thaw cycles (Vliet-Lanoë, 2010). Thus, this combination of anthropogenic (practices) and taphonomic factors led to a mixing and homogenization of originally differentiated deposits and the macroscopic appearance of the fill as a homogeneous black layer.

Based on micromorphological results, the depression infill displays obvious microstratigraphic succession, notably because of the intercalation of MFT_B with MFT C. On the margins of the CK, a succession of layers can be observed that indicate this sequence, whereas in the center the deposits are much more homogeneous. MFT_B is formed by rounded aggregates of material similar to the local sediment that may correspond to colluvial input into the depression related to wall collapse, potentially favoured by weathering [START_REF] Langohr | Creusement, érosion et comblement des fossés ; l'approche des sciences de la terre[END_REF]. This assumption is verified both at the edge of the structure, with material additions characteristic of CK edge collapse backfill, and in the centre of the CK where a similar deposit is present. Those materials are also regularly affected by biopedological agents leading to blending with combustion residues from over-and underlying deposits. This pattern implies, therefore, a potential mixing of charcoal fragments from different charcoal production episodes in SUs 2 and 3. Lastly, the summit of the depression fill corresponds to a humic horizon (MFT_E, SU 2), which represents the development of forest soil [START_REF] Gerasimova | Topsoils -Mollic, Takyric and Yermic Horizons[END_REF] after the last episode of the CK use.

Overall, a succession of different episodes of charcoal production is detected through microstratigraphy where it was otherwise undetectable by macroscopic observation. Some are observable through the presence of in-situ thermally altered sediment (MFT_D), but others are only evidenced through their reworked residues (MFT_C). The infilling of the depression by unaltered sediments points to abandonment phases indicating discontinuous charcoal production activities. This kind of multiperiod use was expected for this type of structure and activity [START_REF] Dupin | The taphonomic characterization of a charcoal production platform. Contribution of an innovative pair of methods: Raman analysis and micromorphology[END_REF], but the assessment of the different types of activities conducted and their precise chronology remains a challenging endeavor because of their destructive nature and their taphonomic context [START_REF] Raab | Architecture of relict charcoal hearths in northwestern Connecticut, USA[END_REF].

The anthracological remains in such a structure are very abundant, but their state of preservation and location are clearly affected by charcoal production practices (e.g. raking out of charcoal, cleaning, reuse of burnt charcoalrich earth) and also bioturbation and weathering processes (e.g. soil mesofauna, root action, frost cycles).

Consequently, these charcoal remains are often subject to further fragmentation [START_REF] Elmer | Active removal of biochar by earthworms (Lumbricus terrestris)[END_REF] and subsequent vertical and lateral displacement. As a result, charcoal fragments sampled from SUs 2 and 3 likely represent a palimpsest of several different charcoal production phases carried out on the terrace (MFT_C).

Accuracy of the CK master chronology

The statistical report for the chronology of Bernadouze's CK has a Student's t-test value of 10.4 and a mean correlation of 0.7. Comparison with similar studies (Supplementary material, Table S3) shows that results of our cross-datation are in the range of those obtained in previous studies (Student's t-test ranging between 3.7 and 12.2 and mean correlation between 0.6 and 0.7), which highlights the high robustness of the chronology of Bernadouze's CK. Indeed, in France, the experimental study conducted by [START_REF] Blondel | Impact of carbonization on growth rings: Dating by dendrochronology experiments on oak charcoals collected from archaeological sites[END_REF] on the site of Saint-Paulien (Haute-Loire) gave a 67-year chronology built using 4 charcoal fragments and validated with 7 regional chronologies with Student's t-test values ranging between 3.7 and 5.5. In the Schwaz/Brixlegg region (Tyrol, Austria), charcoal fragments from the mine of Mauk E allowed the construction of a 149-year master chronology, based on 133 charcoal samples, which was cross-dated with two regional chronologies, giving Student's t-test values of 4.5 and 5 [START_REF] Pichler | Charcoal from a prehistoric copper mine in the Austrian Alps: dendrochronological and dendrological data, demand for wood and forest utilisation[END_REF]. In Patagonia (Argentina), 3 charcoal fragments from the archaeological site of Cerro Pintado give a 50-year chronology with a mean correlation of 0.6 with a regional chronology [START_REF] Fermé | Dendrochronological study of charred wood at the Cerro Pintado archaeological site (Patagonia, Argentina)[END_REF]. In the USA, the 72 charcoal samples from the CKs of Old Greenwood and Gray's Crossing (Truckee, California) allowed the reconstruction of a 122-year master chronology that was cross-dated with one regional chronology giving a Student's t-test value of 9.3 and mean correlation of 0.6. In Mt Hope (Central Nevada), a floating chronology of 100 years based on 12 samples was cross-dated with a local chronology giving a Student's t-test value of 10.6 and a mean correlation of 0.7. Finally, in Central Nevada, the 202-year floating chronology constructed based on 41 charcoal samples from the Cortez site was cross-dated with a Student's t-test value of 12.2 and a mean correlation of 0.6 [START_REF] Strachan | Application of Dendrochronology to Historical Charcoal-Production Sites in the Great Basin, United States[END_REF].

We are aware that the high Student's t-test and mean correlation values we obtained may be influenced by the short length of both the individual charcoal chronologies, for which we got a maximum length of 66 years, and the master chronology, 61 years. However, the master chronology length is comparable to the length of the chronologies of Saint Paulien and Cerro Pintado, which were constructed with less than 5 samples. In addition, the CK charcoal master chronology has been cross-dated (contrary to some others except those of the USA) with a local living tree chronology built using beech trees growing on the same site. The high correlation values resulting from the cross-dating (charcoal vs living tree chronologies) allow us to confirm the hypothesis according to which, in the context of charcoal production, burnt wood is almost always collected close to the site of charcoal production. On the other hand, some samples showed strong correlation with very similar growth patterns suggesting an ecologically identical supply location (Blanc-Garidel and Wicha, 2019) or a common tree origin [START_REF] Pichler | Charcoal from a prehistoric copper mine in the Austrian Alps: dendrochronological and dendrological data, demand for wood and forest utilisation[END_REF]. Moreover, we used only charcoal which retained at least 30 rings to build the master chronology, although other studies suggested using charcoal fragments with at least 50 rings [START_REF] Kuniholm | How to Collect Archaeological Wood and Charcoal for Dendrochronological Analysis[END_REF][START_REF] Blondel | Impact of carbonization on growth rings: Dating by dendrochronology experiments on oak charcoals collected from archaeological sites[END_REF]. The silvicultural treatment implemented in the "Forêt de Bernadouze" was the "coppice selection system" (taillis fureté in French), which was considered as the best suited for European beech, especially in the Pyrenees [START_REF] Chassinat | Forêt domaniale de Suc-et-Sentenac, aménagement[END_REF][START_REF] Rochebrunel | [END_REF][START_REF] Salvador | Simples notes sur l'aménagement et l'exploitation des forêts pyrénéennes françaises[END_REF][START_REF] Nicolescu | Silvicultural Guidelines for European Coppice Forests[END_REF]. This system consists of cutting on each stump only shoots that have reached a target diameter, with 50 years old considered as the fertile age and 130 years the end of their exploitability. In 1892, the shoots were between 20 to 50 years old. In the same year, forest managers decided to convert this coppice into high forest [START_REF] Chassinat | Forêt domaniale de Suc-et-Sentenac, aménagement[END_REF]. They planned a 16-years rotation including biannual cleaning cuts from 1906 to 1921. In 1920, forest managers reported the existence of shoots between 20 to 70 years old and renewed the 16-years rotation from 1922 to 1937 [START_REF] Rochebrunel | [END_REF] and once again from 1938 to 1953 [START_REF] Noiriel | Projet de révision d'aménagement[END_REF]. The conversion into high forest was partly reached in 1965 [START_REF] Noiriel | Projet de révision d'aménagement[END_REF][START_REF] Izard | Historique et dynamique du paysage forestier de la vallée de Suc et Sentenac (Vicdessos, Ariège)[END_REF].

High-temporal chronological reconstruction of charcoal production activity on the target CK terrace

A non-intensive initial phase of charcoal making

The deeper MFT_C was observed in Bern 1_b and Bern 2-3_b layers and confirmed an early charcoal making phase on the terrace. The substratum was excavated by charcoal burners to make a flat area, slightly curved, and accommodate a wood pile of about 5 meters in diameter. This size corresponds with the average diameter of the CKs detected in the Forêt de Bernadouze with 5.1 x 4.7 meters in diameter, and elsewhere in the mountains of the northern Pyrenees [START_REF] Bonhôte | Charcoal kilns and environmental history in the eastern Pyrenees (France). A methodological approach[END_REF]. According to Bonhôte 1998, Annex 2, pp 220-223), in mountain areas this diameter could correspond to CKs consuming up to 15 to 20 "steres" (cubic metres) of wood. The construction technique of the type of kiln used in the Pyrenees is less well known than those found in other regions such as Liguria or Portugal [START_REF] Poggi | Carbonai dell'Antola, video-cassetta 40[END_REF]. In the studied CK, the remains of this former episode were thoroughly scraped, thus there is very little evidence with which to reconstruct the process. However, it is not possible to conclude if this phase involved a single episode of the use of the terrace (only one CK) or different successive uses. In the Mediterranean region, several charcoal making episodes (up to five successive events) can take place on the same CK in a single year [START_REF] Battesti | Une charbonnière expérimentale, Méthodologie de l'anthracologie sur charbonnière pour une approche de la gestion du milieu forestier en Languedoc[END_REF]. It is possible, therefore, that this first phase of exploitation probably hosted several episodes in a short period of time. However, there are no historical records from the Pyrenean area to support this hypothesis.

No charcoal was collected and studied from the Bern1_b and Bern2-3_b layers corresponding to SU 5 and the deepest part of MFT_C. Therefore it is not possible to date precisely this first charcoal production event. Based on our knowledge of the charcoal production chronology in the "Forêt de Bernadouze", we hypothesise that this former phase is potentially attested to elsewhere in this area. In addition, it may be only a few years or decades ahead of the next. Indeed, two CKs were radiocarbon dated at the same time i.e. between the end of 19th and the beginning of the 20th c. with a high probability of dating from the 1811-1926 period (71.5% of probability) (Poz-73024, 65 ± 30 BP, 1691-1921(95.4%) Cal AD and Lyon-14958 (SacA-53114), 80± 30 BP, 1690-1926 (95.4%) Cal AD)). The anthracological study of these two CKs also showed a taxonomic composition similar to the target CK discussed (Sus 1 to 3). While one CK was composed almost exclusively of beech (97.4%) with some trace of fir (1.5%), the other CK composition indicated the combined use of beech (75.5%) and fir (23.6%). It could therefore be envisaged that this first charcoal making event can be dated to the early 20th c. Based on this assumption, this initial phase of charcoal production recorded in the CK could be related to the 16-years rotation set between 1906 to 1921 [START_REF] Chassinat | Forêt domaniale de Suc-et-Sentenac, aménagement[END_REF][START_REF] Rochebrunel | [END_REF].

The MFT_B, observed both at the edge (Bern1_c) and in the middle of the terrace (Bern2-3_c) and corresponding to SUs 3 and 4, contained inputs of eroded sediment from the edge of the CK terrace. This partial filling of the terrace, attributed to the collapse of the edges, suggests an abandonment phase. However, the duration of this phase of abandonment, which probably ranged from one to several years remains difficult to accurately date. In mountain areas, the duration of the infilling process depends mainly on the instability of the environment (i.e erosion, landslides, often accentuated by acute slopes) and on the combined effects of climatic and edaphic factors [START_REF] Izard | Les montagnes du fer : éco-histoire de la métallurgie et des forêts dans les Pyrénées méditerranéennes (de l'Antiquité à nos jours) : pour une histoire de l'environnement[END_REF][START_REF] Langohr | Creusement, érosion et comblement des fossés ; l'approche des sciences de la terre[END_REF]. [START_REF] Dupin | The taphonomic characterization of a charcoal production platform. Contribution of an innovative pair of methods: Raman analysis and micromorphology[END_REF] showed that the reuse of a CK, without a long period of interruption, generally results in a homogeneous stratigraphy where abandonment markers are completely removed by cleaning. This evidence suggests that a more or less prolonged interruption in the use of a CK may result in the creation of an abandonment layer interspersed between the charcoal sequences. This hypothesis presents the best fit with the evidence collected from the present study. Comparable observations on other CKs have been interpreted as an indication of abandonment and multiperiod use of the terrace [START_REF] Dubois | Archéologie de la forêt charbonnée : questions et méthodes, illustrées par I'exemple du site de Lercoul (Ariège)[END_REF][START_REF] Izard | Les montagnes du fer : éco-histoire de la métallurgie et des forêts dans les Pyrénées méditerranéennes (de l'Antiquité à nos jours) : pour une histoire de l'environnement[END_REF][START_REF] Montanari | The study of charcoal-burning sites in the Apennine mountains of Liguria (NW Italy) as a tool for forest history[END_REF][START_REF] Nelle | Woodland history of the last 500 years revealed by anthracological studies of charcoal kiln sites in the Bavarian Forest, Germany[END_REF][START_REF] Gebhardt | Impact of charcoal production activities on soil profiles: the micromorphological point of view[END_REF][START_REF] Rouaud | Les forêts de pente de la haute vallée de la Dordogne : enjeux écologiques et énergétiques d'une ancienne forêt charbonnée[END_REF][START_REF] Raab | Architecture of relict charcoal hearths in northwestern Connecticut, USA[END_REF].

Repeated combustion on one site can lead to relative sterility of the soil, which can take several decades to recover (Acovitsioti-Hameau, 2001b;[START_REF] Carrari | Former charcoal platforms in Mediterranean forest areas: a hostile microhabitat for the recolonization by woody species[END_REF] and lead to its abandonment. The abandonment also may be directly related to the management of the surrounding forest resource. The CK remains represent material evidence for the management of a reduced surface area of the forest (called "parcel management" by [START_REF] Fabre | Première approche du charbonnage et de l'histoire de la forêt méditerranéenne: l'exemple du canton d'Aniane (Hérault)[END_REF]. Availability of wood fuel is not infinite, and this forces the temporary abandonment of a "parcel" for the time coppice/woodland takes to regenerate [START_REF] Battesti | Une charbonnière expérimentale, Méthodologie de l'anthracologie sur charbonnière pour une approche de la gestion du milieu forestier en Languedoc[END_REF]. This refers to the notion of the logging district, based on a rotation of cutting cycles [START_REF] Acovitsioti-Hameau | L'habitat des artisans de la forêt en Moyenne Provence: l'exemple des charbonniers[END_REF], in theory 8, 10, 12 or 15 years, which is a forest management method very often used in montane areas [START_REF] Duhamel Du Monceau | Art du charbonnier, ou Manière de faire le charbon de bois. Description des arts et métiers, faits ou approuvés par Messieurs de l'Académie royale des Sciences[END_REF][START_REF] Salvador | Simples notes sur l'aménagement et l'exploitation des forêts pyrénéennes françaises[END_REF] and also very well suited with charcoal production in order to ensure the sustainability of fuel resources.

The partial preservation of the filling layer also suggests that when the terrace was reoccupied for a new episode of charcoal making, (i) the charcoal burner did not consider it necessary to clean-up and dig-up the flat area, or (ii) the layer of sediment filling the pit was so thick that it was not completely removed. Indeed, when preparing the terrace for a new burning, the charcoal maker would start by thoroughly cleaning the terrace. Due to its quality (fine, smooth and not very porous), this mixture of earth and charcoals was collected in order to cover the new wood pile, allowing a better insulation from the air [START_REF] Battesti | Une charbonnière expérimentale, Méthodologie de l'anthracologie sur charbonnière pour une approche de la gestion du milieu forestier en Languedoc[END_REF][START_REF] Musset | Charbonniers, le métier du diable ? Le Monde alpin et rhodanien[END_REF]Acovitsioti-Hameau, 2001a;[START_REF] Hanus | Charbonnier, fais fumer! Chronique ethnographique d'une fête du Vercors[END_REF][START_REF] Rouaud | Les forêts de pente de la haute vallée de la Dordogne: retour sur le passée d'une forêt charbonnée[END_REF]. The reuse of abandoned terraces is often advised by modern technical treatises as well as by former and current charcoal burners because of specific soil properties (see the overall review in Krebs et al. 2017, Table 2). Indeed, charcoal production on an old terrace allows better control of pyrolysis and thus, increases charcoal productivity by up to 10% [START_REF] Larcher | Les charbonniers à Paimpont : contribution à l'histoire d'une commune[END_REF].

The second phase of use of the CK terrace

Overlapping SUs 3 and 2, the MFT_C (observed in Bern 2-3_d and from Bern1_d to Bern1_f) suggested a new phase of charcoal production. In the center of the terrace (Bern2-3_d), the micromorphological analysis of oxidized sediment highlighted several CK uses, which appeared more regular and intense than in the previous phase as suggested by the high concentration of charcoal fragments [START_REF] Dupin | Anthracology of charcoal kilns in the forest of Chailluz (France) as a tool to understand Franche-Comte forestry from the mid-15th to the early 20th century AD[END_REF]2019). The successive episodes of use led to significant mixing of the materials remaining inside the CK. Individual chronologies from two fragments collected in SU 3 (squares 1 and 2) suggested a restart date of activity in 1924 (Fig. 10 a,b). According to historical documentation, this date could correspond to the second phase of rotation between 1922 and 1937 [START_REF] Rochebrunel | [END_REF]. Following the activity restart, several charcoal making episodes were identified based on individual chronologies of charcoal remains in the early 1930s. Two other individual chronologies constructed from a charcoal fragment extracted from SU 3 in square 1 and from another fragment extracted from SU 3 in square 2 provided two additional probable cutting dates in 1931 and 1932 respectively. Although no dates were obtained between 1924 and 1931, the absence of markers indicating the abandonment of the terrace suggests that the charcoal production was probably continuous and intensive [START_REF] Dupin | The taphonomic characterization of a charcoal production platform. Contribution of an innovative pair of methods: Raman analysis and micromorphology[END_REF]. It is also highly likely that charcoal fragments from this period were completely removed during the systematic preparation of the CK to allow renewed use by the charcoal burner. In any case, three new dates provided by the individual chronologies of three other charcoal specimens extracted from SU 3 in squares 2 and 4 point to date a new event in 1935. This event may be linked to the intensification of charcoal production as revealed by micromorphology in this facies. Subsequently, three separate charcoal chronologies from SU 2 (2 charcoals) and SU3 (1 charcoal) (squares 2 and 3) provided three additional dating points in 1939. From 1939 onwards, the new dates obtained became annual: 2 charcoal fragments extracted from SU 2 and 1 fragment from SU 3 dated to 1940 (squares 2 and 4), 2 fragments from SU 3 dated to 1941 (squares 1 and 4) and finally 1 fragment collected in SU 3 (square 2) dated to 1942. These three consecutive dates characterise several successive, very brief and closely spaced, charcoal making episodes located at the transition from SU 3 to 2 (Fig. 10 a,b). In addition, the presence of reworked, heterogeneous materials, mixing weathered sediments and charcoal fragments, at the top of the CK terrace fill testifies to the occurrence of one or more charcoal making operations. The lumpy microstructure, abundant earthworm droppings and the traces of root action in this part of the facies, are characteristic of the homogenization of this deposit. Once filled in, the terrace appears to have been permanently abandoned since pedogenesis resumed in the form of the establishment of a humus horizon (MFT_E, Bern1_g). The CK was definitively abandoned in 1942 (Fig. 10 a,b). Finally, the chronology demonstrates that the CK stratigraphy was totally mixed since both the most recent and the oldest charcoal making episodes were found in SU3, while the intermediate episodes were evidenced in SU2. This result is consistent with the CK complex history reconstructed through micromorphological analysis.

Some insights about charcoal burner practices and seasonality

The formation progress of the terminal ring provided new information about the seasonality of harvesting practices for supplying wood. Our results suggested that trees felling usually occurred at the end of the growing season, i.e. in autumn. Previous archeological works already highlighted that fuelwood harvesting for charcoal making was generally done from late autumn until the first sap rises at the end of winter [START_REF] Montanari | The study of charcoal-burning sites in the Apennine mountains of Liguria (NW Italy) as a tool for forest history[END_REF]Fouédjeu et al., 2019). The period from late winter up to early spring (before the beginning of growing season) is also considered to be the best time to cut coppices for new shoot growth [START_REF] Nicolescu | Silvicultural Guidelines for European Coppice Forests[END_REF]. As the felled wood was rarely charred green, it was first stacked in the forest for a few months [START_REF] Bond | Medieval Charcoal-burning in England[END_REF], providing a partially dry wood that was easier to burn. Thereafter, wood was charred, very often as early as the following spring, i.e. in the same year of felling or at the latest at the beginning of the following year [START_REF] Duhamel Du Monceau | Art du charbonnier, ou Manière de faire le charbon de bois. Description des arts et métiers, faits ou approuvés par Messieurs de l'Académie royale des Sciences[END_REF][START_REF] Py-Saragaglia | Chapter 10: Long-term forest evolution and woodland uses in an ancient charcoal-production forest of the eastern French Pyrenees: an interdisciplinary approach at high spatio-temporal resolution[END_REF]. According to local modern and contemporary historical records charcoal burning season in the mountains of Ariège lasted generally from September, sometimes late August, to June or some cases to July, with a short stop during the months of January and February (Bonhôte, 1998, annex 2, pp 220-223). In addition, the log calibre reconstruction of charcoal fragments highlights that charcoal burners used mainly very small-calibre wood (0.8 to 3.6 cm), i.e. shoots with their terminal twigs, confirming that wood fuel came from coppice harvesting [START_REF] Py-Saragaglia | Chapter 10: Long-term forest evolution and woodland uses in an ancient charcoal-production forest of the eastern French Pyrenees: an interdisciplinary approach at high spatio-temporal resolution[END_REF].

Various charcoal fragments lacked terminal rings; these 66 charcoals designate terminus post quem dates. However, even though we were unable to date accurately their felling dates, these fragments may highlight one or more additional phases/episodes of charcoal making and/or the continuity of practices on the CK terrace. Indeed, we cannot completely exclude the possibility that the last dated rings correspond to the felling of the shoot. In fact, it is not unusual in an archaeological context for excavation, sampling or sieving to increase the risk of charcoal fragmentation and bark detachment (bark cells are already poorly resistant to carbonisation) [START_REF] Blondel | Impact of carbonization on growth rings: Dating by dendrochronology experiments on oak charcoals collected from archaeological sites[END_REF]. But in the absence of sapwood which is not visible on beech wood, it is impossible to validate that the last dated ring corresponds to the date of felling.

Conclusion

This paper (i) highlighted, through the implementation of an original and integrated approach, the complexity of the CK tool, archive and chronology of charcoal making practices, and (ii) provided new insights about the successive phases of use and abandonment of an ancient CK terrace. By combining micromorphology and dendrochronology, we have also been able to characterise the rate of terrace reuse at high temporal resolution. Micromorphological analysis highlighted the multiperiod nature of the CK use. The dendrochronological study of charcoal fragments allows us to overcome the radiocarbon dating (AMS) limitation by improving the accuracy of CK use dates to an annual resolution. Our study reveals (i) a regular and repetitive use of the CK terrace and (ii) its operation over a very specific period of time. We therefore highlighted two major phases of charcoal production interspersed with a phase of abandonment resulting from the need for (to allow) the beech coppices to grow and to regenerate. The first phase has only been evidenced by the microstratigraphy analysis but cannot be dated due to the lack of charcoal fragments. The second phase was precisely dated between 1924-1942 and was characterised by an intensification of charcoal making production of this CK from 1939 onwards becoming annual.

Our study highlighted the relevance of an integrated approach to better characterise the long-term variation in intensity of charcoal production and of its pressure on local woodland. All results underline the need to develop such integrated approaches for the study of post-1650 CKs in order to refine the chronology of charcoal burner practices and better understand their impact on forest ecosystems. Locally, the construction of the "Forêt de Bernadouze" master chronology has provided a new dendrochronological reference for beech on the northern slope of the Pyrenees and highlights the need to build new reference sets, at least for the dominant species (beech, fir, oak, pine). As emphasized by [START_REF] Deforce | Tree-ring analysis of archaeological charcoal as a tool to identify past woodland management: The case from a 14th century site from Oudenaarde (Belgium)[END_REF], the comparison of charcoal growth patterns with those derived from current forest management models (coppice, coppice with standard, selection coppice system, high forest, etc.) may facilitate the reconstruction of past woodland management practices.

Figures

Figure 1: The Forêt de Bernadouze represented with sampling points of charcoals kilns (CK). Green squares represent four CKs near to CK target and red square CKs dated by radiocarbon between 19th and 20th c. White points represent the totality of 64 dendrochronological samples. White circles, the 46 dendrochronological samples analysed for this study. The red star represents the CK selected for the micromorphological and dendrochronological study.

Supplementary material

Table S1: Detailed micromorphological description of the five soil thin sections Table S2: Description of individual and master chronologies built on felled trees and charcoal fragments Table S3: Synthesis of the different charcoal-based chronologies used to compare the new Bernadouze charcoal chronology This SU is located directly on the mineral layer. It contains small, powdery rocks in some sections and some rare glass fragments.

4

5-12

This SU presents a mix of sediment and charcoals (low density) inserted in the SU 3 at the western edge of the section. Charcoals in this SU have not been studied (Fig. 3D). The limit is clear, but irregular with Bern1_c. The microstructure and constituents are similar to those of Bern1_b. The presence of charcoals and micro-charcoals is significant and aggregates of oxidized sediments are occasionally observed. The biological features frequently intersect these usual elements of superficial horizons affected by successive freeze/thaw cycles.

Bern 1 Bern1_c

The contact is clear with underlying facies and mineral materials are similar to those of Bern1_a. Frequent charcoals and microcharcoals (less than in Bern1_b however) are present in the porosity. Moreover, aggregates of materials similar to those Bern1_b as well as centimetric and sub-rounded aggregates of yellow silty-clay are observed.

and 4

Bern 1 Bern1_b

The contact with Bern1_a is diffuse. The materials are similar to those of Bern1_a, with, however, a markedly higher proportion of organic constituents, notably millimeter-to centimeter-sized charcoal. Aggregates of oxidised sediments are also observed. The biological activity is materialized by numerous fillings formed by faecal pellets and root tissues.

5

Bern 1 Bern1_a Reworked material with a crumby/granular microstructure and rounded aggregates. Mixed with combustion residues (burned sediments, charcoal) through biological activity (mesofauna channels, roots).

Figure 2 :

 2 Figure 2: Protocol illustrating excavation and sampling of the CK. (a) White dotted circle shows the CK boundary after removal of litter. Grid with 1 x 1 m squares. Sampling of a bucket of 10 litres of sediment per tiles and SU. (b) Excavation of a quarter of the CK.

Figure 3 :

 3 Figure 3: (a): North section of the trench. (b) and (c). Detail of the western and the central part of the section after an additional 10 cm, discovery of a fifth SU and the localization of micromorphological samples. (d): Schematic survey of the trench.

Figure 4 :

 4 Figure 4: Scans and microphotography of Bern1 sample. (A) Scan of the indurated sample with schematic organization of micro-facies (a to g). (B, C and D) Scans of the three thin-sections. (E) Bern1_a: main coarse mineral materials; (red arrow: phlogopite; yellow arrows: dolomite) surrounded by carbonated clay-silt. (F) Bern1_b: centimetric charcoal surrounded by a mix of mineral and organic materials. (G) Bern1_b: heterogeneous constituents and organization with both unaltered aggregates (yellow arrow), oxidized aggregates (red arrow) and carbonaceous material (blue arrow). (H) Bern1_c: integration of combustion residues-rich material as mesofauna faecal pellets (yellow arrows) and a large rounded aggregate of unaltered sediment (red arrow). (I) Bern1_d: aggregate of oxidized sediment in secondary position within a heterogeneous matrix including abundant combustion residues and biological features. (J) Bern1_g: Topsoil organic layer with abundant leaf fragments (red arrows).

Figure 5 :

 5 Figure 5: Scans and microphotography of samples Bern2 (A, C) and 3 (B, D). (A, B) Scan of the indurated samples (A: Bern3, B: Bern2) with schematic organization of micro-facies (a to d). (C, D) Scans of the two thin-sections. (E) Bern2-3_a: channel infilling including fragmented charcoal. (F) Bern2-3_b: aggregates, mostly of biological origin, made of combustion residues rich heterogeneous material (red arrow) or of mineral components similar to the local regolith. (G) Bern2-3_c: Oxidized upper part of the micro-facies. (H) Bern2-3_b: General organization and microstructure with mixture of aggregates enriched in combustion residue (red arrow) and others similar to the surrounding regolith (yellow arrow).

Figure 6 :Figure 7 :

 67 Figure 6: A synthesis of the organization of the identified MFTs in samples (images from resin-indurated blocs cut, wetted and taken on a flatbed scanner)

Figure 8 :

 8 Figure 8: Correlation matrix illustrating the very good synchronization of live trees (a) and charcoal (b). Black color indicates the quality of correlation.

Figure 9 :

 9 Figure9: Interdating charcoal chronology with live wood chronology. The table at the bottom indicates the main statistical results of the cross-dating processes.

Figure 10 :

 10 Figure 10: Diagram of the 82 dated charcoal series presented as horizontal bars. Stars represent the number of the stratigraphic units (SU) from which the charcoal fragments come. Fragments without bark are indicated in green and fragments with bark in red. Blue points indicate the last ring produced by the tree only on charcoal with bark and yellow vertical lines indicate dates of felling. (b): Time course of use of the terrace illustrating the trend of charcoal production on the CK. Green color represents charcoal without bark and red color charcoal with bark. We can see an intensification from 1939.

  Figure 7

  Figure 10

  SU, detected only at the micromorphological thin section, is intercalated in the mineral soil. It has very low charcoal content and was not sampled and studied by dendrochronology (Fig.3B-C-D). 6 35-45 Substrate without charcoal, silty-clays colluvium (Fig. 3C, D)

Table 1 :

 1 Synthesis and description of stratigraphic units (SU) frequent fragments of charcoal. Their microstructure is essentially friable, with small, well-separated aggregates. The more massive and larger aggregates seem to correspond essentially to earthworm pellets. Fecal pellets of various sizes (earthworms, enchytraeidae) and fresh root tissues are observed. a few µm to several centimetres) are abundant. The microstructure evolves vertically from a relatively porous and friable assemblage (microstructure of separate crumbs) to a more massive top. The presence of channels their filling with faecal pellets and frequent fresh root tissue indicate marked biological activity.

	Sample Facies	Description

Table S3 :

 S3 Synthesis of the different charcoal-based chronologies used to compare the new Bernadouze charcoal chronology

	Mt. Hope	Central Nevada	Single-leaf	10.6	0.7	100
	Cortez			Pinyon	12.2	0.6	202	Biondi and Strachan (2009)
							6
	Bern 2-Bern2-3_d	Heterogeneous mixture of materials: local sediments, sometines rubefied, 2

Acknowledgements

This paper forms part of L. Fouédjeu's doctoral thesis co-supervised by V. Py-Saragaglia and D. Galop and funded by the University Toulouse 2 Jean Jaurès. Financial support for field investigations was provided by the FODYNA and TRANSYLVE projects directed by V. Py-Saragaglia and S. Burri. These projects are part of the OHM-Haut Vicdessos interdisciplinary program led by D. Galop and (co)funded by the LabEx DRIIHM, French program "Investissements d'Avenir" (ANR-11-LABX-0010) which is managed by the ANR. Special thanks are due to the Regional Archaeological Service of Occitanie (SRA), led by D. Delhoume, which supported the micromorphological study. The authors express their gratitude to the French National Forestry Office (ONF) of Ariège and the loggers who helped collect the wood disks and to V. Pescini for help during fieldwork. Special thanks are due to F. Philippe for his help and knowledge in the field, to the carpenter of the University Toulouse II, H. Come, who helped to prepare pieces of wood for the dendrochronological measurements, and to the Department of Forest Ecology of the Czech university of Life Sciences Prague (Czech Republic), led by M. Svoboda, for their support of the dendrochronological analyses in the framework of a doctoral fellowship funded by TESC PhD School and the GEODE lab. Finally, we thank two anonymous reviewers whose constructive comments led to a significant improvement of this manuscript.

List and legends of tables and figures:

Tables Table 1: Synthesis and description of stratigraphic units (SU)

Bern2-3_c

The contact with the underlying facies is diffuse. Organization is characterised by the accumulation of more or less rounded centimetric aggregates. Biological features are always present, notably through the integration of combustion residues in the faecal pellets filling the porosity and root tissues.

3 Bern 2-3

Bern2-3_b

The constituents are heterogeneous but enriched in charcoals and microcharcoals and brown to greyish in PPL. Crumby and micro aggregated microstructure. The traces of biological activity are abundant, whether in the form of channels, root tissues in the porosity or faecal pellets attributable to soil mesofauna.