Otwarty dostęp

Feasibility study of flood risk monitoring based on optical satellite data


Zacytuj

Introduction

Flood hazard monitoring covers a wide range of activities: from forecasting rainfalls to monitoring flood infrastructures. The research of flood risk monitoring based on the optical satellite data was conducted as part of the SAFEDAM project – “Advanced technologies in the prevention of flood hazard”. The project is financed by the National Centre for Research and Development in Defence and Security Programme. The aim of the SAFEDAM project is to develop an innovative system enabling the risk of river floods in Poland to be monitored using, inter alia, unmanned aerial platform, satellite and aerial data. For this purpose, different photogrammetric and remote sensing techniques were tested and will be used in monitoring the flood risk, the flood phenomenon and the condition of levees. A comprehensive system will enable the multi-source collection of data, automatic data analysis and visualization for hydrological services and crisis management professionals.

Satellite Earth Observation

For the purpose of flood risk monitoring, depending on the scale on which the analysis is performed, a variety of satellite data is used. In order to carry out effective monitoring, several conditions must be met: 1. data should be collected at the appropriate frequency depending on the mode (preventive/ operational); 2. the data should be adequately detailed to provide information on the state of the water and the status of the flood embankments.

Data from meteorological systems is useful to forecast precipitation and its intensity. Very high-resolution satellite images (VHRS) and high-resolution satellite images (HRS) can be used to monitor flood embankments and their condition, and to monitor floodplains. This applies to both radar and optical data. Radar data is weather independent and generally considered to be better suited for water detection (Pluto-Kossakowska et al. 2017), while optical data is easier for photointerpretation and acts as a compaction of measurements from radar data.

The method of performing flood risk maps and flood risk monitoring has been presented by Skakun (2012) for Namibia with the use of multi-temporal LANDSAT images. Shaker et al. (2012) determined one of the best methods of using SPOT PAN satellite data to estimate the risk of flooding. Malinowski et al. (2015) conducted a study in which the total accuracy of twelve maps of flood ranges obtained by the decision tree method ranged between 77% and 95%. Asare-Kyei et al. (2015) used RapidEye satellite data to produce land cover maps in order to obtain flood risk indicator maps. An analogous approach which presents the usage of optical and radar satellite data to produce flood risk maps (where optical data is used for the implementation of land cover maps) and, based on radar data, identified the extent of the flood presented, was introduced by Stancalie et al. (2007). Byun et al. (2015) present a methodology for detection of flood changes based on VHRS KOMPSAT images from two dates without the need for user interaction. Cundill’s (2014) article presents the usage of remote sensing as a tool for the inspection of grass-covered embankments. Summarising the existing studies, various research teams use different methods and indices to identify flood risk areas depending on both the type of the area of interest and the data availability (the issue of costs and availability of data in the archive). Optical data is typically used to obtain information about the current status of land cover, while the extent of flood waters has frequently been obtained on the basis of radar data.

Optical satellite systems

Optical satellite remote sensing systems are characterized by a variety of technical parameters, including very important spatial and spectral resolution. Both of these parameters are crucial for monitoring the condition of embankments and forecasting the occurrence of floods. A detailed analysis of currently available VHRS and VHR satellite data was carried out, both in terms of their spatial resolution and their spectral resolution (Tables 1 and 2). The tables present the currently operating satellite systems as well as the systems operating in 2010, because this was a year in which there was a flood threat in Poland and data from this period can be used as part of the research to test various image processing methods in terms of flood risk monitoring.

VHRS optical satellite data parameters

Satellite/SensorIkonosEros AQuickBirdOrbView-3Eros BKOMP SAT-2CARTO SAT-2World-View-1
Nationality/CompanyGeoEye (previously Space Imaging), USAImageSat International, IsraelDigital Globe, USAOrbImage, USAImageSat International, IsraelKARI, South KoreaISRO, IndiaDigital Globe, USA
Operating timesince 24.09.1999since 05.12.200018.10.2001 – 17.12.201426.06.2003 – 03.2007since 25.04.2006since 28.07.2006since 10.01.2007since 18.09.2007
Orbit altitude [km]681-709475-491450470500685630496
Spectral bands [um]PAN0.45-0.900.50-0.900.45-0.900.45-0.900.50-0.900.50-0.900.50-0.850.45-0.80
MS0.45-0.520.45-0.520.45-0.520.45-0.52
0.52-0.610.52-0.600.52-0.600.52-0.60
0.64-0.720.63-0.690.62-0.690.63-0.69
0.72-0.880.76-0.900.76-0.900.76-0.90
Spatial resolution[m]PAN0.82×0.821.0×1.0 -1.8×1.80.61×0.611×10.7×0.71×10.8×0.80.50×0.50
MS3.28×3.282.44×2.444×41.9×1.94×4
Radiometric resolution [bit]1111111111101111
Footprint [km]11×1113.5×13.516×168×814×1415×159.6×9.617.6×17.6
Maximum angle of incidence+/- 45°+/- 43°+/- 30°+/- 45°+/- 45°+/- 60°+/- 45°+/- 45°
Satellite/ SensorGeoeye-1WorldView-2CartoSat-2BPlèiades-1AKOMPSAT-3Plèiades-1BWorld-View-3KOMPSAT-3A
Nationality/CompanyGeoEye (previously Space Imaging), USADigital Globe, USAISRO, IndiaCNES, AIRBUS Defence & Space, FranceKARI, South KoreaCNES, AIRBUS Defence & Space, FranceDigital Globe, USAKARI, South Korea
Operating timesince 06.09.2008since 08.10.2009since 12.07.2010since 17.12.2011since 17.05.2012since 02.12.2012since 13.08.2014since 25.03.2015
Orbit altitude [km]681770637694685694617528
Spectral bands[um]PAN0.45-0.900.45-0.800.45-0.850.48-0.830.45-0.900.48-0.830.45-0.800.45-0.90
MS0.45-0.520.40-0.450.43-0.550.45-0.520.43-0.550.40-0.450.45-0.52
0.52-0.600.45-0.520.49-0.610.52-0.600.49-0.610.45-0.510.52-0.60
0.63-0.700.51-0.580.60-0.720.63-0.690.60-0.720.51-0.580.63-0.69
0.76-0.900.59-0.630.75-0.950.76-0.900.75-0.950.59-0.630.76-0.90
0.63-0.690.63-0.69
0.71-0.750.71-0.75
0.77-0.900.77-0.90
0.86-1.040.86-1.04
Spatial resolution [m]PAN0.41×0.410.46×0.4610.8×0.80.5×0.50.7×0.70.5×0.50.31×0.3110.5×0.5
MS1.65×1.651.84×1.842×22.8×2.82×21.24×1.242×2
Radiometric resolution [bit]1111101214121114
Footprint [km]15×1516.4×16.49.6×9.620×2015×1520×2013.1×13.112×12
Maximum angle of incidence+/- 60°+/-45°+/-45°+/-45°+/-45°+/-45°+/-45°+/-45°

Source: own elaboration based on: Airbus Defence and Space 2018, Satellite Imaging Corporation 2018, DigitalGlobe 2018, ESA Earth Online 2018.

Table 1 presents all available VHRS optical satellites, along with their characteristics. Among these systems, WorldView-2 data was analysed, mainly due to the availability of data for the area of interest during the flooding and Plèiades data was analysed because of the daily revisit time and future use of data in the SAFEDAM system.

Table 2 presents selected HRS optical satellites, along with their characteristics, as a complement to the VHRS data. Sentinel-2 was chosen for research because of the free data policy and convenient access to the imagery database.

Selected high resolution (HRS) optical satellite data parameters

Satellite/SensorLANDSAT-7/ETM+SPOT5ALOS/PRISM/AVNIR-2SPOT6/7LDCM/OLISentinel-2ASentinel-2B
Nationality/CompanyNASA, USAAIRBUS Defence & Space (previously SpotImage), FranceJAXA, JapanAIRBUS Defence & Space, FranceNASA, USAESA, European UnionESA, European Union
Operating timesince 15.04.199904.05.2002 -31.03.2015since 24.01.2006since 12.09.2012since 11.02.2013since 23.06.2015since 07.03.2017
Orbit altitude [km]705832692832705786786
Spectral bands [um]PAN0.52-0.900.48-0.710.52-0.770.455-0.7450.500-0.6800.430-0.4570.420-0.465
MS0.45-0.520.50-0.590.42-0.500.455-0.5250.433-0.4530.448-0.5460.443-0.541
0.52-0.600.61-0.680.52-0.600.530-0.5900.450-0.5150.538-0.5830.536-0.582
0.63-0.690.78-0.890.61-0.690.625-0.6950.525-0.6000.646-0.6840.646-0.685
0.76-0.901.58-1.750.76-0.890.760-0.8900.630-0.6800.694-0.7130.694-0.714
1.55-1.750.845-0.8850.731-0.7490.730-0.748
10.40-12.501.560-1.6600.769-0.7970.766-0.794
2.08-2.352.100-2.3000.763-0.9080.767-0.900
1.360-1.3900.848-0.8810.848-0.880
0.932-0.9580.930-0.957
1.336-1.4111.339-1.415
1.542-1.6851.540-1.681
2.081-2.3232.067-2.305
Spatial resolution [m]PAN15×152.5×2.5, 5.0×5.02.5×2.51.5×1.515×1510×1010×10
MS30×3010×1010×106.0×6.030×3010×10 20×20 60×6010×10 20×20 60×60
Radiometric resolution [bit]88812121212
Footprint [km]185×18560×6070×7060×60185×185290290
Maximum angle of incidence+/- 31°+/- 27°+/-46°+/-46°

Source: own elaboration based on: Airbus Defence and Space 2018, Satellite Imaging Corporation 2018, ESA Earth Online 2018.

The temporal resolution of the satellite systems means the revisit time and the possibility of acquiring new imageries. For the area of medium latitudes, this ranges from 1-2 days (Plèiades and SPOT constellation) up to several days (e.g. LANDSAT, Sentinel). It depends on the construction of the satellite system (a single satellite or a constellation, a scanner operating only in a nadir or a scanner with the possibility of tilting along and across the direction of flight), and the possibilities of its programming. In the case of individual satellite systems with HVRS, where programming is possible, it is usually a few days (3-4 days). For preventive monitoring this is a sufficient revisit time, but the possibilities of the usage of satellite optical systems are limited by the weather conditions and cloud cover. The experience so far shows that it is generally possible to obtain optical data every 7-10 days. However, there are periods when obtaining optical data is impossible due to unfavourable weather conditions.

Capabilities and restrictions of usage of earth observation

Optical data, due to frequent cloudiness occurring in the period of increased flood hazard, is less commonfor operational use in monitoring damage to levees. It is possible to use such data only in cloudless conditions or in the presence of slight cloud cover. This means that the identification of damage to flood embankments with the use of optical data can be carried out as quasi-continuous monitoring, conducted in order to detect damages early. Examples of optical data with very high spatial resolution (IKONOS, QuickBird) being used to identify damage to flood embankments are shown in the literature - these cover, for example, damage to the sod in the main body of the embankment and on the embankment slopes (Hossain et al. 2004). Optical data is much more frequently used to assess the state of damage after a flood event or to determine the extent of flooding, or for flood risk assessment (Stancalie et al. 2007; Ho et al. 2010; Skakun 2012; Shaker et al. 2012; Asare-Kyei et al. 2015; Byun et al. 2015; Hu et al. 2015; Malinowski et al. 2015).

In order to clarify the nomenclature and build a coherent vocabulary, the following definitions are provided: ‘identification’, as the possibility of indicating the location of the damage, but without specifying its type, and ‘recognition’, as providing the type of damage, which seems quite a limited task in the context of satellite images, especially those with high resolution. As a result of the analysis of the type and size of the potential damage to flood embankments and the available satellite data (Tables 1,2), the possibilities of using individual types of optical data to identify damage to flood embankments have been indicated (Tables 3, 4, 5). Most often with this analysis, it is possible to point out a potential place of damage, but without indicating the cause. Satellite data is only used to monitor and indicate places for more accurate measurements. Supplementary measurements based on more precise technologies are recommended.

Damage identification based on VHRS 0.3 m × 0.3 m pixel optical data

Type of damageDamageDimensions (approximate)Spectral rangesEstimated accuracy of results
Damage/loss of groundDen, cavity, tunnel caused by the activity of burrowing animals0.1 - 1 mVIS, NIR0.3 - 0.6 m
Surface damageTraces in the crown of the levees0.1 - 0.3 mVIS, NIR0.15 - 0.30 m

Source: Own elaboration

Damage identification based on VHRS 0.5 m × 0.5 m optical data

Type of damageDamageDimensions(approximate)Spectral rangesEstimated accuracy of results
Surface damageWild crossings0.3 - 1 mVIS, NIR0.3 - 0.5 m
Surface damageDefects in the body and on the sides of the levees caused by wild and burrowing animalsUnspecifiedVIS, NIR0.3 - 0.5 m
Slope deformationLocal soil runoffUnspecifiedVIS, NIR0.3 - 0.5 m

Source: Own elaboration

Levee monitoring and identification based on VHRS optical data

Using optical data with a maximum spatial resolution of 0.3 m × 0.3 m it is possible to detect damages of size equal to or bigger than 0.3 m × 0.3 m. Table 3 presents the estimated accuracy of results based on the use of airborne imagery or satellite data with the highest available resolution, such as WorldView-3. It is possible to see a smaller artefact on a VHRS image, but it is not feasible to determine details such as type of damage.

Table 4 presents the estimated accuracy of results based on the use of optical VHRS data with a 0.5 m × 0.5 m pixel size. An example of satellite sensors are Plèiades satellites.

Figure 1 presents damaged flood banks caused by animals. The holes were caused by the beavers living in the area. Grazing animals on the levees and in their immediate vicinity violates the construction of levees and causes changes in soil mass.

Figure 1

Presents damaged flood banks caused by animals. The holes were caused by the beavers living in the area. Grazing animals on the levees and in their immediate vicinity violates the construction of levees and causes changes in soil mass. Photo by Janusz Orzepowski

Damage identification based on HRS optical data

Table 5 presents the estimated accuracy of results for optical data with a 2 m × 2 m pixel size. It is possible to use SPOT 6 and 7 data to identify damage with approximate dimensions of 2-5 metres and higher. Large-scale damage, with dimensions of several or several dozen metres can be detected on the basis of optical data with a 10 m × 10 m pixel, such as Sentinel-2.

Damage identification based on HRS 2 m × 2 m optical data

Type of damageDamageDimensions (approximate)Spectral rangesEstimated accuracy of results
Surface defectsTheft of part of the infrastructure2 - 10 mVIS, NIR1 - 2 m
Surface damageTranshumance within the body of the levee2 - 10 mVIS, NIR1 - 2 m
LandslidesPossibility of landslide in the section with very low geotechnical parametersfrom several to dozens of metresVIS, NIR1 - 2 m
LandslidesUndercutting of the base of the levee by the current of the river during floodsfrom several to dozens of metresVIS, NIR1 - 2 m
Cracks in the bodyLongitudinal slits5 - 15 mVIS, NIR1 - 2 m
Blur of the shaftMay occur when the flood wave elevation is greater than the lowest point of the levee crownUnspecifiedNIR1 - 2 m

Source: Own elaboration

Figure 2 presents an example of a landslide in Poland in 2015 and levee surface damage caused by people. Wheel rims are clearly visible, probably from cars, quads or agricultural machinery. The levee must maintain the correct geometry. Changes in the soil mass cause shifts in the geometry and changes in shaft stability. In the case of traces, water is collected in them and the filtration is disturbed, which reduces the stability of the shaft.

Figure 2

Presents an example of a landslide in Poland in 2015 and levee surface damage caused by people. Photo by Janusz Orzepowski

Analysis of the possibilities of using optical data in flood risk monitoring
Purpose of the research

The main goal of the research is to examine the potential and effectiveness of the optical satellite data in flood risk monitoring. As input data, Plèiades (10.06.2010) and WorldView-2 (20.06.2013) data was used. Sentinel-2 imageries had also been taken into account earlier, but they were excluded from the research due to low spatial resolution. As mentioned above, this data is more suitable for large area changes. The concept of SAR and optical Sentinel data use for flood risk monitoring in Poland, prepared as part of the SAFEDAM project, was presented by Pluto-Kossakowska et al. (2017).

Research methodology

The research methodology is based on bare soil and water body detection using optical satellite imagery. Therefore, vegetation and soil indices utilizing some of the image bands were used. Furthermore, the bare soil and water body areas can be converted to vector layers and compared with other data (Airborne Laser Scanning, reference water boundaries, range of water for warning state, range of water for emergency state, etc.), which is useful in the analysis of data acquired with different sensors and over different time periods, as is the case in the SAFEDAM project (Figure 3).

Figure 3

Workflow of optical satellite imagery analysis

Source: Own elaboration

In the presented approach two analyses were conducted. The first analysis consisted of bare soil detection (places of potential levee damage), based on the Soil-adjusted Vegetation Index (SAVI), the Optimized Soil-adjusted Vegetation Index (OSAVI), the Optimized Soil-adjusted Vegetation Index using green band (OSAVI-G), and the Modified Soil-adjusted Vegetation Index (MSAVI).

The second analysis consisted of water range detection, based on the Normalized Differential Vegetation Index (NDVI), the Normalized Difference Water Index (NDWI), the Soil-adjusted Vegetation Index (SAVI), the Simple Ratio (SR), the Wetness (TCA) component, the Differential Vegetation Index (DVI), the Atmospherically Resistant Vegetation Index (ARVI), the Index of Free Water (IFW), the Water Index (WI), and the Water Impoundment Index (WII).

Results and evaluation

The aim of the conducted analysis was to identify indices to delimit the bare soil, as a signal that the structure of the shaft and the water range may have changed unambiguously during a flood, and to determine the threshold values at which the best and reproducible results are obtained. The analysis was conducted for optical satellite imagery which is used in the prevention of flood hazards. In the SAFEDAM project, the detection of bare soil and water is conducted automatically; therefore, it is important to use reliable indices.

Figure 4 illustrates an example of bare soil detection based on Plèiades data. The best results in bare soil detection were obtained using the MSAVI and NDVI indices. The MSAVI threshold was 0.065, while the NDVI values for bare soil ranged from -0.1 (wet soil) to 0.1 (dry soil). As can be seen, some parts of the bare soil are overestimated, which may be caused by the shadows and low greenness level of the vegetation.

Figure 4

Bare soil detection based on Plèiades data (a) using MSAVI (b) and NDVI (c)

Source: Own elaboration based on Pleiades data

Figure 5 shows the results of water body detection. One of the problems which may occur during water body extraction involves elements which are beyond the river, which may be interpreted as water – for example, areas of shadow or building roofs. Therefore, it is necessary to filter the extracted water and include only the real water body in further flood risk analysis. The NDVI, OSAVI and SAVI indices proved to be the most stable when detecting water [Table 6]. Using these indices, it is possible to find a threshold value that will enable the automatic detection of water bodies. In some cases, NDWI and ARVI give very good results, but the threshold for all analysed data (Sentinel-2, WorldView-2 and Plèiades) could not be unambiguously determined. Regarding the WI index, it was determined that it significantly overestimates the water range by including shaded areas. As a result of the conducted research, it was found that the most effective at delimiting water was the NDVI index, for which a threshold value of -0.22 gives reproducible results (for cloudless images). Using a threshold value which is too high results in a larger error of commission (Table 6) – built-up areas (roofs covered with roofing paper) can be incorrectly indicated as water areas. In the presence of clouds (Cirrus, Cirrostratus, Cirrocumulus or Altostratus) on the satellite image, a two-stage operation is necessary: 1. cloud detection (for Sentinel-2 data, the cloud masks provided by ESA can be used), 2. water detection with a NDVI threshold value between -0.13 and -0.10.

Figure 5

Detection of water bodies based on WorldView-2 (a) data using NDVI (b), SAVI (c), SARVI (d), WV-WI (e) and NDWI (f) - results comparison with reference data (a)

Source: Own elaboration based on WorldView-2 data

Assessment of the accuracy of determining the range of water bodies on the basis of index images obtained from WorldView-2 data (partly cloudy)

IndexThreshold valuesCorrectly detected waterError of omissionError of commission
DVI<-3071.61%28.39%2.93%
NDVI<-0.1494.76%5.24%10.04%
NDVI<-0.1594.40%5.60%7.54%
NDVI<-0.1693.99%6.01%5.68%
NDVI<-0.2289.73%10.27%1.40%
NDWI>0.3278.57%21.43%6.24%
SAVI<-0.2094.96%5.04%12.01%
SAVI<-0.2593.66%6.43%4.71%
SAVI<-0.2792.97%7.03%3.49%
OSAVI<-0.2590.22%9.78%1.56%
SARVI<-0.4, 0.0>88.48%11.52%15.01%
ARVI<-0.185.06%4.37%14.94%
WV-WI-0.593.36%6.64%19.88%
WV-WII<2991.26%8.74%16.49%

Source: Own elaboration

Figure 6 presents the results of Plèiades data processing. The extracted water body area covers the water presented on the satellite scene properly.

Figure 6

Detection of water bodies based on Plèiades data using NDVI

Source: Own elaboration based on Pleiades dat

The analysis shows that it is clearly possible to detect damage in turf on levees of not less than 0.3 m. Minor changes of a few or a few dozen centimetres are impossible to recognise. Only significant changes can be detected using satellite imagery with a resolution of 5 m × 5 m or more. Examples of such damage may be those related to human activity, with approximate dimensions of 10-15 metres. Sentinel-2 can be used to identify large area changes. The automatization of satellite data processing (from different satellite systems) gives opportunities for hydrological services and crisis management professionals to access the accurate and actual information about the condition of levees for the early detection of flood risk.

The initial results of the studies (using vegetation, soil and water indices based on satellite data) demonstrate the potential of very high and high-resolution satellite imagery for the operational detection of damage and risk areas:

minor changes of a few and a few dozen cm are virtually impossible to recognise. There must be special conditions (e.g. high contrast with the surrounding terrain) so that they can be seen in satellite imagery;

changes below 1 m can be detected using VHR images such as WorldView-3, WorldView-2, GeoEye, Plèiades;

changes from 1 to several metres can be detected using SPOT 6 and 7 images;

only significant changes with approximate dimensions of 10-15 metres can be detected using satellite imagery with a resolution of 5 m × 5 m or more, such as Sentinel-2 data.

Discussions and conclusions

The main goal of the SAFEDAM project is to create an IT system for managing spatial data related to levees. The system will enable data collection with a non-invasive method of data acquisition (UAV and satellite data), automatic data analysis, visualization for hydrological services and for society, and provide tools for crisis management professionals and rescue teams during the intervention phase.

Earth observations are one of the key elements of the system, providing continuous information about a large area of interest at the same time. Analysis based on satellite data is used in both preventive and intervention mode. Satellite data can be used for large-scale image analysis to pre-select areas dedicated to unmanned platform measurement. Using processing and the observation of time series, the multi-temporal analysis tool may indicate simple analyses based on vegetation indices, to identify landslide areas, the occurrence of illegal objects, anthropopressure and the impact of animals.

As a result of the conducted analysis it was found that NDVI is the most suitable index for the automatic determination of a water range with a threshold value of -0.22 (for cloudless satellite scenes). This threshold value makes it possible to detect water with a slight overestimation error of 1-1.5% of the reference water mask area. In Rokni et al. (2014), where Landsat imagery was used for water surface detection, similar NDVI values were acquired (-0.21 and -0.29). However, the detection accuracy in Rokni et al. (2014) was slightly higher, which may have been caused by the selection of a more homogeneous, water-covered area than the test area selected for this study (Fig. 5). Moreover, whereas the NDWI and modification of this index gives accurate results in many studies (El Asmar et al. 2013; Rokni et al. 2014), in the analysis presented above the NDWI-based water detection was not satisfactory. Additionally, according to the proposed approach, in the event that high-level clouds (e.g. Cirrus) or middle-level clouds (e.g. Altus) occur over the analysed area, it is necessary to apply a two-stage processing workflow and two different threshold values for cloudless and cloudy parts of the scene.

Damage and changes on the surface of the embankments are detectable, but it is difficult to clearly determine the type and cause. The presented research has shown that damage visible as bare soil can be detected using the MSAVI or NDVI indices. However, if the threshold value is too high, the bare soil is overestimated.

This feasibility study will support the SAFEDAM System by providing a clear overview of the existing optical satellite systems and their capabilities, as well as taking into account their limitations.

eISSN:
2084-6118
Język:
Angielski
Częstotliwość wydawania:
4 razy w roku
Dziedziny czasopisma:
Geosciences, Geography, other