The initial settlement of humans in a location, in particular “slash and burn” land clearance practices associated with agricultural development, is habitually detected in paleoecological records as an increase in macroscopic charcoal particles (e.g., Glikson, 2013). Charcoal analysis provides insights into wood-burning practices, but is restricted by the size of identifiable particles (Lehndorff et al., 2014). In addition, this approach has spatial and temporal limitations, is affected by secondary deposition up to several years later, and does not directly detect fire intensity (Denis et al., 2012). A number of molecular markers such as levoglucosan and benzene polycarboxylic acids (BPCA) have been proposed recently for the quantification of pyrogenic carbon in fire-impacted soils and sediments (Elias et al., 2001; D'Anjou et al., 2012; Lehndorff et al., 2014).

The signatures and concentrations of organic compounds in smoke aerosols from the controlled burning of agricultural, silvicultural, and garden vegetation have been shown to vary according to burn temperature, flame structure, pile density, moisture content, and vegetation/biomass type (Freeman and Cattell, 1990; Oros and Simoneit, 2001; García-Falcón et al., 2006; Denis et al., 2012). Since the biomarkers and their combustion alteration products are source-specific, they can be used as tracers for determining fuel type and source contributions from biomass burning. Oros and Simoneit (2001) showed that the major organic components directly emitted to smoke particulate matter from conifers were straight chain aliphatic compounds from vegetation wax and diterpenoid acids from resin. The major natural products altered by combustion included derivatives from phenolic (lignin), monosaccharide (cellulose) biopolymers, and oxygenated and aromatic products from diterpenoids. Finally, other biomarkers present as minor components included phytosterols, unaltered high molecular weight wax esters and polycyclic aromatic hydrocarbons (PAH).

Vane et al. (2013) assessed the transfer pathways of PAH between upland peats subject to managed burning and reservoir sediments in South Yorkshire, England. These authors observed that the major transfer pathway is in the particulate form (i.e., suspended sediment) rather than the dissolved form. In addition, the PAH signatures of the recently formed, shallow peat samples were strongly correlated with those of the sediment core samples, suggesting that the former are a dominant source for the latter. Similarly, Denis et al. (2012) showed that PAH fluxes in annually varved lacustrine sediments from the Yosemite National Park (USA) record the observational fire history within 0.5 km of the lake. Low molecular weight PAHs (e.g., fluoranthene, pyrene, and benz[a]anthracene; Supplementary Figure 1) were found to be the best recorders of fire, whereas high molecular weight PAHs recorded fire intensity.

PAHs can also be derived from vegetation and its decayed remains during early diagenetic processes (Wakeham et al., 1980; Malawska et al., 2002) or from atmospheric deposition of industrial emissions (Sanders et al., 1995). A simultaneous rise in phenanthrene (PAH) in Lake Michigan (USA) and Lake Ontario (USA) around 1900 has for instance been attributed to high temperature combustion of coal for steel production (see Figure 11 by Meyers, 2003). The correspondence between the phenanthrene maxima in the 1950s and use of coal suggests that airborne distribution of industrial stack gases and particles has distributed PAHs across the Great Lakes region with, however, significant contributions from local sources.

Schüpbach et al. (2015) combined macroscopic charcoal and the molecular fire markers levoglucosan, mannosan, and galactosan to investigate Holocene fire history in the Maya Lowlands of Petén (Guatemala). The authors detected three periods of high fire activity during the Holocene: 9500–6000, 3700, and 2700 cal year BP. The first period of high fire activity was attributed to climate conditions, but interestingly, all three molecular fire indicators show an earlier increase compared to charcoal for the first peak. The authors suggest that the earlier increase of the molecular markers might be a result of low temperature and smoldering fires prevailing over the entire region of Mesoamerica during the early Holocene. This signal would have been captured less efficiently by the sedimentary charcoal record, which is mostly a local fire indicator. The fire maximum around 3700 cal year BP might already be influenced by early anthropogenic agricultural activities in the region, based on the fast decrease in the relative abundance of wooden charcoal (Schüpbach et al., 2015). The youngest fire maximum (2700 cal year BP) was linked to the agricultural activity of the Maya at Lake Petén Itzá. It lies well within the “Maya Clay layer,” which corresponds to the interval with the maximum soil erosion rate (Anselmetti et al., 2007).

This biomarker approach might prove useful to distinguish between wild fires that flare during drought years and intentional fires lit independently of climate conditions. Analyses of 11 PAH compounds in the particulate matter produced by burning various Australian woods and leaves under different conditions showed that large scale uncontrolled bush fires produced much higher relative amounts of benzo(a)pyrene, the most toxic of the parent PAH compounds, than other combustion sources (Freeman and Cattell, 1990).

Pollen has long been used to identify early human impacts on the landscape. It is possible to evaluate characteristic vegetation shifts resulting from land clearance and to identify pollen from cultivated plants or weeds associated with human-disturbed landscapes (Behre, 1981). Lignin phenols and plant waxes (long chain n-alkanes and n-alkanoic acids) are molecules characteristic of vascular plants, and therefore can serve as biomarkers for terrestrial plants (Feng et al., 2013). Their resistance to biodegradation makes them excellent biomarkers as they pass into the environment in and on leaf fragments and as aerosols formed by wind abrasion and fires (Eglinton and Eglinton, 2008). Plant waxes can be transported from land to their sedimentary repository in lakes and oceans within days to weeks (Eglinton and Eglinton, 2008). The chain length distribution and their carbon isotopic composition reflect the vegetation present, in particular the C3–C4 plant distribution (Chikaraishi et al., 2004; Rommerskirchen et al., 2006).

n-Alkane distributions in sediments provide an approximate record of the kinds of plants that have populated the watershed (Meyers, 2003). In watersheds where grasses dominate, C31 is the major sedimentary n-alkane, whereas C27 and C29 are more abundant in lake sediments where trees dominate (Cranwell, 1973). The average chain length (ACL) is commonly used to evaluate the relative contribution of trees, shrubs, and herbs (Peltzer and Gagosian, 1989): ACL27–33 = Σ([Ci])×i)/Σ[Ci] where [C_i] is the concentration of the n-alkane with carbon number i, over the range 27–33.

Ternary plots of the relative proportions of the C27, C29, and C31 n-alkanes in sediments can be used to reconstruct past abundances of grasses, shrubs, coniferous, and deciduous trees in the vicinity of a lake (e.g., Schwark et al., 2002).

Fisher et al. (2003), for instance, showed an extremely close correspondence between the mean carbon number of (C27–C31) n-alkanes and the percentage of pollen derived from grass in the sediments of Gormire Lake (UK). Thermochemolysis of lignin confirmed that there was a significant change from woody to grass-dominated material deposited in the lake sediments, contemporaneous with a first deforestation event at ~600 BC. Changes in both the sedimentary C/N ratios and yields of lignin-derived phenols suggest that deforestation events at ~600 BC and AD 1200 led to enhanced aquatic productivity, initially through the development of reed-swamp vegetation.

A less typical example comes from a stalagmite record from Ethiopia. Preserved stalagmite plant-derived n-alkanes and n-alkanols show clear changes in composition over time, related to known land-use changes in the area (Blyth et al., 2007). Clearance of the original shrub and woody vegetation in the 1930s to make way for agriculture is reflected as a decrease in the C27/C31 n-alkane ratio, which is a simple expression of the ACL (Marseille et al., 1999).

More specific molecular markers of vascular plants have been proposed to track changes in terrestrial vegetation. Tricyclic diterpenes (Supplementary Figure 1) are considered to be mainly derived from gymnosperms (Simoneit et al., 1985), whereas pentacyclic triterpenes with oleanane, ursane, or lupane structure (Supplementary Figure 1) are mostly produced by angiosperms (Abe et al., 1993). Other pentacyclic triterpenes, such as those with a fernane or arborane structure (Supplementary Figure 1), are distributed between angiosperms and ferns, and potentially bacteria (Hauke et al., 1992a,b; Jaffé and Hausmann, 1995). Within major vascular groups, even more specific molecular biomarkers, mostly pentacyclic triterpenes, have recently been developed (see next section).

Forest clearings, land-use practices, and grazing by livestock often diminish catchment soil stability, causing detrital events that are reflected in sedimentary basins as magnetic susceptibility maxima, higher amounts of detrital elements (e.g., Dearing, 2008; Debret et al., 2010; Arnaud et al., 2012; Magny et al., 2013) higher sedimentation rates, and changes in grain size distribution. However, higher erosion rates can clearly also be associated with increased precipitation. Thus, the link to human activities needs to be made based on additional indicators.

Recently, Jacob et al. (2009) proposed a new molecular biomarker of soil erosion: Trimethyl tetrahydrochrysene (TTHC) a diagenetic derivative of oleanane-type pentacyclic triterpenes (Wakeham et al., 1980). TTHC profiles reconstructed in Lake Le Bourget (France) sediment revealed increased soil erosion linked to the arrival of new agricultural practices during the Bronze Age in the lake catchment (see also Section Cultivation, Figure 1). This compound was similarly used by Guillemot et al. (in review) to track soil erosion due to recent agropastoral pressure around Lake Igaliku (Greenland). However, the assumption that TTHC arises from degradation of pentacyclic triterpenes in litters and soils (Wakeham et al., 1980) remains to be demonstrated, since we have not yet been able to detect TTHC in several litter and soil samples.

Another type of biomarker, branched glycerol dialkyl glycerol tetraethers (brGDGTs) are regarded as specific biomarkers for soil, as they are allegedly formed by soil bacteria (Hopmans et al., 2004). The Branched and Isoprenoid Tetraether (BIT) index was thus created to reconstruct the contribution of soil organic matter to marine and lacustrine sediments. It is based on the relative abundance of brGDGTs, representing terrestrial organic matter, and crenarchaeol, representing aquatic organic matter (Hopmans et al., 2004). Although brGDGTs were also found to be produced in situ in lakes, means exist to determine the significance of in situ production (Naeher et al., 2014). The BIT index could thus represent a potential marker for increased soil erosion caused by human land-use, but to our knowledge it has never been applied as such.

Organic matter in lacustrine sediments is derived from production within the aquatic system, inputs of terrestrial material from the surrounding lake catchment and bacterial production within the sediments themselves (Prartono and Wolff, 1998). The relative importance of these sources is determined by local factors such as climate, nutrient supply, hydrodynamic conditions, and the biogeochemistry of the lake sediments. Human perturbations in a catchment often result in an increase in primary productivity in the aquatic ecosystem downstream. In sedimentological records, such increases are reflected in the C:N ratio as well as the δ¹³C of organic matter, which trace the relative contributions of terrestrial and aquatic organic matter (Massa et al., 2012). In addition, diatom-based trophic reconstructions have been used to track the impacts of past land use on water quality (e.g., Miras et al., 2015).

Algal and microbial biomarker records can also reflect changes in aquatic ecosystems. Whereas odd long chain n-alkanes (C27, C29, C31) are major constituents of terrestrial plants (see Section Vegetation shifts), short chain n-alkanes (C15, C17, C19) are produced by many phytoplankton species (Cranwell, 1982). The abundance of short chain n-alkanes has thus been used to reconstruct lacustrine paleoproductivity rates (Meyers, 2003), and the ratio of terrigenous to aquatic n-alkanes or TAR_(HC) was defined by Bourbonnière and Meyers (1996) to determine the relative contribution of terrestrial vs. aquatic OM:

Silliman and Schelske (2003), for instance, documented a dramatic drop in the TAR_HC, which indicates a shift from macrophyte to algal dominance of organic matter production, in Lake Apopka (USA). This trophic shift followed agricultural development of the watershed in the mid-1940s and reflects increased nutrient delivery to the lake.

Similar to n-alkanes, sterols can be used to evaluate the proportion of allochtonous to autochtonous organic matter. In this case, the ratio of C₂₉/C_27(ST) is used, where C29 is mostly synthesized by terrestrial higher plants (as β-sitosterol), whereas the C27 sterols are derived from cholesterol, which has essentially an aquatic origin (Nishimura and Koyama, 1977; Muri et al., 2004). Prartono and Wolff (1998), for instance, observe increases in concentrations of unsaturated and branched chain fatty acids, phytol, and C27 and C28 sterols toward the surficial sediments of lake Rostherne Mere (UK), suggesting that the autochthonous contribution to sedimentary organic matter has increased as a result of the recent eutrophication. Isotopic carbon composition and C/N ratios support the increasing inputs of autochthonous organic matter toward surficial sediments. Interestingly, the authors note that elevated levels of relatively reactive autochthonous organic matter in the shallower sediments and the resulting reducing conditions in the water column have led to a change in the diagenetic reactions prevailing in the sediments.

Last but not least, the proportion of allochtonous vs. autochtonous organic matter and the redox potential in sediments has been suggested to control the microbial reduction of sterols in sediments (Nishimura, 1977). The 5α(H)-stanols/Δ⁵-sterols ratio is thus used to monitor changes in redox potential, and has been defined as:

As an example, Bertrand et al. (2012) used a suite of molecular ratios to identify the different developmental steps of a watershed, from a swamp to the progressive settlement of a pond in the Lorraine region (France). Several transitions are clearly recorded, especially drying or high water level phases, as well as the extension of the forest cover in the watershed. Interestingly, draining and dredging of the sediment in 2002 resulted in a redox change (an oxic event) 15 cm below the 2002 horizon.

In palynological records, the appearance of Cerealia pollen (Behre, 1981) or Zea mays pollen in Mesoamerica (e.g., Pohl et al., 2007) is often used to detect the beginning of cultivation practices. The biomarker approach can be especially useful in places were pollen are poorly preserved, or for plants that are difficult to identify in palynological records. For instance, whereas the distinction between Cannabis and Humulus (hop genus) pollen grains requires morphometric analyses due to their resemblance (Whittington and Gordon, 1987; Mercuri et al., 2002), cannabinol allows a clear differentiation, as this organic molecule is specific to hemp (Lavrieux et al., 2013). The palynological record of cultivated taxa may also be skewed in tropical Pacific insular ecosystems where settlers consistently cut the flowers of introduced plants (taro and yam) to simulate tuber growth (Combettes et al., 2015).

The molecule miliacin is a prime example of such an approach. Jacob et al. (2008a,b) demonstrated that sedimentary miliacin (olean-18-en-3β-ol methyl ether) was the only PTME specific to the broomcorn millet. Bossard et al. (2013) later confirmed that seeds of broomcorn millet are sedimentary miliacin's dominant source, thus its detection in lacustrine sediments indicates that this PTME is produced in significant amounts in the catchment. Jacob et al. (2009) used miliacin to track the cultivation of broomcorn millet (Panicum milliaceum) in the catchment of Lake le Bourget (France, Figure 1). As millet is an allochtonous plant introduced by humans in the region, this organic molecule is unambiguously indicative of human activities. Their records show that the onset of millet cultivation coincides with the onset of major soil erosion in the catchment during the Middle Bronze Age (~1800 BC; Figure 1). It is worth noting that the beginning of millet cultivation was not associated with major vegetation changes in the catchment (Gauthier and Richard, 2009; Jacob et al., 2009). Only pollen data from the Late Bronze Age showed significant disturbances. Millet fields were probably restrained in area (the catchment is very steep and colonized by conifers), and probably not larger than those used during the Neolithic. But these highly productive hot spots of miliacin production in the catchment were efficient enough to be recorded in the sediment.

Simonneau et al. (2013) found an almost synchronous first appearance of miliacin in the sediments of Lake Paladru (France, ~1600 BC), revealing that millet cultivation was a regional practice but diachronously imported in the Alps. Miliacin has also been detected in Ukrainian paleosols, where it was likewise interpreted as a tracer of former millet cultivation (Motuzaite-Matuzeviciute et al., 2016). In this particular case, the use of molecular biomarkers helped to disentangle the origin of charred cereal grain remains in pits that were not associated with any local settlement. High levels of miliacin in ancient soils stratigraphically connected to the pits proved that millet was cultivated locally and not being brought from other locations by mobile groups (Motuzaite-Matuzeviciute et al., 2016).

Organic geochemical studies of soils under experimental control and paleosols have shown that compound specific δ¹⁵N of amino acids reveal early land use practices (Simpson et al., 1999b). Despite the rapidly decomposing nature of the amino acid fraction in soils, the signal demonstrates differences between manured grassland, unmanured grassland and continuous cereal cultivation (Simpson et al., 1999b). Amino acids such as valine, alanine, leucine, and isoleucine are isotopically lighter where grassland soils have not been manured in comparison to manured areas. In addition, δ¹⁵N values are lighter under cereals and heavier under grassland. Relict soils from the Mediaeval Period and fossil soils from the Bronze Age from Orkney retain a δ¹⁵N signal typical of heavily manured soils.

In more recent times, Zocatelli et al. (2010) used changes in grass specific biomarkers to trace turf production around Lagoa do Boqueirão (Brazil) for football stadiums and golf courses. The increase in PTME concentrations in the lake sediments during the period 1996–2000 results mainly from high amounts of a compound tentatively assigned as arbor-8-en-3βol ME, a potential diagenetic derivative of cylindrin produced by Zoysia japonica, a grass introduced for turf cultivation (Figure 2).

Novel biomarkers specific to introduced species could thus provide additional indicators for the domestication and spreading of certain cultivation. Of interest could be the iso- and anteiso-monomethyl alkanes reported in plants from the mint family (see Section Vegetation shifts above). In addition, biochemical investigations can provide useful information on past land-use, especially when settlements or macrobotanical remains are absent (Motuzaite-Matuzeviciute et al., 2016).

Hemp, along with millet, is one of the oldest cultivated plants. Ancient societies used hemp for food, medicine, and to produce ropes, paper, and clothes (Russo, 2007). Hemp retting, the process of submersion of hemp stems to facilitate the extraction of its fibers, was tracked in Lake Aydat sediments (France) using cannabinol (Lavrieux et al., 2013; Figure 3). This biomarker allowed the palynological challenge of discriminating hemp cultivation from hemp retting to be overcome. Cannabinol allows not only tracing a type of cultivation, but also the post-harvest processing of cultivated material by humans. Strong environmental after effects are inherent to its use: hemp retting leads to dramatic alteration of water quality (Lavrieux et al., 2013). Hemp retting likely influenced some bacterial communities dynamics (Miras et al., 2015).

The domestication of animals has long been a determinant factor in the modification of landscapes (Diamond, 2002). Spores of the dung fungus Sporormiella recovered from natural archives have been used to trace the extinction of the North American megafauna during the Late Pleistocene, as well as the introduction of livestock (i.e., modern grazing herbivores) during the historic period (Davis and Shafer, 2006; Raper and Bush, 2009). The application of Sporormiella as a proxy for herbivore presence and abundance is however complicated by its spatially sensitive representation: the spores decline sharply with increasing distance from the dung source (Raper and Bush, 2009). Sporormiella spores are released close to the ground and have a gelatinous sheath that facilitates spore attachment on nearby vegetation, thus hampering their dispersion.

Fecal markers such as stanols and bile acid have been used as biomarker indicators for the detection of grazing in a catchment (e.g., D'Anjou et al., 2012). 5β-stanols are organic compounds produced by the microbially mediated alteration of cholesterol in the intestinal tracts of most mammals, making them ideal fecal biomarkers. Bull et al. (2002) proposed the following ratio to distinguish between omnivore (pig/human) and ruminant species: (coprostanol + epicoprostanol/5β-stigmastanol + epi-5β-stigmastanol), where coprostanol (5β-cholestan-3β-ol) is mainly produced by humans and 5β-stigmastanol serves as the biomarker of herbivore fecal matter. Fecal markers have the advantage of dispersing more widely than spores of Sporomiella.

Bile acid products of gut flora can be used to distinguish human contributions from those of pigs: hyodeoxycholic acid is a biomarker of pig feces, whereas human feces are dominated by deoxycholic acid and lithocholic acid (Simpson et al., 1999a).

Guillemot et al. (2015), for instance, reconstructed the evolution of livestock in southwest Greenland over the last two millennia using the concentration of deoxycholic acid in a sedimentary sequence retrieved from Lake Igaliku (Figure 4). The flux of deoxycholic acid correlated quantitatively with that of coprophilous fungal spores, with maximum flux recorded during the two periods of human settlement and grazing activity in the region (i.e., the Norse settlement during the Middle Ages and the recent Danish agricultural phase since 1920).

Another paleolimnological application of fecal sterols comes from the study by D'Anjou et al. (2012), who used these molecular markers as a record of prehistoric human and livestock presence in the sedimentary archive of Lake Liland (Lilandsvatnet, Norway). The arrival of humans in northern Norway is indicated by an abrupt increase in coprostanol (and its epimer epicoprostanol) and an associated increase in 5β-stigmastanol (and 5β-epistigmastanol), which resulted from human and animal feces washing into the lake. Human settlement was accompanied by an abrupt increase in landscape fires (indicated by the rise in pyrolytic PAHs) and a decline in woodland (registered by a change in n-alkane chain lengths).

Recently, Giguet-Covex et al. (2014) used an approach based on DNA metabarcoding of lake sediments to provide a high-resolution reconstruction of plant cover and livestock farming in the catchment of Lake Anterne (French Alps). The presence of cattle and sheep DNA in the lake sediments indicates intense pasturing activity with mixed herds in the catchment during the Late Iron Age and Roman Period, as well as during the second half of the Middle Ages. By comparing these data with a previous reconstruction of erosive event frequency (Giguet-Covex et al., 2011), the authors showed that the most intense erosion period was caused by deforestation and overgrazing by sheep and cowherds during the Late Iron Age and Roman Period.

The ideal biomarker has one single specific source, for instance cannabinol (Section Cultivation). This desire, however, is rarely met. Most organic compounds retrieved from natural archives have multiple sources, some more specific than others (e.g., PTMEs for Graminae vs. leaf waxes for any terrestrial plant). It is thus assumed that a biomarker is specific to a given taxa as long as there is no counter-example. The more specific, the more useful a biomarker is for tracing human-related processes. More detailed studies of the organic compounds in the contemporary biosphere will further “ground-truth” the biomarker approaches.

The most problematic pitfall regarding source material is when there are additional and/or unknown sources of the biomarker of interest, as is the case for branched GDGTs or TTHC (see Section Soil erosion). Changes in the production or proportion of additional sources can influence the preserved record and thus the final interpretation, but these can often be corrected. For instance, for n-alkanes, the potential contribution from bacterial/algal inputs or contamination with petroleum can be assessed using the relative abundance of odd-over-even n-alkanes or the relative abundance of higher molecular-weight components (see Section Vegetation shifts, Pancost and Boot, 2004).

Absolute concentrations in biomarkers are rarely discussed in terms of biomass of its source organisms since they depend on production, transport, sedimentation, and post-deposition processes. Instead, changes in the relative intensity of the signal, or the presence/absence of a biomarker are preferably interpreted and discussed.

The first issue related with quantification is the extent to which environmental factors (e.g., sunlight for plants, fire temperature for fire biomarkers) influence the production of the biomarker over time. The common assumption is that this effect is negligible in comparison to changes in the distribution of the source organism.

The transport mechanisms of biomarkers can also affect its representativeness and should always be assessed in the depositional setting under investigation. Both rivers and wind transport large amounts of terrestrial organic material to the lakes and oceans repositories. The proportion of terrestrial biomarker derived from aeolian as opposed to fluvial processes is geographically variable and variable over time at any given site. In addition, the mode of transport depends on the biomarker source: fire biomarkers will essentially be transported by wind, whereas soil and plant biomarkers dominantly experience fluvial transport. These modes of transport have obviously an impact on the deposition area of the biomarker. The fact that transport processes and catchment areas of biomarkers and microfossils differ can prove very useful in multi-proxy reconstructions. Indeed, when taken into account, the different types of indicators can provide complementary information. For instance, macroscopic charcoal is usually transported only a few 100 m to a few kilometers and is thus mostly a local fire proxy (Whitlock and Larsen, 2001), whereas levoglucosan can be transported as an aerosol hundreds of kilometers and is therefore a regional to hemispheric fire proxy (Zennaro et al., 2014). The opposite is true for pollen and plant biomarker, with pollen having a larger catchment area than biomarker (e.g., Jansen et al., 2013; Garel et al., 2014).

Transport time can also be a source of uncertainties in the biomarker record. For instance, plant waxes have been shown to be transported from continent to ocean in days to weeks (Conte and Weber, 2002), but molecular radiocarbon ages have revealed that some compounds can reside for thousands of years in soils (Eglinton et al., 2002). There is also radiocarbon evidence for long biomarker transit times during fluvial transport (Drenzek, 2007). These residence times in intermediate reservoirs such as soils, river, and littoral sediments can result in lags of the biomarker record compared to more rapidly deposited coarser material, but can be assessed by molecular radiocarbon analyses.

In conclusion, the terrestrial biomarkers transferred from land have to be considered in terms of both the catchment areas that supply them and the selectivity and rates of transfer and deposition by the transport processes. For more detailed reviews of the controls on biomarker transport, we refer the reader to Pancost and Boot (2004) and Eglinton and Eglinton (2008).

Degradation could alter or attenuate the molecular stratigraphic record. However, biomarkers have been shown to be highly persistent compounds. The strong covalent bonds, C–C and C–H, which characterize most of the compounds used as biomarkers, are very difficult to break without high energy input such as heating at high temperatures, radical attack, or through the highly specific activity of enzymes and other catalysts (Eglinton and Eglinton, 2008). Most biomarkers thus survive passage through food web processes (Volkman et al., 1980), consolidation, and the effects of diagenesis (Brassell and Eglinton, 1981). Thus, this type of molecule tends to pass intact through the environment and is conserved in sediments over geological time. We use the following case study from Lake le Bourget (France) to further discuss questions related to the representativeness of molecular biomarkers of human activities, and thus explore the extent to which such a record can be interpreted.

Since the data published by Jacob et al. (2009), the resolution of the miliacin record in core LDB04 has been considerably improved for the Bronze Age (Figure 5). In addition, it has been replicated by another miliacin record on core LDB01, drilled 4 km away in the same lake (Figure 1A). The first striking feature is the comparable evolution of miliacin concentrations over the Bronze Age in the two cores. This implies that miliacin was homogeneously distributed in sediments over the lake (at least at 4 km of distance). Few changes in detritism over the interval exclude varying dilution by the mineral matrix (see Arnaud et al., 2005; Debret et al., 2010). Since miliacin is considered, as with most pentacyclic triterpenes, to be an antifungal and antibacterial agent, changes in early diagenesis can also be discarded. The main control on miliacin concentration in sediments must thus be related to incoming fluxes of this biomarker.

The second interesting feature is the marked transition in miliacin concentrations. For example, the drastic decrease in miliacin concentrations noted around 800 BC is rapid and synchronous in both cores. It coincides with the abandonment of Bronze Age villages around the lake (the youngest dendrochronological date on a wooden pile corresponds to 805 BC; Billaud, 2012), associated with a major lake level rise (Magny, 2004). This implies a limited time lag between villages' abandonment and the decrease in miliacin concentration. This also suggests that there is an immediate (considering the timescale) transport of miliacin from the site of cultivation to the sediment, with limited residence time in soils. This observation is also valid for other marked transitions at 1250, 1200, 1100, 1150, and 850 BC. Our running hypothesis is that most of the miliacin produced in the catchment was directly and quantitatively transported to the sediment. Bossard et al. (2013) showed that miliacin is essentially concentrated in millet seeds. Hence, miliacin transport is related to millet seed transport. In this case, ethnography provides very pertinent information. Poissonnier (1995) reported that millet harvesting was not performed with a sickle but with a dedicated knife, panicle by panicle, to avoid any loss. Therefore, one could expect limited fall of millet seeds, and thus miliacin, during harvesting. Miliacin was therefore directly transported to the villages, sites of millet preparation and consumption. Since villages were located on the shores, miliacin was then directly released to the lake and dispersed above the sediment.

Except for the 850–800 BC interval (with values higher than 350 ng/g), most miliacin concentration values oscillate between ca. 75 and 150 ng/g. As mentioned above, absolute concentrations in biomarkers are rarely discussed. But one could question the origin of such recurrent doublings of miliacin concentrations in Lake le Bourget sediments. The period of maximum miliacin concentration coincides with the maximum density of dendrochronological dates measured on wooden piles (Figure 5) and is considered as the period of maximum population density around the lake (Billaud, 2012). It is thus tempting to relate miliacin concentrations with demography, as millet was staple food for Late Bronze Age populations (the “Golden Age” of millet; Marinval, 1995). However, it is not known whether miliacin production by the plant changes with climate conditions. Varying miliacin concentrations would thus attest to varying millet production due to varying needs from populations of varying size. We currently have no information on whether millet was solely used for human or also cattle consumption and no information on the proportion of millet in diet. Thus, these two explanations for varying miliacin cannot be discarded. However, the pace of climate variability, as depicted by high lake levels (Magny, 2004) and, more distant, sea surface temperatures in Northern Iceland (that proved connected to climate evolution in Europe; Sicre et al., 2008a,b) shows puzzling similarities with that of the miliacin record. Periods of high temperatures coincide roughly with periods of high miliacin concentrations. A causal link between climate variations affecting demography and millet consumption is highly attractive but still speculative and should be confirmed. In summary, the correlation of large shifts in concentration with village abandonment, lake level variations, and overall climate indicates that miliacin concentration could potentially track demographic changes possibly under climate control, but this needs to be confirmed.

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.