Research Article

Ecology

Carbon recovery dynamics following disturbance by selective logging in Amazonian forests

UMR EcoFoG (Agroparistech, CNRS, Inra, Université des Antilles, Cirad), French Guiana
UMR EcoFoG (Agroparistech, CNRS, Inra, Université des Antilles, Université de Guyane), French Guiana
UMR EcoFoG (Agroparistech, Inra, Université des Antilles, Université de Guyane, Cirad), French Guiana
Cirad, UR Forests and Societies, France
Embrapa Amazônia Oriental, Brazil
Wageningen University, Netherlands
University of Florida, United States
CarbonForExpert, Switzerland
University of Oxford, United Kingdom
Instituto Boliviano de Investigación Forestal, Bolivia
Embrapa Amazônia Ocidental, Brazil
Florida International University, United States
Embrapa Amapa, Brazil
Instituto de Investigaciones de la Amazonia Peruana, Peru
Embrapa Acre, Brazil
University of São Paulo, Brazil

Dec 20, 2016

Open access
Copyright information

Abstract
eLife digest
Introduction
Results
Discussion
Materials and methods
Appendix
Data availability
References
Article and author information
Metrics

Abstract

When 2 Mha of Amazonian forests are disturbed by selective logging each year, more than 90 Tg of carbon (C) is emitted to the atmosphere. Emissions are then counterbalanced by forest regrowth. With an original modelling approach, calibrated on a network of 133 permanent forest plots (175 ha total) across Amazonia, we link regional differences in climate, soil and initial biomass with survivors’ and recruits’ C fluxes to provide Amazon-wide predictions of post-logging C recovery. We show that net aboveground C recovery over 10 years is higher in the Guiana Shield and in the west (21 $\pm$ 3 Mg C ha $^{- 1}$ ) than in the south (12 $\pm$ 3 Mg C ha $^{- 1}$ ) where environmental stress is high (low rainfall, high seasonality). We highlight the key role of survivors in the forest regrowth and elaborate a comprehensive map of post-disturbance C recovery potential in Amazonia.

https://doi.org/10.7554/eLife.21394.001

eLife digest

The Amazon rainforest in South America is the largest tropical forest in the world. Along with being home to a huge variety of plants and wildlife, rainforests also play an important role in storing an element called carbon, which is a core component of all life on Earth. Certain forms of carbon, such as the gas carbon dioxide, contribute to climate change so researchers want to understand what factors affect how much carbon is stored in rainforests. Trees and other plants absorb carbon dioxide from the atmosphere and then incorporate the carbon into carbohydrates and other biological molecules. The Amazon rainforest alone holds around 30% of the total carbon stored in land-based ecosystems.

Humans selectively harvest certain species of tree that produce wood with commercial value from the Amazon rainforest. This “selective logging” results in the loss of stored carbon from the rainforest, but the loss can be compensated for in the medium to long term if the forest is left to regrow. New trees and trees that survived the logging grow to fill the gaps left by the felled trees. However, it is not clear how differences in the forest (for example, forest maturity), environmental factors (such as climate or soil) and the degree of the disturbance caused by the logging affect the ability of the forest ecosystem to recover the lost carbon.

Piponiot et al. used computer modeling to analyze data from over a hundred different forest plots across the Amazon rainforest. The models show that the forest’s ability to recover carbon after selective logging greatly differs between regions. For example, the overall amount of carbon recovered in the first ten years is predicted to be higher in a region in the north known as the Guiana Shield than in the south of the Amazonian basin where the climate is less favorable.

The findings of Piponiot et al. highlight the key role the trees that survive selective logging play in carbon recovery. The next step would be to couple this model to historical maps of logging to estimate how the areas of the rainforest that are managed by selective logging shape the overall carbon balance of the Amazon rainforest.

https://doi.org/10.7554/eLife.21394.002

Introduction

With on-going climate change, attention is increasingly drawn to the impacts of human activities on carbon (C) cycles (Griggs and Noguer, 2002), and in particular to the 2.1 $\pm$ 1.1 Pg C yr $^{- 1}$ of C loss caused by various forms and intensities of anthropogenic disturbances in tropical forests (Grace et al., 2014). Among those disturbances, selective logging, i.e. the selective harvest of a few merchantable tree species, is particularly widespread: in the Brazilian Amazon alone, about 2 Mha yr $^{- 1}$ were logged in 1999–2002 (Asner et al., 2005). The extent of selective logging in the Brasilian Amazon was equivalent to annual deforestation in the same period, and resulted in C emissions of 90 Tg C yr $^{- 1}$ (Huang and Asner, 2010) which increased anthropogenic C emissions by almost 25% over deforestation alone (Asner et al., 2005). In contrast to deforested areas that are used for agriculture and grazing, most selectively logged forests remain as forested areas (Asner et al., 2006) and may recover C stocks (West et al., 2014). Previously logged Amazonian forests may thus accumulate large amounts of C (Pan et al., 2011), but this C uptake is difficult to accurately estimate, because while detecting selective logging from space is increasingly feasible (Frolking et al., 2009) (even if very few of the IPCC models effectively account for logging), directly quantifying forest recovery remains challenging (Asner et al., 2009; Houghton et al., 2012; Goetz et al., 2015). Studies based on field measurements (e.g. Sist and Ferreira, 2007; Blanc et al., 2009; West et al., 2014; Vidal et al., 2016), sometimes coupled with modeling approaches (e.g. Gourlet-Fleury et al., 2005; Valle et al., 2007) or airborne light detection and ranging (LiDAR) measurements (e.g. Andersen et al., 2014) have assessed post-logging dynamics at particular sites. Nonetheless, to our knowledge no spatially-explicit investigation of post-logging C dynamics at the Amazon biome scale is available.

C losses from selective logging are determined by harvest intensity (i.e. number of trees felled or volume of wood extracted) plus the care with which harvest operations are conducted, which affects the amount of collateral damage. After logging, C losses continue for several years due to elevated mortality rates of trees injured during harvesting operations (Shenkin et al., 2015). Logged forests may recover their aboveground carbon stocks (ACS) via enhanced growth of survivors and recruited trees (Blanc et al., 2009). Full recovery of pre-disturbance ACS in logged stands reportedly requires up to 125 years, depending primarily on disturbance intensity (Rutishauser et al., 2015). The underlying recovery processes (i.e. tree mortality, growth and recruitment) are likely to vary with the clear geographical patterns in forest structure and dynamics across the Amazon Basin and Guiana Shield. In particular, northeast-southwest gradients have been reported for ACS (Malhi and Wright, 2004), net primary productivity (Aragão et al., 2009), wood density (Baker et al., 2004), and floristic composition (ter Steege et al., 2006). Such gradients coincide with climate and edaphic conditions that range from nearly a seasonal nutrient-limited in the northeast to seasonally dry and nutrient-rich in the southwest (Quesada et al., 2012). These regional differences in biotic and abiotic conditions largely constrain demographic processes that ultimately shape forest C balances.

Here we partition the contributions to post-disturbance ACS gain (from growth and recruitment of trees $\geq$ 20 cm DBH) and ACS loss (from mortality) of survivors and recruited trees to detect the main drivers and patterns of ACS recovery in forests disturbed by selective logging across Amazonia sensu lato (that includes the Amazon Basin and the Guiana Shield). Based on long-term (8–30 year) inventory data from 13 experimentally-disturbed sites (Sist et al., 2015) across Amazonia (Figure 1—figure supplement 1), 133 permanent forest plots (175 ha in total) that cover a large gradient of disturbance intensities (ACS losses ranging from 1% to 71%) were used to model the trajectory of those post-disturbance ACS changes (Figure 1) in a comprehensive Bayesian framework. We quantify the effect of pre-disturbance ecosystem characteristics [the site’s average pre-logging ACS ( $a c s 0$ ) and the relative difference between each plot and $a c s 0$ as a proxy of forest maturity ( $d a c s$ )], disturbance intensity [percentage of pre-logging ACS lost ( $l o s s$ )], and interactions with the environment [annual precipitation ( $p r e c$ ), seasonality of precipitation ( $s e a s$ ), and soil bulk density ( $b d$ )] (Figure 2) on the rates at which post-disturbance ACS changes converge to a theoretical steady state (as in Figure 1, see Materials and methods for more details). With global maps of ACS (Avitabile et al., 2016), climatic conditions (Hijmans et al., 2005) and soil bulk density (Nachtergaele et al., 2008), we up-scale our results to Amazonia (sensu lato) and elaborate predictive maps of potential ACS changes over 10 years under the hypothesis of a 40% ACS loss, which is a common disturbance intensity after conventional logging in Amazonia (Blanc et al., 2009; Martin et al., 2015; West et al., 2014). Summing these ACS changes over time gives the net post-disturbance rate of ACS accumulation. Disentangling ACS recovery into demographic processes and cohorts is essential to reveal mechanisms underlying ACS responses to disturbance and to make more robust predictions of ACS recovery compared to an all-in-one approach (see Appendix).

Figure 1 with 1 supplement see all

Download asset Open asset

Post-disturbance annual ACS changes of survivors and recruits in 133 Amazonian selectively logged plots.

Data is available between the year of minimum ACS ( $t = 0$ ) and $t =$ 30 years. ACS changes are: recruits’ ACS growth (orange), recruits’ ACS loss (gold), new recruits’ ACS (red), survivors’ ACS growth (light green) and survivors’ ACS loss (dark green). Thick solid lines are the maximum-likelihood predictions (for an average plot, when all covariates are null), and dashed lines are the model theoretical behaviour. New recruits’ ACS, recruits’ ACS growth, and recruits’ ACS loss converge over time to constant values. A dynamic equilibrium is then reached: ACS gain from recruitment and recruits’ growth compensate ACS loss from recruits’ mortality. Survivors’ ACS growth and loss. decline over time and tend to zero when all initial survivors have died.

https://doi.org/10.7554/eLife.21394.003

Figure 2 with 1 supplement see all

Download asset Open asset

Effect of covariates on the rate at which post-disturbance ACS changes converge to a theoretical steady state (in yr $^{- 1}$ ).

Covariates are : disturbance intensity ( $l o s s$ ) , i.e. the proportion of initial ACS loss; mean site’s ACS ( $a c s 0$ ), and relative forest maturity, i.e. pre-logging plot ACS as a % of $a c s 0$ ( $d a c s$ ); annual precipitation ( $p r e c$ ); seasonality of precipitation ( $s e a s$ ), soil bulk density ( $b d$ ). Covariates are centred and standardized. Red and black levels are 80% and 95% credible intervals, respectively. The median rate is the prediction of the convergence rate for an average plot (when all covariates are set to zero). Negative covariate values indicate slowing and positive values indicate accelerating rates. (a) Survivors’ ACS growth. (b) New recruits’ ACS. (c) Recruits’ ACS growth. (d) Survivors’ ACS loss. (e) Recruits’ ACS loss.

https://doi.org/10.7554/eLife.21394.005

Figure 2—source data 1 Parameters posterior distribution. Columns are the 2.5%, 10%, 50%, 90% and 97.5% quantiles of the posterior distribution of the model parameters (rows).: https://doi.org/10.7554/eLife.21394.006
Download elife-21394-fig2-data1-v2.xlsx

Results

Local variations of ACS changes

At a given site, variations of post-logging ACS changes are explained with the disturbance intensity ( $l o s s$ ) and the relative forest maturity ( $d a c s$ ). At high disturbance intensity (positive $l o s s$ ) as well as in relatively immature forests (negative $d a c s$ ), ACS gain from recruits is high: recruitment decreases slowly (Figure 2b and Figure 3b) and recruits’ growth increases rapidly (Figure 2c and Figure 3c). In the same conditions of high disturbance intensity, survivors’ ACS growth is lower in the first years following logging than for low disturbance intensities, but declines slowly (Figure 2a and Figure 3a). Disturbance intensity and relative forest maturity have a weak effect on ACS loss from both survivors and recruits (Figures 2d,e and 3d,e). Overall, net ACS change stays high longer at high disturbance intensity (Figure 3f).

Figure 3

Download asset Open asset

Predicted effect of disturbance intensity on ACS changes along time in an Amazonian-average plot.

(a) Survivors’ ACS growth. (b) New recruits’ ACS. (c) Recruits’ ACS growth. (d) Survivors’ ACS loss. (e) Recruits’ ACS loss. (f) Net ACS change. The net ACS change is the sum of all five ACS changes. ACS changes were calculated with all parameters set to their maximum-likelihood value and covariates (except standardized disturbance intensity $l o s s$ ) set to 0. Time since minimum ACS varies from 0 to 30 year (i.e. the calibration interval) and disturbance intensity ranges between 5% and 60% of initial ACS loss.

https://doi.org/10.7554/eLife.21394.008

Regional variations of ACS changes

Variations of post-logging ACS changes between sites are explained with the mean ACS of each site ( $a c s 0$ ), climatic conditions [annual precipitation ( $p r e c$ ), seasonality of precipitation ( $s e a s$ )] and the soil bulk density ( $b d$ ). Contribution of survivors’ growth to ACS recovery declined slowly in sites with low $a c s 0$ and high water stress (low precipitation, high seasonality and high bulk density) (Figure 2a). Survivors’ ACS loss showed the opposite pattern (Figure 2d) except in apparent response to high seasonality of precipitation ( $s e a s$ ) that slowed the post-disturbance rates of decline of both ACS growth and loss. Despite slower recruits’ ACS growth in sites with high pre-logging ACS ( $a c s 0$ ), no other regional covariate had significant effects on recruits’ ACS changes (Figure 2b,c and e).

Prediction maps

While no significant environmental effects were detected for recruits’ ACS changes (Figures 2 and 4), the survivors showed a highly structured regional gradient: (i) ACS gain from survivors’ ACS growth is high in the west and in the Guiana Shield, but low in the south (Figure 4a), whereas (ii) survivors’ ACS loss is low in the south and in the Guiana Shield but high in the west (Figure 4d). To illustrate how these regional differences will be critical for future ACS across Amazonia, we developed a map of net ACS recovery over the first 10 years after a 40% ACS loss by integrating the sum of ACS change predictions through time (Figure 5). Across the region, net ACS recovery over the first ten years after a 40% ACS loss is predicted to be 17 $\pm$ 7 Mg C ha $^{- 1}$ , with higher values in the west and in the Guiana Shield (Figure 5a). The uncertainty in predictions was low to medium (coefficient of variation under 40%) in 82% of the mapped area, and high (coefficient of variation above 50%) in 5% of the mapped area (Figure 5b).

Figure 4

Download asset Open asset

Predicted cumulative ACS changes (Mg C ha $^{- 1}$ ) over the first 10 year after losing 40% of ACS.

Extrapolation was based on global rasters: topsoil bulk density from the Harmonized global soil database (Nachtergaele et al., 2008), Worldclim precipitation data (Hijmans et al., 2005) and biomass stocks from Avitabile et al. map (Avitabile et al., 2016). Cumulative ACS changes are obtained by integrating annual ACS changes through time. We here show the median of each pixel. Top graphs are ACS gain and bottom graphs are ACS loss. (a) ACS gain from survivors’ growth. (b) ACS gain from new recruits. (c) ACS gain from recruits’ growth. (d) ACS loss from survivors’ mortality. (e) ACS loss from recruits’ mortality. Black dots are the location of our experimental sites. Survivors’ ACS changes (a and d) show strong regional variations unlike to recruits’ ACS changes (b,c and e).

https://doi.org/10.7554/eLife.21394.009

Figure 5

Download asset Open asset

Predicted net ACS recovery over the first 10 year after losing 40% of pre-logging ACS.

(a) median predictions. (b) coefficient of variation (per pixel). Four areas were arbitrarily chosen to illustrate four different geographical behaviours: (1) the Guiana Shield and (2) northwestern Amazonia are two areas with high ACS recovery; the Guiana Shield has higher initial ACS and slower ACS dynamics whereas northwestern Amazonia has lower initial ACS and faster ACS dynamics. (3) central Amazonia has intermediate ACS recovery. (4) southern Amazonia has low ACS recovery.

https://doi.org/10.7554/eLife.21394.010

Four areas (Figure 5a) were selected to represent four contrasted cases of net ACS recovery in time (Figure 6): two areas, northwestern Amazonia and the Guiana Shield, with high ACS accumulation (21 $\pm$ 3 Mg C ha $^{- 1}$ over 10 year), one intermediate area, central Amazonia (15 $\pm$ 1 Mg C ha $^{- 1}$ over 10 year) and one area with low ACS accumulation, southern Amazonia (12 $\pm$ 3 Mg C ha $^{- 1}$ over 10 year). Survivors’ contribution to the sum of ACS gains (recruitment and growth) over the first 10 years after disturbance was 71 $\pm$ 4% in the Guiana Shield, 71 $\pm$ 2% in the west; 63 $\pm$ 4% in central Amazonia and 55 $\pm$ 6% in the south. Predicted net ACS recovery (Figure 5) and survivors’ ACS growth (Figure 4a) are highly correlated: $ρ = 0.90$ (Pearson’s correlation coefficient).

Figure 6

Download asset Open asset

Predicted contribution of annual ACS changes in ACS recovery in four regions of Amazonia (Figure 5).

The white line is the net annual ACS recovery, i.e. the sum of all annual ACS changes. Survivors’ (green) and recruits’ (orange) contribution are positive for ACS gains (survivors’ ACS growth, new recruits’ ACS and recruits’ ACS growth) and negative for survivors’ and recruits’ ACS loss. Areas with higher levels of transparency and dotted lines are out of the calibration period (0–30 year). In the Guiana Shield and in nothwestern Amazonia, high levels of net ACS recovery are explained by large ACS gain from survivors’ growth. Extrapolation was based on global rasters: topsoil bulk density from the Harmonized global soil database (Nachtergaele et al., 2008), precipitation data from Worldclim (Hijmans et al., 2005) and biomass stocks from Avitabile et al. (Avitabile et al., 2016) map.

https://doi.org/10.7554/eLife.21394.011

Discussion

Contrasting post-disturbance ACS dynamics were detected among the western Amazon, Guiana Shield, and southern Amazon (Figure 4). (i) In the western Amazon, environmental stress is reduced due to fertile soils and abundant, mostly non-seasonal precipitation, but forests are prone to frequent and sometimes large-scale wind-induced disturbances (Espírito-Santo et al., 2014). Such conditions of low stress and high disturbance tend to favor fast-growing species with rapid life cycles (He et al., 2013), which results in fast ACS gain and loss from survivors even after the logging disturbance (Figures 4a,d and 6). (ii) Forests of the Guiana Shield are generally dense and grow on nutrient-poor soils (Quesada et al., 2012), where wood productivity is highly constrained by competition for key nutrients, especially phosphorus and nitrogen (Santiago, 2015; Mercado et al., 2011). The short duration pulse of nutrients released from readily decomposed stems, twigs and leaves of trees damaged and killed by logging may thus explain the substantial but limited-duration increase in growth of survivors on these nutrient-poor soils (Figure 6). Yet post-disturbance ACS loss from survivors’ mortality decreases slowly in the Guiana Shield (Figure 6). This is consistent with the low mortality rates and the high tree longevity reported in old-growth forests of this region (Phillips et al., 2004). (iii) In the southern Amazon, high seasonal water stress is the main constraint on ACS recovery (Wagner et al., 2016). Stress-tolerant trees are generally poor competitors (He et al., 2013) and this may explain the slow ACS changes of survivors in this region (Figures 4a,d and 6). Finally, Central Amazonia is a transition zone for the main environmental and biotic gradients found in Amazonia: (1) a competition gradient between dense and nutrient-poor northeastern forests and nutrient-rich western forests; (2) an environmental gradient between northern wet forests and southern drier forests (Quesada et al., 2012).

Across Amazonia, survivors contribute most to post-disturbance ACS recovery. In regions where survivors’ ACS gain is high (west and northeast), net ACS recovery is also high: annual ACS recovery is between 1 and 3 Mg C ha $^{- 1}$ yr $^{- 1}$ in the first 10 year after logging (Figure 6), lower than in Amazonian secondary forests (3–5 Mg C ha $^{- 1}$ yr $^{- 1}$ in the first 20 year after abandonment of land use [Poorter et al., 2016]). Recruits, for their part, have very low geographical variations in post-logging ACS changes: 10 years after the disturbance they are predicted to store similar amounts of ACS almost everywhere in Amazonia. Nevertheless, small trees with DBH $<$ 20 cm have not been accounted for in our study and may play an important role in post-logging ACS changes. The 10–20 cm DBH size class contains as much as 14% of total ACS and may be highly dynamic in some Amazonian forests (Vieira et al., 2004). Because of the slow tree growth rates in Amazonia (Vieira et al., 2005; Herault et al., 2010), many trees will not reach the 20 cm DBH threshold 10 years after logging: the effects of the 10–20 cm DBH stratum on post-logging ACS changes are likely to be missed in sites with less than 10 years of measurements (e.g. Peteco, Ecosilva, Iracema, Cumaru) and should be studied, together with the natural regeneration, in the future.

At the stand level, high disturbance intensities reduce survivors’ ACS: survivors’ ACS growth is consequently lower (Figure 3a), resulting in lower net ACS change during the first 10 years of the recovery period (Figure 3f). High disturbance intensities as well as relatively low forest maturity alleviate competition, and this is probably why ACS contributions from recruits remain high for longer (Figure 2b) in such enhanced growth conditions (Herault et al., 2010). In the first years after logging, net ACS recovery depends little on disturbance intensity (Figure 3f), but recovery is predicted to last longer in heavily logged forests. In immature forests, intense self-thinning (Swaine et al., 1987) may explain fast ACS losses from survivors’ mortality (Figure 2d).

In the tropics, reduced-impact logging techniques (RIL; [Putz et al., 2008]) are promoted to reduce collateral damage to residual stands and biodiversity. Our results reveal that lower disturbance intensities, as a direct consequence of the employment of RIL techniques, could increase survivors’ ACS growth and slow down their ACS loss. Given that government specified minimum cutting cycles are short, e.g. 35 year in the Brazilian Amazon (Blaser et al., 2011), and that many commercial species are slow-growing and dense-wooded (Dauber et al., 2005; Wright et al., 2010), available timber stocks for the next cutting cycle will be comprised mostly of survivors. Attention should be taken to high harvest intensities and/or substantial incidental damage due to poor harvesting practices that diminish stocks of survivors, even if they promote recruitment. Most trees that recruit are fast-growing pioneers that are favored by disturbance but are vulnerable to water stress (Bonal et al., 2016) and competition (Valladares and Niinemets, 2008), and because their height is lower than in mature forests (Rutishauser et al., 2016), they might have reduced carbon sequestration potential. With ongoing climate change and increased frequencies and intensities of droughts in Amazonia (Malhi et al., 2008), betting on recruits to store C in forests disturbed by selective logging might thus be a risky gamble.

In this study, we focus on one type of disturbance: selective logging. Because of its economic value and implications for forest management, selective logging is a long-studied human disturbance in tropical forests, and the data gathered by the TmFO network are unique in terms of experiment duration and spatial extent. We nevertheless believe that our study gives clues on the regional differences in Amazonian forests response to large ACS losses induced by other disturbances (e.g. droughts, fire) that are expected to increase in frequency with ongoing global changes (Bonal et al., 2016).

Materials and methods

Site description

Request a detailed protocol

Our study includes data from thirteen long-term (8–30 year) experimental forest sites located in the Amazon Basin and the Guiana Shield (Figure 1—figure supplement 1). Sites meet the following criteria: (i) located in tropical forests with mean annual precipitation above 1000 mm; (ii) a total censused area above 1 ha; (iii) at least one pre-logging census and (iv) at least two post-logging censuses. For each site, we extracted annual precipitation and seasonality of precipitation data from WorldClim (RRID:SCR_010244) (Hijmans et al., 2005), topsoil bulk density data from the Harmonized World Soil database (Nachtergaele et al., 2008), and the synthetic climatic index from Chave et al. (Chave et al., 2014), using in all cases the highest resolution data available (30 arc-seconds). For one of our sites (La Chonta, see Figure 1—figure supplement 1), field measurements of precipitation (mean = 1580 mm yr $^{- 1}$ ) differed substantially from WorldClim data (1032 mm yr $^{- 1}$ ): in this particular case we used the measured value and adjusted the synthetic climatic index (E) in the allometric equation (Chave et al., 2014) accordingly. Sites' data is available at Dryad Digital Repository (Piponiot et al., 2016).

ACS computation

Request a detailed protocol

In all plots, diameter at breast height (DBH) of trees $>$ 20 cm DBH were measured, and trees were identified to the lowest taxonomic level: to the species level (75%) when possible, or to the genus level (15%); 10% of trees were not identified. To get the wood density, we applied the following standardized protocol to all sites: (i) trees identified to the species level were assigned the corresponding wood specific gravity value from the Global Wood Density Database (GWDD, doi:10.5061/dryad.234/1) (Zanne et al., 2009); (ii) trees identified to the genus level were assigned a genus-average wood density; (iii) trees with no botanical identification or that were not in the GWDD were assigned the site-average wood density. The aboveground biomass (AGB) was estimated with the allometric equations from Chave et al. (Chave et al., 2014). Biomass was assumed to be 50% carbon (Penman et al., 2003). The ACS of every tree $i$ was then computed as follows:

\begin{array}{l} \hat{A C S_{i}} = e x p (- 1.803 - 0.976 \times E + 0.976 \times l n (W D_{i}) + \\ + 2.673 \times l n (D B H_{i}) - 0.0299 \times l n (D B H_{i}^{2})) \times 0.5 \end{array}

where $W D_{i}$ and $D B H_{i}$ are the specific wood density and diameter at breast height of the tree $i$ and $E$ is the synthetic climatic index (Chave et al., 2014). The ACS changes data that was generated is available at Dryad Digital Repository (Piponiot et al., 2016).

The recovery period

Request a detailed protocol

After logging, plot ACS decreases rapidly until it reaches its minimum value ( $a c s m i n$ ) a few years later. This transition point determines the beginning $t_{m i n} = t_{0}$ of the recovery period. $a c s m i n$ was estimated as the minimum ACS in the 4 years following logging activities. Because our focus is on post-logging ACS recovery, we did not include in our analysis plots where the minimum ACS value was not reached within the 4 years after logging, either because the logging activity did not affect the plot or because there were other sources of disturbance long after logging (fire, road opening, silvicultural treatments).

ACS changes computation

Request a detailed protocol

For each plot $j$ and census $k$ , with $t_{k}$ the time since the beginning of the recovery period $t_{0}$ , we define 5 ACS changes : new recruits’ ACS ( $R r_{j, k}$ ) is the ACS of all trees $<$ 20 cm DBH at $t_{k - 1}$ and $\geq$ 20 cm DBH at $t_{k}$ ; recruits’ ACS growth ( $R g_{j, k}$ ) is the ACS increment of living recruits between $t_{k - 1}$ and $t_{k}$ ; recruits’ ACS loss ( $R l_{j, k}$ ) is the C in recruits that die between $t_{k - 1}$ and $t_{k}$ ; survivors’ ACS growth ( $S g_{j, k}$ ) is the ACS increment of living survivors between $t_{k - 1}$ and $t_{k}$ ; survivors’ ACS loss ( $S l_{j, k}$ ) is the ACS of survivors that die between $t_{k - 1}$ and $t_{k}$ . ACS gains ( $S g$ , $R r$ , $R g$ ) are positive and ACS losses ( $S l$ , $R l$ ) are negative. Instantaneous ACS changes are subject to stochastic variation over time: because we are less interested in year-to-year variations than in long-term ACS trajectories, we modelled cumulative ACS changes instead of annual ACS changes. Cumulative ACS changes (Mg C ha $^{- 1}$ ) were defined as follows:

c C h a n g e_{j, k} = \sum_{m = 0}^{k} C h a n g e_{j, m}

where $j$ is the plot, $t_{k}$ the time since $t_{0}$ (yr) and $C h a n g e$ is the annual ACS change (Mg C ha $^{- 1}$ yr $^{- 1}$ ), either recruits’ ACS ( $R r$ ), recruits’ ACS growth ( $R g$ ), recruits’ ACS loss ( $R l$ ), survivors’ ACS growth ( $S g$ ), or survivors’ ACS loss ( $S l$ ).

Covariates

Request a detailed protocol

To model ACS changes, we chose six covariates : (1) $l o s s$ disturbance intensity, i.e. percentage of initial ACS loss; (2) $a c s 0$ mean ACS of the site; (3) $d a c s$ relative ACS of the plot, as a % of $a c s 0$ ; (4) $p r e c$ annual precipitation; (5) $s e a s$ precipitation seasonality; (6) $b d$ topsoil bulk density. To give equivalent weight to all covariates, we centred and standardized them in order to have a mean of zero and a standard deviation of one over all observations. The uncertainty associated with ACS covariates ( $l o s s$ , $a c s 0$ , $d a c s$ ) is less than 10% (Chave et al., 2014). Climatic covariates (annual precipitation $p r e c$ and precipitation seasonality $s e a s$ ) were extracted from Worldclim rasters (RRID:SCR_010244). Error in Worldclim precipitation data was estimated to be $<$ 10 mm in Amazonia (Hijmans et al., 2005). There is no information on the uncertainty on topsoil bulk density but we expect it to be higher than the uncertainty on other covariates, due to measurement (De Vos et al., 2005) and interpolation methods (Hendriks et al., 2016).

Survivors’ model

Request a detailed protocol

Survivors’ cumulative ACS changes are null at $t = 0$ (by definition). When all survivors are dead, their ACS changes stop: annual ACS changes become null and cumulative ACS changes reach a constant/finite limit. We decided to model survivors’ cumulative ACS growth $c S g$ and ACS loss $c S l$ as:

c S_{i, j, k} \sim 𝒩 (α_{j}^{S} \times (1 - e x p (- β_{j}^{S} \times t_{k})), {(σ_{E}^{S})}^{2})

where $j$ is the plot, $t_{k}$ is the time since $t_{0},$ $S$ is either $S g$ or $S l . α_{j}^{S}$ is the finite limit of the cumulative ACS change and $β_{j}^{S}$ the rate at which the cumulative ACS change converges to this limit. By choosing an exponential kernel, we assume that survivors’ ACS change at $t_{k}$ is proportional to survivors’ ACS change at $t_{k} - 1$ .

Because $α_{j}^{S}$ values are expected to vary among plots, they are modelled with the following distribution:

α_{j}^{S} \sim 𝒩 (α_{0}^{S}, {(σ_{α}^{S})}^{2})

Parameter $β_{j}^{S}$ is the rate at which survivors’ ACS change (from growth or mortality) on plot $j$ converges to a finite limit after the disturbance: it reflects the response rapidity of survivors’ ACS changes to disturbance. Because we are interested in predicting variations in $β_{j}^{S}$ ( $S$ is either $S g$ or $S l$ ), we expressed $β_{j}^{S}$ as a function of covariates:

β_{j}^{S} = β_{0}^{S} + \sum_{l = 1}^{6} (λ_{l}^{S} \times V_{j, l})

where $\sum_{l = 1}^{6} (λ_{l}^{S} \times V_{j, l})$ , is the effect of covariates ( $V_{j, l}$ ) on the post-logging rate $β_{j}$ . Covariates are centred and standardized and are (1) $l o s s$ : disturbance intensity, i.e. percentage of initial ACS loss; (2) $a c s 0$ : mean ACS of the site; (3) $d a c s$ relative ACS of the plot, as a % of $a c s 0$ ; (4) $p r e c$ annual precipitation; (5) $s e a s$ precipitation seasonality; (6) $b d$ topsoil bulk density.

When all survivors in plot $j$ are dead, all the C gained by their growth ( $c S g_{j, \infty} = α_{j}^{S g}$ ) plus their initial ACS ( $a c s m i n_{j}$ ) will have been lost ( $c S l_{j, \infty} = α_{j}^{S l}$ ). We thus added the following constraint to each plot $j$ :

α_{j}^{S l} + α_{j}^{S g} + a c s m i n_{j} = 0

with $α_{j}^{S g}, α_{j}^{S l}$ the finite limits of survivors’ cumulative ACS growth and ACS loss respectively, and $a c s m i n_{j}$ the ACS of the plot $j$ at $t_{m i n} = t_{0}$ .

Recruits’ model

Request a detailed protocol

When survivors are all dead, recruits will constitute the new forest. We made the assumption that the ACS of this new forest will reach a dynamic equilibrium: recruits’ annual ACS changes are expected to converge to constant values (that are however prone to small inter-annual variations), with ACS gains compensating ACS losses. Because there are no recruits yet at $t_{0}$ , recruits’ annual ACS growth ( $R g$ ) and ACS loss ( $R l$ ) are zero, and progressively increase to reach their asymptotic values. Recruits’ annual ACS growth and ACS loss can be thus modelled with the function:

f (t; α, β) = α \times (1 - e x p (- β \times t))

where $t$ is the time since the beginning of the recovery period. In the same logic as survivors’ cumulative ACS change, $α$ is the asymptotic value of recruits’ annual ACS change (Mg C ha $^{- 1}$ yr $^{- 1}$ ), and $β$ is the rate at which this asymptotic value is reached.

Contrary to recruits’ annual ACS growth and ACS loss, the ACS of new recruits ( $R r$ ) is high at $t_{0}$ because of the competition drop induced by logging, but then progressively decreases to reach its asymptotic value. We modelled it with the following function:

f (t; α, β, η) = α \times (1 + η \times e x p (- β \times t))

where $t$ is the time since logging. The parameter $η$ was added to allow annual recruited ACS to be higher than $α$ at $t_{0}$ .

As stated before, we chose to model cumulative ACS changes instead of annual ACS changes. The general model for recruits’ cumulative ACS changes is deduced by integrating annual ACS changes from $t_{0}$ to $t_{k}$ :

c R_{i, j, k} \sim 𝒩 (α_{i}^{R} \times (t_{k} + η \times \frac{1 - e x p (- β_{j}^{R} \times t_{k})}{β_{j}^{R}}), (σ_{E}^{R})^{2})

where $i$ is the site, $j$ is the plot, $t_{k}$ is the time since $t_{0} R$ is either $R r, R g$ or $R l$ . When $R$ is $R g$ or $R l$ , $η = - 1$ ; when $R$ is $R r$ , $η > 0$ .

Once the forest reaches a new dynamic equilibrium, recruits’ annual ACS changes should depend mostly on each site’s characteristics: we expect there to be more inter-site than intra-site variation in recruits’ asymptotic ACS changes $α^{R}$ . This is why we use one value $α_{i}^{R}$ per site $i$ , and model it as follows:

α_{i}^{R} \sim 𝒩 (α_{0}^{R}; {(σ_{α}^{R})}^{2})

When the dynamic equilibrium is reached, annual ACS gain (growth and recruitment) compensates annual ACS loss (mortality). We thus added the following constraint for every site $i$ :

α_{i}^{R r} + α_{i}^{R g} + α_{i}^{R l} = 0

With the same logic as for survivors, we are interested in predicting variation in $β^{R}$ . Given that we use one value $α_{i}^{R}$ per site $i$ (i.e. all plots in one site $i$ have the same value for $α_{i}^{R}$ ), we chose to take into account the inter-plot variability as follows:

β_{j}^{R} \sim 𝒩 (β_{0}^{R} + \sum_{l = 1}^{6} (λ_{l}^{R} \times V_{j, l}), {(σ_{β}^{R})}^{2})

Inference

Request a detailed protocol

Bayesian hierarchical models were inferred through MCMC methods using an adaptive form of the Hamiltonian Monte Carlo sampling (Carpenter et al., 2015). Each observation was given a weight proportional to the size of the plot. Codes were developed using the R language (RRID:SCR_001905) (R Developement Core Team, 2015) and the Rstan package (Carpenter et al., 2015). A detailed list of priors is provided in Table 1.

Table 1

List of priors used to infer ACS changes in a Bayesian framework. Models are : ( $S g$ ) survivors’ ACS growth, ( $S l$ ) survivors’ ACS loss, ( $R r$ ) new recruits’ ACS, ( $R g$ ) recruits’ ACS growth, ( $R l$ ) recruits’ ACS loss. $λ_{l o s s}$ is the parameter relative to the covariate $l o s s$ (logging intensity).

https://doi.org/10.7554/eLife.21394.012

Model	Parameter	Prior	Justification
$S g$	$α_{j}^{S g}$	$𝒰 [25, 250]$	On average 100 survivors/ha storing 0.25 to 2.5 MgC each
$S g$	$β_{j}^{S g}$	$𝒰 [0.015, 0.04]$	${75 < t_{0.95}^{S g}}^{*} < 200$ yr
$S l$	$β_{j}^{S l}$	$𝒰 [0.006, β^{S g}]$	$t_{0.95}^{S g} < {t_{0.95}^{S l}}^{*} < 500$ yr
$R r$	$α_{i}^{R r}$	$𝒰 [0.1, 1]$	Range of observed values in TmFO control plots
$R r$	$β_{j}^{R r}$	$𝒰 [0.006, 0.6]$	$5 < {t_{0.95}^{R r}}^{*} < 500$ yr
$R r$	$η$	$𝒰 [0, 3]$	$R r (t = 0) < 3 \times R r (t = \infty)$
$R g$	$α_{i}^{R g}$	$𝒰 [0.5, 3]$	Range of observed values in Amazonia (Johnson et al., 2016)
$R g$	$β_{j}^{R g}$	$𝒰 [0.006, 0.15]$	$20 < {t_{0.95}^{R g}}^{*} < 500$ yr
$R l$	$β_{j}^{R l}$	$𝒰 [0.003, 0.06]$	$50 < {t_{0.95}^{R l}}^{*} < 1000$ yr
All models M^†	$λ_{l o s s}^{M}$	$𝒰 [- β^{M}, β^{M}]$	Avoid multicollinearity problems
All models M^†	${(λ_{l}^{M})}_{l \neq l o s s}$	$𝒰 [- \frac{β^{M}}{4}, \frac{β^{M}}{4}]$	Avoid multicollinearity problems

^∗ $t_{0.95} = \frac{l n (20)}{β}$ is the time when the ACS change has reached 95% of its asymptotic value.
^†M is one of the five models: either $S g$ , $S l$ , $R r$ , $R g$ , $R l$ .

Prediction maps

Request a detailed protocol

Maps were obtained with the following steps: (i) spatially-explicit covariates are extracted at the resolution of 30 arc-second from: the pan-tropical carbon map of Avitabile et al. for pre-disturbance aboveground carbon stocks (Avitabile et al., 2016); WorldClim (RRID:SCR_010244) (Hijmans et al., 2005) for annual precipitation and seasonality of precipitation, and the Harmonized World Soil database (Nachtergaele et al., 2008) for topsoil bulk density; (ii) disturbance intensity is set to 40% of pre-logging ACS loss, which is a common value for disturbance intensity after conventional logging in Amazonia (West et al., 2014; Blanc et al., 2009; Martin et al., 2015) , and the relative forest maturity $d a c s$ is set to zero; (iii) parameters are drawn from their previously calibrated distribution; (iv) to simulate random effects, all five parameters ( $α$ ) are taken from their distribution $𝒩 (α_{0}, σ_{α}^{2})$ ; (v) for every pixel, we estimate the five cumulative ACS changes ( $c S g$ , $c S l$ , $c R r$ , $c R g$ , $c R l$ ) 10 years after the 40% ACS loss, given the parameters value and the pixel covariates values extracted from global rasters. Steps (iii) to (v) are repeated 200 times and summary statistics are calculated for every pixel. Because a significant part of our sites have experiment duration lower than 10 years (Figure 1—figure supplement 1), we are less confident in Amazonian-wide predictions after that 10 year period. Maps were elaborated under the R statistical software (RRID:SCR_001905) (R Developement Core Team, 2015).

Appendix

The importance of a process-based approach

We here study C recovery dynamics with a demographic process-based approach, i.e. by segregating ACS changes into cohorts (survivors and recruits) and demographic processes (growth, recruitment, mortality), as opposed to an all-in-one model in which only the ecosystem net ACS change is modelled, without examination of demographic processes.

To compare the goodness of fit of the two approaches (all-in-one and process-based), we calibrated an all-in-one model with our data and compared the accuracy of its predictions with the process-based predictions reported in this study. The all-in-one model was written as follows:

C_{j, k} \sim 𝒩 ((a c s 0_{j} - a c s m i n_{j}) \times (1 - e x p (- β_{j}^{C} \times t_{k})), {(σ_{E}^{C})}^{2})

where $C_{j, k}$ (MgC ha $^{- 1}$ ) is the total C accumulation $t_{k}$ years after the disturbance in plot $j a c s 0_{j} - a c s m i n_{j}$ is the ACS lost by logging and $β_{j}^{C}$ is the rate at which the plot ACS returns to its pre-logging ACS. We took into account the effect of covariates and dependencies for $β_{j}^{C} :$

β_{j}^{C} \sim N (β_{0}^{C} + \sum_{l = 1}^{6} (λ_{l}^{C} \times V_{j, l}); {(σ_{β}^{C})}^{2})

The process-based model made better predictions (RMSE = 0.24) than the all-in-one model (RMSE = 0.31). In some sites, for example Paracou (black diamonds in Appendix—figure 1), there is a clear bias in the all-in-one model predictions: C accumulation is overestimated at the beginning of the recovery period and underestimated towards the end. This bias may be due to the non-adequacy of the negative exponential curve in the classic all-in-one model ( Appendix—figure 2a) to the C recovery observed in experimental plots (Rutishauser et al., 2015). The process-based model does not predict a constant instantaneous C accumulation rate (Appendix—figure 2b), and is thus more accurate.

Appendix 1—figure 1

Download asset Open asset

Observed vs fitted values of net ACS accumulation (MgC ha $^{- 1}$ ).

(a) Fitted values from the all-in-one model. (b) Fitted values from the process-based model (right). Net ACS accumulation is the sum of cumulative ACS changes (gain and loss). Each combination of a colour and shape is specific to a site. The closer the dots are to the x=y line, the better the prediction.

https://doi.org/10.7554/eLife.21394.013

Appendix 1—figure 2

Download asset Open asset

Predicted trajectories of net ACS accumulation (MgC ha $^{- 1}$ ) per site with (a) the all-in-one model and (b) the process-based model.
https://doi.org/10.7554/eLife.21394.014

Data availability

The following data sets were generated

1. Piponiot C
2. Sist P
3. Mazzei L
4. Peña-Claros M
5. Putz F
6. Rutishauser E
7. Shenkin A
8. Ascarrunz N
9. de Azevedo C
10. Baraloto C
11. França M
12. Guedes M
13. Honorio Coronado E
14. d'Oliveira MVN
15. Ruschel AR
16. da Silva KE
17. Doff Sotta E
18. de Souza CR
19. Vidal E
20. West TAP
21. Hérault B
(2016) Data from: Post-disturbance carbon recovery in Amazonian forests
Available at Dryad Digital Repository under a CC0 Public Domain Dedication.

http://dx.doi.org/10.5061/dryad.rc279

References

(2014) Monitoring selective logging in western Amazonia with repeat lidar flights
Remote Sensing of Environment 151:157–165.

https://doi.org/10.1016/j.rse.2013.08.049
- Google Scholar
1. Aragão LEOC
2. Malhi Y
3. Metcalfe DB
4. Silva-Espejo JE
5. Jiménez E
6. Navarrete D
7. Almeida S
8. Costa ACL
9. Salinas N
10. Phillips OL
11. Anderson LO
12. Alvarez E
13. Baker TR
14. Goncalvez PH
15. Huamán-Ovalle J
16. Mamani-Solórzano M
17. Meir P
18. Monteagudo A
19. Patiño S
20. Peñuela MC
21. Prieto A
22. Quesada CA
23. Rozas-Dávila A
24. Rudas A
25. Silva Jr. JA
26. Vásquez R
(2009) Above- and below-ground net primary productivity across ten Amazonian forests on contrasting soils
Biogeosciences 6:2759–2778.

https://doi.org/10.5194/bg-6-2759-2009
- Google Scholar
1. Asner GP
2. Broadbent EN
3. Oliveira PJ
4. Keller M
5. Knapp DE
6. Silva JN
(2006) Condition and fate of logged forests in the brazilian Amazon
PNAS 103:12947–12950.

https://doi.org/10.1073/pnas.0604093103
- PubMed
- Google Scholar
1. Asner GP
2. Knapp DE
3. Broadbent EN
4. Oliveira PJ
5. Keller M
6. Silva JN
(2005) Selective logging in the brazilian Amazon
Science 310:480–482.

https://doi.org/10.1126/science.1118051
- PubMed
- Google Scholar
1. Asner GP
2. Rudel TK
3. Aide TM
4. Defries R
5. Emerson R
(2009) A contemporary assessment of change in humid tropical forests
Conservation Biology 23:1386–1395.

https://doi.org/10.1111/j.1523-1739.2009.01333.x
- PubMed
- Google Scholar
1. Avitabile V
2. Herold M
3. Heuvelink GB
4. Lewis SL
5. Phillips OL
6. Asner GP
7. Armston J
8. Ashton PS
9. Banin L
10. Bayol N
11. Berry NJ
12. Boeckx P
13. de Jong BH
14. DeVries B
15. Girardin CA
16. Kearsley E
17. Lindsell JA
18. Lopez-Gonzalez G
19. Lucas R
20. Malhi Y
21. Morel A
22. Mitchard ET
23. Nagy L
24. Qie L
25. Quinones MJ
26. Ryan CM
27. Ferry SJ
28. Sunderland T
29. Laurin GV
30. Gatti RC
31. Valentini R
32. Verbeeck H
33. Wijaya A
34. Willcock S
(2016) An integrated pan-tropical biomass map using multiple reference datasets
Global Change Biology 22:1406–1420.

https://doi.org/10.1111/gcb.13139
- PubMed
- Google Scholar
1. Baker TR
2. Phillips OL
3. Malhi Y
4. Almeida S
5. Arroyo L
6. Di Fiore A
7. Erwin T
8. Killeen TJ
9. Laurance SG
10. Laurance WF
11. Lewis SL
12. Lloyd J
13. Monteagudo A
14. Neill DA
15. Patino S
16. Pitman NCA
17. M. Silva JN
18. Vasquez Martinez R
(2004) Variation in wood density determines spatial patterns inAmazonian forest biomass
Global Change Biology 10:545–562.

https://doi.org/10.1111/j.1365-2486.2004.00751.x
- Google Scholar
1. Blanc L
2. Echard M
3. Herault B
4. Bonal D
5. Marcon E
6. Chave J
7. Baraloto C
(2009) Dynamics of aboveground carbon stocks in a selectively logged tropical forest
Ecological Applications 19:1397–1404.

https://doi.org/10.1890/08-1572.1
- PubMed
- Google Scholar
1. Blaser J
2. Sarre A
3. Poore D
4. Johnson S
(2011)
Status of Tropical Forest Management

Status of Tropical Forest Management.
- Google Scholar
1. Bonal D
2. Burban B
3. Stahl C
4. Wagner F
5. Hérault B
(2016) The response of tropical rainforests to drought-lessons from recent research and future prospects
Annals of Forest Science 73:27–44.

https://doi.org/10.1007/s13595-015-0522-5
- PubMed
- Google Scholar
1. Carpenter B
2. Gelman A
3. Hoffman M
4. Lee D
5. Goodrich B
6. Betancourt M
7. Brubaker MA
8. Li P
9. Riddell A
(2015)
Stan: A probabilistic programming language

Journal of Statistical Software.
- Google Scholar
(2014) Improved allometric models to estimate the aboveground biomass of tropical trees
Global Change Biology 20:3177–3190.

https://doi.org/10.1111/gcb.12629
- PubMed
- Google Scholar
(2005) Sustainability of timber harvesting in bolivian tropical forests
Forest Ecology and Management 214:294–304.

https://doi.org/10.1016/j.foreco.2005.04.019
- Google Scholar
(2005) Predictive quality of pedotransfer functions for estimating bulk density of forest soils
Soil Science Society of America Journal 69:500.

https://doi.org/10.2136/sssaj2005.0500
- Google Scholar
1. Espírito-Santo FD
2. Gloor M
3. Keller M
4. Malhi Y
5. Saatchi S
6. Nelson B
7. Junior RC
8. Pereira C
9. Lloyd J
10. Frolking S
11. Palace M
12. Shimabukuro YE
13. Duarte V
14. Mendoza AM
15. López-González G
16. Baker TR
17. Feldpausch TR
18. Brienen RJ
19. Asner GP
20. Boyd DS
21. Phillips OL
(2014) Size and frequency of natural forest disturbances and the Amazon forest carbon balance
Nature Communications 5:3434.

https://doi.org/10.1038/ncomms4434
- PubMed
- Google Scholar
1. Frolking S
2. Palace MW
3. Clark DB
4. Chambers JQ
5. Shugart HH
6. Hurtt GC
(2009) Forest disturbance and recovery: A general review in the context of spaceborne remote sensing of impacts on aboveground biomass and canopy structure
Journal of Geophysical Research: Biogeosciences 114:n/a.

https://doi.org/10.1029/2008JG000911
- Google Scholar
1. Goetz SJ
2. Hansen M
3. Houghton RA
4. Walker W
5. Laporte N
6. Busch J
(2015) Measurement and monitoring needs, capabilities and potential for addressing reduced emissions from deforestation and forest degradation under REDD+
Environmental Research Letters 10:123001.

https://doi.org/10.1088/1748-9326/10/12/123001
- Google Scholar
1. Gourlet-Fleury S
2. Cornu G
3. Jésel S
4. Dessard H
5. Jourget J-G
6. Blanc L
7. Picard N
(2005) Using models to predict recovery and assess tree species vulnerability in logged tropical forests: A case study from French Guiana
Forest Ecology and Management 209:69–85.

https://doi.org/10.1016/j.foreco.2005.01.010
- Google Scholar
(2014) Perturbations in the carbon budget of the tropics
Global Change Biology 20:3238–3255.

https://doi.org/10.1111/gcb.12600
- PubMed
- Google Scholar
1. Griggs DJ
2. Noguer M
(2002) Climate change 2001: The scientific basis. contribution of working group I to the third assessment report of the intergovernmental panel on climate change
Weather 57:267–269.

https://doi.org/10.1256/004316502320517344
- Google Scholar
(2013) Global shifts towards positive species interactions with increasing environmental stress
Ecology Letters 16:695–706.

https://doi.org/10.1111/ele.12080
- PubMed
- Google Scholar
(2016) Exploring the challenges with soil data in regional land use analysis
Agricultural Systems 144:9–21.

https://doi.org/10.1016/j.agsy.2016.01.007
- Google Scholar
1. Hijmans RJ
2. Cameron SE
3. Parra JL
4. Jones PG
5. Jarvis A
(2005) Very high resolution interpolated climate surfaces for global land areas
International Journal of Climatology 25:1965–1978.

https://doi.org/10.1002/joc.1276
- Google Scholar
(2012) Carbon emissions from land use and land-cover change
Biogeosciences 9:5125–5142.

https://doi.org/10.5194/bg-9-5125-2012
- Google Scholar
1. Huang M
2. Asner GP
(2010) Long-term carbon loss and recovery following selective logging in Amazon forests
Global Biogeochemical Cycles 24:n/a.

https://doi.org/10.1029/2009GB003727
- Google Scholar
1. Herault B
2. Ouallet J
3. Blanc L
4. Wagner F
5. Baraloto C
(2010) Growth responses of neotropical trees to logging gaps
Journal of Applied Ecology 47:821–831.

https://doi.org/10.1111/j.1365-2664.2010.01826.x
- Google Scholar
1. Johnson MO
2. Galbraith D
3. Gloor E
4. De Deurwaerder H
5. Guimberteau M
6. Rammig A
7. Thonicke K
8. Verbeeck H
9. von Randow C
10. Monteagudo A
11. Phillips OL
12. Brienen RJW
13. Feldpausch TR
14. Lopez Gonzalez G
15. Fauset S
16. Quesada C
17. Christoffersen B
18. Ciais P
19. Gilvan S
20. Kruijt B
21. Meir P
22. Moorcroft P
23. Zhang K
24. Alvarez E
25. Alves de Oliveira A
26. Amaral I
27. Andrade A
28. Aragao L
29. Araujo-Murakami A
30. Arets E
31. Arroyo L
32. Aymard G
33. Baraloto C
34. Barroso J
35. Bonal D
36. Boot R
37. Camargo J
38. Chave J
39. Cogollo A
40. Cornejo FV
41. Costa LD
42. di Fiore A
43. Ferreira L
44. Higuchi N
45. Honorio E
46. Killeen TJ
47. Laurance SG
48. Laurance WF
49. Licona J
50. Lovejoy T
51. Malhi Y
52. Marimon B
53. Marimon BHJ
54. Matos DCL
55. Mendoza C
56. Neill D
57. Pardo G
58. Peña-Claros M
59. Pitman NCa
60. Poorter L
61. Prieto A
62. Ramirez-Angulo H
63. Roopsind A
64. Rudas A
65. Salomao RP
66. Silveira M
67. Stropp J
68. Ter Steege H
69. Terborgh J
70. Thomas R
71. Toledo M
72. Torres-Lezama A
73. van der Heijden GMF
74. Vasquez R
75. Vieira I
76. Vilanova E
77. Vos V
78. Baker TR
(2016) Variation in stem mortality rates determines patterns of aboveground biomaass inAmazonian forests: implications for dynamic global vegetation models
Global Change Biology 44:1–18.

https://doi.org/10.1111/gcb.13315
- Google Scholar
1. Malhi Y
2. Roberts JT
3. Betts RA
4. Killeen TJ
5. Li W
6. Nobre CA
(2008) Climate change, deforestation, and the fate of the Amazon
Science 319:169–172.

https://doi.org/10.1126/science.1146961
- PubMed
- Google Scholar
1. Malhi Y
2. Wright J
(2004) Spatial patterns and recent trends in the climate of tropical rainforest regions
Philosophical Transactions of the Royal Society B: Biological Sciences 359:311–329.

https://doi.org/10.1098/rstb.2003.1433
- PubMed
- Google Scholar
1. Martin PA
2. Newton AC
3. Pfeifer M
4. Khoo M
5. Bullock JM
(2015) Impacts of tropical selective logging on carbon storage and tree species richness: A meta-analysis
Forest Ecology and Management 356:224–233.

https://doi.org/10.1016/j.foreco.2015.07.010
- Google Scholar
1. Mercado LM
2. Patino S
3. Domingues TF
4. Fyllas NM
5. Weedon GP
6. Sitch S
7. Quesada CA
8. Phillips OL
9. Aragao LEOC
10. Malhi Y
11. Dolman AJ
12. Restrepo-Coupe N
13. Saleska SR
14. Baker TR
15. Almeida S
16. Higuchi N
17. Lloyd J
(2011) Variations in Amazon forest productivity correlated with foliar nutrients and modelled rates of photosynthetic carbon supply
Philosophical Transactions of the Royal Society B: Biological Sciences 366:3316–3329.

https://doi.org/10.1098/rstb.2011.0045
- Google Scholar
(2008)
Harmonized world soil database

Technical Report.
- Google Scholar
1. Pan Y
2. Birdsey RA
3. Fang J
4. Houghton R
5. Kauppi PE
6. Kurz WA
7. Phillips OL
8. Shvidenko A
9. Lewis SL
10. Canadell JG
11. Ciais P
12. Jackson RB
13. Pacala SW
14. McGuire AD
15. Piao S
16. Rautiainen A
17. Sitch S
18. Hayes D
(2011) A large and persistent carbon sink in the world's forests
Science 333:988–993.

https://doi.org/10.1126/science.1201609
- PubMed
- Google Scholar
Book
1. Penman J
2. Gytarsky M
3. Hiraishi T
4. Krug T
(2003)
Good Practice Guidance for Land Use, Land-Use Change and Forestry

The Intergovernmental Panel on Climate Change.
- Google Scholar
1. Phillips OL
2. Baker TR
3. Arroyo L
4. Higuchi N
5. Killeen TJ
6. Laurance WF
7. Lewis SL
8. Lloyd J
9. Malhi Y
10. Monteagudo A
11. Neill DA
12. Vargas PN
13. Silva JN
14. Terborgh J
15. Martínez RV
16. Alexiades M
17. Almeida S
18. Brown S
19. Chave J
20. Comiskey JA
21. Czimczik CI
22. Di Fiore A
23. Erwin T
24. Kuebler C
25. Laurance SG
26. Nascimento HE
27. Olivier J
28. Palacios W
29. Patiño S
30. Pitman NC
31. Quesada CA
32. Saldias M
33. Lezama AT
34. Vinceti B
(2004) Pattern and process in Amazon tree turnover, 1976-2001
Philosophical Transactions of the Royal Society B: Biological Sciences 359:381–407.

https://doi.org/10.1098/rstb.2003.1438
- PubMed
- Google Scholar
Data
1. Piponiot C
2. Sist P
3. Mazzei L
4. Peña-Claros M
5. Putz F
6. Rutishauser E
7. Shenkin A
8. Ascarrunz N
9. de Azevedo CP
10. Baraloto C
11. França M
12. Guedes M
13. Honorio Coronado E
14. d'Oliveira MVN
15. Ruschel AR
16. da Silva KE
17. Doff Sotta E
18. de Souza CR
19. Vidal E
20. West TAP
21. Hérault B
(authors) (2016) Data from: Post-disturbance carbon recovery in Amazonian forests
Dryad Digital Repository.

https://doi.org/10.5061/dryad.rc279
(2016) Biomass resilience of neotropical secondary forests
Nature 530:211–214.

https://doi.org/10.1038/nature16512
- PubMed
- Google Scholar
(2008) Reduced-impact logging: Challenges and opportunities
Forest Ecology and Management 256:1427–1433.

https://doi.org/10.1016/j.foreco.2008.03.036
- Google Scholar
1. Quesada CA
2. Phillips OL
3. Schwarz M
4. Czimczik CI
5. Baker TR
6. Patiño S
7. Fyllas NM
8. Hodnett MG
9. Herrera R
10. Almeida S
11. Alvarez Dávila E
12. Arneth A
13. Arroyo L
14. Chao KJ
15. Dezzeo N
16. Erwin T
17. di Fiore A
18. Higuchi N
19. Honorio Coronado E
20. Jimenez EM
21. Killeen T
22. Lezama AT
23. Lloyd G
24. López-González G
25. Luizão FJ
26. Malhi Y
27. Monteagudo A
28. Neill DA
29. Núñez Vargas P
30. Paiva R
31. Peacock J
32. Peñuela MC
33. Peña Cruz A
34. Pitman N
35. Priante Filho N
36. Prieto A
37. Ramírez H
38. Rudas A
39. Salomão R
40. Santos AJB
41. Schmerler J
42. Silva N
43. Silveira M
44. Vásquez R
45. Vieira I
46. Terborgh J
47. Lloyd J
48. Czimczik CI
49. Patiño S
50. Arneth A
51. Chao KJ
52. López-González G
53. Neill DA
54. Núñez Vargas P
55. Prieto A
56. Rudas A
57. Santos AJB
58. Lloyd J
(2012) Basin-wide variations in Amazon forest structure and function are mediated by both soils and climate
Biogeosciences 9:2203–2246.

https://doi.org/10.5194/bg-9-2203-2012
- Google Scholar
Book
1. R Developement Core Team
(2015)
R: A Language and Environment for Statistical Computing

Vienna, Austria: R Foundation for Statistical Computing.
- Google Scholar
1. Rutishauser E
2. Hérault B
3. Baraloto C
4. Blanc L
5. Descroix L
6. Sotta ED
7. Ferreira J
8. Kanashiro M
9. Mazzei L
10. d'Oliveira MV
11. de Oliveira LC
12. Peña-Claros M
13. Putz FE
14. Ruschel AR
15. Rodney K
16. Roopsind A
17. Shenkin A
18. da Silva KE
19. de Souza CR
20. Toledo M
21. Vidal E
22. West TA
23. Wortel V
24. Sist P
(2015) Rapid tree carbon stock recovery in managed Amazonian forests
Current Biology 25:R787–R788.

https://doi.org/10.1016/j.cub.2015.07.034
- PubMed
- Google Scholar
(2016) Tree height reduction after selective logging in a tropical forest
Biotropica 48:285–289.

https://doi.org/10.1111/btp.12326
- Google Scholar
1. Santiago LS
(2015) Nutrient limitation of eco-physiological processes in tropical trees
Trees 29:1291–1300.

https://doi.org/10.1007/s00468-015-1260-x
- Google Scholar
(2015) Fates of trees damaged by logging in amazonian Bolivia
Forest Ecology and Management 357:50–59.

https://doi.org/10.1016/j.foreco.2015.08.009
- Google Scholar
1. Sist P
2. Ferreira FN
(2007) Sustainability of reduced-impact logging in the Eastern Amazon
Forest Ecology and Management 243:199–209.

https://doi.org/10.1016/j.foreco.2007.02.014
- Google Scholar
1. Sist P
2. Rutishauser E
3. Peña-Claros M
4. Shenkin A
5. Hérault B
6. Blanc L
7. Baraloto C
8. Baya F
9. Benedet F
10. da Silva KE
11. Descroix L
12. Ferreira JN
13. Gourlet-Fleury S
14. Guedes MC
15. Bin Harun I
16. Jalonen R
17. Kanashiro M
18. Krisnawati H
19. Kshatriya M
20. Lincoln P
21. Mazzei L
22. Medjibé V
23. Nasi R
24. d'Oliveira MVN
25. de Oliveira LC
26. Picard N
27. Pietsch S
28. Pinard M
29. Priyadi H
30. Putz FE
31. Rodney K
32. Rossi V
33. Roopsind A
34. Ruschel AR
35. Shari NHZ
36. Rodrigues de Souza C
37. Susanty FH
38. Sotta ED
39. Toledo M
40. Vidal E
41. West TAP
42. Wortel V
43. Yamada T
(2015) The tropical managed forests observatory: a research network addressing the future of tropical logged forests
Applied Vegetation Science 18:171–174.

https://doi.org/10.1111/avsc.12125
- Google Scholar
(1987) The dynamics of tree populations in tropical forest: a review
Journal of Tropical Ecology 3:359–366.

https://doi.org/10.1017/S0266467400002339
- Google Scholar
1. ter Steege H
2. Pitman NC
3. Phillips OL
4. Chave J
5. Sabatier D
6. Duque A
7. Molino JF
8. Prévost MF
9. Spichiger R
10. Castellanos H
11. von Hildebrand P
12. Vásquez R
13. Steege HT
(2006) Continental-scale patterns of canopy tree composition and function across Amazonia
Nature 443:444–447.

https://doi.org/10.1038/nature05134
- PubMed
- Google Scholar
1. Valladares F
2. Niinemets Ülo
(2008) Shade tolerance, a Key Plant Feature of Complex Nature and Consequences
Annual Review of Ecology, Evolution, and Systematics 39:237–257.

https://doi.org/10.1146/annurev.ecolsys.39.110707.173506
- Google Scholar
1. Valle D
2. Phillips P
3. Vidal E
4. Schulze M
5. Grogan J
6. Sales M
7. van Gardingen P
(2007) Adaptation of a spatially explicit individual tree-based growth and yield model and long-term comparison between reduced-impact and conventional logging in eastern Amazonia, Brazil
Forest Ecology and Management 243:187–198.

https://doi.org/10.1016/j.foreco.2007.02.023
- Google Scholar
1. Vidal E
2. West TAP
3. Putz FE
(2016) Recovery of biomass and merchantable timber volumes twenty years after conventional and reduced-impact logging in amazonian Brazil
Forest Ecology and Management 376:1–8.

https://doi.org/10.1016/j.foreco.2016.06.003
- Google Scholar
1. Vieira S
2. de Camargo PB
3. Selhorst D
4. da Silva R
5. Hutyra L
6. Chambers JQ
7. Brown IF
8. Higuchi N
9. dos Santos J
10. Wofsy SC
11. Trumbore SE
12. Martinelli LA
(2004) Forest structure and carbon dynamics in amazonian tropical rain forests
Oecologia 140:468–479.

https://doi.org/10.1007/s00442-004-1598-z
- PubMed
- Google Scholar
(2005) Slow growth rates of amazonian trees: Consequences for carbon cycling
PNAS 102:18502–18507.

https://doi.org/10.1073/pnas.0505966102
- Google Scholar
1. Wagner FH
2. Hérault B
3. Bonal D
4. Stahl C
5. Anderson LO
6. Baker TR
7. Becker GS
8. Beeckman H
9. Boanerges Souza D
10. Botosso PC
11. Bowman DMJS
12. Bräuning A
13. Brede B
14. Brown FI
15. Camarero JJ
16. Camargo PB
17. Cardoso FCG
18. Carvalho FA
19. Castro W
20. Chagas RK
21. Chave J
22. Chidumayo EN
23. Clark DA
24. Costa FRC
25. Couralet C
26. da Silva Mauricio PH
27. Dalitz H
28. de Castro VR
29. de Freitas Milani JE
30. de Oliveira EC
31. de Souza Arruda L
32. Devineau J-L
33. Drew DM
34. Dünisch O
35. Durigan G
36. Elifuraha E
37. Fedele M
38. Ferreira Fedele L
39. Figueiredo Filho A
40. Finger CAG
41. Franco AC
42. Freitas Júnior JL
43. Galvão F
44. Gebrekirstos A
45. Gliniars R
46. Graça PMLdeA
47. Griffiths AD
48. Grogan J
49. Guan K
50. Homeier J
51. Kanieski MR
52. Kho LK
53. Koenig J
54. Kohler SV
55. Krepkowski J
56. Lemos-Filho JP
57. Lieberman D
58. Lieberman ME
59. Lisi CS
60. Longhi Santos T
61. López Ayala JL
62. Maeda EE
63. Malhi Y
64. Maria VRB
65. Marques MCM
66. Marques R
67. Maza Chamba H
68. Mbwambo L
69. Melgaço KLL
70. Mendivelso HA
71. Murphy BP
72. O'Brien JJ
73. Oberbauer SF
74. Okada N
75. Pélissier R
76. Prior LD
77. Roig FA
78. Ross M
79. Rossatto DR
80. Rossi V
81. Rowland L
82. Rutishauser E
83. Santana H
84. Schulze M
85. Selhorst D
86. Silva WR
87. Silveira M
88. Spannl S
89. Swaine MD
90. Toledo JJ
91. Toledo MM
92. Toledo M
93. Toma T
94. Tomazello Filho M
95. Valdez Hernández JI
96. Verbesselt J
97. Vieira SA
98. Vincent G
99. Volkmer de Castilho C
100. Volland F
101. Worbes M
102. Zanon MLB
103. Aragão LEOC
(2016) Climate seasonality limits leaf carbon assimilation and wood productivity in tropical forests
Biogeosciences 13:2537–2562.

https://doi.org/10.5194/bg-13-2537-2016
- Google Scholar
1. West TAP
2. Vidal E
3. Putz FE
(2014) Forest biomass recovery after conventional and reduced-impact logging in Amazonian Brazil
Forest Ecology and Management 314:59–63.

https://doi.org/10.1016/j.foreco.2013.11.022
- Google Scholar
1. Wright SJ
2. Kitajima K
3. Kraft NJ
4. Reich PB
5. Wright IJ
6. Bunker DE
7. Condit R
8. Dalling JW
9. Davies SJ
10. Díaz S
11. Engelbrecht BM
12. Harms KE
13. Hubbell SP
14. Marks CO
15. Ruiz-Jaen MC
16. Salvador CM
17. Zanne AE
(2010) Functional traits and the growth-mortality trade-off in tropical trees
Ecology 91:3664–3674.

https://doi.org/10.1890/09-2335.1
- PubMed
- Google Scholar
1. Zanne AE
2. Lopez-Gonzalez G
3. Coomes DAA
4. Ilic J
5. Jansen S
6. Lewis S
7. Miller RBB
8. Swenson NGG
9. Wiemann MCC
10. Chave J
(2009) Global wood density database
Technical Report 235:.

https://doi.org/10.5061/dryad.234/1
- Google Scholar

Article and author information

Author details

Camille Piponiot
1. Université de Guyane, UMR EcoFoG (Agroparistech, CNRS, Inra, Université des Antilles, Cirad), Kourou, French Guiana
2. Cirad, UMR EcoFoG (Agroparistech, CNRS, Inra, Université des Antilles, Université de Guyane), Kourou, French Guiana
3. CNRS, UMR EcoFoG (Agroparistech, Inra, Université des Antilles, Université de Guyane, Cirad), Kourou, French Guiana
4. Cirad, UR Forests and Societies, Montpellier, France
Contribution
CP, Conception and design, Analysis and interpretation of data, Drafting or revising the article

For correspondence
camille.piponiot@gmail.com

Competing interests
The authors declare that no competing interests exist.

"This ORCID iD identifies the author of this article:" 0000-0002-3473-1982
Plinio Sist

Cirad, UR Forests and Societies, Montpellier, France

Contribution
PS, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article

Competing interests
The authors declare that no competing interests exist.
Lucas Mazzei

Embrapa Amazônia Oriental, Belém, Brazil

Contribution
LM, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article

Competing interests
The authors declare that no competing interests exist.
Marielos Peña-Claros

Forest Ecology and Forest Management Group, Wageningen University, Wageningen, Netherlands

Contribution
MP-C, Conception and design, Acquisition of data, Drafting or revising the article

Competing interests
The authors declare that no competing interests exist.
Francis E Putz

Department of Biology, University of Florida, Gainesville, United States

Contribution
FEP, Analysis and interpretation of data, Drafting or revising the article

Competing interests
The authors declare that no competing interests exist.
Ervan Rutishauser

CarbonForExpert, Hermance, Switzerland

Contribution
ER, Conception and design, Analysis and interpretation of data, Drafting or revising the article

Competing interests
The authors declare that no competing interests exist.
Alexander Shenkin

Environmental Change Institute, University of Oxford, Oxford, United Kingdom

Contribution
AS, Analysis and interpretation of data, Drafting or revising the article

Competing interests
The authors declare that no competing interests exist.
Nataly Ascarrunz

Instituto Boliviano de Investigación Forestal, Santa Cruz, Bolivia

Contribution
NA, Acquisition of data, Drafting or revising the article

Competing interests
The authors declare that no competing interests exist.
Celso P de Azevedo

Embrapa Amazônia Ocidental, Belém, Brazil

Contribution
CPdA, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article

Competing interests
The authors declare that no competing interests exist.
Christopher Baraloto

Department of Biological Sciences, International Center for Tropical Botany, Florida International University, Miami, United States

Contribution
CB, Analysis and interpretation of data, Drafting or revising the article

Competing interests
The authors declare that no competing interests exist.
Mabiane França

Embrapa Amazônia Ocidental, Belém, Brazil

Contribution
MF, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article

Competing interests
The authors declare that no competing interests exist.
Marcelino Guedes

Embrapa Amapa, Macapa, Brazil

Contribution
MG, Acquisition of data, Drafting or revising the article

Competing interests
The authors declare that no competing interests exist.
Eurídice N Honorio Coronado

Instituto de Investigaciones de la Amazonia Peruana, Iquitos, Peru

Contribution
ENHC, Acquisition of data, Drafting or revising the article

Competing interests
The authors declare that no competing interests exist.

"This ORCID iD identifies the author of this article:" 0000-0003-2314-590X
Marcus VN d'Oliveira

Embrapa Acre, Rio Branco, Brazil

Contribution
MVNd'O, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article

Competing interests
The authors declare that no competing interests exist.
Ademir R Ruschel

Embrapa Amazônia Oriental, Belém, Brazil

Contribution
ARR, Acquisition of data, Drafting or revising the article

Competing interests
The authors declare that no competing interests exist.
Kátia E da Silva

Embrapa Amazônia Ocidental, Belém, Brazil

Contribution
KEdS, Acquisition of data, Drafting or revising the article

Competing interests
The authors declare that no competing interests exist.
Eleneide Doff Sotta

Embrapa Amapa, Macapa, Brazil

Contribution
EDS, Acquisition of data, Drafting or revising the article

Competing interests
The authors declare that no competing interests exist.
Cintia R de Souza

Embrapa Amazônia Ocidental, Belém, Brazil

Contribution
CRdS, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article

Competing interests
The authors declare that no competing interests exist.
Edson Vidal

Departamento de Ciências Florestais, University of São Paulo, Piracicaba, Brazil

Contribution
EV, Acquisition of data, Drafting or revising the article

Competing interests
The authors declare that no competing interests exist.
Thales AP West

Department of Biology, University of Florida, Gainesville, United States

Contribution
TAPW, Acquisition of data, Drafting or revising the article

Competing interests
The authors declare that no competing interests exist.
Bruno Hérault

Cirad, UMR EcoFoG (Agroparistech, CNRS, Inra, Université des Antilles, Université de Guyane), Kourou, French Guiana

Contribution
BH, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article

For correspondence
bruno.herault@cirad.fr

Competing interests
The authors declare that no competing interests exist.

"This ORCID iD identifies the author of this article:" 0000-0002-6950-7286

Funding

Agence Nationale de la Recherche (ANR-10-LABEX-0025)

Camille Piponiot
Bruno Hérault

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP: 2013/16262-4 and 2013/50718-5)

Edson Vidal

European Regional Development Fund (FEDER 2014-2020, GY0006894)

Camille Piponiot
Bruno Hérault

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

We are in debt with all technicians and colleagues who helped setting up the plots and collecting data over years. Without their precious work, this study would have not been possible and they may be warmly thanked here. We are grateful to CIRAD, the GFclim project (FEDER 2014–2020, Project GY0006894) and the Sao Paulo Research Foundation (FAPESP: 2013/16262–4 and 2013/50718–5) for financial support. This study was partially funded by an Investissement d’Avenir grant of the ANR (CEBA: ANR-10-LABEX-0025) and carried out in the framework of the Tropical managed Forests Observatory (TmFO), supported by the Sentinel Landscape program of CGIAR (Consultative Group on International Agricultural Research) - Forest Tree and Agroforestry Research Program.

Copyright

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

3,061

views
586

downloads
47

citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Mendeley

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

Camille Piponiot
Plinio Sist
Lucas Mazzei
Marielos Peña-Claros
Francis E Putz
Ervan Rutishauser
Alexander Shenkin
Nataly Ascarrunz
Celso P de Azevedo
Christopher Baraloto
Mabiane França
Marcelino Guedes
Eurídice N Honorio Coronado
Marcus VN d'Oliveira
Ademir R Ruschel
Kátia E da Silva
Eleneide Doff Sotta
Cintia R de Souza
Edson Vidal
Thales AP West
Bruno Hérault

(2016)

Carbon recovery dynamics following disturbance by selective logging in Amazonian forests

eLife 5:e21394.

https://doi.org/10.7554/eLife.21394

Share this article

Cite this article

Post-disturbance annual ACS changes of survivors and recruits in 133 Amazonian selectively logged plots.

Effect of covariates on the rate at which post-disturbance ACS changes converge to a theoretical steady state (in yr-1).

Figure 2—source data 1

Predicted effect of disturbance intensity on ACS changes along time in an Amazonian-average plot.

Predicted cumulative ACS changes (Mg C ha−1) over the first 10 year after losing 40% of ACS.

Predicted net ACS recovery over the first 10 year after losing 40% of pre-logging ACS.

Predicted contribution of annual ACS changes in ACS recovery in four regions of Amazonia (Figure 5).

Observed vs fitted values of net ACS accumulation (MgC ha-1).

Predicted trajectories of net ACS accumulation (MgC ha-1) per site with (a) the all-in-one model and (b) the process-based model.

Author details

Camille Piponiot

Contribution

For correspondence

Competing interests

Plinio Sist

Contribution

Competing interests

Lucas Mazzei

Contribution

Competing interests

Marielos Peña-Claros

Contribution

Competing interests

Francis E Putz

Contribution

Competing interests

Ervan Rutishauser

Contribution

Competing interests

Alexander Shenkin

Contribution

Competing interests

Nataly Ascarrunz

Contribution

Competing interests

Celso P de Azevedo

Contribution

Competing interests

Christopher Baraloto

Contribution

Competing interests

Mabiane França

Contribution

Competing interests

Marcelino Guedes

Contribution

Competing interests

Eurídice N Honorio Coronado

Contribution

Competing interests

Marcus VN d'Oliveira

Contribution

Competing interests

Ademir R Ruschel

Contribution

Competing interests

Kátia E da Silva

Contribution

Competing interests

Eleneide Doff Sotta

Contribution

Competing interests

Cintia R de Souza

Contribution

Competing interests

Edson Vidal

Contribution

Competing interests

Thales AP West

Contribution

Competing interests

Bruno Hérault

Contribution

For correspondence

Competing interests

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

Effect of covariates on the rate at which post-disturbance ACS changes converge to a theoretical steady state (in yr $^{- 1}$ ).

Predicted cumulative ACS changes (Mg C ha $^{- 1}$ ) over the first 10 year after losing 40% of ACS.

Observed vs fitted values of net ACS accumulation (MgC ha $^{- 1}$ ).

Predicted trajectories of net ACS accumulation (MgC ha $^{- 1}$ ) per site with (a) the all-in-one model and (b) the process-based model.