- Open Access
Human and environmental controls over aboveground carbon storage in Madagascar
© Asner et al; licensee BioMed Central Ltd. 2012
Received: 2 November 2011
Accepted: 30 January 2012
Published: 30 January 2012
Accurate, high-resolution mapping of aboveground carbon density (ACD, Mg C ha-1) could provide insight into human and environmental controls over ecosystem state and functioning, and could support conservation and climate policy development. However, mapping ACD has proven challenging, particularly in spatially complex regions harboring a mosaic of land use activities, or in remote montane areas that are difficult to access and poorly understood ecologically. Using a combination of field measurements, airborne Light Detection and Ranging (LiDAR) and satellite data, we present the first large-scale, high-resolution estimates of aboveground carbon stocks in Madagascar.
We found that elevation and the fraction of photosynthetic vegetation (PV) cover, analyzed throughout forests of widely varying structure and condition, account for 27-67% of the spatial variation in ACD. This finding facilitated spatial extrapolation of LiDAR-based carbon estimates to a total of 2,372,680 ha using satellite data. Remote, humid sub-montane forests harbored the highest carbon densities, while ACD was suppressed in dry spiny forests and in montane humid ecosystems, as well as in most lowland areas with heightened human activity. Independent of human activity, aboveground carbon stocks were subject to strong physiographic controls expressed through variation in tropical forest canopy structure measured using airborne LiDAR.
High-resolution mapping of carbon stocks is possible in remote regions, with or without human activity, and thus carbon monitoring can be brought to highly endangered Malagasy forests as a climate-change mitigation and biological conservation strategy.
The spatial distribution of carbon stored in the aboveground tissues of vegetation - also known as aboveground carbon density (ACD; units of Mg C ha-1) - is a time-integrated expression of ecological and land-use processes ranging from photosynthesis and nutrient cycling to disturbance and climate change. Spatial variation in ACD is also the largest source of uncertainty in monitoring carbon emissions for voluntary carbon offset markets and for developing international action for Reduced Emissions from Deforestation and Forest Degradation (REDD) at national and sub-national levels [1–4]. Mapping the geographic patterns of carbon storage is thus a high priority in scientific, conservation, and resource-management communities.
Carbon mapping efforts have proven challenging for a variety of reasons. Field inventory plots are critically important at local scales, but they are time consuming, costly, and limited by accessibility. As a result, they usually do not capture the variation in ACD that exists throughout the environment . This is particularly true in human-dominated landscapes, or in remote areas that may harbor remaining forests, and these two scenarios have become ubiquitous in many tropical regions . Madagascar is an important example of these challenges; high rates of deforestation and degradation have transformed millions of hectares of dry, mesic and humid tropical forest, leaving highly fragmented landscapes in most regions [7, 8]. Remaining forests cover up to 15% of Madagascar, and thus the vast majority of landscapes are partially-vegetated, human-dominated systems . The spatial heterogeneity of vegetation cover and structure in these landscapes has resulted in a poor overall knowledge of their role in storing carbon.
The remaining forests of Madagascar are also often remote, whether by distance or by difficult terrain, and they are among the last strongholds of the unique forest flora and fauna remaining on the island. These remote forests appear partially to fully intact in satellite imagery [9, 10], but satellite sensors do not resolve their carbon stocks at a spatial scale sufficient to understand ecological controls over those stocks. As a result, we know little about controls over carbon storage in Malagasy forests - information that is critically needed to support their inclusion in carbon retention strategies to conserve them.
Mapping ACD beyond the reach of field plots is an important step toward resolving regional carbon dynamics in Madagascar and elsewhere. Airborne Light Detection and Ranging (LiDAR), when combined with field plots, has provided spatially contiguous, high-resolution ACD estimates of temperate and tropical forests [e.g., [11–13]]. Airborne systems can readily acquire thousands of hectares of data per day, changing our understanding of environmental controls over carbon stocks. Airborne LiDAR-based mapping of tropical forest ACD has also proven highly precise and accurate, with errors recently becoming indistinguishable from those derived from field measurements .
Despite the carbon mapping accuracy, resolution, and extensive coverage provided by airborne LiDAR, these approaches also reach geographic limitations due to cost, so additional methods are required to extend LiDAR-based carbon estimates to even larger regions. Importantly, the pronounced heterogeneity of carbon stocks in Madagascar requires high spatial resolution when extending to the regional or whole-country scale. One approach to extrapolating outside the LiDAR coverage is to stratify the region into narrow vegetation classes, and then to apply additional analyses to account for vegetation losses and gains from deforestation, degradation and land abandonment [15, 16]. This has proven useful in regions already containing detailed vegetation maps derived from satellite or other sources, resulting in regionally-integrated uncertainties in total carbon storage of just a few percent or less . However, such highly stratified maps of vegetation type do not exist for many areas of the tropics [2, 18], and other means are necessary to extend sampling-based estimates of ACD - whether taken from field, aircraft or satellite measurements - to broader geographic scales. Again, Madagascar serves as an example of the challenges faced when mapping carbon stocks in most tropical regions, where vegetation maps may not exist or may only be based on coarse representations of geology [e.g., ]. A more detailed, top-down approach utilizing satellite imagery is required; however, the most widely available optical satellite images from NASA Landsat or MODIS (Moderate Resolution Imaging Spectrometer) do not directly measure carbon stocks. Rather, optical sensors are most sensitive to vegetation type and cover , and so spatial extrapolation of LiDAR-based biomass estimates using optical satellite data alone is subject to substantial uncertainty. Synthetic aperture radar (SAR), which is sensitive to aspects of canopy texture and structure, can be related to biomass stocks [21, 22], but direct SAR-based metrics of forest biomass are limited by saturation of the signal in vegetation harboring more than 50-70 Mg C ha-1, depending upon radar frequency [23, 24]. Another option is to evaluate how carbon stocks vary with topography and other geologic drivers mapped with SAR, but relationships between topography or geology and ACD have rarely been explored [25, 26].
Here, we present the first large-scale, high-resolution mapping estimates of aboveground carbon stocks in one northern and one southern region of Madagascar. Using airborne LiDAR-based measurements of ACD, calibrated to a network of field plots, we quantitatively assessed which topographic and land-cover factors best predicted ACD across environmental and land-use gradients in each region. We then used this information, combined with spaceborne optical and SAR data, to extend field and LiDAR carbon estimates to the regional level. Our underlying goal was to understand both human and environmental controls over the carbon landscape in some of the most remote portions of Madagascar.
Results and discussion
Field Calibration of Airborne LiDAR
Landscape Controls over LiDAR-derived ACD
In both the N and S regions, aboveground carbon stocks reached their highest median levels at mid-elevation (Figure 4). In the north, we found that ACD peaked in the 1350-1650 m zone, whereas highest median values were closer to 700 m in the south (Appendix 3, Additional file 1). Above these zones of highest carbon storage, both ACD median and variance decreased steadily to much lower and narrower ranges in montane conditions. The sharpest high-elevation decline in ACD occurred above 1950 m in the N region, where forests give way to short-statured Philippia shrublands . In the S region, the decrease in ACD at higher elevations was less pronounced. We note, however, that ACD levels in the S sub-montane (e.g., 1100-1500 m) and montane (> 1500 m) zones were lower and less variable than in the N forests at comparable elevations. Importantly, most ACD distributions were non-normal, suggesting that without very large sample sizes, or carefully distributed sampling, field plots alone may not accurately resolve carbon stock variation.
Deforestation, disturbance and regrowth all imparted differences in standing carbon stocks (Appendix 3, Additional file 1). In the N region, secondary regrowth from deforestation of at least 5 and 10 years of age resulted in mean (± s.d.) ACD of 39.4 (29.4) and 52.9 (49.5) Mg C ha-1, respectively. Forested areas subjected to disturbance (also known as degradation), as mapped with CLASlite, contained 70.9 (± 56.6) Mg C ha-1. In the S region, carbon stocks in forest regrowth varied slightly by elevation (Appendix 3, Additional file 1), but were generally between 16 and 49 Mg C ha-1, depending upon age and whether they were deforested or disturbed prior to assessment. In sum, aboveground carbon stocks in forest regrowth were extremely variable, but consistently much lower than in intact forest.
Regional Sources of Carbon Stock Variation
The relative error between high-resolution LiDAR-scale ACD maps and regional ACD maps was calculated as the percent difference between the square root of the median squared residual value (predicted ACD - observed ACD) and mean observed ACD. We used the median rather than mean squared residual value to compensate for highly right-skewed distributions of pixel level error. Thus, our pixel-level uncertainties reflect median uncertainties. We estimated a pixel-level uncertainty of 35% and 10% in the N and S regions, respectively, at 1 ha resolution. Greater uncertainty in the N region is caused by generally weaker relationships between environmental variables and measured ACD (Figure 5). Nonetheless, a 35% uncertainty at the pixel level is on par with errors inherent to field inventory plots [30, 31], and when averaged over much larger jurisdictional scales of thousands to millions of hectares, the overall uncertainty drops to extremely low values [3, 17].
Additional comparisons between regionally-extrapolated ACD and high-resolution forest cover were made by overlaying regional ACD maps on Google Earth imagery (Figure S2, Additional file 1). In one southern landscape, the largest tree crowns were clearly aligned with areas of high ACD (reds), whereas areas containing small-statured vegetation (naturally or from land use) aligned with low-to-moderate ACD (yellow-green), and no ACD (blues) estimates. In one portion of the southern montane landscape, the abrupt orographically-mediated transition from high biomass windward forests (red) to lower biomass leeward forests (orange) was evident (Figure S2, Additional file 1). This break was so distinct, it was expressed over just a few hundred meters of distance. Similar extreme gradients can be found on other tall oceanic islands where orographic effects are common . Close inspection against Google imagery also revealed that suppressed ACD levels along slopes (greens-blues; Figure S2, Additional file 1) were caused either by natural landslides or human activities.
As LiDAR coverage was statistically reduced, we found that median ACD remained stable for most vegetation classes - within 10% of the true median value - until the coverage declined below ~ 300 ha (Figure S3, Additional file 1). However, for some classes, at least 2000 ha of LiDAR sampling was required to be within 10% of the true median. Conservatively, median ACD values reported for classes with > 500 ha of LiDAR coverage can be considered to have very high confidence. These large classes account for 90% and 88% of the core study areas in the N and S regions, respectively (Appendix 3, Additional file 1). These results point to the importance of obtaining thousands of samples in order to derive the true median and distribution of ACD throughout the landscape.
Terrestrial carbon landscapes depict the integrated effects of many ecological and socio-economic processes occurring at spatial scales ranging from localized human activities (e.g., wood gathering) to synoptic climate-driven processes affecting vegetation growth and carbon storage. As a result, carbon maps are most informative when created over large environmental gradients while maintaining high spatial resolution. With a combination of airborne LiDAR, field calibration plots and satellite data, we mapped aboveground carbon stocks throughout two remote regions of Madagascar covering a total of 2,372,680 ha. Despite the widespread evidence of human activities throughout the regions, we found that large-scale natural controls over aboveground carbon density were driven primarily by physiography and vegetation cover. These results strongly suggest that environmental surrogates can be used to extend direct observations made by airborne LiDAR to much larger geographic scales. Our results support other work on biomass surrogates at coarser resolutions [21, 35], but the high-resolution regional stratification and airborne LiDAR approach presented here yields carbon storage detail that can be validated point by point through a large region. The global and regional approaches are complimentary and could be combined in future work to improve carbon mapping at national and global levels.
Here our regionally mapped carbon stocks are the result of using very high resolution airborne LiDAR to understand the physiographic and vegetation cover controls, and then applying that knowledge outside of the LiDAR coverage using high resolution satellite data. Our previous efforts along these lines have relied on habitat classification maps, often provided by government agencies, but in this case, comparable effectiveness was achieved using freely available SRTM and Landsat imagery. The results thus suggest that verifiable, high-resolution carbon assessments can be accomplished in remote regions of the tropics. The high resolution and accuracy of the results may support new efforts to sustain and expand Malagasy forests using current voluntary and planned compliance carbon emission offset strategies.
Regional Pre-flight Stratification
The pre-flight stratification maps were used to direct airborne LiDAR sampling. Specifically, LiDAR polygons were placed over each study region to maximize coverage of each pre-flight stratum while attending to aircraft logistical constraints imposed by mountainous terrain and winds. A total of 12 and 10 polygons ranging in size from 2,000 ha to 12,000 ha were selected for the N and S regions, respectively (Figure 1). We used the Carnegie Airborne Observatory (CAO) to collect the LiDAR data . The flights were conducted at 2000 m above ground level, providing LiDAR sampling at 1.12 m LiDAR spot spacing with a 28-degree field of view, 50 kHz pulse repetition frequency and 50% overlap between adjacent flight lines within each mapping polygon. The aircraft maintained a ground speed of 95 knots or less. LiDAR spatial errors are less than 0.15 m vertically and 0.36 m horizontally (RMSE). The total LiDAR coverage was 56,921 ha (19.6% of the core study area) and 70,780 ha (14.3%) in the N and S regions, respectively.
Field Plot Measurements
We estimated ground-based ACD in 83 circular plots of up to 30 m radius, or 0.28 ha (Figure 1). Of these, 19 were located in N humid forests, 28 in S humid forests, and 36 in S spiny forests. We used the plots to calibrate the LiDAR data to ACD, and thus we placed the plots in a wide range of conditions, from dry to humid forests and from nearly deforested land to degraded forests to intact, closed canopy forest. This approach reduced the number of plots required overall, while providing the biophysical variation needed to calibrate the LiDAR. We used a nested sampling design in each plot, which differed depending on forest structural complexity and thus climactic zone. For humid forest plots, we included all trees > 50 cm in diameter at breast height (dbh, 1.3 m from the ground) or above buttress inside a radius of 30 m from the plot center, all trees > 20 cm dbh inside a radius of 14 m, and all trees > 5 cm dbh inside a radius of 4 m. In spiny forests, we included all trees > 5 cm dbh inside a radius of 20 m and all trees > 0 cm dbh inside a radius of 5 m. We located each plot center using a global positioning system (GPS) with a survey-grade receiver (GS50+, Leica Geosystems, St. Fallen, Switzerland). The GPS data were later differentially corrected, in most cases yielding < 1 m positional error. For some plots in the S humid forests, Leica GPS reception was poor and thus we utilized a Garmin CS+ receiver (Garmin Ltd., Olathe, Kansas) with the point averaging function enabled for several hours. When comparing Garmin and Leica data, we found that the loss of accuracy with Garmin data was 1-4 m. Allometric methods for estimating plot-level carbon stocks made use of both general  and recent, locally-developed models , as detailed in Appendix 1 (Additional file 1).
where x is MCH, y is ACD, and a, b, and k are model parameters. A non-arithmetic error term (xkε) was used to account for heteroscedasticity that typifies MCH-to-ACD relationships [e.g., ]. This approach is analogous to fitting a linear model to log-transformed variables, but avoids the need for back-transformation . The model was fit using maximum likelihood in the R programming language (v.9.2, R Development Core Team 2009). We conducted additional analyses to determine the dependency of our LiDAR-to-ACD model on the number of field plots (below under Error Analyses). We also considered the dependency of our LiDAR-to-ACD model on the spatial resolution at which the equation was constructed (Appendix 2, Additional file 1).
LiDAR Carbon and Terrain Analysis
We tested the degree to which terrain and forest cover variables were related to LiDAR-derived ACD. Our goal was to develop an improved method to scale up the LiDAR-based data to the regional level using freely available satellite imagery. To that end, we examined relationships between ACD and elevation data, the latter provided by NASA's Shuttle Radar Topography Mission (SRTM; 90 m resolution), as well as CLASlite forest cover data derived from Landsat imagery (30 m resolution). To improve signal via averaging, and to reduce errors due to misalignment among diverse image data types, we convolved all variables to 1 ha spatial resolution prior to analysis using a spatial averaging filter. Within areas under airborne LiDAR coverage, we considered several SRTM-derived variables, including elevation, slope, aspect, and height above nearest stream (i.e., the vertical distance from a point on the landscape to an interpolated water table, also known as a relative elevation model or REM). From CLASlite, we considered the fractional cover of photosynthetic vegetation (PV) and the fractional cover of bare soil (described in detail by Asner et al. 2009). We considered the influence on ACD of all variables-both individually and in combination-at 1 ha resolution using correlation and ordinary linear regression analyses with the R programming language (v.9.2, R Development Core Team 2009).
Final Regional Stratification and Carbon Mapping
Final stratification was intended to utilize information obtained from the LiDAR-based environmental controls analysis, as well as the deforestation and forest disturbance results from CLASlite (Figure 9). We first segmented the region into forest, non-forest, and regrowth, derived from CLASlite analysis of Landsat images from four periods; 1990, 2000, 2005, and 2010. The CLASlite analysis provided a detailed representation of 2010 forest condition, enhanced by the consideration of past land-cover change. The broader classes of forest, non-forest and regrowth were further stratified to capture regional ACD variation. This process was guided by existing scientific knowledge of the ecology and terrain of Madagascar, as well as statistical considerations to limit over-stratification.
We stratified the forest class using the SRTM- and CLASlite-derived inputs that best predicted ACD in forested areas at the LiDAR scale. As presented in the results (Landscape Controls over LiDAR-derived ACD), the best predictors were SRTM elevation and CLASlite PV cover fraction. The nature of the influence of these variables on LiDAR-derived ACD suggested that a stratified approach to carbon mapping was preferable to modeling ACD in continuous fashion (see results). Thus, in the N region, we parsed the landscape into 150 m vertical increments, creating 13 elevation strata. In the S region, we created a total of 11 elevation strata with 100 m elevation increments below 500 m (characterized by dry, spiny forests and flatter terrain), and 200 m increments above 500 m in elevation (steeper, variable terrain containing the region's remaining humid forests). The two regions were further stratified by PV cover fraction with thresholds of 80-90% and 91-100% PV in the N region and 60-70%, 71-80%, 81-90%, and 91-100% in the S Region. Our approach to the S region employed a lower PV threshold to better capture dry spiny forests, which often contain lower PV fractions due to naturally open canopies and space between trees and other woody plants. The intersection of elevation and PV cover fraction strata produced 26 forested classes in the N region and 44 classes in the S region.
In both regions, we used CLASlite  to stratify deforested and disturbed lands (e.g., areas entering a non-forest category following previous mapping as forest cover), as well as those in regrowth following deforestation or disturbance. Here, we define disturbance as the diffuse thinning of forest cover. We split forest regrowth following deforestation into categories of greater than 5 or 10 years since the most recent detection of deforestation. We created another class for regrowth following forest disturbance detected in any year. Areas not in forest and not undergoing regrowth following deforestation or forest disturbance made up the broad non-forest class, which we stratified into five 20% intervals of PV cover fraction. In the S region, all regrowth and non-forest classes were further parsed by a 500 m elevation threshold that generally corresponds to the separation of lowland, dry spiny forests and mid-elevation humid forests. Combined, this yielded 34 total classes in the N region and 52 in the S region. Using this final stratification, we mapped ACD at 1.0 ha resolution by assigning the LiDAR-derived median carbon value for each class to the area occupied by that class across the regional maps.
We determined the sensitivity of our LiDAR ACD model to the number of field calibration plots used by fitting simulated calibration models with randomized multiples of 4 plots (~ 5% of the data), from 4 to 76 plots (up to 90% of 83 total plots) . We repeated this simulated calibration 1000 times to yield a mean regression standard error of the estimate (SEE) at every interval of 4 field plots. For each run, we then assessed the ability of the model to predict (as measured by the RMSE) the carbon density of an additional set of 4 randomized field plots not included as part of the simulated calibration.
To determine the amount of LiDAR coverage needed to accurately characterize inter-class variation in ACD, we simulated reduced LiDAR coverage by successively excluding portions of all LiDAR flight boxes. To do so, we reduced the widths of all boxes beginning at a minimum width of 30 m (a single pixel of our LiDAR-derived carbon map), doubling repeatedly up to 3840 m width (an approximation of the largest width that was divisible by the pixel size of 30 m). With each successively more limited regional LiDAR sample, we recomputed median carbon for each vegetation class in both the N and S study areas, and compared the resulting estimates to those determined using the full LiDAR coverage.
We thank the GoodPlanet Foundation and WWF-Madagascar for vegetation mapping and field measurement assistance, and additional logistical support. We thank G. Powell and N. O'Connor for decision support and logistical help, and C. Grinand and three anonymous reviewers for providing helpful comments on the manuscript. The Carnegie Airborne Observatory flight campaign to the region was supported by the Andrew Mellon Foundation, the Carnegie Institution and Air France. The Carnegie Airborne Observatory is made possible by the Grantham Foundation for the Protection of the Environment, the Gordon and Betty Moore Foundation, the W.M. Keck Foundation, and William Hearst III.
- Angelsen A: Moving Ahead with REDD: issues, options and implications. Bogor, Indonesia: Center for International Forestry Research (CIFOR); 2008.Google Scholar
- Pelletier J, Ramankutty N, Potvin C: Diagnosing the uncertainty and detectability of emission reductions for REDD + under current capabilities: an example for Panama. Environmental Research Letters 2011, 6: 024005. 10.1088/1748-9326/6/2/024005View ArticleGoogle Scholar
- Asner GP: Painting the world REDD: addressing scientific barriers to monitoring emissions from tropical forests. Environmental Research Letters 2011, 6: 021002. 10.1088/1748-9326/6/2/021002View ArticleGoogle Scholar
- Avitabile V, Herold M, Henry M, Schmullius C: Mapping biomass with remote sensing: a comparison of methods for the case study of Uganda. Carbon Balance and Management 2011, 6: 7. 10.1186/1750-0680-6-7View ArticleGoogle Scholar
- Goetz S, Baccini A, Laporte N, Johns T, Walker W, Kellndorfer J, Houghton R, Sun M: Mapping and monitoring carbon stocks with satellite observations: a comparison of methods. Carbon Balance and Management 2009., 4: doi:10.1186/1750–0680–1184–1182 doi:10.1186/1750-0680-1184-1182Google Scholar
- Asner GP, Rudel TK, Aide TM, Defries R, Emerson RE: A Contemporary Assessment of Change in Humid Tropical Forests. Conservation Biology 2009, 23: 1386–1395. 10.1111/j.1523-1739.2009.01333.xView ArticleGoogle Scholar
- Brown KA, Gurevitch J: Long-term impacts of logging on forest diversity in Madagascar. Proceedings of the National Academy of Sciences 2004, 101: 6045–6049. 10.1073/pnas.0401456101View ArticleGoogle Scholar
- Green GM, Sussman RW: Deforestation history of the eastern rain forests of Madagascar from satellite images. Science 1990, 248: 212–215. 10.1126/science.248.4952.212View ArticleGoogle Scholar
- Harper GJ, Steininger MK, Tucker CJ, Juhn D, Hawkins F: Fifty years of deforestation and forest fragmentation in Madagascar. Environmental Conservation 2007, 34: 325–333.View ArticleGoogle Scholar
- Ingram JC, Dawson TP, Whittaker RJ: Mapping tropical forest structure in southeastern Madagascar using remote sensing and artificial neural networks. Remote Sensing of Environment 2005, 94: 491–507. 10.1016/j.rse.2004.12.001View ArticleGoogle Scholar
- Gonzalez P, Asner GP, Battles JJ, Lefsky MA, Waring KM, Palace M: Forest carbon densities and uncertainties from Lidar, QuickBird, and field measurements in California. Remote Sensing of Environment 2010, 114: 1561–1575. 10.1016/j.rse.2010.02.011View ArticleGoogle Scholar
- Lefsky MA, Cohen WB, Harding DJ, Parker GG, Acker SA, Gower ST: Lidar remote sensing of above-ground biomass in three biomes. Global Ecology and Biogeography 2002, 11: 393–399. 10.1046/j.1466-822x.2002.00303.xView ArticleGoogle Scholar
- Drake JB, Dubayah RO, Clark DB, Knox RG, Blair JB, Hofton MA, Chazdon RL, Weishampel JF, Prince SD: Estimation of tropical forest structural characteristics using large-footprint lidar. Remote Sensing of Environment 2002, 79: 305–319. 10.1016/S0034-4257(01)00281-4View ArticleGoogle Scholar
- Mascaro J, Detto M, Asner GP, Muller-Landau HC: Evaluating uncertainty in mapping forest carbon with airborne LiDAR. Remote Sensing of Environment 2011, 115: 3770–3774. 10.1016/j.rse.2011.07.019View ArticleGoogle Scholar
- Asner GP: Tropical forest carbon assessment: integrating satellite and airborne mapping approaches. Environmental Research Letters 2009, 3: 1748–9326.Google Scholar
- Helmer EH, Lefsky MA, Roberts DA: Biomass accumulation rates of Amazonian secondary forest and biomass of old-growth forests from Landsat time series and the Geoscience Laser Altimeter System. Journal of Applied Remote Sensing 2009, 3: 033505. 10.1117/1.3082116View ArticleGoogle Scholar
- 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. Proceedings of the National Academy of Sciences 2010, 107: 16738–16742. 10.1073/pnas.1004875107View ArticleGoogle Scholar
- Herold M, Skutsch M: Monitoring, reporting and verification for national REDD + programmes: two proposals. Environmental Research Letters 2011, 6: 014002. 10.1088/1748-9326/6/1/014002View ArticleGoogle Scholar
- Du Puy D, Moat J: A refined classification of the primary vegetation of Madagascar based on the underlying geology: Using GIS to map its distribution and to assess its conservation status. Edited by: Lourenco WR. Editions de l'ORSTOM - Paris; 1996:205–218. Proceedings of the International Symposium on the Biogeography of Madagascar; September 1995; ParisGoogle Scholar
- Justice CO, Vermote E, Townshend JRG, DeFries R, Roy DP, Hall DK, Salomonson VV, Privette JL, Riggs G, Strahler A, Lucht W, Myneni RB, Knyazikhin Y, Running SW, Nemani RR, Wan Z, Huete AR, Leeuwen Wv, Wolfe RE, Giglio L, Muller J-P, Lewis P, Barnsley MJ: The Moderate Resolution Imaging Spectroradiometer (MODIS): Land remote sensing for global change research. Geoscience and Remote Sensing 1998, 36: 1212–1227. 10.1109/36.701027View ArticleGoogle Scholar
- 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. Proceedings of the National Academy of Sciences 2011, 108: 9899–9904. 10.1073/pnas.1019576108View ArticleGoogle Scholar
- Ranson KJ, Sun G, Weishampel JF, Knox RG: Forest biomass from combined ecosystem and radar backscatter modeling. Remote Sensing of Environment 1997, 59: 118–133. 10.1016/S0034-4257(96)00114-9View ArticleGoogle Scholar
- Luckman A, Baker J, Kuplich TM, Yanasse CDF, Frery AC: A study of the relationship between radar backscatter and regenerating tropical forest biomass for spaceborne SAR instruments. Remote Sensing of Environment 1997, 60: 1–13. 10.1016/S0034-4257(96)00121-6View ArticleGoogle Scholar
- Kasischke ES, Melack JM, Dobson MC: The use of imaging radars for ecological applications - a review. Remote Sensing of Environment 1997, 59: 141–156. 10.1016/S0034-4257(96)00148-4View ArticleGoogle Scholar
- Clark DB, Clark DA: Landscape-scale variation in forest structure and biomass in a tropical rain forest. Forest Ecology and Management 2000, 137: 185–198. 10.1016/S0378-1127(99)00327-8View ArticleGoogle Scholar
- Saatchi SS, Houghton RA, Dos Santos Alvala RC, Soares JV, Yu Y: Distribution of aboveground live biomass in the Amazon basin. Global Change Biology 2007, 13: 816–837. 10.1111/j.1365-2486.2007.01323.xView ArticleGoogle Scholar
- Asner GP, Mascaro J, Muller-Landau HC, Vieilledent G, Vaudry R, Rasamoelina M, Hall J, van Breugal M: A universal airborne LiDAR approach for tropical forest carbon mapping. Oecologia 2011. doi:10.1007/s00442–011–2165-z doi:10.1007/s00442-011-2165-zGoogle Scholar
- Mascaro J, Asner GP, Muller-Landau HC, van Breugal M, Hall J, Dahlin K: Controls over aboveground forest carbon density on Barro Colorado Island, Panama. Biogeosciences 2011, 8: 1615–1629. 10.5194/bg-8-1615-2011View ArticleGoogle Scholar
- Drake JB, Knox RG, Dubayah RO, Clark DB, Condit R, Blair JB, Hofton M: Aboveground biomass estimation in closed canopy neotropical forests using lidar remote sensing: Factors affecting the generality of relationships. Global Ecology and Biogeography 2003, 12: 147–159. 10.1046/j.1466-822X.2003.00010.xView ArticleGoogle Scholar
- Keller M, Palace M, Hurtt G: Biomass estimation in the Tapajos National Forest, Brazil: examination of sampling and allometric uncertainties. Forest Ecology and Management 2001, 154: 371–382. 10.1016/S0378-1127(01)00509-6View ArticleGoogle Scholar
- Chave J, Chust G, Condit R, Aguilar S, Hernandez A, Lao S, Perez R: Error propagation and scaling for tropical forest biomass estimates. In Tropical forests and global atmospheric change. Edited by: Malhi Y, Phillips O. London: Oxford University Press; 2004:155–166.Google Scholar
- Mueller-Dombois D, Fosberg RF: Vegetation of the tropical Pacific islands. New York: Springer-verlag; 1998.View ArticleGoogle Scholar
- Allnutt TF, Ferrier S, Manion G, Powell GVN, Ricketts TH, Fisher BL, Harper GJ, Irwin ME, Kremen C, Labat J-N, Lees DC, Pearce TA, Rakotondrainibe F: A method for quantifying biodiversity loss and its application to a 50-year record of deforestation across Madagascar. Conservation Letters 2008, 1: 173–181. 10.1111/j.1755-263X.2008.00027.xView ArticleGoogle Scholar
- Gade DW: Deforestation and Its Effects in Highland Madagascar. Mountain Research and Development 1996, 16: 101–116. 10.2307/3674005View ArticleGoogle Scholar
- Baccini A, Laporte N, Goetz SJ, Sun M, Dong H: A first map of tropical Africa's above-ground biomass derived from satellite imagery. Environmental Research Letters 2008, 3: 045011. 10.1088/1748-9326/3/4/045011View ArticleGoogle Scholar
- Asner GP, Knapp DE, Balaji A, Paez-Acosta G: Automated mapping of tropical deforestation and forest degradation: CLASlite. Journal of Applied Remote Sensing 2009, 3: 033543. 10.1117/1.3223675View ArticleGoogle Scholar
- Asner GP, Knapp DE, Kennedy-Bowdoin T, Jones MO, Martin RE, Boardman J, Field CB: Carnegie Airborne Observatory: In-flight fusion of hyperspectral imaging and waveform light detection and ranging (LiDAR) for three-dimensional studies of ecosystems. Journal of Applied Remote Sensing 2007., 1: DOI:10-.1117/1111.2794018 DOI:10-.1117/1111.2794018Google Scholar
- Chave J, Andalo C, Brown S, Cairns MA, Chambers JQ, Eamus D, Fölster H, Fromard F, Higuchi N, 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-xView ArticleGoogle Scholar
- Vieilledent G, Vaudry R, Andriamanohisoa SFD, Rakotonarivo OS, Randrianasolo HZ, Razafindrabe HN, Bidaud Rakotoarivony C, Ebeling J, Rasamoelina M: A universal approach to estimating biomass and carbon stock in tropical forests using generic allometric models. Ecological Applications 2011. doi: . http://dx.doi.org/10.1890/11–0039-.1 doi:.Google Scholar
- Popescu SC, Wynne RH, Scrivani JA: Fusion of small-footprint lidar and multispectral data to estimate plot-level volume and biomass in deciduous and pine forests in Virginia, USA. Forest Science 2004, 50: 551–565.Google Scholar
- Drake JB, Dubayah RO, Knox RG, Clark DB, Blair JB: Sensitivity of large-footprint lidar to canopy structure and biomass in a neotropical rainforest. Remote Sensing of Environment 2002, 81: 378–392. 10.1016/S0034-4257(02)00013-5View ArticleGoogle Scholar
- Lefsky MA, Cohen WB, Acker SA, Parker GG, Spies TA, Harding D: Lidar remote sensing of the canopy structure and biophysical properties of Douglas-fir western hemlock forests. Remote Sensing of Environment 1999, 70: 339–361. 10.1016/S0034-4257(99)00052-8View ArticleGoogle Scholar
- Lefsky MA, Cohen WB, Parker GG, Harding DJ: Lidar remote sensing for ecosystem studies. BioScience 2002, 52: 19–30. 10.1641/0006-3568(2002)052[0019:LRSFES]2.0.CO;2View ArticleGoogle Scholar
- Asner GP, Hughes RF, Varga TA, Knapp DE, 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-5View ArticleGoogle Scholar
- Baskerville G: Use of logarithmic regression in the estimation of plant biomass. Canadian Journal of Forest Research-Revue Canadienne De Recherche Forestiere 1972, 2: 49–53. 10.1139/x72-009View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.