The ESA-CCI AGB product was only released in late 2019 and was in its early phases of development at the time of use. However, given that its spatial resolution was high enough to separate areas of only secondary forest and its recent acquisition warranted its use for this research. Only areas of secondary forest greater than 9000 m2 were considered for further analysis, an area approximately equal to 1 pixel of the ESA-CCI product. Despite limiting the study to these larger secondary forest polygons, we were still left with just under 2.5 million polygons of secondary forest to analyse. The secondary forest map was laid over the AGC data and the modal AGC was extracted for each secondary forest polygon using the “zonal_stats” function available in the “rasterstats” module for the programming language “Python” (v3.6). We then aggregated the AGC values by the age of secondary forest and used the median AGC value for each age in further analysis. We applied a bias correction to the median AGC values, subtracting the smallest median value from all values to shift the data to begin at or near 0 Mg C ha−1 AGC for a 1-year-old secondary forest.
Following this, we used six remote sensing products of driving variables widely accepted to influence regrowth of forests. The data products included four environmental drivers (1–4) and two anthropogenic disturbance drivers (5–6): (1) Mean annual downward shortwave radiation (for the period 1985–2017)26 (link), (2) Mean annual precipitation (for the period 1985–2017)27 , (3) the mean Maximum Cumulative Water Deficit (MCWD) (for the period 1985–2017)65 (link),66 , (4) Soil Cation Concentration30 (link), (5) Annual burned areas (between 2001 and 2017)31 and (6) Number of times a secondary forest area was deforested between 1987 and 2017 (repeated deforestations) (this study). These products all have different spatial resolutions (Supplementary Table
The drivers were then grouped according to numerical limits, such as the 25, 50 and 75th percentiles. We then modelled the AGC for the age of secondary forest under these groupings using the commonly used Chapman-Richard model for regrowth67 (link):
where Yt refers to the AGC at age t; A is the AGC asymptote or the AGC of the old-growth forest; k is a growth rate coefficient of Y as a function of age; c is a coefficient that determines the shape of the growth curve; and ε is an error term. We applied the “nls” function available in the “nlstools” package for the statistical software R (v4.0.2)68 ,69 . We assumed that after a given amount of time, the AGC could return to levels equivalent to old-growth forests, and reach a pre-calculated asymptote. As such, we extracted the median, bias-corrected AGC value of old-growth forests under each variable condition from the ESA-CCI AGC product to represent the value of the asymptote (Supplementary Fig.