Next Article in Journal
Modeling Irrigation Water Requirement of Mixed Crop with Coupled Smart Irrigation System and System Dynamic Model
Previous Article in Journal
Evaluation of Operational Projects Supported by Cohesion Funds for the National Forest Parks of Greece
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Proceeding Paper

Estimation of Tree Height in Burned Areas with GEDI Laser Data in Northern Portugal and Galicia (Spain) †

by
Saulo Folharini
*,
Sarah Moura Batista dos Santos
,
António Bento-Gonçalves
and
António Vieira
Communication and Society Research Centre (CECS), The University of Minho, 4800-058 Guimarães, Portugal
*
Author to whom correspondence should be addressed.
Presented at the 3rd International Electronic Conference on Forests — Exploring New Discoveries and New Directions in Forests, 15–31 October 2022; Available online: https://iecf2022.sciforum.net/.
Environ. Sci. Proc. 2022, 22(1), 50; https://doi.org/10.3390/IECF2022-13050
Published: 15 October 2022

Abstract

:
We analyzed the tree height of four areas affected by wildfires in northern Portugal and Galicia (Spain) using the GEE platform, random forest regression model, and GEDI data. Before the fires, the heights varied from 5.21 to 20.16 m, with r2 values of 0.82 and 0.67. After the fires, heights of 5.55 to 9.12 m were recorded, with values of r2 0.47 and 1. These r2 values after fires indicate the absence or limitation of sample data. We concluded that the GEDI data have great potential to assist in the mapping of areas affected by wildfires.

1. Introduction

Large wildfires constitute a very significant disturbance to fire-prone landscapes. In different locations around the world, such as California (USA), Portugal, Spain, and Australia, the intensity of these natural events is enhanced by the effects of climate change and the overexploitation of natural resources by the population [1,2,3,4,5].
Wildfires are events that result in a great loss of biodiversity, in addition to important changes to the surface where they occur [6,7]. Portugal is a country where the recurrence of forest fires is increasing due to a combination of environmental and human factors, such as drier years, and climatic conditions favorable to their propagation, such as lack of precipitation, low humidity, and wind. Considering the human factor, in the last three decades, the unequal distribution of the population throughout the Portuguese territory is related to the regional distribution of ignition sources [8]. The same scenario is observed in the northwest region of Spain.
The consequences of the fires are of social and economic order, affecting rural properties and their agricultural production, resulting in a significant loss of capital, and of environmental order, degrading the landscape’s natural values and decreasing the provision of ecosystem services. The forest is responsible for providing different ecosystem services, such as carbon storage in vegetation and soil, food, water filtration, and climate regulation at local and regional levels. With the increase in the occurrence of wildfires, the availability of ecosystem services is altered [9].
In the context of climate change, quantifying carbon stocks, an important ecosystem service to mitigate global warming, is essential. To quantify the carbon stock, biophysical parameters, such as the height of the trees are necessary. To measure this parameter, there is the possibility of using a worldwide database called the Global Ecosystem Dynamics Investigation (GEDI), a LIDAR instrument on board the International Space Station that has made data available since April 2019, with a scan diameter of 25 m. This instrument makes high-resolution observations of the vertical structure of the forest, providing metrics needed to understand processes and patterns of biodiversity [10]. The studies of [11,12,13] used GEDI data to estimate the height of trees in forests. We used the Google Earth Engine platform to analyze satellite imagery, digital elevation model, and GEDI data to measure vegetation height before and after fires.

2. Material and Method

2.1. Study Area

The study was performed in 4 areas affected by wildfires in the summer of 2020 in northern Portugal and Galicia (Spain) (Figure 1).
The areas cover 1294.8 hectares (ICNF, 2020), with one of the areas affecting the Peneda-Gerês National Park (PNPG), a protected area together with the Spanish Natural Park of Baixa Limia–Serra do Xurés, the Gerês-Xurés Transfrontier Park, and the Biosphere reserve of the same name.

2.2. Data Acquisition and Processing

The research steps were adapted from the Kamusoko [13] and are available in Figure 2.
Sentinel data for the period of September to December 2020 were used to model vegetation structural characteristics. The Sentinel-2 multispectral sensor data at level 2A (Surface Reflectance) were processed together with C-band Synthetic Aperture Radar (SAR) data from Sentinel-1 (ground range detected). Both Sentinel datasets have a geometric resolution of 10 m, which will be the target spatial resolution of the models [14]. Sentinel-1 data were filtered for single copolarization of vertical transmission (VV) and dual-band cross-polarization in vertical transmission and horizontal reacceptance (VH). Sentinel-2 has been filtered for a cloud percentage of less than 20%.
Shuttle Radar Topography Mission (SRTM) [15] data were used to obtain elevation and slope values. This SRTM V3 product (SRTM Plus) was provided by NASA at a resolution of 1 arc-second (approximately 30 m). The World Cover 2020 product of the European Space Agency (ESA) was used for the delimitation of the forest mask; this image had a 10m resolution based on Sentinel-1 and Sentinel-2 data.
We used the Global Ecosystem Dynamics Investigation (GEDI) level 2A data [11,16] for forest area selection. GEDI data were stored as point geometries representing the centroid of the 25 m sample footprint. Vegetation structural attribute that was modelled was canopy height (98th percentile of the relative height metrics from GEDI L2A).
All data were processed in Google Earth Engine (GEE) using the Python API. The predictor variables for the model were determined from Sentinel-1, Sentinel-2, SRTM, and World Cover data. A machine learning approach was used with the application of a random forest regression model. A random forest (RF) model is a machine learning ensemble technique that has proven its popularity in remote sensing applications due to its high accuracy and performance [14,17].
In total, about 10,000 quality-filtered GEDI samples were used for model training and validation of each vegetation structure attribute. The GEDI sample set was randomly split before modeling into a collection of training samples (70%) and test samples (30%) for model validation. The bands included in the classification were (‘B2’, ‘B3’, ‘B4’, ‘B5’, ‘B6’, ‘B7’, ‘B8’, ‘B11’, ‘B12’, ‘VV_iqr’, ‘VH_iqr’, ‘elevation’, ‘slope’, and ‘Map’). The test samples, model-independent, were used to calculate several model efficiency criteria, such as the coefficient of determination (R2) and root mean square error (RMSE).

3. Results and Discussion

The importance of the variables used in the regression using random forest can be seen in Figure 3.
B4 (red) was the one that underwent the greatest change, as it is the most sensitive to vegetation variations in the visible spectrum. Therefore, the lack of or decrease in the amount of vegetation after fires resulted in a greater variation in its importance in the regression.
The wavelengths of the Red-edge range (B5, B6, and B7) and B8 (NIR) are more sensitive in identifying the severity of fires. The bands of the Red-edge range can be used in the calculation of the NDVI705 (Difference Vegetation Red-Edge Index) most suitable for forest monitoring and wildfires [18,19].
Sentinel’s SWIR waveform is sensitive to changes in the structural crust, deformation of soil minerals, and soil salinity. These changes are identified in soil degradation processes resulting from wildfires that result in greater soil exposure [20], hence the greater variation in the importance of these variables in the regression. More details on soil degradation using hyperspectral data and remote sensing can be found in Goldshleger et al. [21].
Sentinel-1 collects C-band synthetic aperture radar (SAR) imagery to provide VV and VH with dual-band cross-polarization and vertically transmitted/horizontally received data applied to forestry studies to identify fire scars. The results show that the VV polarization is more important than the VH, maintaining the same proportion before and after fires. Regarding terrain data, elevation is no longer an important variable in the postfire regression, and the slope becomes one of the most important variables. This is because the slope becomes more easily identified in a digital elevation model (DEM) in an area with fewer trees. The slope values help to identify regions where the trees did not burn.
Concerning the height of the trees, the results indicate that before the fires, the heights varied from 5.21 to 20.16 m and between 5.55 and 9.12 m after the fires. The regression results of the training (70%) and validation (30%) datasets before the fires are shown in the Figure 4.
In the training data set before the fires, a good result of r2 0.82 and RMSE 3.47 was recorded, with a greater concentration of data close to the fit line and in the validation data set, r2 0.67 and RMSE 5.23, with greater dispersion of the data to the surroundings of the model fit line.
After the fires, the training and validation data are shown in the Figure 5.
Concerning the regression of the training data observed after the fires, low values of r2 were identified in the training data of r2 0.47 and RMSE 3.36 with a large data dispersion and in the validation data of r2 1 and RMSE 3.34, constituting bad statistical results. However, these results are acceptable in the context of a work that intends to measure the height of trees in burned areas before and after fires. The result is explained by the drastic reduction in trees for model validation after the fires, dramatically decreasing the number of predictor variable data (GEDI) needed for random forest regression.

4. Conclusions

The use of the GEE platform to analyze large amounts of data has shown potential in recent years, due to its flexibility in code construction and data crossing. In the present study, we sought to analyze the GEDI dataset, an important LIDAR instrument that since 2019 has provided data from planet Earth with great potential for analysis related to climate change, such as measuring carbon stocks and vegetation structures.
We can conclude that the results presented in the present research are important and have great potential for the application in the study of wildfires using remote sensing. It was possible to identify with relative precision the variations observed in the height of trees in burned areas in northern Portugal and Spain.
The regression results corroborate the decrease in tree height due to the occurrence of fires, demonstrating great potential for the application of the random forest model to predict tree height using LIDAR data and other predictor variables, such as SAR and hyperspectral images from the Sentinel satellite.

Author Contributions

Conceptualization, S.F., S.M.B.d.S., A.B.-G. and A.V.; methodology, software, validation, and formal analysis, S.F. and S.M.B.d.S.; investigation and resources, S.F., S.M.B.d.S., A.B.-G. and A.V.; data curation and writing—original draft preparation, S.F. and S.M.B.d.S.; writing—review and editing, A.B.-G. and A.V.; visualization and supervision, A.B.-G. and A.V.; project administration, A.V.; funding acquisition, A.V. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the European Regional Development Fund. Climate Change Resilient Tourism in Protected Areas of Northern Portugal (CLICTOUR-Project NORTE-01-0145-FEDER-000079).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Sentinel-1 SAR GRD: C-band Synthetic Aperture Radar Ground Range Detected, available at: https://developers.google.com/earth-engine/datasets/catalog/COPERNICUS_S1_GRD; NASA SRTM Digital Elevation 30m, available at: https://developers.google.com/earth-engine/datasets/catalog/USGS_SRTMGL1_003; GEDI L2A Raster Canopy Top Height (Version 2), available at: https://developers.google.com/earth-engine/datasets/catalog/LARSE_GEDI_GEDI02_A_002_MONTHLY.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Nauslar, N.; Abatzoglou, J.; Marsh, P. The 2017 North Bay and Southern California Fires: A Case Study. Fire 2018, 1, 18. [Google Scholar] [CrossRef] [Green Version]
  2. Russell-Smith, J.; McCaw, L.; Leavesley, A. Adaptive prescribed burning in Australia for the early 21st Century—Context, status, challenges. Int. J. Wildl. Fire 2020, 29, 305. [Google Scholar] [CrossRef] [Green Version]
  3. Pereira, M.; Calado, T.; DaCamara, C.; Calheiros, T. Effects of regional climate change on rural fires in Portugal. Clim. Res. 2013, 57, 187–200. [Google Scholar] [CrossRef] [Green Version]
  4. Tedim, F.; Remelgado, R.; Borges, C.; Carvalho, S.; Martins, J. Exploring the occurrence of mega-fires in Portugal. For. Ecol. Manag. 2013, 294, 86–96. [Google Scholar] [CrossRef]
  5. Iriarte-Goñi, I.; Ayuda, M.-I. Should Forest Transition Theory include effects on forest fires? The case of Spain in the second half of the twentieth century. Land Use Policy 2018, 76, 789–797. [Google Scholar] [CrossRef] [Green Version]
  6. dos Santos, S.M.B.; Bento-Gonçalves, A.; Vieira, A. Research on Wildfires and Remote Sensing in the Last Three Decades: A Bibliometric Analysis. Forests 2021, 12, 604. [Google Scholar] [CrossRef]
  7. dos Santos, S.M.B.; Bento-Gonçalves, A.; Franca-Rocha, W.; Bapt ista, G. Assessment of Burned Forest Area Severity and Postfire Regrowth in Chapada Diamantina National Park (Bahia, Brazil) Using dNBR and RdNBR Spectral Indices. Geosciences 2020, 10, 106. [Google Scholar] [CrossRef] [Green Version]
  8. Nunes, A.; Lourenço, L.; Bento-Gonçalves, A.; Vieira, A. Três décadas de incêndios florestais em Portugal: Incidência regional e principais fatores responsáveis. Cad. Geogr. 2013, 32, 133–143. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Haines-Young, R.; Potschin, M.B. Common International Classification of Ecosystem Services (CICES) V5.1: Guidance on the Application of the Revised Structure, 1st ed.; Barton in Fabis: Nottingham, UK, 2018; Available online: https://cices.eu/content/uploads/sites/8/2018/01/Guidance-V51-01012018.pdf (accessed on 1 June 2022).
  10. Dubayah, R.; Blair, J.B.; Goetz, S.; Fatoyinbo, L.; Hansen, M.; Healey, S.; Hofton, M.; Hurtt, G.; Kellner, J.; Luthcke, S.; et al. The Global Ecosystem Dynamics Investigation: High-resolution laser ranging of the Earth’s forests and topography. Sci. Remote Sens. 2020, 1, 100002. [Google Scholar] [CrossRef]
  11. Potapov, P.; Li, X.; Hernandez-Serna, A.; Tyukavina, A.; Hansen, M.C.; Kommareddy, A.; Pickens, A.; Turubanova, S.; Tang, H.; Silva, C.E.; et al. Mapping global forest canopy height through integration of GEDI and Landsat data. Remote Sens. Environ. 2021, 253, 112165. [Google Scholar] [CrossRef]
  12. Rishmawi, K.; Huang, C.; Zhan, X. Monitoring Key Forest Structure Attributes across the Conterminous United States by Integrating GEDI LiDAR Measurements and VIIRS Data. Remote Sens. 2021, 13, 442. [Google Scholar] [CrossRef]
  13. Kamusoko, C. Vegetation and Forest Cover Applications. Ai Geolabs Publish. 2022. Available online: https://aigeolabs.com/courses/course-6/ (accessed on 1 June 2022).
  14. Kacic, P.; Hirner, A.; Ponte, E. Da Fusing Sentinel-1 and -2 to Model GEDI-Derived Vegetation Structure Characteristics in GEE for the Paraguayan Chaco. Remote Sens. 2021, 13, 5105. [Google Scholar] [CrossRef]
  15. Farr, T.G.; Rosen, P.A.; Caro, E.; Crippen, R.; Duren, R.; Hensley, S.; Kobrick, M.; Paller, M.; Rodriguez, E.; Roth, L.; et al. The Shuttle Radar Topography Mission. Rev. Geophys. 2007, 45, RG2004. [Google Scholar] [CrossRef] [Green Version]
  16. Dubayah, R.; Blair, J.B.; Beck, J.; Wirt, B.; Armston, J.; Hofton, M.; Luthcke, S.; Tang, H. GLOBAL Ecosystem Dynamics Investigation (GEDI) Level 2 User Guide; 2021; Volume 3. Available online: https://lpdaac.usgs.gov/documents/986/GEDI02_UserGuide_V2.pdf (accessed on 1 June 2022).
  17. Belgiu, M.; Drăguţ, L. Random forest in remote sensing: A review of applications and future directions. ISPRS J. Photogramm. Remote Sens. 2016, 114, 24–31. [Google Scholar] [CrossRef]
  18. Evangelides, C.; Nobajas, A. Red-Edge Normalised Difference Vegetation Index (NDVI705) from Sentinel-2 imagery to assess post-fire regeneration. Remote Sens. Appl. Soc. Environ. 2020, 17, 100283. [Google Scholar] [CrossRef]
  19. Fernández-Manso, A.; Fernández-Manso, O.; Quintano, C. SENTINEL-2A red-edge spectral indices suitability for discriminating burn severity. Int. J. Appl. Earth Obs. Geoinf. 2016, 50, 170–175. [Google Scholar] [CrossRef]
  20. White, J.D.; Ryan, K.C.; Key, C.C.; Running, S.W. Remote sensing of forest fire severity and vegetation recovery. Int. J. Wildl. Fire 1996, 6, 125–136. [Google Scholar] [CrossRef] [Green Version]
  21. Goldshleger, N.; Ben-Dor, E.; Lugassi, R.; Eshel, G. Soil Degradation Monitoring by Remote Sensing: Examples with Three Degradation Processes. Soil Sci. Soc. Am. J. 2010, 74, 1433–1445. [Google Scholar] [CrossRef]
Figure 1. Study area. Source: Organized by the authors.
Figure 1. Study area. Source: Organized by the authors.
Environsciproc 22 00050 g001
Figure 2. Workflow diagram depicting data preprocessing steps. Source: Organized by the authors.
Figure 2. Workflow diagram depicting data preprocessing steps. Source: Organized by the authors.
Environsciproc 22 00050 g002
Figure 3. Importance of the variables before and after fires. Source: Organized by the authors.
Figure 3. Importance of the variables before and after fires. Source: Organized by the authors.
Environsciproc 22 00050 g003
Figure 4. Prediction and observation are split into prefire training and validation data. Source: Organized by the authors.
Figure 4. Prediction and observation are split into prefire training and validation data. Source: Organized by the authors.
Environsciproc 22 00050 g004
Figure 5. Prediction and observation are split into training and validation data after fires. Source: Organized by the authors.
Figure 5. Prediction and observation are split into training and validation data after fires. Source: Organized by the authors.
Environsciproc 22 00050 g005
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Folharini, S.; dos Santos, S.M.B.; Bento-Gonçalves, A.; Vieira, A. Estimation of Tree Height in Burned Areas with GEDI Laser Data in Northern Portugal and Galicia (Spain). Environ. Sci. Proc. 2022, 22, 50. https://doi.org/10.3390/IECF2022-13050

AMA Style

Folharini S, dos Santos SMB, Bento-Gonçalves A, Vieira A. Estimation of Tree Height in Burned Areas with GEDI Laser Data in Northern Portugal and Galicia (Spain). Environmental Sciences Proceedings. 2022; 22(1):50. https://doi.org/10.3390/IECF2022-13050

Chicago/Turabian Style

Folharini, Saulo, Sarah Moura Batista dos Santos, António Bento-Gonçalves, and António Vieira. 2022. "Estimation of Tree Height in Burned Areas with GEDI Laser Data in Northern Portugal and Galicia (Spain)" Environmental Sciences Proceedings 22, no. 1: 50. https://doi.org/10.3390/IECF2022-13050

APA Style

Folharini, S., dos Santos, S. M. B., Bento-Gonçalves, A., & Vieira, A. (2022). Estimation of Tree Height in Burned Areas with GEDI Laser Data in Northern Portugal and Galicia (Spain). Environmental Sciences Proceedings, 22(1), 50. https://doi.org/10.3390/IECF2022-13050

Article Metrics

Back to TopTop