Assessment of Fire Effects on Surface Runoff Erosion Susceptibility: The Case of the Summer 2021 Forest Fires in Greece

: The wildﬁres of summer 2021 in Greece were among the most severe forest ﬁre events that have occurred in the country over the past decade. The conﬂagration period lasted for 20 days (i.e., from 27 July to 16 August 2021) and resulted in the devastation of an area of more than 3600 Km 2 . Forest ﬁre events of similar severity also struck other Mediterranean countries during this period. Apart from their direct impacts, forest ﬁres also render an area more susceptible to runoff erosion by massively removing its vegetation, among other factors. It is clear that immediately after a forest ﬁre, most areas are much more susceptible to erosion. In this paper, we evaluate the erosion hazard of Attica, Northern Euboea, and the Peloponnese that were devastated by forest ﬁres during the summer of 2021 in Greece, in comparison with their geological and geomorphological structures, as well as land cover and management. Given that a very signiﬁcant part of these areas were burnt during the major conﬂagrations of this summer, erosion risk, as well as ﬂood risk, are expected to be very high, especially for the coming autumn and winter. For the evaluation of erosion risk, the burnt areas were mapped, and the ﬁnal erosion-risk maps were constructed through GIS software. The ﬁnal maps suggest that most of the burnt areas are highly susceptible to future surface runoff erosion events.


Introduction
The wildfires of summer 2021 in Greece were among the most severe forest fire events that have occurred in the country over the past decade. The conflagration period lasted for 20 days (i.e., from 27 July to 16 August 2021) and resulted in the devastation of an area of more than 3600 Km 2 [1].
Erosion refers to the removal of soil and rock particles (sediments) and their transportation away from their initial source [2][3][4]. In most cases, particles are mainly removed from their initial position via water.
Surface runoff erosion-risk maps have been developed for certain areas. The most common mapping methods include either mathematical models [5][6][7][8], or in situ observations and bibliographical research [6][7][8][9][10]. In the first method, erosion risk is evaluated and determined based on certain criteria posed by researchers. Therefore, in many cases, the results are subjective. In the second method, the main drawback is the fact that it is usually non-functional when it comes to areas at a regional scale; many errors might arise which would not be of high significance at a small scale, and the necessary field measurements might be very challenging [6] or very rough [8]. Field studies have been improved during the last decades, due to the development of remote sensing techniques [11][12][13]. Among the most important advances in field measurements and observations are the LiDAR [14][15][16][17][18] and unmanned aerial vehicles (UAVs or drones) [19,20]. against erosion than a bush-covered area. The tree type also influences soil erosion [48]. Wide-leaved trees (such as walnut trees, fig trees and sycamores) infiltrate a greater amount of water than needle-leaved ones (e.g., pines) or thin-leaved ones (e.g., kermes).
In addition, land cover refers to both land use and land management. Land cover primarily regards vegetation characteristics. For example, the denser the vegetation cover, the higher the protection it offers against soil erosion. This means that forest areas are less vulnerable to erosion than shrubs. Land management refers to the applied land techniques. An area covered by forests is much more protected than an agricultural area, which is, in turn, less susceptible to erosion than an urban area. Even in the case of agricultural areas, there are many subcases. A cultivated area is more easily eroded than a non-cultivated agricultural area [49]. Furthermore, different agricultural practices bear a different impact on the soil's vulnerability to runoff erosion [50]. For example, tilling techniques generally lead to an increase in the sand's proportion in the soil, as well as a decrease in the soil's permeability [51,52]. In addition, overgrazing reduces an area's protection against soil erosion. On the contrary, the creation of artificial terraces is very helpful in reducing the negative effects of flowing water. Increasing the soil's roughness can also be protective. Finally, another human intervention that severely affects an area's susceptibility to runoff erosion are extensive fires, although they may also occasionally be owed to natural causes. It has been proven that almost any region is far more vulnerable to erosion immediately after a fire. A conflagration results in the reduction of both vegetation and biomass protective cover. This means that all types of protection against erosion are diminished. Therefore, all areas are far more susceptible to erosion, as well as floods, after a conflagration.
Generally speaking, runoff erosion is increased immediately after a fire [53][54][55][56][57][58][59]. However, certain researchers [60][61][62] have studied erosion risk in certain areas after fires, and did not spot any significant increase in erosion risk in comparison with the erosion risk of periods where no fires had occurred. In most cases, in areas affected by forest fires, runoff erosion, as well as flooding problems arise. They usually result in many infrastructure damages, as well as fatalities, not to mention ecological and environmental issues [63]. Additionally, erosion results in the denudation of cultivable areas, thus creating agricultural problems, as well as problems regarding both the quality and quantity of potable water [64,65]. Soil is the component that enables plants to grow and, therefore, allows agricultural activities to take place. Consequently, its removal hinders said activities. Moreover, through the removal of soil, several nutrients (either naturally present in soils or artificially through fertilizers, pesticides etc.) are removed as well. They follow the course of water, thus usually ending up in water bodies (e.g., lakes), or in the sea. One of the most typical consequences of this regime is eutrophication.
Fire affects an area's susceptibility to erosion by impacting the soil properties, among other factors [66,67]. For example, physical and other properties can be altered due to the heat caused by the fire, including the soil's content in organic material, its porosity and permeability, and its shear strength, as well as its infiltration capacity [62,[68][69][70][71][72][73][74][75]. Another direct effect of fires is the removal of vegetation [76], which acts as a protective cover that prevents the soil particles from being removed under normal conditions [75,77,78]. Fire also reduces soil roughness [41]. Of course, the intensity of the fire itself plays a very important role, as a more intense fire and/or a fire that has incinerated a large area, increases erosion, as well as flood risk, to a much greater effect than a milder one [72,79].
In this paper, we examine the vulnerability to surface runoff erosion of several Greek regions that were recently affected by forest fires, and try to identify areas that are more likely to suffer from future erosion events. The areas studied in this research do not share the same geological and geomorphological characteristics, meaning that the effects of forest fires will be examined in connection with the overall structure. The selected areas have suffered by fire events during the summer of 2021, namely Attica, Euboea, and the Peloponnese. The goal of this research is to assess their vulnerability to erosion immediately after the devastating conflagrations, through the development of corresponding maps. For this purpose, the burnt areas were mapped through GIS, and they were combined to other data of the study areas, such as geology, geomorphology, land cover, etc. The final erosionrisk maps illustrate which parts of the burnt areas are more prone to imminent extreme weather phenomena regarding surface runoff erosion.

Study Areas
Three Greek areas were studied: Attica, Northern Euboea, and the Peloponnese (Figure 1). Attica is one of the regional units (formerly prefectures) in the geographical department of Central Greece, covering an area of 3023 Km 2 . Euboea is Greece's second largest island (after Crete) and one of Central Greece's regional units (formerly prefectures), covering an area of 3658 Km 2 . The Peloponnese is one of Greece's geographical departments, covering an area of 21650 Km 2 . All areas are characterized by Mediterranean climate, i.e., by hot and dry summers and mild winters. Rainfalls mainly occur during the wet season, during autumn and winter. According to the climatic data provided by the Hellenic National Meteorological Service [80], in Attica, the mean annual precipitation reaches 402 mm, as measured from the stations of Elliniko, Elefsina, Tatoi and Nea Philadelphia. The minimum precipitation is during the summer months (8 mm in June, 7 mm in July and 6 mm in August), whereas maximum precipitation occurs during November (61 mm), December (69 mm) and January (54 mm). In Euboea, mean annual precipitation reaches 506 mm, according to measurements from adjacent meteorological stations (Skyros, Tanagra, Aliartos and Lamia). Maximum and minimum precipitation height occurs during the same months as in Attica (maximum precipitation reaches 66 mm in November, 82 in December and 70 in January, whereas minimum precipitation reaches 13 mm in June, 10 in July and 13 in August). In the Peloponnese, measurements from the stations of Tripolis, Kalamata and Velo indicate that the mean annual precipitation reaches 665 mm, with minimum values during the summer (8 mm in June, 6 in July and 11 in August) and maximum values during November (108 mm), December (115 mm) and January (94 mm). Rapid and intense weather phenomena, such as thunderstorms, and strong winds are not uncommon in the study areas.

Attica
The areas burnt during the devastating fires of August 2021 include three individual areas: the broader regions of Varympompi, Villia and Sounion national forest (Lavrion The increased vulnerability of these areas to forest fires during this period is totally justified. Initially, the largest part of these areas was, as will be shown, covered by softleaved vegetation (forests, cultivations, etc.) that renders them highly vulnerable to forest fires. Moreover, several severe heat waves have previously struck the Balkan region. Temperatures reached and even exceeded 40 • C for several successive days. As a result, the air temperature was very high during the summer of 2021, whereas soil moisture was very low. These conditions, in connection to the lack of precipitation during this period, rendered the Balkan Peninsula, and Greece in particular, extremely vulnerable to forest fires [1].

Attica
The areas burnt during the devastating fires of August 2021 include three individual areas: the broader regions of Varympompi, Villia and Sounion national forest (Lavrion area) ( Figure 2). In Varympompi area, the burnt land is mainly composed of lacustrine, fluvial and terrestrial Neogene sediments found in its northwestern part, and, to a lesser extent, limestones (locally dolomites), covering its southeastern part. A small part is covered by metamorphic rocks (clastic and, secondarily, marbles), mainly located in the far northeastern part [1]. With regard to land cover, the largest part of the burnt land was covered by forest vegetation (mainly broad-leaved, coniferous and mixed forests) and transitional zones between forest and shrub vegetation, whereas a smaller part was covered by cultivable and/or cultivated areas and urban fabric [77]. The burnt area was 85 Km 2 . Settlements within the burnt areas include Varympompi, Ippokrateios Politeia, Tatoi, Afidnes, Pefkofyton, Aghios Stephanos and Kryonerion. The urban fabric in the area, especially in the settlements of Ippokrateios Politeia, Tatoi and Afidnes, is well-developed. According to the aforementioned, land cover did facilitate the spreading of the fo fires, in association with the heat waves that had previously affected Greece. Regar the fire-affected areas of Attica Prefecture, they host well-developed drainage syst partly draining urban areas. There are two main hydrographic systems in the wider of Varympompi area ( Figure 2). The main drainage networks are the Oinois river, flow towards Marathon, and Kifissos river, flowing towards the western suburbs of Attic the fire-affected area in Villia, a drainage system crosses it for about 10 Km in an E West direction, and finally ends up in the area of Elefsina. There are also smaller ca In Sounion National Forest (Lavrion area), the burnt area was mainly covered by schists and, to a lesser extent, marbles. As its name states, its largest part was covered by forests of several types (e.g., sclerophyllous, coniferous) and shrubs, although there were also cultivated areas. More specifically, five main vegetation units were found, each covering an almost equal area to the others. These units included cultivations, coniferous forests, sclerophyllous forests, mixed forests and shrubs. The total burnt area in Lavrion was 5.5 Km 2 [1]. As far as Villia area is concerned, almost all of the burnt lands were covered by natural vegetation, i.e., forests, mainly sclerophyllous and coniferous, as well as shrubs and grasslands. The total burnt area in Villia reached 66 Km 2 .
According to the aforementioned, land cover did facilitate the spreading of the forest fires, in association with the heat waves that had previously affected Greece. Regarding the fire-affected areas of Attica Prefecture, they host well-developed drainage systems, partly draining urban areas. There are two main hydrographic systems in the wider area of Varympompi area ( Figure 2). The main drainage networks are the Oinois river, flowing towards Marathon, and Kifissos river, flowing towards the western suburbs of Attica. In the fire-affected area in Villia, a drainage system crosses it for about 10 Km in an East-West direction, and finally ends up in the area of Elefsina. There are also smaller catchments passing through the fire-affected area, flowing towards Nea Peramos (Figure 3).  In Attica, the areas burnt during the devastating fires of August 2021 inc broader regions of Varympompi, Villia and Sounion national forest, with the tot area reaching approximately 155 Km 2 [1].
The morphological slope distribution of Attica is shown in Figure 4. Larg cover more than half of the fire-affected areas. More specifically, 53% of the area have an inclination of more than 10°, 1% of which exceeds the inclination of 30 tionally, 29% of the total area is of medium slope, with inclinations fluctuating 10° and 4°, whereas 18% is characterized by an inclination of less than 4°. In Attica, the areas burnt during the devastating fires of August 2021 include the broader regions of Varympompi, Villia and Sounion national forest, with the total burnt area reaching approximately 155 Km 2 [1].
The morphological slope distribution of Attica is shown in Figure 4. Large slopes cover more than half of the fire-affected areas. More specifically, 53% of the area's slopes have an inclination of more than 10 • , 1% of which exceeds the inclination of 30 • . Additionally, 29% of the total area is of medium slope, with inclinations fluctuating between 10 • and 4 • , whereas 18% is characterized by an inclination of less than 4 • .

Northern Euboea
In Euboea island, a significant area in its northern part was burnt, specifically the area of Pefki-Aghia Anna-Rovies, i.e., an area of more than 505 Km 2 [1] ( Figure 5). As far as geology is concerned, the largest part of the burnt area is composed of lacustrine, terrestrial and fluvio-lacustrine sediments of Neogene age, with ultramafic rocks (peridotites and amphibolites) covering small, individual parts. Carbonates and mica schists cover only a very small part. According to CORINE [77], almost all the studied area was covered by coniferous and broad-leaved forests, whereas only a small part included cultivable land, meaning that any fire of a small scale could easily obtain the severity of the August 2021 fires. Among the most significant settlements within the conflagrated area are Aghia Anna, Rovies, Lymni, Pappades, Kerasia, Vassilika and Pefki. The rural fabric is not welldeveloped, as most of the area is woodland and all the aforementioned settlements are small villages.

Northern Euboea
In Euboea island, a significant area in its northern part was burnt, specifically the area of Pefki-Aghia Anna-Rovies, i.e., an area of more than 505 Km 2 [1] ( Figure 5). As far as geology is concerned, the largest part of the burnt area is composed of lacustrine, terrestrial and fluvio-lacustrine sediments of Neogene age, with ultramafic rocks (peridotites and amphibolites) covering small, individual parts. Carbonates and mica schists cover only a very small part. According to CORINE [77], almost all the studied area was covered by coniferous and broad-leaved forests, whereas only a small part included cultivable land, meaning that any fire of a small scale could easily obtain the severity of the August 2021 fires. Among the most significant settlements within the conflagrated area are Aghia Anna, Rovies, Lymni, Pappades, Kerasia, Vassilika and Pefki. The rural fabric is not welldeveloped, as most of the area is woodland and all the aforementioned settlements are small villages.
The fire-affected area is characterized by abrupt slopes, downcutting erosion and a dense drainage system. Most of the drainage basins are small with small hydrographic networks, yet there are few large drainage basins covering a significant part of the fireaffected area, ending up at Gouves, Rovies, Neochori and Krya Vrysi.
The morphological slope distribution of Euboea is shown in Figure 6. Large slopes cover more than half of the fire-affected areas. Moreover, 63.1% of the slopes have an inclination between 10 • and 30 • (medium slopes), and 14.5% of the total area is characterized by low slope, with inclinations fluctuating between 10 • and 4 • , whereas 5.8% is characterized by an inclination of less than 4 • (very low slopes). The remaining 16.6% of the area's slopes have an inclination of more than 30 • (high slopes). Land 2022, 11, x FOR PEER REVIEW 9 of 20 The fire-affected area is characterized by abrupt slopes, downcutting erosion and a dense drainage system. Most of the drainage basins are small with small hydrographic networks, yet there are few large drainage basins covering a significant part of the fireaffected area, ending up at Gouves, Rovies, Neochori and Krya Vrysi.
The morphological slope distribution of Euboea is shown in Figure 6. Large slopes cover more than half of the fire-affected areas. Moreover, 63.1% of the slopes have an inclination between 10° and 30° (medium slopes), and 14.5% of the total area is characterized by low slope, with inclinations fluctuating between 10° and 4°, whereas 5.8% is characterized by an inclination of less than 4° (very low slopes). The remaining 16.6% of the area's slopes have an inclination of more than 30° (high slopes).

The Peloponnese
In the Peloponnese, four main individual fire outbreaks took place, namely in Aigialeia, Ancient Olympia, Diavolitsi and East Mani. The burnt area reached approximately 291 Km 2 (Figure 7) [1].
In the Aigialeia region, post-alpine sediments are exclusively found, mainly marls of the Pliocene-Pleistocene age, sandy clays and Quaternary conglomerates [1]. Land cover in the burnt area includes coniferous and shrub vegetation and cultivations, their proportion being almost equal [77]. Artificial surfaces, primarily railways and motorways, mainly occupy the coastal part of the region. A total of 3 Km 2 were burnt, including the settlements of Lampirion, Ziria, Ano Ziria and Kamari [78].  Land cover in these four areas favored the spreading of the devastating fire significant part of them was covered by low vegetation (e.g., cultivations), even t forests covered a smaller part in comparison with the previous regions (Euboea a tica). The fire-affected areas are characterized by high elevations (Figure 7) and stee phological slopes (Figure 8). In the areas of Ancient Olympia and East Mani, a welloped drainage system is observed, contrary to the area of Aigialeia, which is main In all cases but Aigialeia, more than 80% of the burnt areas are characterized by whose inclination exceeds 10°. More specifically, in Aigialeia, 22% of the slopes h inclination of more than 30°, and 32% are between 10° and 30°. The remaining 20 27% include slopes with an inclination of 4° to 10°, and less than 4°, respectively. cient Olympia, the slopes with an inclination of 10° to 30° reach 52%, whereas tho ceeding 30° reach 28%. Only 5% is characterized by an inclination of less than 4°. Sim in the Diavolitsi area, 62% and a 24% of the slopes have an inclination of 10° to 30 more than 30°, respectively; the flat areas (less than 4° in inclination) only cover 3% total area. Finally, in Mani, the areas with more than 30°, and between 10° to 30° s are almost equal, covering 48% and 35% of the burnt areas, respectively. Ultimatel of the area is characterized by slopes between 4° and 10°. In the Aigialeia region, post-alpine sediments are exclusively found, mainly marls of the Pliocene-Pleistocene age, sandy clays and Quaternary conglomerates [1]. Land cover in the burnt area includes coniferous and shrub vegetation and cultivations, their proportion being almost equal [77]. Artificial surfaces, primarily railways and motorways, mainly occupy the coastal part of the region. A total of 3 Km 2 were burnt, including the settlements of Lampirion, Ziria, Ano Ziria and Kamari [78].
In Ancient Olympia, post-alpine sediments are also exclusively found, specifically of Pliocene to Quaternary age, and of several phases [81]. Land cover mainly includes cultivations, and only forest to a small extent. A total area of 135 Km 2 was devastated during the wildfires, the affected settlements including Lampeti, Lalas and Ancient Olympia [1].
In Diavolitsi area, almost 47 Km 2 were burnt. The region was half occupied by cultivations and half by forests. In Messinia (Diabolitsi) area, the burnt land is mainly covered by Radiolarites and pelagic limestones with intercalations of silex, marls and sandstones. A small part is covered by flysch and limestones. Finally, in Eastern Mani, almost 105 Km 2 were burnt, including the settlements Kastania, Melissa, Sidirokastron, Myrsini and Gytheion. This area too was occupied by forests and cultivations, with the former covering a larger part than the latter. In all cases, the settlements are small, mainly rural villages with no significant urban fabric.
Land cover in these four areas favored the spreading of the devastating fires, as a significant part of them was covered by low vegetation (e.g., cultivations), even though forests covered a smaller part in comparison with the previous regions (Euboea and Attica). The fire-affected areas are characterized by high elevations (Figure 7) and steep morphological slopes (Figure 8). In the areas of Ancient Olympia and East Mani, a welldeveloped drainage system is observed, contrary to the area of Aigialeia, which is mainly flat. In all cases but Aigialeia, more than 80% of the burnt areas are characterized by slopes whose inclination exceeds 10 • . More specifically, in Aigialeia, 22% of the slopes have an inclination of more than 30 • , and 32% are between 10 • and 30 • . The remaining 20% and 27% include slopes with an inclination of 4 • to 10 • , and less than 4 • , respectively. In Ancient Olympia, the slopes with an inclination of 10 • to 30 • reach 52%, whereas those exceeding 30 • reach 28%. Only 5% is characterized by an inclination of less than 4 • . Similarly, in the Diavolitsi area, 62% and a 24% of the slopes have an inclination of 10 • to 30 • , and more than 30 • , respectively; the flat areas (less than 4 • in inclination) only cover 3% of the total area. Finally, in Mani, the areas with more than 30 • , and between 10 • to 30 • slopes, are almost equal, covering 48% and 35% of the burnt areas, respectively. Ultimately, 11% of the area is characterized by slopes between 4 • and 10 • .

Materials and Methods
For this research, the overall geological and geomorphological structures of the areas were studied. The data used were derived from the bibliography, from several urements through GIS programs (namely ArcGIS Pro and MapInfo Pro) and from work conducted for the needs of previous research in the study areas. Additionally, in each of the burnt areas was conducted immediately after the fires. Therefore, mor observations were made regarding the morphology of these areas after the fire in co ison with the previous morphology, as well as considering the direct effects of the such as the destruction of vegetation.
All the collected data were imported into GIS. (ArcGIS Pro version 2.8.3 and Ma Pro version 12.5) and the erosion risk was evaluated using a Boolean logic-based m The main parameters used in this stage were land uses, slope, drainage density and ogy. These logical associations are expressed through a set of Boolean logic rules S1). Boolean rules assign four values to lithology and slope parameters (high, me low, very low) and take into account high drainage density and land-use type value final estimation of erosion risk is categorized into six values (very high, high, me low, very low and deposition).

Results and Discussion
The surface runoff erosion hazard for Attica Prefecture is shown in Figure  broader area of Vilia is not characterized by a high risk, despite the medium to high gradient, which would facilitate the removal of soil and rock material (i.e., accelerat sion). Yet, the burnt areas were slightly more vulnerable to erosion than the surrou area. The former was characterized by a low risk, whereas the latter by a very low The same situation was observed in Varympompi. The morphological slopes wer tively high in the broader area, meaning that they already facilitated runoff erosion were generally uniform with less fluctuations. Yet, the burnt area was characterize low risk, and the surrounding area was generally of very low risk. The area of So showed a similar regime. In this case too, the slopes were generally uniform and in diate, compared with the other two regions. The erosion risk was low to medium f burnt area, and very low for the rest of the region.

Materials and Methods
For this research, the overall geological and geomorphological structures of the burnt areas were studied. The data used were derived from the bibliography, from several measurements through GIS programs (namely ArcGIS Pro and MapInfo Pro) and from field work conducted for the needs of previous research in the study areas. Additionally, a visit in each of the burnt areas was conducted immediately after the fires. Therefore, more field observations were made regarding the morphology of these areas after the fire in comparison with the previous morphology, as well as considering the direct effects of the fires, such as the destruction of vegetation.
All the collected data were imported into GIS. (ArcGIS Pro version 2.8.3 and MapInfo Pro version 12.5) and the erosion risk was evaluated using a Boolean logic-based model. The main parameters used in this stage were land uses, slope, drainage density and lithology. These logical associations are expressed through a set of Boolean logic rules (Table S1). Boolean rules assign four values to lithology and slope parameters (high, medium, low, very low) and take into account high drainage density and land-use type values. The final estimation of erosion risk is categorized into six values (very high, high, medium, low, very low and deposition).

Results and Discussion
The surface runoff erosion hazard for Attica Prefecture is shown in Figure 9. The broader area of Vilia is not characterized by a high risk, despite the medium to high slope gradient, which would facilitate the removal of soil and rock material (i.e., accelerates erosion). Yet, the burnt areas were slightly more vulnerable to erosion than the surrounding area. The former was characterized by a low risk, whereas the latter by a very low risk. The same situation was observed in Varympompi. The morphological slopes were relatively high in the broader area, meaning that they already facilitated runoff erosion: they were generally uniform with less fluctuations. Yet, the burnt area was characterized by a low risk, and the surrounding area was generally of very low risk. The area of Sounion showed a similar regime. In this case too, the slopes were generally uniform and intermediate, compared with the other two regions. The erosion risk was low to medium for the burnt area, and very low for the rest of the region. The erosion risk map of Northern Euboea is shown in Figure 10. It is clear that whole area is particularly prone to runoff erosion. This is due to the combination of geological structure and the morphological slope distribution. Neogene sediments p vailed in the study area, which are very erodible. Additionally, the slope gradient w characterized by high to very high values. Yet, the burnt areas were more susceptible one factor. More specifically, the burnt areas were generally characterized by a very hi risk, which was high to medium in some individual parts; the non-burnt areas were high to medium risk, as a total. The erosion risk map of Northern Euboea is shown in Figure 10. It is clear that the whole area is particularly prone to runoff erosion. This is due to the combination of its geological structure and the morphological slope distribution. Neogene sediments prevailed in the study area, which are very erodible. Additionally, the slope gradient was characterized by high to very high values. Yet, the burnt areas were more susceptible by one factor. More specifically, the burnt areas were generally characterized by a very high risk, which was high to medium in some individual parts; the non-burnt areas were of high to medium risk, as a total. Figure 11 shows the erosion risk distribution for the fire-affected areas of the Peloponnese. In the area of Aigialeia, erosion risk is well-distributed; it ranges from very low to very high depending on the area. The part near the sea was of very low to non-existent risk. The rest of the area shows a similar regime to northern Euboea; the burnt area was of medium to very high risk, whereas the surrounding region was of low to medium risk. High risk only characterized a few individual parts. This regime was also facilitated by the slope, which ranges from very high in the inland part to very low in the coastal part. Regarding eastern Mani, the situation was similar; the burnt region was of low, and in certain areas, medium erosion risk, whereas the erosion risk was very low in the surrounding area. The morphological slope had an evident impact on the distribution of the high and low hazard areas, as the slopes generally decrease eastwards, as do the erosion hazards in general. In Diavolitsi, correspondingly, the burnt areas were characterized by a low to medium risk, whilst the surrounding areas were of low to very low risk. Finally, in ancient Olympia, the difference in the erosion risk was more intense, with the burnt areas showing high to very high erosion risk and the rest of the area low to non-existent risk.  Figure 11 shows the erosion risk distribution for the fire-affected areas of the Peloponnese. In the area of Aigialeia, erosion risk is well-distributed; it ranges from very low to very high depending on the area. The part near the sea was of very low to non-existent risk. The rest of the area shows a similar regime to northern Euboea; the burnt area was of medium to very high risk, whereas the surrounding region was of low to medium risk. High risk only characterized a few individual parts. This regime was also facilitated by the slope, which ranges from very high in the inland part to very low in the coastal part. Regarding eastern Mani, the situation was similar; the burnt region was of low, and in certain areas, medium erosion risk, whereas the erosion risk was very low in the surrounding area. The morphological slope had an evident impact on the distribution of the high and low hazard areas, as the slopes generally decrease eastwards, as do the erosion hazards in general. In Diavolitsi, correspondingly, the burnt areas were characterized by a low to medium risk, whilst the surrounding areas were of low to very low risk. Finally, in ancient Olympia, the difference in the erosion risk was more intense, with the burnt areas showing high to very high erosion risk and the rest of the area low to non-existent risk. Besides land degradation, other indirect impacts of soil erosion include increased flood risk, problems in the quality of potable water, disruptions in the ecosystems and others [82]. Runoff erosion poses a significant geological problem of Europe. Yet, the in- Besides land degradation, other indirect impacts of soil erosion include increased flood risk, problems in the quality of potable water, disruptions in the ecosystems and others [82]. Runoff erosion poses a significant geological problem of Europe. Yet, the intensity of runoff erosion is not uniform all over the continent. In the Mediterranean, runoff erosion is very intense, as the dry period lasts for a significant amount of time and heavy rain and storm events are very common during the wet season, whereas in many parts, the slopes are steep and/or covered with erodible material [6,83]. On the contrary, the northern parts of Europe are much less susceptible to runoff erosion [6]. In most cases, runoff erosion is also facilitated by inappropriate cultivation methods, as well as other human interventions, such as deforestation and overgrazing, which are very common in Mediterranean countries [84][85][86][87].
According to CORINE [88], the evaluation of erosion risk in Greece, Italy, southern France, Spain and Portugal, and Southern France, showed that a part of Italy and northern Greece are of low risk, whereas the rest of Greece is characterized by moderate risk. The central part of Italy and almost all the Iberian Peninsula are characterized by moderate to high risk. Another research by RIVM [89], considering almost all Europe, has classified most parts of the Balkan peninsula, Italy, Spain, Turkey, northern Switzerland, Germany and western Ukraine as highly susceptible to surface runoff erosion, whereas northern France and eastern Ukraine and Belarus are of medium risk.
In Greece, according to van der Knijff et al. [90], who applied the USLE model in Europe, the annual soil loss generally fluctuates between 0 and 20 t/ha. The largest part of northern and northwestern Greece and the far northwestern part of the Peloponnese are characterized by a low annual soil loss (generally less than 1 t/ha). An exception would be the Strymonas basin, Chalkidiki and Olympus, where annual soil loss reaches, or possibly exceeds, 10 t/ha [90]. Central Greece, the Peloponnese and the Ionian islands are characterized by the highest soil rate, which is less than 5 t/ha/year in only a few areas. Most of these areas undergo a soil loss of 5 to 10 t/ha/year. Higher rates (10-20 t/ha/year) are mainly found in the northern Peloponnese.
Soil erosion, as well as other natural disasters (e.g., landslides, floods) are in most cases hindered or reduced when vegetation is abundant. Plants intercept the falling water through their leaves. In this way, water during a rainfall has a decreased velocity, i.e., a lower erodibility, whereas plants also absorb soil water through their roots, meaning that the soil becomes less saturated, and more water can infiltrate into the soil; hence, the surface runoff is significantly reduced [91][92][93][94]. Moreover, plants restrain soil particles, thus hindering them from being removed; they also create holes in the soil through their roots, further enabling water infiltration. Furthermore, a strong presence of plants means that soil roughness is increased, both directly (through the plants themselves) and indirectly (through dead leaves, trunks, etc.). In this way, soil particles are prevented from being transported downstream [95][96][97]. It is clear that in most cases, surface runoff erosion, especially when it comes to soils, is highly increased after a conflagration, as a massive amount of vegetation is usually rapidly dismantled.
The areas studied in this paper have undergone a significant vegetation loss due to the devastating fires of summer 2021. Therefore, surface runoff erosion was expected to be higher in the burnt areas than the areas that remained unaffected. In addition, the morphological slope of the study areas clearly affects the distribution of runoff erosion hazard. The results of the geomorphological and erosion-risk mapping have confirmed this expectation. All the burnt parts of the 10 individual areas that were severely affected by the forest fires, which were studied in this research (Vilia, Sounion, Varympompi, Northern Euboea, Aigialeia, Diavolitsi, Ancient Olympia, and Eastern Mani) are characterized by an increased susceptibility to erosion than their surrounding regions, despite the fact that, for each of these areas, the geomorphological and geological structure of the burnt and the non-burnt parts are similar.

Conclusions
Vegetation plays a very important role in reducing surface runoff and hindering erosion, especially as far as soils are concerned. Areas affected by wildfires are more susceptible to surface runoff erosion, particularly when the fire event is very recent, and the dismantled vegetation has not regenerated. The more intense the conflagrations and/or the larger the area affected, the higher the erosion risk. The wildfires of summer 2021 were among the most devastating conflagration events that have struck Greece during recent decades, not only regarding duration and area, but also intensity. Therefore, surface runoff erosion risk is very high in the burnt areas. Our research has shown that all the burnt parts of the studied areas (Vilia, Sounion and Varympompi in Attica, Northern Euboea, Diavolitsi, Ancient Olympia, Aigialeia and Eastern Mani in the Peloponnese) are more susceptible to near-future erosion events than their surrounding non-burnt areas.