How much organic carbon have UK lakes stored in the Holocene? A preliminary estimate

Abstract

Temperate lake sediments store a substantial amount of organic carbon (OC) over millennia. Despite the importance of quantifying terrestrial carbon budgets for Nature-based Solutions, the long-term accumulation of OC in European temperate lakes is poorly constrained. In this study, we analyzed 30 lake sediment records to generate a preliminary first-order estimate of Holocene OC accumulation rate (OCAR) and OC storage in UK lakes. We also examined the environmental variables that influence OCAR and produced synthesized Holocene records of %OC and z-scores of log-transformed OCAR and sediment accumulation rate (SAR) at 500-year resolution. Based on our estimation, we report an average Holocene OCAR of 7.4 ± 5.5 g C m⁻² yr⁻¹ and a Holocene total OC storage of 0.24 ± 0.18 Pg C in UK lakes. Apart from latitude, no relationship was found between the average Holocene OCAR and the various environmental variables (i.e. temperature, precipitation, surface area, catchment area, depth, altitude, and geology type). During the Holocene, OCAR closely resembles variations in SAR, whereas the increase in %OC is likely explained by the warming climate. Early Holocene variations in OCAR were primarily climate-driven. In contrast, the anthropogenic impact on the landscape exerted a predominant influence on OC burial during the middle-late Holocene. Our results improve the current understanding of terrestrial carbon budgets in the UK and demonstrate the under-appreciated importance of lakes as long-term OC stores.

Keywords

Holocene Lake sediments Organic carbon burial Paleolimnology Temperate lakes Terrestrial carbon budgets UK

Introduction

In the context of anthropogenic climate change, there is an unparalleled need to achieve carbon neutrality by this century. Nature-based Solutions (NbS) comprise a crucial component of climate change mitigation due to the large amount of carbon stored in natural and managed ecosystems. NbS has the potential to mitigate CO₂ emission at a rate of 10 Gt CO₂ yr⁻¹ and, in the long-term, reduce peak warming by 0.1–0.3°C depending on the different warming scenarios (Girardin et al., 2021). Quantifying terrestrial carbon budgets is conducive to implementing NbS and could help guide more targeted and informed land management decisions (Cole et al., 2007; Gregg et al., 2021; Keith et al., 2021).

Despite only constituting 3.7% of the Earth’s non-glaciated land surface area (Verpoorter et al., 2014), lakes are highly spatially variable and play a critical role in the global carbon cycle (Buffam et al., 2011; Einsele, 2001; Prairie, 2008; Radbourne et al., 2017; Tranvik et al., 2009). Lake ecosystems are sites where active processing of carbon (e.g. carbon production, transformation, burial, and evasion) takes place (Cole et al., 2007; Prairie, 2008; Tranvik et al., 2009). Allochthonous input of organic matter from catchments contributes to lake’s dissolved organic carbon (OC) or dissolved organic matter pool, which can be assimilated by heterotrophs as an energy source (Hanson et al., 2015; Prairie, 2008). Dissolved inorganic carbon, mainly in the form of CO₂, is consumed during photosynthetic assimilation by heterotrophs and oxidized into organic carbon in both dissolved and particulate forms (McGowan et al., 2016). OC may be reduced back into CO₂ and CH₄ through respiration. Particulate carbon suspended in the water column may sink to the lake bottom and contribute to OC burial in lake sediments (Alin and Johnson, 2007).

Although lakes mostly behave as net carbon sources and are super-saturated in CO₂ in the surface waters (Cole et al., 1994; Tranvik et al., 2009), lake sediments can store a large amount of OC in the long-term due to their high OC burial efficiency, rapid sediment accumulation rates and high aquatic productivity (Cole et al., 2007; Sobek et al., 2009). Some of the deeper and larger lakes of tectonic origin could persist in the landscape for millions of years, and lake sediments in glaciated terrains could accumulate for more than 10,000 years (Alin and Johnson, 2007; Cole, 2013). In contrast, many smaller lakes are short-lived or artificially constructed (i.e. reservoirs). Globally, inland waters (including lakes, rivers, streams, and reservoirs) receive carbon at the magnitude estimated at 2.9 Pg C yr⁻¹, amongst which 1.4 Pg C yr⁻¹ are evaded into the atmosphere, 0.9 Pg C yr⁻¹ are transported to the sea, and the remaining 0.6 Pg C yr⁻¹ are buried to sediments (Tranvik et al., 2009). The total amount of carbon storage in lake sediments during the Holocene is estimated at 400–800 Pg C (Cole et al., 2007), and the long-term OC accumulation rate (OCAR) ranges between 4.5 and 14 g C m⁻² yr⁻¹ (Cole et al., 2007; Dean and Gorham, 1998; Stallard, 1998). Recently, OCAR in lake sediments has increased due to cultural eutrophication (Anderson et al., 2014, 2020). Carbon burial is expected to continue to increase, driven by climate change, increased runoff and eutrophication (Tranvik et al., 2009).

Lake sediments document spatially and temporally integrated information. Hence, paleolimnological approaches are well-suited for quantifying OC accumulation and storage in lakes and may offer critical insights for understanding future changes in the terrestrial carbon cycle (McGowan et al., 2016). One widely adopted method uses loss-on-ignition 550°C (%LOI₅₅₀) data to estimate %OC in lake sediments (Anderson et al., 2014; Pajunen, 2000) or directly makes use of downcore total organic carbon (TOC) data measured using an elemental analyzer (McGowan et al., 2016). To date, many studies have adopted a paleolimnological approach for quantifying the amount of OC storage and OCAR in lakes at different spatial and temporal scales (e.g. Anderson et al., 2020; Dean and Gorham, 1998; Kastowski et al., 2011). However, there is a bias toward Holocene scale estimates on boreal and arctic lakes within Europe (e.g. Anderson et al., 2009; Chmiel et al., 2015; Kortelainen et al., 2004), whereas temperate lakes, which are large organic matter stores, receive less attention (Scott, 2014). Temperate lakes cover ~25% of the European lake surface area but contribute to ~35% of European lakes’ total carbon accumulation rate (Scott, 2014).

The UK is home to 43,738 water bodies covering 213,911 ha in Great Britain (Hughes et al., 2004) and 1,668 water bodies covering 62,600 ha in Northern Ireland (Gibson et al., 1994), altogether constituting around 1% of the total UK surface area. There are few systematic inventories of OCAR in the UK. Scott (2014) examined the OCAR of two lakes in the Shropshire–Cheshire meres region, and Anderson et al. (2014) included a number of UK lakes in their OCAR estimation for European lakes. Both studies focused on the last 100–150 years. Only a few UK lakes were included in the estimation of Holocene OCAR in European lakes (Kastowski et al., 2011). Thus, there is no clear overview of Holocene OCAR and OC storage estimates in the UK. There are extensive lake sediment records in the UK covering part or the complete sequence of the Holocene (e.g. Bennett et al., 1992; Fossitt, 1996; Watkins et al., 2007), offering an excellent opportunity to conduct a Holocene OC inventory.

Here, we present a first-order preliminary estimate for OCAR and OC storage in UK lakes across the Holocene by adopting a paleolimnological approach. In addition, we investigate the temporal variations in OCAR, sediment accumulation rate (SAR) and %OC during the Holocene and examine the environmental variables influencing OCAR in UK lakes (i.e. temperature, precipitation, surface area, catchment area, depth, altitude, latitude, and geology type). By comparing the OCAR in UK lakes to estimates from lakes in different regions and various habitat types in the UK, we evaluate the relative importance of lake sediments as an OC store in the UK.

Methods

Study sites and data collection

The selection of lake sediment cores was based on four main criteria: (1) the sediment record covers at least 3 kyr during the Holocene; (2) a %LOI₅₅₀, TOC, or organic matter (OM) record exists; (3) a minimum of three radiocarbon dates to ensure a relatively reliable age control and (4) extant lakes.

A wide range of search terms was employed to identify relevant studies. These include “lake,” “loch,” “llyn,” “lough,” “Holocene,” “Britain,” “England,” “Scotland,” “Wales,” “Northern Ireland,” “organic matter,” “LOI,” “TOC,” “vegetation,” and “pollen.” We are aware that this search was not exhaustive, and additional sediment records must exist in earlier publications or remain unpublished. However, in contrast to other Holocene estimates for lakes (e.g. Anderson et al., 2009, n = 11; Chmiel et al., 2015, n = 7), fjords (e.g. Smeaton et al., 2017, n = 5), and peatlands (e.g. Ratcliffe et al., 2018, n = 12), the number of lakes included is sufficient for a preliminary first-order estimate.

A total of 30 lake sediment records across the UK were selected (Figure 1). This includes three new lake sediment records (Loweswater, Loch na Claise, and Loch an Aigeil) that were analyzed for %LOI₅₅₀ following Heiri et al. (2001), which entails burning at 550°C for 4 hours and calculating where DW₁₀₅ refers to the dry weight of sediments after drying at 105°C. %LOI₅₅₀ data for the 27 published records were generated using a temperature range between 450°C and 550°C and furnace times between 1 and 12 hours (references are in the Supplementary Information, available online). Scottish lakes (n = 19) comprise most of the selected records, followed by lakes in England (n = 6), Wales (n = 3), and Northern Ireland (n = 2). The proportional representation of lakes from each country is roughly consistent with the spatial distribution of lakes in the UK, dominated by lakes in Scotland (57%) and England (36%) in terms of absolute numbers (Gibson et al., 1994; Hughes et al., 2004).

% L O I_{550} = 100 (\frac{(D W_{105} - D W_{550})}{D W_{105}})

(1)

Figure 1.

Overview map showing the 30 study sites. Green triangles (n = 11) represent lake sediment records with reported dry bulk density data (required to calculate organic carbon accumulation rate); red circles (n = 19) represent lake sediment records lacking dry bulk density data. List of lake names, published reference sources, and further information associated with each lake are given in the Supplementary Information, available online.

Chronology

To obtain Holocene temporal trends in OCAR, SAR, and %OC, radiocarbon dates for all 30 lake records were (re)calibrated with the IntCal20 curve (Reimer et al., 2020) using Bayesian age-depth modeling in the rbacon package (version 2.5.8) (Blaauw and Christen, 2011) in R (see Supplemental Figures S1–S30, available online for the new age models and Supplemental Tables S1–S3, available online for radiocarbon ages from the three new records). The 95% confidence interval from the Bacon age-depth models were used as conservative age uncertainties. The maximum age uncertainties across all interpolated depth layers ranged from 289 (Lochan Uaine) to 1,705 years (Loch na Claise), with 77% of the lake records having maximum age uncertainties below 1,000 years.

Calculation of average Holocene OCAR and OC storage

%LOI₅₅₀ was multiplied by a factor of 0.469 to convert it into %OC based on previous comparisons of %LOI₅₅₀ and %OC data (Dean, 1974). Most of the published literature (n = 19) did not report data on dry bulk density (DBD) (Figure 1), but DBD and SAR are required to calculate dry mass accumulation rate (DMAR), which is essential in the calculation of OCAR (see below). Following Dean and Gorham (1998), an empirical relationship between %LOI₅₅₀ and DBD was established using data from the 11 lakes with DBD reported (Figure 1) (equation (2)), which has r2 = 0.6 and p < 2.2 × 10⁻¹⁶. This %LOI₅₅₀-DBD relationship is highly similar to existing empirical relationships (Dean and Gorham, 1998; Moyle et al., 2021) (see Supplemental Figure S31, available online).

D B D = 2.4355 % L O I_{550}^{- 0.827}

(2)

This relationship was applied to the estimation of DBD for the 19 lakes. OCAR (g C m⁻² yr⁻¹) was determined by multiplying DMAR (g cm⁻² yr⁻¹) by %OC and a conversion factor of 10,000 (Anderson et al., 2014).

To account for the higher-than-average SAR occurring at the deepest part of the lake (Anderson et al., 2014; Engstrom and Rose, 2013), a sediment focusing factor based on the ratio of sediment ²¹⁰Pb flux to regional atmospheric ²¹⁰Pb flux was applied. Due to the paucity of sediment ²¹⁰Pb flux data for most study sites, we employ a regional average focusing factor of 1.78 based on data from 90 UK lakes since 1850 CE (Handong Yang, unpublished data). Sediment focusing varies over time as the pattern of deposition evolves and lake morphometry or size changes (Engstrom and Rose, 2013). However, we assume that the focusing factor remains constant over the Holocene for ease of estimation. The focusing-corrected average Holocene OCAR in each lake sediment record was then averaged to estimate the average focusing-corrected Holocene OCAR in UK lakes. The total Holocene OC storage in UK lakes was calculated by multiplying the average focusing-corrected Holocene OCAR by 11,650 years and a total UK lake area of 276,511 ha (Gibson et al., 1994; Hughes et al., 2004). Standard deviations were calculated for the OCAR in each sediment record and propagated to obtain uncertainty of the average Holocene OCAR in UK lakes. We recognize that there are compounding error sources associated with each step of the estimation, but the reported error values represent a conservative estimate of the uncertainty, which means that the true value lies within the range provided.

Examining controlling factors on OCAR

Pearson’s product-moment correlation was performed to determine the relationship between the average Holocene OCAR from each sediment record and various modern environmental variables at each site, assuming that OCAR and these variables have remained constant over time. Data for surface area, catchment area, altitude, mean depth, maximum depth, geology type (categorized into low, moderate and high alkalinity groups) and latitude were obtained from the UK Lakes Portal (UK Centre for Hydrology and Ecology, 2016) and the original publications; the maximum annual temperature and mean annual precipitation between 1961 and 2020 at the nearest weather station for each site were obtained from the Met Office (2022). Spearman’s rank correlation coefficient was performed to analyze the relationship between OCAR and the three geology types.

Average UK Holocene temporal trends

The OCAR, SAR and %OC data for each lake were linearly interpolated between neighboring data points to the nearest 500-year since the sampling resolution of all sediment cores lies below 500-year. The OCAR and SAR data were log-transformed for normality (hereafter referred to as OCAR_log and SAR_log). For each sediment record, z-scores of OCAR_log and SAR_log were calculated by subtracting the mean and dividing by the standard deviation. Synthesized OCAR_log and SAR_log z-scores were produced by taking the average z-scores for all lakes at each resampled age (i.e. one average value every 500 years). Pearson’s product-moment correlation was performed to assess the relationship between the synthesized OCAR_log z-score and the synthesized SAR_log z-score and %OC.

Results

Average Holocene OCAR and OC storage estimates

The average focusing-corrected Holocene OCAR in UK lakes is 7.4 ± 5.5 g C m⁻² yr⁻¹ (uncorrected OCAR: 13.2 ± 9.9 g C m⁻² yr⁻¹), with values ranging from 2.5 ± 0.9 g C m⁻² yr⁻¹ for Loch Buailaval Beag to 20.3 ± 8.9 g C m⁻² yr⁻¹ for Lough Catherine (Figure 2). The area-weighted OCAR is 21 ± 15 Gg C yr⁻¹ (rounded to two significant figures). During the Holocene, the total OC storage in UK lakes is approximately 0.24 ± 0.18 Pg C with areal OC storage of 86.5 ± 64.6 kg C m⁻².

Figure 2.

Focusing-corrected average Holocene organic carbon accumulation rate (OCAR) (g C m⁻² yr⁻¹) in each lake sediment record.

Relationship between OCAR and environmental variables

The only variable that significantly correlates with the average Holocene OCAR from each sediment record is latitude (r = −0.5, p = 3.2 × 10⁻³). No significant correlations were found between OCAR and any other of the variables tested (see Supplemental Information, available online).

Average UK Holocene temporal trends

Since 12 cal kyr BP, %OC has increased rapidly from ca. 6.8% to 17.3%, stabilizing after 8 cal kyr BP around 16–17% with a gradual increase (Figure 3a) to present-day values of 18.1%. The OCAR_log z-score remains stable between 12 and 10 cal kyr BP and subsequently declines between 10 and 8 cal kyr BP (Figure 3b and 3d). In contrast, the SAR_log z-score displays a more extended decline between 12 and 8 cal kyr BP (Figure 3c). Both records show a long-term increase from 8 cal kyr BP to the present, interrupted by a temporary decline between 6 and 5 cal kyr BP. After 3 cal kyr BP, OCAR_log z-score lies predominantly above the average (i.e. when z-score = 0), which is greater than the early and middle Holocene.

Figure 3.

Average Holocene temporal trends for 30 UK lakes with LOESS smoothing (95% confidence interval, span = 0.4) at 500-year resolution: (a) percentage organic carbon (%OC); (b) raw log-transformed organic carbon accumulation rate (OCAR_log) z-score showing data points from all sediment records; (c) synthesized log-transformed sediment accumulation rate (SAR_log) z-score; (d) synthesized OCAR_log z-score.

A very strong correlation exists between the OCAR_log and SAR_log z-scores (r = 0.9, p = 1.8 × 10⁻⁸). An r² of 0.8 suggests that SAR_log can explain up to 80% of the variation in OCAR_log. On the contrary, %OC displays no relationship with the OCAR_log z-score (r = −0.2, p = 0.5).

Discussion

Average Holocene OCAR and OC storage estimates

The total Holocene OC storage for UK lakes is estimated to be 0.24 ± 0.18 Pg C (areal OC storage: 86.5 ± 64.6 kg C m⁻²), equivalent to 0.88 ± 0.66 Pg CO₂. The magnitude of the average Holocene OCAR of UK lakes at 7.4 ± 5.5 g C m⁻² yr⁻¹ is in good agreement with the global range for lakes between 4.5 and 14 g C m⁻² yr⁻¹ (Cole et al., 2007; Dean and Gorham, 1998; Stallard, 1998), and is higher than the average Holocene OCAR estimates for lakes in Europe (5.2 g C m⁻² yr⁻¹) (Kastowski et al., 2011), southwest Greenland (6 g C m⁻² yr⁻¹) (Anderson et al., 2009), Finland (1.6 g C m⁻² yr⁻¹) (Kortelainen et al., 2004), and Sweden (3.2 g C m⁻² yr⁻¹) (Chmiel et al., 2015). Variations in OCAR in Europe are likely due to region-specific lake characteristics. For example, lower temperatures in Finland may suppress primary production, thus leading to lower OCAR in boreal lakes (Kortelainen et al., 2004). In contrast, the higher OCAR in southwest Greenland may be explained by the selection of very small lakes in the study (range: 0.05–0.73 km²) (Anderson et al., 2009) compared to the selection of larger Finnish lakes (average area: 122 km²) (Kortelainen et al., 2004) because OCAR is often found to be inversely related with lake area (Mulholland and Elwood, 1982; Pajunen, 2000).

Within habitats of the UK, the Holocene OCAR of lakes (7.4 ± 5.5 g C m⁻² yr⁻¹) is higher than the modern estimates for broadleaf and mixed woodlands (4.2 g C m⁻² yr⁻¹, range: 0.7–12 g C m⁻² yr⁻¹), coniferous woodlands (4.9 g C m⁻² yr⁻¹, range: 0.8–11.7 g C m⁻² yr⁻¹) and grasslands (2.2 g C m⁻² yr⁻¹) (Downing et al., 2008), but lower than the Holocene estimates for fjords (8.3 ± 2.0 g C m⁻² yr⁻¹) (Smeaton et al., 2017), peatlands (17.8 ± 4.9 g C m⁻² yr⁻¹) (Ratcliffe et al., 2018), the modern estimate for ponds aged between 18 and 20 years (142 ± 19 g C m⁻² yr⁻¹) (Taylor et al., 2019) and millennial-scale estimates for seagrasses (138 ± 38 g C m⁻² yr⁻¹) and salt marshes (218 ± 24 g C m⁻² yr⁻¹) (Mcleod et al., 2011) (Figure 4a). The area-weighted OCAR for UK lakes (21 ± 15 Gg C yr⁻¹) is lower than all UK habitats apart from seagrasses (Figure 4b) due to the small total surface area of lakes in the UK compared with other habitats.

Figure 4.

Comparison of (a) average Holocene organic carbon accumulation rate (OCAR) and (b) area-weighted Holocene OCAR in a range of habitats in the UK. Sources of OCAR data include Mcleod et al. (2011) for seagrasses and salt marshes, Smeaton et al. (2017) for fjords, Taylor et al. (2019) for ponds, Downing et al. (2008) for woodlands and grasslands and Ratcliffe et al. (2018) for peatlands. Data for habitat area used for calculating area-weighted Holocene OCAR was obtained from Gibson et al. (1994) and Hughes et al. (2004) for lakes, Taylor et al. (2019) for ponds, Carey et al. (2008) for woodlands, grasslands and peatlands, Smeaton et al. (2017) for fjords, Phelan et al. (2011) for salt marshes and Green et al. (2021) for seagrasses.

In comparison to most terrestrial habitats, vegetated coastal habitats exhibit very high OCARs owing to their high efficiency in trapping sediments from both riverine and oceanic sources, facilitated by their complex ecosystem structure (Mcleod et al., 2011). Peatlands and ponds also exhibit higher OCARs than lakes. Peatlands are highly efficient OC sinks due to active plant growth (Loisel et al., 2014). Ponds are generally very productive and are effective in trapping and preserving carbon in sediments (Cole et al., 2007; Downing et al., 2006). The higher OCAR in lakes compared to woodlands and grasslands is due to the higher OC preservation rate of lakes (Cole et al., 2007).

The OC turnover rates in ponds or forest biomass last for decades to centuries (Cole, 2013). The long-term storage of OC can only occur in soils once forests reach maturity, and yet, the long-term OC burial in soil is much limited compared to its instantaneous carbon assimilation rate due to periodic burning and/or high decay rates (Prairie, 2008). Moreover, despite having very high OCAR, ponds are ephemeral features of the landscape that is more susceptible to in-filling and desiccation compared to lakes. In contrast, lake sediments, peatlands and coastal habitats can store OC for millennia (Cole et al., 2007; Mcleod et al., 2011; Prairie, 2008; Ratcliffe et al., 2018; Taylor et al., 2019). Consequently, although the average Holocene OCAR of lakes is modest in contrast to various habitat types in the UK, the capacity of lakes to act as long-term carbon stores cannot be overlooked.

The estimation of Holocene OC storage and accumulation rate in UK lakes provides an important constraint on the regional terrestrial carbon budget. Using the past as an analog, our results are highly relevant for understanding how carbon budgets in lakes will change in the future (Hanson et al., 2015; McGowan et al., 2016). Moreover, despite the promising potential of NbS for climate change mitigation, there is currently a bias toward carbon flows (e.g. sequestration) and an underappreciation of carbon storage in carbon accounting (Keith et al., 2021). While both processes are essential for the ecosystem service of climate regulation, the overlooking of carbon storage hampers current NbS efforts by providing incomprehensive information about the mitigation benefits, which may result in perverse outcomes in NbS implementation (Keith et al., 2021). Hence, by quantifying the long-term carbon storage in lake sediments in the UK, we provide critical information that may contribute to a more comprehensive carbon accounting of terrestrial aquatic systems beneficial for facilitating better informed NbS practices and associated land management decisions. Despite the long-term mitigation benefits NbS may offer, we stress that these benefits are limited compared to decarbonization efforts (Girardin et al., 2021). Therefore, NbS should not act as a distraction from emission reduction activities in the near term.

Controlling factors on OCAR

Previous studies have identified a wide range of environmental variables controlling OC burial in lakes. The main variables include depth (Anderson et al., 2014), surface area (Downing et al., 2008; Mulholland and Elwood, 1982; Pajunen, 2000), lake morphometry (Ferland et al., 2014), precipitation (Tranvik et al., 2009), temperature (Gudasz et al., 2015; Rantala et al., 2016), nutrients (Anderson et al., 2020), oxygen exposure time, and mineralization (Sobek et al., 2009). Here, we found a moderate relationship between OCAR and latitude, which may suggest an indirect influence of cooling as latitude increases (Kastowski et al., 2011). The absence of a relationship between OCAR and the environmental variables we tested may suggest the relationships found in previous studies may not apply universally. Moreover, none of the variables alone shape variations in OCAR between lakes, but rather, a combination of many site-specific factors. However, it is beyond the scope of this study to explain variations between lakes. Besides, with over 67% of the studied lakes ⩽10 ha, the absence of a relationship between OCAR and the surface area might be due to the bias toward very small lakes with variable OCAR yet narrow environmental gradients (Kastowski et al., 2011).

Average UK Holocene temporal trends

The strong positive correlation between the synthesized OCAR_log and SAR_log z-scores indicates that variations in OCAR_log are largely driven by SAR_log, and %OC of the sediment has not changed markedly during the Holocene. The increase in %OC during the early Holocene (12–8 cal kyr BP) likely resulted from the increasing annual temperature in Europe (Figure 5a and 5b) (Davis et al., 2003; Lang et al., 2010), which facilitated the switch from late glacial minerogenic sediment to organic gyttja as catchment soil developed (Bennett et al., 1990; Edwards and Whittington, 1998; Snowball and Thompson, 2008). The decrease in the synthesized OCAR_log and SAR_log z-scores prior to 8 cal kyr BP can be explained by the progressive establishment of woodlands fostered by the increasingly warm climate, which stabilizes the soils and reduces runoff and erosion (Figure 5a, 5f, and 5g) (Bennett et al., 1992; Fossitt, 1996; Watkins et al., 2007).

Figure 5.

Dominant drivers of Holocene organic carbon accumulation rate (OCAR) and sediment accumulation rate (SAR) variations in 30 UK lakes (see Figure 1). (a) Pollen-reconstructed annual temperature record for Europe (Davis et al., 2003); (b) percentage organic carbon (%OC) at 500-year resolution; (c) radiocarbon (¹⁴C) date probability density functions from archeological sites between 9 and 3.4 cal kyr BP (Woodbridge et al., 2014); (d) pollen-inferred deciduous woodland change since 9 cal kyr BP (Woodbridge et al., 2014); (e) C:N ratios from three study sites (blue: Lochnagar; green: Loch nan Corr; red: Loch na Claise) with LOESS-smoothing (span = 0.2) (Dalton et al., 2005; Mackie et al., 2007; Matthews, 2022); (f and g) synthesized log-transformed SAR (SAR_log) and OCAR (OCAR_log) z-scores at 500-year resolution. The yellow bar above (g) represents the time range for the elm decline between 6,347 and 5,281 cal kyr BP based on statistical analyses of radiocarbon dates (Parker et al., 2002). The period-based chronology for pre-historic Britain was obtained from Pollard (2008). The green shading separates out intervals where climate exerts more influence on OCAR and SAR variations from periods where anthropogenic impact on the landscape dominates (see text for details).

Since 8 cal kyr BP, the long-term increase in the synthesized OCAR_log and SAR_log z-scores suggests increasing anthropogenic disturbance on the landscape (Figure 5f and 5g). The initial population increase in Britain started from ca. 7.6 cal kyr BP during the Late Mesolithic, followed by a long-term population growth trend, indicated by the increase in ¹⁴C date densities from archeological evidence and the sustained decline in deciduous woodlands (Figure 5c and 5d) (Woodbridge et al., 2014). On the contrary, %OC has remained relatively stable since 8 cal kyr BP, resembling the stabilized trend in temperature (Figure 5a and 5b).

The local maxima in the synthesized OCAR_log and SAR_log z-scores between 6.5 and 5.5 cal kyr BP coincides with the widely documented elm decline on the British Isles (6,347–5,281 cal yr BP) (Figure 5f and 5g) for which human activity is one of the hypothesized causes (Parker et al., 2002). The elevated sedimentation and OC burial may have been driven by the transition from Mesolithic hunter-gatherers to Neolithic agriculture in the British Isles between 6.4 and 6 cal kyr BP (Woodbridge et al., 2014) and the substantial population growth in Britain between 6.1 and 5.4 cal kyr BP (Collard et al., 2010) (Figure 5c), both increasing anthropogenic pressure on the landscape. The increase in the synthesized z-scores is interrupted by a decline toward lower values in both records between 5.5 and 4.5 cal kyr BP. This interval coincides with reduced ¹⁴C date density between 5.3 and 4.4 cal kyr, indicating reduced Neolithic impact on the landscape, which leads to the re-establishment of woodlands (Figure 5c, 5d, 5f, and 5g) (Woodbridge et al., 2014).

Since 5–4 cal kyr BP, the continued increase in the synthesized OCAR_log and SAR_log z-scores suggests enhanced soil erosion, which might be attributed to anthropogenic activities (Edwards and Whittington, 2001; Van Vliet-Lanoë et al., 1992). The gradual increase in both z-scores between 4 and 3 cal kyr BP coincides with the British Bronze Age (2,200–800 BC or 4,150–2,750 BP), whilst the accelerated increase since ca. 3–2.5 cal kyr BP coincides with the British Iron Age (800 BC–43 CE/500 CE or 2,750–1,907/1,450 BP) (Pollard, 2008). Changes in OCAR and SAR are consistent with the long-term increase in C:N at Loch na Claise, Loch nan Corr, and Lochnagar, with an accelerated increase over the last 2 kyr (Figure 5e) (Dalton et al., 2005; Mackie et al., 2007; Matthews, 2022). These suggest an increase in allochthonous sediment input linked to human activities. However, there are limited records of C:N from the studied lakes to help further discriminate between allochthonous and autochthonous sources of OC. Compared to natural conditions, anthropogenic sediment delivery can be 5–10 times higher (Dearing and Jones, 2003). Increased land-use pressure, including deforestation, agriculture, and grazing, have all exacerbated soil erosion and increased allochthonous carbon input to lakes (Beales, 1980; Bennett et al., 1992; Edwards and Whittington, 1998; Jones et al., 1985; Mulder, 1999). Soil erosion would also elevate nutrient supply into lakes, increasing autochthonous OC accumulation due to increased aquatic productivity.

Conclusion

This study represents a preliminary first-order estimate of Holocene OCAR and OC storage in UK lakes. We acknowledge a number of limitations, namely, the assumption of a non-varying %LOI₅₅₀-DBD relationship, spatial and temporal extrapolations of OCAR, uncertainties of radiocarbon chronologies and issues of linear interpolation, which could be improved with analyses of new dated sediment cores. Although the OCAR of lakes is modest compared to various coastal and terrestrial habitats in the UK, our results demonstrate the importance of lakes as long-term OC stores. While more lake sediment records could be incorporated into the estimation for better spatiotemporal representation of UK lakes, our estimates based on available data improve the current understanding of terrestrial carbon storage in the UK and OC accumulation in European temperate lakes. Our results may contribute important information for understanding how lake carbon budgets would change in the future, as well as for implementing a more comprehensive carbon accounting for better-informed NbS and land management practices.

Moreover, our results indicate that the first-order response of lake OCAR was to climate amelioration during the early Holocene, initiating organic matter production in terrestrial landscapes and its burial in lake sediments. The second major driver of OCAR since the middle Holocene has been the anthropogenic impact on the landscape, which has resulted in enhanced transport and erosion of OC from the deforested landscapes into lakes and also increased aquatic productivity due to a greater nutrient supply. Further work using geochemical analyses is required to distinguish the type or source of OC that reaches lakes during the middle-late Holocene.

Supplemental Material

sj-pdf-2-hol-10.1177_09596836231157062 – Supplemental material for How much organic carbon have UK lakes stored in the Holocene? A preliminary estimate

Supplemental material, sj-pdf-2-hol-10.1177_09596836231157062 for How much organic carbon have UK lakes stored in the Holocene? A preliminary estimate by Mengyao Du, Vivienne J Jones, Maarten van Hardenbroek, Louisa Matthews and Kejia Ji in The Holocene

Supplemental Material

sj-xlsx-1-hol-10.1177_09596836231157062 – Supplemental material for How much organic carbon have UK lakes stored in the Holocene? A preliminary estimate

Supplemental material, sj-xlsx-1-hol-10.1177_09596836231157062 for How much organic carbon have UK lakes stored in the Holocene? A preliminary estimate by Mengyao Du, Vivienne J Jones, Maarten van Hardenbroek, Louisa Matthews and Kejia Ji in The Holocene

Footnotes

Acknowledgements

We would like to sincerely thank Neil Rose and the two anonymous reviewers for their helpful comments and suggestions. We are also grateful to Handong Yang for providing unpublished ²¹⁰Pb data for calculating the average sediment focusing factor.

Correction (May 2023):

Figures 2 and 4a have been updated to correct the y-axis labels.

Correction (January 2024):

Figure 5 has been updated to correct the x-axis label.

Funding

The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: Radiocarbon chronologies for three new sediment records analyzed in this study were provided by NEIF (NF/2019/1/11 and NEIF 2364.0321), and the PhD studentship for Louisa Matthews was provided by NERC-UKRI (IAP-17-91).

Data availability

The OCAR dataset is available in supplementary materials.

ORCID iD

Mengyao Du

Supplemental material

Supplemental material for this article is available online.

References

Alin

Johnson

(2007) Carbon cycling in large lakes of the world: A synthesis of production, burial, and lake-atmosphere exchange estimates. Global Biogeochemical Cycles 21(3): GB3002.

Anderson

Bennion

Lotter

(2014) Lake eutrophication and its implications for organic carbon sequestration in Europe. Global Change Biology 20(9): 2741–2751.

Anderson

D’Andrea

Fritz

(2009) Holocene carbon burial by lakes in SW Greenland. Global Change Biology 15(11): 2590–2598.

Anderson

Heathcote

Engstrom

(2020) Anthropogenic alteration of nutrient supply increases the global freshwater carbon sink. Science Advances 6(16): eaaw2145.

Beales

(1980) The late Devensian and Flandrian vegetational history of Crose Mere, Shropshire. New Phytologist 85(1): 133–161.

Bennett

Boreham

Sharp

, et al. (1992) Holocene history of environment, vegetation and human settlement on Catta ness, lunnasting, Shetland. Journal of Ecology 80(2): 241–273.

Bennett

Fossitt

Sharp

, et al. (1990) Holocene vegetational and environmental history at Loch Lang, South Uist, Western Isles, Scotland. New Phytologist 114(2): 281–298.

Blaauw

Christen

(2011) Flexible paleoclimate age-depth models using an autoregressive gamma process. Bayesian Analysis 6(3): 457–474.

Buffam

ISHI

Turner

Desai

, et al. (2011) Integrating aquatic and terrestrial components to construct a complete carbon budget for a north temperate lake district: Carbon budget in lake-rich region. Global Change Biology 17(2): 1193–1211.

10.

Carey

Wallis

Emmett

, et al. (2008) Countryside survey: UK results from 2007. Report for the Natural Environment Research Council, 10 December. Swindon: Natural Environment Research Council.

11.

Chmiel

Niggemann

Kokic

, et al. (2015) Uncoupled organic matter burial and quality in boreal lake sediments over the holocene. Journal of Geophysical Research Biogeosciences 120(9): 1751–1763.

12.

Cole

(2013) Freshwater in flux. Nature Geoscience 6(1): 13–14.

13.

Cole

Caraco

Kling

, et al. (1994) Carbon dioxide supersaturation in the surface waters of lakes. Science 265(5178): 1568–1570.

14.

Cole

Prairie

Caraco

, et al. (2007) Plumbing the global carbon cycle: Integrating inland waters into the terrestrial carbon budget. Ecosystems 10(1): 172–185.

15.

Collard

Edinborough

Shennan

, et al. (2010) Radiocarbon evidence indicates that migrants introduced farming to Britain. Journal of Archaeological Science 37(4): 866–870.

16.

Dalton

Birks

HJB

Brooks

, et al. (2005) A multi-proxy study of lake-development in response to catchment changes during the holocene at Lochnagar, north-east Scotland. Palaeogeography Palaeoclimatology Palaeoecology 221(3–4): 175–201.

17.

Davis

BAS

Brewer

Stevenson

, et al. (2003) The temperature of Europe during the holocene reconstructed from pollen data. Quaternary Science Reviews 22(15–17): 1701–1716.

18.

Dean

Jr (1974) Determination of carbonate and organic matter in calcareous sediments and sedimentary rocks by loss on ignition: Comparison with other methods. SEPM Journal of Sedimentary Research 44(1): 242–248.

19.

Dean

Gorham

(1998) Magnitude and significance of carbon burial in lakes, reservoirs, and peatlands. Geology 26(6): 535–538.

20.

Dearing

Jones

(2003) Coupling temporal and spatial dimensions of global sediment flux through lake and marine sediment records. Global and Planetary Change 39(1–2): 147–168.

21.

Downing

Cole

Middelburg

, et al. (2008) Sediment organic carbon burial in agriculturally eutrophic impoundments over the last century: Sediment organic carbon burial. Global Biogeochemical Cycles 22(1): GB1018.

22.

Downing

Prairie

Cole

, et al. (2006) The global abundance and size distribution of lakes, ponds, and impoundments. Limnology and Oceanography 51(5): 2388–2397.

23.

Edwards

Whittington

(1998) The palaeoenvironmental background: Pollen studies at Rae Loch. In: Barclay

Maxwell

(eds) The Cleaven Dyke and Littleour: Monuments in the Neolithic of Tayside, vol. 13. Edinburgh: Society Antiquaries Scotland Monograph, pp.5–12.

24.

Edwards

Whittington

(2001) Lake sediments, erosion and landscape change during the holocene in Britain and Ireland. CATENA 42(2–4): 143–173.

25.

Einsele

(2001) Atmospheric carbon burial in modern lake basins and its significance for the global carbon budget. Global and Planetary Change 30(3–4): 167–195.

26.

Engstrom

Rose

(2013) A whole-basin, mass-balance approach to paleolimnology. Journal of Paleolimnology 49(3): 333–347.

27.

Ferland

Prairie

Teodoru

, et al. (2014) Linking organic carbon sedimentation, burial efficiency, and long-term accumulation in boreal lakes. Journal of Geophysical Research Biogeosciences 119(5): 836–847.

28.

Fossitt

(1996) Late quaternary vegetation history of the Western Isles of Scotland. New Phytologist 132(1): 171–196.

29.

Gibson

Smith

Wolfe-Murphy

(1994) The Northern Ireland lake survey. In: Brian

(ed.) Irish Botanical News. Belfast: The Committee for Ireland Botanical Society of the British Isles, pp.5–11.

30.

Girardin

CAJ

Jenkins

Seddon

, et al. (2021) Nature-based solutions can help cool the planet - if we act now. Nature 593(7858): 191–194.

31.

Green

Unsworth

RKF

Chadwick

, et al. (2021) Historical analysis exposes catastrophic seagrass loss for the United Kingdom. Frontiers in Plant Science 12: 629962.

32.

Gregg

Elias

Alonso

, et al. (2021) Carbon storage and sequestration by habitat: a review of the evidence (second edition). Natural England Research Report no. NERR094, 20 April. York: Natural England.

33.

Gudasz

Sobek

Bastviken

, et al. (2015) Temperature sensitivity of organic carbon mineralization in contrasting lake sediments. Journal of Geophysical Research Biogeosciences 120(7): 1215–1225.

34.

Hanson

Pace

Carpenter

, et al. (2015) Integrating landscape carbon cycling: Research needs for resolving organic carbon budgets of lakes. Ecosystems 18(3): 363–375.

35.

Heiri

Lotter

Lemcke

(2001) Loss on ignition as a method for estimating organic and carbonate content in sediments: Reproducibility and comparability of results. Journal of Paleolimnology 25(1): 101–110.

36.

Hughes

Hornby

Bennion

, et al. (2004) The development of a GIS-based inventory of standing waters in Great Britain together with a risk-based prioritisation protocol. Water Air and Soil Pollution Focus 4(2–3): 73–84.

37.

Jones

Benson-Evans

Chambers

(1985) Human influence upon sedimentation in Llangorse Lake, Wales. Earth Surface Processes and Landforms 10(3): 227–235.

38.

Kastowski

Hinderer

Vecsei

(2011) Long-term carbon burial in European lakes: Analysis and estimate. Global Biogeochemical Cycles 25(3): GB3019.

39.

Keith

Vardon

Obst

, et al. (2021) Evaluating nature-based solutions for climate mitigation and conservation requires comprehensive carbon accounting. The Science of the Total Environment 769: 144341.

40.

Kortelainen

Pajunen

Rantakari

, et al. (2004) A large carbon pool and small sink in boreal holocene lake sediments: carbon stock in boreal lakes. Global Change Biology 10(10): 1648–1653.

41.

Lang

Bedford

Brooks

, et al. (2010) Early-Holocene temperature variability inferred from chironomid assemblages at Hawes Water, northwest England. The Holocene 20(6): 943–954.

42.

Loisel

Beilman

, et al. (2014) A database and synthesis of northern peatland soil properties and Holocene carbon and nitrogen accumulation. The Holocene 24(9): 1028–1042.

43.

Mackie

EAV

Lloyd

Leng

, et al. (2007) Assessment of δ¹³C and C/N ratios in bulk organic matter as palaeosalinity indicators in holocene and Lateglacial isolation basin sediments, northwest Scotland. Journal of Quaternary Science 22(6): 579–591.

44.

Matthews

(2022) Iron age palaeoenvironments of Northwest Scotland. Unpublished PhD Thesis, Newcastle University, UK.

45.

McGowan

Anderson

Edwards

, et al. (2016) Long-term perspectives on terrestrial and aquatic carbon cycling from palaeolimnology. WIREs Water 3(2): 211–234.

46.

Mcleod

Chmura

Bouillon

, et al. (2011) A blueprint for blue carbon: Toward an improved understanding of the role of vegetated coastal habitats in sequestering CO₂. Frontiers in Ecology and the Environment 9(10): 552–560.

47.

Met Office (2022) UK climate averages. Available at: https://www.metoffice.gov.uk/research/climate/maps-and-data/uk-climate-averages (accessed 12 June 2022).

48.

Moyle

Boyle

Chiverrell

(2021) Towards a history of holocene P dynamics for the Northern Hemisphere using lake sediment geochemical records. Biogeosciences 18(20): 5609–5638.

49.

Mulder

YLA

(1999) Aspects of vegetation and settlement history in the Outer Hebrides, Scotland. PhD Thesis, University of Sheffield, UK.

50.

Mulholland

Elwood

(1982) The role of lake and reservoir sediments as sinks in the perturbed global carbon cycle. Tellus 34(5): 490–499.

51.

Pajunen

(2000) Lake sediments: Their carbon store and related accumulation rates. Geological Survey of Finland, Special Paper 29: 39–69.

52.

Parker

Goudie

Anderson

, et al. (2002) A review of the mid-holocene elm decline in the British Isles. Progress in Physical Geography: Earth and Environment 26(1): 1–45.

53.

Phelan

Shaw

Baylis

(2011) The extent of saltmarsh in England and Wales: 2006–2009. Report for the Environment Agency, 1 August. Bristol: Environment Agency.

54.

Pollard

(2008) The construction of prehistoric Britain. In: Pollard

(ed.) Prehistoric Britain. Oxford: Blackwell, pp.1–17.

55.

Prairie

(2008) Carbocentric limnology: Looking back, looking forward. Canadian Journal of Fisheries and Aquatic Sciences 65(3): 543–548.

56.

Radbourne

Ryves

Anderson

, et al. (2017) The historical dependency of organic carbon burial efficiency. Limnology and Oceanography 62(4): 1480–1497.

57.

Rantala

Luoto

Nevalainen

(2016) Temperature controls organic carbon sequestration in a subarctic lake. Scientific Reports 6(1): 34780.

58.

Ratcliffe

Payne

Sloan

, et al. (2018) Holocene carbon accumulation in the peatlands of northern Scotland. Mires and Peat 23: 1–30.

59.

Reimer

Austin

WEN

Bard

, et al. (2020) The IntCal20 Northern Hemisphere radiocarbon age calibration curve (0–55 cal kBP). Radiocarbon 62(4): 725–757.

60.

Scott

(2014) Carbon fixation, flux and burial efficiency in two contrasting eutrophic lakes in the UK (Rostherne Mere & Tatton Mere). PhD Thesis, Loughborough University, UK.

61.

Smeaton

Austin

WEN

Davies

, et al. (2017) Scotland’s forgotten carbon: A national assessment of mid-latitude fjord sedimentary carbon stocks. Biogeosciences 14(24): 5663–5674.

62.

Snowball

IAN

Thompson

ROY

(2008) A mineral magnetic study of holocene sedimentation in Lough Catherine, Northern Ireland. Boreas 19(2): 127–146.

63.

Sobek

Durisch-Kaiser

Zurbrügg

, et al. (2009) Organic carbon burial efficiency in lake sediments controlled by oxygen exposure time and sediment source. Limnology and Oceanography 54(6): 2243–2254.

64.

Stallard

(1998) Terrestrial sedimentation and the carbon cycle: Coupling weathering and erosion to carbon burial. Global Biogeochemical Cycles 12(2): 231–257.

65.

Taylor

Gilbert

Cooke

, et al. (2019) High carbon burial rates by small ponds in the landscape. Frontiers in Ecology and the Environment 17(1): 25–31.

66.

Tranvik

Downing

Cotner

, et al. (2009) Lakes and reservoirs as regulators of carbon cycling and climate. Limnology and Oceanography 54(6part2): 2298–2314.

67.

UK Centre for Hydrology and Ecology (2016) UK lakes portal. Available at: https://eip.ceh.ac.uk/apps/lakes/index.html (accessed 2 November 2021).

68.

Van Vliet-Lanoë

Helluin

Pellerin

, et al. (1992) Soil erosion in western Europe: From the last interglacial to the present. In: Bell

Boardman

(eds) Past and Present Soil Erosion, vol. 22. Oxford: Oxbow Monograph, pp.101–114.

69.

Verpoorter

Kutser

Seekell

, et al. (2014) A global inventory of lakes based on high-resolution satellite imagery. Geophysical Research Letters 41(18): 6396–6402.

70.

Watkins

RUTH

Scourse

Allen

JRM

(2007) The holocene vegetation history of the Arfon platform, North Wales, UK. Boreas 36(2): 170–181.

71.

Woodbridge

Fyfe

Roberts

, et al. (2014) The impact of the Neolithic agricultural transition in Britain: A comparison of pollen-based land-cover and archaeological 14C date-inferred population change. Journal of Archaeological Science 51: 216–224.

Supplementary Material

Please find the following supplemental material available below.

For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.

For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.

0.00 MB

5.84 MB

0.03 MB