Estimation of forest stem volume using ALOS-2 PALSAR-2 satellite images

A first evaluation of ALOS-2 PALSAR-2 data for forest stem volume estimation has been performed at a coniferous dominated test site in southern Sweden. Both the Fine Beam Dual (FBD) polarization and the Quad-polarimetric mode were investigated. Forest plots with stem volume reaching up to a maximum of about 620 m3 ha-1 (corresponding to 370 tons ha-1) were analyzed by relating backscatter intensity to field data using an exponential model derived from the Water Cloud Model. The estimation accuracy of stem volume at plot level (0.5 ha) was calculated in terms of Root Mean Square Error (RMSE). For the best case investigated an RMSE of 43.1% was obtained using one of the FBD HV-polarized images. The corresponding RMSE for the FBD HH-polarized images was 43.9%. In the Quadpolarimetric mode the lowest RMSE at HV- and HHpolarization was found to be 39.8% and 47.4%, respectively

Estimation of change in forest variables using synthetic aperture radar

Large scale mapping of changes in forest variables is needed for both environmental monitoring, planning of climate actions and sustainable forest management. Remote sensing can be used in conjunction with field data to produce wall-to-wall estimates that are practically impossible to produce using traditional field surveys. Synthetic aperture radar (SAR) can observe the forest independent of sunlight, clouds, snow, or rain, providing reliable high frequency coverage. Its wavelength determines the interaction with the forest, where longer wavelengths interact with larger structures of the trees, and shorter wavelengths interact mainly with the top part of the canopy, meaning that it can be chosen to fit specific applications. This thesis contains five studies conducted on the Remningstorp test site in southern Sweden. Studies I – III predicted above ground biomass (AGB) change using long wavelength polarimetric P- (in I) and L-band (in I – III) SAR data. The differences between the bands were small in terms of prediction quality, and the HV polarization, just as for AGB state prediction, was the polarization channel most correlated with AGB change. A moisture correction for L-band data was proposed and evaluated, and it was found that certain polarimetric measures were better for predicting AGB change than all of the polarization channels together. Study IV assessed the detectability of silvicultural treatments in short wavelength TanDEM-X interferometric phase heights. In line with earlier studies, only clear cuts were unambiguously distinguishable. Study V predicted site index and stand age by fitting height development curves to time series of TanDEM-X data. Site index and age were unbiasedly predicted for untreated plots, and the RMSE would likely decrease with longer time series. When stand age was known, SI was predicted with an RMSE comparable to that of the field based measurements. In conclusion, this thesis underscores SAR data's potential for generalizable methods for estimation of forest variable changes

Timber production assessment of a plantation forest: An integrated framework with field-based inventory, multi-source remote sensing data and forest management history

Timber production is the purpose for managing plantation forests, and its spatial and quantitative information is critical for advising management strategies. Previous studies have focused on growing stock volume (GSV), which represents the current potential of timber production, yet few studies have investigated historical process-harvested timber. This resulted in a gap in a synthetical ecosystem service assessment of timber production. In this paper, we established a Management Process-based Timber production (MPT) framework to integrate the current GSV and the harvested timber derived from historical logging regimes, trying to synthetically assess timber production for a historical period. In the MPT framework, age-class and current GSV determine the times of historical thinning and the corresponding harvested timber, by using a "space-for-time" substitution. The total timber production can be estimated by the historical harvested timber in each thinning and the current GSV. To test this MPT framework, an empirical study on a larch plantation (LP) with area of 43,946 ha was conducted in North China for a period from 1962 to 2010. Field-based inventory data was integrated with ALOS PALSAR (Advanced Land-Observing Satellite Phased Array L-band Synthetic Aperture Radar) and Landsat-8 OLI (Operational Land Imager) data for estimating the age-class and current GSV of LP. The random forest model with PALSAR backscatter intensity channels and OLI bands as input predictive variables yielded an accuracy of 67.9% with a Kappa coefficient of 0.59 for age-class classification. The regression model using PALSAR data produced a root mean square error (RMSE) of 36.5 m(3) ha(-1). The total timber production of LP was estimated to be 7.27 x 10(6) m(3), with 4.87 x 10(6) m(3) in current GSV and 2.40 x 10(6) m(3) in harvested timber through historical thinning. The historical process-harvested timber accounts to 33.0% of the total timber production, which component has been neglected in the assessments for current status of plantation forests. Synthetically considering the RMSE for predictive GSV and misclassification of age-class, the error in timber production were supposed to range from -55.2 to 56.3 m(3) ha(-1). The MPT framework can be used to assess timber production of other tree species at a larger spatial scale, providing crucial information for a better understanding of forest ecosystem service. (C) 2016 Elsevier B.V. All rights reserved.ArticleINTERNATIONAL JOURNAL OF APPLIED EARTH OBSERVATION AND GEOINFORMATION.52:155-165(2016)journal articl

Predictions of Biomass Change in a Hemi-Boreal Forest Based on Multi-Polarization L- and P-Band SAR Backscatter

Above-ground biomass change accumulated during four growth seasons in a hemi-boreal forest was predicted using airborne L- and P-band synthetic aperture radar (SAR) backscatter. The radar data were collected in the BioSAR 2007 and BioSAR 2010 campaigns over the Remningstorp test site in southern Sweden. Regression models for biomass change were developed from biomass maps created using airborne LiDAR data and field measurements. To facilitate training and prediction on image pairs acquired at different dates, a backscatter offset correction method for L-band data was developed and evaluated. The correction, based on the HV/VV backscatter ratio, facilitated predictions across image pairs almost identical to those obtained using data from the same image pair for both training and prediction. For P-band, previous positive results using an offset correction based on the HH/VV ratio were validated. The best L-band model achieved a root mean square error (RMSE) of 21 t/ha, and the best P-band model achieved an RMSE of 19 t/ha. Those accuracies are similar to that of the LiDAR-based biomass change of 18 t/ha. The limitation of using LiDAR-based data for training was considered. The findings demonstrate potential for improved biomass change predictions from L-band backscatter despite varying environmental conditions and calibration uncertainties

Remote Sensing of Savannas and Woodlands

Savannas and woodlands are one of the most challenging targets for remote sensing. This book provides a current snapshot of the geographical focus and application of the latest sensors and sensor combinations in savannas and woodlands. It includes feature articles on terrestrial laser scanning and on the application of remote sensing to characterization of vegetation dynamics in the Mato Grosso, Cerrado and Caatinga of Brazil. It also contains studies focussed on savannas in Europe, North America, Africa and Australia. It should be important reading for environmental practitioners and scientists globally who are concerned with the sustainability of the global savanna and woodland biome

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

Calibration of DART Radiative Transfer Model with Satellite Images for Simulating Albedo and Thermal Irradiance Images and 3D Radiative Budget of Urban Environment

Remote sensing is increasingly used for managing urban environment. In this context, the H2020 project URBANFLUXES aims to improve our knowledge on urban anthropogenic heat fluxes, with the specific study of three cities: London, Basel and Heraklion. Usually, one expects to derive directly 2 major urban parameters from remote sensing: the albedo and thermal irradiance. However, the determination of these two parameters is seriously hampered by complexity of urban architecture. For example, urban reflectance and brightness temperature are far from isotropic and are spatially heterogeneous. Hence, radiative transfer models that consider the complexity of urban architecture when simulating remote sensing signals are essential tools. Even for these sophisticated models, there is a major constraint for an operational use of remote sensing: the complex 3D distribution of optical properties and temperatures in urban environments. Here, the work is conducted with the DART (Discrete Anisotropic Radiative Transfer) model. It is a comprehensive physically based 3D radiative transfer model that simulates optical signals at the entrance of imaging spectro-radiometers and LiDAR scanners on board of satellites and airplanes, as well as the 3D radiative budget, of urban and natural landscapes for any experimental (atmosphere, topography,…) and instrumental (sensor altitude, spatial resolution, UV to thermal infrared,…) configuration. Paul Sabatier University distributes free licenses for research activities. This paper presents the calibration of DART model with high spatial resolution satellite images (Landsat 8, Sentinel 2, etc.) that are acquired in the visible (VIS) / near infrared (NIR) domain and in the thermal infrared (TIR) domain. Here, the work is conducted with an atmospherically corrected Landsat 8 image and Bale city, with its urban database. The calibration approach in the VIS/IR domain encompasses 5 steps for computing the 2D distribution (image) of urban albedo at satellite spatial resolution. (1) DART simulation of satellite image at very high spatial resolution (e.g., 50cm) per satellite spectral band. Atmosphere conditions are specific to the satellite image acquisition. (2) Spatial resampling of DART image at the coarser spatial resolution of the available satellite image, per spectral band. (3) Iterative derivation of the urban surfaces (roofs, walls, streets, vegetation,…) optical properties as derived from pixel-wise comparison of DART and satellite images, independently per spectral band. (4) Computation of the band albedo image of the city, per spectral band. (5) Computation of the image of the city albedo and VIS/NIR exitance, as an integral over all satellite spectral bands. In order to get a time series of albedo and VIS/NIR exitance, even in the absence of satellite images, ECMWF information about local irradiance and atmosphere conditions are used. A similar approach is used for calculating the city thermal exitance using satellite images acquired in the thermal infrared domain. Finally, DART simulations that are conducted with the optical properties derived from remote sensing images give also the 3D radiative budget of the city at any date including the date of the satellite image acquisition