5 p.Nowadays Earth observation satellites, in particular Landsat, provide a valuable help to forest managers in
post-fire operations; being the base of post-fire damage maps that enable to analyze fire impacts and to develop
vegetation recovery plans. Sentinel-2A MultiSpectral Instrument (MSI) records data in similar spectral wavelengths
that Landsat 8 Operational Land Imager (OLI), and has higher spatial and temporal resolutions. This
work compares two types of satellite-based maps for evaluating fire damage in a large wildfire (around 8000 ha)
located in Sierra de Gata (central-western Spain) on 6–11 August 2015. 1) burn severity maps based exclusively
on Landsat data; specifically, on differenced Normalized Burn Ratio (dNBR) and on its relative versions (Relative
dNBR, RdNBR, and Relativized Burn Ratio, RBR) and 2) burn severity maps based on the same indexes but
combining pre-fire data from Landsat 8 OLI with post-fire data from Sentinel-2A MSI data. Combination of both
Landsat and Sentinel-2 data might reduce the time elapsed since forest fire to the availability of an initial fire
damage map. Interpretation of ortho-photograph Pléiades 1 B data (1:10,000) provided us the ground reference
data to measure the accuracy of both burn severity maps. Results showed that Landsat based burn severity maps
presented an adequate assessment of the damage grade (κ statistic = 0.80) and its spatial distribution in wildfire
emergency response. Further using both Landsat and Sentinel-2 MSI data the accuracy of burn severity maps,
though slightly lower (κ statistic = 0.70) showed an adequate level for be used by forest managersS

Fernández Manso, Alfonso

Fernández-Manso, Óscar

Quintano Pastor, Carmen

English

Leon University (Spain)

UNCORRECTEDPROOFInternational Journal of Applied Earth Observations and Geoinformation xxx (2017) xxx-xxxContents lists available at ScienceDirectInternational Journal of Applied Earth Observations andGeoinformationjournal homepage: www.elsevier.comCombination of Landsat and Sentinel-2 MSI data for initial assessing of burn severityC. Quintano ⁠a⁠, ⁠⁎, A. Fernández-Manso ⁠b, O. Fernández-Manso ⁠ca Electronic Technology Department, Sustainable Forest Management Research Institute, University of Valladolid-INIA, C/Francisco Mendizábal, s/n, 47014-Valladolid, Spainb Agrarian Science and Engineering Department, University of León,Av. Astorga s/n, 24400 Ponferrada, Spainc Civil Protection Agency, Castilla y León Government, Valladolid, SpainA R T I C L E I N F OKeywords:Sentinel-2LandsatdNBRBurn severityMediterranean ecosystemsA B S T R A C TNowadays Earth observation satellites, in particular Landsat, provide a valuable help to forest managers inpost-fire operations; being the base of post-fire damage maps that enable to analyze fire impacts and to developvegetation recovery plans. Sentinel-2A MultiSpectral Instrument (MSI) records data in similar spectral wave-lengths that Landsat 8 Operational Land Imager (OLI), and has higher spatial and temporal resolutions. Thiswork compares two types of satellite-based maps for evaluating fire damage in a large wildfire (around 8000 ha)located in Sierra de Gata (central-western Spain) on 6–11 August 2015. 1) burn severity maps based exclusivelyon Landsat data; specifically, on differenced Normalized Burn Ratio (dNBR) and on its relative versions (Rela-tive dNBR, RdNBR, and Relativized Burn Ratio, RBR) and 2) burn severity maps based on the same indexes butcombining pre-fire data from Landsat 8 OLI with post-fire data from Sentinel-2A MSI data. Combination of bothLandsat and Sentinel-2 data might reduce the time elapsed since forest fire to the availability of an initial firedamage map. Interpretation of ortho-photograph Pléiades 1 B data (1:10,000) provided us the ground referencedata to measure the accuracy of both burn severity maps. Results showed that Landsat based burn severity mapspresented an adequate assessment of the damage grade (κ statistic = 0.80) and its spatial distribution in wildfireemergency response. Further using both Landsat and Sentinel-2 MSI data the accuracy of burn severity maps,though slightly lower (κ statistic = 0.70) showed an adequate level for be used by forest managers.1. IntroductionForest Mediterranean ecosystems are greatly affected by fires(San-Miguel-Ayanz et al., 2016). Accurate knowledge of fire-damagedareas is fundamental for fire management, planning and monitoringvegetation restoration (Brewer et al., 2005). Satellite-based data, partic-ularly Landsat data, is becoming key information to map damaged area(both burned area and burn severity level) accurately and quasi-imme-diately after fire (Chen et al., 2015; Fang and Yang, 2014; Quintanoet al., 2013). Sentinel-2 MultiSpectral Instrument (MSI) records datain similar spectral wavelengths that Landsat 8 Operational Land Im-ager (OLI), has a higher spatial resolution (10/20 m vs 30 m) and ahigher temporal resolution (5 days revisited −considering both Sen-tinel-2A and Sentinel-2B- vs 16 days). Fernandez-Manso et al. (2016)and Navarro et al. (2017), among others, have evaluated burn sever-ity based on Sentinel-2 MSI data successfully. Additionally, recentstudies (Shoko and Mutanga, 2017; van der Werff and van der Meer,2016) have showed the suitability and even superiority of Sentinel-2MSI data in natural resources applications.Threshold-based classification of Normalized Burn Ratio difference(dNBR) (Key and Benson, 2006) has turned into a methodological ref-erence to obtain burn severity maps (Soverel et al., 2010). However,scene selection for the multitemporal approach is not always easy. Ide-ally, scene pairs should represent similar phenology and moisture andthey should not exhibit between-scene land-cover changes. Addition-ally, they should be cloud-free, and they should not contain the firescar near the scene edge (Key and Benson, 2006). The aim of thisstudy is to validate the combined used of Sentinel-2 MSI and Landsat8 OLI data to obtain a burn severity map based on dNBR, or on itsrelative versions. In this way the process of finding the scenes couldbe greatly reduced compared to using just Landsat data, and an initialfire damage map would be available almost immediately after fire. We⁎ Corresponding author.Email addresses: menchu@tele.uva.es (C. Quintano); alfonso.manso@unileon.es (A. Fernández-Manso); out-fermanos@jcyl.es (O. Fernández-Manso)https://doi.org/10.1016/j.jag.2017.09.014Received 10 July 2017; Received in revised form 22 September 2017; Accepted 25 September 2017Available online xxx0303-2434/ © 2017.Short communicationUNCORRECTEDPROOFC. Quintano et al. International Journal of Applied Earth Observations and Geoinformation xxx (2017) xxx-xxxwill compare the accuracy of two types of burn severity maps: basedexclusively on Landsat data (as usual) and based on both Landsat andSentinel-2 MSI data. This is the first Landsat/Sentinel-2 MSI based studythat assesses burn damage in Mediterranean ecosystems.2. MaterialsThe study area was located in Sierra de Gata (central-western Spain)where a wildfire burned 79.50 km⁠2 from 6 to 11 August 2015 (Fig. 1).Climate is Mediterranean. Affected vegetation was mainly a combina-tion of shrubland and forest dominated by Pinus pinaster Ait. and Quer-cus pyrenaica Wild.Burn severity maps based exclusively on Landsat data were ob-tained from two USGS Landsat 8 OLI surface reflectance images (path/row 203/32, higher-level data product on demand), acquired on 6August 2015 and on 7 September 2015 respectively. Additionally, apost-fire Sentinel-2A MSI image acquired on 21 August 2015 (tile num-ber T29TPE, level-1C) was downloaded from www.earthex-plorer.usgs.gov.As in most research studies and official burn severity maps, threeburn severity levels were considered: high, moderate and low. Inter-pretation of ortho-photograph Pléiades 1 B data (1:10,000) acquiredon 15 August 2015 enabled us to obtain 305 photo-interpreted plotsused to measure the accuracy of the burn severity maps: 90 unburned, 60 low burnseverity, 75 moderate burn severity and 90 high burn severity. The ac-quisition times of Pléiades 1 B (11:39 UTC), Landsat 8 OLI (11:07 UTC)and Sentinel-2A MSI data (11:21 UTC) are quite similar, and the smalldifferences did not impact on the accuracy measurement process of burnseverity maps.3. MethodsBoth Landsat 8 OLI surface reflectance images were subset to theselected forest fire and pre- and post-fire NBR images were calculated(NBR_pre_L, NBR_post_L, respectively). Next, we computed dNBR, Rela-tive dNBR (RdNBR) (Miler and Thode, 2007) and Relativized Burn Ratio(RBR) (Parks et al., 2014) and denoted them respectively by dNBR_L,RdNBR_L, and RBR_LThe level-1C Sentinel-2A MSI image (Top of Atmosphere, TOA, Re-flectance) was converted to surface reflectance image (or Botton of At-mospehre, BOA, Reflectance) by using Sen2Cor tool provided by the Eu-ropean Space Agency (ESA). The surface reflectance image was subset tothe selected fire and NBR was computed. We denoted it by NBR_post_S.Regridding was applied to match the spatial resolution of Sen-tinel-2A MSI image to Landsat 8 OLI data. Storey et al. (2016) alertedof a temporary misalignment between Sentinel-2A MSI image to Land-sat 8 OLI data, thatFig. 1. Location of study area.2UNCORRECTEDPROOFC. Quintano et al. International Journal of Applied Earth Observations and Geoinformation xxx (2017) xxx-xxxvary geographically. However, no misalignment was observed betweenthese two types of remote sensed data in our study area. Finally, com-bining NBR_pre_L and NBR_post_S, dNBR, RdNBR, and RBR images wereobtained as well. We denoted them respectively by dNBR_LS, RdNBR_LS,and RBR_LS.Burn severity maps from dNBR, RdNBR and RBR images were ob-tained by applying the threshold levels recommended by Botella andFernández-Manso (in press), who adapted the methodology proposed byKey and Benson (2006) to Spain. Error matrixes of burn severity esti-mates were calculated from the 305 photo-interpreted sampling plots.Overall accuracy (OA), producer’s accuracy (PA), user’s accuracy (UA)and κ statistic were computed as well (Congalton and Green, 2009). AZ-test enabled us to check whether κ statistic differences between thetwo burn severity maps had statistical significance.Fig. 2. Upper: difference normalized burned ratio, dNBR; left: from Landat data (dNBR_L8), right: from Landsat and Sentinel data (dNBR_L8_S2A); lower: burn severity estimates; left:from dNBR_L8; right: from dNBR_L8_S2A.3UNCORRECTEDPROOFTable 1Accuracy measurement of burn severity maps based on dNBR_L, dNBR_LS, RdNBR_L, RdNBR_LS, NBR_L, and NBR_LS.dNBR_L dNBR_LS RdNBR_L RdNBR_LS RBR_L RBR_LSPA Total 0.77 0.68 0.81 0.74 0.81 0.73U 0.93 1.00 0.95 0.93 0.94 0.94L 0.66 0.50 1.00 0.43 0.82 0.40M 0.68 0.48 0.76 0.58 0.76 0.58H 0.82 0.82 0.77 1.00 0.77 1.00UA Total 0.77 0.70 0.85 0.74 0.82 0.73U 0.77 0.77 1.00 0.77 0.90 0.75L 0.77 0.40 0.77 0.30 0.77 0.30M 0.69 0.69 0.69 1.00 0.69 1.00H 0.87 0.87 0.87 0.87 0.87 0.87OA 0.79 0.72 0.84 0.77 0.83 0.76κ statistic 0.71 0.62 0.80 0.69 0.76 0.70σκ 0.001 0.001 0.001 0.001 0.001 0.001Z 1.990 2.562 1.455PA: producer’s accuracy; UA: user’s accuracy; OA: overall accuracy; σκ: variance of κ statistic; dNBR_L: Normalized Burn Ratio difference based exclusively on Landsat 8 OLI data; dNBR_LS: Normalized Burn Ratio difference based on both Landsat 8 OLI andSentinel-2 MSI data; RdNBR_L: Relative Normalized Burn Ratio difference based exclusively on Landsat 8 OLI data; RdNBR_LS: Relative Normalized Burn Ratio difference based on both Landsat 8 OLI and Sentinel-2 MSI data; RBR_L: Relativized Burn Ratiobased exclusively on Landsat 8 OLI data; RBR_LS: Relativized Burn Ratio based on both Landsat 8 OLI and Sentinel-2 MSI data; U: unburned; L: low bur severity; M: moderate burn severity; H: high burn severity.Bold: burn severity maps based on Landsat data exclusively and based on Landsat and Sentinel-2 MSI data have significant differences.UNCORRECTEDPROOFC. Quintano et al. International Journal of Applied Earth Observations and Geoinformation xxx (2017) xxx-xxx4. ResultsFig. 2 (upper) displays dNBR_L8 and dNBR_L8_LS. Visually both im-ages are quite similar. However, we can appreciate some differences be-tween the burn severity maps obtained from them (Fig. 2, lower). Themap based exclusively on Landsat data shows some unburned areas in-side the fire perimeter that were labeled as low burn severity level onthe burn severity map based on dNBR_L8_LS image.Accuracy measurements are shown on Table 1. Accuracy of burnseverity maps based on the combination of Landsat and Sentinel-2 datareached the values of 0.62, 0.69 and 0.70, respectively, when dNBR_LS,RdNBR_LS and RBR_LS were classified. Burn severity maps based exclu-sively on Landsat data had a κ statistic higher than burn severity mapsbased on both Landsat and Sentinel-2 data (0.71, 0.80 and 0.76). How-ever, there were not significant differences between the bun severitymaps based on RBR_L and RBR_LS (Z <1.96). Burn severity maps basedon the relative versions of dNBR had higher accuracy than maps basedon dNBR. Particularly, the highest accuracy was reached when RdNBR_Limage were classified (κ = 0.80), followed by RBR indexes (0.76, 0.70,there was not significant difference between these two κ values). PA val-ues for low and moderate burn severity levels were higher when burnseverity maps based exclusively on Landsat data were taken into ac-count. Similarly, UA values for low burn severity level are lower forburn severity maps based on both Landsat 8 OLI and Sentinel-2A MSIdata. These differences show a higher degree of misclassification at lowthan at high burn severity levels.5. DiscussionThe differences between the burn severity maps displayed in Fig. 2,particularly at low severity and unburned classes, point to some misclas-sification problems at low severity levels, that are showed in Table 1 aswell. Misclassification problems at low severity pixels are common andhave been reported by different studies (Cocke et al., 2005; Sundermanand Weisberg, 2011; Stambaugh et al., 2015). Among the main possiblecauses of the visual differences found between the burn severity mapson this study we might mention: 1) reflectance calibration differencesbetween Landsat 8 OLI and Sentinel-2 MSI data, and 2) the thresholdlevels used to obtain the burn severity estimation based on both Land-sat and Sentinel-2 data should be carefully tested, though they showedvery promising results (Botella and Fernández-Manso 2017)(Botella andFernández-Manso, in press). Additional research is necessary to performsuch threshold revision. We leave it for future work.From a forest manager's perspective the most common reason for es-timating burn severity is, however, to identify areas that need to be re-covered (Cocke et al., 2005). For that reason, discriminating accuratelythe high burn severity level from the rest of levels may supply the keyinformation. Both burn severity maps are quite similar when taking intoaccount the high burn severity level. However, the burn severity mapbased on the combination of Landsat and Sentinel-2 data may be oper-ative as a preliminary emergency map, as it allows an accurate initialapproximation of fire damage level. Emergency maps help forest man-agers to minimize the potentially harmful effects of excessive post-fireerosion (Miller and Yool, 2002).6. ConclusionWe developed an operative method based on the combination ofLandsat 8 OLI and Sentinel-2 MSI data to obtain emergency burn sever-ity maps thatallow accurate initial approximation of fire damage level. This workproved that in case of unavailability of a pre-fire or a post-fire Land-sat image, one of them could be substituted for a Sentinel-2 MSI imageto obtain an initial burn severity map. In our study, the burn severitymap based exclusively on Landsat data showed higher accuracy than theburn severity map based on both Landsat and Sentinel-2 data, thoughmaps based on both sensors reached a more than acceptable level of ac-curacy (κ = 0.70). Future research should adjust the used classificationthresholds.AcknowledgmentsThe research work was financially supported by the Spanish Ministryof Economy and Competitiveness, and the European Regional Develop-ment Fund (ERDF) in the frame of the GES-FIRE project “Multi-scaletools for the post-fire management of fire prone ecosystems in the con-text of global change” (AGL201348189-C2-1-R), and by the RegionalGovernment of Castilla y León within the frame of the FIRECYL project(LE033U14).ReferencesBotella, M.A., Fernández-Manso, A., 2017. Study of post-fire severity in the Valencia re-gion comparing the NBR, RdNBR and RBR indexes derived from Landsat 8 images.Revista de Teledetección)in press.Brewer, C.K., Winne, J.C., Redmond, R.L., Opitz, D.W., Mangrich, M.V., 2005. Classifyingand mapping wildfire severity: a comparison of methods. Photogramm. Eng. RemoteSens. 71, 1311–1320.Chen, G., Metz, M.R., Rizzo, D.M., Meentemeyer, R.K., 2015. Mapping burn severity in adisease-impacted forest landscape using Landsat and MASTER imagery. Int. J. Appl.Earth Obs. Geoinf. 40, 91–99.Cocke, A.E., Fulé, P.Z., Crouse, J.E., 2005. Comparison of burn severity assessments usingdifferenced normalized burn ratio and ground data. Int. J. Wildland Fire 14, 189–198.Congalton, R.G., Green, K., 2009. Assessing the Accuracy of Remotely Sensed Data: Prin-ciples and Practice. CRC press, Boca Raton.Fang, L., Yang, J., 2014. Atmospheric effects on the performance and threshold extrapo-lation of multi-temporal Landsat derived dNBR for burn severity assessment. Int. J.Appl. Earth Obs. Geoinf. 33, 10–20.Fernandez-Manso, A., Fernandez-Manso, O., Quintano, C., 2016. SENTINEL-2A red-edgespectral indices suitability for discrimination burn severity. Int. J. Appl. Earth Obs.Geoinf. 50, 170–175.Miler, J.D., Thode, A.E., 2007. Quantifying burn severity in a heterogeneous landscapewith a relative version of the delta Normalized Burn Ratio (dNBR). Remote Sens. En-viron. 109, 66–80.Miller, J.D., Yool, S.R., 2002. Mapping forest post-fire canopy consumption in severaloverstory types using multi-temporal Landsat TM and ETM data. Remote Sens. Envi-ron. 82, 481–496.Navarro, G., Caballero, I., Silva, G., Parra, P., Vázquez, A., Caldeira, R., 2017. Evaluationof forest fire on Madeira Island using Sentinel-2A MSI imagery. Int. J. Appl. Earth Obs.Geoinf. 58, 97–106.Parks, S.A., Dillon, G.K., Miller, C., 2014. A new metric for quantifying burn severity: therelativized burn ratio. Remote Sensing 6, 1827–1844.Quintano, C., Fernández-Manso, A., Roberts, D.A., 2013. Multiple Endmember SpectralMixture Analysis (MESMA) to map burn severity levels from Landsat images inMediterranean countries. Remote Sens. Environ. 136, 76–88.San-Miguel-Ayanz, J., Durrant, T., Boca, R., Libertá, G., Boccacci, F., Di-Leo, M.,López-Pérez, J., Schulte, E., 2016. Forest Fires in Europe, Middle East and NorthAfrica 2015 . Joint Research Centre Technical Reports, EUR 28148 EN (ISBN978–92-79–62959-4).Shoko, C., Mutanga, O., 2017. Seasonal discrimination of C3 and C4 grasses functionaltypes: an evaluation of the prospects of varying spectral configurations of new gener-ation sensors. Int. J. Appl. Earth Obs. Geoinf. 62, 47–56.Soverel, N.O., Perrakis, D.B., Coops, N.C., 2010. Estimating burn severity from Land-sat dNBR and RdNBR indices across western Canada. Remote Sens. Environ. 114,1896–1909.Stambaugh, M.C., Hammer, L.D., Godfrey, R., 2015. Performance of burn-Severity metricsand classification in oak Woodlands and grasslands. Remote Sens. 7, 10501–10522.Storey, J., Roy, D.P., Masek, J., Gascond, F., Dwyer, J., Choate, M., 2016. A note on thetemporary misregistration of landsat-8 operational land imager (OLI) and sentinel-2multi spectral instrument (MSI) imagery. Remote Sens. Environ. 186, 121–122.Sunderman, S., Weisberg, P., 2011. Remote sensing approaches for reconstructing fireperimeters and burs severity mosaics in desert spring ecosystems. Remote Sens. Envi-ron. 115, 2384–2389.van der Werff, H., van der Meer, F., 2016. Sentinel-2A MSI and landsat 8 OLI provide datacontinuity for geological remote sensing. Remote Sens. 8, 883.5

Combination of Landsat and Sentinel-2 MSI data for initial assessing of burn severity

https://buleria.unileon.es/bitstream/10612/6810/1/International%20Journal%20of%20Applied%20Earth.pdf

Combination of Landsat and Sentinel-2 MSI data for initial assessing of burn severity

Abstract

Similar works

Full text

Available Versions

Leon University (Spain)