applied sciences
Article Geostatistical Analysis of the Spatial Correlation between Territorial Anthropization and Flooding Vulnerability: Application to the DANA Phenomenon in a Mediterranean Watershed
Salvador Garcia-Ayllon 1,2,* and John Radke 2,3
1 Department of Civil Engineering, Technical University of Cartagena, 30009 Cartagena, Spain 2 Department of Landscape Architecture and Environmental Planning, University of California, Berkeley, CA 94720, USA; [email protected] 3 Department of City and Regional Planning, University of California, Berkeley, CA 94720, USA * Correspondence: [email protected]
Abstract: Climate change is making intense DANA (depresión aislada en niveles altos) type rains a more frequent phenomenon in Mediterranean basins. This trend, combined with the transformation of the territory derived from diffuse anthropization processes, has created an explosive cocktail for many coastal towns due to flooding events. To evaluate this problem and the impact of its main guiding parameters, a geostatistical analysis of the territory based on GIS indicators and an NDVI (Normalized Difference Vegetation Index) analysis is developed. The assessment of the validity of a proposed methodology is applied to the case study of the Campo de Cartagena watershed located
around the Mar Menor, a Mediterranean coastal lagoon in Southeastern Spain. This area has suffered three catastrophic floods derived from the DANA phenomenon between 2016 and 2019. The results
Citation: Garcia-Ayllon, S.; Radke, J. show that apart from the effects derived from climate change, the real issue that amplifies the damage Geostatistical Analysis of the Spatial caused by floods is the diffuse anthropization process in the area, which has caused the loss of the Correlation between Territorial natural hydrographic network that traditionally existed in the basin. Anthropization and Flooding Vulnerability: Application to the Keywords: flooding vulnerability; DANA; Mar Menor; geostatistical analysis; Campo de Cartagena; DANA Phenomenon in a anthropization Mediterranean Watershed. Appl. Sci. 2021, 11, 809. https://doi.org/ 10.3390/app11020809 1. Introduction Received: 23 December 2020 The evaluation of the problems associated with floods has typically been fundamen- Accepted: 14 January 2021 tally approached from a hydrological perspective [1,2], although several points of view can Published: 16 January 2021 be found in the scientific literature. However, in the field of flood risk management, apart from traditional methodological approaches based on the simulation of flood models [3,4], Publisher’s Note: MDPI stays neutral there are other perspectives based on multicriteria socioeconomic assessments [5], statistical with regard to jurisdictional claims in published maps and institutional affil- rainfall indicators [6], or several different natural resilience strategies [7–10], for example. iations. Currently, the increasing effect of climate change on the phenomena associated with torrential rains is forcing a rethink in those methodologies and approaches in many places [4,11,12]. One of the locations where this change has been particularly virulent is in Spanish Mediterranean basins. There, the appearance of the DANA phenomenon (Spanish acronym for depresión aislada en niveles altos, meaning upper-level isolated at- Copyright: © 2021 by the authors. mospheric depression) has replaced the traditional flash floods associated with cold fall Licensee MDPI, Basel, Switzerland. phenomena [13]. Technically, it is a depression isolated at high levels, where a cold front This article is an open access article distributed under the terms and clashes with warm air and produces heavy storms and torrential rains. In practice, DANA conditions of the Creative Commons is a weather phenomenon quite similar to a cold drop, but its intensity and effects are Attribution (CC BY) license (https:// proving to be much more devastating. Some researchers have already begun to warn creativecommons.org/licenses/by/ Europe about the growing phenomenon of “medicanes” [14–16]. These meteorological 4.0/).
Appl. Sci. 2021, 11, 809. https://doi.org/10.3390/app11020809 https://www.mdpi.com/journal/applsci Appl. Sci. 2021, 11, x FOR PEER REVIEW 2 of 22
Appl.Appl. Sci. Sci.2021 2021, 11, 11, 809, x FOR PEER REVIEW 22 of of 21 22
after the analogy in their effects with the usual hurricanes in the USA, are increasingly phenomena,afterpresent the analogyevery named year in after theirin the the effects Mediterranean analogy with in the their dueusual effects to hurricanes climate with the change usualin the [17] hurricanes USA,. are increasingly in the USA, arepresent increasinglyMediterranean every year present in thecoastal every Mediterranean areas year inhave the traditionally Mediterraneandue to climate been change due characterized to climate[17]. change as being [17 subject]. to frequentMediterraneanMediterranean flooding, coastal coastalsome of areas areas which have have are traditionally traditionally even quite been relevant.been characterized characterized However, as asduring being being subjectthe subject last to twoto frequentfrequentdecades, flooding,flooding, relevant somesome flooding ofof whichwhich events areare and even even damage quite quite relevant. relevant.have increased However,However, exponentially duringduring thethe last( Figlasture two two 1) . decades,decades,This phenomenon, relevantrelevant floodingflooding which events eventsoccurs andinand most damagedamage of the havehave Mediterranean increasedincreased exponentially exponentiallybasins, is likely (Figure(Fig influencedure1 ).1). ThisThisby phenomenon, phenomenon,the effects of whichclimate which occurs occurschange, in in most asmost some of of the therecent Mediterranean Mediterranean studies have basins, basins, begun is is to likely likely notice influenced influenced [18,19]. In bybygeneral, thethe effectseffects the of ofexisting climate climate predictions change, change, asas in some some the recent meteorologicalrecent studies studies have havestudies begun begun carried to to notice notice out in [ 18[18,19]. Mediterra-,19]. In In general,general,nean basins the the existing existing such as predictions predictionsthis one [20,21] in in the the show meteorological meteorological a global trend studies studies towards carried carried desertification out inout Mediterranean in Mediterra- in the area basinsneandue basins to such the as annualsuch this as one lossthis [20 ofone, 21rainfall. [20,21]] show This show a global is ain global contrast trend trend towards to an towards increasing desertification desertification value in of the the in area maximum the duearea toduedaily the to annual the rains annual that loss denote loss of rainfall. of anrainfall. increasing This This is in is trend contrastin contrast of these to to an extr an increasing increasingeme episodes value value (Fig of of theure the maximum2). maximum However, dailydailythese rains rains meteorological thatthat denotedenote an ananalyses, increasingincreasing which trendtrend are ofof carried thesethese extremeextr outeme by episodestheepisodes administrations (Figure(Figure2 ).2). However, However,responsible thesethesefor the meteorologicalmeteorological matter and which analyses,analyses, also whichunderlinewhich areare carried thecarried statistical out out by by increases the the administrations administrations in frequency responsible responsibleand intensity forfordetected the the matter matter in and rainfalland which which for also also recent underlineunderline years, thethe do statisticalstatisticalnot justify increasesincreases the increase in in frequencyfrequency in damage and and intensity fromintensity cata- detecteddetectedstrophic in in episodes rainfall rainfall for that for recent appeared recent years, years, in do the not do historical justifynot justify the series increase the. increase in damage in damage from catastrophic from cata- episodesstrophic thatepisodes appeared that appeared in the historical in the historical series. series.
Figure 1. Impact level of flooding events between 1950 and 2020 in the Campo de Cartagena basin. Source: AEMET and Figure 1. Impact level of flooding events between 1950 and 2020 in the Campo de Cartagena basin. Source: AEMET and FigConsorcioure 1. Impact de Compensación level of flooding de segurosevents between de España 1950 [22,23]. and 2020 in the Campo de Cartagena basin. Source: AEMET and ConsorcioConsorcio de de Compensaci Compensaciónón de de seguros seguros de de España España [ 22[22,23].,23].
(a) (b) (a) (b) FigureFigure 2. Climate2. Climate change change scenarios scenarios for for 2020–2100 2020–2100 regarding regarding precipitation precipitation in in the the Region Region of Murciaof Murcia area area using using multiple multiple regressionFigregressionure 2. models:Climate models: annualchange annual rainfallscenarios rainfall (a) (for anda) and2020 maximum maximum–2100 regarding daily daily rainfall rainfall precipitation (b). (b Source:). Source: in AEMETthe AEMET Region [22 [22].]. of Murcia area using multiple regression models: annual rainfall (a) and maximum daily rainfall (b). Source: AEMET [22]. A commonA common variable variable to to all all these these Mediterranean Mediterranean basins basins is theis the intense intense anthropization anthropization processprocessA common that that these these variable territories territories to all have have these undergone undergone Mediterranean [24 [24,25],25]. basins. The The Spanish Spanishis the intense Mediterranean Mediterranean anthropization coastcoast is isprocess undoubtedly that these the territories European have region undergone in which [24,25] activities. The such such Spanish as as tourism, tourism, Mediterranean agriculture, agriculture, coast infra- in- is frastructureundoubtedly expansion, the European etc., region have experiencedin which activities the greatest such as degree tourism, of growthagriculture, in recent infra- decades [26]. Nevertheless, it has traditionally been difficult to analyze the extent to which Appl. Sci. 2021, 11, x FOR PEER REVIEW 3 of 22
Appl. Sci. 2021, 11, 809 structure expansion, etc., have experienced the greatest degree of growth in recent3 dec- of 21 ades [26]. Nevertheless, it has traditionally been difficult to analyze the extent to which phenomena such as diffuse anthropization processes are responsible for the increasing vulnerability to the flooding of an area. Complex phenomena, such as the effect of soil phenomena such as diffuse anthropization processes are responsible for the increasing sealingvulnerability due to to the the artificialization flooding of an area.of theComplex territory, phenomena, the alteration such of theas the natural effect hydro- of soil graphicsealing duenetwork to the of artificialization the basins due ofto thechanges territory, in land the alterationuse, or the of impact the natural of the hydrographicbarrier effect generatednetwork of by the linear basins transport due to changesinfrastructures, in land use,have or no the wide impact-ranging of the methodological barrier effect gener- pro- posalsated by in linear the current transport scientific infrastructures, literature. haveSeveral no wide-rangingapproaches for methodological specific case studies proposals can bein thefound current [27–29] scientific, but comprehensive literature. Several large- approachesscale evaluatio forn specific methods case from studies a territorial can be pointfound of [27 view–29], are but still comprehensive lacking. large-scale evaluation methods from a territorial point of viewIn this are context, still lacking. this research proposes an analytical method based on the evaluation of GISIn indicators this context, using this geostatistical research proposes and Normalized an analytical Difference method basedVegetation on the Index evaluation (NDVI of) algorithms.GIS indicators This using innovative geostatistical approach and enables Normalized the statistical Difference and Vegetationspatial correlation Index (NDVI of the) degreealgorithms. of territorial This innovative transformation approach derived enables from the anthropic statistical processes and spatial with correlation the increasing of the vulnerabilitydegree of territorial of a territory transformation to floods. derived To assess from its anthropicvalidity, the processes methodology with the is applied increasing to thevulnerability Campo de of Cartagena a territory tocase floods. study, To a assess large its basin validity, located the methodologyin the southeast is applied part of to the SpanishCampo deMediterranean Cartagena case coast study, that a has large been basin suffering located infrom the catas southeasttrophic part floods of the in Spanish recent years.Mediterranean coast that has been suffering from catastrophic floods in recent years.
2.2. Materials Materials and and Methods Methods 2.1.2.1. Area Area of of Study Study TheThe study study area area is the Campo de Cartagena watershed. This This basin basin is is located in the RegionRegion of of Murcia, Murcia, a a semiarid semiarid Mediterranean Mediterranean territory territory in in Southeastern Southeastern Spain. Spain. The The area area is characterizedis characterized by scarce by scarce precipitation precipitation (<300 (<300 mm per mm year) per year)that mainly that mainly occurs occurs during during storm eventsstorm eventsin autumn in autumn and winter. and winter. This traditional This traditional dryland dryland agriculture agriculture territory territory has been has been ex- tensivelyextensively transformed transformed from from a relatively a relatively natural natural landscape landscape to to one one of of intense intense human human activity activity inin the the last last 60 60 years years despite despite having having an important important environmental environmental value [30]. [30]. A Att present, present, we can find a varied catalogue of human activities and anthropic elements, such as high-density can find a varied catalogue of human activities and anthropic elements, such as high-den- urban sprawl in coastal towns, isolated tourist resorts, intensive irrigated agriculture, sity urban sprawl in coastal towns, isolated tourist resorts, intensive irrigated agriculture, important infrastructures such as highways, railways, ports or airports, etc. (Figure3). A important infrastructures such as highways, railways, ports or airports, etc. (Figure 3). A coastal lagoon called Mar Menor is situated at the end of this large watershed, delimited coastal lagoon called Mar Menor is situated at the end of this large watershed, delimited by a group of mountains configuring an extensive plain of about 1440 km2. by a group of mountains configuring an extensive plain of about 1440 km2.
Figure 3. Area of study: Campo de Cartagena watershed surrounding the Mar Menor coastal lagoon. Figure 3. Area of study: Campo de Cartagena watershed surrounding the Mar Menor coastal lagoon.
Freshwater that inputs into the lagoon were traditionally restricted to six main ephemeral watercourses called wadis. These wadis can reach lengths of over 50 km. The Appl. Sci. 2021, 11, x FOR PEER REVIEW 4 of 22 Appl. Sci. 2021, 11, x FOR PEER REVIEW 4 of 22
Appl. Sci. 2021, 11, 809 4 of 21 Freshwater that inputs into the lagoon were traditionally restricted to six main Freshwater that inputs into the lagoon were traditionally restricted to six main ephemeral watercourses called wadis. These wadis can reach lengths of over 50 km. The ephemeral watercourses called wadis. These wadis can reach lengths of over 50 km. The surface area of the watershed nourishing some of these wadis, such as the Albujon, even surfacesurface area area of of the watershed nourishing some of t thesehese wadis wadis,, such as the Albujon, even exceeds that of the Mar Menor lagoon itself (Figure 4). These wide and shallow gullies are exceedsexceeds that that of of the the Mar Mar Menor Menor lagoon lagoon itself itself (Fig (Figureure 4).4 ).These These wide wide and and shallow shallow gullies gullies are generally inactive but can carry great quantities of water and sediment during the tradi- generallyare generally inactive inactive but can but carry can carry great greatquantities quantities of water of water and sediment and sediment during during the tradi- the tional Mediterranean flood episodes in autumn. The torrential nature of the water supply tionaltraditional Mediterranean Mediterranean flood floodepisodes episodes in autumn. in autumn. The torrential The torrential nature natureof the water of the supply water is now aggravated by the impermeable soils and scarce vegetation cover of the watershed issupply now aggravated is now aggravated by the impermeable by the impermeable soils and scarce soils and vegetation scarce vegetation cover of the cover watershed of the areas, going directly to the Mar Menor. As stated in Figure 2, during the last two decades, areas,watershed going areas, directly going to directlythe Mar toMenor. the Mar As Menor. stated in As Fig statedure 2, in during Figure2 the, during last two the decades, last two damage from relevant flooding events has been catastrophic and has increased by 318% damagedecades, from damage relevant from relevantfloodingflooding events has events been has catastrophic been catastrophic and has and increased has increased by 318% by [22]. Among them, the DANAs from 2016 and 2019 are especially serious (Figure 5). [22].318% Among [22]. Among them, them,the DANAs the DANAs from 2016 from and 2016 2019 and are 2019 especially are especially serious serious (Figure (Figure 5). 5).
FigFigureure 4. 4. ((aa)) Wadis Wadis going going to the Mar Menor lagoon lagoon.. ( (b) Watershed nourishing Albujon wadi. Figure 4. (a) Wadis going to the Mar Menor lagoon. (b) Watershed nourishing Albujon wadi.
(a) (b) (a) (b) Figure 5. (a) Flooding of the coastal town of Los Alcazares after the DANA (depresión aislada en niveles altos) in 2016 (b) Figure 5. (a) Flooding of the coastal town of Los Alcazares after the DANA (depresión aislada en niveles altos) in 2016 (b) SentinelFigure 5. 2 (satellitea) Flooding image of of the the coastal effects town of the of 2019 Los AlcazaresDANA in afterthe Mar the Menor. DANA Source: (depresi óEmergencyn aislada en service niveles altosof the) inRegion 2016 (ofb) Sentinel 2 satellite image of the effects of the 2019 DANA in the Mar Menor. Source: Emergency service of the Region of Murcia.Sentinel 2 satellite image of the effects of the 2019 DANA in the Mar Menor. Source: Emergency service of the Region Murcia. of Murcia. 2.2. Methodological Overview 2.2.2.2. Methodo Methodologicallogical Overview Overview The process of analyzing the impact of diffuse territorial anthropization in increasing The process of analyzing the impact of diffuse territorial anthropization in increasing the vulnerabilityThe process of a analyzing territory theof floods impact caused of diffuse by DANA territorial is carried anthropization out in applied in increasing stages the vulnerability of a territory of floods caused by DANA is carried out in applied stages inthe the vulnerability following framework. of a territory of floods caused by DANA is carried out in applied stages in the following framework. in theIn following the first stage, framework. an analysis of contrasting the expected results from the application In the first stage, an analysis of contrasting the expected results from the application of simulationsIn the first following stage, an analysisthe models of contrasting collected from the expected the official results regulations from the for application a highly ofof simulations following following the the models models collected collected from from the the official official regulations for a highly anthropized study area at the territorial level is initiated. In this case, simulations are carried out in the Mar Menor and Campo de Cartagena areas for different return periods using Appl. Sci. 2021, 11, 809 5 of 21
the national system of cartography of flooded areas by the Ministry of the Environment in Spain. The theoretical results obtained for the study area are contrasted with the results observed through an NDVI analysis of the latest relevant extreme events (in this case the DANAs of 2016, 2018, and 2019). Second, in this case study, a relevant damage situation is examined, as an unexpected dissociation between the theoretical floodplain and those observed in the last events is detected. As a result, this situation is analyzed to determine what the implications of the anthropic territorial transformations are in the area. To accomplish this, a territorial analysis of the study area is performed, as executed in the previous subsection, to determine what the possible anthropic phenomena may be that are influenced by the global phenomenon of the territorial vulnerability to floods. Third, based on this territorial analysis, a series of GIS indicators of anthropic phe- nomena and a GIS indicator of the damages derived from floods are elaborated. The first indicators must spatially represent the different levels of damage caused. The latter must be able to model, in a significant way, the spatiotemporal evolution of the anthropic phenomena estimated as more related to the problem of flooding in the chosen study area. Fourth, and lastly, a geostatistical analysis is carried out to assess the spatial correlation of the evolution of the selected indicators of territorial anthropization to the generation of damages derived from floods. To accomplish this, we first calculate the global Moran’s I autocorrelation statistic to verify that the spatial distribution patterns obtained for the said correlation have statistical significance. The verification would show that we are facing a real phenomenon and not a random distribution without physical significance. Secondly, a bivariate analysis is performed using the Anselin Local Moran I statistic to assess numerically the individual correlation level of each of the different anthropization indicators with the damage indicator. Finally, we apply a hot spot analysis where spatial statistics indicate the behavior patterns of the spatial distribution of the previous statistic. We next describe in detail the methodological process of these stages.
2.3. Elaboration of GIS Indicators Two georeferenced databases are implemented in this model to generate the indicators that serve as evaluation elements in the geostatistical analysis: (a) spatial values of the damage that occurred in the last three DANAs and (b) a spatiotemporal distribution of the most relevant GIS transformation indicators found in the bibliography.
2.3.1. Spatial Damage Values Database To analyze the spatial distribution of damage caused by DANAs, we created an index called the flood damage severity index (IFDS). From the information provided by the emergency services from various local and regional administrations that are in charge of processing the damage files of those affected by the flooding events, we generated a spatial qualitative punctual database of damage (Figure6). These data were obtained on aggregate units such as tourist resorts, residential buildings, industrial estates, shopping centers, etc., in order to make georeferenced treatment possible whilst preserving the legal requirements of anonymity (technical details of the database configuration can be found in AppendixA ). To obtain a uniform sequence of discrete values, the alphanumeric data are classified into three categories based on the level of significance of the damage: minor, relevant, and catastrophic damage, following the criteria from Table1. Appl. Sci. 2021, 11, 809 6 of 21 Appl. Sci. 2021, 11, x FOR PEER REVIEW 6 of 22
FigureFigure 6.6. DetailDetail ofof spatial spatial distribution distribution of of the the georeferenced georeferenced database database generated generated for for damage. damage. Source Source of data:of data: Own Own elaboration elabora- fromtion from several several emergency emergency services services of various of various municipalities, municipalities, regional regional and state and emergencystate emergency service service [31,32 [31]. ,32].
Table 1. 퐼 criteria for the level of significance of the damage. Table 1. IFDS퐹퐷푆criteria for the level of significance of the damage.
MinorMinor DamageDamage Relevant Relevant DamageDamage Catastrophic Catastrophic Damage Damage Block on highways or critical Block on highways or critical ImpactImpact on communicationscommunications TemporaryTemporary pathpath blocks blocks Temporary Temporary roadblocksroadblocks roadroad infrastructure Structural destruction or Structural destruction or DamageDamage to buildingsbuildings and and FloodingFlooding of garages of garages in build- FloodFlood inside inside homes homes and and main completecomplete loss ofof the the value value of of materialmaterial goodsgoods in buildingsings mainareas areas of buildingsof buildings materialmaterial goodsgoods TreeTree uprooting uprooting or or drowning drowning AgricultureAgriculture andand livestock livestock PartialPartial crop damagedamage Complete Complete lossloss ofof the the crop crop of animalsanimals TemporaryTemporary outages outages of ofessen- PermanentPermanent or orprolonged prolonged out- MalfunctionMalfunction ofof nonessentialnonessential OperationOperation of publicpublic services services tialessential supply supply services services (water, agesoutages of essential of essential supply supply ser- servicesservices (e.g.,(e.g., streetlamps)streetlamps) (water,electricity) electricity) servicesvices (water, (water, electricity) electricity) Damage to unprotected areas Recoverable damage to Damage to unprotected areas Recoverable damage to pro- Irreversible damage in Environmental damage with no environmental, protected areas (e.g., loss of Irreversible damage in pro- Environmental damage with no environmental, his- tected areas (e.g., loss of the protected areas historical, or cultural value the beach line) tected areas torical, or cultural value beach line)
2.3.2.2.3.2. GISGIS IndicatorsIndicators ofof AnthropizationAnthropization TheThe phenomenon of of anthropization anthropization is is spatially spatially evaluated evaluated using using GIS GIS indicators indicators of ter- of territorialritorial transformation. transformation. This This spatial spatial transformation transformation generates generates subphenomena subphenomena such asas thethe “dam“dam effect” effect” and and “soil “soil sealing” sealing” or or alterations alterations in in the the orography orography of of the the land land which, which according, accord- toing the to the bibliography bibliography consulted, consulted, may may have have important important implications implications on on the the risk risk of of flooding flooding ofof aa territory.territory. ByBy usingusing historichistoric GISGIS cartography,cartography, thethe evolutionevolution ofof variousvarious dimensionlessdimensionless indicatorsindicators associatedassociated withwith thesethese subphenomenasubphenomena isis evaluatedevaluated numericallynumerically overover time.time. TheThe indicatorsindicators selectedselected toto evaluateevaluate thethe patternspatterns ofof territorialterritorial transformationtransformation inin thethe studystudy areaarea areare detaileddetailed below:below: • LinearLinear infrastructureinfrastructure densitydensity ((LID)LID) index:index: LinearLinear communicationcommunication infrastructuresinfrastructures areare largelylargely associatedassociated with with the the generation generation of of “ “damdam effects effects”” during during episodes episodes of oftorrential torren- tialrain rain that that generate generate floods floods [33,34] [33,.34 This]. This index index is calculated is calculated over over time time to show to show the theim- importance of the impact produced by these infrastructures within the whole of the portance of the impact produced by these infrastructures within the whole of the phenomenon of the anthropization of a territory. This dimensionless index evaluates phenomenon of the anthropization of a territory. This dimensionless index evaluates Appl. Sci. 2021, 11, 809 7 of 21
the spatial presence of these infrastructures in each analysis spatial cell by taking the size and relevance of the linear infrastructure into consideration.
∑ h L 2 LID = i i (1) Str
Li = length of existing linear infrastructures (m) hi = weighting coefficient (highway = 1, normal road = 0.75, country road = 0.5) 2 Str = territorial surface in reference (m ) • Index of soil artificialization (SA): The urbanization of the territory is, together with erosion, the main cause of the “sealing effect” of the soil [35,36]. This phenomenon notably reduces the soil’s ability to absorb rainwater. To analyze this phenomenon, the artificial transformation process of the territory is evaluated over the last decades by applying the CORINE Land Cover criteria [37]. The indicator differentiates areas of low permeability (such as urban areas without green areas) from artificial areas with absorption capacity (such as golf resorts). This dimensionless index evaluates land artificial transformation in each analysis cell by taking the size of the linear infrastructure into consideration. ∑ h A SA = i i (2) Str
2 Ai = artificial surface area (m ) hi = weighting coefficient (highly waterproof artificial surface = 1, medium waterproof artificial surface = 0.75, low waterproof artificial surface = 0.5) 2 Str = territorial surface in reference (m ) • Index of alteration of the terrain orography (ATO): The transformations in the orog- raphy of the land are also a parameter strongly implicated in the flood risks of the territory [38]. Apart from the cases mentioned in the two previous phenomena (infras- tructural and building construction works), agricultural activity is the most common and relevant cause, in terms of the affected surface area, of these land alterations. In this section, aggressive transformations such as the replacement of trees or herbaceous crops by fruit and vegetable crops are differentiated from intermediate actions, such as relevant pending changes or elimination of the terraced structures of the crops, and minor ones, such as simple changes in the crop direction or land use without three-dimensional alteration of its orography.
∑ h A ATO = i i (3) Str
2 Ai = orographically altered surface (m ) hi = weighting coefficient (relevant orographic alterations = 1, intermediate alterations = 0.75, simple land use change = 0.5) 2 Str = territorial surface in reference (m )
2.4. Computation of Indicators First, to determine the exact scope of the analysis, the current configuration at the official level of the structure of natural channels in the area was taken into account. Unregu- lated channels with a level greater than 3 according to the Horton–Strahler criteria [39] were selected, in order to generate a representative model of the current hydromorphometric situation of the theoretical natural hydrographic network in the Campo de Cartagena area (Figure7). Appl. Sci. 20202121, 11, x 809 FOR PEER REVIEW 8 of 22 21
FigFigureure 7. Area of study selected following a watershed model based on Horton–StrahlerHorton–Strahler criteria established by national regulations.
ToTo contrast the distribution of the real floodplain floodplain produced by the latest DANA in comparisoncomparison to thethe situationsituation expectedexpected by by the the current current theoretical theoretical models models of of the the regulations, regulations, an anNDVI NDVIanalysis analysis was was carried carried out out to identify to identify the changesthe changes in land in land cover cover during during the DANA the DANA days. days.The Normalized The Normalized Difference Difference Vegetation Vegetation Index (IndexNDVI )(NDVI is a land) is covera land analysis cover analysis index used index to usedestimate to estimate the quantity, the quantity, quality, quality, and development and development of vegetation of vegetation based on based the measurement on the meas- urementof the radiation of the radiation intensity ofintensity certain of bands certain of thebands electromagnetic of the electromagnetic spectrum thatspectrum vegetation that vegetationemits or reflects. emits Foror reflects. the calculation For the ofcalculation this vegetation of this index, vegetation the information index, the foundinformation in the foundred and in infraredthe red and bands infrared of this bands electromagnetic of this electromagnetic spectrum was spectrum used. The was calculation used. The of cal- the NDVI culationwas of the done NDVI using was the done following using formula: the following formula: 푁퐼푅 − 푅푒푑NIR − Red 푁퐷푉퐼 =NDVI = (4(4)) 푁퐼푅 + 푅푒푑NIR + Red where Red and NIR stand for the spectral reflectance reflectance measurements acquired in the red (visible)(visible) and and near near-infrared-infrared regions, respectively. This method allows us to identify all the areasareas in which which the the rains rains have have produced produced a flood flood plate on the surface surface of the the basin basin regardless regardless ofof whether whether or or not not damage damage has has occurred occurred in in an an area. area. In In this this way, way, the areas subject to runoff becausebecause of of the the rains rains were identified. identified. The data data were obtained through the Sentinel 2A and 2B2B satellites satellites with with an an accuracy level of 10 m for the days in which the DANAs o off the year yearss 2016,2016, 2018 2018,, and 2019 occurred. Then,Then, the the GIS analysis wa wass undertaken using a constructed 5 mm22 digital elevationelevation mmodelodel (DEM) (DEM) of of the the Campo Campo de deCartagena Cartagena area area to ass toess assess the statistical the statistical correlation correlation between be- thetween anthropizing the anthropizing process process of the ofterritory the territory from the from case the study case studyand the and increase the increase in its invul- its nerabilityvulnerability to flooding. to flooding. Laser Laser imaging imaging detection detection and andranging ranging (LiDAR) (LiDAR) tools tools from from the IM- the IDAIMIDA LiDAR LiDAR SDI SDI [32] [32 are] are used used to toimplement implement relevant relevant elements elements such such as as urban urban areas, areas, wadi orographyorography,, and the slopes of the agriculturalagricultural crops inin thethe DEMDEM (Figure(Figure8 8).). Appl. Sci. 2021, 11, 809 9 of 21 Appl. Sci. 2021, 11, x FOR PEER REVIEW 9 of 22
FigFigureure 8 8.. DigitalDigital Elevation Elevation Model Model of of Campo Campo de de Cartagena Cartagena area: area: comparison comparison of of wadis wadis’’ (blue (blue lines) lines) origin origin and and finish. finish. Source Source ofof the the data: data: L LiDARiDAR SDI SDI tools tools from from IMIDA IMIDA Institute Institute of of the the Region Region of of Murcia Murcia [32] [32]..
OnceOnce the the distributions distributions at at the the spatial spatial level level of of the the flooding flooding impact impact and and territorial territorial an- an- thropizationthropization indices werewere generated,generated, we we evaluated evaluate thed the possible possible spatial spatial correlation correlation between be- tweenthem usingthem geostatisticalusing geostatistical tools. For tools. this, For the this, territory the territory needs to needs be “discretized” to be “discretize as a continu-d” as a continuousous element element in 100 × in100 100 meter × 100 cells meter in order cells toin numericallyorder to numerically evaluate the evaluate spatial the correlation spatial correlationbetween indicators between with indicators specific with values specific on a continuous values on a space. continuous These cells space. were These assigned cells werevalues assigned according values to the according results obtainedto the results in the obtainedNDVI analysisin the NDVI and analy the GISsis cartographyand the GIS cartographydatabase for database the damage for ofthe the damage DANAs. of Onthe theDANAs. other hand,On the anthropic other hand, values anthropic were obtained values weforre the obtained indicators for the of territorial indicators transformation of territorial transformation of historical GIS of historical cartography GIS forcartography the years for1956, the 1981, years and 1956, 2016. 1981, The and cartography 2016. The cartography detailed data detailed are summarized data are summarized in Table2. in Table 2. Table 2. Technical characteristics of the georeferenced data used. Table 2. Technical characteristics of the georeferenced data used. Pixel Size Projected on the Planimetric Accuracy (X,Y) Altimetric Accuracy (Z) Mapping Data GSD Ground (cm) Pixel Size Projected on the Planimetric Ac- Altimetric Ac-Mesh Step Mean Squared Error (m) Mean Squared Error (m) FlightMapping Orthophoto GSD Ground (cm) curacy (X,Y) curacy (Z) Mesh Data Mean Squared Mean Squared ×Step 1956–1981 45 50Flight <1.00Orthophoto <2.00 5 5 1982–2016 22 25 <0.50Error (m) <1.00 Error (m) 5 × 5 1956–1981 45 50 <1.00 <2.00 5 × 5 2.5.1982 Geostatistical–2016 Analysis22 25 <0.50 <1.00 5 × 5 This analysis enables us to assess the extent to which the transformations made by 2.5. Geostatistical Analysis human activity in the scope of action selected for the territory have influenced the current vulnerabilityThis analysis to flooding. enables us The to assess spatial the relationships extent to which were the parameterized transformations and made assessed by humanthrough activity the use in of the global scope Moran’s of actionI [40 selected] and Anselin for the Local territory Moran’s haveI influenced[41] bivariate the statistics; current vulnerabilityboth are geoprocessing to flooding. tools The from spatial GvSIG relationships desktop 2.5.1 we (AGS,re parameterized Madrid, Spain) and and asse ArcGISssed through10.8 (Esri, the Sacramento, use of global CA, Moran USA).’s I [40] and Anselin Local Moran’s I [41] bivariate statis- tics; bothBivariate are geoprocessing global spatial tools autocorrelation from GvSIG allows desktop us to2.5.1 assess (AGS, the Spain) statistical and correlation ArcGIS 10.8 of (Esri,a set ofUSA). geolocated data obtained spatially and whether the autocorrelation is positive or negative.Bivariate The global bivariate spatial global autocorrelation Moran’s I statistic allows formula us to asses is givens the as statisticalI (9): correlation of a set of geolocated data obtained spatially and whether the autocorrelation is positive n n or negative. The bivariate global Morann ’∑s iI= statistic1 ∑j=1 w formulai,jzizj is given as 퐼 (9): = I n 2 (5) S0 푛 ∑ 푛 z 푛 ∑푖=1 ∑i=푗=11 푤i 푖,푗푧푖푧푗 퐼 = 푛 2 (5) 푆0 ∑푖=1 푧푖 Appl. Sci. 2021, 11, 809 10 of 21