Abstract
The nutrient-poor grasslands of Western Europe are of major conservation concern because land use changes threaten their high biodiversity. Studies assessing their characteristics show that their past and on-going dynamics are strongly related to human activities. Yet, the initial development patterns of this specific ecosystem remain unclear. Here, we examine findings from previous paleoecological investigations performed at local level on European grassland areas ranging from several hundred square meters to several square kilometers. Comparing data from these locally relevant studies at a regional scale, we investigate these grasslands’ spatiotemporal patterns of origin and long-term dynamics. The study is based on taxonomic identification and radiocarbon AMS dating of charcoal pieces from soil/soil sediment archives of nutrient-poor grasslands in Mediterranean and temperate Western Europe (La Crau plain, Mont Lozère, Grands Causses, Vosges Mountains, Franconian Alb, and Upper-Normandy region). We address the following questions: (1) What are the key determinants of the establishment of these nutrient-poor grasslands? (2) What temporal synchronicities might there be? and (3) What is the spatial scale of these grasslands’ past dynamics? The nutrient-poor grasslands in temperate Western Europe are found to result from the first anthropogenic woodland clearings during the late Neolithic, revealed by fire events in mesophilious mature forests. In contrast, the sites with Mediterranean affinities appear to have developed at earlier plant successional stages (pine forest, matorral), established before the first human impacts in the same period. However, no general pattern of establishment and dynamics of the nutrient-poor grasslands could be identified. Local mechanisms appear to be the key determinants of the dynamics of these ecosystems. Nevertheless, this paleoecological synthesis provides insights into past climate or human impacts on present-day vegetation.
Keywords
Introduction
One of the major lessons from paleoecology is that no fixed reference state or baseline exists, and that vegetation is dynamic in time (Birks et al., 2016). Nevertheless, defining targets for environmental management means identifying trajectories, a key issue for conservation, restoration, or sustainability (Bakker et al., 2000; Egan and Howell, 2001). This involves examining the origins of the ecosystem to determine what its trajectory was before a disturbance shifted the ecosystem onto another trajectory (Balaguer et al., 2014).
Identifying the origin of the ecosystem is particularly important in the case of nutrient-poor grasslands. In Europe, these grasslands, also widely known as semi-natural grasslands, are of major conservation importance, presenting high levels of biodiversity (Poschlod and WallisDeVries, 2002) and being threatened by decreased traditional land use (Kahmen and Poschlod, 2004). Human activities have been identified as the key disturbances triggering the presence of nutrient-poor grasslands for at least several millennia (Goepp, 2007; Poschlod and Baumann, 2010; Poschlod et al., 2009). It is thus assumed that without human deforestation and agro-pastoral activities (Bredenkamp et al., 2002), nutrient-poor grasslands would be restricted to locally small, extreme habitats like rock outcrops or extremely steep south-facing slopes with edaphic conditions too dry for tree or shrub growth (Ellenberg, 1996).
Over the last decades, the vegetation in most nutrient-poor grasslands (and in other types of grasslands) has been changing due to shrub encroachment. The species typical of the grassland are being excluded, thereby inducing a loss of species richness and diversity (Bonanomi et al., 2013; Habel et al., 2013; Kahmen et al., 2002; Vassilev et al., 2011). Once the human constraint is removed, natural secondary succession takes place (Kahmen and Poschlod, 2004; WallisDeVries et al., 2002). As a consequence of this on-going habitat change and biodiversity loss, stakeholders have made substantial efforts to conserve and/or restore this type of ecosystem (Dolnik et al., 2009; Schrautzer et al., 2009). Usually, however, not much is known about the primary succession at the spatial scale of present-day grassland distribution, nor about the mechanisms behind the development of such grassland ecosystems. In particular, the respective impacts of recurrent disturbances such as fires, windfalls, and grazing by wild herds remain relatively unexplored (Blakesley and Buckley, 2016; Svenning, 2002; Vera, 1999).
Whether nutrient-poor grasslands are a strictly man-made ecosystem or are of natural origin has been debated for decades (Baumann and Poschlod, 2008; Bush, 1988, 1989, 1993; Thomas, 1989; Wilmans, 1997). This lack of consensus on the origin of nutrient-poor grasslands has so far prevented reference states from being clearly defined (Blakesley and Buckley, 2016). However, recent studies have provided new insights into the local origin and history of various nutrient-poor grasslands in Atlantic, Continental, and Mediterranean Western Europe (Bal et al., 2014; Dutoit et al., 2009; Goepp, 2007; Henry et al., 2010; Poschlod and Baumann, 2010; Schwartz et al., 2005; Vernet, 2006).
Here, we review these studies to investigate whether the origins of the nutrient-poor grasslands show any similarities in chronology and spatial patterns through a large-scale comparison of locally relevant data. Our results are then compared with previous regional pollen data to extend local scale conclusions to the regional scale.
Methods
Studying the long-term history of nutrient-poor grassland communities (e.g. grasslands established on nutrient-poor soils) is challenging because they are generally restricted to specific small areas (e.g. hilltops). Paleo-indicators significant at a relevant spatial scale need to be collected from areas that rarely offer suitable archives recording and preserving the standard paleo-indicators, such as pollen grains. Analysis of charcoal assemblages from soil/sediment has been shown to provide suitable paleo-signals, especially relevant at local scale (Nelle et al., 2013; Thinon, 1978). Therefore, the studies hereunder relied on the analysis of macroscopic charcoal pieces accumulated and preserved locally in the soil/sediment (pedoanthracological approach). The charcoal assemblages are composed of macroscopic pieces of 400 µm or larger, recorded and preserved locally (Dutoit et al., 2009; Talon, 2010; Touflan et al., 2010). Moreover, the macroscopic charcoal assemblages are extracted either from soil horizons formed in-situ, or from soil sediment eroded and deposited in a catchment area. Thus, the spatial relevance of the charcoal data is related to the size of the catchment area (Robin and Nelle, 2014; Robin et al., 2014). Moreover, the presence of fossil assemblages of charcoal is directly related to fire disturbances that could directly explain past changes in communities.
Following Carcaillet and Thinon (1996), a similar protocol was applied to each site, both to sample the macroscopic charcoal assemblages and to analyze them qualitatively and quantitatively. For each study site considered here, the assemblages of macroscopic charcoal pieces were extracted from natural soil or soil sediment archives. Sampling exposures were dug manually to bedrock and are therefore described in terms of soil horizon and/or levels of soil sediment accumulation (soil stratigraphy; Robin et al., 2014). Between 5 and 10 L of material was sampled for each vertical unit described. These vertical sampling units were defined either according to soil horizon or stratigraphy (Goepp, 2007) or by systematic sampling every 10 cm of depth (Dutoit et al., 2009). The number of exposures per study site varied according to the availability of suitable archive site for sampling and the local sampling strategy to relevantly fit to the surface of the study area (Table 1).
General biotic and abiotic characteristics of the various nutrient-poor grasslands compared here.
includes single subsurface sampling (Vernet, 2006).
only the dated taxa.
not specifically compatibilized (Goepp, 2007).
Once sampled, each bag of material was laboratory treated to extract the macroscopic charcoal assemblages. First, the contents of the bag were dried and then weighed. Material larger than 5 mm was excluded by sieving. Next, the dry material was carefully wet-sieved with a deflocculant (Na4P2O7) through a column of sieves of three mesh sizes: 2 mm, 800 µm, and 400 µm. The charcoal pieces were then extracted from mineral and organic sections, for each size class of each sample, by hand-sorting under a binocular microscope.
Charcoal pieces were then randomly selected from the extracted charcoal assemblages for taxonomic identification, in balanced proportions from the three size classes of each sample. The charcoal pieces were identified according to wood anatomy criteria, since carbonization of wood fixes the wood anatomy so that it remains visible with only minor deformation and shrinking (Braadbaart and Poole, 2008; Théry-Parisot, 2001). The charcoal pieces were identified under an incident light microscope, mainly at magnifications of ×200 and ×500, and compared with charred wood reference collections and wood anatomy atlases (Jacquiot, 1955; Jacquiot et al., 1973; Schweingruber, 1990).
The macroscopic charcoal signals were actually further quantified per identified taxa, and/or in total amount per vertical unit of sampling, and/or in total amount per profile of sampling, and/or in charcoal concentration per kilo of sampled material smaller than 5 mm. However, the present comparative study only considered the taxonomically identified charcoal pieces that were directly dated. For each study area, chronological information was gathered from the radiocarbon dating of various quantities of single charcoal pieces, after their formal taxonomic identification. Only these latter data provide insights without reference to the complex and unclear taphonomic processes. These taphonomic processes may significantly influence the recording and preservation of soil charcoal assemblages. In particular, bioturbation causes mixing of assemblages, and fragmentation of the charcoal pieces causes quantitative bias in records (Carcaillet, 2001a, 2001b; Robin et al., 2012). Poorly documented such processes mean that extensive extrapolation is required for taxonomic and chronological interpretation of any unidentified or undated pieces in charcoal assemblages, often with questionable results (Carcaillet and Vernet, 2001; Giesecke, 2016; Robin et al., 2016; Willis et al., 2001).
Since the data on radiocarbon ages and the corresponding taxa identification for single charcoal pieces are less likely to involve extrapolation, the database for this comparative approach uses the C14 measurements of each study site and the corresponding identified taxa. Our database moreover eliminates potential comparative biases due to site-specific variations in sampling strategy at local scale. Nevertheless, of the large quantity of charcoal pieces extracted overall from the different study sites, only a small proportion were identified and dated. Thus, the level of representativeness of the identified and dated taxa at the scale of each investigated site was also taken into account (Table 2).
Radiocarbon ages of the charcoal pieces from the various nutrient-poor grasslands compared here, with information about charcoal assemblages.
: taxa with low frequency; + +: taxa with medium frequency; + + +: taxa with high frequency.
Study sites
We considered six Western European nutrient-poor grasslands, ranging widely from nutrient-poor Mediterranean grasslands to nutrient-poor subalpine grasslands in very wet areas, where soil charcoal data were available (Figure 1 and Supplementary material, available online):
La Crau (Henry et al., 2010): This plain area is situated in South-Eastern France, ca. 50 km north-west of Marseille. The plain is a fossilized river bed formed by Durance river deposits between 2,000,000 and 70,000 years ago (Molliex et al., 2013). In this flat topography, stone cover amounts to ca. 50%, and average soil depth is 40 cm, covering a 1–5 m thick layer of impermeable conglomerate. Soils of the study area (Crau de Miramas) are Truncated Leptic Cambisol originating from the dissolution and removal of carbonates after soil deposition by the Durance river between 70,000 and 125,000 years ago (Riss-Würm interglacial period). It should also be noted that soil abrasion by eolian deflation occurred in this area between 10,000 and 75,000 years ago, resulting in the erosion of the upper silty layer (Molliex et al., 2013). Thus, the alluvial water table is inaccessible to plant roots. Moreover, dryness is also induced by the Mediterranean climate, with a mean temperature of 14.5°C and a yearly precipitation of 500–600 mm (Buisson and Dutoit, 2006). Previous palynological studies realized at the regional scale (Triat-Laval, 1982) showed that during the late Neolithic and Bronze Age, the landscape was already covered by matorral vegetation and Pinus pinea/halepensis forests but not by forest dominated by the late successional species Quercus ilex. Today, the Mediterranean nutrient-poor grassland plant community, also called ‘steppe of La Crau’ or ‘Coussoul’, is characterized by the dominance of a perennial Poaceae tussock species, Brachypodium retusum, while in the western part Thymus vulgaris is the dominant plant species. Since antiquity, from 100 BC to fourth/fifth century AD in this area, the plain has chiefly been used for sheep grazing (Badan et al., 1995). No traces of shrub encroachment have therefore been observed in its central part since the 19th century (Buisson and Dutoit, 2006). A total of 14 soil profiles were sampled at different locations in the remaining area of the plain of La Crau, from which 23 single charcoal pieces were dated.
Mont Lozère (Bal et al., 2014): This study site is situated in South-Central France, in the southern part of the mountainous area of the Massif Central. Mont Lozère is a 30-km-long mountain range, oriented northwest–southeast, ranging from 900 to 1700 m a.s.l. It is a granitic area, with marked relief. De Beaulieu and Gilot (1972) put forward the hypothesis of a human origin for the grassland, consistent with pollen data showing the opening of the beech-forest landscape by agro-pastoral activities during the same period as in the Causses, in the middle Neolithic–Chalcolithic. Soils are influenced by the leucogranite substratum (silicate rock) and belong to Hyperdystric Cambisol. Mean annual temperature is 6.5°C and yearly precipitation is 1500 mm. Old and traditional agricultural systems, where nutrients were translocated through sheep grazing to arable fields, led to this nutrient-poor vegetation (Bal et al., 2014). Vegetation distribution along the slopes of Mont Lozère can be divided as follows: from 400 to 900 m a.s.l. forests are dominated by Quercus pubescens and Quercus petraea; from 900 to about 1400 m a.s.l. forests are dominated by Fagus sylvatica, with a few Abies alba; and above 1400 m a.s.l. lie Nardus stricta grasslands with Vaccinium myrtillus, Calluna vulgaris, Cytisus oromediterraneus, and Genista purgans. Sporadically, early successional species such as Pinus sylvestris, Betula pendula, and Juniperus nana are encountered. Four soil profiles were sampled in the Nardus stricta grasslands, along an altitudinal transect from 1450 m to 1670 m a.s.l. A total of 14 single charcoal pieces was selected and dated.
Grands Causses (Causse Méjean & Causse Larzac; Vernet, 2006): The site of Grands Causses is situated in the same region as Mont Lozère and includes two areas, Causse Méjean and Causse Larzac. Both are located on a Jurassic calcareous-dolomitic plateau from 700 to 1000 m a.s.l., with an active Karstic system (Corbel, 1954). Soils belonging to Dolomitic Cambisol occurred on Jurassic calcareo-dolomitic substrate. Mean annual temperature is 8–9°C and yearly precipitation 1000–1200 mm (Météo France online data). Supramediterranean forest species grow in close proximity to oromediterranean and mediterraneo-montane species (Lecompte and Cohen, 1985). Both areas have distinct zones of forest, dominated by Pinus sylvestris, and grassland traditionally grazed by sheep. Because of the lack of wetlands and peat bogs on the calcareous-dolomitic plateau of the Causses, only one reference was found concerning previous palynological investigations (Van Campo and Aymonin, 1962), simply suggesting the existence of a forest of Pinus sylvestris in the Mesolithic. Moreover, prehistoric traces have already shown that the spread of human settlement in the Causse Méjean began during the Neolithic period, with a peak during the late Neolithic–Chalcolithic. This assumption is based on data from the distribution of archeological remains and is confirmed by the results of Vernet (2006), which support the existence of several phases of fire affecting pine forests in the northern areas of the Causse during the period 4800–3000 BP cal. yr. Soil and sediment samples in the Grands Causses were taken from natural small topographical depressions, such as dolines, at various locations on the plateau. To ensure the comparability of the data on grassland history, only samples from open areas were taken into consideration here. A total of 216 samples were collected from 85 subsurface locations and yielded a total of 28 dates (26 from Causse Méjean and 2 from Causse Larzac).
Vosges Mountains (Rossberg and Hohneck; Goepp, 2007:; Schwartz et al., 2005): The Vosges are a range of low mountains in North-Eastern France, ca. 120 km long, oriented north–south and formed of sandstone and crystalline stones. Their altitude ranges from ca. 300–1424 m a.s.l. Average annual temperature is 8°C (Gerardmer station, Météo France), and yearly precipitation is from 1800 to 2000 mm (Météo France on-line data). Janssen and Janssen-Kettlitz (1972) showed a slight increase of pollen from non-forested herbaceous species in the Bronze Age around 2500 yr BP, but no signs of the existence of former open areas have been detected. Soils from the Vosges sites belong to Hyperdystric Cambisol and Umbric Leptosol originating from the degradation of biotite porphyroid granites under mountainous climate conditions. The grasslands are mainly in the southern half of the massif, where the elevation of the relief leads to what are known as ‘chaumes’, above ca. 1000 m a.s.l. The vegetation communities of these grasslands include mainly Nardus stricta, Festuca ovina, and Festuca rubra. They have been grazed by cattle for centuries. Some dwarf shrubs also occur, such as Vaccinium myrtillus and Calluna vulgaris. Sites were investigated on two chaumes: the Rossberg site, with two soil profiles providing 18 radiocarbon dates, and the Hohneck site, with 11 soil profiles providing 22 radiocarbon dates (Schwartz et al., 2014).
Franconian Alb (Baumann, 2006; Poschlod and Baumann, 2010): This study site includes several sampling areas in the central and southern Franconian Alb (or Franconian Jura), in South-Eastern Germany. It is a relatively low mountain range, part of the Jura, reaching ca. 600 m a.s.l. and mainly composed of limestone from the Jurassic period. Poschlod and Baumann (2010) in the same survey compared two different methodological approaches (soil charcoal and pollen analyses) on these sites. Their results show that the peaks of Pinus pollen fully correspond to the periods to which the charcoal found in the calcareous grasslands were dated. Soils of the Rendzic Leptosol type occurred on Jurassic Malm. Mean annual temperature is around 7–8°C, and yearly precipitation is 600–700 mm (Baumann, 2006). The main phytosociological association is Gentiano-Koelerietum. A total of 17 profiles from four different sampling areas in formerly grazed grasslands were investigated, and eight radiocarbon dates were obtained.
Upper-Normandy (France; Dutoit et al., 2009): The study site is the chalk hillsides of Saint-Adrien situated ca.100 km north-west of Paris, in the Lower Seine Valley of Upper-Normandy. The Seine River and its tributaries, such as the Bequet River, shaped the geomorphology of cliffs and rocky outcrops on chalk from the Cretaceous period, with hills ranging from 50 to 110 m a.s.l. Previous palynological studies in the Seine valley showed that forest openings began to occur in the Pre-Boreal Period but that this extensively forested landscape dominated by Fagus sylvatica was probably not conducive to large areas of chalk grassland until the late Neolithic (Huault, 1986). The soils of the studied site are Rendzic Leptosol originating from the degradation of the chalk outcrop since the beginning of the post-Glacial period 9,000 years ago. The local topography of steep slopes induces shallow soil, with low water-availability for vegetation. Thus, vegetation is nutrient-poor dry grassland dominated by Sesleria albicans and Bromus erectus, grazed by sheep, goats, and cattle for centuries (Dutoit et al., 2003). Mean annual temperature is around 10°C, and yearly precipitation is 800 mm. Five profiles were sampled, one at the base of the hillside and four at the top. The latter four profiles are in an area comparable to the minimum sampling area required to provide a good representation of the floristic composition, richness, and diversity of chalk grassland communities (around 500 m²; Dutoit et al., 2003). In total, 14 single charcoal pieces were dated.

Map showing location of the six compared nutrient-poor grasslands in Western Europe (map under GNU-free documentation license).
Results
Soil charcoal synthesis
A total of 127 single charcoal pieces from the six study sites were radiocarbon dated (Bal et al., 2014; Dutoit et al., 2009; Goepp, 2007; Henry et al., 2010; Poschlod and Baumann, 2010; Schwartz et al., 2005; Vernet, 2006; Table 2). When considered together, the 14C-AMS data compiled show an uneven distribution, both within and between study areas (Figure 2). Obviously, it can be postulated that some soil charcoal pieces could have been lost through charcoal fragmentation and/or soil erosion (Robin and Nelle, 2014). However, the absence of charcoal dates for some periods, for instance, from La Crau area before 5000 cal. BP, or before 2000 cal. BP in the Franconian Alb, cannot be attributed to the absence of soils, as they were forming or already developed during the early Holocene. While it is true that the soil of La Crau area has been truncated (Molliex et al., 2013), that occurred between 10,000 to 70,000 years ago and not during the period for which dates are lacking. To avoid any over-interpretation, the chronological data are not extrapolated here to the total quantities of charcoal per soil layer; we interpret the chronological data exclusively with reference to the single charcoal pieces dated from each study site.

Age distribution of the dated charcoal pieces from the six compared nutrient-poor grasslands in Western Europe, related to the chrono-cultural phases with which the dated charcoal pieces are contemporaneous.
A clear clustering of dates is only observed during three periods at La Crau: the end of the Neolithic/Bronze Age, the Roman period, and the Medieval period/early Modern period. In Mont Lozère, clustering is less marked: some dates are found loosely grouped in the early Neolithic, others in the late Neolithic, then in the Iron Age, with a few isolated dates during the Roman and Medieval periods. The radiocarbon dataset for Grands Causses shows more dates during the Neolithic, the Bronze Age, and the Modern period, with only a few dates during the Medieval period. The two Vosges areas show differences in the temporal distribution of dates: Hohneck has a few dates occurring together during the early Iron Age, with others during the Medieval period and the Bronze Age, while the Rossberg dates are more concentrated in the Iron Age and the Roman period, with fewer Medieval period dates. Only the Franconian Alb has no dates older than Roman times. Seven dates relate either to the Roman period or to the late Medieval period/early Modern period. The Upper-Normandy dates are more or less grouped in the late Neolithic/Bronze Age and Medieval period/early Modern period.
The taxonomic spectrum of the dated charcoal pieces from the various study sites show 30 different taxa in total. The large number of taxa revealed by charcoal analysis in these study areas, compared with single-site or regional studies (e.g. Robin et al., 2013, 2014), is due to their very diverse biogeographic contexts. Therefore, to ensure comparability between study areas, the dated taxa are grouped here by functional type, according to their life-history traits with respect to plant community succession as described in the specialist literature (Denslow, 1980; Rameau et al., 1999, 2005, 2008). Some are early successional taxa, which could indicate a possibly short-lived expansion of shrubs and trees: sclerophyllous Quercus spp., Pinus pinea/halepensis and Pinus type sylvestris (including Pinus nigra, Pinus sylvestris and Pinus uncinata), Cupressus spp., Prunus spp. Phillyrea spp., Acer spp., Genista spp., Populus spp., Juniperus spp., Fraxinus spp., Vaccinium spp., Calluna sp., Corylus sp., Buxus sp., Betula spp., Salix spp., while some characterize strictly open habitats (Lavandula sp., Thymus sp.). A second group consists of taxa tolerant to competition, with Taxus sp., Sorbus spp., and Tilia spp. The last group consists of late successional taxa, with Fagus sp., deciduous Quercus and Abies sp., which are trees that could close the canopy and restrict plants of open vegetation.
The greater number of taxa in the early successional group can be explained by the selection process. Charcoal pieces were usually selected, so as to provide dating as accurate as possible with respect to the fire event. Thus, short-lived species were chosen, and these species are often early successional. Therefore, the dating of the successional groups is not interpreted quantitatively (i.e. by frequency of taxa occurrence) but only qualitatively (i.e. by site-to-site comparison of the chronological occurrence of successional groups). The taxonomic results interpreted according to chronological data are shown in Figure 3.

Number of dated charcoal pieces on a time axis per investigated grassland in Western Europe, and per successional group of the identified taxon, related to the chrono-cultural phases with which the dated charcoal pieces are contemporaneous.
The plain of La Crau mainly contained early successional species like Pinus pinea/halepensis, Phillyrea sp. (cf. angustifolia), and Prunus sp. at the end of the Neolithic or during the Bronze Age. For these periods, there is no evidence of late successional species. The only dated late successional taxon is from the Roman period. Mont Lozère had early successional taxa before La Crau, notably Prunus sp., as well as late successional taxa such as Fagus sylvatica. Both were already present during the early Neolithic. Then, both groups are dated at various periods during the middle and late Holocene. On the Grands Causses, several dates for the early successional taxon Pinus type sylvestris are from the late Mesolithic/early Neolithic (as early as the dates for the early successional taxa on Mont Lozère). The dates for late successional taxa are much later, starting with Fagus sylvatica in the late Bronze Age. The subsequent late successional taxa all date from the Modern period. Early successional taxa also show several dates from the Modern period; however, there are more prior occurrences, during the Bronze and Iron Ages. The Vosges contains the oldest dated taxa of the entire dataset, with early successional taxa such as Pinus type sylvestris, and Salix spp. from the Mesolithic. It also contains both a tolerant taxon (Tilia spp.) and a late successional taxon (Fagus sylvatica) dated to the late Mesolithic. Then, from the late Neolithic onward, various taxa from both late and early successional groups occur regularly, up to the Modern period. In the Franconian Alb, the only taxon from the early successional group dated was Pinus type sylvestris, from the Roman period to the Medieval period and the Modern period. However, observations of the charcoal assemblages reveal late successional taxa such as Fagus sylvatica, Abies sp., and deciduous Quercus, along with early successional taxa such as Juniperus and Ligustrum. In Upper-Normandy, no early successional taxa are dated. Thus, although they were most likely present locally at some stages of the forest dynamic, the absence of dating and their weak occurrence within the local soil charcoal assemblages suggest that this vegetation stage did not burn locally. The oldest dating is for late successional taxa from the late Mesolithic; they subsequently occur sporadically in more recent periods, along with tolerant taxa such as Taxus baccata.
Discussion
The presence of macroscopic charcoal pieces in natural soils indicates the importance of fire as a disturbance in the opening of the ecosystem, as pointed out by the local scale investigations at each site (Bal et al., 2014; Dutoit et al., 2009; Goepp, 2007; Henry et al., 2010; Poschlod and Baumann, 2010; Schwartz et al., 2005; Vernet, 2006). From Prehistoric times until the beginning of the Modern period, our synthesis on soil charcoal data shows that nutrient-poor grasslands existed throughout the study region. However, caution should be used in interpreting these data, since the charcoal proxy here comes from soil/soil sediment archives. While it was possible to identify the geographical origin of the charcoal signal here by the catchment area of each sampling, taphonomical processes make recording and preservation of such charcoal assemblages inadequate for reliable interpretation. This is an issue particularly regarding transfer of the results to a larger geographical range, such as the regional scale. For instance, there is no relation between the depth of sampling and the measured radiocarbon ages (Carcaillet, 2001a, 2001b). This is why our study sites simply represent a few areas in different landscapes from Northern to Southwestern Europe. Moreover, the spatial distribution of this soil charcoal sampling was designed to capture reliable local information: taxonomic identifications and radiocarbon dating, independently of the archive, to permit their comparison at regional scale. Actually, taxonomic identification and radiocarbon dating yield only partial information, since we have no guarantee that we dated the oldest charcoal pieces of a specific taxa or that we identified all the taxa historically present in the study areas (probably overlooking taxa that did not burn locally for certain chronological periods). Nevertheless, they provide some evidence of the local presence of formally identified taxa for a dated period on the site sampled.
The 14C age distribution points to fire activity occurring locally during the Mesolithic, and from the Neolithic at various sites, but still sporadically. This shows the early role of fire in the dynamics of nutrient-poor grasslands. From the late Neolithic to the Bronze Age, fire activity progressively spread throughout the study sites, consistent with fire history patterns at continental scale in Europe (Carcaillet et al., 2002; Power et al., 2008) and in the Mediterranean (Vannière et al., 2008). However, we were unable to identify clear synchronicities between fires events dated at local scale in different regions. Thus, the fire events in nutrient-poor grasslands at local scale are highly heterogeneous, probably due to the existence of local semi-open mosaics of fuel and to specific local scale fire regime drivers, like local vegetation dynamics, human activities, and potential effects from recurrent grazing by wild herds in the Mesolithic (Robin et al., 2012). However, during this period, the impacts of wild herds on the degree of vegetation opening (in space and time) could not be reliably assessed from natural bone archeological traces (Svenning, 2002). Thus, the key impact of wild herds on vegetation openness and maintenance of nutrient-poor grasslands during the Mesolithic periods in Europe is still purely a hypothesis (Vera, 1999), as the precise dynamics of wild herds have not yet been determined.
However, identifying the ignition sources of the dated charcoal is challenging. The soil charcoal database does not enable us to assess directly whether fire activity is related to human activities or to climate forcing, since the fire responsible for the presence of charcoal assemblages can be either natural or human-made. In addition, the charcoal recorded might result from the domestic use of wood for cooking, even out of archeological archives or be related to clearly identified human occupation sites such as La Crau plain (Henry et al., 2010). However, the taxonomic identification and radiocarbon dating of these charcoal pieces possibly from ‘hearths’ still provide direct evidence of the local presence of some taxa for some periods. From the local charcoal assemblages, therefore, two groups of sites can be distinguished: one at northernmost latitudes (i.e. Upper-Normandy, Vosges, Franconian Alb), and one with Mediterranean affinities (i.e. Grands Causses, Mont Lozère, La Crau).
For the first group, the impact of man-made fire is plausible considering that the burnt vegetation is identified as originating from low flammable forest systems since at least the Neolithic for the Vosges and Normandy sites (i.e. presence of Fagus sylvatica, Abies sp., and/or deciduous Quercus; Dutoit et al., 2009; Goepp, 2007; Poschlod and Baumann, 2010; Schwartz et al., 2005; Figure 3). The role played by human disturbance in the emergence of grassland ecosystems has already been pointed out by paleo-investigations of other European grassland sites (e.g. Bush, 1988, 1989). These studies feed into the debate on whether human disturbance should be seen as causing the development of grassland or as maintaining existing grassland inherited from cold periods (Bush, 1993; Thomas, 1989). The soil charcoal dataset for the sites considered here supports the hypothesis that nutrient-poor grasslands in Western Europe are the result of woodland opening since the Neolithic, which became even more marked during the Bronze Age, Roman Empire, and Medieval Age. This is consistent with the Holocene vegetation history of Western Europe: reestablishment of closed canopy forest during the early/middle Holocene (ca. 8–4 ka yr BP), followed by progressive afforestation related to human activities and use of natural resources (e.g. Berglund et al., 1996; Kalis et al., 2003; Robin et al., 2014). This is especially marked for the late Holocene (ca. the last 4 ka yr BP), which saw the onset of human use of fire as a tool to open wooded areas (Marlon et al., 2008), related to the development of cultural landscapes (Gaillard et al., 2007).
For the second group (Grands Causses, La Crau), with the sites at lower latitudes, fires may have occurred naturally due to the presence of highly flammable vegetation (e.g. pine woodland; Bal et al., 2014; Henry et al., 2010; Vernet, 2006; Figure 3). The importance of wildfire related to summer drought as a key disturbance during the Holocene in the Mediterranean region has already been shown (Vannière et al., 2011), supporting climate determinism (Combourieu-Nebout et al., 2013; Jalut et al., 2009; Pons and Quezel, 1998; Thiebault, 1997). In La Crau, the charcoal assemblages show a degradation of pre-existing matorral vegetation and pine forests by fire and grazing during the Neolithic, before any advanced forest succession established locally (Bazile-Robert, 1987; Henry et al., 2010). In the Grands Causses sites, the same pattern emerges even earlier, with fire occurrences in pioneer forests from the beginning of the Mesolithic and the Neolithic (Bal et al., 2014; Vernet, 2006). This early successional vegetation type remained dominant during the Neolithic, probably because the Mediterranean climate and soil erosion slowed primary plant successions (Carrión et al., 2003; Grove and Rackham, 2003). Human activities intensified the effect of climate constraint, keeping open ecosystems dominated by herbaceous plant communities, notably through grazing (Henry et al., 2010; Vernet, 2006).
A relatively similar phenomenon has been observed in Central Europe, where natural grasslands occurred even before the start of agricultural activities in the early Neolithic (5500 BC). However, their area was estimated to be relatively small and grassland patches were probably rather fragmented in a dominantly forested landscape matrix. Thus, large-scale enlargement of grasslands is not expected to have occurred before the Iron Age (Hájková et al., 2011; Hejcman et al., 2012; Kuneš et al., 2015; Pokorny et al., 2015).
Finally, on Mont Lozère, the finding of some Fagus sylvatica dates from the Neolithic points to the presence of mature stands before the establishment of grassland. However, for the same period and from the same charcoal assemblages, pioneer forest vegetation was also found (i.e. Betula spp., Salix spp., Corylus avellana; Bal et al., 2014). Thus, it can be postulated that the Fagus dates do not result from dense closed forests. The high altitude of the Mont Lozère grasslands sampled, with their wet mountain climate, explain the differences from La Crau and even from Grands Causses, where climate conditions are much drier. However, for those three study sites with Mediterranean affinities, it is clearly shown that the first fire events occurred in non-mature forest stages, as revealed by charcoal assemblages. This is an important distinction, marking them out at large scale from the study sites at northernmost latitudes, where the first fire occurrences were in late successional forest stages.
Concerning the last few hundred years, the absence of charcoal from this period in La Crau, Mont Lozère, and the Franconian Alb may indicate a continuous openness of these grasslands. This may have led to higher floristic diversity for sites not plowed during the Modern period (Dutoit et al., 2003; Forey and Dutoit, 2012; Gibson and Brown, 1992; Karlík and Poschlod, 2009; Römermann et al., 2005). Grassland diversity has already been linked to late Iron Age human population density (Pärtel et al., 2007). Conversely, the abundance of young charcoal from Fagus should indicate that the other grasslands are younger and would quickly disappear if human land use stopped. These grasslands would certainly differ floristically. However, since our sites are situated in different biogeographic areas, testing this hypothesis would require discriminating climatic factors from those linked to history and landscape organization at the regional and local scale across Europe (Cousins and Eriksson, 2002).
Conclusion
Comparing the locally relevant chronological and taxonomic data for the six study areas revealed no clear pattern of nutrient-poor grassland development and dynamics between areas in contrasting parts of Europe. Local-scale mechanisms appear to be the dominant drivers of these ecosystems. It can be concluded that the origin and history of European nutrient-poor grasslands are related to a combination of climatic constraint and/or human influence impacting vegetation dynamics at local or at least ‘infra-regional’ scale. Nevertheless, we observed a clear distinction between the sites under Mediterranean influence (i.e. Grands Causses, Mont Lozère, and La Crau), where matorral and pioneer forest stages were already established when the first human impacts occurred during the Neolithic, and the sites at central Europe latitudes (i.e. Upper-Normandy, Vosges, and Franconian Alb), where the grassland appears to result more from anthropogenic woodland clearing from the Neolithic onwards. Moreover, comparisons between local pedoanthracological data and previous regional pollen data clearly show how the results can be additional for a better understanding of the spatial areas really affected by vegetation changes. Here, charcoal-based data suggest a Prehistoric origin for nutrient-poor grasslands (from the Neolithic to the Bronze age) but their extension in space appears to have occurred no earlier than the Middle-Ages according to most of the regional pollen analyses previously realized in our study regions.
