Land-use changes have been shown to have large effects on climate and biogeochemical cycles, but so far most studies have focused on the effects of conversion of natural vegetation to croplands and pastures. By contrast, relatively little is known about the long-term influence of past agriculture on vegetation regrowth and carbon sequestration following land abandonment. We used the LPJ-GUESS dynamic vegetation model to study the legacy effects of different land-use histories (in terms of type and duration) across a range of ecosystems. To this end, we performed six idealized simulations for Europe and Africa in which we made a transition from natural vegetation to either pasture or cropland, followed by a transition back to natural vegetation after 20, 60 or 100 years. The simulations identified substantial differences in recovery trajectories of four key variables (vegetation composition, vegetation carbon, soil carbon, net biome productivity) after agricultural cessation. Vegetation carbon and composition typically recovered faster than soil carbon in subtropical, temperate and boreal regions, and vice versa in the tropics. While the effects of different land-use histories on recovery periods of soil carbon stocks often differed by centuries across our simulations, differences in recovery times across simulations were typically small for net biome productivity (a few decades) and modest for vegetation carbon and composition (several decades). Spatially, we found the greatest sensitivity of recovery times to prior land use in boreal forests and subtropical grasslands, where post-agricultural productivity was strongly affected by prior land management. Our results suggest that land-use history is a relevant factor affecting ecosystems long after agricultural cessation, and it should be considered not only when assessing historical or future changes in simulations of the terrestrial carbon cycle but also when establishing long-term monitoring networks and interpreting data derived therefrom, including analysis of a broad range of ecosystem properties or local climate effects related to land cover changes.
Historically, many natural forests or grasslands on Earth have been cleared
or cultivated for grazing, timber, food production, mining or settlements.
However, land-use change (LUC) in these areas has rarely been continuous,
and land cover and management have often changed for a variety of reasons
(Burgi and Turner, 2002). Based on the HYDE dataset, Campbell et al. (2008)
estimated that 269 Mha of cropland and 479 Mha of pasture
have been abandoned between 1700 and 2000. Recently, agricultural cessation
rates have risen globally, especially in the temperate region. For example,
during the last decades, large areas in Europe previously used for pasture
or crop cultivation have been abandoned (e.g., Schierhorn et al., 2013; Smith
et al., 2005). Following agricultural abandonment, and in the absence of
further anthropogenic influence, natural vegetation recolonizes in a typical
succession from herbaceous vegetation to shrubland and forests, if
environmental conditions are suitable for tree growth. These secondary
forests act as an important carbon (C) sink during the years of regrowth,
thereby reducing the growth rate of global atmospheric CO
The immediate effects of land-use (LU) practices on C fluxes and nutrient cycles have been studied in some detail over recent decades. Generally, agriculture significantly reduces C and, in the absence of supplementary sources, nitrogen (N) pools due to initial deforestation, reduced soil litter input, and accelerated soil decomposition and erosion (Davidson and Ackerman, 1993; Fujisaki et al., 2015; Guo and Gifford, 2002; McLauchlan, 2006; Murty et al., 2002). Pasture soils can be an exception as they have been found to accumulate C, depending on location and management (McSherry and Ritchie, 2013; Milchunas and Lauenroth, 1993). The long-term importance of past LU on ecosystems, however, was recognized only recently, and much less effort has been put so far into the investigation of legacy effects of LU history on ecosystem processes, how long these effects persist, or whether they may even be irreversible (Chazdon, 2014; Compton and Boone, 2000; Cramer et al., 2008; Hobbs et al., 2009; McLauchlan, 2006). This is important not only for understanding present-day ecological systems but also because, due to demographical, social, technological, economic and environmental changes, LUC and land abandonment will continue to occur in the future (Hurtt et al., 2011).
Most observational studies that looked at the recovery of ecosystems after agricultural cessation focused on the first years of succession. Analyses of the long-term effects of historical LU are often limited by the availability of adequate LU information and the absence of undisturbed ecosystems, and usually rely on chronosequences (Chazdon, 2003; Knops and Tilman, 2000). Only a few long-term observational study plots like the one maintained at the Rothamsted Experimental Station (e.g., Poulton et al., 2003) exist. Differences between (near) pristine and post-agricultural forests or grasslands have been reported to persist for decades or centuries after agricultural abandonment for various variables, including soil pH (Falkengren-Grerup et al., 2006); microbial communities (Fichtner et al., 2014); soil C, N and phosphorus (Compton and Boone, 2000); and other nutrients (Wall and Hytonen, 2005). Furthermore, aboveground (ag) biomass (Wandelli and Fearnside, 2015), percentage vegetation cover (Lesschen et al., 2008), biodiversity (Vellend, 2004), species composition (Aide et al., 2000) and structure (Bellemare et al., 2002) remained affected for years to decades, or even longer. These effects have consequences, not only for the C sink capacity of the ecosystem but also for water and energy exchange between the land and the atmosphere (Foley et al., 2003), which also has important, albeit still highly uncertain, implications for regional climate change (e.g., Arora and Montenegro, 2011; Brovkin et al., 2013; de Noblet-Ducoudre et al., 2012). Some studies have detected an influence of ancient agriculture on forest composition and diversity even thousands of years later (Dambrine et al., 2007; Dupouey et al., 2002; Willis et al., 2004). However, the persistence of legacy effects varies considerably with former LU, geographical location, sampling methods and examined variables, making recovery trajectories often hard to predict (Cramer et al., 2008; Foster et al., 2003; Guariguata and Ostertag, 2001; Norden et al., 2015; Post and Kwon, 2000; Suding et al., 2004).
In this study, we performed idealized simulations with the LPJ-GUESS dynamic global vegetation model (DGVM) to explore the importance of agricultural LU history in terms of type and duration for the regeneration of ecosystems and C stocks and fluxes under a range of environmental conditions. We converted natural vegetation to either pasture or cropland, followed by a re-transition to natural vegetation after time periods of 20, 60 and 100 years. While there are numerous variables suitable to measure recovery (Chazdon, 2003; Martin et al., 2016), we analyzed recovery times for vegetation composition (represented here by the dominant plant functional type), vegetation C, soil C, and net biome productivity to evaluate the longevity of the effects of LU history on the C cycle component of ecosystems and to ascertain whether the system eventually recovers to its pre-disturbance state.
LPJ-GUESS is a process-based DGVM that is driven by climate, atmospheric
CO
Conversion of natural to managed land in LPJ-GUESS is characterized by the
initial killing of all living vegetation in the affected area. The
corresponding woody biomass is partly oxidized immediately (67–76 %) and
partly transferred to the product (21 %) or litter (3–12 %) pool.
Ten percent of the leaves are oxidized, while the rest of the leaves and the fine
roots enter the litter pool. Only the litter thus remains in the ecosystem
subsequent to land conversion. Pastures are represented by preventing tree
establishment and wildfires and by splitting the aboveground biomass of the
grasses equally between atmosphere (harvest) and litter at the end of each
year. Crops were represented by grass PFTs modified to mimic aspects of
cropland important for the C and N cycles. Settings for croplands and
pastures were as follows:
For transitions from natural vegetation to cropland, we transferred only
3 % of the cleared woody biomass to the litter instead of 12 % for
natural vegetation–pasture transitions. This accounts for the practice
that farmers would try to remove as many coarse roots as possible before
planting of crops. Harvest efficiency (in this study: fraction of aboveground biomass that is
oxidized) was 0.5 yr While we removed 100 % of harvested N biomass for croplands, we changed
this value to 65 % for pastures. That accounts for significant urine N
regain from animals fed on pastures (Dean et al., 1975; Lauenroth and Milchunas, 1992). Root turnover rate was 0.7 yr In croplands we estimated tillage effects by increasing heterotrophic
respiration by a factor of 1.94 (Pugh et al., 2015). We simulated N fertilization in croplands by applying 75 kg ha
After patch-destroying disturbances or managed land converting back to
natural vegetation, there is a typical succession from grasses to
light-demanding pioneer trees, eventually followed in many ecosystems by the
establishment of shade-tolerant PFTs. It has been shown that LPJ-GUESS is
able to realistically simulate observed succession pathways and species
variations (Hickler et al., 2004; Smith et al., 2014).
During spin-up (500 years) and the simulation period (900 years), we forced
LPJ-GUESS with temperature-detrended, repeated 1981–2000 climate from the
University of East Anglia Climate Research Unit 3.21 dataset
(CRU, 2013), 1990s mean N deposition (Lamarque et al.,
2013) and a fixed atmospheric CO
To facilitate the interpretation, we classified each grid cell to one biome.
We used the same classification rules as Smith et al. (2014), aggregated
to eight biomes as in Bayer et al. (2015). Afterwards, we excluded grid cells
from the analyses which were classified as desert or tundra, had a mean net
primary productivity (NPP) below 0.1 kg C m
We studied the influence of LU history on ecosystems by analyzing four key
variables: dominant PFT, vegetation C, soil C (excluding litter) and net
biome productivity (NBP). NBP is the net atmosphere–land carbon flux after
C losses associated with respiratory fluxes, fire, harvest, land clearing and
decomposition of LUC product pools are subtracted from gross primary
productivity. We investigated the legacy effects of LU history by
calculating a recovery time for each variable, simulation and grid cell
after the conversion back to natural vegetation. For vegetation C, soil C
and NBP, recovery time was defined as the year in which the 20-year running
mean of the variable exceeded the threshold of one standard deviation (
For the dominant PFT recovery, we first identified which PFT dominates each
grid cell in the reference simulation based on the annual maximum leaf area
index (LAI) amongst PFTs. We then checked for dominant PFT recovery in the same
way as we did for vegetation C, soil C and NBP (i.e., whether its LAI exceeded the
threshold of 1
Vegetation C (kg C m
Maps of simulated vegetation and soil C, as well as dominant PFT and biomes
derived from PFT composition for the reference simulation, are shown in Fig. 1.
The salient features of biome and C storage distribution at the regional
scale are captured (Haxeltine and Prentice, 1996; Scharlemann et al.,
2014). Vegetation C reaches its highest values in tropical forests of
central Africa and decreases towards the deserts of southern and northern
Africa. Patterns are more homogeneous in Europe, where most areas store
5–10 kg C m
The LAI of the dominant PFT recovers on average within around one century
for all LU histories (Fig. 2). Maps of the recovery time (Fig. 3) show
distinct geographical patterns which occur in all simulations. Most
subtropical grasslands and savannas, and parts of the temperate and boreal
forests, recover within several decades, some grasslands even within 5
years. In contrast, recovery times are clearly longer (
Overall, differences across simulations of different LU histories are moderate, with generally only small differences in temperate forests, savannas and shrublands (Fig. 3; see also biome averages in Table 1 and the histogram in Fig. A2). Areas of major differences are central Africa, where P20 recovers faster than other simulations because post-agricultural net mineralization rates are higher in this region for P20 than for the other simulations (Fig. 4), thereby relatively increasing post-agricultural N availability compared to the other simulations (Fig. 5), and the African Mediterranean coast, where croplands recover much faster because the reduced C : N ratio in the soil (not shown) enhances N mineralization and thus plant N availability compared to pastures. Furthermore, in parts of the boreal zone recovery takes several hundred years for C100 instead of a few decades for the other simulations because lower available N levels relatively reduce the growth of IBS (the dominant PFT in this region) compared to other woody PFTs. Figure 6 shows the maximum differences between recovery times across all simulations per biome (black dots), as well as across a subset of simulations (colored squares and triangles). The differences were first calculated for each grid cell and only then averaged over biomes, thereby providing an estimate of the relative importance of former LU duration versus former LU type on recovery times. While substantial differences occur across the pasture simulations (P20, P60, P100) in tropical forests, savannas and grasslands, and across cropland simulations (C20, C60, C100) in boreal forests (emphasizing the importance of LU duration in these regions), major differences between P100 and C100 occur in boreal forests and grasslands (emphasizing the importance of LU type if agricultural duration was long). On the other hand, in our simulations, dominant PFT recovery in temperate forests is hardly influenced by the type of former LU or, conversely, pasture duration has negligible effects on boreal forest recovery. Interestingly, temperate forests recover faster for P100 and C100 then for P20 and C20. This pattern is generally restricted to areas where the TeBS PFT dominates. We interpret this behavior as reduced soil N favoring TeBS in the competition with other tree PFTs, thereby reaching its background LAI levels earlier.
Compared to dominant PFT, recovery occurs slightly later for vegetation C (Fig. 2, Table 1). Spatial patterns look more homogeneous than for the dominant PFT (Fig. 3). While most grasslands recover within a few decades for all simulations, in particular so for post-cropland recovery, recovery occurs only after several decades or centuries in forest ecosystems. Lower standard deviations for the mean differences in vegetation C recovery times compared to the standard deviations for the mean differences in dominant PFT recovery times for most biomes (Fig. 6a and b) reflect the spatially more uniform response of vegetation C recovery. Exceptions are tropical forests and grasslands, where the standard deviation is higher for vegetation C recovery compared to dominant PFT recovery.
Time series (20-year running mean) of dominant PFT, vegetation C,
soil C and NBP for the different experiments, starting from the time of
reconversion to natural vegetation and area-averaged over all grid cells.
Dominant PFT, vegetation C and soil C are shown in relative values compared
to reference simulation mean, while NBP is shown as absolute values
(kg C m
Maps of recovery times in years for the dominant PFT, vegetation C, soil C, and NBP for the P20, P100, C20, and C100 simulations.
Average net N mineralization rates (kg N ha
Significant differences in recovery times occur between simulations of different LU types that have the same duration, and between simulations of the same LU type but with different duration. For example, in the grasslands and savannas of southern, eastern and northern Africa, former croplands recover much faster than former pastures (see also Table 1 and Fig. A2) because post-agricultural N availability is enhanced in these regions (Fig. 5). In former croplands in these environments, the combined effect of fertilizing and harvest is a net N flux to the ecosystem (not shown) and mineralization rates are enhanced after cropland abandonment (Fig. 4). This net N flux can partially be explained by high levels of water stress in these savannas and grasslands, resulting in greater C and N allocation to roots relative to leaves and thereby decreased harvest removal in this region (Fig. A3). Conversely, recovery in northern European forests is delayed for C60 and, to an even greater extent, C100 because in this region N removal by annual harvest exceeds N addition through fertilization during the agricultural period (not shown) and post-agricultural N mineralization rates in this region are substantially reduced compared to the other simulations many decades or even a few centuries after abandonment (Fig. 4). Differences in vegetation recovery times resulting from agricultural duration are mostly found in temperate and boreal forests for the cropland simulations (here longer durations result in longer recovery times due to reduced N availability, Fig. 5) and in tropical forests and shrublands for the pasture simulations, emphasizing the importance of agricultural duration in these regions (see also Fig. 6b).
Relative depletion of soil C content under crop and pasture LU is not as large (loss of 0–11 % compared to the reference simulation) as for vegetation C (Fig. 2). However, regeneration proceeds over longer timescales due to slower C accumulation in soils than in vegetation. C depletion is generally more pronounced for former crops than for pastures due to the greater harvest efficiency, which leads to more biomass removed each year, and the effect of tillage enhancing soil respiration (Sect. 2.1). Upon re-conversion, soil C accumulation is delayed for the pasture simulations compared to the cropland simulations, especially for P20, where the residual roots and other litter left after the original deforestation event continue to decay and soil C decreases for some decades. The general delay for pastures is associated with larger heterotrophic respiration rates (not shown) compared to rates calculated in recovering croplands.
Soil C recovery rates are highly latitude-dependent (Fig. 3), being much
slower in temperate (
Soil C recovery times differ substantially between simulations in many areas. LU type is particularly important in grasslands and non-tropical forests. While croplands tend to recover faster than pastures in grasslands of southern and northern Africa, the opposite occurs in most temperate and boreal forests but also the northern Sahel, where soil C after re-conversion from croplands does not recover at all. Post-agricultural N availability is enhanced in parts of the Sahel for the cropland simulations due to increased N mineralization rates (Figs. 4 and 5), and trees benefit more than grasses, leading to a shift in the equilibrium vegetation state towards woody species (not shown), which results in an overall lower soil C pool size. It should be noted that even though some regions do not recover within 800 years, a large fraction of the original C loss is already replenished after a few centuries, thereby limiting implications for the C cycle. Counter to a priori expectations, for tropical and temperate forests and for shrublands, the difference between P20 and C20 is usually higher than between P60 and C60 or P100 and C100 (Fig. 6c). Pasture duration is relevant for speed of soil C recovery in most ecosystems and, apart from in the tropics, a longer duration usually delays recovery, mainly due to substantial initial depletions after long pasture durations (Fig. 2). For croplands, longer durations tend to delay recovery in temperate and boreal forests but accelerate soil C recovery in the (sub)tropics. This is somewhat unexpected for the tropical forest biome, where longer cropland durations usually do not increase N availability upon abandonment in our simulations (Fig. 5). However, while tropical soils lose large amounts of C during the first decades of cropland use, slow C accumulation takes place thereafter, resulting in higher soil C values at the end of the agricultural period for C100 than for C20 in large parts of eastern Africa. This occurs because tillage-driven C losses in more labile soil pools, which dominate the system's response during the first decades, are eventually supplanted as the dominant process by accumulation in more stable pools. This is different to temperate and boreal forest, where soil C decreases throughout the entire cropland period. Overall, the greatest sensitivity of soil C recovery times to different LU histories is found in boreal forests and grasslands, where maximum differences across simulations are often several centuries (Fig. 6c). The maximum differences across all simulations (P20/P60/P100/C20/C60/C100) in boreal forests are mainly due to differences across simulations of same LU type but different duration (e.g., P20/P60/P100), whereas the sensitivity of grasslands mainly reflects differences across simulations of different LU type but same duration (e.g., P100/C100), emphasizing the importance of duration and type of agriculture in a range of biomes.
Average N limitation on vegetation RuBisCO capacity (and thus on gross primary production) for the reference simulation (full simulation period) and during the first 100 years of regrowth for the P20, P100, C20, and C100 simulations. N limitation is a number scaling from 0 (completely N-limited) to 1 (no N limitation) (Smith et al., 2014).
Maximum difference in recovery time (longest recovery time minus shortest recovery time of all selected simulations) for the dominant PFT, vegetation C, soil C, and NBP. Black dots show maximum differences across all six simulations (P20, P60, P100, C20, C60, C100), green squares differences across 20-year pasture and cropland simulations (P20, C20), blue squares differences across 60-year pasture and cropland simulations (P60, C60), red squares differences across 100-year pasture and cropland simulations (P100, C100), orange triangles differences across pasture simulations (P20, P60, P100), and purple triangles differences across cropland simulations (C20, C60, C100). Background colors indicate associated biomes, arrows one standard deviation, and the dashed line 0 years' difference. Thus, the black dots show the sensitivity of recovery times to LU history across all simulations for each biome. The red, blue and green squares indicate the relative contribution of LU type for a specific LU duration to this sensitivity, and the orange and purple squares indicate the relative contributions of pasture and of cropland duration. For example, if recovery times for one variable in one grid cell were to be 50, 60, 65, 90, 100, 110 years (for P20, P60, P100, C20, C60, C100), the maximum difference in recovery time across all simulations (black) would be 60 years, across the 20-year simulations (green) 40 years, across the 60-year simulations (blue) 40 years, across the 100-year simulations (red) 45 years, across the pasture simulations (orange) 15 years and across the cropland simulations (purple) 20 years.
NBP switches from being a C source to the atmosphere during the period of land management to a C sink after reconversion to natural vegetation (Fig. 2). The sink capacity of the recovering ecosystem is greatest during the first decades and then gradually returns to the NBP levels of the reference simulation. P20 and, to a lesser extent, C20 act as a smaller sink than the other simulations at least during the first 100 years of regrowth. Recovery generally occurs slower in temperate and boreal regions than in the tropics for all simulations (Fig. 3). Apart from boreal forests, standard deviations of mean differences in recovery times are very small in all biomes compared to the other variables (Fig. 6d). Recovery times are often somewhat lower than those which would be expected from vegetation and soil C recovery times. This is because the greater standard deviation of NBP in our reference simulation (Fig. 2) reduces the threshold value in our recovery definition, thereby making it easier to reach recovery levels for NBP. We discuss the implications of this further in Sect. 4.2.
Differences in NBP recovery times between simulations are relatively small (typically a few years to decades; see Table 1). The largest differences in recovery times are found in the boreal forests between the cropland simulations, and, as for soil C, the differences are often greater between P20 and C20 than between P100 and C100 (Fig. 6d).
Average recovery times and standard deviations per biome and for each simulation. Recovery times are depicted in Fig. 3.
The effects of forest conversion to croplands or pastures are relatively
well studied. Tilled croplands typically show large depletions of soil C
compared to natural forest vegetation, but the picture for pasture is more
diverse (Davidson and Ackerman, 1993; Don et al., 2011; Guo and Gifford,
2002). Table 2 summarizes recent reviews about observed soil C changes in
agriculture compared to our results. LPJ-GUESS tends to simulate lower
C loss in croplands than commonly reported in observations. We attribute this
to a combination of the observation's focus on the top soil (while in
LPJ-GUESS soil C is implicitly averaged over the whole soil column) and our
relatively high fertilizer rates increasing productivity and thereby C input
to the soil. Pugh et al. (2015) studied the C dynamics of soils in
managed lands in LPJ-GUESS and found C accumulation even after 100 years of
grazed pasture at some locations, especially for low atmospheric CO
In contrast to studies of LU effects compared to previously natural ecosystems, the regeneration of ecosystems after agricultural abandonment has been studied less, and a direct comparison to our simulations is challenging, either because limited information about former LU or reference conditions was provided in these studies or because there are important differences from our setup in terms of management and LU duration or other site-specific characteristics. Additionally, most of the available studies were conducted in Amazonia or North America (Don et al., 2011) and there is large variability in physical and biotic characteristics as well as in land management (Kauffman et al., 2009). Many studies focus on the recovery of biodiversity or species richness (Cramer et al., 2008; Queiroz et al., 2014), but these variables cannot be adequately captured by our large-scale PFT approach. It is often assumed that the ecosystem will gradually return to its previous state and that intensive LU delays recovery but the timescales are widely unknown and differ across variables and regions, e.g., tropical species composition recovers much slower than forest structure and soil nutrients (Chazdon, 2003). Different recovery processes are strongly interlinked, e.g., vegetation accumulation and turnover are key factors in the replenishment of soil quality and nutrients which in turn determine plant productivity, and post-agricultural soil C and N dynamics have been shown to correlate during the regeneration of ecosystems (Knops and Tilman, 2000; Li et al., 2012).
Table 2 includes several studies about ecosystem vegetation and soil recovery after agricultural abandonment. Overall, the studies that looked at vegetation recovery upon abandonment indicate that biomass accumulation slows down after some decades and that accumulation rates correlate negatively with agricultural duration. Our simulations show that the rate of vegetation C sequestration indeed declines over time and that longer LU durations delay recovery in each of the analyzed biomes. Observations also indicate that use of land for pasture delays recovery in the tropics upon pasture abandonment compared to cropping, but in our simulations this seems to be the case only after long agricultural durations. For studies about soil C dynamics after agricultural abandonment, interpretation is often hindered by combining different soil layers or aggregating different LU types (Li et al., 2012) and by large variations observed across studies (Post and Kwon, 2000). Nevertheless, most of the observed patterns are reproduced in our simulations, suggesting that LPJ-GUESS captures the salient processes: after abandonment, croplands accumulate C faster than pastures, and recovery often takes more than a century. The impact of LU duration has rarely been studied; however, our results suggest that even though longer agricultural durations mostly result in greater initial soil C depletions, recovery can occur at similar or even faster speed in the subtropics and tropics. In temperate and boreal forests long LU durations tend to delay recovery.
Observations and LPJ-GUESS results of soil C changes during agriculture (cropland and/or pasture) and vegetation and soil C recovery after abandonment.
Observations and LPJ-GUESS results of soil C changes during agriculture (cropland and/or pasture) and vegetation and soil C recovery after abandonment.
The LPJ-GUESS model has been successfully tested against a range of observations and observation-based products, including vegetation distribution and dynamics and soil C response to changes in vegetation cover (Hickler et al., 2004; Miller et al., 2008; Pugh et al., 2015; Smith et al., 2014). In our simulations, we used only two different agricultural land cover types (intensive grazing and fertilized, tilled crops). Our analysis would therefore not identify effects of, for instance, clearing technique (e.g., burning compared to mechanical removal) or different land management practices (e.g., repeated burning or irrigation) within one land cover type. For example, recovery of species richness and maximum tree height of secondary forests occurs faster under no tillage compared to heavy tillage (Randriamalala et al., 2012).
Our study is intended as an idealized experiment to highlight the importance of LU history on ecosystem state and fluxes across biomes. Still, some processes with the potential to affect post-agricultural ecosystem recovery, at least regionally, are not currently included in LPJ-GUESS. One aspect is the phosphorus cycle, which is not implemented in LPJ-GUESS, even though it can be significantly altered by LUC (MacDonald et al., 2012; McLauchlan, 2006). Moreover, while C and N cycles interact in LPJ-GUESS (Smith et al., 2014), the uniform annual fertilizer rate we applied in this study might be realistic in some regions, such as parts of Europe, but exceeds present-day fertilizer use in Africa (Potter et al., 2010). Seed availability, remnant trees and resprouting from surviving roots are important factors during initial stages of tree colonization following agricultural cessation (Bellemare et al., 2002; Cramer et al., 2008). While LPJ-GUESS does not account for these effects explicitly, seedling establishment is limited by a suitable growth environment, such that effects like re-sprouting or remnant trees as seed sources are mimicked. The model has been shown to, for example, reproduce vegetation recolonization in northern Europe during the Holocene well (Miller et al., 2008), as well as canopy structural changes as a function of forest age (Smith et al., 2014). What is more, by using a prescribed climate in our simulations, hydrological biosphere–atmosphere interactions and feedbacks are not captured (Eltahir and Bras, 1996; Giambelluca, 2002), which could alter regional climate in response to land cover change, potentially affecting recovery rates, especially in tropical regions. Biophysical effects are not restricted to modifications of the water cycle but also include changes in surface albedo and roughness length as a function of ecosystem structure and composition, thereby affecting air mixing and heat transfer. While forests generally absorb more sunlight than grasslands (e.g., Culf et al., 1995), differences amongst tree species and age classes exist as well. Substantial impacts related to realistic land-use have been found on local-to-regional scales (Alkama and Cescatti, 2016; Peng et al., 2014). Whether or not the locally observed changes translate to a significant global radiative forcing is still debated as the direction of change differs across regions in some climate models, which may cancel when integrated globally (Pielke et al., 2011). Additionally, while we focus on C sequestration rates in our analysis, there might be biogeochemical implications beyond C. For instance, the emissions of biogenic volatile organic compounds (BVOCs) to the atmosphere vary greatly amongst plant species (Kesselmeier and Staudt, 1999). BVOCs affect atmospheric composition and climate via ozone production, lengthening the lifetime of atmospheric methane, and contributing to secondary organic aerosol formation (Penuelas and Staudt, 2010; Wu et al., 2012). BVOC emission factors might also be drastically influenced by wildfires (Ciccioli et al., 2014), which in turn are driven by species composition and vegetation density. Thus, different successional trajectories of ecosystem structure and composition recovery have the potential to directly modify air quality and climatic conditions under which regrowth occurs, potentially creating positive or negative climate system feedbacks.
The term recovery is subjective and, in the absence of a universal definition amongst ecologists, several definitions currently exist. The definition used in this study examines recovery from a C sequestration perspective which does not capture situations, for example, where the system approaches a new equilibrium (as soil C did in some regions in the cropland simulations). In order to obtain a better understanding of the uncertainties related to our definition we therefore explored four alternative plausible recovery definitions.
When applying a mean minus 2
Recovery based on percentage change (Fig. A4) results in more heterogeneous
patterns across variables when compared to our standard recovery definition.
Applying a threshold of 95 % of the mean, instead of a mean minus
1
By expanding our standard recovery definition by an upper threshold
(reference mean plus 1
Finally, when using the mean
Altogether, the alternative recovery definitions agree on the general findings when applying our standard definition, especially in terms of relative recovery rates. For all definitions, vegetation C and dominant PFT recover faster in grasslands than in forest-dominated ecosystems, and soil C recovery takes much longer in higher latitudes. However, some areas, especially in the subtropics, “recover” from values higher than in the reference simulation, and these cases are not captured by our standard definition. Additionally, in the tropics, soil C accumulation sometimes does not stop once background values are reached and soil C leaves the reference range. When recovery is defined based on standard deviation, NBP recovery is often quicker than recovery of the C pools. This inconsistency emphasizes the importance of both recovery definition and selected variables when studying the recovery of ecosystems (Jones and Schmitz, 2009). This is particularly relevant for flux tower measurements, where an underlying long-term trend caused by the recovery from previous, often unquantified or unknown LU change, might be overlooked due to a large interannual variability in net ecosystem exchange.
Most studies which have explored the effects of distant human activities on
present-day ecosystems were restricted by sampling difficulties, small
spatial scales, short time periods since abandonment, and little information
about background conditions or the specific LU history of the site. Here, we
use a model-based approach to study the legacy effects of agricultural LU
history (type and duration) on ecosystem regeneration and C sink capacity
after the cessation of agriculture in a range of biomes across Europe and
Africa. The model reproduces qualitatively the response found at study
locations, including distinct differences in recovery between different
variables of the terrestrial carbon cycle. Long-lasting legacy effects of
former agricultural intensity emerge as important for present-day ecosystem
functions. These findings have implications for various scientific
applications:
Long-term monitoring sites (e.g., FLUXNET) and Earth observation systems need
to collect and maintain detailed information about past and present land
cover and land management to adequately interpret their data. Assessments of trends in data from sites that seek to identify impacts of
climate change and/or increasing atmospheric CO Simulation experiments need to move beyond deforestation but also
represent, in a more detailed manner, re-growth dynamics following agricultural
abandonment at the sub-grid level. At the moment a few DGVMs have started to
do so (Shevliakova et al., 2009; Stocker et al., 2014; Wilkenskjeld et
al., 2014) based on model products of tropical shifting cultivation
(Hurtt et al., 2011), but accounting for gross land cover changes is also
important in other regions like Europe (Fuchs et al., 2015). Failure to
consider LU history may lead to errors in the simulation of vegetation
properties, potentially resulting in biases in carbon sequestration or
energy balance calculations, with subsequent implications for simulations of
regional and global climate. Our study suggests that, for vegetation and soil C
studies, accounting for LUC over the last 100–150 years is sufficient in
the tropics, while more than 200 years
might be necessary in the temperate and boreal zone; studies restricted to vegetation should not have to
account for LUC more than 150 years ago in any major climatic zone.
Assessing the efficiency of climate mitigation through large-scale
reforestation or afforestation projects will require knowledge about the
type and duration of previous LU. Our simulations suggest that the potential
to rapidly sequester C in biomass and soil is greatest in tropical forests
following short periods of cropland, while boreal forests accumulate C
slowest, especially when previously used for pasture. Special attention
should be paid to monitoring changes in belowground C, as in most places
the accumulation of soil C is much more sensitive to LU history than
C accumulation in re-growing trees. In terms of soil C, our results suggest that some subtropical regions might
not recover at all on timescales relevant for humans. However, given the low
absolute amounts of C “missing” in these soils, implications for the global
C cycle are expected to be small.
Researchers interested in the LPJ-GUESS source code can contact the model
developers (
Plant functional types used in this study.
Soil C for the six simulations after conversion to natural
vegetation at one single example site to illustrate how recovery time was
calculated according to our definition. The cyan-shaded area corresponds to
reference simulation mean
Histograms of recovery times for the dominant PFT, vegetation C, soil C, and NBP for the six experiments. Colors indicate different biomes.
Annual ratio of C removed by harvest and C stored in vegetation, averaged over the whole agricultural period and for P60. As only aboveground biomass is harvested, lower values indicate increased C allocation to roots compared to leaves due to limited water supply.
Maps of recovery time for the dominant PFT, vegetation C, soil C and
NBP with an alternative recovery definition for the P60 and C60
simulations. The definition is the same as our standard definition but with a
mean
Maps of recovery time for the dominant PFT, vegetation C, soil C and
NBP with an alternative recovery definition for the P60 and C60
simulations. The definition is the same as our standard definition but with a
mean
This work was funded by the Helmholtz Association through the International Research Group CLUCIE and by the European Commission's 7th Framework Programme, under grant agreement number 603542 (LUC4C). A. D. Bayer acknowledges support by the European Commission's 7th Framework Programme, under grant agreement number 308393 (OPERAs). This work was supported, in part, by the German Federal Ministry of Education and Research (BMBF), through the Helmholtz Association and its research program ATMO. It also represents paper number 20 of the Birmingham Institute of Forest Research. The article processing charges for this open-access publication were covered by a Research Centre of the Helmholtz Association. Edited by: V. Arora two anonymous referees