Measuring biomass changes due to woody encroachment and deforestation/degradation in a forest-savanna boundary region of central Africa using multi-temporal L-band radar backscatter

Satellite L-band synthetic aperture radar backscatter data from 1996 and 2007 (from JERS-1 and ALOS PALSAR respectively), were used with field data collected in 2007 and a back-calibration method to produce biomass maps of a 15 000 km² forest-savanna ecotone region of central Cameroon. The relationship between the radar backscatter and aboveground biomass (AGB) was strong (r² = 0.86 for ALOS HV to biomass plots, r² = 0.95 relating ALOS-derived biomass for 40 suspected unchanged regions to JERS-1 HH). The root mean square error (RMSE) associated with AGB estimation varied from ~ 25% for AGB < 100 Mg ha^− 1 to ~ 40% for AGB > 100 Mg ha^− 1 for the ALOS HV data. Change detection showed a significant loss of AGB over high biomass forests, due to suspected deforestation and degradation, and significant biomass gains along the forest-savanna boundary, particularly in areas of low population density. Analysis of the errors involved showed that radar data can detect changes in broad AGB class in forest-savanna transition areas with an accuracy > 95%. However, quantitative assessment of changes in AGB in Mg ha^− 1 at a pixel level will require radar images from sensors with similar characteristics collecting data from the same season over multiple years.     

Empirical relationships between AIRSAR backscatter and LiDAR-derived forest biomass, Queensland, Australia

To evaluate the use of multi-frequency, polarimetric Synthetic Aperture Radar (SAR) data for quantifying the above ground biomass (AGB) of open forests and woodlands, NASA JPL AIRSAR (POLSAR) data were acquired over a 37 × 60 km area west of Injune, central Queensland, Australia. From field measurements recorded within 32 50 × 50 m plots, AGB was estimated by applying species-specific allometric equations to stand measurements. AGB was then scaled-up to the larger area using relationships established with Light Detection and Ranging (LiDAR) data acquired over 150 (10 columns, 15 rows) 500 × 150 m cells (or Primary Sampling Units, PSUs) spaced 4 × 4 km apart in the north- and east-west directions. Large-scale (1 : 4000) stereo aerial photographs were also acquired for each PSU to assess species composition. Based on the LiDAR extrapolations, the median AGB for the PSU grid was 82 Mg ha^− 1 (maximum 164 Mg ha^− 1), with the higher levels associated with forests containing a high proportion of Angophora and Callitris species. Empirical relationships between AGB and SAR backscatter confirmed that C-, L- and P-band saturated at different levels and revealed a greater strength in the relationship at higher incidence angles and a larger dynamic range and consistency of relationships at HV polarizations. A higher level of saturation (above ∼50 Mg ha^− 1) was observed at C-band HV compared to that reported for closed forests which was attributable to a link between foliage projected cover (FPC) and AGB. The study concludes that L-band HV backscatter data acquired at incidence angles approaching or exceeding 45° are best suited for estimating the AGB up to the saturation level of ∼80-85 Mg ha^− 1. For regional mapping of biomass below the level of saturation, the use of the Japanese Space Exploration Agency (JAXA) Advanced Land Observing Satellite (ALOS) Phase Arrayed L-band SAR (PALSAR) is advocated.     

L- and P-band backscatter intensity for biomass retrieval in hemiboreal forest

At present, the greatest source of uncertainty in the global carbon cycle is in the terrestrial ecosystems. In order to reduce these uncertainties it is necessary to provide consistent and accurate global estimates of the world forest biomass. One of the most promising methods for obtaining such estimates is through polarimetric SAR backscatter measurements at low frequencies. In this paper, the relation between polarimetric SAR backscatter at L- and P-bands and forest biomass is investigated using data acquired within the BioSAR-I campaign in southern Sweden during 2007. Methods for estimating biomass on stand level using these data are developed and evaluated, and the results for the two frequency bands are compared. For L-band data, the best results were obtained using HV-polarized backscatter only, giving estimation errors in terms of root mean square errors (RMSE) between 31% and 46% of the mean biomass for stands with biomass ranging from 10 to 290 t/ha, and an (adjusted) coefficient of determination (R²) between 0.4 and 0.6. For P-band data, the results are better than for L-band. Models using HV- or HH-polarized P-band backscatter give similar results, as does a model including both HV and HH. The RMSEs were between 18 and 27%, and the R² values were between 0.7 and 0.8.     

Integration of radar and Landsat-derived foliage projected cover for woody regrowth mapping, Queensland, Australia

In Queensland, Australia, forest areas are discriminated from non-forest by applying a threshold (∼ 12%) to Landsat-derived Foliage Projected Cover (FPC) layers (equating to ∼ 20% canopy cover), which are produced routinely for the State. However, separation of woody regrowth following agricultural clearing cannot be undertaken with confidence, and is therefore not mapped routinely by State Agencies. Using fully polarimetric C-, L- and P-band NASA AIRSAR and Landsat FPC data for forests and agricultural land near Injune, central Queensland, we corroborate that woody regrowth dominated by Brigalow (Acacia harpophylla) cannot be discriminated using either FPC or indeed C-band data alone, because the rapid attainment of a canopy cover leads to similarities in both reflectance and backscatter with remnant forest. We also show that regrowth cannot be discriminated from non-forest areas using either L-band or P-band data alone. However, mapping can be achieved by thresholding and intersecting these layers, as regrowth is unique in supporting both a high FPC (> ∼ 12%) and C-band SAR backscatter (> ~ − 18 dB at HV polarisation) and low L-band and P-band SAR backscatter (e.g. < =∼ 14 dB at L-band HH polarisation). To provide a theoretical explanation, a wave scattering model based on that of Durden et al. [Durden, S.L., Van Zyl, J.J. & Zebker, H.A. (1989). Modelling and observation of radar polarization signature of forested areas. IEEE Trans. Geoscience and Remote Sensing, 27, 290-301.] was used to demonstrate that volume scattering from leaves and small branches in the upper canopy leads to increases in C-band backscattering (particularly HV polarisations) from regrowth, which increases proportionally with FPC. By contrast, low L-band and P-band backscatter occurs because of the lack of double bounce interactions at co-polarisations (particularly HH) and volume scattering at HV polarisation from the stems and branches, respectively, when their dimensions are smaller than the wavelength. Regrowth maps generated by applying simple thresholds to both FPC and AIRSAR L-band data showed a very close correspondence with those mapped using same-date 2.5 m Hymap data and an average 73.7% overlap with those mapped through time-series comparison of Landsat-derived land cover classifications. Regrowth mapped using Landsat-derived FPC from 1995 and JER-1 SAR data from 1994-1995 also corresponded with areas identified within the time-series classification and true colour stereo photographs for the same period. The integration of Landsat FPC and L-band SAR data is therefore expected to facilitate regrowth mapping across Queensland and other regions of Australia, particularly as Japan's Advanced Land Observing System (ALOS) Phase Arrayed L-band SAR (PALSAR), to be launched in 2006, will observe at both L-band HH and HV polarisations.     

Aboveground biomass retrieval in tropical forests — The potential of combined X- and L-band SAR data use

In the context of reducing emissions from deforestation and forest degradation (REDD) and the international effort to reduce anthropogenic greenhouse gas emissions, a reliable assessment of aboveground forest biomass is a major requirement. Especially in tropical forests which store huge amounts of carbon, a precise quantification of aboveground biomass is of high relevance for REDD activities. This study investigates the potential of X- and L-band SAR data to estimate aboveground biomass (AGB) in intact and degraded tropical forests in Central Kalimantan, Borneo, Indonesia. Based on forest inventory data, aboveground biomass was first estimated using LiDAR data. These results were then used to calibrate SAR backscatter images and to upscale the biomass estimates across large areas and ecosystems. This upscaling approach not only provided aboveground biomass estimates over the whole biomass range from woody regrowth to mature pristine forest but also revealed a spatial variation due to varying growth condition within specific forest types. Single and combined frequencies, as well as mono- and multi-temporal TerraSAR-X and ALOS PALSAR biomass estimation models were analyzed for the development of accurate biomass estimations. Regarding the single frequency analysis overall ALOS PALSAR backscatter is more sensitive to AGB than TerraSAR-X, especially in the higher biomass range (> 100 t/ha). However, ALOS PALSAR results were less accurate in low biomass ranges due to a higher variance. The multi-temporal L- and X-band combined model achieved the best result and was therefore tested for its temporal and spatial transferability. The achieved accuracy for this model using nearly 400 independent validation points was r² = 0.53 with an RMSE of 79 t/ha. The model is valid up to 307 t/ha with an accuracy requirement of 50 t/ha and up to 614 t/ha with an accuracy requirement of 100 t/ha in flat terrain. The results demonstrate that direct biomass measurements based on the synergistic use of L- and X-band SAR can provide large-scale AGB estimations for tropical forests. In the context of REDD monitoring the results can be used for the assessment of the spatial distribution of the biomass, also indicating trends in high biomass ranges and the characterization of the spatial patterns in different forest types.     

Exploring LiDAR-RaDAR synergy—predicting aboveground biomass in a southwestern ponderosa pine forest using LiDAR, SAR and InSAR

Scanning Light Detecting and Ranging (LiDAR), Synthetic Aperture Radar (SAR) and Interferometric SAR (InSAR) were analyzed to determine (1) which of the three sensor systems most accurately predicted forest biomass, and (2) if LiDAR and SAR/InSAR data sets, jointly considered, produced more accurate, precise results relative to those same data sets considered separately. LiDAR ranging measurements, VHF-SAR cross-sectional returns, and X- and P-band cross-sectional returns and interferometric ranges were regressed with ground-estimated (from dbh) forest biomass in ponderosa pine forests in the southwestern United States. All models were cross-validated. Results indicated that the average canopy height measured by the scanning LiDAR produced the best predictive equation. The simple linear LiDAR equation explained 83% of the biomass variability (n = 52 plots) with a cross-validated root mean square error of 26.0 t/ha. Additional LiDAR metrics were not significant to the model. The GeoSAR P-band (λ = 86 cm) cross-sectional return and the GeoSAR/InSAR canopy height (X-P) captured 30% of the forest biomass variation with an average predictive error of 52.5 t/ha. A second RaDAR-FOPEN collected VHF (λ ∼ 7.8 m) and cross-polarized P-band (λ = 88 cm) cross-sectional returns, none of which proved useful for forest biomass estimation (cross-validated R² = 0.09, RMSE = 63.7 t/ha). Joint consideration of LiDAR and RaDAR measurements produced a statistically significant, albeit small improvement in biomass estimation precision. The cross-validated R² increased from 83% to 84% and the prediction error decreased from 26.0 t/ha to 24.9 t/ha when the GeoSAR X-P interferometric height is considered along with the average LiDAR canopy height. Inclusion of a third LiDAR metric, the 60th decile height, further increased the R² to 85% and decreased the RMSE to 24.1 t/ha. On this 11 km² ponderosa pine study area, LiDAR data proved most useful for predicting forest biomass. RaDAR ranging measurements did not improve the LiDAR estimates.     

Forest carbon densities and uncertainties from Lidar, QuickBird, and field measurements in California

Greenhouse gas inventories and emissions reduction programs require robust methods to quantify carbon sequestration in forests. We compare forest carbon estimates from Light Detection and Ranging (Lidar) data and QuickBird high-resolution satellite images, calibrated and validated by field measurements of individual trees. We conducted the tests at two sites in California: (1) 59 km² of secondary and old-growth coast redwood (Sequoia sempervirens) forest (Garcia-Mailliard area) and (2) 58 km² of old-growth Sierra Nevada forest (North Yuba area). Regression of aboveground live tree carbon density, calculated from field measurements, against Lidar height metrics and against QuickBird-derived tree crown diameter generated equations of carbon density as a function of the remote sensing parameters. Employing Monte Carlo methods, we quantified uncertainties of forest carbon estimates from uncertainties in field measurements, remote sensing accuracy, biomass regression equations, and spatial autocorrelation. Validation of QuickBird crown diameters against field measurements of the same trees showed significant correlation (r = 0.82, P < 0.05). Comparison of stand-level Lidar height metrics with field-derived Lorey's mean height showed significant correlation (Garcia-Mailliard r = 0.94, P < 0.0001; North Yuba R = 0.89, P < 0.0001). Field measurements of five aboveground carbon pools (live trees, dead trees, shrubs, coarse woody debris, and litter) yielded aboveground carbon densities (mean ± standard error without Monte Carlo) as high as 320 ± 35 Mg ha^− 1 (old-growth coast redwood) and 510 ± 120 Mg ha^− 1 (red fir [Abies magnifica] forest), as great or greater than tropical rainforest. Lidar and QuickBird detected aboveground carbon in live trees, 70-97% of the total. Large sample sizes in the Monte Carlo analyses of remote sensing data generated low estimates of uncertainty. Lidar showed lower uncertainty and higher accuracy than QuickBird, due to high correlation of biomass to height and undercounting of trees by the crown detection algorithm. Lidar achieved uncertainties of < 1%, providing estimates of aboveground live tree carbon density (mean ± 95% confidence interval with Monte Carlo) of 82 ± 0.7 Mg ha^− 1 in Garcia-Mailliard and 140 ± 0.9 Mg ha^− 1 in North Yuba. The method that we tested, combining field measurements, Lidar, and Monte Carlo, can produce robust wall-to-wall spatial data on forest carbon.     

Model-assisted regional forest biomass estimation using LiDAR and InSAR as auxiliary data: A case study from a boreal forest area

There is a need for accurate inventory methods that produce relevant and timely information on the forest resources and carbon stocks for forest management planning and for implementation of national strategies under the United Nations Collaborative Program on Reduced Emissions from Deforestation and Forest Degradation in Developing Countries (REDD). Such methods should produce information that is consistent across various geographical scales. Airborne scanning Light Detection and Ranging (LiDAR) is among the most promising remote sensing technologies for estimation of forest resource information such as timber volume and biomass, while acquisition of three dimensional data with Interferometric Synthetic Aperture Radar (InSAR) from space is seen as a relevant option for inventory in the tropics because of its ability to “see through the clouds” and its potential for frequent updates at low costs. Based on a stratified probability sample of 201 field survey plots collected in a 960 km² boreal forest area in Norway, we demonstrate how total above-ground biomass (AGB) can be estimated at three distinct geographical levels in such a way that the estimates at a smaller level always sum up to the estimate at a larger level. The three levels are (1) a district (the entire study area), (2) a village, local community or estate level, and (3) a stand or patch level. The LiDAR and InSAR data were treated as auxiliary information in the estimation. At the two largest geographical levels model-assisted estimators were employed. A model-based estimation was conducted at the smallest level. Estimates of AGB and corresponding error estimates based on (1) the field sample survey were compared with estimates obtained by using (2) LiDAR and (3) InSAR data as auxiliary information. For the entire study area, the estimates of AGB were 116.0, 101.2, and 111.3 Mg ha⁻¹, respectively. Corresponding standard error estimates were 3.7, 1.6, and 3.2 Mg ha⁻¹. At the smallest geographical level (stand) an independent validation on 35 large field plots was carried out. RMSE values of 17.1-17.3 Mg ha⁻¹ and 42.6-53.2 Mg ha⁻¹ were found for LiDAR and InSAR, respectively. A time lag of six years between acquisition of InSAR data and field inventory has introduced some errors. Significant differences between estimates and reference values were found, illustrating the risk of using pure model-based methods in the estimation when there is a lack of fit in the models. We conclude that the examined remote sensing techniques can provide biomass estimates with smaller estimated errors than a field-based sample survey. The improvement can be highly significant, especially for LiDAR.     

Distinguishing between live and dead standing tree biomass on the North Rim of Grand Canyon National Park, USA using small-footprint lidar data

Accurate estimation of live and dead biomass in forested ecosystems is important for studies of carbon dynamics, biodiversity, wildfire behavior, and for forest management. Lidar remote sensing has been used successfully to estimate live biomass, but studies focusing on dead biomass are rare. We used lidar data, in conjunction with field measurements from 58 plots to distinguish between and map standing live and dead tree biomass in the mixed coniferous forest of the North Rim of Grand Canyon National Park, USA. Lidar intensity and canopy volume were key variables for estimating live biomass, whereas for dead biomass, lidar intensity alone was critical for accurate estimation. Regression estimates of both live and dead biomass ranged between 0 and 600 Mg ha^− 1, with means of 195.08 Mg ha^− 1 and 65.73 Mg ha^− 1, respectively. Cross validation with field data resulted in correlation coefficients for predicted vs. observed of 0.85 for live biomass (RMSE = 50 Mg ha^− 1 and %RMSE (RMSE as a percent of the mean) = 26). For dead biomass, correlation was 0.79, RMSE was 42 Mg ha^− 1, and %RMSE was 63. Biomass maps revealed interesting patterns of live and dead standing tree biomass. Live biomass was highest in the ponderosa pine zone, and decreased from south to north through the mixed conifer and spruce-fir forest zones. Dead biomass exhibited a background range of values in these mature forests from zero to 100 Mg ha^− 1, with lower values in locations having higher live biomass. In areas with high dead biomass values, live biomass was near zero. These areas were associated with recent wildfires, as indicated by fire maps derived from the Monitoring Trends in Burn Severity Project (MTBS). Combining our dead biomass maps with the MTBS maps, we demonstrated the complementary power of these two datasets, revealing that MTBS burn intensity class can be described quantitatively in terms of dead biomass. Assuming a background range of dead biomass up to 100 Mg ha^− 1, it is possible to estimate and map the contribution to the standing dead tree biomass pool associated with recent wildfire.     

Improved forest biomass estimates using ALOS AVNIR-2 texture indices

Optical remote sensing is still one of the most attractive choices for obtaining biomass information, as new sensors are available with fine spatial and spectral resolutions. Better biomass estimates may be possible if suitable processing techniques for these sensors can be demonstrated. This research investigates the potential of high resolution optical data from the ALOS AVNIR-2 sensor for biomass estimation in a mountainous, subtropical forested region using four different types of image processing techniques including i) spectral reflectance and simple spectral band ratio, ii) commonly used vegetation indices, iii) texture parameters and iv) ratio of texture parameters. Simple linear and stepwise multiple regression models were developed between biomass data from 50 field plots, and image parameters derived from these techniques.Results indicate that spectral reflectance, the simple band ratio, and commonly used vegetation indices have relatively low potential for biomass estimation, as only about 58% of the variability in the field data was explained by the model (adjusted r² = 0.58 and RMSE = 64 t/ha). However, the texture parameters of spectral bands were found to be effective for biomass estimation with an explained variability of ca. 76% (adjusted r² = 0.76 and RMSE = 46 t/ha). The result was further improved to adjusted r² = 0.88 (RMSE = 32 t/ha) using the simple ratio of texture parameters. The results suggest that the performance of biomass estimation can be improved significantly using the texture parameters of high resolution optical data, and further improvement can be obtained using the ratio of texture parameters, as this combines the advantages of both texture and ratio.     

The BIOMASS mission: Mapping global forest biomass to better understand the terrestrial carbon cycle

In response to the urgent need for improved mapping of global biomass and the lack of any current space systems capable of addressing this need, the BIOMASS mission was proposed to the European Space Agency for the third cycle of Earth Explorer Core missions and was selected for Feasibility Study (Phase A) in March 2009. The objectives of the mission are 1) to quantify the magnitude and distribution of forest biomass globally to improve resource assessment, carbon accounting and carbon models, and 2) to monitor and quantify changes in terrestrial forest biomass globally, on an annual basis or better, leading to improved estimates of terrestrial carbon sources (primarily from deforestation); and terrestrial carbon sinks due to forest regrowth and afforestation. These science objectives require the mission to measure above-ground forest biomass from 70° N to 56° S at spatial scale of 100-200 m, with error not exceeding ± 20% or ± 10 t ha^− 1 and forest height with error of ± 4 m. To meet the measurement requirements, the mission will carry a P-Band polarimetric SAR (centre frequency 435 MHz with 6 MHz bandwidth) with interferometric capability, operating in a dawn-dusk orbit with a constant incidence angle (in the range of 25°-35°) and a 25-45 day repeat cycle. During its 5-year lifetime, the mission will be capable of providing both direct measurements of biomass derived from intensity data and measurements of forest height derived from polarimetric interferometry. The design of the BIOMASS mission spins together two main observational strands: (1) the long heritage of airborne observations in tropical, temperate and boreal forest that have demonstrated the capabilities of P-band SAR for measuring forest biomass; (2) new developments in recovery of forest structure including forest height from Pol-InSAR, and, crucially, the resistance of P-band to temporal decorrelation, which makes this frequency uniquely suitable for biomass measurements with a single repeat-pass satellite. These two complementary measurement approaches are combined in the single BIOMASS sensor, and have the satisfying property that increasing biomass reduces the sensitivity of the former approach while increasing the sensitivity of the latter. This paper surveys the body of evidence built up over the last decade, from a wide range of airborne experiments, which illustrates the ability of such a sensor to provide the required measurements.At present, the BIOMASS P-band radar appears to be the only sensor capable of providing the necessary global knowledge about the world's forest biomass and its changes. In addition, this first chance to explore the Earth's environment with a long wavelength satellite SAR is expected to make yield new information in a range of geoscience areas, including subsurface structure in arid lands and polar ice, and forest inundation dynamics.     

Biomass estimation over a large area based on standwise forest inventory data and ASTER and MODIS satellite data: A possibility to verify carbon inventories

According to the IPCC GPG (Intergovernmental Panel on Climate Change, Good Practice Guidance), remote sensing methods are especially suitable for independent verification of the national LULUCF (Land Use, Land-Use Change, and Forestry) carbon pool estimates, particularly the aboveground biomass. In the present study, we demonstrate the potential of standwise (forest stand is a homogenous forest unit with average size of 1-3 ha) forest inventory data, and ASTER and MODIS satellite data for estimating stand volume (m³ ha^− 1) and aboveground biomass (t ha^− 1) over a large area of boreal forests in southern Finland. The regression models, developed using standwise forest inventory data and standwise averages of moderate spatial resolution ASTER data (15 m × 15 m), were utilized to estimate stand volume for coarse resolution MODIS pixels (250 m × 250 m). The MODIS datasets for three 8-day periods produced slightly different predictions, but the averaged MODIS data produced the most accurate estimates. The inaccuracy in radiometric calibration between the datasets, the effect of gridding and compositing artifacts and phenological variability are the most probable reasons for this variability. Averaging of the several MODIS datasets seems to be one possibility to reduce bias. The estimates obtained were significantly close to the district-level mean values provided by the Finnish National Forest Inventory; the relative RMSE was 9.9%. The use of finer spatial resolution data is an essential step to integrate ground measurements with coarse spatial resolution data. Furthermore, the use of standwise forest inventory data reduces co-registration errors and helps in solving the scaling problem between the datasets. The approach employed here can be used for estimating the stand volume and biomass, and as required independent verification data.     

A LiDAR-derived canopy density model for tree stem and crown mapping in Australian forests

The retrieval of tree and forest structural attributes from Light Detection and Ranging (LiDAR) data has focused largely on utilising canopy height models, but these have proved only partially useful for mapping and attributing stems in complex, multi-layered forests. As a complementary approach, this paper presents a new index, termed the Height-Scaled Crown Openness Index (HSCOI), which provides a quantitative measure of the relative penetration of LiDAR pulses into the canopy. The HSCOI was developed from small footprint discrete return LiDAR data acquired over mixed species woodlands and open forests near Injune, Queensland, Australia, and allowed individual trees to be located (including those in the sub-canopy) and attributed with height using relationships (r² = 0.81, RMSE = 1.85 m, n = 115; 4 outliers removed) established with field data. A threshold contour of the HSCOI surface that encompassed ∼ 90% of LiDAR vegetation returns also facilitated mapping of forest areas, delineation of tree crowns and clusters, and estimation of canopy cover. At a stand level, tree density compared well with field measurements (r² = 0.82, RMSE = 133 stems ha^− 1, n = 30), with the most consistent results observed for stem densities ≤ 700 stems ha^− 1. By combining information extracted from both the HSCOI and the canopy height model, predominant stem height (r² = 0.91, RMSE = 0.77 m, n = 30), crown cover (r² = 0.78, RMSE = 9.25%, n = 30), and Foliage & Branch Projective Cover (FBPC; r² = 0.89, RMSE = 5.49%, n = 30) were estimated to levels sufficient for inventory of woodland and open forest structural types. When the approach was applied to forests in north east Victoria, stem density and crown cover were reliably estimated for forests with a structure similar to those observed in Queensland, but less so for forests of greater height and canopy closure.     

The use of waveform lidar to measure northern temperate mixed conifer and deciduous forest structure in New Hampshire

The direct retrieval of canopy height and the estimation of aboveground biomass are two important measures of forest structure that can be quantified by airborne laser scanning at landscape scales. These and other metrics are central to studies attempting to quantify global carbon cycles and to improve understanding of the spatial variation in forest structure evident within differing biomes. Data acquired using NASA's Laser Vegetation Imaging Sensor (LVIS) over the Bartlett Experimental Forest (BEF) in central New Hampshire (USA) was used to assess the performance of waveform lidar in a northern temperate mixed conifer and deciduous forest.Using coincident plots established for this study, we found strong agreement between field and lidar measurements of height (r² = 0.80, p < 0.000) at the footprint level. Allometric calculations of aboveground biomass (AGBM) and LVIS metrics (AGBM: r² = 0.61, PRESS RMSE = 58.0 Mg ha^− 1, p < 0.000) and quadratic mean stem diameter (QMSD) and LVIS metrics (r² = 0.54, p = 0.002) also showed good agreement at the footprint level. Application of a generalized equation for determining AGBM proposed by Lefsky et al. (2002a) to footprint-level field data from Bartlett resulted in a coefficient of determination of 0.55; RMSE = 64.4 Mg ha^− 1; p = 0.002. This is slightly weaker than the strongest relationship found with the best-fit single term regression model.Relationships between a permanent grid of USDA Forest Service inventory plots and the mean values of aggregated LVIS metrics, however, were not as strong. This discrepancy suggests that validation efforts must be cautious in using pre-existing field data networks as a sole means of calibrating and verifying such remote sensing data. Stratification based on land-use or species composition, however, did provide the means to improve regression relationships at this scale. Regression models established at the footprint level for AGBM and QMSD were applied to LVIS data to generate predicted values for the whole of Bartlett. The accuracy of these models was assessed using varying subsets of the USFS NERS plot data. Coefficient of determinations ranged from fair to strong with aspects of land-use history and species composition influencing both the fit and the level of error seen in the predicted relationships.     

Estimating biomass carbon stocks for a Mediterranean forest in central Spain using LiDAR height and intensity data

Biomass fractions (total aboveground, branches and foliage) were estimated from a small footprint discrete-return LiDAR system in an unmanaged Mediterranean forest in central Spain. Several biomass estimation models based on LiDAR height, intensity or height combined with intensity data were explored. Raw intensity data were normalized to a standard range in order to remove the range dependence of the intensity signal. In general terms, intensity-based models provided more accurate predictions of the biomass fractions. Height models selected were mainly based on a percentile of the height distribution. Intensity models selected included variables that consider the percentage of the intensity accumulated at different height percentiles, which implicitly take into account the height distribution. The general models derived considering all species together were based on height combined with intensity data. These models yielded R² values greater than 0.58 for the different biomass fractions considered and RMSE values of 28.89, 18.28 and 1.51 Mg ha⁻¹ for aboveground, branch and foliage biomass, respectively. Results greatly improved for species-specific models using the main species present in each plot, with R² values greater than 0.85, 0.70 and 0.90 for black pine, Spanish juniper and Holm oak, respectively, and with lower RMSE for the biomass fractions. Reductions in LiDAR point density had only a small effect on the results obtained, except for those models based on a variation of the Canopy Reflection Sum, which was weighted by the mean point density. Based on the species-specific equations derived, Holm oak dominated plots showed the highest average carbon contained by aboveground biomass and branch biomass 44.66 and 31.42 Mg ha^− 1 respectively, while for foliage biomass carbon, Spanish juniper showed the highest average value (3.04 Mg ha^− 1).     

Mapping biomass and stress in the Sierra Nevada using lidar and hyperspectral data fusion

In this paper, we explored fusion of structural metrics from the Laser Vegetation Imaging Sensor (LVIS) and spectral characteristics from the Airborne Visible Infrared Imaging Spectrometer (AVIRIS) for biomass estimation in the Sierra Nevada. In addition, we combined the two sensors to map species-specific biomass and stress at landscape scale. Multiple endmember spectral mixture analysis (MESMA) was used to classify vegetation from AVIRIS images and obtain sub-pixel fractions of green vegetation, non-photosynthetic vegetation, soil, and shade. LVIS metrics, AVIRIS spectral indices, and MESMA fractions were compared with field measures of biomass using linear and stepwise regressions at stand (1 ha) level. AVIRIS metrics such as water band indices and shade fractions showed strong correlation with LVIS canopy height (r² = 0.69, RMSE = 5.2 m) and explained around 60% variability in biomass. LVIS variables were found to be consistently good predictors of total and species specific biomass (r² = 0.77, RMSE = 70.12 Mg/ha). Prediction by LVIS after species stratification of field data reduced errors by 12% (r² = 0.84, RMSE = 58.78 Mg/ha) over using LVIS metrics alone. Species-specific biomass maps and associated errors created from fusion were different from those produced without fusion, particularly for hardwoods and pines, although mean biomass differences between the two techniques were not statistically significant. A combined analysis of spatial maps from LVIS and AVIRIS showed increased water and chlorophyll stress in several high biomass stands in the study area. This study provides further evidence that lidar is better suited for biomass estimation, per se, while the best use of hyperspectral data may be to refine biomass predictions through a priori species stratification, while also providing information on canopy state, such as stress. Together, the two sensors have many potential applications in carbon dynamics, ecological and habitat studies.     

Investigating RaDAR-LiDAR synergy in a North Carolina pine forest

A low frequency (80-120 MHz) VHF RaDAR, BioSAR, specifically designed for forest biomass estimation and a profiling LiDAR, PALS, were flown over loblolly pine plantations in the southeastern United States. LiDAR-only, RaDAR-only, and joint LiDAR-RaDAR linear models were developed to determine if returns from two sensors could be used to estimate pine biomass more accurately and precisely than returns from either sensor alone. The best five-variable RaDAR model explained 81.8% (R²) of the stem green biomass variability, with a regression RMSE of 57.5 t/ha. The best one-variable LiDAR model explained 93.3% of the biomass variation (RMSE = 33.9 t/ha). Combining the RaDAR normalized volumetric returns with the profiling LiDAR ranging measurements did little to improve the best LiDAR-only model. The best LiDAR-RaDAR model explained 93.8% of the biomass variation (RSME = 32.7 t/ha). Cross-validation and training/test validation procedures demonstrated (1) that all models are unbiased and (2) the increased precision of the LiDAR-only and LiDAR-RaDAR models. The results of this investigation and a companion study indicate that there is little to be gained combining VHF-RaDAR volumetric returns and profiling LiDAR ranging measurements in pine forests; a LiDAR ranging system is sufficient for accurate, precise biomass estimation.     

Sensitivity of SAR data to post-fire forest regrowth in Mediterranean and boreal forests

Disturbed forests may need decades to reach a mature stage and optically-based vegetation indices are usually poorly suited for monitoring purposes due to the rapid saturation of the signal with increasing canopy cover. Spaceborne synthetic aperture radar (SAR) data provide an alternate monitoring approach since the backscattered microwave energy is sensitive to the vegetation structure. Images from two regions in Spain and Alaska were used to analyze SAR metrics (cross-polarized backscatter and co-polarized interferometric coherence) from regrowing forests previously affected by fire. TerraSAR-X X-band backscatter showed the lowest sensitivity to forest regrowth, with the average backscatter increasing by 1-2 dB between the most recent fire scar and the unburned forest. Increased sensitivity (around 3-4 dB) was observed for C-band Envisat Advanced Synthetic Aperture (ASAR) backscatter. The Advanced Land Observing Satellite (ALOS) Phased Array-type L-band Synthetic Aperture Radar (PALSAR) L-band backscatter presented the highest dynamic range from unburned to recently burned forests (approximately 8 dB). The interferometric coherence showed low sensitivity to forest regrowth at all SAR frequencies. For Mediterranean forests, five phases of forest regrowth were discerned whereas for boreal forest, up to four different regrowth phases could be discerned with L-band SAR data. In comparison, the Normalized Difference Vegetation Index (NDVI) provided reliable differentiation only for the most recent development stages. The results obtained were consistent in both environments.     

Evaluating uncertainty in mapping forest carbon with airborne LiDAR

Airborne LiDAR is increasingly used to map carbon stocks in tropical forests, but our understanding of mapping errors is constrained by the spatial resolution (i.e., plot size) used to calibrate LiDAR with field data (typically 0.1-0.36 ha). Reported LiDAR errors range from 17 to 40 Mg C ha^− 1, but should be lower at coarser resolutions because relative errors are expected to scale with (plot area)^-1/2. We tested this prediction empirically using a 50-ha plot with mapped trees, allowing an assessment of LiDAR prediction errors at multiple spatial resolutions. We found that errors scaled approximately as expected, declining by 38% (compared to 40% predicted from theory) from 0.36- to 1-ha resolution. We further reduced errors at all spatial resolutions by accounting for tree crowns that are bisected by plot edges (not typically done in forestry), and collectively show that airborne LiDAR can map carbon stocks with 10% error at 1-ha resolution — a level comparable to the use of field plots alone.     

Tropical forest backscatter anomaly evident in SeaWinds scatterometer morning overpass data during 2005 drought in Amazonia

A severe drought occurred in southwestern Amazonia in the dry season (June-September) of 2005. We analyzed 10 years (7/99-10/09) of SeaWinds active microwave Ku-band backscatter data collected over the Amazon Basin, developing monthly means and anomalies from those means in an effort to detect landscape responses to this drought. We compared these to seasonal accumulating water deficit anomalies generated using Tropical Rainfall Monitoring Mission (TRMM) precipitation data (1999-2009) and 100 mm mo^− 1 evapotranspiration demand as a water deficit threshold. There was significant interannual variability in dry-season monthly mean backscatter only for morning (c. 06:00 LST) overpass data, and little interannual variability in dry-season monthly mean backscatter for afternoon (c. 18:00 LST) overpass data. Strong negative anomalies in both morning-overpass backscatter and accumulating water deficit developed during July-October 2005, centered on the southwestern Amazon Basin, with a strong spatial correlation between morning-overpass backscatter anomalies and water deficit anomalies in September. This is the first reporting of tropical forest seasonal drought detection by active microwave scatterometry. Based on the differences between early-morning and late-afternoon backscatter variability, we hypothesize that as the drought persisted over several months, the forest canopy was increasingly unable to recover full leaf moisture content over night, resulting in anomalously low early-morning overpass backscatter.