Mapping and monitoring of forest carbon stocks across large areas in the tropics will necessarily rely on remote sensing approaches, which in turn depend on field estimates of biomass for calibration and validation purposes. Here, we used field plot data collected in a tropical moist forest in the central Amazon to gain a better understanding of the uncertainty associated with plot-level biomass estimates obtained specifically for the calibration of remote sensing measurements. In addition to accounting for sources of error that would be normally expected in conventional biomass estimates (e.g., measurement and allometric errors), we examined two sources of uncertainty that are specific to the calibration process and should be taken into account in most remote sensing studies: the error resulting from spatial disagreement between field and remote sensing measurements (i.e., co-location error), and the error introduced when accounting for temporal differences in data acquisition. We found that the overall uncertainty in the field biomass was typically 25% for both secondary and primary forests, but ranged from 16 to 53%. Co-location and temporal errors accounted for a large fraction of the total variance (>65%) and were identified as important targets for reducing uncertainty in studies relating tropical forest biomass to remotely sensed data. Although measurement and allometric errors were relatively unimportant when considered alone, combined they accounted for roughly 30% of the total variance on average and should not be ignored. Our results suggest that a thorough understanding of the sources of error associated with field-measured plot-level biomass estimates in tropical forests is critical to determine confidence in remote sensing estimates of carbon stocks and fluxes, and to develop strategies for reducing the overall uncertainty of remote sensing approaches.

We investigated the potential of hyperspectral remote sensing to estimate aboveground biomass (AGB) over the Brazilian savannas (Cerrado), the second-largest source of carbon emissions in Brazil. For this purpose, a Hyperion/Earth Observing-1 (EO-1) image was collected in the dry season at the Ecological Station of Águas Emendadas (ESAE). In order to estimate the AGB, we evaluated the performance of five machine learning models (Classification and Regression Trees – CART; Cubist – CB, Partial Least Squares Regression – PLS; Random Forest – RF; and Support Vector Machine – SVM) and four sets of metrics (reflectance, narrowband vegetation indices – VIs; absorption band parameters; and the combination of these attributes). The lowest root mean square error (RMSE) was obtained for RF using VIs (29%) and a combination of metrics (28%). For VIs, RF differed from CUB, PLS and SVM at 5% significance level. From cross-validation results, the RMSE was 26.36% for grasslands, 35.04% for open savannas, and 24.85% for dense savannas. The RF model with VIs had the most stable predictive performance across the models, as indicated by small variations in RMSE from CART to SVM. The five most important ranked VIs in the RF model were the Normalized Difference Vegetation Index (NDVI), Pigment Specific Simple Ratio (PSSR), Enhanced Vegetation Index (EVI), Red Edge Normalized Difference Vegetation Index (RENDVI) and Structure Insensitive Pigment Index (SIPI). Most of their relationships with AGB were non-linear. The resultant AGB estimates showed consistent results with a vegetation cover map of the ESAE. Areas of the ESAE with AGB lower than 10 Mg.ha⁻¹ were coincident with the occurrence of grassland physiognomies (savanna grasslands and shrub savannas), while areas with AGB higher than 25 Mg.ha⁻¹ matched the occurrence of dense savanna physiognomies (woodland savanna and dense woodland savanna). Grassland areas showed larger values of coefficient of variation (CV) than areas of dense savannas. These first-hand results set a baseline of models and metrics for AGB modeling of savannas during the future transition from current sampling-type hyperspectral missions (< 10 km of swath) to large-coverage hyperspectral satellites (> 100 km of swath).

Changes in tropical-forest structure and aboveground biomass (AGB) contribute directly to atmospheric changes in CO${}_2$ , which, in turn, bear on global climate. This paper demonstrates the capability of radar-interferometric phase-height time series at X-band (wavelength = 3 cm) to monitor changes in vertical structure and AGB, with sub-hectare and monthly spatial and temporal resolution, respectively. The phase-height observation is described, with a focus on how it is related to vegetation-density, radar-power vertical profiles, and mean canopy heights, which are, in turn, related to AGB. The study site covers $18\times 60$ km in the Tapajós National Forest in the Brazilian Amazon. Phase-heights over Tapajós were measured by DLR’s TanDEM-X radar interferometer 32 times in a 3.2 year period from 2011–2014. Fieldwork was done on 78 secondary and primary forest plots. In the absence of disturbance, rates of change of phase-height for the 78 plots were estimated by fitting the phase-heights to time with a linear model. Phase-height time series for the disturbed plots were fit to the logistic function to track jumps in phase-height. The epochs of clearing for the disturbed plots were identified with ≈1-month accuracy. The size of the phase-height change due to disturbance was estimated with ≈2-m accuracy. The monthly time resolution will facilitate REDD+ monitoring. Phase-height rates of change were shown to correlate with LiDAR RH90 height rates taken over a subset of the TanDEM-X data’s time span (2012–2013). The average rate of change of phase-height across all 78 plots was 0.5 m-yr ${}^{-1}$ with a standard deviation of 0.6 m-yr ${}^{-1}$ . For 42 secondary forest plots, the average rate of change of phase-height was 0.8 m-yr ${}^{-1}$ with a standard deviation of 0.6 m-yr ${}^{-1}$ . For 36 primary forest plots, the average phase-height rate was 0.1 m-yr ${}^{-1}$ with a standard deviation of 0.5 m-yr ${}^{-1}$ . A method for converting phase-height rates to AGB-rates of change was developed using previously measured phase-heights and field-estimated AGB. For all 78 plots, the average AGB-rate was 1.7 Mg-ha ${}^{-1}$ -yr ${}^{-1}$ with a standard deviation of 4.0 Mg-ha ${}^{-1}$ -yr ${}^{-1}$ . The secondary-plot average AGB-rate was 2.1 Mg-ha ${}^{-1}$ -yr ${}^{-1}$ , with a standard deviation of 2.4 Mg-ha ${}^{-1}$ -yr ${}^{-1}$ . For primary plots, the AGB average rate was 1.1 Mg-ha ${}^{-1}$ -yr ${}^{-1}$ with a standard deviation of 5.2 Mg-ha ${}^{-1}$ -yr ${}^{-1}$ . Given the standard deviations and the number of plots in each category, rates in secondary forests and all forests were significantly different from zero; rates in primary forests were consistent with zero. AGB-rates were compared to change models for Tapajós and to LiDAR-based change measurements in other tropical forests. Strategies for improving AGB dynamical monitoring with X-band interferometry are discussed.

In northern Brazilian Amazon, the crops, savannahs and rainforests form a complex landscape where land use and land cover (LULC) mapping is difficult. Here, data from the Operational Land Imager (OLI)/Landsat-8 and Phased Array type L-band Synthetic Aperture Radar (PALSAR-2)/ALOS-2 were combined for mapping 17 LULC classes using Random Forest (RF) during the dry season. The potential thematic accuracy of each dataset was assessed and compared with results of the hybrid classification from both datasets. The results showed that the combination of PALSAR-2 HH/HV amplitudes with the reflectance of the six OLI bands produced an overall accuracy of 83% and a Kappa of 0.81, which represented an improvement of 6% in relation to the RF classification derived solely from OLI data. The RF models using OLI multispectral metrics performed better than RF models using PALSAR-2 L-band dual polarization attributes. However, the major contribution of PALSAR-2 in the savannahs was to discriminate low biomass classes such as savannah grassland and wooded savannah.

We discriminated different successional forest stages, forest degradation, and land use classes in the Tapajós National Forest (TNF), located in the Central Brazilian Amazon. We used full polarimetric images from ALOS/PALSAR-2 that have not yet been tested for land use and land cover (LULC) classification, neither for forest degradation classification in the TNF. Our specific objectives were: (1) to test the potential of ALOS/PALSAR-2 full polarimetric images to discriminate LULC classes and forest degradation; (2) to determine the optimum subset of attributes to be used in LULC classification and forest degradation studies; and (3) to evaluate the performance of Random Forest (RF) and Support Vector Machine (SVM) supervised classifications to discriminate LULC classes and forest degradation. PALSAR-2 images from 2015 and 2016 were processed to generate Radar Vegetation Index, Canopy Structure Index, Volume Scattering Index, Biomass Index, and Cloude–Pottier, van Zyl, Freeman–Durden, and Yamaguchi polarimetric decompositions. To determine the optimum subset, we used principal component analysis in order to select the best attributes to discriminate the LULC classes and forest degradation, which were classified by RF. Based on the variable importance score, we selected the four first attributes for 2015, alpha, anisotropy, volumetric scattering, and double-bounce, and for 2016, entropy, anisotropy, surface scattering, and biomass index, subsequently classified by SVM. Individual backscattering indexes and polarimetric decompositions were also considered in both RF and SVM classifiers. Yamaguchi decomposition performed by RF presented the best results, with an overall accuracy (OA) of 76.9% and 83.3%, and Kappa index of 0.70 and 0.80 for 2015 and 2016, respectively. The optimum subset classified by RF showed an OA of 75.4% and 79.9%, and Kappa index of 0.68 and 0.76 for 2015 and 2016, respectively. RF exhibited superior performance in relation to SVM in both years. Polarimetric attributes exhibited an adequate capability to discriminate forest degradation and classes of different ecological succession from the ones with less vegetation cover.

This study aims to analyze the capability of the target decomposition techniques and the polarimetric ratios applied to the ALOS/PALSAR-2 satellite polarimetric images to discriminate the land use and land cover classes in the Tapajós National Forest region, Pará State. Three full polarimetric ALOS/PALSAR-2, level 1 single look complex scenes were selected to generate the coherence and the covariance matrices to derive the Cloude-Pottier and the Freeman-Durden target decomposition attributes. From the radiometrically calibrated PALSAR-2 images, we generated the backscatter coefficients, the cross polarized ratio (RC; HV/HH), the parallel polarized ratio (RP; VV/HH) and the Radar Forest Degradation Index (RFDI). The images resulting from these polarimetric attributes were processed by the Maximum Likelihood (MAXVER) classifier coupled with the Iterated Conditional Modes (ICM) contextual algorithm. We found that the classifications derived from the target decomposition attributes, mainly from the Cloude-Pottier technique, with a Kappa index of 0.75, presented a significant higher performance than those derived from the RC ratio, RP ratio, and RFDI.