Quantifying and understanding carbon storage and sequestration within the Eastern Arc Mountains of Tanzania, a tropical biodiversity hotspot
Carbon Balance and Management volume 9, Article number: 2 (2014)
The carbon stored in vegetation varies across tropical landscapes due to a complex mix of climatic and edaphic variables, as well as direct human interventions such as deforestation and forest degradation. Mapping and monitoring this variation is essential if policy developments such as REDD+ (Reducing Emissions from Deforestation and Forest Degradation) are to be known to have succeeded or failed.
We produce a map of carbon storage across the watershed of the Tanzanian Eastern Arc Mountains (33.9 million ha) using 1,611 forest inventory plots, and correlations with associated climate, soil and disturbance data. As expected, tropical forest stores more carbon per hectare (182 Mg C ha-1) than woody savanna (51 Mg C ha-1). However, woody savanna is the largest aggregate carbon store, with 0.49 Pg C over 9.6 million ha. We estimate the whole landscape stores 1.3 Pg C, significantly higher than most previous estimates for the region. The 95% Confidence Interval for this method (0.9 to 3.2 Pg C) is larger than simpler look-up table methods (1.5 to 1.6 Pg C), suggesting simpler methods may underestimate uncertainty. Using a small number of inventory plots with two censuses (n = 43) to assess changes in carbon storage, and applying the same mapping procedures, we found that carbon storage in the tree-dominated ecosystems has decreased, though not significantly, at a mean rate of 1.47 Mg C ha-1 yr-1 (c. 2% of the stocks of carbon per year).
The most influential variables on carbon storage in the region are anthropogenic, particularly historical logging, as noted by the largest coefficient of explanatory variable on the response variable. Of the non-anthropogenic factors, a negative correlation with air temperature and a positive correlation with water availability dominate, having smaller p-values than historical logging but also smaller influence. High carbon storage is typically found far from the commercial capital, in locations with a low monthly temperature range, without a strong dry season, and in areas that have not suffered from historical logging. The results imply that policy interventions could retain carbon stored in vegetation and likely successfully slow or reverse carbon emissions.
Tropical forests are globally significant ecosystems; accounting for ~50% of global forest area , storing ~ 45% of all carbon in terrestrial vegetation [2–4], maintaining high biodiversity , and providing ecosystem services, such as timber, non-timber forest products , and climate change mitigation [7, 8]. However, within the last few decades, vast areas of tropical forests have been converted to other land-uses or degraded. For example, between 1990 and 1997, 4.4-7.2 million hectares of humid tropical forest were converted each year and an additional 1.6-3.0 million hectares of forest were visibly degraded . This process increased in the early 2000s, with an estimated 5.1-5.7 million hectares of humid tropical forest (and 3.5-4.7 million hectares of dry tropical forest) deforested per year between 2000 and 2005 . The gradual and sustained reduction in forest quality and quantity has resulted in substantial emissions of CO2 . Globally, deforestation and forest degradation accounted for 6-20% of anthropogenic GHG emissions in the 1990s and early 2000s [12–14]. Tropical regions make a substantial contribution to this, emitting 0.7-1.5 Pg C yr-1 between 1990 and 1999 [9, 15–17] and 0.71.5 Pg C yr-1 between 2000 and 2007 [13, 16–18]. These processes also impact the future potential of forests to remove carbon from the atmosphere [7, 19, 20].
Recently, attempts to mitigate increasing anthropogenic CO2 emissions through reducing emissions from degradation and deforestation (REDD+) have been instigated . The REDD+ programme is aimed at contributing to a reduction in greenhouse emissions whilst providing economic incentives for better management and protection of forests. This policy has been widely welcomed and may provide a financial incentive to significantly reduce carbon emissions [22, 23], although the equity and justice issues surrounding the impact on local livelihoods are actively debated [24, 25]. Key technical issues for the successful implementation of REDD+ include (but are not limited to) the accuracy of monitoring systems, preventing leakage and establishing accurate historical baselines. Thus, the success of REDD+, in part, rests on robust scientific information on the magnitude and extent of carbon storage in tropical regions and how it changes over time .
The Intergovernmental Panel on Climate Change (IPCC) provide a three “Tier” system through which carbon stocks and emissions can be reported, each with a different level of methodological complexity and accuracy. Tier 1 is the simplest method, using global default values obtained from the IPCC literature [27, 28]. The intermediate Tier 2 level improves on Tier 1 by using country specific data. Tier 3 is the most rigorous approach, using local forest inventory data, focusing on the direct measurement of trees, repeated over a time series [27–29]. Here we develop a Tier 3 methodology for the Eastern Arc Mountains (EAM) watershed area.
The estimates become progressively more robust from Tier 1 to 3 due to changes in two main systematic errors . The first, completeness, refers to the number of IPCC carbon pools that are included, with studies including all five pools (aboveground live, litter, coarse wood debris [CWD], belowground and soil carbon) considered complete. The second, representativeness, derives from the substantial natural variability in the carbon stored across landscapes, even within a biome or country . The aboveground biomass of a forest within a landscape may differ considerably from global default (Tier 1) values or even from country-specific (Tier 2) values. For example, in the Peruvian Amazon, data from the Los Amigos Conservation Concession  were shown not to be representative of forests nationally. Nearby forests situated to the north and south of this local study are estimated to contain 20-35% less carbon per unit area , suggesting that Los Amigos Conservation Concession is an area of locally high biomass. Since Tier 3 methods account for variation observed within biomes and countries, the representativeness of the carbon estimates is higher than those associated with Tier 1 and 2 methodologies [32, 33].
However, Tier 3 methods are more expensive [34, 35] and some nations may lack the capacity to adopt such methods . Whilst, in some cases, the capability to apply Tier 3 guidelines is being rapidly developed, multi-temporal inventory data and data on historical carbon stock changes can take several decades to accrue [37, 38]. It is expected that REDD+ requirements will allow data provisions from several tiers in a single report. Highly variable and/or substantial carbon pools should be estimated using Tier 3 methodology (e.g. forest aboveground live carbon [ALC]), whilst Tier 1 or Tier 2 methodology may be sufficient for smaller carbon pools (e.g. CWD) or carbon poor land cover categories (e.g. bare ground).
In Tier 3 methods, in order to extrapolate from plot data, it is necessary to develop correlations with remotely sensed data to scale to the study area or country-wide estimates. Generally, carbon storage is either estimated via statistical correlation with electromagnetic properties, ground-truthed by volumetric measurements, such as diameter at breast height (DBH), which are converted to biomass estimates using allometric equations. A variety of remotely sensed data sources have been employed for carbon mapping and these can be aggregated into four groups: photographic imagery, RADAR, LiDAR, and ancillary geographic information systems (GIS) data (see Additional file 1: SI1 for an evaluation of each method). Here, we use ancillary GIS data as such data have three main advantages: 1) wide availability, often free of charge; 2) a suitable resolution (e.g. 90 m ); and 3) correlations with these ancillary GIS data may indicate which variables directly affect carbon storage. Developing an understanding of how these variables influence carbon storage is vital for accurate scenarios of future emissions.
Here, we correlate carbon storage estimates from tree inventory plots (n = 1,611, median size = 0.1 ha) with data on climatic (e.g. temperature, precipitation, and solar radiation), edaphic (e.g. soil water holding capacity and soil fertility) and proxy variables for direct human interventions (e.g. governance type, distance from the main economic demand centres, population pressure, and historical logging), and variables that derive from climate-human interactions (e.g. burnt area index) for the Tanzanian watershed of the Eastern Arc Mountains (hereafter, EAM ), which covers 33.9 million ha (Figure 1; see Swetnam et al (2011)  for further details). We develop Tier 3 type correlation equations to estimate the total ALC stored across the forested and wooded land cover categories, an advancement on previous Tier 2 estimates for the region presented in Willcock et al (2012) . Additionally, we investigate the most influential correlates of spatial differences in carbon storage and how these result from changes in either species composition affecting wood density (specific gravity) or the number of large trees present. Lastly, a smaller number of inventory plots (n = 43, median size 0.1 ha) have two censuses, and by applying the same mapping procedures, we assess changes in carbon storage over time, providing a first-order estimate of sequestration across the region.
Utilising 1,611 plots and scaling to the 33.9 million ha study area we estimate that 1.32 (95% confidence interval [CI] ranges from 0.89 to 3.16) Pg C was stored in the aboveground live vegetation in the year 2000 (Figure 2; Table 1). Woodland and bushland contributed most to the amount of stored aboveground live carbon (ALC) in the study region, with open woodland storing the most ALC (0.49 [0.47 to 1.60] Pg C over 9.6 million ha); followed by bushland (0.29 [0.15 to 0.51] Pg C over 5.0 million ha) and closed woodland (0.18 [0.13 to 0.61] Pg C over 1.8 million ha).
Best estimate values from our methodology, per unit area, in each land cover class, are given in Table 2. Forest contained the greatest ALC per unit area, with highest values in sub-montane forest (189 [95 to 588] Mg ha-1), followed by lowland (182 [152- to 360] Mg ha-1), upper montane (166 [69 to 533] Mg ha-1), montane (130 [62 to 702] Mg ha-1), and forest mosaic (121 [55 to 485] Mg ha-1). Woodlands held less ALC than forests, with closed woodland storing 100 (70 to 331) Mg ha-1 and open woodland storing 51 (38 to 165) Mg ha-1 (Table 2), but more than the landscape average of 39 (26 to 93) Mg ha-1.
Our sequestration model suggests that the landscape may be losing 0.05 (-0.07 to 0.26) Pg C yr-1 (mean net flux to atmosphere of 1.47 [-2.13 to 7.75] Mg C ha-1 yr-1). Of the 12.3 million ha of tree-dominated land in our study area, only 1.4% (0.17 million ha) shows a carbon decrease over the entire 95% CI range and only 0.8% (0.10 million ha) a definite carbon increase (Figure 3). The locations showing net carbon uptake are in the Udzungwa mountains, while the locations with net reductions in carbon storage are mainly in the Pare and Usambara mountains.
Links between carbon stock and influential variables
The variables that influence carbon storage and sequestration may be inferred from relationships within the correlation models. Forward selection results are presented in the following paragraphs as these best indicate causal relationships [43–45]. In general, backward models were in close agreement with forward models (Tables 3 and 4; Additional file 1: Tables S1-S3).
Carbon storage (adjusted R-squared [Adj R-sq] = 0.18) is correlated positively with the natural logarithm of the population pressure with decay constant of 12.5 km (p-value < 0.001) and increased by 1 Mg ha-1 for every 8700 km from a road (p-value < 0.010), and every 30,000 units in the cost distance to Dar es Salaam (p-value < 0.010). Carbon storage decreased by 1 Mg ha-1 for every 1°C increase in mean annual monthly temperature range (p-value < 0.001), every 2.7% rise in the total available water capacity of the soil (p-value < 0.001), and every 4.4 month increase in the mean number of dry months annually (p-value < 0.050). Carbon storage was 2.1 Mg ha-1 lower in areas where historical logging was present (p-value < 0.010), and 4.2 Mg ha-1 higher in areas under the control of local communities/governments (p-value < 0.010). Thus, carbon storage is high in areas far from the commercial capital, with a low monthly temperature range, without a dry season, that have not suffered from historical logging and are under local community/government control (Figure 4; Table 3).
The rate of carbon sequestration correlated with three principal component (PC) axes (presented in order of influence; Adj R-sq = 0.41). Carbon sequestration was negatively correlated with the soil fertility axis (PC5; p-value < 0.050), warmer temperatures and longer dry seasons (PC3; p-value < 0.050), and with increased anthropogenic disturbance (PC1; p-value < 0.010). Thus, carbon sequestration was highest in less fertile areas with little or no drought and little anthropogenic disturbance (Table 4).
Wood specific gravity (WSG; Adj R-sq = 0.28; see Additional file 1: SI2) was most strongly affected by the annual mean burned area probability (increasing by 1 g cm-3 for every 0.04 increase; p-value < 0.001) and the total available water capacity of the soil (decreasing by 1 g cm-3 for every 82.0% increase; p-value < 0.001). Thus, WSG is higher in burnt areas with little available water (Additional file 2: Figure S1; Additional file 3: Figure S2; Additional file 1: Table S1).
The intercept of the power law relationship (an indication of potential stem density [see Additional file 1: SI3]; Adj R-sq = 0.30) was most affected by the natural logarithm of the population pressure with decay constant of 12.5 km (positive correlation; p-value < 0.001) and the mean annual monthly temperature range (increasing by 1.0 for every 1.2°C increase; p-value < 0.001). Thus, the density of smaller stems increases in areas with a high population pressure and large temperature fluctuations (Additional file 3: Figure S2; Additional file 4: Figure S3; Additional file 1: Table S2).
Correlations identified for the gradient of the power law relationship (an indication of the proportion of larger stems; see Additional file 1: SI3) were broadly the inverse of those identified for the intercept. The gradient of the power law relationship was most affected by the natural logarithm of the population pressure with decay constant of 20.8 km (negative correlation; p-value < 0.001) and the mean burned area probability in the fourth quarter (decreasing by 1.0 for every 0.2 increase; p-value < 0.001). Thus, the proportion of large stems was greater in areas experiencing few disturbances from people or fire (Additional file 3: Figure S2; Additional file 5: Figure S4; Additional file 1: Table S3).
When investigating the most influential correlates of spatial differences in carbon storage and how these result from changes in either species composition affecting wood density (specific gravity) or the number of large trees present, we found that the final Tier 3 carbon storage estimates were positively correlated with both size-frequency distribution estimates (both intercept and gradient [p-values < 0.001]), and negatively correlated with WSG estimates (p- value < 0.001) and maximum height estimates (p-value < 0.001; Additional file 1: see SI4). All possible interactions were investigated and were significant (Adj R-sq = 0.35; p- values < 0.001), however, the majority of the explanatory power lay within the second order interactions (Adj R-sq = 0.33; p-values < 0.001; Additional file 1: Table S5). Broadly, WSG and the proportion of larger stems had largest influence over the carbon storage estimate. Considering only second order interactions, in areas of low potential stem density, carbon storage is positively correlated with maximum canopy height (Additional file 6: Figure S5). However, the opposite correlation is observed in areas of higher stem density. Although similar interactions are observed between both size-frequency distribution estimates (gradient and intercept), the interaction between WSG and maximum canopy height is inverse, with carbon storage only showing positive correlations with maximum canopy height in areas of high WSG. Both size-frequency distribution estimates also interacted similarly with WSG, with both showing positive correlations with carbon storage in areas of low WSG, but negative correlations in areas of high WSG (Additional file 6: Figure S5). Finally carbon sequestration correlation values were positively correlated with carbon storage estimates (p-value < 0.001), indicating that areas storing the most carbon are also those that are increasing in stock at the fastest rate.
Tier 3 correlation-based method vs. Tier 1 and 2 methods
Our estimates of 1.3 Pg C stored across the 33.9 million hectares is larger than most previous Tier 1 estimates [46–48], although below the most recently produced estimate  (Table 1). Underestimation of the amount of carbon stored in the EAM region in global analyses can be a result of their poor resolution and/or application of data from other regions which may differ systematically compared to East African forests, woodlands and savannas . When separated by land cover category, our locally derived carbon estimates are comparable to those presented in other local [49–52] and global studies, the latter often containing little or no data from East Africa [3, 4, 46, 47, 53]. This suggests differences between our estimates and other studies have arisen because many previous studies mapped carbon storage at lower resolution [3, 4, 46, 47, 53]. When considering homogenous landscapes, scale effects are unlikely to cause a dramatic difference in carbon estimates. However, in highly fragmented and heterogeneous landscapes, such as East Africa, the effects of scale are likely to be substantial. Forest fragments, typically of high carbon storage, may be omitted at lower resolutions, being ‘replaced’ by more dominant, but low carbon, land cover categories (e.g. open woodland), resulting in underestimation of carbon storage.
It must be noted that, the landscape-scale confidence intervals surrounding our Tier 3 estimates are considerably wider than those around previous estimates [3, 4, 42, 47, 53]. This result is consistent with Hill et al (2013), who also showed increasing methodological sophistication does not necessarily result in reduced uncertainty, as is often assumed . Confidence intervals derived from look-up table values may show a systematic bias. The ranges provided are an artefact of the study area, the number of land cover categories and the resolution, as when summed across a large number of pixels, pixel error is mostly negated as underestimates in one part of the landscape are counterbalanced by overestimates in other parts. The 95% CI developed from correlation equations are effectively based on numerous continuous variables, containing the uncertainty relating to anthropogenic, climatic and edaphic variables, thus have many thousands of possible combinations, severely limiting the ability of the ‘law of averages’ to act. Hence, the 95% CI presented in this investigation may better reflect that of the actual landscape, containing more variables that make-up the complex landscape heterogeneity (i.e. improved representativeness), although this is only true for those pixels estimated using the correlation equations (86% of the EAM but only 52% of the study area). Therefore, the look-up table 95% CI presented in Willcock et al (2012), and used in this study, may underestimate uncertainty . Future studies should expand the existing plot network (Figure 1), enabling the correlation equations (and improved 95% CI) to be applied to the entire study area. This process has already begun under a new WWF-REDD+ project (which focusses on better sampling the data-deficient land cover categories identified in this study ) and the National Forest Monitoring and Assessment (NAFORMA) project [56, 57].
Links between carbon stock and influential variables
The results presented here indicate that ALC storage in tree-dominated ecosystems is correlated with anthropogenic, climatic and edaphic variables. However, in all our models there is a large amount of unexplained variation (R-squared values for our correlation models vary between 0.18 and 0.41). This is likely to be due to three main reasons (Additional file 1: SI6). Firstly, although we used the highest resolution datasets that are freely available, several of the associated variables are of relatively poor resolution across the EAM (including; wind, light and soil nutrient variables [Additional file 1: Table S6]). This is particularly important here as low resolution GIS data is unlikely to correlate well with the response variables from our plot network as many plots (with high variance ) may fall within a single cell . Thus, our study may be biased against retaining low resolution explanatory variables in our models. Secondly, contemporary forest characteristics are the result of growth, recruitment and mortality over many years. It is difficult to obtain data on historical variables and yet these could have had a significant impact on present day carbon storage and other forest characteristics . Thirdly, present day information is also lacking, for example datasets describing physical soil properties in the study area are unavailable. Thus, future work is needed to develop additional high resolution GIS data, particularly for historic time periods.
Of the variance explained in our forward and backward models, direct anthropogenic factors are the most influential explanatory variables (as noted by the largest coefficient of explanatory variables on the response variable, in contrast to those [e.g. temperature] with smaller p-values but also smaller influence [Table 3]) and so are the focus of our remaining discussion (see Additional file 1: SI5 for discussion of climatic and edaphic variables).
Within our study area, people are clustered around high carbon areas (Figure 4). We suggest this could be due to these areas having favourable climatic conditions with more moisture for plant (and thus crop) growth. Further, the incidence of malaria is lower at high elevations , making these locations more habitable for human populations. Thus there is a peak in population density near the base of high-carbon montane forests . Our interpretation that it is the landscape suitability driving human population density is consistent with the observation that when individual localities are followed over time, degradation at the local level caused by the population is evident [62, 63]. This emphasises that our results are not proof of causation and that the drivers may be a correlate of the explanatory variables retained in our models (Additional file 1: SI6). Our results also show a decrease in carbon storage in previously logged areas and in areas nearer the commercial capital, Dar es Salaam. This confirms previous reports that areas near the capital have lower biomass due to the local demand of low grade timber by the city, as well as international demand for high grade timber via the city’s port ; emphasising the connections between the rural and urban landscape, and how the sphere of urban influence drives change in rural ecosystems. Future investigations should use simulation modelling and direct experimentation to identify if the influential variables highlighted here can be confirmed as drivers of carbon storage and sequestration, providing a deeper understanding of the process-based relationships.
The decrease in carbon storage as a result of logging (51-77% of the ALC is retained) is of similar magnitude to other reported estimates . However, the historical logging data we utilised was based on expert opinion (Additional file 1: Table S6) so, given its importance, further work developing and evaluating historical variables is needed (Additional file 1: Table S7). We observe a comparable decrease due to differing governance. Land under national control holds between 40% and 65% of the ALC stored in areas under decentralised governance. This perhaps indicates that decentralisation of management (e.g. participatory and community led forestry) is successful in our study area [37, 65]. However, it is not possible to prove causation within the framework of this study. Many locally managed forests are located in the south-east of our study area within an area of naturally high carbon storage, whereas land under national control covers much larger areas, including the dry, carbon-poor east. Hence, our finding that carbon storage is higher in areas under decentralised control may be an artefact of the differing areas where this type of land management occurs. Further studies monitoring change in carbon storage over time under the two different governance regimes would enable the effect of land management to be determined.
The overall effects on carbon storage are a result of many changes in forest characteristics. Both WSG and the proportion of larger stems decrease with increasing anthropogenic disturbance, however, stem density (≥ = 10 cm DBH) increases. Anthropogenic disturbance, for example logging, is often a commercial activity and results in the preferential removal of the largest, most valuable stems . The more open canopy, following stem removal, would result in increased recruitment from young forest trees , leading to the high numbers of small stems observed. However, the opposite would be expected in woodlands and savannas, with more open canopies resulting in more grass, high fire intensity and so less recruitment [67, 68]. Our results highlight how influential the negative effect of people on tropical forest carbon storage can be. This assertion is supported by data from across the tropics [69–71]. The significant impact of anthropogenic activities implies that REDD+ could, at the local scale, have significant positive impacts on carbon storage. However, careful policy designs to limit leakage of deforestation and encourage the involvement of the local population are needed to ensure REDD+ schemes achieve their carbon storage and sequestration aims .
Like carbon storage and its components, carbon sequestration is also correlated with anthropogenic, climatic and edaphic variables. We estimate that some localities (for example the Udzungwa Mountains National Park; Figure 4) provide a carbon sink of comparable per-area magnitude to modelled estimates in East Africa  and to that observed over recent decades in structurally intact African forest . However, many areas of forest and woodland within the study area experience a high level of degradation and disturbance, and so are net sources. Here, we have shown that anthropogenic disturbance is a key determinant of the trend in carbon storage over time in eastern Tanzania. Important locations of high carbon losses are the Pare and Usambara mountains (Table 5), which historically have seen the highest rates of degradation and disturbance . The national population of Tanzania is increasing  and this may increase the pressure on tree-dominated ecosystems which could result in the study area becoming a significant source of carbon in the future. Furthermore, the effect of increase in anthropogenic pressures could be compounded by potential decrease in carbon storage as a result of increasing temperatures [76, 77] and changes in soil nutrients (see Additional file 1: SI5). However, these future effects could be complicated by increasing levels of atmospheric CO2, varying effectiveness of legally protected areas and shifting consumption patterns.
Our results show that the amount of carbon stored in forests across 33.9 million ha of the Eastern Arc Mountains of Tanzania is considerable: 1.32 (0.89 to 3.16) Pg. Our estimate is significantly higher than most previous estimates. However, our more sophisticated method also has higher uncertainty, implying that other methods may substantially underestimate the uncertainty involved. Within the tree-dominated land cover categories, historical logging is the most influential direct anthropogenic factor, while the mean number of dry months is the most influential environmental factor, with an order of magnitude less impact on carbon storage. We show that WSG, size-frequency distribution variables and height variables are all important in determining carbon storage. Our estimates indicate that, between 2004 and 2008, tree-dominated communities across the study areas showed no significant change, however some areas were identified as large sinks (0.8% of the study area) and others large sources (1.4% of the study area), showing the importance of taking a landscape scale approach. The carbon maps produced and statistical relationships documented can assist policy-makers in designing policies to maintain and enhance carbon storage for climate mitigation and other ecosystem services.
We collated data from 2,462 tree inventory plots within our study area (see Additional file 1: SI3), then applied a quality control and standardisation protocol. This consists of two main steps: (1) Metadata quality control; and (2) Measurement bias detection.
Firstly, all plots lacking a recorded spatial location and a fixed area were discarded (770 plots). Plots where one or more diameter at breast height (DBH) data were known to be missing were also excluded (7 plots). Furthermore, plots smaller than 0.025 ha (16 plots) were deemed to produce unreliable carbon estimates so also removed from the dataset.
Secondly, to assess possible measurement bias, i.e. not measuring over buttresses and so overestimating biomass , the remaining plots were grouped by the lead field researcher. Size-frequency distributions, using 10 cm size classes, were created for each of these groups. Forest size-frequency distributions are suggested to conform to the -2 power law based on metabolic scaling . Although it has been argued that this rule is not globally applicable , many studies accept this as a theoretical maximum value for the abundance of large stems . Thus, researchers with many plots above this maximum value likely measured stems around buttresses and so were removed (1 researcher, 100 Plots).
The quality control and standardisation procedure resulted in a dataset of 1,611 tree inventory plots (median 0.1 ha, mean 0.1 ha, mode 0.1 ha [43 plots with multiple censuses; median 0.1 ha, mean 0.5 ha, mode 1.0 ha]; Figure 1; see Additional file 1: SI3 for a further information) from which we calculated plot-level stand structure indices and aboveground carbon storage per unit area (see Additional file 1: SI2 for full details). We obtained the exponent and intercept of the population size-frequency distribution using the power law fit for each plot using the log-log transformation method. Whereby, for each plot, we created 10 cm bin size-frequency distributions based on DBH, and a linear model of the logarithm of frequency against the logarithm of the size class was fitted. Whilst not as accurate as the maximum likelihood estimation method, our simpler method is more stable for many of our plots, providing both the intercept and slope indicators of population structure .
We obtained WSG data via the phylogenetic information provided by our tree inventory plots. We used a global wood density database to extract species average WSG . This procedure provided over 32,000 trees with WSG data. When this was not possible we adopted a hierarchical approach, first applying the appropriate genus average if available (~14,000 trees) before considering family average (~9,500 trees), plot average (~4,500 trees) and dataset average (~80 trees) in turn . Including WSG as an additional parameter in allometric equations reduces the biomass estimation error [49, 85, 86].
In addition, we estimated plot biomass using moist forest tree allometry  based on measurements of DBH from our tree inventory plots, WSG (as described above) and height data (derived from our dataset using the best fit DBH-height equation form [Equation 5.1; see Additional file 1: SI4], if not measured in the tree inventory plots). Finally, carbon was assumed to be 50% of biomass .
For a smaller number of plots, multiple measurements were available over time (n = 43; mean plot size = 0.5 ha; mean measurement period = 3.9 years). We calculated changes in carbon storage rates by dividing the difference in carbon storage estimates between censuses by the number of years separating them.
For our 1,611 geo-referenced tree inventory plots, we obtained further information on variables falling into five broad categories; anthropogenic, climatic, geographic, edaphic, and pyrologic (median resolution 1.0 ha, mean resolution 22.0 ha, mode resolution 1.0 ha; Additional file 1: Table S6). Anthropogenic data, further divided into six subcategories, were obtained: (1) population pressure variables (n = 14 related variables) were obtained from Platts (2012)  (see Additional file 1: SI7); (2) Dar es Salaam related variables (n = 3; e.g. distance to Dar es Salaam), (3) market town related variables (n = 3; e.g. distance to market towns), and (4) infrastructure related variables (n = 2; e.g. distance to roads) were derived from available topographic maps; (5) historical logging (n = 1) from Swetnam et al (2011) ; and (6) governance (n = 1) from the World Database on Protected Areas . Climate data were divided into three subcategories (precipitation [n = 2; maximum mean cumulative water deficit and mean number of dry months annually], temperature [n = 4; mean annual temperature, mean annual minimum monthly temperature, mean annual monthly maximum temperature, and mean annual monthly temperature range] and wind speed [n = 1]) and were derived from the Tropical Rainfall Measuring Mission [90, 91], WorldClim [92, 93], and United States National Aeronautics and Space Administration Surface meteorology and Solar Energy  datasets. Similarly, geographic data have two variables (aspect [n = 1] and incoming solar radiation [n = 1]) derived from Shuttle Radar Topography Mission  and National Renewable Energy Laboratory [95, 96] datasets respectively. Lastly, we extracted edaphic data (n = 6) from the International Soil Reference and Information Centre database [97, 98] and fire-related variables (n = 5) derived from MODIS images .
We then correlated these variables with carbon storage, and following this, its components: WSG, the intercept of the power law relationship, and the gradient of the power law relationship, in each case using general linear models (see Additional file 1: SI2-5). No transformations were required to ensure a normal distribution when correlating either WSG, the intercept of the power law relationship or the gradient of the power law relationship with the individual variables. However, carbon storage estimates required a square root transformation to ensure a normal distribution within the general linear models (normality was confirmed using the Shapiro-Wilk test; p-value > 0.05). In all models, plots were weighted by the square root of their area as confidence in biomass estimation increases with the area surveyed [100, 101]. Landscape scale spatial autocorrelation was accounted for by including spatial terms (latitude, longitude and the interactions between them) in the model (Additional file 1: Table S6) . The numerous possible interactions were excluded from the models, as these were found to add very little explanatory power to the models, only increasing R-squared values by ~0.001 with the addition of each interaction term. All analyses were performed using R 2.12.1  and mapped in ArcGIS v9.3.1 .
When assessing carbon sequestration (n = 43) fewer degrees of freedom were available, therefore explanatory variables need to be grouped. Therefore, we conducted a principle components (PC) analysis, obtaining five PC which explained >90% of the cumulative variance of the individual influential variables (Additional file 1: Table S4). Then, covariation of PC with carbon sequestration was assessed instead of the individual influential variables. Carbon sequestration estimates required a cube-root transformation to ensure a normal distribution within the general linear models (confirmed using the Shapiro-Wilk test; p-value > 0.05). This enabled the effect of multiple variables to be examined even with this limited dataset. PC analysis of the variables was performed on the scaled data using the prcomp package  within R 2.12.1 . All other aspects of the model (weighting and spatial autocorrelation) were performed identically to the models for carbon storage and its components.
The most appropriate model was chosen using forward and backward stepwise selection. Forward models are more useful for inferring causal relationships  and so were preferentially used to infer the influential variables of carbon storage and sequestration. However, averaging forward–backwards and backward–forwards predictions outperforms conventional selection procedures  and so both methods were used when estimating the spatial distributions within the study area. Akaike information criterion (AIC) was used to reduce/expand the models, with variable selection occurring when the variable reduced the mean squared error (MSE) under ten-fold cross validation . Unlike model selection using R-squared, which neglects the principles of parsimony, AIC considers both model fit and complexity, resulting in better predictions and allowing inferences to be made from multiple models . Model selection continued until the addition/removal of further variables able to reduce cross validation MSE no longer increased AIC, thereby producing the best-fit model with the lowest prediction error .
Within each category (anthropogenic, climatic, geographic, edaphic, and pyrologic), some variables were highly correlated (Additional file 1: Table S7) and this may confound the stepwise procedure as each variable does not carry enough distinct information . For example, all temperature related variables (Additional file 1: Table S7) were correlated (R-squared > 0.6). However, it is unclear which correlated best with the variables of interest, e.g. carbon storage and sequestration. Many studies include mean annual temperature in biomass models [77, 109], but theory suggests that it may be the temperature range driving this relationship as photosynthesis correlates with maximum temperatures, but respiration with minimum temperatures [76, 110, 111]. We found that, if we removed correlated variables prior to model selection, the final models were artefacts of the variables we had selected. For example, if we included mean annual temperature in the model, but not temperature range, then the significant correlations between mean annual temperature and ALC storage were found. However, these correlations were insignificant if temperature range was added to the model, with the newly added variable showing a significant effect instead. In short, the resultant models were automatically biased towards a priori expectations. To avoid this bias, we devised a procedure by which the influential variables included in model selection were selected by their ability to explain variation within the data of interest (e.g. carbon storage). All variables (describe above) were included in model selection. Once this had run to completion the model was assessed. The subcategory with the most correlated variables retained within the model was selected and all but the most influential, significant variable were removed. For example, if all four temperature-related variables were included in the initial model and this was the largest group of variables then this group would be selected. Then, if mean annual temperature was the most influential and significant temperature-related variable, all other temperature-related variables would be excluded in the next round of model selection. Thus, stepwise model selection was then repeated for all remaining variables. This process was repeated until no highly correlated variables remained within the model produced.
Since only landscape-scale variation was accounted for by the spatial terms already included in the model (latitude, longitude and the interactions between them; Table 1; Additional file 1: Table S6), it was necessary to investigate the effect of local-scale (<10 km2) spatial autocorrelation . To do this, the separate forward and backward models, containing no highly correlated variables (produced above), were mapped. Then, the sum of the model estimates within the maps were extracted at 1, 3, 5, 7 and 10 km2 resolutions, and included as additional variables (representing local spatial autocorrelation terms) into the stepwise model selection process, which was re-run a final time . However, in all cases, local spatial autocorrelation terms were rejected as they did not reduce cross validated MSE.
Since it was not necessary to include local spatial autocorrelation terms in the models, the preliminary maps produced above could be regarded as final spatial representations of the ten best fit models, two (forward and backward) for each of the five variables of interest (carbon storage, carbon sequestration, WSG, the intercept of the power law relationship and the gradient of the power law relationship). Each pair of maps (forward and backward) were then combined into a single, final weighted mean estimate. The ratio of the relevant cross validated MSE of the forward and backward models was used to create the weighted mean, with the model showing lowest error receiving the highest weighting . Thus, we ultimately produced five maps (from ten best fit models); one each for carbon storage, carbon sequestration, WSG, the intercept of the power law relationship, and the gradient of the power law relationship. As our carbon storage estimates were derived from data representing trees with a DBH greater than or equal to 10 cm, regionally estimates of ratios from Willcock et al (2012) were used to estimate the unmeasured component of ALC storage , this was summed with our modelled carbon storage estimate, providing an estimate of total ALC storage.
Although the five maps produced covered the entire study area, we were concerned that extrapolating predictions beyond the range of observed predictor variables from our dataset could result in large, unquantifiable errors. Thus, we limited the models to localities where all the associate variables were within the range of that shown in our dataset, thus only interpolating within our correlation models for tree-dominated land cover categories. For any pixels outside the data range, look-up table methods were used in preference to the correlation model estimates. Thus, for every land cover in our study area containing trees (open woodland; closed woodland; forest mosaic; lowland forest; sub-montane forest; montane forest; and upper montane forest ) that fell within the limits of our dataset, the estimate of carbon storage derived from the correlation equations was used. For all other land cover categories, and for those localities for which predictor variables fell outside the ranges of values used in model construction, land cover based look-up table values from Willcock et al (2012) were used to estimate ALC storage . In total, look-up table values were applied to 52% of the landscape, although this was predominantly to low carbon land cover categories, with 86% of the EAM (which hold the majority of the regions tropical forest ) estimated using the correlation approach described above. Estimates of WSG and population structure were only made for wooded land cover categories, with estimates for areas within our dataset range being derived from the relevant correlation equations and estimates for other areas coming from land cover based look-up table values derived from the median value of our WSG and population structure data (weighted by the square root of plot size and derived via sampling with replacement 10,000 times) for each land cover category (Additional file 1: Table S8). For carbon sequestration, again, estimates were only made for wooded land cover categories for those areas inside the range of our dataset estimates derived from the correlation equations were used. However, unlike carbon storage, WSG and population structure, for areas outside the range of our dataset, a land cover based look-up table was not used as several land cover categories were poorly represented due to the small sample size available (n = 43). Instead, for pixels outside the range of the correlation-derived carbon sequestration model (16% of pixels with wooded land cover), the median value of data from our recensused plots (again weighted by the square root of plot size and derived via sampling with replacement 10,000 times) was utilised.
For every 1 ha pixel of each map derived from correlation equations, we produced 95% confidence intervals (CI). If the pixel estimate was derived from the general linear models, then the pixel 95% CI was calculated by adding and subtracting the square root of the cross validation MSE. For look-up table pixels the look up table 95% CI were used. The pixel 95% CI describes, for every pixel, the range we would expect each of our estimates to lie within. However, as we are also interested in estimating carbon storage and sequestration on a landscape scale, indications of uncertainty are also required at landscape-scale. Simply summing the pixel 95% CI to derive 95% CI of the overall landscape-scale estimates would incorrectly treat random error as a region-wide systematic bias. Thus, to derive 95% CI for landscape-scale estimates, we randomly allocated each pixel an estimate within the range dictated by its 95% pixel CI, and summed these values across the entire landscape. This process was performed 10,000 times and the median value and 95% CI (the 250th and 9,750th ranked values, which may not be equally distributed around the median) for aboveground carbon storage and sequestration in the study area were obtained.
For the final model of carbon storage estimates, we investigated how the components of carbon storage (population structure, WSG and tree height) interacted to ultimately produce the ecosystem service of carbon storage. We obtained estimates of maximum canopy height from the best fit DBH-height equation [Equation 5.1; see Additional file 1: SI4 and Additional file 7: Figure S6], and combined this spatially with our correlation model derived estimates of WSG, the intercept of the power law relationship and gradient of the power law relationship. We then correlated these against our estimates of carbon storage, allowing all possible interactions, and selected the best-fit model (via AIC) using both forwards and backwards stepwise regression.
Ethical approval for the above study was obtained from the Faculty of Environment Research Ethics Committee, in accordance with the University of Leeds research ethics policy.
Akaike information criteria
Aboveground live carbon
Coarse woody debris
Diameter at breast height
Eastern Arc Mountains
Effective cation exchange capacity
Geographic information systems
History Database of the Global Environment
Intergovernmental Panel on Climate Change
International Union for Conservation of Nature
York Institute for Tropical Ecosystems
Mean annual temperature
Mean squared error
Reducing Emissions from Deforestation and Forest degradation
Centre for Sustainability and Global Environment
Wood specific gravity.
Malhi Y, Grace J: Tropical forests and atmospheric carbon dioxide. Trends Ecol Evol 2000, 15: 332–337. 10.1016/S0169-5347(00)01906-6
IPCC: Land use, land-use change, and forestry. Cambridge, UK: Cambridge University Press; 2000.
Baccini A, Goetz SJ, Walker WS, Laporte NT, Sun M, Sulla-Menashe D, Hackler J, Beck PSA, Dubayah R, Friedl MA, Samanta S, Houghton RA: Estimated carbon dioxide emissions from tropical deforestation improved by carbon-density maps. Nat Clim Change 2012, 2: 182–185. 10.1038/nclimate1354
Saatchi SS, Harris NL, Brown S, Lefsky M, Mitchard ETA, Salas W, Zutta BR, Buermann W, Lewis SL, Hagen S, Petrova S, White L, Silman M, Morel A: Benchmark map of forest carbon stocks in tropical regions across three continents. Proc Natl Acad Sci 2011, 108: 9899–9904. 10.1073/pnas.1019576108
Myers N, Mittermeier RA, Mittermeier CG, da Fonseca GAB, Kent J: Biodiversity hotspots for conservation priorities. Nature 2000, 403: 853–858. 10.1038/35002501
Timko JA, Waeber PO, Kozak RA: The socio-economic contribution of non-timber forest products to rural livelihoods in Sub-Saharan Africa: knowledge gaps and new directions. Int For Rev 2010, 12: 284–294.
Lewis SL, Lopez-Gonzalez G, Sonke B, Affum-Baffoe K, Baker TR, Ojo LO, Phillips OL, Reitsma JM, White L, Comiskey JA, Djuikouo M-N, Ewango CEN, Feldpausch TR, Hamilton AC, Gloor M, Hart T, Hladik A, Lloyd J, Lovett JC, Makana J-R, Peacock J, Peh KS-H, Sheil D, Sunderland T, Swaine M, Taplin J, Taylor D, Thomas C, Votere R, Wöll H: Increasing carbon storage in intact African tropical forests. Nature 2009, 457: 1003–1006. 10.1038/nature07771
Phillips OL, Malhi Y, Higuchi N, Laurance WF, Nunez PV, Vasquez RM, Laurance SG, Ferreira LV, Stern M, Brown S, Grace J: Changes in the carbon balance of tropical forests: evidence from long-term plots. Science 1998, 282: 439–442.
Achard F, Eva HD, Stibig H-J, Mayaux P, Gallego J, Richards T, Malingreau J-P: Determination of deforestation rates of the world's humid tropical forests. Science 2002, 297: 999–1002. 10.1126/science.1070656
Hansen MC, Stehman SV, Potapov PV: Quantification of global gross forest cover loss. Proc Natl Acad Sci 2010, 107: 8650–8655. 10.1073/pnas.0912668107
Putz FE, Zuidema PA, Pinard MA, Boot RGA, Sayer JA, Sheil D, Sist P, Elias M, Vanclay JK: Improved Tropical Forest Management for Carbon Retention. PLoS Biol 2008, 6: e166. 10.1371/journal.pbio.0060166
IPCC: Climate Change 2007: The Physical Science Basis. Available at [Accessed 05/01/12]. Agenda 2007, 6 http://www.ipcc.ch/pdf/assessment-report/ar4/wg1/ar4-wg1-spm.pdf Available at [Accessed 05/01/12]. Agenda 2007, 6
van der Werf GR, Morton DC, DeFries RS, Olivier JGJ, Kasibhatla PS, Jackson RB, Collatz GJ, Randerson JT: CO2 emissions from forest loss. Nat Geosci 2009, 2: 737–738. 10.1038/ngeo671
Dixon RK, Borwn S, Houghton RA, Solomon AM, Trexler MC, Wisniewski J: Carbon pools and fluxes of global forest ecosystems. Science 1994, 263: 185–190. 10.1126/science.263.5144.185
DeFries RS, Houghton RA, Hansen MC, Field CB, Skole D, Townshend J: Carbon emissions from tropical deforestation and regrowth based on satellite observations for the 1980s and 1990s. Proc Natl Acad Sci U S A 2002, 99: 14256–14261. 10.1073/pnas.182560099
Houghton RA: TRENDS: A Compendium of Data on Global Change. Available at: . An update of estimated carbon emissions from land use change, typically updated every year. Oak Ridge National Laboratory, US Department of Energy; 2008. http://cdiac.ornl.gov/trends/landuse/houghton/houghton.html
Pan Y, Birdsey RA, Fang J, Houghton R, Kauppi PE, Kurz WA, Phillips OL, Shvidenko A, Lewis SL, Canadell JG, Ciais P, Jackson RB, Pacala SW, McGuire AD, Piao S, Rautiainen A, Sitch S, Hayes D: A large and persistent carbon sink in the world’s forests. Science 2011, 333: 988–993. 10.1126/science.1201609
Hansen MC, Stehman SV, Potapov PV, Loveland TR, Townshend JRG, DeFries RS, Pittman KW, Arunarwati B, Stolle F, Steininger MK, Carroll M, DiMiceli C: Humid tropical forest clearing from 2000 to 2005 quantified by using multitemporal and multiresolution remotely sensed data. Proc Natl Acad Sci 2008, 105: 9439–9444. 10.1073/pnas.0804042105
Chave J, Condit R, Muller-Landau HC, Thomas SC, Ashton PS, Bunyavejchewin S, Co LL, Dattaraja HS, Davies SJ, Esufali S: Assessing evidence for a pervasive alteration in tropical tree communities. PLoS Biol 2008, 6: e45. 10.1371/journal.pbio.0060045
Field CB, Behrenfeld MJ, Randerson JT, Falkowski P: Primary production of the biosphere: integrating terrestrial and oceanic components. Science 1998, 281: 237–240.
REDD overview: REDD overview. [http://www.un-redd.org/] 
Strassburg B, Turner RK, Fisher B, Schaeffer R, Lovett A: Reducing emissions from deforestation—The “combined incentives” mechanism and empirical simulations. Glob Environ Chang 2009, 19: 265–278. 10.1016/j.gloenvcha.2008.11.004
Kindermann G, Obersteiner M, Sohngen B, Sathaye J, Andrasko K, Rametsteiner E, Schlamadinger B, Wunder S, Beach R: Global cost estimates of reducing carbon emissions through avoided deforestation. Proc Natl Acad Sci U S A 2008, 105: 10302–10307. 10.1073/pnas.0710616105
Cattaneo A, Lubowski R, Busch J, Creed A, Strassburg B, Boltz F, Ashton R: On international equity in reducing emissions from deforestation. Environ Sci Pol 2010, 13: 742–753. 10.1016/j.envsci.2010.08.009
Farris M: The sound of falling trees: integrating environmental justice principles into the climate change framework for Reducing Emissions from Deforestation and Degradation (REDD). Fordham Environ Law Rev 2010, 20: 515.
Birdsey R, Pan Y, Houghton R: Sustainable landscapes in a world of change: tropical forests, land use and implementation of REDD+: Part II. Carbon Manag 2013, 4: 567–569. 10.4155/cmt.13.67
IPCC: Guidelines for National Greenhouse Gas Inventories. In National Greenhouse Gas Inventories Programme. Edited by: Eggleston HS, Buendia L, Miwa K, Ngara T, Tanabe K. Kanagawa, Japan: Institute for Global Environmental Strategies; 2006.
IPCC: Good practice guidance for land use, land-use change and forestry. In IPCC National Greenhouse Gas Inventories Programme. Edited by: Penman J, Gytarsky M, Hiraishi T, Krug T, Kruger D, Pipatti R, Buendia L, Miwa K, Ngara T, Tanabe K. Kanagawa, Japan: Institute for Global Environmental Strategies; 2003.
GOFC-GOLD A: A Sourcebook of Methods and Procedures for Monitoring and Reporting Anthropogenic Greenhouse Gas Emissions and Removals Caused by Deforestation, Gains and Losses of Carbon Stocks in Forests Remaining Forests, and Forestation. AB, Canada: GOFC-GOLD GOFC-GOLD Project Office, Natural Resources Canada Calgary; 2010.
Asner GP, Powell GVN, Mascaro J, Knapp DE, Clark JK, Jacobson J, Kennedy-Bowdoin T, Balaji A, Paez-Acosta G, Victoria E, Secada L, Valqui M, Hughes RF: High-resolution forest carbon stocks and emissions in the Amazon. Proc Natl Acad Sci 2010, 107: 16738–16742. 10.1073/pnas.1004875107
International W: Carbon Storage in the Los Amigos Conservation Concession, Madre de Dios, Perú. Washington DC, USA: Winrock International; 2006.
Asner GP, Mascaro J, Clark JK, Powell G: Reply to Stoke et al.: Regarding high-resolution carbon stocks and emissions in the Amazon. Proc Natl Acad Sci 2011, 108: E13-E14. 10.1073/pnas.1017675108
Cláudia Dias A, Louro M, Arroja L, Capela I: Comparison of methods for estimating carbon in harvested wood products. Biomass Bioenergy 2009, 33: 213–222. 10.1016/j.biombioe.2008.07.004
Pedroni L, Dutschke M, Streck C, Porrúa ME: Creating incentives for avoiding further deforestation: the nested approach. Clim Pol 2009, 9: 207–220. 10.3763/cpol.2008.0522
Hardcastle P, Baird D, Harden V, Abbot PG, O’Hara P, Palmer JR, Roby A, Haüsler T, Ambia V, Branthomme A, Wilkie M, Arends E, González C: Capability and Cost Assessment of the Major Forest Nations to Measure and Monitor their Forest Carbon. Edinburg, Scotland: LTS International Ltd; 2008.
Romijn E, Herold M, Kooistra L, Murdiyarso D, Verchot L: Assessing capacities of non-Annex I countries for national forest monitoring in the context of REDD+. Environ Sci Pol 2012, 19–20: 33–48.
Burgess ND, Bahane B, Clairs T, Danielsen F, Dalsgaard S, Funder M, Hagelberg N, Harrison P, Haule C, Kabalimu K, Kilahama F, Kilawe E, Lewis SL, Lovett LC, Lyatuu G, Marshall AR, Meshack C, Miles L, Milledge SAH, Munishi PKT, Nashanda E, Shirima D, Swetnam R, Willcock S, Williams A, Zahabu E: Getting ready for REDD+ in Tanzania: a case study of progress and challenges. Oryx 2010, 44: 339–351. 10.1017/S0030605310000554
Maniatis D, Malhi Y, Andr LS, Mollicone D, Barbier N, Saatchi S, Henry M, Tellier L, Schwartzenberg M, White L: Evaluating the Potential of Commercial Forest Inventory Data to Report on Forest Carbon Stock and Forest Carbon Stock Changes for REDD+ under the UNFCCC. Int J Forest Res 2011, 2011: 1–13.
Science for a changing world http://www.usgs.gov/default.asp
Platts PJ, Burgess ND, Gereau RE, Lovett JC, Marshall AR, Mcclean CJ, Pellikka PKE, Swetnam RD, Marchant R: Delimiting tropical mountain ecoregions for conservation. Environ Conserv 2011, 38: 312–324. 10.1017/S0376892911000191
Swetnam RD, Fisher B, Mbilinyi BP, Munishi PKT, Willcock S, Ricketts T, Mwakalila S, Balmford A, Burgess ND, Marshall AR, Lewis SL: Mapping socio-economic scenarios of land cover change: a GIS method to enable ecosystem service modelling. J Environ Manag 2011, 92: 563–574. 10.1016/j.jenvman.2010.09.007
Willcock S, Phillips OL, Platts PJ, Balmford A, Burgess ND, Lovett JC, Ahrends A, Bayliss J, Doggart N, Doody K, Fanning E, Green JMH, Hall J, Howell KL, Marchant R, Marchant R, Marshall A, Mbilinyi B, Munishi PKT, Owen N, Swetnam RD, Topp-Jorgensen EJ, Lewis SL: Towards regional, error-bounded landscape carbon storage estimates for data-deficient areas of the world. PLoS ONE 2012, 7: e44795. 10.1371/journal.pone.0044795
Platts PJ, McClean CJ, Lovett JC, Marchant R: Predicting tree distributions in an East African biodiversity hotspot: model selection, data bias and envelope uncertainty. Ecol Model 2008, 218: 121–134. 10.1016/j.ecolmodel.2008.06.028
Stuart EA: Matching methods for causal inference: A review and a look forward. Stat Sci Rev J Inst Math Stat 2010, 25: 1.
Wiegand RE: Performance of using multiple stepwise algorithms for variable selection. Stat Med 2010, 29: 1647–1659.
Hurtt GC, Frolking S, Fearon MG, Moore B, Shevliakova E, Malyshev S, Pacala SW, Houghton RA: The underpinnings of land-use history: three centuries of global gridded land-use transitions, wood-harvest activity, and resulting secondary lands. Glob Chang Biol 2006, 12: 1208–1229. 10.1111/j.1365-2486.2006.01150.x
Baccini A, Laporte N, Goetz SJ, Sun M, Dong H: A first map of tropical Africa's above-ground biomass derived from satellite imagery. Environ Res Lett 2008, 3: 045011. 10.1088/1748-9326/3/4/045011
Saatchi SS, Houghton RA, Dos Santos AlvalÁ RC, Soares JV, Yu Y: Distribution of aboveground live biomass in the Amazon basin. Glob Chang Biol 2007, 13: 816–837. 10.1111/j.1365-2486.2007.01323.x
Marshall AR, Willcock S, Lovett JC, Balmford A, Burgess ND, Latham JE, Munishi PKT, Platts PJ, Salter R, Shirima DD, Lewis SL: Measuring and modelling above-ground carbon storage and tree allometry along an elevation gradient. Biol Conserv 2012, 154: 20–33.
Munishi PKT, Shear TH: Carbon storage in afromontane rain forests of the Eastern Arc Mountains of Tanzania: their net contribution to atmospheric carbon. J Trop For Sci 2004, 16: 78–93.
Shirima DD, Munishi PKT, Lewis SL, Burgess ND, Marshall AR, Balmford A, Swetnam RD, Zahabu EM: Carbon storage, structure and composition of miombo woodlands in Tanzania’s Eastern Arc Mountains. Afr J Ecol 2011, 49: 332–342. 10.1111/j.1365-2028.2011.01269.x
Pfeifer M, Platts PJ, Burgess ND, Swetnam RD, Willcock S, Lewis SL, Marchant R: Land use change and carbon fluxes in East Africa quantified using earth observation data and field measurements. Environ Conserv 2013, 40: 1–12. FirstView FirstView 10.1017/S0376892912000276
Ruesch A, Gibbs HK: New IPCC Tier1 Global Biomass Carbon Map For the Year 2000. Oak Ridge, Tennessee: Oak Ridge National Laboratory; 2008.http://cdiac.ornl.gov/ Available online from the Carbon Dioxide Information Analysis Center [Accessed 15/01/12]
Hill TC, Williams M, Bloom AA, Mitchard ETA, Ryan CM: Are inventory based and remotely sensed above-ground biomass estimates consistent? PLoS ONE 2013, 8: e74170. 10.1371/journal.pone.0074170
Burgess ND, Mwakalila S, Munishi P, Pfeifer M, Willcock S, Shirima D, Hamidu S, Bulenga GB, Rubens J, Machano H, Marchant R: REDD herrings or REDD menace: response to Beymer-Farris and Bassett. Glob Environ Chang 2013, 23: 1349–1354. 10.1016/j.gloenvcha.2013.05.013
Kessy JF, Anderson K, Dalsgaar S: NAFORMA Field Manual: Socioeconomic survey. Dar es Salaam, Tanzania: Forestry and Beekeeping Division, Ministry of Natural Resources and Tourism; 2010:96.
Vesa L, Malimbwi RE, Tomppo E, Zahabu E, Maliondo S, Chamuya N, Nssoko E, Otieno J, Miceli G, Kaaya AK, Dalsgaar S: NAFORMA Field Manual: Biophysical survey. Dar es Salaam, Tanzania: Forestry and Beekeeping Division, Ministry of Natural Resources and Tourism; 2010:96.
Chave J, Condit R, Lao S, Caspersen JP, Foster RB, Hubbell SP: Spatial and temporal variation of biomass in a tropical forest: results from a large census plot in Panama. J Ecol 2003, 91: 240–252. 10.1046/j.1365-2745.2003.00757.x
Baccini A, Friedl M, Woodcock C, Zhu Z: Scaling field data to calibrate and validate moderate spatial resolution remote sensing models. Photogramm Eng Remote Sens 2007, 73: 945. 10.14358/PERS.73.8.945
Ramankutty N, Gibbs HK, Achard F, Defries R, Foley JA, Houghton RA: Challenges to estimating carbon emissions from tropical deforestation. Glob Chang Biol 2007, 13: 51–66. 10.1111/j.1365-2486.2006.01272.x
Balls MJ, Bødker R, Thomas CJ, Kisinza W, Msangeni HA, Lindsay SW: Effect of topography on the risk of malaria infection in the Usambara Mountains, Tanzania. Trans R Soc Trop Med Hyg 2004, 98: 400–408. 10.1016/j.trstmh.2003.11.005
Ahrends A, Burgess ND, Milledge SAH, Bulling MT, Fisher B, Smart JCR, Clarke GP, Mhoro BE, Lewis SL: Predictable waves of sequential forest degradation and biodiversity loss spreading from an African city. Proc Natl Acad Sci 2010, 107: 14556–14561. 10.1073/pnas.0914471107
Bayon G, Dennielou B, Etoubleau J, Ponzevera E, Toucanne S, Bermell S: Intensifying weathering and land use in iron age Central Africa. Science 2012, 335: 1219–1222. 10.1126/science.1215400
Putz FE, Zuidema PA, Synnott T, Peña-Claros M, Pinard MA, Sheil D, Vanclay JK, Sist P, Gourlet-Fleury S, Griscom B, Palmer J, Zagt R: Sustaining conservation values in selectively logged tropical forests: the attained and the attainable. Conserv Lett 2012, 5: 296–303. 10.1111/j.1755-263X.2012.00242.x
Topp-Jørgensen E, Poulsen MK, Lund JF, Massao JF: Community-based monitoring of natural resource use and forest quality in Montane Forests and Miombo Woodlands of Tanzania. Biodivers Conserv 2005, 14: 2653–2677. 10.1007/s10531-005-8399-5
Silva JNM, de Carvalho JOP, Lopes JCA, de Almeida BF, Costa DHM, de Oliveira LC, Vanclay JK, Skovsgaard JP: Growth and yield of a tropical rain forest in the Brazilian Amazon 13 years after logging. Forest Ecol Manag 1995, 71: 267–274. 10.1016/0378-1127(94)06106-S
Govender N, Trollope WSW, Van Wilgen BW: The effect of fire season, fire frequency, rainfall and management on fire intensity in savanna vegetation in South Africa. J Appl Ecol 2006, 43: 748–758. 10.1111/j.1365-2664.2006.01184.x
Hoffmann WA: The effects of fire and cover on seedling establishment in a Neotropical Savanna. J Ecol 1996, 84: 383–393. 10.2307/2261200
Chhatre A, Agrawal A: Trade-offs and synergies between carbon storage and livelihood benefits from forest commons. Proc Natl Acad Sci 2009, 106: 17667–17670. 10.1073/pnas.0905308106
Chhatre A, Agrawal A: Forest commons and local enforcement. Proc Natl Acad Sci 2008, 105: 13286–13291. 10.1073/pnas.0803399105
Mbwambo L, Eid T, Malimbwi RE, Zahabu E, Kajembe GC, Luoga E: Impact of decentralised forest management on forest resource conditions in Tanzania. Forests, Trees and Livelihoods 2012, 1–17.
Fisher B, Lewis SL, Burgess ND, Malimbwi RE, Munishi PK, Swetnam RD, Kerry Turner R, Willcock S, Balmford A: Implementation and opportunity costs of reducing deforestation and forest degradation in Tanzania. Nat Clim Chang 2011, 1: 161–164. 10.1038/nclimate1119
Doherty RM, Sitch S, Smith B, Lewis SL, Thornton PK: Implications of future climate and atmospheric CO2 content for regional biogeochemistry, biogeography and ecosystem services across East Africa. Glob Chang Biol 2009, 16: 617–640.
Willcock S: Long-term Changes in Land Cover and Carbon Storage in Tanzania, East Africa. University of Leeds, School of Geography 2012.
NBS: Tanzania Census 2002 - Analytical Report Volume X. Dar es Salaam: National Bureau of Statistics, Ministry of Planning, Economy and Empowerment; 2006:193.
Clark DA, Piper SC, Keeling CD, Clark DB: Tropical rain forest tree growth and atmospheric carbon dynamics linked to interannual temperature variation during 1984–2000. Proc Natl Acad Sci U S A 2003, 100: 5852–5857. 10.1073/pnas.0935903100
Raich JW, Russell AE, Kitayama K, Parton WJ, Vitousek PM: Temperature influences carbon accumulation in moist tropical forests. Ecology 2006, 87: 76–87. 10.1890/05-0023
Phillips OL, Malhi Y, Vinceti B, Baker T, Lewis SL, Higuchi N, Laurance WF, Vargas PN, Martinez RV, Laurance S, Ferreira LV, Stern M, Brown S, Grace J: Changes in growth of tropical forests: evaluating potential biases. Ecol Appl 2002, 12: 576–587. 10.1890/1051-0761(2002)012[0576:CIGOTF]2.0.CO;2
Enquist BJ, Niklas KJ: Invariant scaling relations across tree-dominated communities. Nature 2001, 410: 655–660. 10.1038/35070500
Li H-T, Han X-G, Wu J-G: Lack of evidence for 3/4 scaling of metabolism in terrestrial plants. J Integr Plant Biol 2005, 47: 1173–1183. 10.1111/j.1744-7909.2005.00167.x
Enquist BJ, West GB, Brown JH: Extensions and evaluations of a general quantitative theory of forest structure and dynamics. Proc Natl Acad Sci 2009, 106: 7046–7051. 10.1073/pnas.0812303106
Goldstein ML, Morris SA, Yen GG: Problems with fitting to the power-law distribution. Eur Phys J B Condens Matter Complex Syst 2004, 41: 255–258. 10.1140/epjb/e2004-00316-5
Zanne AE, Lopez-Gonzalez G, Coomes DA, Ilic J, Jansen S, Lewis SL, Miller RB, Swenson NG, Wiemann MC, Chave J: Global wood density database. 2009.http://hdl.handle.net/10255/dryad.235 [Accessed 5/12/2008]
Baker TR, Phillips OL, Malhi Y, Almeida S, Arroyo L, Di Fiore A, Erwin T, Killeen TJ, Laurance SG, Laurance WF, Lewis SL, Lloyd J, Monteagudo A, Neill DA, Patiño S, Pitman NCA, Silva JNM, Vásquez Martínez R: Variation in wood density determines spatial patterns in Amazonian forest biomass. Glob Chang Biol 2004, 10: 545–562. 10.1111/j.1365-2486.2004.00751.x
Djomo AN, Ibrahima A, Saborowski J, Gravenhorst G: Allometric equations for biomass estimations in Cameroon and pan moist tropical equations including biomass data from Africa. Forest Ecol Manag 2010, 260: 1873–1885. 10.1016/j.foreco.2010.08.034
Chave J, Andalo C, Brown S, Cairns MA, Chambers JQ, Eamus D, Folster H, Fromard F, Higuchi N, Kira T, Lescure J-P, Nelson BW, Ogawa H, Puig H, Riéra B, Yamakura T: Tree allometry and improved estimation of carbon stocks and balance in tropical forests. Oecologia 2005, 145: 87–99. 10.1007/s00442-005-0100-x
Platts PJ: Spatial Modelling, Phytogeography and Conservation in the Eastern Arc Mountains of Tanzania and Kenya. University of York, Environment Department 2012.
Swetnam RD: Historical logging in protected areas, E. Tanzania. V2. Cambridge, UK: Zoology Department, Cambridge University; 2011.
IUCN, UNEP-WCMC: The World Database on Protected Areas (WDPA). Cambridge, UK: UNEP- WCMC; 2010.http://www.protectedplanet.net Available at: [Accessed 03/05/2010)]
Zomer RJ, Trabucco A, Bossio DA, Verchot LV: Climate change mitigation: A spatial analysis of global land suitability for clean development mechanism afforestation and reforestation. Agr Ecosyst Environ 2008, 126: 67–80. 10.1016/j.agee.2008.01.014
Tropical Rainfall Measuring Mission http://trmm.gsfc.nasa.gov/
Hijmans RJ, Cameron SE, Parra JL, Jones PG, Jarvis A: Very high resolution interpolated climate surfaces for global land areas. Int J Climatol 2005, 25: 1965–1978. 10.1002/joc.1276
Hole-filled SRTM for the globe Version 4, from the CGIAR-CSI SRTM 90m Database http://srtm.csi.cgiar.org
Wind Speed At 50 m Above The Surface Of The Earth http://eosweb.larc.nasa.gov/sse/
Perez R, Ineichen P, Moore K, Kmiecik M, Chain C, George R, Vignola F: A new operational model for satellite-derived irradiances: description and validation. Sol Energy 2002, 73: 307–317. 10.1016/S0038-092X(02)00122-6
Low Resolution Solar Data http://www.nrel.gov/gis/
Batjes NH: SOTER-based soil parameter estimates for Southern Africa. Volume 4. Wageningen: ISRIC - World Soil Information; 2004:27.
ISRIC: SOTER and WISE-based soil property estimates for Southern Africa. 2010.http://www.isric.org/ Available at [Accessed 17/2/2010]
Roy DP, Jin Y, Lewis PE, Justice CO: Prototyping a global algorithm for systematic fire-affected area mapping using MODIS time series data. Remote Sens Environ 2005, 97: 137–162. 10.1016/j.rse.2005.04.007
Houghton RA, Lawrence KT, Hackler JL, Brown S: The spatial distribution of forest biomass in the Brazilian Amazon: a comparison of estimates. Glob Chang Biol 2001, 7: 731–746. 10.1046/j.1365-2486.2001.00426.x
Clark DB, Clark DA: Landscape-scale variation in forest structure and biomass in a tropical rain forest. Forest Ecol Manag 2000, 137: 185–198. 10.1016/S0378-1127(99)00327-8
Dormann CF, McPherson JM, Araújo MB, Bivand R, Bolliger J, Carl G, Davies RG, Hirzel A, Jetz W, Daniel Kissling W, Kühn I, Ohlemüller R, Peres-Neto PR, Reineking B, Schröder B, Schurr FM, Wilson R: Methods to account for spatial autocorrelation in the analysis of species distributional data: a review. Ecography 2007, 30: 609–628. 10.1111/j.2007.0906-7590.05171.x
R Development Core Team: R: A Language and Environment for Statistical Computing. In R Foundation for Statistical Computing. Vienna, Austria: R Foundation for Statistical Computing; 2010.
ESRI: ArcMap 9.3.1. 1999–2009.
Venables WN, Ripley BD: Modern Applied Statistics with S. Oxford, UK: Springer; 2002.
Varma S, Simon R: Bias in error estimation when using cross-validation for model selection. BMC Bioinforma 2006, 7: 91. 10.1186/1471-2105-7-91
Johnson JB, Omland KS: Model selection in ecology and evolution. Trends Ecol Evol 2004, 19: 101–108. 10.1016/j.tree.2003.10.013
Chong I-G, Jun C-H: Performance of some variable selection methods when multicollinearity is present. Chemom Intell Lab Syst 2005, 78: 103–112. 10.1016/j.chemolab.2004.12.011
Asner G, Flint Hughes R, Varga T, Knapp D, Kennedy-Bowdoin T: Environmental and Biotic Controls over Aboveground Biomass Throughout a Tropical Rain Forest. Ecosystems 2009, 12: 261–278. 10.1007/s10021-008-9221-5
Lloyd J, Farquhar GD: Effects of rising temperatures and CO 2 on the physiology of tropical forest trees. Phil Trans Roy Soc B Biol Sci 2008, 363: 1811–1817. 10.1098/rstb.2007.0032
Graham EA, Mulkey SS, Kitajima K, Phillips NG, Wright SJ: Cloud cover limits net CO2 uptake and growth of a rainforest tree during tropical rainy seasons. Proc Natl Acad Sci 2003, 100: 572–576. 10.1073/pnas.0133045100
Maggini R, Lehmann A, Zimmermann NE, Guisan A: Improving generalized regression analysis for the spatial prediction of forest communities. J Biogeogr 2006, 33: 1729–1749. 10.1111/j.1365-2699.2006.01465.x
Burgess ND, Butynski TM, Cordeiro NJ, Doggart NH, Fjeldsa J, Howell KM, Kilahama FB, Loader SP, Lovett JC, Mbilinyi B, Menegon M, Moyer DC, Nashanda E, Perkin A, Rovero F, Stanley WT, Stuart SN: The biological importance of the Eastern Arc Mountains of Tanzania and Kenya. Biol Conserv 2007, 134: 209–231. 10.1016/j.biocon.2006.08.015
This study is part of the Valuing the Arc research programme (http://valuingthearc.org/) funded by the Leverhulme Trust (http://www.leverhulme.ac.uk/). Manuscript preparation and later analyses took place under the ‘Which Ecosystem Service Models Best Capture the Needs of the Rural Poor?’ project (WISER; NE/L001322/1), funded with support from the United Kingdom’s Ecosystem Services for Poverty Alleviation program (ESPA; http://www.espa.ac.uk). ESPA receives its funding from the Department for International Development (DFID), the Economic and Social Research Council (ESRC) and the Natural Environment Research Council (NERC). SLL was funded by a Royal Society University Research Fellowship; SW additionally by the Stokenchurch Charity; OLP was supported by an Advanced Grant from the European Research Council and is a Royal Society-Wolfson Research Award holder. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. We thank the Tanzanian Commission for Science and Technology (COSTECH), the Tanzanian Wildlife Institute (TAWIRI) and the Sokoine University of Agriculture for their support of this work, as well as all the field assistants involved. Furthermore, we would like to thank the two anonymous reviewers, whose comments and insight vastly improved the manuscript.
The authors declare that they have no competing interests.
The majority of this jointly-authored publication was led by SW. Contributions to the collaborative dataset came from PJP, AA, ND, KD, EF, JG, JH, KH, ARM, BM, PKTM, NO, EJTJ, AM, SW and RDS. The analysis was performed by SW, supervised by OLP and SLL. The manuscript was prepared by SW, with assistance from OLP, SLL, AB, PP, NDD and RM. All authors read and approved the final manuscript.
An erratum to this article can be found online at http://dx.doi.org/10.1186/s13021-017-0088-7.
Electronic supplementary material
Additional file 2: Figure S1: The spatial variation of WSG in tree-dominated land cover categories within the study area (a), with upper (b) and lower (c) pixel based 95% CI. See text for details on methods. (DOC 410 KB)
Additional file 3: Figure S2: The most influential, significant influential variables on WSG (a and b), the intercept of the power law relationship (c and d), and the gradient of the power law relationship (e and f). Dashed red lines indicate 95% CI. (DOC 77 KB)
Additional file 4: Figure S3: The spatial variation in the intercept of the power law relationship (a proxy measure for potential stem density) in tree dominated land cover categories within the study area (a), with upper (b) and lower (c) pixel based 95% CI. See text for details on methods. (DOC 273 KB)
Additional file 5: Figure S4: The spatial variation in the gradient of the power law relationship (a proxy measure for the proportion of larger stems) in tree-dominated land cover categories within the study area (a), with upper (b) and lower (c) pixel based 95% CI. See text for details on methods. (DOC 254 KB)
Additional file 6: Figure S5: The 2nd order interactions relating my carbon storage derivatives (wood specific gravity, maximum canopy height, the intercept of the power law relationship, and the gradient of the power law relationship [shown here as WSG, height, intercept, and gradient respectively]) to aboveground live carbon storage. Dashed red lines indicate 95% CI. (DOC 91 KB)
Additional file 7: Figure S6: The effect of MAT on tree height for a range of DBH. The data (points) correspond to DBH ranges whereas the Gompertz model fits (solid lines) illustrate the relationship for mid-point of this range only. Dotted lines represent the 95CI of the model fits. (DOC 122 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Willcock, S., Phillips, O.L., Platts, P.J. et al. Quantifying and understanding carbon storage and sequestration within the Eastern Arc Mountains of Tanzania, a tropical biodiversity hotspot. Carbon Balance Manage 9, 2 (2014). https://doi.org/10.1186/1750-0680-9-2