Rapid flood inundation mapping using social media, remote sensing and topographic data

Flood events cause substantial damage to urban and rural areas. Monitoring water extent during large-scale flooding is crucial in order to identify the area affected and to evaluate damage. During such events, spatial assessments of floodwater may be derived from satellite or airborne sensing platforms. Meanwhile, an increasing availability of smartphones is leading to documentation of flood events directly by individuals, with information shared in real-time using social media. Topographic data, which can be used to determine where floodwater can accumulate, are now often available from national mapping or governmental repositories. In this work, we present and evaluate a method for rapidly estimating flood inundation extent based on a model that fuses remote sensing, social media and topographic data sources. Using geotagged photographs sourced from social media, optical remote sensing and high-resolution terrain mapping, we develop a Bayesian statistical model to estimate the probability of flood inundation through weights-of-evidence analysis. Our experiments were conducted using data collected during the 2014 UK flood event and focus on the Oxford city and surrounding areas. Using the proposed technique, predictions of inundation were evaluated against ground-truth flood extent. The results report on the quantitative accuracy of the multisource mapping process, which obtained area under receiver operating curve values of 0.95 and 0.93 for model fitting and testing, respectively.


Introduction
Flooding causes considerable damage to people, infrastructure and economies of many countries of the world.During a flooding event, rapid estimation of inundated areas is critical to effectively manage response operations.Emergency managers require timely and accurate information on areas affected by floodwater to prioritise relief efforts and plan mitigation measures against damage.
Nowadays, large amounts of geospatial data can be derived from a multitude of sources.In a time-critical disaster situation, utilisation of multiple data sources is particularly desirable.Remote sensing from satellite and aircraft platforms can provide detailed snapshots of a situation.However, these frequently do not provide the spatio-temporal coverage required for flood assessment and may not be validated in the field.For example, Landsat 8 provides multispectral data at 30-m spatial resolution, yet urban and semi-urban environments can contain many features affecting floodwater movement and presence within this cell size.Remote assessments can also have quality and uncertainty issues resulting from the choice of processing algorithm, mixed spectral responses and background noise (Jiang et al. 2014).
Furthermore, the increasing proliferation of consumer devices, such as smartphones, tablets and wearable logging sensors is facilitating the creation of large streams of data (Goodchild 2007).However, despite much interest in using such sources in natural hazard assessment (Goodchild and Glennon 2010), exploiting this information is not trivial.The data usually have no validation or assessment of quality (Goodchild and Li 2012) and may contain deliberate or unintended bias (Xiao et al. 2015).
To address these issues, we describe a method to integrate information derived from multiple sources to aid in the estimation of flood inundation extent.Observations made by in situ crowd users are combined with medium-scale optical remote sensing and a highresolution terrain model.Our main contribution is a probabilistic model, which quantitatively evaluates the contribution of each data source, and exploits contemporaneous remote sensing, social media and high-resolution topographic models to estimate flood inundation extent.Although our sample of social media is relatively small (n = 205), the importance of this information can be quantified by our method and use this ranking within an overall flood presence probability prediction.This enables even relatively small amounts of crowdsourced information to add value to hazard assessments.Validation and assessment of the proposed fusion-conflation model is undertaken demonstrating quantitative improvements over mapping flood extent using a single data source.

Crowdsourcing and disaster assessment
In this section, we review approaches and research undertaken to exploit crowdsourced data during natural disasters.We separate this work into two separate areas.Section 2.1 deals with citizen science type approaches where the ''crowd'' is actively engaged or directed.Section 2.2 reviews studies using web-harvested geospatial data.

Volunteer-based assessment
The use of crowdsourced data has been explored in various papers and research programmes in relation to natural disaster assessment.Coordinated deployment of crowd contributors is relatively well explored following earthquakes (Zook et al. 2010;Barrington et al. 2011;Ghosh et al. 2011;Kerle and Hoffman 2013).These approaches typically rely on citizen-based satellite or aerial image interpretation to assess disaster damage.Crowdsourcing software platforms may also be used to assess ground-collected images captured using smartphones.For example, the UNITAR-UNOSAT project GeoTag-X aims at using volunteers to categorise images captured during disaster situations and providing assessments, such as floodwater extent, damage to buildings or the safety of temporary shelters for refugees (GeoTag-X 2015).Poser and Dransch (2010) evaluated the use of volunteered geographic information for flood damage assessment.The authors surveyed individuals affected by a river flooding event to understand their perception of inundation and compared this to hydraulic modelling predictions, identifying that non-professionals can make useful assessments of water extent and depth.Hung et al. (2016) describe the use of geo-referenced citizen reports in relation to flood events.They describe a method to assess the credibility of the reports based on a probabilistic model trained on a previous flood event.The model was then used to categorise credibility on a subsequent event, highlighting that volunteered information may be automatically ranked to enable more rapid and adoption within emergency management contexts.
Le Coz et al. (2016) highlight the use of volunteer-based hydrological monitoring and prediction.They describe the adoption of citizen-captured video footage to estimate hydraulic data (e.g.water velocity and discharge) in flash flooding situations.Furthermore, the authors describe development and deployment of an interpolation approach using citizen's photographs of maximum flood levels in conjunction with LiDAR surface models, and supplementary field assessment by professionals, in order to generate water depth maps.
In our work, we choose to investigate use of social media data as that follows events in real time, enabling capture of unpredictable situations and with the potential of a large sample of well-positioned data.For these reasons, we focus the remaining review of the literature on this topic.

Social media
Using information from social media for aiding decision-making is a problematic task.This is due to multiple reasons including the challenge of extracting relevant information (e.g.identifying topicality) within unstructured or semi-structured web-harvested data, unknown quality (there may be little or no relevant metadata), and difficulty integrating it with other sources (which may in turn have their own issues of quality and uncertainty).When using social media, biases may also be present, for example, from a lack of digital engagement within certain demographics of the populations of particular areas (Xiao et al. 2015).Panteras et al. (2014) cross-reference geotagged points of images mined from Flickr and tweets to estimate spatial footprints events.The authors extract toponyms expressed in tweets to aid in estimating the viewing direction of Flickr images depicting a wildfire event, leading to a more accurate delineation of the event.However, the authors do not utilise remotely sensed earth observation data within their methodology, which is often available in disaster situations.
Several research projects have analysed social media specifically during flood disasters.Albuquerque et al. (2015) provide an example of the use of social media during flood events.They built a statistical model of tweets as they relate to authoritative river levels.In particular, the model tested the association between tweet locations and the water level in flood-affected catchments.Schnebele et al. (2014b) describe an interpolation-based approach that incorporates nonauthoritative data such as geotagged airborne (oblique) photographs and terrestrial videos to provide a spatio-temporal flood damage assessment.Furthermore, integrating social media data with remote sensing, river gauge measurements and digital elevation models (DEM) has been shown to be effective at producing near-real-time flood hazard maps (Schnebele and Cervone 2013).Fusion-based approaches can also aid extent estimation when the data sources have incomplete coverage (Schnebele et al. 2014a).However, these approaches rely on subjective user determination of weights for the model variables.

Study area
From December 2013 to February 2014, many parts of the UK were affected by flooding.These events resulted from a series of major storms generating intense rainfall, severe coastal gales and exceptional river flow rates (Slingo et al. 2014).Furthermore, a lack of investment in river dredging, flood defence schemes, intensive agricultural practice and inappropriate development of floodplains are thought to have exacerbated levels of flooding (Thorne 2014).The impact of the inundation was sufficient to invoke activation of the International Charter (International Charter 2000), a mechanism for making satellite data available for disaster response, on 4/12/13, 6/1/14 and 6/2/14.
The study area is centred on the city of Oxford, shown in Fig. 1.The area encompasses 154.35 km 2 .Multiple watercourses flow through and converge in the city and its surroundings.The Thames (Isis) and Cherwell form the major river systems with many smaller streams and brooks as tributaries.The landscape is predominantly made up of urban and agricultural uses with small tracts of mixed forest cover.Oxford city is thought to be primarily at risk of fluvial flood events (OCC 2011).Groundwater flooding is also known to be a risk factor (OCC 2011;Macdonald et al. 2012).In the city, 4500 properties

Methodology
The main processing steps designed for estimating and integrating evidence regarding flood extent are illustrated in Fig. 2. The methodology describes integration of near-realtime Landsat 8 imagery and geotagged photographs from Flickr, together with ancillary topographic data (derived from a LIDAR DTM).These datasets are processed and classified (see Sects.4.1, 4.2 and 4.3), and the results are then integrated using the fusion model (see Sect. 4.4).
Figure 3 shows a 1-month period detailing the acquisition times for the different data sources used.The period details the number of geotagged images captured each day and the Landsat 8 acquisition date (7 January).The figure also shows the river level for a gauge station (Kings Mill) close to the centre of the city of Oxford, which notably reaches a maximum height on the day of the Landsat 8 overpass.As discussed below, our experiments focus on the 7-day period delineated by the shaded area (5 Jan 2014-11 Jan 2014).

Crowdsourced flood extent estimation
In this work, the photograph-sharing site Flickr was used as a source of volunteered photographs.Flickr data were retrieved through the public API.For this study, the query was restricted to geotagged imagery falling within the study area and time frame (5 Jan 2014-11 Jan 2014).Flickr has fields for title, description and tags-each of which may contain information regarding the image content.Therefore, the keyword ''flood'' was used to detect relevant records across any of these fields (this term functions as a wildcard, detecting words as, for example, ''flooding'', ''floods'').This retrieved a total of 205 images taken by 25 different users.Figure 4 (left) shows the distribution of image locations across the study area.
Several issues regarding the quality of information are present when using geotagged web-harvested images.To further analyse the sample, the exchangeable image file (EXIF) metadata associated with each image were retrieved.The EXIF provides metadata regarding each individual image which includes information populated by the camera sensor or device used for capture.This includes potentially valuable geographic data such as the direction the device was facing during capture (GPSImageDirection).However, as noted in previous work, such as Panteras et al. (2014), these EXIF tags are not often populated.For the sample of 205 geotagged images that were collected, only 4 were found to contain EXIF data specifying the image direction.
To address the lack of metadata and derive a 2D representation of flood extent, an appropriate strategy for modelling these images is required.An analysis of the spatial point pattern of geotagged positions can be undertaken using kernel density estimation (KDE), as  adopted in Schnebele and Cervone (2013).These unsupervised methods learn a probabilistic surface based on the 2D point pattern.Rather than use a kernel-based method, we take a similar approach to Panteras et al. (2014) and adopt a viewshed-based model.Specifically, the cumulative viewshed is used, first adopted for archaeological analysis (Wheatley 1995).Here a viewshed for each photograph is calculated with each pixel labelled as visible (value = 1) or not visible (value = 0).A summation over all viewsheds therefore gives the number of times a location is visible from any photograph.This avoids the necessity of selecting a bandwidth parameter, as required by kernel density estimation, and takes account of topography-an important consideration in flood scenarios.The viewshed methodology provides a measure of corroboration between individual geotagged photographs at the pixel level and thus provides a model of agreement on where possible flooded areas may be.The viewshed methodology implemented here assumes a viewer height of 1.6 m, 360°viewing angle and a maximum viewing distance of 500 m.Although the focal length of the camera image might be utilised, this information was not present in the EXIF data.Thus, we chose a maximum viewing distance of 500 m based upon the range that would be effectively observable (e.g. a photograph taken at edge of a flooded field could extend to this range).The LIDAR DTM (described in further detail below) is used by viewshed algorithm for defining topography.The result of the cumulative viewshed analysis is shown in Fig. 4 (right).

Landsat 8 water detection
Both multispectral and SAR satellite remote sensing is commonly used for assessing inundation extent during flood events.The operational land imager (OLI) instrument of Landsat 8 provides multispectral data at 30-m spatial resolution.A single Landsat 8-OLI image was obtained from the US Geological Survey (USGS) Global Visualization Viewer as Level 1 Terrain Correct (L1T) product.One scene acquired on 7th January 2014 (path/ row: 203, 24) was obtained with cloud-free coverage of the study area.A small section in the south-eastern corner of the study area fell outside the acquisition swath and thus is missing data.Digital number (DN) values for the scene were converted to surface reflectance (see Fig. 5, left) using the dark object subtraction (DOS) method (Chavez Jr 1996) available in QGIS.Further data pre-processing steps undertaken here involved stacking Landsat bands 1-7 (30 m), re-projecting to British National Grid, cropping to the study area undertaken using GDAL.
One frequently applied approach for detecting water given spectral response at appropriate infrared wavelengths is to compute the Normalised Difference Water Index (NDWI) (McFeeters 1996).However, using the Modified NDWI (MNDWI) has the advantage of suppressing the response of both vegetation and built-up areas leading to enhanced water detection for these areas (Xu 2006).This is applicable to the landscape found in the Oxford study area and as such, we adopt this method here (see Fig. 5, right).MNDWI can be expressed as: where MIR corresponds to band 7 (SWIR2, 2.11-2.29 l) of the L8-OLI instrument.As an index method, a threshold for discriminating between water and non-water areas is required.The threshold value for identifying areas of water can be assessed by expert or automatically determined using methods such as Otsu's algorithm (Otsu 1979).The method uses the maximum between-class variance to determine the threshold.For further information, a detailed discussion of using Otsu's method for delineating water features is in Li et al. (2013).

Topographic flood factors
The topography of an area can greatly affect the flow and accumulation of floodwater.In general, the use of ancillary topographic data for aiding flood monitoring is widely applied for aiding land and flood mapping procedures.For example, Pierdicca et al. (2008) developed a method for estimating flooded areas from SAR with additional data from DEM and land cover mapping.
Here we build a digital terrain model (DTM) from elevation data sourced from both the UK Government Environment Agency and National Mapping Agency of Great Britain (Ordnance Survey).Environment Agency LIDAR DTM tiles at 2-m spatial resolution covering the study area were collated (Environment Agency 2015).As the LIDAR data did not cover the edges of the study area, Ordnance Survey Terrain 5 grid tiles were also acquired providing DTM data at a 5-m spatial resolution.These sources were merged and resampled using cubic convolution resampling to form a single elevation model at 5-m spatial resolution (see Fig. 6, right).The slope variable was calculated from the DTM using the GDAL (GDAL 2016) slope function (see Fig. 6, left).

Data fusion
A dataset for training and testing the statistical model, described in detail below, was obtained using a randomly distributed sample of 120 points drawn from a ground-truth inundation extent of 11th January 2014.The ground-truth flood inundation was created by digitisation of the extent according to a visual interpretation of colour infrared (CIR) aerial imagery undertaken by experts in the Environment Agency Geomatics group.This sample was partitioned into sets of 100 and 20 points for training and testing, respectively (Pradhan et al. 2010;Tehrany et al. 2014b).See Fig. 1 for the locations of training and testing points.

Weights of evidence
A Bayesian probabilistic method is adopted to establish the relation between the evidence variables.Weights-of-evidence (WoE) analysis first gained widespread usage within the geosciences community for geological prospecting.The method has since been widely used for analysis of natural hazards including landslides (Regmi et al. 2010;Sterlacchini et al. 2011) and flood susceptibility (Tehrany et al. 2014b).Bonham-Carter (1994) provides detailed description of WoE, the main parts of which are described here.
The weight terms are calculated as where W þ i and W À i denote the positive and negative weights for evidence B. The terms PfBjAg=PfBjAg and PfBjAg=PfBjAg denote, respectively, the sufficiency and necessity ratios (also known as likelihood ratios) respectively.These are calculated from where A denotes the area of known occurrences of the phenomena based on the training sample, B denotes the area of an evidence variable, and N{} denotes the count of the grid cells.Intuitively, within this work, these equations calculate proportions for the presence and absence of flood given presence and absence of the binary evidence.
In WoE, continuous variables must be discretised to form multiclass evidence data.Once weights for the evidence data have been calculated, the maps may be integrated together.The standard approach to applying WoE is to then convert multiclass evidence to binary evidence variables.While an index of susceptibility can be calculated based on a summation of contrast C (Oh and Lee 2010; Regmi et al. 2010), a map of posterior probability may be calculated through identification of the appropriate cut-off value determined from an analysis of contrast C (Bonham-Carter 1994).The ratio of the contrast to standard deviation, C/s(C), provides a statistical test for determination of the optimal cut-off point for each variable.Specifically, the threshold can be determined as the point where C/s(C) is greater than or equal to a defined confidence level.The two final weights for this binary map are assigned from the W? and W-values at this threshold point.
Assuming conditional independence between evidence variables (i.e.P(B 1 B 2 |-A) = P(B 1 |A)P(B 2 |A)), logarithmic posterior probability maps may be computed using where logit{A} denotes the prior logarithmic odds of a flooded location given its prior probability P(A).Therefore The prior probability p(A) can be calculated as the proportion of known flood locations to the study area.
where N{T} is the total area.Weight W i corresponds to either W þ i or W À i based on the presence or absence of the evidence, respectively.Then converting to a probability where P flood is the probability of a flooded location.
Computing posterior probability maps using weights of evidence requires an assumption of conditional independence between evidence variables.In geospatial problems, this assumption may be violated to some degree as spatial data can often exhibit correlations and therefore examination of the conditional independence should be undertaken (Bonham-Carter 1994;Agterberg and Cheng 2002) For discretisation in this work, a six-class quantile reclassification scheme was applied to variables with continuous values, i.e.Flickr cumulative viewshed, slope and elevation variables.Quantile classification was chosen due to its previous application to hazard mapping using evidence variables (Dickson et al. 2006;Tehrany et al. 2014b).As binary evidence, the MNDWI variable does not require reclassification.

Results
In this section, we present experimental results of applying the proposed workflow for estimating flood extent.First, we analyse the calculated weight and contrast values which quantify positive and negative correlations with flooded locations provided by each individual data source (Sect.5.1).Secondly, evaluation of the predictive power of the individual data sources and their integration together as a posterior probability model is made against ground-truth flood extent data (Sect.5.2).Thirdly, we generate a flood extent map based on the posterior probability (Sect.5.3).

Evidence weights
Table 1 details the weight calculations (W?, W-) for each evidence variable class.Figure 7 plots the weight contrast (C) for each of the multiclass data sources.Recall from Sect.4.4.1 that the contrast can be used to help assess the importance of each class and determine a binary cut-off threshold for the influence of a variable.
As described in Sect.4.4.1, the cumulative viewshed Flickr extent was classified using quantiles (Viewshed-01 to Viewshed-06).In Table 1 and Fig. 7 (left), the influence (W?) and weight contrast (C) of the viewshed extent remain negative for the first three classes (Viewshed-01 to Viewshed-03), with a strong positive correlation for classes [4][5][6].This corresponds to what we would expect given an overestimation of flooded extent provided by the viewshed methodology described in Sect.4.1.However, classes 4-6 relate to areas where the viewshed extent is corroborated by multiple images, correlating with flooded locations and lead to an increased positive weight (weight values of 1.72, 1.60, and 1.54).
Notably, the weight values for the Flickr cumulative viewshed are close to that found in Landsat-8 MNDWI water detection class, which has a strong positive weight (W?) of 2.11.Recall from Sect.4.2 that the water detection is based on a binary threshold derived using the Otsu method, a commonly used method for detecting flood presence.Given this, we can expect this detection to be influential in detecting flood presence.
The results also demonstrate correlations of topographic evidence variables with flood locations.As seen in Table 1, lower elevations exhibited a high positive correlation with flood locations, ranging from 1.42 to 0 for positive weight values (W?).The positive weight value for the slope variable is seen to provide a slightly weaker correlation ranging from 0.63 to 0.0 for gentle to steep slopes respectively.As discussed in Sect.4, analysis of the contrast (C) and C/s(C) values enables determination of an appropriate cut-off to determine the final evidence variable class weight (visual assessment of Fig. 7 (centre) and Fig. 7 (right) does not show a clear distinction for this cut-off).Here, a value of 2 is chosen for the confidence level (Bonham-Carter 1994).This calculation forms a Studentised measure of the certainty with which the contrast is known-i.e. a large value indicates an effect that is more likely to be real.Slope values covering values of 0°-0.86°andelevation covering altitudes of 28.63-58.81m defined this threshold, and the final evidence weight W fin , shown in Table 1.Here, ROC is used to assess both model fitting (Fig. 8a) and predictive power (Fig. 8b).The cumulative viewshed computed from geotagged Flickr imagery shows results of 0.6 and 0.55 for model fitting and testing, respectively.These results indicate that these data exhibit relatively poor sensitivity and specificity.This is partly to be expected as the assumption of a 360°viewing angle for each image in the viewshed calculation results in a maximum spatial extent, and this increases false positives in flood presence.Furthermore, the ROC model results are based on a relatively small split of randomly spatially distributed test data (20 points) which directly influence the results.Testing of the resulting probability map incorporating all data sources was carried out indicating no significant violations of conditional independence.A conditional independence ratio of 6.38 was calculated, well above the recommended upper threshold of 1 suggested by Bonham-Carter (1994).

Flood extent map generation
Using a cut-off point determined from the ROC analysis of the training data, binary flood extent maps may be generated from the probabilistic flood map.The cut-off point that equally maximised both sensitivity and specificity was chosen using the ROCR package (Sing et al. 2005).Figure 9 (right) shows the binary extent flood map against the CIR ground-truth extent.

Discussion
The results demonstrate that incorporation of multiple sources of data can aid prediction of flood extents using the proposed methodology.Despite this, several important assumptions must be noted.Firstly, this work does not consider temporal aspects of the data within the Fig. 9 Flood extent maps based on posterior probability (left) and probability cut-off calculated using ROC analysis modelling process as this case study involved a prolonged flooding event.Consideration of the temporal aspects would be important for generalising the work to flash flood mapping scenarios.Secondly, one drawback of the weights-of-evidence approach is the need for variable discretisation.This can result in information loss if sufficient consideration is not given to the modelling problem and evidence data used.This means that the proposed methodology cannot be applied in a completely automated fashion unless an automatic choice of number of classes is performed.
The work described here exploited a data-driven fusion model to integrate variables which inform the desired analytical output-in this case a map of flood extent.However, in other disaster scenarios alternative outputs might be more valuable.For example, identification of the areas where there is high uncertainty might be useful, allowing for further remote interrogation or field investigation.The weights-of-evidence approach supports this through enabling analysis of the weight variance (Dickson et al. 2006).Extension of the analysis of either the weights or uncertainty with ensemble learning methods also offers the potential of improving the results.
Although the experiments here exploit geotagged images sourced from social media, the rapid inundation mapping methodology employed could incorporate other crowdsourced assessments such as Twitter (Panteras et al. 2014;Smith et al. 2015) or volunteered by citizens (Poser et al. 2009).Similarly, other authoritative data sources may be relevant depending on the context.Within the UK, for example, many large watercourses are gauged.Where this is the case, both river flow measurements and river level data are usually available and this might be incorporated through linear interpolation of depth levels as an indication of maximal water extent (Apel et al. 2009) or inclusion of an appropriate 2D hydrological inundation model (Liu et al. 2015).Other remote sensing platforms may also be utilised as independent sources of evidence, such as SAR (Mason et al. 2010;Matgen et al. 2011;Schumann et al. 2013).
Data quality is often discussed in relation to use of crowdsourced information.In terms of rapid inundation mapping, the data sources expose different levels of uncertainty.On one side with authoritative data, this uncertainty comes from a controlled quality assurance but which needs to be modified due to the assessment time difference.On the other side, various crowdsourcing data (e.g.Flickr, citizen science, Twitter) come with less controlled quality assurance, but are timely.This reflects a recognised issue in the use of volunteered data for natural hazard assessment where data quality is a concern, but its adoption is particularly desirable (Goodchild and Glennon 2010;Haworth 2016).Our work evaluated a sample of social media data taken for a large and serious flooding event retrieved via keyword search.Despite this, the small size of the retrieved social media, and relatively poor performance of this information when used on its own, highlights that crowdsourced data currently must be used in combination with traditional flood mapping methods.However, novel additional techniques for identifying relevant social media reports and increasing samples sizes of digital reporting from citizens may relax this requirement (as long as its quality can be assured).The fusion method proposed did not take into account the quality of the information, but the workflow proposed here could integrate a quality assurance step (e.g.such as a probabilistic credibility assessment (Hung et al. 2016)), even if the algorithms used in quality assurance and in data fusion present some entanglement (Leibovici et al. 2015).This paper presented the methodology and experimental results of a rapid flood mapping workflow using social media, remote sensing and topographic map data.Using a Bayesian statistical model based on weights-of-evidence analysis, a data-driven approach was developed which quantifies the contribution of each information source.Combining these sources together enables generation of a probability map denoting likelihood of the presence of floodwater at the pixel level.A case study using data relating to the 2014 UK flood event was completed.ROC assessment and AUROC values demonstrated that the multisource flood mapping method effectively predicts spatial extent of water.Furthermore, the methodology shows that even relatively simple voting-based procedures for controlling for the quality of non-authoritative social media data may be integrated in a flood assessment workflow.In future work, we will investigate the effects of uncertainty within the data, such as positional accuracy of the geotagged photographs and the imprecision of viewing direction, on the predicted flood extent.

Fig. 1
Fig. 1 Study area location within Great Britain, ground-truth flood extent, training and testing points

Fig. 2
Fig. 2 Processing methodology illustrating data sources, analytical steps and validation

Fig. 3
Fig. 3 Data sources and availability including number of geotagged images (Flickr), Landsat-8 acquisition time, CIR aerial photograph (used for ground-truth) acquisition time and river level height (at the Kings Mill gauge).The shaded area corresponds to the experimental time frame for data analysis

5. 2
Model testing and validationReceiver operating characteristic (ROC) and area under curve ROC (AUROC) analysis was used to assess the performance of the individual data sources derived from Landsat-8 MNDWI and Flickr analysis, and the posterior probability model generated from combining all sources using weights of evidence.In brief, ROC assessments are typically used to measure the predictive power of a model to detect a binary outcome and are often used for validating susceptibility models and mapping tasks (such asDevkota et al. (2013) andTehrany et al. (2014a)).The ROC curve plots the false positive rate (1-specificity) against the true positive rate (sensitivity).The sensitivity and specificity provide an evaluation of the model at correctly identifying flooded locations (true positive rate) and non-flooded locations (true negative rate).