Flood Disaster Mapping Using Geospatial Techniques: A Case Study of the 2022 Pakistan Floods †

: Remote sensing images are an essential tool for mapping the amount of ﬂood inundation after ﬂood events. For early ﬂood estimation, ﬂood mapping is a crucial component. This study used an integration of geospatial techniques to evaluate the ﬂood extent in District Dera Ghazi Khan, Pakistan. The modiﬁed normalized difference water index (MNDWI) was utilized to estimate the ﬂood extent using Landsat data. For a thorough ﬂood investigation, pre-ﬂood, during, and post-ﬂood images were obtained. The analysis enabled us to delineate ﬂood extent as well as ﬂood duration. The result showed that the ﬂood continued for nearly 5 weeks in the study area. This proposed geospatial technique provides a framework for the identiﬁcation of inundated areas, which allows emergency responses to be focused on newly ﬂooded areas. Hence, the current study offers a novel perspective on ﬂood mapping and signiﬁcantly contributes to ﬂood monitoring.


Introduction
Flood disasters are frequently occurring with more severity worldwide [1,2], and they account for more than 50% of all natural disasters that happen [3,4]. In the recent decade, flood disasters have caused over 0.25 million causalities and also registered a total of about USD 4.8 billion in economic losses [5,6]. Flooding is the primary global cause of catastrophic loss and destruction, and the number of people exposed to floods is increasing faster than the world's population [7,8].
South Asian countries, such as Bangladesh, India, and Pakistan, claim the highest number of flood-related causalities and economic losses [7,[9][10][11]. In the recent decade, flood disasters have been common in Pakistan and have caused severe damage to all socio-economic sectors [12][13][14].
Floods in Pakistan generally occur in the monsoon season because of an excessive amount of rainfall in the upstream watershed of the main rivers, the Indus and the Chenab [15][16][17]. As a result, floods with high magnitudes have occurred in floodplain areas, which ultimately resulted in a huge impact on people's lives, property, and the overall economy of the country [18,19]. Despite the fact that this catastrophic flood menace cannot be entirely prevented, the effects can be lessened by employing effective remote sensingbased flood risk mapping [17,20]. The remote sensing technique is the most affordable way to provide extensive data on all stages of a flood disaster, which may be used for flood mapping [9,13].
Therefore, the main goal of the current research is to delineate flooded areas using suitable satellite-derived MNDWI water indices. Results from the water indices propose a step-wise approach for accurately performing flood mapping. Such an approach gives 2 of 6 an insight into the flood situation in this flooded area with minimum inclusion of other non-water classes. In addition, in case of flood disaster occurrence, such insights can assist in proper decision making for emergency flood management in District Dera Ghazi Khan.

Study Area
Tehsil Taunsa, District DG Khan was chosen for this study. Its geographic range is between 30.4078 • N and 70.5265 • E (Figure 1). It encompasses the south-western portion of Punjab and is bordered by the Punjab districts of Muzaffargarh, Layyah, and Rajanpur. The entire study area covered an area of 11,294 km 2 [21]. The highest and lowest rainfalls ever recorded were 50 mm in July and 2 mm in October, respectively. The average annual rainfall is 22.18 mm [22]. The locals' primary economic activity is agriculture. In the study area, cotton, sugarcane, wheat, rice, and sunflowers are the primary crops grown. a step-wise approach for accurately performing flood mapping. Such an approach gives an insight into the flood situation in this flooded area with minimum inclusion of other non-water classes. In addition, in case of flood disaster occurrence, such insights can assist in proper decision making for emergency flood management in District Dera Ghazi Khan.

Study Area
Tehsil Taunsa, District DG Khan was chosen for this study. Its geographic range is between 30.4078° N and 70.5265° E (Figure 1). It encompasses the south-western portion of Punjab and is bordered by the Punjab districts of Muzaffargarh, Layyah, and Rajanpur. The entire study area covered an area of 11,294 km² [21]. The highest and lowest rainfalls ever recorded were 50 mm in July and 2 mm in October, respectively. The average annual rainfall is 22.18 mm [22]. The locals' primary economic activity is agriculture. In the study area, cotton, sugarcane, wheat, rice, and sunflowers are the primary crops grown.

Data and Methodology
The Landsat 8-9 OLI/TIRS data were obtained from the United States Geological Survey (USGS) [23][24][25]. We obtained 6 Landsat images from 21 July and 21 September 2022. Furthermore, these images were utilized to estimate MNDWI water indices for flood mapping using ArcGIS 10.8 [26,27]. This index, first presented by Xu [28], effectively calculates inundated areas. In essence, it employs the green and shortwave infrared bands to extract the flood extent, as illustrated in Equation (1):

Data and Methodology
The Landsat 8-9 OLI/TIRS data were obtained from the United States Geological Survey (USGS) [23][24][25]. We obtained 6 Landsat images from 21 July and 21 September 2022. Furthermore, these images were utilized to estimate MNDWI water indices for flood mapping using ArcGIS 10.8 [26,27]. This index, first presented by Xu [28], effectively calculates inundated areas. In essence, it employs the green and shortwave infrared bands to extract the flood extent, as illustrated in Equation (1): The MNDWI index reduces accumulated noise while highlighting the water surface reflectance using the green band; the resulting values range from 1 to +1. The generated Environ. Sci. Proc. 2023, 25, 78 3 of 6 images show negative values for built-up regions and positive values for water areas, respectively, depending on how much water area reflectance is high and how little built-up area reflectance is low in the shortwave infrared band.
Optical Landsat data, which provide up-to-date information with comparatively high temporal resolution have been widely used in the detection of flooded areas [7,17,23]. On the other hand, the presence of clouds lowers the availability of Landsat images during flooding. Certainly, synthetic aperture radar (SAR) and radar satellites can pierce clouds and obtain images in all weather situations and are commonly employed for flood mapping under these circumstances [29,30]. SAR data are normally very expensive and, in our case, the free SAR data were not available and the study area was entirely cloud-free during the 2022 flood.

Spatio-Temporal Flooded Area Mapping
The 2022 flood extent was identified to delineate the most inundated areas in Taunsa, Dera Ghazi Khan, Pakistan. The results show that the cumulative flood inundation from pre-flood to post-flood instances is as shown in Figures 2 and 3. The results also show that the flood inundation spans an area of approximately 526 km 2 in the 21 July image, whereas the highest flood inundation of 1462.24 km 2 can be seen in the 31 August image. Additionally, the peak flood scenario remained constant until 7 September as rain and hill torrent flow persisted from the Suleiman range and reduced flood inundation to an area of 1196 km 2 at a pace of 38 km 2 /day. On 14 September, the inundation receded to about 702 km 2 , with a decreasing rate of 72 km 2 /day. Lastly, the floodwaters continued to recede and the normal stage was seen as a reducing trend persisting throughout the month of September, to the point that on September 21st, the extent of the inundation covered only 491 km 2 and returned to its pre-flood stage. The trend of the flood extent during three flood instances is presented in Figures 2 and 3. The MNDWI index reduces accumulated noise while highlighting the water surface reflectance using the green band; the resulting values range from 1 to +1. The generated images show negative values for built-up regions and positive values for water areas, respectively, depending on how much water area reflectance is high and how little built-up area reflectance is low in the shortwave infrared band.
Optical Landsat data, which provide up-to-date information with comparatively high temporal resolution have been widely used in the detection of flooded areas [7,17,23]. On the other hand, the presence of clouds lowers the availability of Landsat images during flooding. Certainly, synthetic aperture radar (SAR) and radar satellites can pierce clouds and obtain images in all weather situations and are commonly employed for flood mapping under these circumstances [29,30]. SAR data are normally very expensive and, in our case, the free SAR data were not available and the study area was entirely cloud-free during the 2022 flood.

Spatio-Temporal Flooded Area Mapping
The 2022 flood extent was identified to delineate the most inundated areas in Taunsa, Dera Ghazi Khan, Pakistan. The results show that the cumulative flood inundation from pre-flood to post-flood instances is as shown in Figures 2 and 3. The results also show that the flood inundation spans an area of approximately 526 km² in the 21 July image, whereas the highest flood inundation of 1462.24 km² can be seen in the 31 August image. Additionally, the peak flood scenario remained constant until 7 September as rain and hill torrent flow persisted from the Suleiman range and reduced flood inundation to an area of 1196 km² at a pace of 38 km 2 /day. On 14 September, the inundation receded to about 702 km², with a decreasing rate of 72 km²/day. Lastly, the floodwaters continued to recede and the normal stage was seen as a reducing trend persisting throughout the month of September, to the point that on September 21st, the extent of the inundation covered only 491 km² and returned to its pre-flood stage. The trend of the flood extent during three flood instances is presented in Figures 2 and 3.  As can be seen in Figures 2 and 3, the water remained in the study area for about 5 weeks after the maximum flood peak was recorded on 31 August 2022. The floodwaters subsided in two stages; during the peak stage, the water decreased gradually until the 7th of September at a rate of about 38 km² per day. Up to 21 September, the floodwaters in the moderate stage significantly decreased at a rate of 72 km² each day. The assessment of the results, however, showed that the agricultural and built-up areas suffered a great deal due to the massive inundation and duration of approximately 5 weeks. Our results reveal that Landsat images, together with satellite-derived MNDWI indices permit detailed flooded area delineation with reliable accuracy. A recent review paper, similarly, analyzed the same indices used to depict areas under flood water and regarded MNDWI to be the most suitable in terms of its ability to differentiate between turbid water and mixed pixels [31]. Suitable satellite data are critical for flood mapping [6,9]. The temporal relationship between satellite characteristics with flood occurrence is a vital parameter in flood mapping. For example, a low-resolution moderate resolution imaging spectroradiometer (MODIS) satellite (~250 m), with its daily revisit time, has been utilized during many floods to acquire flood mapping, but with questionable accuracy [9]. However, due to the lack of during-flood SAR satellite data, such as Sentinel-1, we employed Landsat data with a resolution of 30 m for flood mapping [2,17,30]. Albertini et al. [31], in their recent review article, also concluded Landsat is the most commonly used satellite for floodwater spatial coverage detection. Despite the better spatial resolution, Landsat satellite data are constrained and cannot obtain geospatial data on time, which usually reduces their applicability to flooded area mapping [27,29]. Furthermore, the MNDWI index used achieved reliable accuracy in terms of flooded area extraction and has been efficiently used in other studies as well for detecting areas under flood water [32].

Conclusions
Tehsil Taunsa is particularly vulnerable to frequent riverine and hillside torrent floods. Locals living along nullahs and the Indus River deal with this flood threat almost every year. In the study region in 2022, floods produced by a two-week-long persistent Areas that have been flooded for several days can be detected using Landsat data [15,17,25]. Landsat data were utilized because they are also capable of accurately detecting flooded areas in built-up and agricultural areas [23,25,30]. However, to properly estimate the extent of flooding, field research and high-resolution data are required. However, due to their high cost, radar and SAR data are not used in the research area. To really observe the flooding scenario for this study, a field survey was also carried out.
As can be seen in Figures 2 and 3, the water remained in the study area for about 5 weeks after the maximum flood peak was recorded on 31 August 2022. The floodwaters subsided in two stages; during the peak stage, the water decreased gradually until the 7th of September at a rate of about 38 km 2 per day. Up to 21 September, the floodwaters in the moderate stage significantly decreased at a rate of 72 km 2 each day. The assessment of the results, however, showed that the agricultural and built-up areas suffered a great deal due to the massive inundation and duration of approximately 5 weeks.
Our results reveal that Landsat images, together with satellite-derived MNDWI indices permit detailed flooded area delineation with reliable accuracy. A recent review paper, similarly, analyzed the same indices used to depict areas under flood water and regarded MNDWI to be the most suitable in terms of its ability to differentiate between turbid water and mixed pixels [31]. Suitable satellite data are critical for flood mapping [6,9]. The temporal relationship between satellite characteristics with flood occurrence is a vital parameter in flood mapping. For example, a low-resolution moderate resolution imaging spectroradiometer (MODIS) satellite (~250 m), with its daily revisit time, has been utilized during many floods to acquire flood mapping, but with questionable accuracy [9]. However, due to the lack of during-flood SAR satellite data, such as Sentinel-1, we employed Landsat data with a resolution of 30 m for flood mapping [2,17,30]. Albertini et al. [31], in their recent review article, also concluded Landsat is the most commonly used satellite for floodwater spatial coverage detection. Despite the better spatial resolution, Landsat satellite data are constrained and cannot obtain geospatial data on time, which usually reduces their applicability to flooded area mapping [27,29]. Furthermore, the MNDWI index used achieved reliable accuracy in terms of flooded area extraction and has been efficiently used in other studies as well for detecting areas under flood water [32].

Conclusions
Tehsil Taunsa is particularly vulnerable to frequent riverine and hillside torrent floods. Locals living along nullahs and the Indus River deal with this flood threat almost every year. In the study region in 2022, floods produced by a two-week-long persistent wet spell occurred upstream of the Indus and in the foothills of the Suleiman range, which resulted in flooding and had a significant impact. Standing crops that were ready to be harvested, homes, animals, and all kinds of infrastructures were affected as a result.
These frequent floods are a severe problem that calls for effective preparation and impact mitigation measures, particularly through effective post-flood monitoring procedures, and particularly in Tehsil Taunsa, District D. G. Khan.
This study used remote sensing data integration and jointly employed appropriate methods to demonstrate the flood inundation severity in Taunsa. In order to estimate the inundated areas, Landsat images of flood instances were analyzed using the GIS-based MNDWI index.
The flood mapping also revealed that the flood water persisted in the study region for a month, which increased and exacerbated flood damage in the study areas. This study found that geospatial techniques can be used to carry out advanced flood mapping, which is important for flood management. As a result, the current study offers an alternative perspective on mapping flood inundation utilizing free satellite data and techniques. Flood risk mapping is the first step toward analyzing flood risk. This study can be used as a basis for further evaluating flood risk assessment and management in the study area. The methodology employed can be integrated with other radar datasets and the Analytic hierarchy process (AHP) techniques to further identify the severity of future flood risk while developing flood risk zonation.