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).     

Assessment of forest structure with airborne LiDAR and the effects of platform altitude

Airborne scanning LiDAR is a spatial technology increasingly used for forestry and environmental applications. However, the accuracy and coverage of LiDAR observations is highly dependent on both the extrinsic specifications of the LiDAR survey as well as the intrinsic effects such as the underlying forest structure. Extrinsic parameters which are set as part of the LiDAR survey include platform altitude, scan angle (half max. angle off nadir), and beam cross sectional diameter at the reflecting surface (referred to as footprint size). In this paper we investigate the effect of a number of these extrinsic parameters, including three different platform altitudes (1000, 2000, and 3000 m), two scan angles at 1000 m (10° and 15° half max. angle off nadir), and three footprint sizes (0.2, 0.4, and 0.6 m). The comparison was undertaken in eucalypt forests at three sites, varying in vegetation structure and topography within the Wedding Bells State Forest, Coffs Harbour, Australia. Results at the plot scale (40 × 90 m areas) indicate that tree heights computed from the 1000 m LiDAR data set (10° half max. angle off nadir) are well correlated with maximum plot heights (difference < 3 m) and field measured canopy volume (r² > 0.75, p < 0.001). Using normalised canopy height profiles (CHP) derived for sites, from data recorded at each altitude, we observed no significant difference between the relative distribution of LiDAR returns, indicating that platform altitude and footprint size have not had a major influence on CHP estimation. Interestingly, comparisons of first and last returns for individual pulses at increasing altitudes identified progressively fewer discrete first/last pulse combinations with more than 70% of pulses recorded as a single return at the highest altitude (3000 m). A possible hypothesis is that greater platform altitude and footprint size reduces the intensity of laser beam incident on a given surface area thus decreasing the probability of recording a last return above the noise threshold. Furthermore, tree scale analysis found a positive relationship between platform altitude and the underestimation of crown area and crown volume. The implications of this work for forest management are: (i) platform altitudes as high as 3000 m can be used to quantify the vertical distribution of phyto-elements, (ii) higher platform altitudes record a lower proportion of first/last return combinations that will further reduce the number of points available for forest structural assessment and development of digital elevation models, and (iii) for discrete LiDAR data, increasing platform altitude will record a lower frequency of returns per crown, resulting in larger underestimates of individual tree crown area and volume if standard algorithms are applied.     

Extracting LiDAR indices to characterise multilayered forest structure using mixture distribution functions

Discrete Light Detection and Ranging (LiDAR) data is used to stratify a multilayered eucalyptus forest and characterise the structure of the vertical profile. We present a methodology that may prove useful for a very broad range of forest management applications, particularly for timber inventory evaluation and forest growth modelling. In this study, we use LiDAR data to stratify a multilayered eucalyptus forest and characterise the structure of specific vegetation layers for forest hydrology research, as vegetation dynamics influence a catchment's streamflow yield. A forest stand's crown height, density, depth, and closure, influence aerodynamic properties of the forest structure and the amount of transpiring leaf area, which in turn determine evapotranspiration rates. We present a methodology that produces canopy profile indices of understorey and overstorey vegetation using mixture models with a wide range of theoretical distribution functions. Mixture models provide a mechanism to summarise complex canopy attributes into a short list of parameters that can be empirically analysed against stand characteristics.Few studies have explored theoretical distribution functions to represent the vertical profile of vegetation structure in LiDAR data. All prior studies have focused on a Weibull distribution function, which is unimodal. In a complex native forest ecosystem, the form of the distribution of LiDAR points may be highly variable between forest types and age classes. We compared 44 probability distributions within a two component mixture model to determine the most suitable bimodal distributions for representing LiDAR density estimates of Mountain Ash forests in south-eastern Australia. An elimination procedure identified eleven candidate distributions for representing the eucalyptus component of the mixture model.We demonstrate the methodology on a sample of plots to predict overstorey stand volumes and basal area, and understorey basal area of 18-, 37-, and 70-year old Mountain Ash forest with variable density classes. The 70-year old forest has been subjected to a range of treatments including: thinning of the eucalyptus layer with two distinct retention rates, removal of the understorey, and clear felling of patches that have 37 year old regenerating forest. We demonstrate that the methodology has clear potential, as observed versus predicted values of eucalyptus basal area and stand volume were highly correlated, with bootstrap based r² ranging from 0.61 to 0.89 and 0.67 to 0.88 respectively. Non-eucalyptus basal area r² ranged from 0.5 to 0.91.     

A temporal analysis of urban forest carbon storage using remote sensing

Quantifying the carbon storage, distribution, and change of urban trees is vital to understanding the role of vegetation in the urban environment. At present, this is mostly achieved through ground study. This paper presents a method based on the satellite image time series, which can save time and money and greatly speed the process of urban forest carbon storage mapping, and possibly of regional forest mapping. Satellite imagery collected in different decades was used to develop a regression equation to predict the urban forest carbon storage from the Normalized Difference Vegetation Index (NDVI) computed from a time sequence (1985-1999) of Landsat image data. This regression was developed from the 1999 field-based model estimates of carbon storage in Syracuse, NY. The total carbon storage estimates based on the NDVI data agree closely with the field-based model estimates. Changes in total carbon storage by trees in Syracuse were estimated using the image data from 1985, 1992, and 1999. Radiometric correction was accomplished by normalizing the imagery to the 1999 image data. After the radiometric image correction, the carbon storage by urban trees in Syracuse was estimated to be 146,800 tons, 149,430 tons, and 148,660 tons of carbon for 1985, 1992, and 1999, respectively. The results demonstrate the rapid and cost-effective capability of remote sensing-based quantitative change detection in monitoring the carbon storage change and the impact of urban forest management over wide areas.     

Lidar remote sensing of forest biomass: A scale-invariant estimation approach using airborne lasers

Researchers in lidar (Light Detection And Ranging) strive to search for the most appropriate laser-based metrics as predictors in regression models for estimating forest structural variables. Many previously developed models are scale-dependent that need to be fitted and then applied both at the same scale or pixel size. The objective of this paper is to develop methods for scale-invariant estimation of forest biomass using lidar data. We proposed two scale-invariant models for biomass: a linear functional model and an equivalent nonlinear model that use lidar-derived canopy height distributions (CHD) and canopy height quantile functions (CHQ) as predictors, respectively. The two models are called functional regression models because the predictors CHD and CHQ are themselves functions or functional data. The model formulation was justified mathematically under moderate assumptions. We also created a fine-resolution biomass map by mapping individual tree component biomass in a temperate forest of eastern Texas with a lidar tree-delineation approach. The map was used as reference data to synthesize training and test datasets at multiple scales for validating the two scale-invariant models. Results suggest that the models can accurately predict biomass and yield consistent predictive performances across a variety of scales with an R² ranging from 0.80 to 0.95 (RMSE: from 14. 3 Mg/ha to 33.7 Mg/ha) among all the fitted models. Results also show that a training data size of around 50 plots or less was enough to guarantee a good fitting of the linear functional model. Our findings demonstrate the effectiveness of CHD and CHQ as lidar metrics for estimating biomass as well as the capability of lidar for mapping biomass at a range of scales. The functional regression models of this study are useful for lidar-based forest inventory tasks where the analysis units vary in size and shape. They also hold promise for estimating other forest characteristics such as below-ground biomass, timber volume, crown fuel weight, and Leaf Area Index.     

Assessing forest structural and physiological information content of multi-spectral LiDAR waveforms by radiative transfer modelling

The concept for a multi-spectral, full-waveform canopy LiDAR instrument was tested by simulating return waveforms using a model providing ecological sound tree structure (TREEGROW) and a model of leaf optical properties (PROSPECT). The proposed instrument will take measurements at four different wavelengths, which were chosen according to physiological processes altering leaf reflectance and transmittance. The modelling was used to assess both the structural and physiological information content such an instrument could provide, especially whether the normally structure-dominated return waveform would pick up small changes in reflectance at the leaf level. Multi-spectral waveforms were simulated for models of single Scots pine trees of different ages and at different stages of the growing season, including chlorophyll concentration induced changes in leaf optical properties. It was shown that the LiDAR waveforms would not only capture the tree height information, but would also pick up the seasonal and vertical variation of NDVI computed from two of the four MSCL wavelengths inside the tree canopy. The instrument concept was further tested in a simulation of a virtual forest stand constructed of 74 trees of different ages according to measurements taken on a field site being 20 by 20 meter in size. A total of 1521 NDVI profiles were computed and mean NDVI corrected backscatter was compared to the actual canopy profile of the virtual stand. The profiles picked up the seasonal variation of chlorophyll within the canopy, while the return of ground remained unchanged from June to September. Thus, it was shown that a MSCL instrument would be able to separately pick up the physiology of canopy and understorey and/or soil. It was found that occlusion would mask the lower parts of the canopy volume within the stand and the seasonal variation of this occlusion effect was quantified, being larger in September, when the absorption of canopy elements is higher. In addition, it could be demonstrated that a new multi-wavelength LiDAR predictor variable was able to significantly improve the retrieval accuracy of photosynthetically active biomass opposed to using a single-wavelength LiDAR alone.     

Strengths and limitations of assessing forest density and spatial configuration with aerial LiDAR

Changes in the structural state of forests of the semi-arid U.S.A., such as an increase in tree density, are widely believed to be leading to an ecological crisis, but accurate methods of quantifying forest density and configuration are lacking at landscape scales. An individual tree canopy (ITC) method based on aerial LiDAR has been developed to assess forest structure by estimating the density and spatial configuration of trees in four different height classes. The method has been tested against field measured forest inventory data from two geographically distinct forests with independent LiDAR acquisitions. The results show two distinct patterns: accurate, unbiased density estimates for trees taller than 20 m, and underestimation of density in trees less than 20 m tall. The underestimation of smaller trees is suggested to be a limitation of LiDAR remote sensing. Ecological applications of the method are demonstrated through landscape metrics analysis of density and configuration rasters.     

Remote sensing of species mixtures in conifer plantations using LiDAR height and intensity data

In even-aged, single species conifer plantations LiDAR height data can be modelled to provide accurate estimates of tree height and volume. However, it is apparent that growth models developed for single species stands are not directly transferable to a more general situation of mixed species plantations. This paper evaluates the ability of small footprint, dual-return, pulsed airborne LiDAR data to estimate the proportion of the productive species when mixed with a nurse crop in closed canopy plantations. A study area located in Galloway Forest District in Scotland is used as an example of Lodgepole pine and Sitka spruce mixed plantation; this area contains good examples of a wide range of pure and mixed species plantation types. Three species groups are studied: areas of pure Sitka spruce, areas of pure Lodgepole pine and areas where the two species have been planted together. Two approaches are assessed for detection of plantation mixtures: the first uses LiDAR intensity data to separate spruce and pine species and the second uses LiDAR-derived canopy density measures, coefficient of variation, skewness, percent of ground returns (which provides a measure of canopy openness) and the mean canopy height, which enables areas with height variations to be identified. From analysis of LiDAR data extracted from 54 study plots using logistic regression, the coefficient of variation and LiDAR intensity data provide the most useful predictors of the proportion of spruce in a pine/spruce mixture with coefficients of determination (R²) of 0.914 and 0.930 respectively. The method could be developed as a mapping tool, which in combination with existing inventory data should help to improve timber volume forecasting for mixed species even-aged plantations.     

Mapping tropical forest structure in southeastern Madagascar using remote sensing and artificial neural networks

Tropical forest condition has important implications for biodiversity, climate change and human needs. Structural features of forests can serve as useful indicators of forest condition and have the potential to be assessed with remotely sensed imagery, which can provide quantitative information on forest ecosystems at high temporal and spatial resolutions. Herein, we investigate the utility of remote sensing for assessing, predicting and mapping two important forest structural features, stem density and basal area, in tropical, littoral forests in southeastern Madagascar. We analysed the relationships of basal area and stem density measurements to the Normalised Difference Vegetation Index (NDVI) and radiance measurements in bands 3, 4, 5 and 7 from the Landsat Enhanced Thematic Mapper Plus (ETM+). Strong relationships were identified among all of the individual bands and field based measurements of basal area (p<0.01) while there were weak and insignificant relationships among spectral response and stem density measurements. NDVI was not significantly correlated with basal area but was strongly and significantly correlated with stem density (r=−0.69, p<0.01) when using a subset of the data, which represented extreme values. We used an artificial neural network (ANN) to predict basal area from radiance values in bands 3, 4, 5 and 7 and to produce a predictive map of basal area for the entire forest landscape. The ANNs produced strong and significant relationships between predicted and actual measures of basal area using a jackknife method (r=0.79, p<0.01) and when using a larger data set (r=0.82, p<0.01). The map of predicted basal area produced by the ANN was assessed in relation to a pre-existing map of forest condition derived from a semi-quantitative field assessment. The predictive map of basal area provided finer detail on stand structural heterogeneity, captured known climatic influences on forest structure and displayed trends of basal area associated with degree of human accessibility. These findings demonstrate the utility of ANNs for integrating satellite data from the Landsat ETM+ spectral bands 3, 4, 5 and 7 with limited field survey data to assess patterns in basal area at the landscape scale.     

Development of a simulation model to predict LiDAR interception in forested environments

Airborne scanning LiDAR systems are used to predict a range of forest attributes. However, the accuracy with which this can be achieved is highly dependent on the sensor configuration and the structural characteristics of the forest examined. As a result, there is a need to understand laser light interactions with forest canopies so that LiDAR sensor configurations can be optimised to assess particular forest types. Such optimisation will not only ensure the targeted forest attributes can be accurately and consistently quantified, but may also minimise the cost of data acquisition and indicate when a survey configuration will not deliver information needs.In this paper, we detail the development and application of a model to simulate laser interactions within forested environments. The developed model, known as the LiDAR Interception and Tree Environment (LITE) model, utilises a range of structural configurations to simulate trees with variable heights, crown dimensions and foliage clumping. We developed and validated the LITE model using field data obtained from three forested sites covering a range of structural classes. Model simulations were then compared to coincident airborne LiDAR data collected over the same sites. Results indicate that the LITE model can be used to produce comparable estimates of maximum height of trees within plots (differences < 2.42 m), mean heights of first return data (differences < 2.27 m), and canopy height percentiles (r² = 0.94, p < 0.001) when compared to airborne LiDAR data. In addition, the distribution of airborne LiDAR hits through the canopy profile was closely matched by model predictions across the range of sites. Importantly, this demonstrates that the structural differences between forest stands can be characterised by LITE. Models that are capable of interpreting the response of small-footprint LiDAR waveforms can facilitate algorithm development, the generation of corrections for actual LiDAR data, and the optimisation of sensor configurations for differing forest types, benefiting a range of experimental and commercial LiDAR applications. As a result, we also performed a scenario analysis to demonstrate how differences in forest structure, terrain, and sensor configuration can influence the interception of LiDAR beams.     

Three-dimensional canopy fuel loading predicted using upward and downward sensing LiDAR systems

We calibrated upward sensing profiling and downward sensing scanning LiDAR systems to estimates of canopy fuel loading developed from field plots and allometric equations, and then used the LiDAR datasets to predict canopy bulk density (CBD) and crown fuel weight (CFW) in wildfire prone stands in the New Jersey Pinelands. LiDAR-derived height profiles were also generated in 1-m layers, and regressed on CBD estimates calculated for 1-m layers from field plots to predict three-dimensional canopy fuel loading. We then produced maps of canopy fuel metrics for three 9 km² forested areas in the Pinelands.Correlations for standard LiDAR-derived parameters between the two LiDAR systems were all highly significant, with correlation coefficients ranging between 0.82 and 0.98. Stepwise linear regression models developed from the profiling LiDAR data predicted maximum CBD and CFW (r² = 0.94 and 0.92) better than those developed from the scanning LiDAR data (r² = 0.83 and 0.71, respectively). A single regression for the prediction of CBD at all canopy layers had r² values of 0.93 and 0.82 for the profiling and scanning datasets, respectively. Individual bin regressions for predicting CBD at each canopy height layer were also highly significant at most canopy heights, with r² values for each layer ranging between 0.36 and 0.89, and 0.44 and 0.99 for the profiling and scanning datasets, respectively. Relationships were poorest mid-canopy, where highest average values and highest variability in fuel loading occurred. Fit of data to Gaussian distributions of canopy height profiles facilitated a simpler expression of these parameters for analysis and mapping purposes, with overall r² values of 0.86 and 0.92 for the profiling and scanning LiDAR datasets, respectively. Our research demonstrates that LiDAR data can be used to generate accurate, three-dimensional representations of canopy structure and fuel loading at high spatial resolution by linking 1-m return height profiles to biometric estimates from field plots.     

Monitoring of cropland practices for carbon sequestration purposes in north central Montana by Landsat remote sensing

We used an object-oriented approach in conjunction with the Random Forest algorithm to classify agricultural practices, including tillage (till or no-till (NT)), crop intensity, and grassland-based conservation reserve (CR). The object-oriented approach allowed for per-field classifications and the incorporation of contextual elements in addition to spectral features. Random Forest is a classification tree-based advanced classifier that avoids data over-fitting associated with many tree-based models and incorporates an unbiased internal classification accuracy assessment. Landsat satellite imagery was chosen for its continuous coverage, cost effectiveness, and image accessibility. Classification results for 2007 included producer's accuracies of 91% for NT and 31% for tillage when applying Random Forest to image objects generated from a May Landsat image. Low classification accuracies likely were attributed to the misclassification of conservation-based tillage practices as NT. Results showed that the binary separation of tillage from NT management is likely not appropriate due to surface spectral and textural similarities between NT and conservation-type tillage practices. Crop and CR lands resulted in producer's accuracies of 100% and 90%, respectively. Crop and fallow producer's accuracies were 95% and 82% in the 2007 classification, despite post-senesced vegetation; misclassification within the fallow class was attributed to pixel-mixing problems in areas of narrow (< 100 m) strip management. A between-date normalized difference vegetation index approach was successfully used to detect areas having “changed” in vegetation status between the 2007 and prior image dates; classified “changed” objects were then merged with “unchanged” objects to produce crop status maps. Field crop intensity was then determined from the multi-year analysis of generated crop status maps.     

Fusion of LiDAR and imagery for estimating forest canopy fuels

Due to increased fuel loading as a result of fire suppression, land managers in the American west are in need of precise information about the fuels they manage, including canopy fuels. Canopy fuel metrics such as canopy height (CH), canopy base height (CBH), canopy bulk density (CBD) and available canopy fuel (ACF) are specific inputs for wildfire behavior models such as FARSITE and emission models such as FOFEM. With finer spatial resolution data, accurate quantification of these metrics with detailed spatial heterogeneity can be accomplished. Light Detection and Ranging (LiDAR) and color near-infrared imagery are active and passive systems, respectively, that have been utilized for measuring a range of forest structure characteristics at high resolution. The objective of this research was to determine which remote sensing dataset can estimate canopy fuels more accurately and whether a fusion of these datasets produces more accurate estimates. Regression models were developed for ponderosa pine (Pinus ponderosa) stand representative of eastern Washington State using field data collected in the Ahtanum State Forest and metrics derived from LiDAR and imagery. Strong relationships were found with LiDAR alone and LiDAR was found to increase canopy fuel accuracy compared to imagery. Fusing LiDAR with imagery and/or LiDAR intensity led to small increases in estimation accuracy over LiDAR alone. By improving the ability to estimate canopy fuels at higher spatial resolutions, spatially explicit fuel layers can be created and used in wildfire behavior and smoke emission models leading to more accurate estimations of crown fire risk and smoke related emissions.     

Modeling forest songbird species richness using LiDAR-derived measures of forest structure

Conservation of biodiversity requires information at many spatial scales in order to detect and preserve habitat for many species, often simultaneously. Vegetation structure information is particularly important for avian habitat models and has largely been unavailable for large areas at the desired resolution. Airborne LiDAR, with its combination of relatively broad coverage and fine resolution provides existing new opportunities to map vegetation structure and hence avian habitat. Our goal was to model the richness of forest songbirds using forest structure information obtained from LiDAR data. In deciduous forests of southern Wisconsin, USA, we used discrete-return airborne LiDAR to derive forest structure metrics related to the height and density of vegetation returns, as well as composite variables that captured major forest structural elements. We conducted point counts to determine total forest songbird richness and the richness of foraging, nesting, and forest edge-related habitat guilds. A suite of 35 LiDAR variables were used to model bird species richness using best-subsets regression and we used hierarchical partitioning analysis to quantify the explanatory power of each variable in the multivariate models. Songbird species richness was correlated most strongly with LiDAR variables related to canopy and midstory height and midstory density (R² = 0.204, p < 0.001). Richness of species that nest in the midstory was best explained by canopy height variables (R² = 0.197, p < 0.001). Species that forage on the ground responded to mean canopy height and the height of the lower canopy (R² = 0.149, p < 0.005) while aerial foragers had higher richness where the canopy was tall and dense and the midstory more sparse (R² = 0.216, p < 0.001). Richness of edge-preferring species was greater where there were fewer vegetation returns but higher density in the understory (R² = 0.153, p < 0.005). Forest interior specialists responded positively to a tall canopy, developed midstory, and a higher proportion of vegetation returns (R² = 0.195, p < 0.001). LiDAR forest structure metrics explained between 15 and 20% of the variability in richness within deciduous forest songbird communities. This variability was associated with vertical structure alone and shows how LiDAR can provide a source of complementary predictive data that can be incorporated in models of wildlife habitat associations across broad geographical extents.     

Quantifying the spatial properties of forest canopy gaps using LiDAR imagery and GIS

The spatial properties of gaps have an important influence upon the regeneration dynamics and species composition of forests. However, such properties can be difficult to quantify over large spatial areas using field measurements. This research considers how we conceptualize and define forest canopy gaps from a remote sensing point of view and highlights the inadequacies of passive optical remotely sensed data for delineating gaps. The study employs the analytical functions of a geographical information system to extract gap spatial characteristics from imagery acquired by an active remote sensing device, an airborne light detection and ranging instrument (LiDAR). These techniques were applied to an area of semi-natural broadleaved deciduous forest, in order to map gap size, shape complexity, vegetation height diversity and gap connectivity. A vegetation cover map derived from imagery from an airborne multispectral scanner was used in combination with the LiDAR data to characterize the dominant vegetation types within gaps. Although the quantification of these gap characteristics alone is insufficient to provide conclusive evidence on specific processes, the paper demonstrates how such information can be indicative of the general status of a forest and can provide new perspectives and possibilities or further ecological research and forest monitoring activities.     

Tropical forest typo mapping and monitoring using remote sensing

Abstract

The Andaman Nicobar group of islands in the Andaman Sea are blessed with luxuriant tropical vegetation. During recent years, vegetation in these islands has been under tremendous pressure due to increased settlement and commercial extraction, Landsat TM data of the Baratang forest division of the Andaman group of islands has been used to prepare forest type maps using visual and digital methods. Digital enhancement techniques have been evaluated to discriminate forest types in a typical environmental set-up. The mapping techniques have been compared with respect to classification and accuracy levels. Finally, the land transformation in the forest division has been studied using past aerial photographs. The study highlights the appropriate methodology required to map forest types.     

Mapping U.S. forest biomass using nationwide forest inventory data and moderate resolution information

A spatially explicit dataset of aboveground live forest biomass was made from ground measured inventory plots for the conterminous U.S., Alaska and Puerto Rico. The plot data are from the USDA Forest Service Forest Inventory and Analysis (FIA) program. To scale these plot data to maps, we developed models relating field-measured response variables to plot attributes serving as the predictor variables. The plot attributes came from intersecting plot coordinates with geospatial datasets. Consequently, these models serve as mapping models. The geospatial predictor variables included Moderate Resolution Imaging Spectrometer (MODIS)-derived image composites and percent tree cover; land cover proportions and other data from the National Land Cover Dataset (NLCD); topographic variables; monthly and annual climate parameters; and other ancillary variables. We segmented the mapping models for the U.S. into 65 ecologically similar mapping zones, plus Alaska and Puerto Rico. First, we developed a forest mask by modeling the forest vs. nonforest assignment of field plots as functions of the predictor layers using classification trees in See5©. Secondly, forest biomass models were built within the predicted forest areas using tree-based algorithms in Cubist©. To validate the models, we compared field-measured with model-predicted forest/nonforest classification and biomass from an independent test set, randomly selected from available plot data for each mapping zone. The estimated proportion of correctly classified pixels for the forest mask ranged from 0.79 in Puerto Rico to 0.94 in Alaska. For biomass, model correlation coefficients ranged from a high of 0.73 in the Pacific Northwest, to a low of 0.31 in the Southern region. There was a tendency in all regions for these models to over-predict areas of small biomass and under-predict areas of large biomass, not capturing the full range in variability. Map-based estimates of forest area and forest biomass compared well with traditional plot-based estimates for individual states and for four scales of spatial aggregation. Variable importance analyses revealed that MODIS-derived information could contribute more predictive power than other classes of information when used in isolation. However, the true contribution of each variable is confounded by high correlations. Consequently, excluding any one class of variables resulted in only small effects on overall map accuracy. An estimate of total C pools in live forest biomass of U.S. forests, derived from the nationwide biomass map, also compared well with previously published estimates.     

Combining remote sensing and ancillary data to monitor the gross productivity of water-limited forest ecosystems

This paper describes the development and testing of a procedure which combines remotely sensed and ancillary data to monitor forest productivity in Italy. The procedure is based on a straightforward parametric model (C-Fix) that uses the relationship between the fraction of photosynthetically active radiation absorbed by plant canopies (fAPAR) and relevant gross primary productivity (GPP). Estimates of forest fAPAR are derived from Spot-VGT NDVI images and are combined with spatially consistent data layers obtained by the elaboration of ground meteorological measurements. The original version of C-Fix is first applied to estimate monthly GPP of Italian forests during eight years (1999-2006). Next, a modification of the model is proposed in order to simulate the short-term effect of summer water stress more efficiently. The accuracy of the original and modified C-Fix versions is evaluated by comparison with GPP data taken at eight Italian eddy covariance flux tower sites. The experimental results confirm the capacity of C-Fix to monitor national forest GPP patterns and indicate the utility of considering the short-term effect of water stress during Mediterranean dry months.     

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.