Operational Monitoring and Damage Assessment of Riverine Flood-2014 in The Lower Chenab Plain, Punjab, Pakistan, Using Remote Sensing and GIS Techniques

: In flood-prone areas, the delineation of the spatial pattern of historical flood extents, damage assessment, and flood durations allow planners to anticipate potential threats from floods and to formulate strategies to mitigate or abate these events. The Chenab plain in the Punjab region of Pakistan is particularly prone to flooding but is understudied. It experienced its worst riverine flood in recorded history in September 2014. The present study applies Remote Sensing (RS) and Geographical Information System (GIS) techniques to estimate the riverine flood extent and duration and assess the resulting damage using Landsat-8 data. The Landsat-8 images were acquired for the pre-flooding, co-flooding, and post-flooding periods for the comprehensive analysis and delineation of flood extent, damage assessment, and duration. We used supervised classification to determine land use/cover changes, and the satellite-derived modified normalized difference water index (MNDWI) to detect flooded areas and duration. The analysis permitted us to calculate flood inundation, damages to built-up areas, and agriculture, as well as the flood duration and recession. The results also reveal that the floodwaters remained in the study area for almost two months, which further affected cultivation and increased the financial cost. Our study provides an empirical basis for flood response assessment and rehabilitation efforts in future events. Thus, the integrated RS and GIS techniques with supporting datasets make substantial contributions to flood monitoring and damage assessment in Pakistan. and duration, and evaluate flood damages on different land uses, such as agricultural land and built-up areas. We used the identified relationship between land uses and flood instances to estimate the overall damage, flood extent, and propose emergency flood management for future floods on the lower Chenab plain.


Introduction
Flood disasters are among the most frequent and destructive of all-natural disasters, posing a potential threat to life and property. Every year, human lives, agricultural activities, and infrastructures are seriously affected by shattering flood disasters around the globe [1][2][3]. In the past three decades, flood disasters affected nearly 2.8 billion people and resulted in over 200,000 causalities with substantial damages to property and economy [4]. Floods account for ~47% of all weather-related disasters that occur across the world [5]. Climate change and rapid population increases in floodplains enforce and boost the frequency and magnitude of riverine flood damages [6,7]. As flood disasters increase in magnitude, the Asian region continues to face a large number of flood hazards and associated losses in lives and all kinds of infrastructure and economic progress [4,7]. In recent years, the south Asia region has been experiencing intense flooding with increased frequency, especially in Pakistan [8][9][10].
In Pakistan, flood events are recognized as a major natural hazard that historically originates from the Indus and Chenab rivers [11][12][13]. Since the creation of Pakistan, 1955Pakistan, , 1959Pakistan, , 1973Pakistan, , 1976Pakistan, , 1988Pakistan, , 1992Pakistan, , 1995Pakistan, , 1996Pakistan, , 1997Pakistan, , 2006Pakistan, , 2010Pakistan, , and 2014 are recorded as the years for destructive floods that resulted in adverse impacts on human lives, property, and the country's economy [9,14,15]. Floods in Pakistan mostly occur in the monsoon months of July and September, due to heavy rains and the melting of snow upstream in the Himalaya region. This results in riverine floods that produce tremendous detrimental impacts on human lives, agriculture, and infrastructure [16][17][18]. For instance, the 2010 riverine flood caused a large inundation that covered an area of 70,238 km², with 884,715 affected houses [8,12,18]. Likewise, the 2011 flood affected 5.88 million people and damaged standing crops and infrastructure covering an area of 16,440 km² and 1882 km², respectively [19]. In 2014, heavy rains coupled with the melting of glaciers in the upstream part of Chenab caused flash flooding in mountainous areas and riverine floods in the upper and lower Chenab plain [20][21][22]. It was considered one of the worst riverine flood catastrophes in terms of damages to standing crops, housing, and infrastructure along the floodplain areas of the Chenab river [13,16,22,23]. This indicates that flooding is a serious problem in Pakistan that requires significant efforts to reduce its effects, mainly through effective post-disaster monitoring and management, especially in the Chenab plain [22][23][24]. As such, we provide a scientific basis that provides a rapid flood 2014 mapping and monitoring using Landsat-8 data, with a focus on the lower Chenab plain. We first provide the limitation of previous studies on the riverine flood of 2014 that focused less on post-flood mapping, monitoring, and damage assessment. These previous studies mainly focused on the upper Chenab plain [13,16,21,22] without focusing on the lower Chenab plain. This is despite the fact that the lower Chenab plain is a fertile flood-prone plain and considered an economically underdeveloped region in the Punjab province of Pakistan. Despite being an understudied area, negligence adds to the problems with a lack of localized flood management strategies. Therefore, accurate post-flood mapping, monitoring, and evaluating damages are the utmost requirements for rapid flood risk assessment in lower Chenab plain. Secondly, to a basis for rapid flood risk mapping, we make use of multispectral remote sensing open data in conjunction with Geographical Information System (GIS) techniques for monitoring and evaluating damage assessment in the severely flood-affected areas of the lower Chenab plain, Punjab, Pakistan.
The use of remote sensing and GIS techniques is chosen as it contributes to an exploration of flood causes and additionally provides accurate mapping of flood extents that enables detailed investigations of flood instances [18,19,25]. Moreover, these techniques can be applied in damage assessment to standing crops and infrastructure. Effective flood modeling and mapping are important for flood assessment, loss estimation, and sustainable land use planning along flood plains to mitigate flood risk effectively [25,26]. Multi-temporal remote sensing images provide a wide source of low-cost information with a reliable accuracy that can be beneficially utilized for flood mapping [27,28]. Furthermore, multispectral remote sensing-derived indexes and GIS-based classifications can be exploited to detect flooded areas [18,29,30]. A wide range of open-source remote sensing data has permitted valuable and accurate historical data, which is essential for a comprehensive study on flood disaster mapping and damage assessment [31]. For example, moderate resolution Landsat data have provided free, up-to-date satellite images across the globe since 1972, which can be used for flood disaster monitoring and damage assessment [32,33]. High-resolution Google Earth (GE) images provide historical and recent data that can also be used for flood monitoring and temporal flood mapping [34]. GE images can also be used as input datasets to digitize land uses, particularly useful in small areas [5], and for the validation of supervised classified images [35]. Thus, the current study uses remote sensing and GIS techniques to detect flood inundation and duration, and evaluate flood damages on different land uses, such as agricultural land and built-up areas. We used the identified relationship between land uses and flood instances to estimate the overall damage, flood extent, and propose emergency flood management for future floods on the lower Chenab plain.

Study Area
The study area is the lower Chenab plain, which is located in the downstream part of river Chenab, central Pakistan. It is located between 70°41'13.2"E and 29°6'0"N, and 71°37'58.8"E and 30°31'33.6"N, within the three main districts of Punjab province, namely the Multan, Muzaffar Garh, and Bahawalpur districts ( Figure 1). It is one of the fertile plains of Punjab province, where floods occur almost every year. The Chenab river originates from the Indian state of Himachal Pradesh. Then, it flows through Indian-occupied Kashmir and enters at Marala Headwork into the province of Punjab, Pakistan. The southern part of the study area includes the Panjnad Barrage where the Sutlej River joins the Chenab, as shown in Figure 1. The total length of the River Chenab is about 974 km, of which 729 km fall within Pakistan. Of the 729 km, we selected approximately 215 km as our study area. We then divided the 215 km into three equal sections, each approximately 71 km long (Figure 1), to enable the visualization and detailed investigation of the results. We also created a 6-km-wide buffer around the Chenab River that covers a distance of about 215 km from north to south as the final extent of our study area. The mean annual Chenab discharge is 1.52 million m³ [36,37]. The highest temperature documented in the study area is 50°C in the month of June with the lowest recorded temperature of 2°C in the month of January [20]. The average yearly rainfall is 157 mm, with July receiving a maximum of 45 mm and October receiving a minimum of only 2mm [20]. The main source of livelihood in the study area is agriculture with main crops being sugarcane, cotton, maize, rice, and fodder. Figure 1. Location of the lower Chenab plain in Pakistan, and the zoomed area represents the river Chenab system, and also shown is the 6-km buffer zone and the extent of our study area.

Materials
The main data source of this study is Landsat 8 imagery (Table 1), acquired from the United States Geological Survey (USGS). Landsat images have a temporal resolution of ~16 days that has provided free up-to-date images across the globe since 1972 [22,33,38,39]. However, the study area provides a unique opportunity with a high temporal resolution of ~8 days since the area lies within two adjacent Landsat paths of 150 and 151 ( Figure 2; Table 1), which enables us to make an in-depth analysis of the flood duration and inundation. In total, we used 9 temporal images acquired between 7 August and 27 November 2014 (Table 1). We used 8 temporal images for flood monitoring. Furthermore, we used 3 temporal images for flood damage assessment using land use/cover classification during flood instances in lower Chenab plain (Table 1).
We also used GPS and the Google Earth (GE) platform to collect spatial training datasets, which were used to support land use/cover classification. The GPS data were collected after the flood waters had ceased, as it was challenging to do field survey during flooding time, due to flood inundation. These points were also used to validate classification and inundation results. Due to its high resolution, the GE platform provides a better visualization of the real-world scenario of land use/cover change in the study area [34]. Furthermore, GE images of three flood instances (and with similar acquisition dates as the Landsat-8 images) were taken for cross-comparison with the Landsat-based land use/cover classification. The Pre-flood image was acquired on 4 August, 2014, the Co-flood image was taken on 13 September, 2014, and the closest available Post-flood image was acquired on 11 November, 2014. To further understand post-flood situation in the study area, a comprehensive field survey was conducted to identify damages and affected areas and observe the Post-flood effects and rehabilitation process.

Methods
The overview of the used methodology is presented in Figure 3. Firstly, we applied radiometric correction, layer stacking and mosaicking and resampling to the Landsat-8 images before use. Then, we converted the digital number values on the Landsat images to reflectance using ENVI's Radiometric Calibration Tool [18]. Finally, we used Arc GISs' spatial analysis tools for Layer stacking and mosaicking, and resampling. The pre-processed Landsat-8 pre-flood, co-flood, and post-flood images were subjected to supervised classification for land use/cover mapping in ArcGIS 10.5, as shown in Figure 3. Landsat data are most common and is often used for different land use/cover mapping and water extraction mapping [13,40]. The Supervised classification method is used to extract information from Landsat data [29,41]. In addition, a maximum likelihood (ML) approach was applied for land use/land cover classification of flood instances. The ML approach is widely used and is easy to apply for land use/cover classification. The ML approach has also been applied to land use/cover mapping of pre-flooding, co-flooding, and post-flooding [11,18,22,[39][40][41][42]. We identified water, agriculture land, and vegetation, built up, sand, barren land, and deposited material as the six land-use classes in the study region. To facilitate the ML approach for classification, a total of 360 spatial training samples were collected from a GPS field survey and comparatively high-resolution GE images, and were further assigned to the pixel values of the most probable land use/cover class. All spatial training datasets were compiled and prepared in Microsoft Excel and were imported into the Arc GIS 10.5 environment. The classified images were converted into a shapefile format and used as input for the GIS-based spatial overlay analysis. This allows us to spatially compare and intersect all classified land uses to provide a clear picture of flood inundated and damaged land uses in lower Chenab.
For damage assessment, we defined a damaged area as an area that experienced a change from 'agricultural' and 'built up' in pre-flood to 'deposited material' in post-flood. A pre-flood image was used as a reference and a change detection technique was used to evaluate damages to agriculture and built-up areas in the study area.

MNDWI Index
The Modified Normalized Difference in Water Index (MNDWI) was used to delineate the spatial pattern of flood-2014 inundation along the lower Chenab plain. The MNDWI index is widely used for the rapid delineation of floodwater required for flood monitoring and assessment, and has been compared to other indices, i.e., the NDWI and Water Ratio Index (WRI) [43,44]. The Normalized Difference Water Index (NDWI) presented by McFeeters [45], is very efficient at delineating water information but has difficulties in case a built-up area exists in water environment. Therefore, Xu [46] presented an effective MNDWI index, which is much better at distinguishing between water and built-up areas. The MNDWI is calculated using green (Band 3) and shortwave infrared (Band 6) wavelengths to delineate water, as given by Equation (1): where GREEN (Band 3) = Green wavelength (0.53-0.59mm) and SWIR (Band 6) = Short-wave infrared Wavelength (1.57-1.65 mm) This index is utilized for the removal of built-up area noise, and it uses Band 3 wavelength to maximize water reflectance. The resulting value ranges from −1 to +1. The low water reflectance and high reflectance of built-up in Band 6 result in positive values of water and negative values of built up in the MNDWI image. The limitation of the MNDWI index is that it does not efficiently distinguish between hill shadow and water body [47,48]. From literature, the MNDWI index has been used for the extraction of water on a flat plain with scattered built-up areas using Landsat 8 OLI and SAR data, achieving high accuracy and better performance. Thus, it was also strongly recommended, in comparison with other indexes, to be applied on Landsat images [49][50][51][52].

Classification and Inundation Validation
The classified images and index-derived inundation were validated using GPS and Google earth (GE) points. For this, the random sample points' tool within the Arc GIS spatial analyst was used to extract random points on the classified land uses of inundated and non-inundated images. These points were then converted into a kml format and overlaid on GE. The values of these points were evaluated using visual interpretation and expert knowledge. The accuracy of classified images was assessed using 557 Arc GIS random points and 150 GPS points. We have used 30 Arc GIS random points of class, namely: water, agriculture, built-up, barren, sand, and deposited material, to each flood instance image. In addition to this, we have further used 25 GPS points of above-mentioned classes to only post-flood instance for the validation process. The accuracy of inundation maps was assessed using 400 ArcGIS random points and 200 GPS points. Furthermore, for each inundation map, we have used 50 ArcGIS random points and 25 GPS points for the validation of water and non-water classes. Lastly, the results from the comparison of random samples and GPS with GE images were used to create the confusion matrix. The confusion matrix was used to evaluate the accuracy of the ML land use/cover classification and MNDWI inundation. We also used the Kappa Coefficient (KC) as an indicator to validate the qualitative agreement, either positive or negative, between classified samples and ground-truth points. It is normally calculated from a statistical assessment to evaluate the proportional improvement by the classifier over ground-truth samples to land-use classes.

Accuracy Assessment
The overall inundation accuracy obtained from all images is above 90%, as shown in Table 2. The highest overall accuracy of 92% is obtained from the peak-flood image (17 September) and the post-flood image (27 November), while the least overall accuracy of 88% is acquired from the postflood (11 November) image. Similarly, the obtained average user and producer accuracy of both classes is nearly 90%. The highest user accuracy of the water class is attained from the peak-flood image (17 September), which is 96%, and the least is 85%, which is obtained from the post-flood (11 November) image. On the other hand, the highest user accuracy of the non-water class is achieved from the post-flood (27 November) image, which is 93%, and the least is 86%, which is from the peak-flood (17 September) image. Furthermore, the highest attained producer accuracy of the water class is 94% from the post-flood (26 October) image, and the least is 87% from the peak-flood (17 September) image. The highest obtained non-water class producer accuracy is 96% from the peak-flood (17 Sepember) image and the least is 86%, which is obtained from the post-flood (11 November) image. The highest KC accuracy of 85% is achieved from the post-flood (27 November) image and the least KC accuracy is obtained from the post-flood (11 November) image. The average KC accuracy of all the images is above 80% ( Table 2).
The highest overall supervised classification accuracy obtained is 92% from the post-flood image. Whereas, the lowest accuracy of 85% was obtained from the co-flood image, as shown in Table  3. The KC accuracy shows that the highest accuracy of 90% is obtained from the post-flood (27 November) image and that the lowest of 81% is obtained from the co-flood (17 September) image. Similarly, the user and producer accuracies were also calculated, which are listed in Table 3. The highest user accuracy was attained by water, barren land, and sand classes of 92%, 96%, and 99%, respectively. Likewise, the highest producer accuracy was obtained by the water, built-up area, and deposited material classes at 99%, 98%, and 95%, respectively. The overall land-use classification accuracy suggests that classified images are reliable for further analysis.

Flood Mapping and Monitoring
The temporal flood 2014 extent maps were prepared and used to determine the most inundated areas in the study area, as a tool for flood monitoring. Figure 4a,b shows the flood inundation with a flood peak on 17 September, which remained stable until the 24 September, and then the inundation gradually decreased till post-flood 27 November. Figure 4a and Table 4 show that northwestern Muzaffargarh and northeastern Multan Saddar part was the most inundated and affected region along with the southeastern ShuJabad and Jalalpur Pirwala. The flooded area receded in three phases: In the first phase from 24 September to 10 October, the southern part gradually receded from west to east. The central part began to recede from 10-26 October 2014 in the second phase, and in the third phase, almost all the flood water receded until 27 November 2014. Hence, the flood duration result shows that Muzaffargarh and Multan Saddar remained inundated for almost two months and reported the most affected areas by flood-2014 in lower Chenab plain.
The peak-flood (17 September) extent is compared in Figures 5a and b. In Figure 5a, the MNDWI index shows the accumulated flood extent that covers an area of 1380km². Whereas, in Figure 5b, the classified image shows that about 1330.03 km² of the area was inundated, which represents a 2% deviation from the MNDWI results. The result shows that both the MNDWI index and supervised classification produce almost similar inundation areas in lower Chenab plain. Furthermore, these inundation results, when incorporated with the GE images acquired at the same time as the Landsat-8 images, also showed the severe spatial pattern of inundation, and they also validate our estimated flood inundation extent, confirming that flood waters remained for almost two months and receded very slowly (Section 3.3).      (Figure 6). The water body showed only 8% in the pre-flood situation. However, after the flood occurrence in the month of September, 2014, an abrupt change appeared in the water body that represents a massive increase to 50% of the area, as shown in Figure 7. The massive water had not only inundated but also severely affected all other land uses. The built-up area considerably decreased from 17.39% to 11.85% in the study area. Likewise, we noticed a massive increase in water that affected the agriculture/vegetation covered area, which decreased from the initial total area of 2536.11 km², from 52.37% to 32.40%. The noticeable decrease of sand was from 8.3% to 1.1%. The slight decrease is noticed in barren land from 4.8% to 4%. Thus, the massive increase in water caused a large part of the agricultural/vegetation area to be flooded and destroyed most of the standing crops. In addition, the built-up area was also severely affected and flooded.

Land Use and Land Cover Changes
In the post-flood period (Figure 8), water receded to its original pre-flood stage and decreased to only 6% from 50% in the co-flood image. Similarly, built up also regained its original pre-flood stage and increased to 17.79%. After the flood occurred, water receded back but large amounts of sediment and other materials remained, representing a dramatic increase from 5.92% to 30%. Similarly, an increase in the vegetation and agriculture area can be noticed from 32.40% to 37.50%. The barren land area did not change much as it remained constant in each instance. Sand slightly increased to 4.4% from 1.17% of total 2536.11 km² area, as shown in Figures 8 and 9. The total area in km² of each class is statistically represented in Figure 9. Hence, all classified images show a similar pattern of land use/cover classes compared to high resolution GE images ( Figure 10).
The classification result and GE images show the same pattern of land use changes in the study area. The agriculture and built-up areas along the river are clearly visible in the pre-flood image ( Figure 10A). Conversely, the co-flood image ( Figure 10B) depicts the abnormal change in water area that has increased drastically and caused an inundation to all the surrounding agricultural and builtup areas. As shown in Figure 10C the post-flood image shows the huge amount of deposited material along the river and its surrounding areas that dramatically increased and destroyed standing crops, which caused late sowing of agricultural crops, and further leads to low productivity in the study area.

Damage Assessment
The results revealed that agriculture/vegetation areas were the most severely inundated, ~495 km², and while inundation in the built-up area covered 229 km², from the total inundated area of 1330 km², as shown in Figures 11a,b. The damage result further showed that a total of 361.43km² of agricultural land and 187.36km² of built-up areas were damaged due to deposited material. This result shows that the Multan district is the most affected/damaged district from flood-2014, as also shown in Figure 12a. Further details on the damaged agriculture land and built-up areas of the involved districts are listed in Figure 12b.

Discussion
Pakistan is a flood-prone country with historical records of various magnitudes of flood events [37,53,54]. In the past decade, flood disasters in Pakistan have surpassed all other disasters in terms of the frequency of occurrences, and also killed over 5700 people coupled with severe damages to the country's economy [5,12,15,17]. In the year of 2014, the flood started from late monsoon rainy season when the river Chenab inflow significantly increased from its upstream tributaries and resulted in high discharge, which exceeded the limits of the river flowing capacity and subsequently caused a huge inundation in the study area, as shown in Figures 7and10b.
Our results show that open satellite data coupled with an ML-supervised classification approach and the MNDWI allow the delineation of flood inundated and damaged areas with high overall accuracy. To obtain high accuracy levels, appropriate satellite data is critical. Firstly, temporal relationship between flood occurrence and satellite characteristics is an important parameter in flood modeling. For instance, a low-resolution (~250 m), multispectral Moderate Resolution Imaging Spectroradiometer (MODIS) satellite with its daily revisit time has already been used to obtain a coflood image for certain flood events in Pakistan, but with questionable accuracy [19,55]. However, due to absence of any co-flood open access high resolution, e.g., Sentinel-1, for the study area, we opted to use moderate spatial resolution Landsat data (~30 m) for flood monitoring and assessment [9,18,22]. Despite a relatively better resolution, Landsat data are limited in such a way that they lack the timely acquisition of geospatial data, which sometimes reduces their suitability for flood monitoring and inundation mapping [26,39,56]. In the study area however, two adjacent Landsat satellite paths (150 and 151) made it possible to acquire high temporal resolution Landsat data (~7 days), against a single temporal resolution of 15 days for an area [9,18,57]. As such, the high temporal resolution of Landsat data enables a detailed investigation of operational flood mapping and monitoring in a study area. Secondly, existences of cloud cover can restrain the availability of flood instances images. Certainly, SAR (Sentinel 1,2,3) and RADAR satellites can easily penetrate clouds and acquire images in comparison to optical satellites (Landsat, MODIS) [26,58]. In our cases, the free SAR data were not available and the study area was entirely cloud free during the 2014 flood, since the causative meteorological events happened in upper Chenab plain. These events generated the flood peak that happened almost 10 days prior to the flooding in lower Chenab plain. Thirdly, consideration of land use/cover change of the flooded area is also an important factor in choosing the satellite data. The inundated agriculture and built-up area can easily be detected by most Radar SAT and optical satellites [59]. With multispectral Landsat images, we can recognize flooded areas even after several days of flood occurrence [41,42,56]. Nevertheless, we used Landsat-8 data as it is also possible to detect sediment material over agriculture fields and built-up areas with reliable accuracy [26,59]. Despite that, high resolution SAR data and a field survey are necessary for accurate and reliable flood mapping and damage assessment [59,60]. Finally, we have to consider the accuracy of the applied methods to the satellite data in order to ascertain the reliability of the flood monitoring results. This study produces an overall accuracy for the MNDWI index of almost 90% (Table 2), while classified images had an average overall accuracy of about 88% ( Table 3). The water class attained the highest overall accuracy in this study, which has the ability to obtain a highest possible accuracy of 100% [49]. Despite that sand also shows reliable accuracy, in some occasions, wet sand was confused with water and vegetation/agriculture land. This is also noted in some instances where agriculture and built-up areas are also confused in transition areas due to mixed pixels.
During the flood-2014 event, the peak-flood water arrived on 17 September, 2014, and then remained constant for a while. The water slightly receded at a rate of 7km²/day, until 24 September, while a recharge of flood water continues from upstream areas. Here, the flood water started to decrease at a rate of 23km²/day until 26 October, exposing a once flooded area of 636.35 km². On 27 November, the river regained its pre-flood stage and flowed normally. However, inundation remained for almost two months in the study area, as shown in Figure 4. As shown in Figures 11 and  12, the flood-2014 caused huge inundation and damages to agriculture land and built up in the study area.
The field survey reveals that most of the local people live in mud houses, which are highly vulnerable to flood disaster and have no flood resilience capacity [25]. As a result, already vulnerable mud houses were badly destroyed by the flood, besides infrastructure being destroyed by the gushing flood-2014. Moreover, the survey reveals that the month of November is considered agriculture land preparation and sowing season for new Rabi season crops and simultaneously harvesting season of certain kharif crops, such as sugarcane and rice in the study area. As a result, vegetation/agriculture land slightly increased to 37% in the post-flooding instance. Likewise, the built-up area also increased to 17.7%, which regained its pre-flood situation, probably due to the receded flood water and reconstruction. Furthermore, the survey shows that the Chenab river carried a huge amount of sediment material from its upstream mountainous areas, and it was ultimately deposited along the lower Chenab plain in 2014. However, deposited material immensely increased to 30% in the post-flooding instance ( Figure 10). In our case, the field survey provided real-time reliable information, which is immensely required for comprehensive flood disaster assessment.
The moderate 30m resolution Landsat data that we used indicated that tracking inundated and damaged agriculture areas and built-up areas is reliable and acceptable. This mostly applies to damage assessment, which could be used for post-flood rehabilitation and a relief operation. However, for mapping the of detailed built up and agriculture area with a high overall accuracy, we must utilize high-resolution satellite data (RADARSAT, SAR), aerial photography, and extensive field survey, which are very often time and cost intensive and the required high-resolution satellite images are not open access. In this work, we did a field survey in order to evaluate real time damages, geo-location information, and the validation of Remote Sensing (RS) results. Thus, we correlated our survey findings for the validation of RS results and formulated an emergency damage assessment. Our damage analysis aims to provide a rapid and low-cost assessment of damaged agricultural and built-up areas but does not provide direct information on monetary losses and indirect damages, which may occur after flooding.
Landsat-8 data are affected by clouds, spectral sensitivity, and moderate spatial resolution and can only detect surface reflectance changes. Conversely, Radar satellite data can easily penetrate clouds and will perform properly in all weather conditions and detect changes in vegetation structure and moisture [59,60]. Radar data is more reliable when extracting flood inundation and damages in agricultural and built-up areas [25]. This well-known advantage allows Radar data to find the inundated and damaged areas in the post-flood stage in order to carry out an emergency flood damage assessment [58][59][60]. However, Radar data is costly and is not applicable for historical data analysis. However, despite such limitations, the integration of field and ancillary data allowed us to extract flood inundation and damage assessment with reliable accuracy. We noticed a number of factors that affected the tracking of damage assessment in this study. The lower Chenab floodplain contains numerous mangoes and other Bela forests. The co-flood Landsat image reveals that flood inundation beneath thin mangoes and other tree canopy covers was easily detected but thick or dense canopy covers' pixels were not identified due to the moderate resolution. This resulted in addition to the overall omission error and also caused the least identification of inundated areas in change detection analysis during flood instances. Secondly, a misclassification error was found in the classified images as the open land classified as the deposited material class in the post-flood image due to the moderate resolution and spectral sensitivity. Then, this error was removed and corrected through field-survey findings. Finally, the moderate 30-m spatial resolution of the Landsat images contributed to mix-up and misclassification in built-up areas where floodwaters pass through housing structures, roads, and other features. This increases the brightness of otherwise dark pixels (water), resulting in an increase in omission errors, which is a limitation for built-up damage assessments.

Conclusion
In this study, we presented a low cost and user-friendly flood monitoring and damage assessment with the integration of open access optical remote sensing and appropriate processing methods jointly exploited with field data. In particular, we used Landsat-8 data and processed them with open-source, GIS-based, supervised classification for damage assessment and a satellite-derived MNDWI index for flood inundation and monitoring. The classified images and inundation mapping produced an excellent overall accuracy of about 90%, which is validated and shows reliable results. we conclude that these methods have been proven to be useful for estimating and understanding a future flood phenomenon with its diverse impacts. These methods do not require more time and provide near real-time information using the user input with indigenous knowledge and expertise. Furthermore, the classified result reveal that the agriculture sector has been the most affected land use/cover in the study area. In a large context, people are mostly engaged with agricultural activities, and this pattern of flood impacts are a major concern as it directly influences the livelihood of the floodplain community. The index result shows that about 75% of the area experienced severe flooding, which lies mostly in the southern and central part of study area, and further revealed that flood inundation remained for almost two months. Finally, despite the fact that flood disasters are recurrent phenomena, our study proves that the used datasets and methods can be useful for emergency, real-time, automated flood monitoring and damage assessment in order to formulate emergency flood disaster management, particularly for relief and response operations. In contrast, traditional flood monitoring and damage assessment with paid on-demand data provides comparatively accurate results, but it is a more resource-and time-consuming process.
Despite that our study is local to Pakistan, a few points from our results can be applied in other areas and thereby contribute to science and global perspectives on flooding. In this paper, we have shown that the existing supervised classification methods, combined with the MNDWI and groundbased validation point, can solve a problem of flooding and provide flood management solutions to local authorities. The combination of supervised classification and the MNDWI has contributed to the accurate mapping of flood extents using Landsat imagery. Despite that these two methods are mostly applied separately, the combined application in flood mapping can help to produce reliable results, with complimentary properties. Furthermore, in the unique agricultural environment along the river, the approach of combining two classification methods can provide insights into flood instances and inundation in both build up and agricultural areas, enabling a rapid check on damages to both built up and agricultural crops and thereby providing an estimation of loss that can support ground-based research. We have also shown that despite a low temporal resolution of Landsat images, the exploitation of two adjacent Landsat paths can provide a high temporal resolution, which can be applied not only to flood management and monitoring but also other fields that require high temporal resolution data, especially where another dataset is not available. This exploitation increased the temporal resolution by two-fold, from ~15 days to ~ 7 days, which can also be factored into operational flood monitoring for emergency responses, such as early response and relief operations. Thus, our study has given another perspective to flood monitoring and management using the available and free satellite imagery.