The SAR Handbook: Comprehensive Methodologies for Forest Monitoring and Biomass Estimation

This Synthetic Aperture Radar (SAR) handbook of applied methods for forest monitoring and biomass estimation has been developed by SERVIR in collaboration with SilvaCarbon to address pressing needs in the development of operational forest monitoring services. Despite the existence of SAR technology with all-weather capability for over 30 years, the applied use of this technology for operational purposes has proven difficult. This handbook seeks to provide understandable, easy-to-assimilate technical material to remote sensing specialists that may not have expertise on SAR but are interested in leveraging SAR technology in the forestry sector

Using satellite estimates of aboveground biomass to assess carbon stocks in a mixed-management, semi-deciduous tropical forest in the Yucatan Peninsula

Information on the spatial distribution of forest aboveground biomass (AGB) and its uncertainty is important to evaluate management and conservation policies in tropical forests. However, the scarcity of field data and robust protocols to propagate uncertainty prevent a robust estimation through remote sensing. We upscaled AGB from field data to LiDAR, and to landscape scale using Sentinel-2 and ALOS-PALSAR through machine learning, propagated uncertainty using a Monte Carlo framework and explored the relative contributions of each sensor. Sentinel-2 outperformed ALOS-PALSAR (R2 = 0.66, vs 0.50), however, the combination provided the best fit (R2 = 0.70). The combined model explained 49% of the variation comparing against plots within the calibration area, and 17% outside, however, 94% of observations outside calibration area fell within the 95% confidence intervals. Finally, we partitioned the distribution of AGB in different management and conservation categories for evaluating the potential of different strategies for conserving carbon stock

Forest Biomass and Land Cover Change Assessment of the Margalla Hills National Park in Pakistan Using a Remote Sensing Based Approach

Climate change is one of the greatest threats recently, of which the developing countries are facing most of the brunt. In the fight against climate change, forests can play an important role, since they hold a substantial amount of terrestrial carbon and can therefore affect the global carbon cycle. Forests are also an essential source of livelihood for a remarkably high proportion of people worldwide and a harbor for rich global biodiversity. Forests are however facing high deforestation rates. Deforestation is regarded as the most widespread process of land cover change (LCC), which is the conversion of one land cover type to the other land cover type. Most of this deforestation occurs in developing countries. Agricultural expansion has been reported as the most significant widespread driver of deforestation in Asia, Africa, and Latin America. This deforestation is altering the balance of forest carbon stocks and threatening biodiversity. Pakistan is also a low forest cover country and faces high deforestation rates at the same time, due to the high reliance of local communities on forests. Moreover, it is also the most adversely affected by climate change. Agricultural expansion and population growth have been regarded as the most common drivers of deforestation in Pakistan. Financial incentives such as ‘Reducing Emissions from Deforestation and Forest Degradation, and the Role of Conservation of Forest Carbon, Sustainable Management of Forests and Enhancement of Forest Carbon Stocks’ (REDD+) offer hope for developing countries for not only halting deforestation but also alleviating poverty. However, such initiatives require the estimation of biomass and carbon stocks of the forest ecosystems. Therefore, it becomes necessary that the biomass and carbon potentials of the forests are explored, as well as the LCCs are investigated for identifying the deforestation and forest degradation hit areas. Based on the aforementioned, the following research objectives/sub-objectives were investigated in the MHNP, which is adjoined with the capital city of Pakistan, Islamabad; A) Forest Biomass and Carbon Stock Assessment of Margalla Hills National Park (MHNP) A.1) Aboveground Biomass (AGB) and Aboveground Carbon (AGC) assessment of the Subtropical Chir Pine Forest (SCPF) and Subtropical Broadleaved Evergreen Forest (SBEF) using Field Inventorying Techniques A.2) Exploring linear regression relationship between Sentinel-1 (S1) and Sentinel-2 (S2) satellite data with the AGB of SCPF and SBEF A.3) AGB estimation combining remote sensing and machine learning approach B) LC Classification and Land Cover Change Detection (LCCD) of MHNP for the time-period between 1999 and 2019 B.1) LC Classification for the years 1999, 2009 and 2019 using Machine Learning Algorithm B.2) LCCD of MHNP between 1999 to 2019

Modelling and Predicting the Growing Stock Volume in Small-Scale Plantation Forests of Tanzania Using Multi-Sensor Image Synergy

Remotely sensed assisted forest inventory has emerged in the past decade as a robust and cost efficient method for generating accurate information on forest biophysical parameters. The launching and public access of ALOS PALSAR-2, Sentinel-1 (SAR), and Sentinel-2 together with the associated open-source software, has further increased the opportunity for application of remotely sensed data in forest inventories. In this study, we evaluated the ability of ALOS PALSAR-2, Sentinel-1 (SAR) and Sentinel-2 and their combinations to predict growing stock volume in small-scale forest plantations of Tanzania. The effects of two variable extraction approaches (i.e., centroid and weighted mean), seasonality (i.e., rainy and dry), and tree species on the prediction accuracy of growing stock volume when using each of the three remotely sensed data were also investigated. Statistical models relating growing stock volume and remotely sensed predictor variables at the plot-level were fitted using multiple linear regression. The models were evaluated using the k-fold cross validation and judged based on the relative root mean square error values (RMSEr). The results showed that: Sentinel-2 (RMSEr = 42.03% and pseudo − R 2 = 0.63) and the combination of Sentinel-1 and Sentinel-2 (RMSEr = 46.98% and pseudo − R 2 = 0.52), had better performance in predicting growing stock volume, as compared to Sentinel-1 (RMSEr = 59.48% and pseudo − R 2 = 0.18) alone. Models fitted with variables extracted from the weighted mean approach, turned out to have relatively lower RMSEr % values, as compared to centroid approaches. Sentinel-2 rainy season based models had slightly smaller RMSEr values, as compared to dry season based models. Dense time series (i.e., annual) data resulted to the models with relatively lower RMSEr values, as compared to seasonal based models when using variables extracted from the weighted mean approach. For the centroid approach there was no notable difference between the models fitted using dense time series versus rain season based predictor variables. Stratifications based on tree species resulted into lower RMSEr values for Pinus patula tree species, as compared to other tree species. Finally, our study concluded that combination of Sentinel-1&2 as well as the use Sentinel-2 alone can be considered for remote-sensing assisted forest inventory in the small-scale plantation forests of Tanzania. Further studies on the effect of field plot size, stratification and statistical methods on the prediction accuracy are recommended. </p

Biomass Estimation Using Satellite-Based Data

Comprehensive measurements of global forest aboveground biomass (AGB) are crucial information to promote the sustainable management of forests to mitigate climate change and preserve the multiple ecosystem services provided by forests. Optical and radar sensors are available at different spatial, spectral, and temporal scales. The integration of multi-sources sensor data with field measurements, using appropriated algorithms to identify the relationship between remote sensing predictors and reference measurements, is important to improve forest AGB estimation. This chapter aims to present different types of predicted variables derived from multi-sources sensors, such as original spectral bands, transformed images, vegetation indices, textural features, and different regression algorithms used (parametric and non-parametric) that contribute to a more robust, practical, and cost-effective approach for forest AGB estimation at different levels

Characterisation and monitoring of forest disturbances in Ireland using active microwave satellite platforms

Forests are one of the major carbon sinks that significantly contribute towards achieving targets of the Kyoto Protocol, and its successors, in reducing greenhouse (GHG) emissions. In order to contribute to regular National Inventory Reporting, and as part of the on-going development of the Irish national GHG reporting system (CARBWARE), improvements in characterisation of changes in forest carbon stocks have been recommended to provide a comprehensive information flow into CARBWARE. The Irish National Forest Inventory (NFI) is updated once every six years, thus there is a need for an enhanced forest monitoring system to obtain annual forest updates to support government agencies and forest management companies in their strategic decision making and to comply with international GHG reporting standards. Sustainable forest management is imperative to promote net carbon absorption from forests. Based on the NFI data, Irish forests have removed or sequestered an average of 3.8 Mt of atmospheric CO2 per year between 2007 and 2016. However, unmanaged and degraded forests become a net emitter of carbon. Disturbances from human induced activities such as clear felling, thinning and deforestation results in carbon emissions back into the atmosphere. Funded by the Department of Agriculture, Food and the Marine (DAFM, Ireland), this PhD study focuses on exploring the potential of data from L-band Synthetic Aperture Radar (SAR) satellite based sensors for monitoring changes in the small stand forests of Ireland. Historic data from ALOS PALSAR in the late 2000s and more recent data from ALOS-2 PALSAR-2 sensors have been used to map forest areas and characterise the different disturbances observed within three different regions of Ireland. Forest mapping and disturbance characterisation was achieved by combining the machine learning supervised Random Forests (RF) and unsupervised Iterative Self-Organizing Data Analysis (ISODATA) classification techniques. The lack of availability of ground truth data supported use of this unsupervised approach which forms natural clusters based on their multi-temporal signatures, with divergence statistics used to select the optimal number of clusters to represent different forest classes. This approach to forest monitoring using SAR imagery has not been reported in the peer-review literature and is particularly beneficial where there is a dearth of ground-based information. When applied to the forests, mapped with an accuracy of up to 97% by RF, the ISODATA technique successfully identified the unique multi-temporal pattern associated with clear-fells which exhibited a decrease of 4 to 5 decibels (dB) between the images acquired before and after the event. The clustering algorithm effectively highlighted the occurrence of other disturbance events within forests with a decrease of 2±0.5dB between two consecutive years, as well as areas of tree growth and afforestation. A highlight of the work is the successful transferability of the algorithm, developed using ALOS PALSAR, to ALOS-2 PALSAR-2 data thereby demonstrating the potential continuity of annual forest monitoring. The higher spatial and radiometric resolutions of ALOS-2 PALSAR-2 data have shown improvements in forest mapping compared to ALOS PALSAR data. From mapping a minimum forest size of 1.8 ha with ALOS PALSAR, a minimum area of 1.1 ha was achieved with the ALOS-2 PALSAR-2 images. Moreover, even with some different backscatter characteristics of images acquired in different seasons, similar signature patterns between the sensors were retrieved that helped to define the cluster groups, thus demonstrating the robustness of the algorithm and its successful transferability. Having proven the potential to monitor forest disturbances, the results from both the sensors were used to detect deforestation over the time period 2007-2016. Permanent land-use changes pertaining to conversion of forests to agricultural lands and windfarms were identified which are important with respect to forest monitoring and carbon reporting in Ireland. Overall, this work has presented a viable approach to support forest monitoring operations in Ireland. By providing disturbance information from SAR, it can supplement projects working with optical images which are generally limited by cloud cover, particularly in parts of northern, western and upland Ireland. This approach adds value to ground based forest monitoring by mapping distinct forests over large areas on an annual basis. This study has demonstrated the ability to apply the algorithm to three different study areas, with a vision to operationalise the algorithm on a national scale. The main limitations experienced in this study were the lack of L-band SAR data availability and reference datasets. With typically only one image acquired per year, and discrepancies and omissions existing within reference datasets, understanding the behaviour of certain cluster groups representing disturbances was challenging. However, this approach has addressed some issues within the reference datasets, for example locating areas for which a felling licence was granted but where trees were never cut, by providing detailed systematic mapping of forests. Future satellites such as Tandem-L, SAOCOM-2A and 2B, P-band BIOMASS mission and ALOS-4 PALSAR-3 may overcome the issue of limited SAR image acquisitions provided more images per year are available, especially during the summer months

Modelling Above Ground Biomass Using Sentinel 2 and Planet Scope Data in Dense Tropical Montane Forests of Tanzania

Forest biomass estimation using field -based inventories at a large scale is challenging and generally entails large uncertainty in tropical regions. In this study, we investigated the performance of Sentinel 2 and Planet Scope data for above ground biomass (AGB) modelling, in the tropical rainforest of Tanzania. A total of 296 field inventory plots were measured across the west Usambara mountain forests. The results showed that, Sentinel 2-based model fitted using GLMs had better performance (cvRMSEr = 67.00 % and pseudo-R2= 20%) as compared to Planet Scope-based models (cvRMSEr = 72.1 % and pseudo-R2= 5.2%). Overall GLMs resulted into models with less prediction errors in contrast to random forest when using Sentinel 2 data. However, for the Planet Scope, there was marginal improvement when using random forest (cvRMSEr = 72.0%). Models that incorporated texture variables produced better prediction accuracy as compared to those with band values and indices only. The study has shown that, Sentinel 2 and Planet Scope remotely sensed data can be used to develop cost-effective method for AGB estimation in tropical rainforests of Tanzania

TerraSAR-X and Wetlands: A Review

Since its launch in 2007, TerraSAR-X observations have been widely used in a broad range of scientific applications. Particularly in wetland research, TerraSAR-X\u27s shortwave X-band synthetic aperture radar (SAR) possesses unique capabilities, such as high spatial and temporal resolution, for delineating and characterizing the inherent spatially and temporally complex and heterogeneous structure of wetland ecosystems and their dynamics. As transitional areas, wetlands comprise characteristics of both terrestrial and aquatic features, forming a large diversity of wetland types. This study reviews all published articles incorporating TerraSAR-X information into wetland research to provide a comprehensive study of how this sensor has been used with regard to polarization, and the function of the data, time-series analyses, or the assessment of specific wetland ecosystem types. What is evident throughout this literature review is the synergistic fusion of multi-frequency and multi-polarization SAR sensors, sometimes optical sensors, in almost all investigated studies to attain improved wetland classification results. Due to the short revisiting time of the TerraSAR-X sensor, it is possible to compute dense SAR time-series, allowing for a more precise observation of the seasonality in dynamic wetland areas as demonstrated in many of the reviewed studies

Biomass Representation in Synthetic Aperture Radar Interferometry Data Sets

This work makes an attempt to explain the origin, features and potential applications of the elevation bias of the synthetic aperture radar interferometry (InSAR) datasets over areas covered by vegetation. The rapid development of radar-based remote sensing methods, such as synthetic aperture radar (SAR) and InSAR, has provided an alternative to the photogrammetry and LiDAR for determining the third dimension of topographic surfaces. The InSAR method has proved to be so effective and productive that it allowed, within eleven days of the space shuttle mission, for acquisition of data to develop a three-dimensional model of almost the entire land surface of our planet. This mission is known as the Shuttle Radar Topography Mission (SRTM). Scientists across the geosciences were able to access the great benefits of uniformity, high resolution and the most precise digital elevation model (DEM) of the Earth like never before for their a wide variety of scientific and practical inquiries. Unfortunately, InSAR elevations misrepresent the surface of the Earth in places where there is substantial vegetation cover. This is a systematic error of unknown, yet limited (by the vertical extension of vegetation) magnitude. Up to now, only a limited number of attempts to model this error source have been made. However, none offer a robust remedy, but rather partial or case-based solutions. More work in this area of research is needed as the number of airborne and space-based InSAR elevation models has been steadily increasing over the last few years, despite strong competition from LiDAR and optical methods. From another perspective, however, this elevation bias, termed here as the “biomass impenetrability”, creates a great opportunity to learn about the biomass. This may be achieved due to the fact that the impenetrability can be considered a collective response to a few factors originating in 3D space that encompass the outermost boundaries of vegetation. The biomass, presence in InSAR datasets or simply the biomass impenetrability, is the focus of this research. The report, presented in a sequence of sections, gradually introduces terminology, physical and mathematical fundamentals commonly used in describing the propagation of electromagnetic waves, including the Maxwell equations. The synthetic aperture radar (SAR) and InSAR as active remote sensing methods are summarised. In subsequent steps, the major InSAR data sources and data acquisition systems, past and present, are outlined. Various examples of the InSAR datasets, including the SRTM C- and X-band elevation products and INTERMAP Inc. IFSAR digital terrain/surface models (DTM/DSM), representing diverse test sites in the world are used to demonstrate the presence and/or magnitude of the biomass impenetrability in the context of different types of vegetation – usually forest. Also, results of investigations carried out by selected researchers on the elevation bias in InSAR datasets and their attempts at mathematical modelling are reviewed. In recent years, a few researchers have suggested that the magnitude of the biomass impenetrability is linked to gaps in the vegetation cover. Based on these hints, a mathematical model of the tree and the forest has been developed. Three types of gaps were identified; gaps in the landscape-scale forest areas (Type 1), e.g. forest fire scares and logging areas; a gap between three trees forming a triangle (Type 2), e.g. depending on the shape of tree crowns; and gaps within a tree itself (Type 3). Experiments have demonstrated that Type 1 gaps follow the power-law density distribution function. One of the most useful features of the power-law distributed phenomena is their scale-independent property. This property was also used to model Type 3 gaps (within the tree crown) by assuming that these gaps follow the same distribution as the Type 1 gaps. A hypothesis was formulated regarding the penetration depth of the radar waves within the canopy. It claims that the depth of penetration is simply related to the quantisation level of the radar backscattered signal. A higher level of bits per pixels allows for capturing weaker signals arriving from the lower levels of the tree crown. Assuming certain generic and simplified shapes of tree crowns including cone, paraboloid, sphere and spherical cap, it was possible to model analytically Type 2 gaps. The Monte Carlo simulation method was used to investigate relationships between the impenetrability and various configurations of a modelled forest. One of the most important findings is that impenetrability is largely explainable by the gaps between trees. A much less important role is played by the penetrability into the crown cover. Another important finding is that the impenetrability strongly correlates with the vegetation density. Using this feature, a method for vegetation density mapping called the mean maximum impenetrability (MMI) method is proposed. Unlike the traditional methods of forest inventories, the MMI method allows for a much more realistic inventory of vegetation cover, because it is able to capture an in situ or current situation on the ground, but not for areas that are nominally classified as a “forest-to-be”. The MMI method also allows for the mapping of landscape variation in the forest or vegetation density, which is a novel and exciting feature of the new 3D remote sensing (3DRS) technique. Besides the inventory-type applications, the MMI method can be used as a forest change detection method. For maximum effectiveness of the MMI method, an object-based change detection approach is preferred. A minimum requirement for the MMI method is a time-lapsed reference dataset in the form, for example, of an existing forest map of the area of interest, or a vegetation density map prepared using InSAR datasets. Preliminary tests aimed at finding a degree of correlation between the impenetrability and other types of passive and active remote sensing data sources, including TerraSAR-X, NDVI and PALSAR, proved that the method most sensitive to vegetation density was the Japanese PALSAR - L-band SAR system. Unfortunately, PALSAR backscattered signals become very noisy for impenetrability below 15 m. This means that PALSAR has severe limitations for low loadings of the biomass per unit area. The proposed applications of the InSAR data will remain indispensable wherever cloud cover obscures the sky in a persistent manner, which makes suitable optical data acquisition extremely time-consuming or nearly impossible. A limitation of the MMI method is due to the fact that the impenetrability is calculated using a reference DTM, which must be available beforehand. In many countries around the world, appropriate quality DTMs are still unavailable. A possible solution to this obstacle is to use a DEM that was derived using P-band InSAR elevations or LiDAR. It must be noted, however, that in many cases, two InSAR datasets separated by time of the same area are sufficient for forest change detection or similar applications