“…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 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 ; (v) for every pixel, we estimate the five cumulative ACS changes (, , ,,) 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.…”