A Structural Classification of Australian Vegetation Using ICESat/GLAS, ALOS PALSAR and Landsat Sensor Data

Australia has historically used structural descriptors of height and cover to characterize, differentiate, and map the distribution of woody vegetation across the continent but no national satellite-based structural classification has been available. In this study, we present a new 30-m spatial resolution reference map of Australian forest and woodland structure (height and cover), with this generated by integrating Landsat Thematic Mapper (TM) and Enhanced TM, Advanced Land Observing Satellite (ALOS) Phased Arrayed L-band Synthetic Aperture Radar (PALSAR) and Ice, Cloud, and land Elevation (ICESat),and Geoscience Laser Altimeter System (GLAS) data. ALOS PALSAR and Landsat-derived Foliage Projective Cover (FPC) were used to segment and classify the Australian landscape. Then, from intersecting ICESat waveform data, vertical foliage profiles and height metrics (e.g., 95% percentile height, mean height and the height to maximum vegetation density) were extracted for each of the classes generated. Within each class, and for selected areas, the variability in ICESat profiles was found to be similar with differences between segments of the same class attributed largely to clearance or disturbance events. ICESat metrics and profiles were then assigned to all remaining segments across Australia with the same class allocation. Validation against airborne LiDAR for a range of forest structural types indicated a high degree of correspondence in estimated height measures. On this basis, a map of vegetation height was generated at a national level and was combined with estimates of cover to produce a revised structural classification based on the scheme of the Australian National Vegetation Information System (NVIS). The benefits of integrating the three datasets for segmenting and classifying the landscape and retrieving biophysical attributes was highlighted with this leading the way for future mapping using ALOS-2 PALSAR-2, Landsat/Sentinel-2, Global Ecosystem Dynamics Investigation (GEDI), and ICESat-2 LiDAR data. The ability to map across large areas provides considerable benefits for quantifying carbon dynamics and informing on biodiversity metrics

An Approach to Mapping Forest Growth Stages in Queensland, Australia through Integration of ALOS PALSAR and Landsat Sensor Data

Whilst extensive clearance of forests in the eastern Australian Brigalow Belt Bioregion (BBB) has occurred since European settlement, appropriate management of those that are regenerating can facilitate restoration of biomass (carbon) and biodiversity to levels typical of relatively undisturbed or remnant formations. However, maps of forests are different stages of regeneration are needed to facilitate restoration planning, including prevention of further re-clearing. Focusing on the Tara Downs subregion of the BBB and on forests with brigalow (Acacia harpophylla) as a component, this research establishes a method for differentiating and mapping early, intermediate and remnant growth stages from Japan Aerospace Exploration Agency (JAXA) Advanced Land Observing Satellite (ALOS) Phased-Array L-band Synthetic Aperture Radar (PALSAR) Fine Beam Dual (FBD) L-band HH- and HV-polarisation backscatter and Landsat-derived Foliage Projective Cover (FPC). Using inventory data collected from 74 plots, located in the Tara Downs subregion, forests were assigned to one of three regrowth stages based on their height and cover relative to that of undisturbed stands. The image data were then segmented into objects with each assigned to a growth stage by comparing the distributions of L-band HV and HH polarisation backscatter and FPC to that of reference distributions using a z-test. Comparison with independent assessments of growth stage, based on time-series analysis of aerial photography and SPOT images, established an overall accuracy of > 70%, with this increasing to 90% when intermediate regrowth was excluded and only early-stage regrowth and remnant classes were considered. The proposed method can be adapted to respond to amendments to user-definitions of growth stage and, as regional mosaics of ALOS PALSAR and Landsat FPC are available for Queensland, has application across the state

Detecting Change in Forest Structure with Simulated GEDI Lidar Waveforms: A Case Study of the Hemlock Woolly Adelgid (HWA; Adelges tsugae) Infestation

The hemlock woolly adelgid (HWA; Adelges tsugae) is an invasive insect infestation that is spreading into the forests of the northeastern United States, driven by the warmer winter temperatures associated with climate change. The initial stages of this disturbance are difficult to detect with passive optical remote sensing, since the insect often causes its host species, eastern hemlock trees (Tsuga canadensis), to defoliate in the midstory and understory before showing impacts in the overstory. New active remote sensing technologies-such as the recently launched NASA Global Ecosystem Dynamics Investigation (GEDI) spaceborne lidar-can address this limitation by penetrating canopy gaps and recording lower canopy structural changes. This study explores new opportunities for monitoring the HWA infestation with airborne lidar scanning (ALS) and GEDI spaceborne lidar data. GEDI waveforms were simulated using airborne lidar datasets from an HWA-infested forest plot at the Harvard Forest ForestGEO site in central Massachusetts. Two airborne lidar instruments, the NASA G-LiHT and the NEON AOP, overflew the site in 2012 and 2016. GEDI waveforms were simulated from each airborne lidar dataset, and the change in waveform metrics from 2012 to 2016 was compared to field-derived hemlock mortality at the ForestGEO site. Hemlock plots were shown to be undergoing dynamic changes as a result of the HWA infestation, losing substantial plant area in the middle canopy, while still growing in the upper canopy. Changes in midstory plant area (PAI 11-12 m above ground) and overall canopy permeability (indicated by RH10) accounted for 60% of the variation in hemlock mortality in a logistic regression model. The robustness of these structure-condition relationships held even when simulated waveforms were treated as real GEDI data with added noise and sparse spatial coverage. These results show promise for future disturbance monitoring studies with ALS and GEDI lidar data

Biomass estimation from simulated GEDI, ICESat-2 and NISAR across environmental gradients in Sonoma County, California

Estimates of the magnitude and distribution of aboveground carbon in Earth's forests remain uncertain, yet knowledge of forest carbon content at a global scale is critical for forest management in support of climate mitigation. In light of this knowledge gap, several upcoming spaceborne missions aim to map forest aboveground biomass, and many new biomass products are expected from these datasets. As these new missions host different technologies, each with relative strengths and weaknesses for biomass retrieval, as well as different spatial resolutions, consistently comparing or combining biomass estimates from these new datasets will be challenging. This paper presents a demonstration of an inter-comparison of biomass estimates from simulations of three NASA missions (GEDI, ICESat-2 and NISAR) over Sonoma county in California, USA. We use a high resolution, locally calibrated airborne lidar map as our reference dataset, and emphasize the importance of considering uncertainties in both reference maps and spaceborne estimates when conducting biomass product validation. GEDI and ICESat-2 were simulated from airborne lidar point clouds, while UAVSAR's L-band backscatter was used as a proxy for NISAR. To estimate biomass for the lidar missions we used GEDI's footprint-level biomass algorithms, and also adapted these for application to ICESat-2. For UAVSAR, we developed a locally trained biomass model, calibrated against the ALS reference map. Each mission simulation was evaluated in comparison to the local reference map at its native product resolution (25 m, 100 m transect, and 1 ha) yielding RMSEs of 57%, 75%, and 89% for GEDI, NISAR, and ICESat-2 respectively. RMSE values increased for GEDI's power beam during simulated daytime conditions (64%), coverage beam during nighttime conditions (72%), and coverage beam daytime conditions (87%). We also test the application of GEDI's biomass modeling framework for estimation of biomass from ICESat-2, and find that ICESat-2 yields reasonable biomass estimates, particularly in relatively short, open canopies. Results suggest that while all three missions will produce datasets useful for biomass mapping, tall, dense canopies such as those found in Sonoma County present the greatest challenges for all three missions, while steep slopes also prove challenging for single-date SAR-based biomass retrievals. Our methods provide guidance for the inter-comparison and validation of spaceborne biomass estimates through the use of airborne lidar reference maps, and could be repeated with on-orbit estimates in any area with high quality field plot and ALS data. These methods allow for regional interpretations and filtering of multi-mission biomass estimates toward improved wall-to-wall biomass maps through data fusion.</p

Exploring the relation between remotely sensed vertical canopy structure and tree species diversity in Gabon

Mapping tree species diversity is increasingly important in the face of environmental change and biodiversity conservation. We explore a potential way of mapping this diversity by relating forest structure to tree species diversity in Gabon. First, we test the relation between canopy height, as a proxy for niche volume, and tree species diversity. Then, we test the relation between vertical canopy structure, as a proxy for vertical niche occupation, and tree species diversity. We use large footprint full-waveform airborne lidar data collected across four study sites in Gabon (Lopé, Mabounié, Mondah, and Rabi) in combination with in situ estimates of species richness (S) and Shannon diversity (H′). Linear models using canopy height explained 44% and 43% of the variation in S and H′ at the 0.25 ha resolution. Linear models using canopy height and the plant area volume density profile explained 71% of this variation. We demonstrate applications of these models by mapping S and H′ in Mondah using a simulated GEDI-TanDEM-X fusion height product, across the four sites using wall-to-wall airborne lidar data products, and across and between the study sites using ICESat lidar waveforms. The modeling results are encouraging in the context of developing pan-tropical structure diversity models applicable to data from current and upcoming spaceborne remote sensing missions

Aboveground biomass density models for NASA's Global Ecosystem Dynamics Investigation (GEDI) lidar mission

NASA's Global Ecosystem Dynamics Investigation (GEDI) is collecting spaceborne full waveform lidar data with a primary science goal of producing accurate estimates of forest aboveground biomass density (AGBD). This paper presents the development of the models used to create GEDI's footprint-level (similar to 25 m) AGBD (GEDI04_A) product, including a description of the datasets used and the procedure for final model selection. The data used to fit our models are from a compilation of globally distributed spatially and temporally coincident field and airborne lidar datasets, whereby we simulated GEDI-like waveforms from airborne lidar to build a calibration database. We used this database to expand the geographic extent of past waveform lidar studies, and divided the globe into four broad strata by Plant Functional Type (PFT) and six geographic regions. GEDI's waveform-to-biomass models take the form of parametric Ordinary Least Squares (OLS) models with simulated Relative Height (RH) metrics as predictor variables. From an exhaustive set of candidate models, we selected the best input predictor variables, and data transformations for each geographic stratum in the GEDI domain to produce a set of comprehensive predictive footprint-level models. We found that model selection frequently favored combinations of RH metrics at the 98th, 90th, 50th, and 10th height above ground-level percentiles (RH98, RH90, RH50, and RH10, respectively), but that inclusion of lower RH metrics (e.g. RH10) did not markedly improve model performance. Second, forced inclusion of RH98 in all models was important and did not degrade model performance, and the best performing models were parsimonious, typically having only 1-3 predictors. Third, stratification by geographic domain (PFT, geographic region) improved model performance in comparison to global models without stratification. Fourth, for the vast majority of strata, the best performing models were fit using square root transformation of field AGBD and/or height metrics. There was considerable variability in model performance across geographic strata, and areas with sparse training data and/or high AGBD values had the poorest performance. These models are used to produce global predictions of AGBD, but will be improved in the future as more and better training data become available

Mapping riparian condition indicators in a sub-tropical savanna environment from discrete return LiDAR data using object-based image analysis

Mapping, monitoring and managing the environmental condition of riparian zones are major focus areas for local and state governments in Australia. New remotely sensed data techniques that can provide the required mapping accuracies, complete spatial coverage and processing and mapping transferability are currently being developed for use over large spatial extents. The research objective was to develop and apply an approach for mapping riparian condition indicators using object-based image analysis of airborne Light Detection and Ranging (LiDAR) data. The indicators assessed were: streambed width; riparian zone width; plant projective cover (PPC); longitudinal continuity; coverage of large trees; vegetation overhang; and stream bank stability. LiDAR data were captured on 15 July 2007 for two 5 km stretches along Mimosa Creek in Central Queensland, Australia. Field measurements of riparian vegetation structural and land form parameters were obtained between 28 May and 5 June 2007. Object-based approaches were developed for mapping each riparian condition indicator from the LiDAR data. The validation and empirical modelling results showed that the object-based approach could be used to accurately map the riparian condition indicators (R-2 = 0.99 for streambed width, R-2 = 0.82 for riparian zone width, R-2 = 0.89 for PPC, R-2 = 0.40 for bank stability). These research findings will be used in a 26,000 km mapping project assessing riparian vegetation and physical form indicators from LiDAR data in Victoria, Australia. (C) 2010 Elsevier Ltd. All rights reserved

Comparison of Discrete Return and Waveform Airborne Lidar Derived Estimates of Fractional Cover in an Australian Savanna

The advance of commercial airborne lidar systems from discrete-return to waveform recording instruments has made repeatable estimates of biophysical variables from these different methods questionable. Using an experimental airborne waveform lidar dataset acquired in an Australian savanna, this study presents a method for the derivation of canopy/ground backscatter coefficients from waveform lidar and a comparison of discrete return and waveform approaches to the estimation of fractional cover. Despite limited validation, the results indicate that waveform estimates of fractional cover can provide consistently higher accuracy than discrete return estimates under varying survey properties. Ongoing work using raw waveform data across larger areas and 3D radiative transfer simulations aims to develop a quantitative understanding of the impact of disparate sensor and survey properties on the detection of change in vegetation structure using commercial lidar instruments.Non peer reviewe