Regional Scale Dryland Vegetation Classification with an Integrated Lidar-Hyperspectral Approach

The sparse canopy cover and large contribution of bright background soil, along with the heterogeneous vegetation types in close proximity, are common challenges for mapping dryland vegetation with remote sensing. Consequently, the results of a single classification algorithm or one type of sensor to characterize dryland vegetation typically show low accuracy and lack robustness. In our study, we improved classification accuracy in a semi-arid ecosystem based on the use of vegetation optical (hyperspectral) and structural (lidar) information combined with the environmental characteristics of the landscape. To accomplish this goal, we used both spectral angle mapper (SAM) and multiple endmember spectral mixture analysis (MESMA) for optical vegetation classification. Lidar-derived maximum vegetation height and delineated riparian zones were then used to modify the optical classification. Incorporating the lidar information into the classification scheme increased the overall accuracy from 60% to 89%. Canopy structure can have a strong influence on spectral variability and the lidar provided complementary information for SAM’s sensitivity to shape but not magnitude of the spectra. Similar approaches to map large regions of drylands with low uncertainty may be readily implemented with unmixing algorithms applied to upcoming space-based imaging spectroscopy and lidar. This study advances our understanding of the nuances associated with mapping xeric and mesic regions, and highlights the importance of incorporating complementary algorithms and sensors to accurately characterize the heterogeneity of dryland ecosystems

Regional scale dryland vegetation classification with an integrated lidar-hyperspectral approach

The sparse canopy cover and large contribution of bright background soil, along with the heterogeneous vegetation types in close proximity, are common challenges for mapping dryland vegetation with remote sensing. Consequently, the results of a single classification algorithm or one type of sensor to characterize dryland vegetation typically show low accuracy and lack robustness. In our study, we improved classification accuracy in a semi-arid ecosystem based on the use of vegetation optical (hyperspectral) and structural (lidar) information combined with the environmental characteristics of the landscape. To accomplish this goal, we used both spectral angle mapper (SAM) and multiple endmember spectral mixture analysis (MESMA) for optical vegetation classification. Lidar-derived maximum vegetation height and delineated riparian zones were then used to modify the optical classification. Incorporating the lidar information into the classification scheme increased the overall accuracy from 60% to 89%. Canopy structure can have a strong influence on spectral variability and the lidar provided complementary information for SAM’s sensitivity to shape but not magnitude of the spectra. Similar approaches to map large regions of drylands with low uncertainty may be readily implemented with unmixing algorithms applied to upcoming space-based imaging spectroscopy and lidar. This study advances our understanding of the nuances associated with mapping xeric and mesic regions, and highlights the importance of incorporating complementary algorithms and sensors to accurately characterize the heterogeneity of dryland ecosystems

Characterizing Dryland Ecosystems Using Remote Sensing and Dynamic Global Vegetation Modeling

Drylands include all terrestrial regions where the production of crops, forage, wood and other ecosystem services are limited by water. These ecosystems cover approximately 40% of the earth terrestrial surface and accommodate more than 2 billion people (Millennium Ecosystem Assessment, 2005). Moreover, the interannual variability of the global carbon budget is strongly regulated by vegetation dynamics in drylands. Understanding the dynamics of such ecosystems is significant for assessing the potential for and impacts of natural or anthropogenic disturbances and mitigation planning, and a necessary step toward enhancing the economic and social well-being of dryland communities in a sustainable manner (Global Drylands: A UN system-wide response, 2011). In this research, a combination of remote sensing, field data collection, and ecosystem modeling were used to establish an integrated framework for semi-arid ecosystems dynamics monitoring. Foliar nitrogen (N) plays an important role in vegetation processes such as photosynthesis and there is wide interest in retrieving this variable from hyperspectral remote sensing data. In this study, I used the theory of canopy spectral invariants (AKA p-theory) to understand the role of canopy structure and soil in the retrieval of foliar N from hyperspectral data and machine learning techniques. The results of this study showed the inconsistencies among different machine learning techniques used for estimating N. Using p-theory, I demonstrated that soil can contribute up to 95% to the total radiation budget of the canopy. I suggested an alternative approach to study photosynthesis is the use of dynamic global vegetation models (DGVMs). Gross primary production (GPP) is the apparent ecosystem scale photosynthesis that can be estimated using DGVMs. In this study, I performed a thorough sensitivity analysis and calibrated the Ecosystem Demography (EDv2.2) model along an elevation gradient in a dryland study area. I investigated the GPP capacity and activity by comparing the EDv2.2 GPP with flux towers and remote sensing products. The overall results showed that EDv2.2 performed well in capturing GPP capacity and its long term trend at lower elevation sites within the study area; whereas the model performed worse at higher elevations likely due to the change in vegetation community. I discussed that adding more heterogeneity and modifying ecosystem processes such as phenology and plant hydraulics in ED.v2.2 will improve its application to higher elevation ecosystems where there is more vegetation production. And finally, I developed an integrated hyperspectral-lidar framework for regional mapping of xeric and mesic vegetation in the study area. I showed that by considering spectral shape and magnitude, canopy structure and landscape features (riparian zone), we can develop a straightforward algorithm for vegetation mapping in drylands. This framework is simple, easy to interpret and consistent with our ecological understanding of vegetation distribution in drylands over large areas. Collectively, the results I present in this dissertation demonstrate the potential for advanced remote sensing and modeling to help us better understand ecosystem processes in drylands

Lidar Aboveground Vegetation Biomass Estimates in Shrublands: Prediction, Uncertainties and Application to Coarser Scales

Our study objectives were to model the aboveground biomass in a xeric shrub-steppe landscape with airborne light detection and ranging (Lidar) and explore the uncertainty associated with the models we created. We incorporated vegetation vertical structure information obtained from Lidar with ground-measured biomass data, allowing us to scale shrub biomass from small field sites (1 m subplots and 1 ha plots) to a larger landscape. A series of airborne Lidar-derived vegetation metrics were trained and linked with the field-measured biomass in Random Forests (RF) regression models. A Stepwise Multiple Regression (SMR) model was also explored as a comparison. Our results demonstrated that the important predictors from Lidar-derived metrics had a strong correlation with field-measured biomass in the RF regression models with a pseudo R2 of 0.76 and RMSE of 125 g/m2 for shrub biomass and a pseudo R2 of 0.74 and RMSE of 141 g/m2 for total biomass, and a weak correlation with field-measured herbaceous biomass. The SMR results were similar but slightly better than RF, explaining 77–79% of the variance, with RMSE ranging from 120 to 129 g/m2 for shrub and total biomass, respectively. We further explored the computational efficiency and relative accuracies of using point cloud and raster Lidar metrics at different resolutions (1 m to 1 ha). Metrics derived from the Lidar point cloud processing led to improved biomass estimates at nearly all resolutions in comparison to raster-derived Lidar metrics. Only at 1 m were the results from the point cloud and raster products nearly equivalent. The best Lidar prediction models of biomass at the plot-level (1 ha) were achieved when Lidar metrics were derived from an average of fine resolution (1 m) metrics to minimize boundary effects and to smooth variability. Overall, both RF and SMR methods explained more than 74% of the variance in biomass, with the most important Lidar variables being associated with vegetation structure and statistical measures of this structure (e.g., standard deviation of height was a strong predictor of biomass). Using our model results, we developed spatially-explicit Lidar estimates of total and shrub biomass across our study site in the Great Basin, U.S.A., for monitoring and planning in this imperiled ecosystem

An assessment of tropical dryland forest ecosystem biomass and climate change impacts in the Kavango-Zambezi (KAZA) region of Southern Africa

The dryland forests of the Kavango-Zambezi (KAZA) region in Southern Africa are highly susceptible to disturbances from an increase in human population, wildlife pressures and the impacts of climate change. In this environment, reliable forest extent and structure estimates are difficult to obtain because of the size and remoteness of KAZA (519,912 km²). Whilst satellite remote sensing is generally well-suited to monitoring forest characteristics, there remain large uncertainties about its application for assessing changes at a regional scale to quantify forest structure and biomass in dry forest environments. This thesis presents research that combines Synthetic Aperture Radar, multispectral satellite imagery and climatological data with an inventory from a ground survey of woodland in Botswana and Namibia in 2019. The research utilised a multi-method approach including parametric and non-parametric algorithms and change detection models to address the following objectives: (1) To assess the feasibility of using openly accessible remote sensing data to estimate the dryland forest above ground biomass (2) to quantify the detail of vegetation dynamics using extensive archives of time series satellite data; (3) to investigate the relationship between fire, soil moisture, and drought on dryland vegetation as a means of characterising spatiotemporal changes in aridity. The results establish that a combination of radar and multispectral imagery produced the best fit to the ground observations for estimating forest above ground biomass. Modelling of the time-series shows that it is possible to identify abrupt changes, longer-term trends and seasonality in forest dynamics. The time series analysis of fire shows that about 75% of the study area burned at least once within the 17-year monitoring period, with the national parks more frequently affected than other protected areas. The results presented show a significant increase in dryness over the past 2 decades, with arid and semi-arid regions encroaching at the expense of dry sub-humid, particularly in the south of the region, notably between 2011-2019

Complexity and Dynamics of Semi-Arid Vegetation Structure, Function and Diversity Across Spatial Scales from Full Waveform Lidar

Semi-arid ecosystems cover approximately 40% of the earth’s terrestrial landscape and show high dynamicity in ecosystem structure and function. These ecosystems play a critical role in global carbon dynamics, productivity, and habitat quality. Semi-arid ecosystems experience a high degree of disturbance that can severely alter ecosystem services and processes. Understanding the structure-function relationships across spatial extents are critical in order to assess their demography, response to disturbance, and for conservation management. In this research, using state-of-the-art full waveform lidar (airborne and spaceborne) and field observations, I developed a framework to assess the complexity and dynamics of vegetation structure, function and diversity across spatial scales in a semi-arid ecosystem. Difficulty in differentiating low stature vegetation from bare ground is the key remote sensing challenge in semi-arid ecosystems. In this study, I developed a workflow to differentiate key plant functional types (PFTs) using both structural and biophysical variables derived from the full waveform lidar and an ensemble random forest technique. The results revealed that waveform lidar pulse width can clearly distinguish shrubs from bare ground. The models showed PFT classification accuracy of 0.81–0.86% and 0.60–0.70% at 10 m and 1 m spatial resolutions, respectively. I found that structural variables were more important than the biophysical variables to differentiate the PFTs in this study area. The study further revealed an overlap between the structural features of different PFTs (e.g. shrubs from trees). Using structural features, I derived three main functional traits (canopy height, plant area index and foliage height diversity) of shrubs and trees that describe canopy architecture and light use efficiency of the ecosystem. I evaluated the trends and patterns of functional diversity and their relationship with non-climatic abiotic factors and fire disturbance. In addition to the fine resolution airborne lidar, I used simulated large footprint spaceborne lidar representing the newly launched Global Ecosystem Dynamics Investigation system (GEDI, a lidar sensor on the International Space Station) to evaluate the potential of capturing functional diversity trends of semi-arid ecosystems at global scales. The consistency of diversity trends between the airborne lidar and GEDI confirmed GEDI’s potential to capture functional diversity. I found that the functional diversity in this ecosystem is mainly governed by the local elevation gradient, soil type, and slope. All three functional diversity indices (functional richness, functional evenness and functional divergence) showed a diversity breakpoint near elevations of 1500 m – 1700 m. Functional diversity of fire-disturbed areas revealed that the fires in our study area resulted in a more even and less divergent ecosystem state. Finally, I quantified aboveground biomass using the structural features derived from both the airborne lidar and GEDI data. Regional estimates of biomass can indicate whether an ecosystem is a net carbon sink or source as well as the ecosystem’s health (e.g. biodiversity). Further, the potential of large footprint lidar data to estimate biomass in semi-arid ecosystems are not yet fully explored due to the inherent overlapping vegetation responses in the ground signals that can be affected by the ground slope. With a correction to the slope effect, I found that large footprint lidar can explain 42% of variance of biomass with a RMSE of 351 kg/ha (16% RMSE). The model estimated 82% of the study area with less than 50% uncertainty in biomass estimates. The cultivated areas and the areas with high functional richness showed the highest uncertainties. Overall, this dissertation establishes a novel framework to assess the complexity and dynamics of vegetation structure and function of a semi-arid ecosystem from space. This work enhances our understanding of the present state of an ecosystem and provides a foundation for using full waveform lidar to understand the impact of these changes to ecosystem productivity, biodiversity and habitat quality in the coming decades. The methods and algorithms in this dissertation can be directly applied to similar ecosystems with relevant corrections for the appropriate sensor. In addition, this study provides insights to related NASA missions such as ICESat-2 and future NASA missions such as NISAR for deriving vegetation structure and dynamics related to disturbance

Peatland leaf-area index and biomass estimation with ultra-high resolution remote sensing

There is fine-scale spatial heterogeneity in key vegetation properties including leaf-area index (LAI) and biomass in treeless northern peatlands, and hyperspectral drone data with high spatial and spectral resolution could detect the spatial patterns with high accuracy. However, the advantage of hyperspectral drone data has not been tested in a multi-source remote sensing approach (i.e. inclusion of multiple different remote sensing datatypes); and overall, sub-meter-level leaf-area index (LAI) and biomass maps have largely been absent. We evaluated the detectability of LAI and biomass patterns at a northern boreal fen (Halssiaapa) in northern Finland with multi-temporal and multi-source remote sensing data and assessed the benefit of hyperspectral drone data. We measured vascular plant percentage cover and height as well as moss cover in 140 field plots and connected the structural information to measured aboveground vascular LAI and biomass and moss biomass with linear regressions. We predicted both total and plant functional type (PFT) specific LAI and biomass patterns with random forests regressions with predictors including RGB and hyperspectral drone (28 bands in a spectral range of 500-900 nm), aerial and satellite imagery as well as topography and vegetation height information derived from structure-from-motion drone photogrammetry and aerial lidar data. The modeling performance was between moderate and good for total LAI and biomass (mean explained variance between 49.8 and 66.5%) and variable for PFTs (0.3-61.6%). Hyperspectral data increased model performance in most of the regressions, usually relatively little, but in some of the regressions, the inclusion of hyperspectral data even decreased model performance (change in mean explained variance between -14.5 and 9.1%-points). The most important features in regressions included drone topography, vegetation height, hyperspectral and RGB features. The spatial patterns and landscape estimates of LAI and biomass were quite similar in regressions with or without hyperspectral data, in particular for moss and total biomass. The results suggest that the fine-scale spatial patterns of peatland LAI and biomass can be detected with multi-source remote sensing data, vegetation mapping should include both spectral and topographic predictors at sub-meter-level spatial resolution and that hyperspectral imagery gives only slight benefits.Peer reviewe

Using Sentinel-2 and canopy height models to derive a landscape-level biomass map covering multiple vegetation types

Vegetation biomass is a globally important climate-relevant terrestrial carbon pool and also drives local hydrological systems via evapotranspiration. Vegetation biomass of individual vegetation types has been successfully estimated from active and passive remote sensing data. However, for many tasks, landscape-level biomass maps across several vegetation types are more suitable than biomass maps of individual vegetation types. For example, the validation of ecohydrological models and carbon budgeting typically requires spatially continuous biomass estimates, independent from vegetation type. Studies that derive biomass estimates across multiple vegetation or land-cover types to merge them into a single landscape-level biomass map are still scarce, and corresponding workflows must be developed. Here, we present a workflow to derive biomass estimates on landscape-level for a large watershed in central Chile. Our workflow has three steps: First, we combine field plotbased biomass estimates with spectral and structural information collected from Sentinel-2, TanDEM-X and airborne LiDAR data to map grassland, shrubland, native forests and pine plantation biomass using random forest regressions with an automatic feature selection. Second, we predict all models to the entire landscape. Third, we derive a land-cover map including the four considered vegetation types. We then use this land-cover map to assign the correct vegetation type-specific biomass estimate to each pixel according to one of the four considered vegetation types. Using a single repeatable workflow, we obtained biomass predictions comparable to earlier studies focusing on only one of the four vegetation types (Spearman correlation between 0.80 and 0.84; normalized-RMSE below 16 % for all vegetation types). For all woody vegetation types, height metrics were amongst the selected predictors, while for grasslands, only Sentinel-2 bands were selected. The land-cover was also mapped with high accuracy (OA = 83.1 %). The final landscape-level biomass map spatially agrees well with the known biomass distribution patterns in the watershed. Progressing from vegetation-type specific maps towards landscape-level biomass maps is an essential step towards integrating remote-sensing based biomass estimates into models for water and carbon management