SAR-Based Estimation of Above-Ground Biomass and Its Changes in Tropical Forests of Kalimantan Using Land C-Band

Kalimantan poses one of the highest carbon emissions worldwide since its landscape is strongly endangered by deforestation and degradation and, thus, carbon release. The goal of this study is to conduct large-scale monitoring of above-ground biomass (AGB) from space and create more accurate biomass maps of Kalimantan than currently available. AGB was estimated for 2007, 2009, and 2016 in order to give an overview of ongoing forest loss and to estimate changes between the three time steps in a more precise manner. Extensive field inventory and LiDAR data were used as reference AGB. A multivariate linear regression model (MLR) based on backscatter values, ratios, and Haralick textures derived from Sentinel-1 (C-band), ALOS PALSAR (Advanced Land Observing Satellite’s Phased Array-type L-band Synthetic Aperture Radar), and ALOS-2 PALSAR-2 polarizations was used to estimate AGB across the country. The selection of the most suitable model parameters was accomplished considering VIF (variable inflation factor), p-value, R2, and RMSE (root mean square error). The final AGB maps were validated by calculating bias, RMSE, R2, and NSE (Nash-Sutcliffe efficiency). The results show a correlation (R2) between the reference biomass and the estimated biomass varying from 0.69 in 2016 to 0.77 in 2007, and a model performance (NSE) in a range of 0.70 in 2016 to 0.76 in 2007. Modelling three different years with a consistent method allows a more accurate estimation of the change than using available biomass maps based on different models. All final biomass products have a resolution of 100 m, which is much finer than other existing maps of this region (>500 m). These high-resolution maps enable identification of even small-scaled biomass variability and changes and can be used for more precise carbon modelling, as well as forest monitoring or risk managing systems under REDD+ (Reducing Emissions from Deforestation, forest Degradation, and the role of conservation, sustainable management of forests, and enhancement of forest carbon stocks) and other programs, protecting forests and analyzing carbon release.


Introduction
The Earth's land surface spans approximately 149.4 million km 2 , of which nearly 30% is characterized by forested areas [1].Tropical forest ecosystems (forests and soils) alone hold about 40% of terrestrial carbon [2,3].However, due to unsustainable use and deforestation, the stored carbon can be released into the atmosphere as CO 2 (carbon dioxide) and will contribute significantly to global climate change.According to the Intergovernmental Panel on Climate Change (IPCC), the resulting greenhouse gas emissions account for approximately 11% of all anthropogenic emissions worldwide [4,5].
The forests of Indonesia are considered to be one of the oldest and most species-rich tropical rainforests on the planet [6].Indonesia's forests alone store 18.6 Gt of carbon [7].Additionally, the country has some of the largest known tropical peat reservoirs on Earth, storing 55-58 Gt of carbon in belowground peatlands ( [8,9]).At the same time the deforestation and degradation of these tropical ecosystems lead to a considerable release of carbon, making Indonesia, especially Kalimantan and Sumatra, one of the largest carbon emitters worldwide [10].Digging canals for means of transport is a detrimental form of degradation, as it causes drying in peatlands and encourages the spread of fires and, thus, increases the potential of forest and peatland loss.As a result of the significant emission levels and immense loss of forests and peatlands, Indonesia has become a prime target for REDD+ projects (Reducing Emissions from Deforestation, forest Degradation, and the role of conservation, sustainable management of forests, and enhancement of forest carbon stocks).One of the intentions of REDD+ includes conditional payments to developing countries for reducing their emissions [11].Within this framework, the most active projects and initiatives worldwide are taking place in Indonesia and its provinces [12].Nevertheless, Indonesia, especially Kalimantan, the largest isle of Indonesia, is affected by ongoing deforestation and degradation processes.
To perform current REDD+ policies, accurate forest monitoring systems, consistent measurements, and information about carbon emissions at national and subnational scales are necessary for participating countries.Especially, subnational projects require accurate high-resolution maps capturing small-scaled variability and changes in forested areas.Forest carbon stocks are usually calculated using above-ground biomass (AGB) by assuming that typically 50% of AGB is carbon [13].Biomass, as the fundamental biophysical parameter quantifying the Earth's living vegetation [14], describes the amount of woody matter within a forest.It is defined by the Global Climate Observing System (GCOS) as an essential climate variable (ECV) [15].Well-known methods for mapping AGB include field-data, airborne and space-borne LiDAR (light detection and ranging) scanning, satellite optical remote sensing, and imaging radar [14].In contrast to ground-based inventories and LiDAR surveys, Earth observation approaches are able to cover larger areas and in a more cost-effective manner.Additionally, forest inventories are not always comparable, since the definition of national forests and sampling strategies vary between countries [16].
An appropriate compromise represents the upscaling of accurate forest inventories or regional LiDAR-derived biomass estimations with large-scale satellite imagery [17].Since optical data are limited by clouds, smoke, and lacking illumination, as well as not being able to capture the vertical structures of trees [18], synthetic aperture radar (SAR) is more suitable in this scope of application, as it is daylight-and almost weather-independent [19].The SAR data-based backscatter approach is well known for forest cover and biomass mapping.This method uses the energy that is received by the sensor after transmission, the so-called backscatter, and relates it to field biomass measurements.The backscatter typically increases with an increasing amount of biomass until a certain value, at which the sensitivity of the backscatter to the AGB stagnates.This biomass saturation level is dependent on the wavelength of the sensor [20].With regard to the sensitivity of vegetation, wavelengths underlay different physical characteristics.While the C-band is able to penetrate through leaves, but is scattered by small branches, the L-band, with a wavelength of up to 30 cm, is scattered mainly by trunks and tall branches.Since P-band SAR data, which is able to penetrate deeply into the canopy cover and is backscattered by trunks and the ground, has been unavailable to date, the L-band represents the most suitable operational data for biomass estimation [21].AGB estimation based on Advanced Land Observing Satellite's Phased Array-type L-band Synthetic Aperture Radar (ALOS PALSAR) data has already been successfully performed by [22,23].AGB studies in tropical forests were also mostly conducted on the basis of L-band SAR data [24][25][26][27][28]. Reported saturation levels using L-band in tropical forests in Indonesia are at approximately 50 t/ha to 200 t/ha [25,[29][30][31][32].A combination of C-and L-bands has shown better results in tropical forests of Colombia than using a single band [33,34] found a correlation of C-band ENVISAT ASAR backscatter and AGB up to 250 t/ha.
Existing biomass maps of pan-tropical ecosystems by [7,30] exhibit resolutions of 500 m or 1 km, respectively.In Southeast Asia the biomass map of [7] overestimates AGB in the lower biomass ranges and does not identify heterogeneity in forests in detail.In contrast, the map of [30] slightly underestimates high biomass ranges, but captures disturbances in forested areas [35].Avitabile et al. [36] fuse those two maps in combination with additional data to create an improved pan-tropical biomass map.The results show a smaller RMSE and bias for all continents.Nevertheless, the spatial resolution is about 1 km, which is why the AGB heterogeneity and small-scaled changes cannot be detected as accurate as necessary for most REDD+ projects.Freely available, large-scale biomass maps at national or subnational scales for Indonesia are not available to date.
The aim of the ESA DUE GlobBiomass project is to improve existing AGB estimation products and reduce uncertainties in different ecosystems by developing an innovative synergistic mapping approach.Within the project above-ground biomass is estimated for five regional sites (Sweden, Poland, Mexico, Eastern South Africa, Kalimantan) for the epochs 2005 (2005 ± 2 years), 2010 (2010 ± 2 years), and 2015 (2015 ± 2 years), as well as a global map for the epoch 2010, using different data and methods.Additionally, the change of biomass during the three epochs is estimated per region.As part of the GlobBiomass project, the present study focuses on large-scale monitoring of biomass in Kalimantan, Indonesia, from space.The aim of the study is to estimate AGB with a finer resolution and better accuracy than other existing AGB maps.Therefore, ALOS PALSAR/ALOS-2 PALSAR-2 L-band, and Sentinel-1 C-band data from 2007, 2009, and 2016 were used.A multivariate linear regression model (MLR) based on SAR backscatter values, polarization ratios, and textures was set up in order to increase the biomass saturation level.As the reference biomass for the model calibration and validation from SAR data, a combination of field inventory data and LiDAR data was used in order to provide a more accurate base.Modelling three different years with a consistent method allows a more accurate estimation of changes between the three time steps than using the available biomass maps based on different models.A higher spatial resolution is important in order to make them a promising alternative building a forest monitoring or risk managing system, but also to achieve the objectives of REDD+, the Global Canopy Programme, UNEP-WCMC, and other programs protecting forests or analyzing carbon release at national and subnational levels.

Study Area and Data
Kalimantan, which is the Indonesian part of Borneo island, has a size of about 544,000 km 2 and lies within the geographic coordinates 4 • 15 41 N to 3 • 45 44 S latitude and 108 • 48 0 E to 118 • 49 41 E longitude (Figure 1).The island's climate is mainly conditioned by the dry southeast monsoon from May to October and the wet northwest monsoon from November to April, and is influenced by frequent rainfall and high temperatures throughout the year.Those conditions are ideal for plant growth, which is why Kalimantan's land cover is characterized mainly by tropical forests covering 301,750 km 2 and, thus, more than 55% of the country.The forests of Borneo are considered to be one of the oldest and most species-rich tropical rainforests on Earth.The dominating forest ecosystems are mangrove forests, peat swamp and freshwater swamp forests, riparian forests, heath forests, lowland dipterocarp forests, ironwood forests, forests on limestone and ultrabasic soils, hill dipterocarp forests, and various montane formations [6].In general lowland dipterocarp and peat swamp forests can be well discriminated in the field by means of average tree height, tree crown diameter, canopy closure, and species composition.Lowland dipterocarp forests are more diverse with taller trees and a more closed canopy [6].Dipterocarps can reach a height of 45-60 m and are a valuable tree species prone to logging.All of the forest types store an extensive amount of carbon [8].Nevertheless, the most significant carbon sinks in this area are belowground peatlands, which can store up to ten times more carbon than the forests growing on top them, since they were formed over the past millennia, as plant debris accumulated under waterlogged conditions [37].The last decades witnessed a decrease of forest due to illegal logging, deforestation for agricultural development but also due to natural and manmade fires.Additionally, the ability of peatlands to store carbon is reduced due to anthropogenic activities like logging and drainage.In particular, draining the normally waterlogged peatlands makes these ecosystems vulnerable to fires.For modelling the reference AGB layer a universal AGB model is used for all different forest types.normally waterlogged peatlands makes these ecosystems vulnerable to fires.For modelling the reference AGB layer a universal AGB model is used for all different forest types.

Above-Ground Data
For the creation of the AGB reference layer existing field inventory data, collected across different forest types using the nested plots method, based on the guidelines provided by [38], were utilized.The methodology varies slightly for the various study sites, since they were processed as part of different research projects.The data was gathered in four test sites across Kalimantan (Figure 1) during 2007, 2009, and 2016 (Table 1).In order to estimate AGB using allometric models after [39], information about forest type, tree species, the diameter at breast height (DBH), and tree height are collected within nested plots based on three or four subplots.Inside the different subplots, trees of a certain diameter at breast height (DBH) were measured, for example, DBH < 10 cm (within the 3 m radius), ≥10 cm to <20 cm (10 m radius), ≥20 cm to <50 cm (20 m radius), and ≥50 cm (30 m radius).The sum of the measured parameters of the subplots was multiplied by an expansion factor in order to obtain the final values for 1 ha.In Pulang

Above-Ground Data
For the creation of the AGB reference layer existing field inventory data, collected across different forest types using the nested plots method, based on the guidelines provided by [38], were utilized.The methodology varies slightly for the various study sites, since they were processed as part of different research projects.The data was gathered in four test sites across Kalimantan (Figure 1) during 2007, 2009, and 2016 (Table 1).In order to estimate AGB using allometric models after [39], information about forest type, tree species, the diameter at breast height (DBH), and tree height are collected within nested plots based on three or four subplots.Inside the different subplots, trees of a certain diameter at breast height (DBH) were measured, for example, DBH < 10 cm (within the 3 m radius), ≥10 cm to <20 cm (10 m radius), ≥20 cm to <50 cm (20 m radius), and ≥50 cm (30 m radius).The sum of the measured parameters of the subplots was multiplied by an expansion factor in order to obtain the final values for 1 ha.In Pulang Airborne LiDAR measurements were acquired during the dry seasons in 2007, 2011, and 2012 in the same areas as the field data was collected.In 2007 (May-October) a Riegl LMS-Q560 2D laser scanner by RIEGL Laser Measurement Systems GmbH (Horn, Austria) was flown at a height of 500 m above-ground and a half scan angle of ±30 • to collect the full-wave LiDAR data.The average point density of the final data 2007 was 1.5 points per m 2 .For the years 2011 and 2012 (August-October) the measurements were acquired using Optech Orion M200 and Optech ALTM 3100 airborne laser scanners by Teledyne Optech (Vaughan, Ontario, Canada) at an altitude of 800 m above-ground.A half scan angle of ±11 • was used and the point density amounted to 10.7 points per m 2 .Since the accuracy of biomass estimations derived from LiDAR metrics increase with a higher point density, a weighting of the plots accordingly to their point density was applied [40].In total, 8300 km 2 were surveyed in different regions in East, West, and Central Kalimantan (Table 1) representing different ecosystems.Since there was a lack of LiDAR data for 2016 due to vast fires in Kalimantan, adjusted data of 2011 and 2012 was used.

SRTM
The Digital Elevation Model (DEM) from the Shuttle Radar Topography mission (SRTM) with a vertical accuracy of ±10 m and a spatial resolution of 30 m is used for topographic analyses.Slope is used to clean up the final AGB map, since extreme overestimation of biomass occurs in steep terrain.Besides, steep terrain can cause layover and shadow effects, decreasing the accuracy of the AGB estimation.

TRMM
Since SAR backscatter is highly sensitive to water content of the surface due to its dielectric properties, daily TRMM (Tropical Rainfall Measuring Mission) of the National Aeronautics and Space Administration (NASA) and JAXA precipitation data with a spatial resolution of 0.25 • were incorporated in order to select satellite acquisition dates with dry and comparable conditions [41].

MODIS Active Fire Data (MDC14DL)
MODIS hotspot information (product MCD14DL, provided by NASA) were used as an additional layer to detect thermal anomalies and active fires in Kalimantan.In this layer active fires are represented in the center of a 1 km × 1 km pixel that is identified by the MODIS MOD14/MYD14 Fire and Thermal Anomalies algorithm as containing one or more fires within the pixel [42,43].

Water Body Mask
ESRI World Water Bodies was used for delineating water bodies within the study area.It provides a base map layer for lakes, seas, oceans and large rivers and as generated on data with a spatial resolution of 100 m [44].

Urban Areas
The ESA CCI land cover map provides three epoch series (2000,2005,2010) of global land cover maps at 300 m spatial resolution which were used to evaluate settlement areas.ESA CCI land cover maps were produced used a multi-sensor and multi-temporal strategy based on MERIS Full and Reduced Resolution (FR and RR) archive.The 10-year product has been served as a baseline to derive the 2000, 2005 and 2010 maps using MERIS and SPOT-Vegetation time series specific to each epoch [45].

Above-Ground Biomass Data
Inventory data and LiDAR data are combined to create more accurate biomass predictions for an area within the SAR images.This upscaling from field inventory to LiDAR transects allows creating a more precise basis for AGB model calibration and validation from SAR backscatter data [40].An overview of the whole methodology of this study can be found in Figure 2.
Space Administration (NASA) and JAXA precipitation data with a spatial resolution of 0.25° were incorporated in order to select satellite acquisition dates with dry and comparable conditions [41].

MODIS Active Fire Data (MDC14DL)
MODIS hotspot information (product MCD14DL, provided by NASA) were used as an additional layer to detect thermal anomalies and active fires in Kalimantan.In this layer active fires are represented in the center of a 1 km × 1 km pixel that is identified by the MODIS MOD14/MYD14 Fire and Thermal Anomalies algorithm as containing one or more fires within the pixel [42,43].

Water Body Mask
ESRI World Water Bodies was used for delineating water bodies within the study area.It provides a base map layer for lakes, seas, oceans and large rivers and as generated on data with a spatial resolution of 100 m [44].

Urban Areas
The ESA CCI land cover map provides three epoch series (2000,2005,2010) of global land cover maps at 300 m spatial resolution which were used to evaluate settlement areas.ESA CCI land cover maps were produced used a multi-sensor and multi-temporal strategy based on MERIS Full and Reduced Resolution (FR and RR) archive.The 10-year product has been served as a baseline to derive the 2000, 2005 and 2010 maps using MERIS and SPOT-Vegetation time series specific to each epoch [45].

Above-Ground Biomass Data
Inventory data and LiDAR data are combined to create more accurate biomass predictions for an area within the SAR images.This upscaling from field inventory to LiDAR transects allows creating a more precise basis for AGB model calibration and validation from SAR backscatter data [40].An overview of the whole methodology of this study can be found in Figure 2. The estimation of AGB in t/ha from field inventory data was achieved by using the tree height, diameter at breast height (DBH), and wood specific density of each tree as the input for a combination of different allometric models.Therefore, allometric models from [46] for saplings (if DBH < 5 cm and height ≤ 1.3 m) or trees (if DBH < 5 cm and height > 1.3 m) and from [39] for moist tropical forest stands (if DBH ≥ 5 cm and height > 1.3 m) were applied.Those ground-based AGB values were related to LiDAR transects in order to estimate biomass reference data using previously-established regression models for the different study sites, as described in [40].
LiDAR height histograms were calculated by normalizing all points within a grid of, e.g., 30 m (similar to the size of the largest nest of the field inventory plots) to the ground using a DTM as reference, i.e., the height of each LiDAR return was calculated relative to the DTM [47,48].The number of points within each 0.5 m interval was stored as a histogram.The first (lowest) interval was considered as the ground return and excluded from further processing.The LiDAR data enabled calculating the two LiDAR point cloud metrics quadratic mean canopy height (QMCH) and centroid height (CH), as described more in detail in [40].Previous studies have shown that those two parameters are able to estimate tropical forests AGB, as they take into account the point distribution of the LiDAR measurements [47,49].This is why the QMCH and the CH were estimated based on height histograms by weighting each 0.5 m height interval by the fraction of points within the interval.CH and QMCH were related to AGB, estimated by the field inventories using regression models for each particular study site in order to obtain large-scale biomass estimations [29,48].Furthermore, the point density was included in the regression model, since [47] showed that the accuracy of AGB estimations derived from LiDAR height histograms increased with higher point densities.The coefficient of determination (R 2 ) of the AGB regression models vary from R 2 = 0.7 (Berau and Kapuas Hulu) to R 2 = 0.8 (Malinau and Central Kalimantan).
The upscaling from field inventory data to LiDAR transects allows creating numerous biomass reference data for the calibration of SAR images and to upscale AGB across large areas and different ecosystems.It provides AGB estimates over the whole biomass range from woody regrowth to pristine forest, and is able to disclose a spatial variation due to varying growth conditions.In order to estimate AGB based on SAR data the reference AGB was rescaled to a spatial resolution of 100 m.
Since there was a lack of LiDAR data for 2016 due to vast fires in Kalimantan, the data of 2011 and 2012 was used.Differences between the reference AGB and the SAR images over time were detected and excluded using MODIS fire hotspots, as fire is the main reason for deforestation in this area.

SAR Data
Since SAR backscatter is highly sensitive to water content of the surface due to its dielectric properties, daily TRMM (Tropical Rainfall Measuring Mission) precipitation data with a spatial resolution of 0.25 • were used to verify the moisture conditions within the SAR imageries.In order to have dry and comparable conditions images with a high influence of precipitation were excluded.During the pre-processing of the Sentinel-1 (10 m) and ALOS PALSAR (25 m) data a co-registration based on an ALOS PALSAR mosaic with a spatial resolution of 25 m, a radiometric calibration estimating γ 0 backscatter coefficients in dB and a geometric correction were accomplished.Furthermore, a multi-temporal speckle filtering using an enhanced Lee filter with a 7 × 7 window was applied [50].The processed data was resampled to a resolution of 100 m resulting in pixels with a spatial resolution of 1 ha.
In order to examine the potential for AGB estimation, ratio images were prepared, using the following equations: R hvhh = HV/HH (1) where HH, HV, VV, and VH indicate the polarization of the γ 0 backscattering coefficients, depending on available polarizations of ALOS PALSAR and Sentinel-1.An evaluation of ratios and textures can be found in [51].
To improve the model accuracy, Haralick textures [52] and their relationship to AGB were also investigated.Texture describes the properties of objects, such as regularity, smoothness, and tonal variation [51].Ten simple Haralick textures applying a gray level co-occurrence matrix (GLCM) and eleven higher-order Haralick textures applying a gray level run-length matrix (GLRM) were calculated within the open-source software Orfeo ToolBox by CNES (Centre National D'Études Spatiales) over a moving window with user-defined radius based on single polarized images and the calculated ratios.The radius used for single polarized images had a size of five, while a radius of three was used for calculating textures based on ratios.The relationship between these generated SAR variables and the LiDAR AGB was analyzed and showed an inversely proportional correlation.The variables were linearized in order to use the variables as an input for a multivariate linear regression model (MLR).

Multivariate Linear Regression Model (MLR)
In Englhart et al. [53] it was found that artificial neural network (ANN) and support vector regression (SVR) were suitable to predict AGB in tropical forests.Nevertheless, in terms of biomass variability and saturation in tropical forest ecosystems, multivariate linear regression models (MLR) are superior to ANN and SVR models [53].Since very high biomass values are expected in the study area, MLR is used to model the above-ground biomass.Multiple regressions are often affected by overfitting and co-linearity among variables.In order to find a model with the greatest explanatory power, a backward stepwise multiple linear regression was performed to automate the selection of the best explanatory variables.The MLR was first set up with all linearized ratios and textures and run iteratively.To decrease the number of inputs, the p-value and the variable inflation factor (VIF) were investigated for each variable to identify their significance and co-linearity.Regarding literature, parameters with a p-value <0.05 and a VIF >5 were excluded from the model [26,51,54].After eliminating redundant information, three variables were finally used for 2007 and 2009 and four variables for 2016.However, the extreme fire events in 2015 burned vast forest areas in Kalimantan.Resulting burned and carbonized trunks enhance double bounce effects in this region, which caused higher backscatter values and, thus, an overestimation of the biomass model.To minimize the overestimation in this areas, a second model was set up for 2016 simply based on two variables that were less sensitive to high backscatter.This model was applied only in burned areas, captured by a mask based on high backscatter values.In the remaining parts of the scene the first model, using four inputs, was used.
Settlements and areas with steep terrains cannot be captured correctly by the model.In order to reduce errors due to radar shadow and layover effects, regions with a terrain steeper than 10 • were excluded from the biomass estimation using the Shuttle Radar Topography Mission (SRTM) DEM.Mountainous areas in Kalimantan are less influenced by human activity since they are difficult to access and due to a lack of transport routes, like canals.For this reason, the forests in mountainous terrain are mostly unaffected by degradation.Altogether 27.2% of the study area of Kalimantan are influenced by slopes steeper than 10 • the main forest type in these regions is the hill-and sub-montane forest which was determined during field inventory campaigns, as well as through detailed forest classification mapping covering the LiDAR areas.The mean value of LiDAR AGB was calculated for the area with a terrain steeper than 10 • and resulted in 350 t/ha.Additionally, urban areas and water bodies were excluded from the final results, using an ESRI World Water Body layer and the ESA CCI land cover map of each particular year (2005, 2010, 2015).

Temporal Change Estimation
The temporal change between 2007-2009, 2009-2016, and 2007-2016 was assessed using the RMSE in order to define a possible biomass range for each estimate at pixel level.The RMSE is calculated for five classes 0-50 t/ha, 50-100 t/ha, 100-150 t/ha, 150-200 t/ha, and >200 t/ha.The specific RMSE of the corresponding AGB range was subtracted from and added to the estimated AGB value at pixel level.The results are two layers consisting of the highest (H) and lowest value (L) of a possible AGB range per pixel for two comparing time-steps (T1, T2).A change/no-change mask, with discrimination of the increase and decrease, was generated by using threshold values based on the following decisions.If HT1 > LT2 or LT1 < HT2, an overlap between the ranges exist and, thus, no-change is assumed.In cases where HT1 < LT2, we suppose an increase of AGB and if LT1 > HT2 a decrease of AGB between the two years, both indicates a change.In order to remove single separated pixel, a minimum change unit of 300 m × 300 m was applied.

Validation and Uncertainty
For the validation of the AGB estimation, field inventory and LiDAR data, was used.The data was randomly split into data used for training (70%) and validation (30%) of the SAR-AGB models.For each map approximately 500 points were randomly set within the reference layer extent.To measure the accuracy of the models different parameters like bias, RMSE, standard deviation (SD) and the R 2 were calculated for each map within an AGB range of 50 t/ha starting from 0 up to >200 t/ha and the overall range.Additionally, the NSE was estimated for the overall biomass range.Regarding [55], the NSE is a dimensionless index measuring the efficiency of a model in a range from −∞ to 1.The closer the NSE is to 1, the more accurate is the model.The NSE is calculated using the following equation: where N is the number of observations, P i is the predicted value, O i is the observed value, and O is the mean of the observed value [55].
An uncertainty map at the pixel level is important for the interpretation of AGB maps, as shown in [56,57].The total uncertainty at the pixel level is composed of different sources of errors which are assumed to be random and independent.These are propagated for each map using the equation proposed by [30], taking into account the errors of measurement, allometry, sampling size, and prediction.Similar to [26], a measurement error supposed to be 10%.In [39] the authors found an error for the estimation of a tree's biomass of approximately ±5%.As we mainly used this allometry, an error of 5% is assumed.According to [30,58], a sampling size error of 20% is supposed.The prediction error includes the sampling error associated with the representativeness of the training data of the actual spatial distribution of AGB and the model predictions and is estimated per pixel.

Modelling Results
Using a backward stepwise approach allows to reduce the parameters within in the MLR, resulting in the finally-used variables listed in Table 2.In addition, coefficient of determination and residual standard error are displayed for each model calibration.Model 1 and model 2 of 2016 are used in combination, whereas model 2 is only applied in burned areas.The overall R 2 of the combined models is 0.69.The RP and the Low Grey-Level Run Emphasis (LGRE) showed the best relationship to LiDAR AGB compared to all Haralick textures.Additionally, cross-polarized-based parameters perform significantly better than co-polarized ones.

Biomass Maps and Temporal Change
The final biomass maps with a resolution of 100 m are presented in Figure 3a-c.Kalimantan is dominated in all years by forests with biomass varying in a range of 50-350 t/ha.However, its land cover is changing in time.Areas close to the coasts and along rivers show low biomass values between 0-50 t/ha.In 2007 15% of Kalimantan was covered by non-forested areas.Certainly, the portion of this class is growing over the years, caused by a decrease of forest and, thus, a loss of biomass.In 2009 the amount of non-forested extents is about 20%, though it reaches the maximum (25%) in 2016.The highest biomass values are reached in areas of mountainous terrain, which can be found in the north and center of Kalimantan.Nevertheless, the extent of forests containing high biomass values is significantly shrinking.In 2007 the model found a percentage of 55% of Kalimantan with biomass values >200 t/ha, while it is 38% in 2016.The biomass variability due to different degradation stages in the forest, as well as different disturbances in contrast to non-disturbed areas or clear-cuts, can be captured in the final maps.
A subset of Kalimantan, showing the change in the south of the island for the period 2009-2016, is displayed in Figure 3d.Red colors show a decrease of biomass, while green colors indicate an increase of AGB.The region is mainly dominated by a loss of biomass with values about − 300 to −100 t/ha.In contrast, only few areas show increased biomass.A similar distribution can be observed across the entire coastline of Kalimantan.Mountainous regions in the center of Kalimantan are less influenced by change.Figure 4 is summarizing the percentage of the forest degradation level per year.Highly-degraded areas increase from 15% in 2007 to 20% in 2009 to 25% in 2016.Accordingly, the area covered by natural forest (AGB > 200 t/ha) decreases.The coverage of forested areas containing an AGB from 50 to 200 t/ha are rising since areas with natural forest are affected by illegal logging, where single trees are felled.Thus the class of natural forest is slowly converted into degraded forest.Ongoing activities are further converting the degraded forest to highly-degraded areas.

Comparison with Pan-Tropical Biomass Maps
Figure 5 displays a subset of a degraded peat swamp forest in Central Kalimantan near the capital Palangka Raya.A visual comparison of the estimated AGB map, the LiDAR-derived reference AGB map and other pan-tropical biomass maps [7,30,36] with a resolution of 500 m and 1 km, respectively, points out that the developed MLR model correctly estimates the variability of biomass.The biomass map of [7] overestimates AGB in the lower biomass ranges and is not identifying heterogeneity in forests in detail.In contrast, the map of [30] underestimates high biomass ranges, but captures disturbances in forested areas.Avitabile et al. [36] combine the pantropical biomass

Comparison with Pan-Tropical Biomass Maps
Figure 5 displays a subset of a degraded peat swamp forest in Central Kalimantan near the capital Palangka Raya.A visual comparison of the estimated AGB map, the LiDAR-derived reference AGB map and other pan-tropical biomass maps [7,30,36] with a resolution of 500 m and 1 km, respectively, points out that the developed MLR model correctly estimates the variability of biomass.The biomass map of [7] overestimates AGB in the lower biomass ranges and is not identifying heterogeneity in forests in detail.In contrast, the map of [30] underestimates high biomass ranges, but captures disturbances in forested areas.Avitabile et al. [36] combine the pantropical biomass

Comparison with Pan-Tropical Biomass Maps
Figure 5 displays a subset of a degraded peat swamp forest in Central Kalimantan near the capital Palangka Raya.A visual comparison of the estimated AGB map, the LiDAR-derived reference AGB map and other pan-tropical biomass maps [7,30,36] with a resolution of 500 m and 1 km, respectively, points out that the developed MLR model correctly estimates the variability of biomass.The biomass map of [7] overestimates AGB in the lower biomass ranges and is not identifying heterogeneity in forests in detail.In contrast, the map of [30] underestimates high biomass ranges, but captures disturbances in forested areas.Avitabile et al. [36] combine the pantropical biomass maps of [7] and [30] and additional data in order to create an improved pan-tropical biomass map.The final product of [36] with a resolution of 1 km has a lower RMSE and bias than the two previous studies.Comparing maps (a) and (e), a similar trend of the biomass distribution is visible.The estimated AGB map based on the MLR model has a much finer resolution (100 m) which allows visualizing AGB variability more precisely and detects even small-scaled changes for carbon modelling, as well as forest monitoring or risk managing systems.Moreover, the modelled maps show a better accordance to the LiDAR-derived AGB in the subset.
Remote Sens. 2018, 10, x FOR PEER REVIEW of 22 maps of [7] and [30] and additional data in order to create an improved pan-tropical biomass map.
The final product of [36] with a resolution of 1 km has a lower RMSE and bias than the two previous studies.Comparing maps (a) and (e), a similar trend of the biomass distribution is visible.The estimated AGB map based on the MLR model has a much finer resolution (100 m) which allows visualizing AGB variability more precisely and detects even small-scaled changes for carbon modelling, as well as forest monitoring or risk managing systems.Moreover, the modelled maps show a better accordance to the LiDAR-derived AGB in the subset.

Validation and Uncertainty
Various validation statistics of the estimated AGB in Kalimantan, calculated using approximately 500 points per map, are listed in Table 3.The sample points were randomly distributed over areas, where reference data was available (training sites, Figure 1).Average AGB estimates are consistently higher than those obtained from the reference biomass, except to the AGB range >200 t/ha.Higher averages indicate a positive bias, while the magnitude of the bias is variable across the AGB ranges and points to the highest values in the biomass classes of 50-100 t/ha and 100-150 t/ha.In contrast, AGB ranges >200 t/ha display a negative bias.The root mean square error (RMSE) is similar in all years, with the highest relative errors exceeding 100% in the lower AGB ranges and low relative errors of around 22% in the highest AGB ranges.The distribution of the relative RMSE is similar in all three years.The relative overall RMSE ranges between 31% (2016), 36% (2009), and 38% (2007).The scatterplots of AGB estimates against the reference AGB display a similar distribution in all years (Figure 6).AGB values up to 250 t/ha show an overestimation of the models, especially in 2007, while AGB ranges higher than 250 t/ha indicate an underestimation.The coefficient of determination (R 2 ) varies in a range of 0.69 in 2016 to 0.77 in 2007 and the standard deviation (SD) is

Validation and Uncertainty
Various validation statistics of the estimated AGB in Kalimantan, calculated using approximately 500 points per map, are listed in Table 3.The sample points were randomly distributed over areas, where reference data was available (training sites, Figure 1).Average AGB estimates are consistently higher than those obtained from the reference biomass, except to the AGB range >200 t/ha.Higher averages indicate a positive bias, while the magnitude of the bias is variable across the AGB ranges and points to the highest values in the biomass classes of 50-100 t/ha and 100-150 t/ha.In contrast, AGB ranges >200 t/ha display a negative bias.The root mean square error (RMSE) is similar in all years, with the highest relative errors exceeding 100% in the lower AGB ranges and low relative errors of around 22% in the highest AGB ranges.The distribution of the relative RMSE is similar in all three years.The relative overall RMSE ranges between 31% (2016), 36% (2009), and 38% (2007).
The scatterplots of AGB estimates against the reference AGB display a similar distribution in all years (Figure 6).AGB values up to 250 t/ha show an overestimation of the models, especially in 2007, while AGB ranges higher than 250 t/ha indicate an underestimation.The coefficient of determination (R 2 ) varies in a range of 0.69 in 2016 to 0.77 in 2007 and the standard deviation (SD) is 53 t/ha in 2009 and 2016 to 56 t/ha in 2007.The NSE, indicating the efficiency of a model in a range from −∞ to 1, shows a good model performance for all years, reaching values between 0.70 and 0.76.   Figure 7 shows the number of sample points (frequency) of the reference AGB in contrast to the estimated AGB in each biomass range.The distribution of the frequency for the reference and the modelled biomass per range are similar, showing just small discrepancies in each year.The histogram of 2007 is affected by a smaller frequency of estimated biomass in each range, except the range >200 t/ha with a relative error of about 25%.In contrast, the histogram of 2009 is dominated by a lower sum of observations per range for all ranges except the smallest biomass range (0-50 t/ha) with a relative error of 42%.2016 shows only little differences in each biomass range.Figure 7 shows the number of sample points (frequency) of the reference AGB in contrast to the estimated AGB in each biomass range.The distribution of the frequency for the reference and the modelled biomass per range are similar, showing just small discrepancies in each year.The histogram of 2007 is affected by a smaller frequency of estimated biomass in each range, except the range >200 t/ha with a relative error of about 25%.In contrast, the histogram of 2009 is dominated by a lower sum of observations per range for all ranges except the smallest biomass range (0-50 t/ha) with a relative error of 42%.2016 shows only little differences in each biomass range.The uncertainty per AGB range (Figure 8), computed as a percentage of the AGB estimates, shows the highest relative errors in areas with low AGB estimates (Table 4).The relative errors in the biomass range 0-50 t/ha vary between 61.3% to 117.9%, in which the lowest values were obtained in 2009 and the highest in 2007.Areas with high biomass estimates show low relative errors of about 6% in each year.Additionally, the overall relative error per year show similar values varying in a range from 7.8% to 9.1%.The uncertainty per AGB range (Figure 8), computed as a percentage of the AGB estimates, shows the highest relative errors in areas with low AGB estimates (Table 4).The relative errors in the biomass range 0-50 t/ha vary between 61.3% to 117.9%, in which the lowest values were obtained in 2009 and the highest in 2007.Areas with high biomass estimates show low relative errors of about 6% in each year.Additionally, the overall relative error per year show similar values varying in a range from 7.8% to 9.1%.The uncertainty per AGB range (Figure 8), computed as a percentage of the AGB estimates, shows the highest relative errors in areas with low AGB estimates (Table 4).The relative errors in the biomass range 0-50 t/ha vary between 61.3% to 117.9%, in which the lowest values were obtained in 2009 and the highest in 2007.Areas with high biomass estimates show low relative errors of about 6% in each year.Additionally, the overall relative error per year show similar values varying in a range from 7.8% to 9.1%.

Biomass Estimation
The results suggest that a MLR using backscatter values in combination with ratios and textures of ALOS PALSAR-1/2 L-band and Sentinel-1 C-band data is eligible to model AGB in tropical forests.Englhart et al. [53], the authors found that in terms of biomass variability and saturation in tropical forest ecosystems MLR models are superior to ANN and SVR models.Since Kalimantan is dominated by high biomass values, a MLR was, therefore, used to model the AGB.
The use of highly accurate field inventory plots in combination with extensive LiDAR surveys as reference datasets improved the SAR-based biomass modelling.These datasets were collected in different ecosystems across Kalimantan, covering a wide range of ecosystems, vegetation types, forest structures and, thus, biomass ranges, for which reason they provide more precise biomass estimates than other sensors or the exclusive use of field inventory data [40].Nevertheless, uncertainties between field data and the modelled reference AGB can originate from different factors.First of all, the lag time between field data and LiDAR data acquisitions can introduce uncertainties due to regrowing or deforestation.For the reason that LiDAR data was not available for 2016, LiDAR data from 2011/2012 was adjusted as the reference layer for 2016.In order to minimize the effect of this temporal shift additional MODIS hotspot data (MCD14DL) was used to eliminate burned areas between 2011 and 2016.Certainly, forest cover changes resulting from logging or regrowing is not captured and can cause inaccuracies.Second, extrapolating field inventory data can introduce errors, since the spatial variability cannot be covered.In addition, the accuracy and precision of AGB extrapolated from field plots are affected by plot size and shape [59].Previous studies showed that errors in LiDAR-estimated AGB decrease exponentially with increasing field plot size, resulting from a smaller effect of co-registration errors and a spatial averaging of errors [58][59][60][61].In contrast, smaller plots have less overlap with the LiDAR data, cannot capture the variability of a forest, and are more sensitive to individual trees [13,62].For lack of data, differing field plot sizes and shapes were used to estimate the reference AGB.However, most of the used field plots exceed an area of 1000 m 2 , which is sufficiently large to be more robust against boundary effects [59,63,64].In contrast, Mauya et al. [65] show improvements of the model accuracy with increasing plot sizes in a range from 200-3000 m 2 .In addition to the size, the shape affects the precision and accuracy of the extrapolated biomass.Circular plots are less influenced by the circumference to area ratio than rectangular ones [65].As the number of rectangular plots in this study is limited and only within regrowing forest areas, their influence in AGB estimation is marginal.The uncertainty resulting from different sampling strategies are considered as the error of sampling size in the overall uncertainty of the SAR model with 20% [58].Finally, uncertainties may also be introduced by reason that a universal AGB model was applied for different forest types [66].Since tropical forests consist of hundreds of different tree species species-specific regression models cannot be applied [39].
After creating a reference biomass layer from LiDAR and field inventory data, AGB can be estimated for large-scaled areas using remote sensing data.The applied SAR backscatter approach is well known for radar-based forest cover and biomass mapping.It is computationally less intensive than alternative approaches and transferable to other regions, but also limited by some factors like backscatter saturation and backscatter variations due to terrain and wetness [19,67].The resulting maps point that the biomass variability due to different degradation stages in the forest, as well as disturbances in contrast to non-disturbed areas or clear-cuts that can be captured very well.The R 2 varies in a range of 0.69 in 2016 to 0.77 in 2007 and the NSE shows good model performance for all years, reaching values between 0.70 and 0.76.Studies modelling biomass and carbon based on ALOS PALSAR HV backscatter values in tropical forests found similar correlations, varying from 0.407 to 0.76 [68][69][70].Nevertheless, all maps show an underestimation at higher AGB levels compared to reference AGB.In addition to the underestimation in higher biomass ranges, biomass in lower AGB ranges is overestimated.Similar distributions of estimated biomass are shown in other regional studies estimating biomass from SAR backscatter values for different forest ecosystems, as a characteristic of SAR data for high biomass levels [17,66,71,72].
One of the general limitations of the applied approach is that SAR-based AGB retrieval suffers from saturation of the backscatter signal in the higher biomass range.The saturation level varies amongst others with the sensor wavelength and polarization, as well as the forest structure [20].AGB studies in tropical forests, such as in Kalimantan, were mostly conducted on the basis of L-band SAR data, being the most suitable operational data for biomass estimation [24][25][26][27][28].The saturation level in tropical forests, using the L-band, varies around 50 t/ha to 200 t/ha [25,30,72].Comparable to Thapa et al. [51], the present study uses backscatter values and, additionally, backscatter ratios and textures, which increases the saturation level to approximately 200-250 t/ha.A higher amount of ratios and textures based on backscatter values of different polarizations could further improve the model results, as shown in [51].With an increasing model complexity, the correlation and the NSE also increase in a range from 0.62 to 0.84 and from 0.54 to 0.83, respectively.
Another limiting factor of the backscatter approach are the moisture conditions of soil and vegetation [73].Especially in tropical forests, located in areas with a high amount of annual precipitation, the estimation of biomass based on backscatter can causes errors.To reduce humidity effects, scenes acquired during different periods of the year were averaged.Furthermore, images with a high influence of precipitation were excluded from the mosaics by incorporating TRMM data as selection criteria.Analogous to literature, variables based on HV backscatter were found as less influenced by changes in moisture and topography conditions at longer wavelengths and more sensitive to biomass than co-polarized data [21,25,30,34,70,72].However, even with the use of variables based on HV polarization, moisture effects can cause differences in the biomass estimation for the different years [74].The very dry conditions in 2007 may be an explanation for the higher backscatter saturation in 2007 and, thus, a better overall correlation.The lower R 2 of 2016 may not only result from the higher amount of precipitation in 2016, but also from the distribution of the validation samples.In contrast to the other years, the biomass range >200 t/ha contains 45% of all samples, which makes this class the most influential.The distinct underestimation of this AGB range is reflected in the overall R 2 of 2016.The given distribution of sample points originated because of the fire catastrophe in 2015, where a lot of areas within the reference biomass layer were burned.In order to reduce errors in the later modelling, those areas were excluded from the reference layer with the use of a MODIS hotspot product.The fires primarily occurred in areas that were already stressed by former burning or logging activities, accordingly in the lower biomass ranges, which leads to a lack of validation data in this class.
Furthermore, topography influences biomass estimation based on SAR data.To reduce possible errors caused by the effect of slope on radar backscatter, Mitchard et al. [26] suggested the use of a DEM, which is why we excluded AGB estimations in steep terrain with slopes >10 • , an area of 27% of the study area.Since the mountainous areas are barely affected by degradation and human activity, those regions were manually set to a biomass level of a healthy natural forest, 350 t/ha, as derived from LiDAR data.Additionally, settlement areas are not reliable in the final AGB maps.The high backscatter values due to double bounce effects result in backscatter values similar to those of forests.The use of an additional urban layer allows to flag settlements in a quality assurance layer, which we provided for each map.
Since models are minimizing the bias and overall error in order to get the best fit, the overestimation in lower biomass ranges results from the model adjustment due to the inability of estimating higher biomass ranges correctly [75].
In contrast to [33,76] a combination of C-and L-bands in the model of 2016 shows only minor improvement in comparison to a single-band approach.The C-band is more sensitive to variabilities in surface roughness resulting in improved modelling in burned areas or grass cover, so the combination of the C-band with the L-band slightly improves the correlation of estimated and predicted biomass compared to the use of the L-band alone.
The final biomass maps have a spatial resolution of 100 m, which is much finer than other existing large-scale biomass maps from Indonesia by [30], yet with a comparable accuracy.Additionally, the model is transferable to other countries and regions, if precise reference data is available.The results have important implications for carbon-related projects, such as UN-REDD, and can help to achieve the objectives through supporting monitoring or risk managing systems.The fine-scaled resolution and the sensitivity to AGB variability of our products allow capturing even small-scale logging activities, and the results can be used as an early warning for more extensive changes of forests due to logging and fire vulnerability.Furthermore, the more accurate biomass maps can lead to more precise carbon estimations when used as input data for carbon-related models.
The results of the biomass estimation in tropical forests could be further improved using biomass estimation approaches based on coherence or phase instead of backscatter, but these methods are limited by data availability.A combined data approach using optical data/vegetation indices and SAR data could also enhance biomass estimation.The launch of new P-band satellites, like ESA's Earth Explorer Biomass in 2021, will allow a more accurate estimation of biomass in tropical forests.

Temporal Change Estimation
To perform a change analysis a comparability of the three biomass maps is required.This is best accomplished using the same processing steps, model, reference, and radar data type for each map in order to avoid differences resulting from the methods and data.The three biomass maps of this study were generated using the same processing steps and model, but included Sentinel-1 as an additional sensor in the calculations for 2016.Differences between the three maps could also occur since different reference layers for each year were used.Nevertheless, the use of contemporaneous reference layers is essential, because of rapid changes in the landscape (e.g., fires and logging) and allows a more accurate calibration of the model.The use of two MLR models in 2016 may also cause discrepancies, since two variables were excluded from the second model.Using a confidence interval of 95% to identify changes allows a reduction of errors during the change mapping.
The resulting change maps show a significant loss of forest and, thus, biomass during 2007-2009, 2009-2016, and for the overall period 2007-2016.In addition to fires, changes can result due to illegal logging activities.The maps provide a helpful tool for REDD+, as well as national projects, since small-scaled deforestation is detected in an accurate and low-cost way for the entirety of Kalimantan over a period of approximately ten years using three time steps for the first time.

Conclusions
It was shown that a multivariate linear regression model using C-and L-band SAR-derived ratios and textures is able to model biomass more precisely and at a better resolution than existing models in this region.Nevertheless, the applied backscatter approach is limited by the fact that SAR-based AGB estimation is defined by saturation effects of the backscatter signal in higher biomass ranges.Due to the use of textures and the large amount of reference data, the saturation level could be increased to approximately 250 t/ha.Accordingly, the estimated AGB maps show an underestimation in higher AGB ranges, but also an overestimation in lower biomass ranges.The correlation of field biomass and estimated biomass varies in a span between 0.69 and 0.77.The model is able to capture biomass variability due to different degradation stages in forested areas.Additionally, different disturbances in contrast to non-disturbed area or clear-cuts can be identified.A biomass overestimation in urban areas and a reduced accuracy because of relief effects (layover, radar shadow) in steep slopes are known issues using SAR data and have been corrected via additional data.Sentinel-1 data, which was additionally used in 2016, did slightly improve the results.Modelling AGB for three different time steps allowed the estimation of change products.The change maps detected, for the first time, deforestation in an accurate and low-cost way for the entirety of Kalimantan over a period of approximately ten years based on three different time steps.The much higher spatial resolution of all products (100 m) allows capturing even small-scaled variabilities in forested areas, with the layer providing important

Figure 1 .
Figure 1.ALOS PALSAR false-color composite of Borneo Island using cross-and co-polarization as RGB (HH; HV; HH).To achieve the full coverage of Borneo 21 images of the year 2009 are mosaicked.

Figure 1 .
Figure 1.ALOS PALSAR false-color composite of Borneo Island using cross-and co-polarization as RGB (HH; HV; HH).To achieve the full coverage of Borneo 21 images of the year 2009 are mosaicked.
Pisau and Kapuas three circular subplots with radii of 20 m, 14 m, and 4 m, and 16 m, 8 m, and 4 m, and in Berau, Malinau and Kapuas Hulu four subplots with radii of 30/35 m, 20/25 m, 10 m, and 3 m were applied.In addition to circular nested plots, information was collected in rectangular plots of 20 m × 50 m to record saplings and trees in regrowing areas in Pulang Pisau and Kapuas.Moreover, nested rectangular plots with three subplots with sizes of 10 m × 10 m, 20 m × 20 m, and 20 m × 50 m were applied in Kapuas Hulu.

2. 2 .
Earth Observing Datasets 2.2.1.SAR Data SAR data was acquired in 2007, 2009, and 2016, using years with preferably dry conditions and less active fires.The study is based on L-band radar data of the Phased Array Type L-band Synthetic Aperture Radar (PALSAR), onboard the Advanced Land Observing Satellite (ALOS) of the Japan Aerospace Exploration Agency (JAXA) and ALOS-2 in combination with C-band radar data based on ENVISAT ASAR and Sentinel-1.ALOS PALSAR and ALOS-2 PALSAR-2 mosaics with a spatial resolution of 25 m were provided by the Kyoto and Carbon Initiative.ENVISAT ASAR was launched in 2002 by the European Space Agency (ESA).Sentinel-1a and Sentinel-1b C-band radars were launched within the Copernicus program by the ESA in 2014 and 2016, respectively.ALOS PALSAR (HH, HV) images were used to estimate the biomass in Kalimantan for the years 2007 and 2009.As C-Band radar data for the years 2007 and 2009 imageries of ENVISAT ASAR (VV) data were tested, certainly the data had to be excluded as the study area is not fully covered by scenes acquired in one sensor mode and because of strong moisture effects.In 2016, data of ALOS-2, launched in 2014, in combination with Sentinel-1 GRD data, acquired in interferometric wide (IW) swath and a resolution of 10 m × 10 m, was used.

Figure 2 .
Figure 2. Flow chart of the AGB map processing steps.Blue refers to the results, and green refers to the main result.

Figure 3 .
Figure 3. Final biomass maps for Kalimantan, per year and subset of the quantity of change layer for 2009-2016.Biomass ranges from 0-350 t/ha, and the quantity of change ranges from −300-100 t/ha.

Figure 3 .
Figure 3. Final biomass maps for Kalimantan, per year and subset of the quantity of change layer for 2009-2016.Biomass ranges from 0-350 t/ha, and the quantity of change ranges from −300-100 t/ha.

Figure 3 .
Figure 3. Final biomass maps for Kalimantan, per year and subset of the quantity of change layer for 2009-2016.Biomass ranges from 0-350 t/ha, and the quantity of change ranges from −300-100 t/ha.

Figure 5 .
Figure 5.Comparison between different AGB maps (a) AGB MLR model result for 2009; (b) LiDAR; (c) Saatchi et al. (2011); (d) Baccini et al. (2012); and (e) Avitabile et al. (2016) in a range of 0-350 t/ha for a subset in the south of Kalimantan.The resolution of AGB maps of the MLR model and the LiDAR data is 100 m; the map of Saatchi et al., 1 km; the map of Baccini et al., 500 m; and the map of Avitabile et al., 1 km.

Figure 5 .
Figure 5.Comparison between different AGB maps (a) AGB MLR model result for 2009; (b) LiDAR; (c) Saatchi et al. (2011); (d) Baccini et al. (2012); and (e) Avitabile et al. (2016) in a range of 0-350 t/ha for a subset in the south of Kalimantan.The resolution of AGB maps of the MLR model and the LiDAR data is 100 m; the map of Saatchi et al., 1 km; the map of Baccini et al., 500 m; and the map of Avitabile et al., 1 km.

Figure 6 .
Figure 6.Linear regression of estimated above-ground biomass and reference above-ground biomass using approximately 500 randomly-generated points across the reference layer for each year (red dashed line = 1:1 line; black line = linear trend including confidence bounds).

Figure 6 .
Figure 6.Linear regression of estimated above-ground biomass and reference above-ground biomass using approximately 500 randomly-generated points across the reference layer for each year (red dashed line = 1:1 line; black line = linear trend including confidence bounds).

Figure 7 .
Figure 7. Histograms of reference above-ground biomass and estimated above-ground biomass per biomass range for each year.Frequency refers to the number of observations per range.

Figure 8 .
Figure 8. Uncertainty maps per AGB range for each year showing highest errors in low biomass ranges.

Figure 7 .
Figure 7. Histograms of reference above-ground biomass and estimated above-ground biomass per biomass range for each year.Frequency refers to the number of observations per range.

Figure 7 .
Figure 7. Histograms of reference above-ground biomass and estimated above-ground biomass per biomass range for each year.Frequency refers to the number of observations per range.

Figure 8 .
Figure 8. Uncertainty maps per AGB range for each year showing highest errors in low biomass ranges.Figure 8. Uncertainty maps per AGB range for each year showing highest errors in low biomass ranges.

Figure 8 .
Figure 8. Uncertainty maps per AGB range for each year showing highest errors in low biomass ranges.Figure 8. Uncertainty maps per AGB range for each year showing highest errors in low biomass ranges.

Table 1 .
Overview of the measured field inventory plots (473) and the acquired LiDAR data (8300 km 2 ) in Kalimantan.

Table 2 .
Overview of the used variables for the different MLR models per year and the R 2 and residual standard error (RSE) for the calibration of the model (RP = higher texture run percentage).

Table 3 .
Overview of validation statistics per AGB ranges and year.

Table 3 .
Overview of validation statistics per AGB ranges and year.

Table 4 .
Uncertainty per AGB range and year in percentage.

Table 4 .
Uncertainty per AGB range and year in percentage.

Table 4 .
Uncertainty per AGB range and year in percentage.