Evaluation of Green and Grey Flood Mitigation Measures in Rural Watersheds

: Floods cause considerable damages worldwide and mitigation of their adverse e ﬀ ects through e ﬀ ective protection measures is needed. Along with the commonly applied “grey” infrastructure, “green” measures that can o ﬀ er additional beneﬁts, such as ecosystem services, are increasingly being considered lately. While the recent research tendencies are focused on the e ﬀ ectiveness and the value of green measures in urban areas, this paper presents a comprehensive ﬁnancial evaluation of green and grey ﬂood mitigation scenarios for a smaller rural watershed. A micro-scale damage model that builds on the hydrodynamic modeling of hazard, detailed asset identiﬁcation, and damage assessment is presented and applied for evaluation of beneﬁts from various ﬂood mitigation measures in the Tamnava watershed in Serbia. Four scenarios are considered: (1) existing ﬂood protection system; (2) green scenario involving new detention basins; (3) grey infrastructure enhancement by rising of the existing levees and diverting ﬂood discharges; and (4) green-grey scenario that combines scenarios (2) and (3). The beneﬁts (loss reduction) are the greatest with the green scenario and marginally higher with the combined green-grey scenario. The results suggest that for small rural watersheds, a holistic, integrative approach that includes both types of infrastructure can provide the most e ﬀ ective ﬂood risk mitigation.


Introduction
State and municipal governments are facing depleting natural resources, adverse climate change impacts, and many socio-economic challenges. Continuation of industrial, transportation, and utility services is critically dependent on continued construction and maintenance of roads, drainage infrastructure, culverts, etc., commonly known as "grey infrastructure".
Despite being essential for economic growth, these infrastructure investments are significant and continue throughout the life of any such project, through annual maintenance costs [1]. This point is of particular importance in infrastructure designed to provide protection from natural hazards, and specifically flood hazards. Large drainage systems and structural flood protection are being exposed to natural decay and require constant investment. They are characterized by a limited life span and declining benefits and are under additional stress by the changing climate [2,3]. Moreover, the human and the material assets protected by grey infrastructure are becoming more valuable over time, so the exposure to flood hazard grows even faster. Grey infrastructure also has a very limited role in ecosystem functions and can lead to degradation of local and regional systems [3]. 2 of 25 In the process of looking for an alternative solution, attention is turning towards nature-based solutions ("green infrastructure") that could be used for mitigation of flood risks [4]. Green infrastructure (GI) is defined as "strategically planned network of natural and semi-natural areas that . . . deliver a wide range of ecosystem services" [5]. Examples of applicable GI for flood risk mitigation are reclaimed wetlands, riverine floodplains, natural detentions, etc. Because of the growing challenges related to urban areas, such as population growth, rapid urbanization, and climate change effects, there has been abundant emergence of projects and research on coping with flood risk in urban areas using GI [6][7][8][9][10][11].
It has been shown in many cases that nature-based GI solutions are less expensive than grey infrastructure (e.g., [12]) and provide a wide array of co-benefits for local economies, the social fabric, and the broader environment [13]. This should be of particular interest to decision-makers, as GI can provide comparable benefits to grey infrastructure at reduced costs in the long term.
On the other hand, there has been very little research on the quantitative and cost effectiveness of GI for flood mitigation in settings other than urban, such as flat valleys of large rivers or small hilly watersheds prone to flash floods. One of the reasons for this is that improved understanding of trade-offs between flood risk reduction and economic consequences in rural settings is needed [2] because the assets in rural settings considerably differ from the ones in urban areas. Recently, the term "large-scale nature-based solutions" is commonly used for the flood risk mitigation measures in rural settings on either watershed or regional scale [4]. In small to medium size watersheds, flood risk reduction measures can vary from conventional grey infrastructure to various scales and levels of GI. In some cases, depending on the scale of the watershed, solutions gravitate towards larger scale nature-based solutions, or some combination of the green and grey infrastructure approach.
Both green and grey infrastructure as flood mitigation measures have their specific advantages and drawbacks and either may not be a universal solution to be applied indiscriminately to all types of watersheds and local conditions. Therefore, the two types of flood mitigation measures should not be a replacement for each other but should rather be considered as complementary measures [12,14,15]. For rural watersheds with some flood mitigation measures already present, the combination of green and grey infrastructure may be the most suitable solution.
Quantitative evaluation of different flood mitigation measures is necessary to select the most effective strategy for flood risk management or for adaptation to climate change [16]. The flood risk assessment, which is generally aimed at evaluating the potential consequences of floods [17], provides a basis for evaluating different mitigation options and making optimal decisions on flood mitigation measures. The flood risk assessment can also serve other various purposes such as identifying vulnerability of communities to floods, developing flood risk maps, or financial appraisals for the insurance sector [18].
The goal of this research is to identify a method of evaluating general flood mitigation strategy for smaller, predominately rural watersheds. The paper presents an evaluation of different flood mitigation scenarios (grey, green, and combined green and grey) for a small rural, agricultural watershed in Serbia through a comprehensive flood risk assessment of each option and their comparison. The evaluation is made with a spatial damage model based on combination of flood hazard and vulnerability assessments. The model can be characterized as the micro-scale flood damage model because it uses high-resolution hazard map outputs from the hydrodynamic modelling and fine-detail, multi-category assets data. The damage assessment is validated against the official post-disaster financial loss appraisals following a major flood event in the region in 2014, and therefore produces estimates of tangible direct and indirect losses with relatively low uncertainties. The flood mitigation options are compared in a traditional way by calculating benefits as avoided damages and losses.

Description of the Study Area
The study area comprises upper watershed of river Tamnava and its tributaries Gračica and Ub, in central Serbia, approximately 60 km southwest of Belgrade, the capital. This area is a part of the larger Kolubara watershed, depicted in Figure 1. The Tamnava watershed covers 726 km 2 and is primarily rural, with 79.3% of the area being cultivated land [19] (Figure 2). Urbanized and industrial land constitutes only 1.2% of the area and is concentrated in two small population centers, the towns of Ub and Koceljeva. The terrain elevation ranges from 470 m a.s.l. in upper reaches of the watershed to 64.4 m a.s.l. at its mouth ( Figure 2), with arable land located in wide floodplains of the three rivers. The agricultural landscape consists primarily of small farms, with emphasis on several main crops. The rest of the economy is based on service industry located in towns of Ub and Koceljeva.

Description of the Study Area
The study area comprises upper watershed of river Tamnava and its tributaries Gračica and Ub, in central Serbia, approximately 60 km southwest of Belgrade, the capital. This area is a part of the larger Kolubara watershed, depicted in Figure 1. The Tamnava watershed covers 726 km 2 and is primarily rural, with 79.3% of the area being cultivated land [19] (Figure 2). Urbanized and industrial land constitutes only 1.2% of the area and is concentrated in two small population centers, the towns of Ub and Koceljeva. The terrain elevation ranges from 470 m a.s.l. in upper reaches of the watershed to 64.4 m a.s.l. at its mouth ( Figure 2), with arable land located in wide floodplains of the three rivers. The agricultural landscape consists primarily of small farms, with emphasis on several main crops. The rest of the economy is based on service industry located in towns of Ub and Koceljeva.   Hydroclimatic regime of the Tamnava watershed is characterized by pronounced seasonal pattern, with the highest flows in early spring due to combined rainfall and snowmelt, and in early summer due to high-intensity convective rainfall events. The latter mainly occur in June, which is the month with the greatest precipitation amount. Mean annual precipitation amounts to 787.7 mm, while mean annual temperature is 11.4 °C for 1981-2010. Mean annual flow of the Tamnava River at Koceljeva is 0.99 m 3 /s, and of the Ub River at the Ub gauge is 1.01 m 3 /s. The highest flows at these gauges are 178 m 3 /s and 146 m 3 /s, respectively, and were observed during a great flood in May 2014.
The study area has been affected by a number of historic flood events (1999,2006,2009,2020, and earlier), including the catastrophic flood in May 2014. Flooding in the Tamnava watershed has often been caused by flooding of the river Kolubara, whose high-water levels would prevent efficient drainage from the Tamnava watershed and cause excessive backflow flooding in its interior. Historically, the floods along the Tamnava River have caused considerable damages, particularly in the vicinity of its confluence with the Kolubara River.
The May 2014 flood was an unprecedented event in the hydrological record of the whole West Balkans region [20]. It was triggered by heavy rain falling during several days over near-saturated soils after several similar antecedent events. The extent of flooding exceeded historical records, with the Kolubara watershed being among the most affected ones. High water levels remained in the river valleys for weeks; the damages were substantial, and the casualties were reported.
A special study on the 2014 catastrophic flood [21], herein referred to as the UNDP study, was intended for reconstruction of this flood event and a comprehensive evaluation of various proposed flood mitigation measures in the Kolubara watershed. This extensive study comprises the following: • A detailed description of the entire Kolubara watershed, including topography, stream network, soil types, geology, and other relevant features of the watershed, such as population, infrastructure, and economy (particularly agriculture and coal mining). Hydroclimatic regime of the Tamnava watershed is characterized by pronounced seasonal pattern, with the highest flows in early spring due to combined rainfall and snowmelt, and in early summer due to high-intensity convective rainfall events. The latter mainly occur in June, which is the month with the greatest precipitation amount. Mean annual precipitation amounts to 787.7 mm, while mean annual temperature is 11.4 • C for 1981-2010. Mean annual flow of the Tamnava River at Koceljeva is 0.99 m 3 /s, and of the Ub River at the Ub gauge is 1.01 m 3 /s. The highest flows at these gauges are 178 m 3 /s and 146 m 3 /s, respectively, and were observed during a great flood in May 2014.
The study area has been affected by a number of historic flood events (1999,2006,2009,2020, and earlier), including the catastrophic flood in May 2014. Flooding in the Tamnava watershed has often been caused by flooding of the river Kolubara, whose high-water levels would prevent efficient drainage from the Tamnava watershed and cause excessive backflow flooding in its interior. Historically, the floods along the Tamnava River have caused considerable damages, particularly in the vicinity of its confluence with the Kolubara River.
The May 2014 flood was an unprecedented event in the hydrological record of the whole West Balkans region [20]. It was triggered by heavy rain falling during several days over near-saturated soils after several similar antecedent events. The extent of flooding exceeded historical records, with the Kolubara watershed being among the most affected ones. High water levels remained in the river valleys for weeks; the damages were substantial, and the casualties were reported.
A special study on the 2014 catastrophic flood [21], herein referred to as the UNDP study, was intended for reconstruction of this flood event and a comprehensive evaluation of various proposed flood mitigation measures in the Kolubara watershed. This extensive study comprises the following: • A detailed description of the entire Kolubara watershed, including topography, stream network, soil types, geology, and other relevant features of the watershed, such as population, infrastructure, and economy (particularly agriculture and coal mining).

•
Thorough analyses of various aspects of flooding in the Kolubara watershed, including descriptions of the major observed floods in the study area. This research uses the part of the UNDP study related to the Tamnava watershed as the starting point. The hydrodynamic model developed in the UNDP study is also used in this research with necessary adjustments as described in Section 2.3. For the sake of consistency with the UNDP study, we adopted its notation for stream sections (sectors) along the Tamnava, Gračica, and Ub rivers. The sectors are distinct in their prevalent land use in the river valleys and proposed options for reducing flood risk. The short description of the sectors is presented in Table 1, and the sector boundaries corresponding to the areas of maximum potential flooding extent are shown in Figure 3. • Thorough analyses of various aspects of flooding in the Kolubara watershed, including descriptions of the major observed floods in the study area. This research uses the part of the UNDP study related to the Tamnava watershed as the starting point. The hydrodynamic model developed in the UNDP study is also used in this research with necessary adjustments as described in Section 2.3. For the sake of consistency with the UNDP study, we adopted its notation for stream sections (sectors) along the Tamnava, Gračica, and Ub rivers. The sectors are distinct in their prevalent land use in the river valleys and proposed options for reducing flood risk. The short description of the sectors is presented in Table 1, and the sector boundaries corresponding to the areas of maximum potential flooding extent are shown in Figure 3.

Flood Mitigation in the Tamnava Watershed
To reduce flood hazard and consequent damages, the flood protection system in the Tamnava watershed has been set up since 1950s. Flood protection system in 2014, which is the existing system used as reference in this study, is primarily based on the levees along several sections of the Tamnava Appl. Sci. 2020, 10, 6913 7 of 25 and Ub rivers ( Figure 3). The levees that protect the populated areas along the Tamnava River are designed for 50-year or 100-year floods, while the remaining levees provide protection of agricultural land from 25-year floods.
The mitigation measures that are proposed in the UNDP study [21] and analyzed in this study include: • Three detention basins with earthen dams at the upper reaches of Tamnava, Gračica, and Ub that should serve as flow regulating and erosion prevention structures during floods and would otherwise retain a minimal amount of water. Considering that detention basins can provide various ecosystem services [22,23], they are considered as green measures in this study. Additionally, natural anti-erosion measures are planned in the UNDP study to complement the detention basins.

•
Reinforcement and crown elevation of the existing levees along the Tamnava and Ub, in sectors 7 , 9, and 12, respectively ( Figure 3). Additionally, construction of a bypass canal connecting Gračica and Ub rivers (sector 13 ) is proposed. All of the above measures are referred to as grey flood protection measures in this study. Locations of the proposed mitigation measures are depicted in Figure 3. The existing levees that already provide certain level of flood protection are planned to be elevated as follows: Sector 7 , bank protection and levee protecting downtown Koceljeva: increasing flood protection level from return period of 50 years to 100 years; Sector 9, agricultural levee on both banks of river Tamnava: increasing flood protection level from return period of 50 years to 100 years; Sector 12, agricultural levee on both banks of River Ub: increasing flood protection level from return period of 25 years to 100 years.
Technical details of the levee system in the watershed and the proposed detention basins are given in the Supplementary material.
To enable comparison of benefits from implementing green and grey measures, and, consequently, identification of optimal flood protection setup for small, predominately rural watersheds, four alternative scenarios of flood protection in the Tamnava watershed are considered in this study:

1.
Baseline scenario (existing conditions, "no action" scenario). This scenario includes only existing flood protection measures and assumes that flood protection level in the watershed would remain the same in the future.

2.
Green scenario includes the three proposed detention basins along with the existing grey infrastructure (i.e., levees) with the flood protection level that they currently provide.

3.
Grey scenario implies raising the existing levees and constructing the bypass canal, as proposed in the UNDP study.

4.
Green-grey scenario includes implementation of both green (detention basins) and grey measures (heightened levees and the bypass canal) in the watershed.
It should be noted that the distinction between green and grey measures is not a clear, straightforward one. Concerning the building material, some elements of green infrastructure (e.g., the dams of the detention basins) may be made of (reinforced) concrete or steel, while the levees (grey measures) are mainly made of local earthen material. In this paper, the distinction between the two groups of measures is not made upon the building material, but rather according to the ecosystem services that the measures can provide. This categorization is in line with the definition provided by European Commission, stating that "green infrastructure provides great benefits for both citizens and biodiversity" [23]. Detention basins can provide habitat for wildlife and keep a fish stock, and their application is deemed to "contribute to meeting the objectives of the 2020 Biodiversity strategy" [23]. In addition, detention basins can offer some recreational opportunities. For the above reasons, the detention basins are considered a green flood mitigation measures in this study, as opposed to the conventional levees and a bypass canal.

Outline of the Methodology
The methodology for evaluating flood mitigation scenarios presented in this paper is based on assessment of flood risk for each scenario. The flood risk assessment follows a process well elaborated in the literature, with risk being quantified as a product of hazard and vulnerability [24,25], and shown in Figure 4. The flood hazard is here represented by maps of flood depth for different hazard scenarios (i.e., flood return periods or probabilities of exceedance) and different mitigation scenarios. The flood hazard maps are obtained as a result of hydrodynamic modelling of flood flows along the Tamnava, Ub, and Gračica rivers, as described in Section 2.3. Therefore, combining the natural hazard and mitigation scenarios produces multiple scenarios that should assist in identifying the optimal solution for the watershed.
European Commission, stating that "green infrastructure provides great benefits for both citizens and biodiversity" [23]. Detention basins can provide habitat for wildlife and keep a fish stock, and their application is deemed to "contribute to meeting the objectives of the 2020 Biodiversity strategy" [23]. In addition, detention basins can offer some recreational opportunities. For the above reasons, the detention basins are considered a green flood mitigation measures in this study, as opposed to the conventional levees and a bypass canal.

Outline of the Methodology
The methodology for evaluating flood mitigation scenarios presented in this paper is based on assessment of flood risk for each scenario. The flood risk assessment follows a process well elaborated in the literature, with risk being quantified as a product of hazard and vulnerability [24,25], and shown in Figure 4. The flood hazard is here represented by maps of flood depth for different hazard scenarios (i.e., flood return periods or probabilities of exceedance) and different mitigation scenarios. The flood hazard maps are obtained as a result of hydrodynamic modelling of flood flows along the Tamnava, Ub, and Gračica rivers, as described in Section 2.3. Therefore, combining the natural hazard and mitigation scenarios produces multiple scenarios that should assist in identifying the optimal solution for the watershed. Assessment of the vulnerability to floods is based on identifying the assets in the focus area that are at risk of flooding. The vulnerability assessment is the basis for flood damage assessment. Vulnerabilities from natural hazards are usually grouped as physical, social, economic, and environmental [26]. However, estimating the potential damages from floods in terms of these four categories can be a complex process associated with high uncertainties. In engineering applications, it is much more common to divide the damages into direct and indirect losses [17,18]. In this study, four categories of vulnerabilities, which were identified as the most important ones and for which sufficient data was available, are considered. These include three categories of direct losses and one category of indirect losses: • Physical damage to the buildings (residential/public).

•
Physical damage to building contents and equipment.

•
Damage to agricultural crops.

•
Population emergency displacement and relocation.
There is an extensive literature on quantifying vulnerability for the considered categories [25,27]. The direct losses are estimated in this study by applying depth-damage relationships, i.e., by calculating the damage in function of the flood depth for each scenario, which is a common approach [16,28] particularly useful for physical assets. The indirect losses of population displacement are Assessment of the vulnerability to floods is based on identifying the assets in the focus area that are at risk of flooding. The vulnerability assessment is the basis for flood damage assessment. Vulnerabilities from natural hazards are usually grouped as physical, social, economic, and environmental [26]. However, estimating the potential damages from floods in terms of these four categories can be a complex process associated with high uncertainties. In engineering applications, it is much more common to divide the damages into direct and indirect losses [17,18]. In this study, four categories of vulnerabilities, which were identified as the most important ones and for which sufficient data was available, are considered. These include three categories of direct losses and one category of indirect losses: • Physical damage to the buildings (residential/public).

•
Physical damage to building contents and equipment.

•
Damage to agricultural crops.

•
Population emergency displacement and relocation.
There is an extensive literature on quantifying vulnerability for the considered categories [25,27]. The direct losses are estimated in this study by applying depth-damage relationships, i.e., by calculating the damage in function of the flood depth for each scenario, which is a common approach [16,28] particularly useful for physical assets. The indirect losses of population displacement are estimated on the basis of daily unit costs per person displaced. Details on the vulnerability and damage assessment are given in Section 2.4. Assessment of risk is performed in function of asset vulnerability and damages for a series of synthetic (design) flood events under each mitigation scenario. Additionally, the same methodology for damage assessment is also applied to the flood event of May 2014 and compared to official damage estimates for the selected losses categories to validate the methodology. The effectiveness of a particular mitigation approach under the specified hazard scenario is measured as the damage avoided by implementing that particular option, i.e., as the difference between the losses of the baseline scenario and the considered mitigation scenarios outlined in Section 2.1.2. This step is elaborated in Section 2.5.

Hazard Assessment
The hydrodynamic model used for flood hazard assessment in this paper was originally developed for reconstruction of the extreme flood event in May 2014 in the UNDP study [21], and it will be referred to as MODEL 2014. It is a 1D model specifically developed to simulate all levee breaches and overtopping, and backwater effects that occurred during flood event in May 2014. This model version is calibrated within the UNDP study according to the observed water stages and surveyed flooding extents. An overview of the characteristics of this model and its calibration are presented in the Supplementary material. Details on the reconstruction of this flood event can be found in the literature [29].
The MODEL 2014 was further modified to enable hydrodynamic simulations and flood hazard assessments under different mitigation scenarios described in Section 2.1. For the existing flood protection, levee breaches and overtopping options were excluded from the MODEL 2014. Further modifications of the hydrodynamic model were aimed at simulations under the future mitigation scenarios by including the three proposed detention basins (green scenario), elevated levees, and the bypass canal (grey scenario) and combination of both types of measures (green-grey scenario).
The hydrodynamic models of the Tamnava, Ub, and Gračica rivers were fed with the hydrographs as the upstream boundary conditions, and stages of the Kolubara River at the confluence with Tamnava as the downstream boundary condition [21]. Only MODEL 2014 was used for non-stationary simulations. The hydrographs used to run MODEL 2014 were either observed at the stream gauges or simulated by hydrologic models [29]. The remaining model versions were used for stationary simulations and were driven with design flood flows of various return periods. The approach of using only the peak flows instead of entire hydrographs was adopted as a conservative approach that maximizes the flooding extents and the flood damage.
The design flows for hydrodynamic simulations were estimated in the UNDP study [21] at the locations of the stream gauges in the watershed by applying the flood frequency analyses of the observed annual maximum floods. These design floods were regionalized to estimate corresponding flows at ungauged locations in the watershed that were needed for hydrodynamic simulations.
All hydrodynamic model versions were developed with the HEC-RAS software [30]. The models included all related structures, such as levees, reservoirs, detention ponds, and bridges. Terrain data were obtained by combining several Digital Terrain Models (DTMs) available for the Kolubara watershed. The combined DTM was further refined according to the structural lines, such as the levees and road network, river cross-sections placed at 500 m distance from additional geodetic surveys, as well as precise data on the bridge structures. The data exchange between HEC-RAS and GIS software was established via the HEC-GeoRAS tool [31].
Simulations with the MODEL 2014 were performed in this paper to enable comparison of losses during the flood event in May 2014 to the official post-disaster estimates of flood damages in the Tamnava watershed, and, thus, to enable validation of the damage assessment approach used in this paper. The other model versions were used for flood hazard and damage assessments under the four analyzed flood mitigation scenarios. Flood hazards were assessed for 2-, 10-, 20-, 50-, 100-, 200-and 1000-year design floods. In all simulations performed, it was assumed that the protection system fulfills its design role completely (no levee breaches).

Asset Identification and Vulnerability Assessment
Three principal types of assets, and four corresponding types of vulnerabilities, were identified for the study: • Building stock and its inventory (contents). • Agricultural production.

•
Population affected by flooding and cost associated with its temporary displacement.
Selection of the vulnerability categories was based on availability of the data, prevalence of assets in the area, and compatibility with categories used in the UNDP study [21] for the purpose of comparison of flood damage estimates.

Building Stock and Its Inventory
In this paper, the micro-level flood damage assessment is conducted, meaning that evaluation of vulnerabilities is performed on a level of individual buildings for a series of hazard events. This approach is more detailed in comparison to the aggregated approach applied in the UNDP study [21]. A total of 5984 permanent structures were identified from aerial photogrammetric imagery within the river sectors' boundaries (the maximum potential extent of flooding). The number of structures and their occupational classes were verified in the field. Of these, 63.5% are residential buildings (primarily single family, low-density and farm houses), 27.8% are agricultural buildings, 6.4% are commercial, educational, and government structures, with the remaining 2.3% being light industrial buildings. The characteristics of the above types of structures (i.e., occupancy classes) are typical for this region and this part of the country [32].
The damages to structures caused by floods are estimated using the depth-damage function (DDF). The damage as a function of water depth in the structure is typically expressed not in absolute monetary values, but as percentage of the "maximum damage value" [33], which is the value of replacing the completely destroyed structure with a new one.
In this paper, the DDFs for residential, commercial, industrial, and transportation occupancy classes of buildings were adopted from the Joint Research Center's global database [33], with modifications applied for Serbia and for inflation. It was conservatively assumed that buildings have no basements. All the modifications applied to DDFs from [33] are made in accordance with field investigations.
The shape of DDF for buildings was also adopted for the corresponding building content. The maximum damage value for the building content ranged between 50% and 150% of the maximum damage value for the corresponding building structure, as shown in Table 2. An example of the DDF for residential buildings used in the study is presented in Figure 5.

Agricultural Production
Agricultural production and its vulnerability are analyzed through three principal crops productions (corn, potatoes, apples) and unbaled hay [21]. The seasonality of the flood event and resulting various levels of vulnerability of each crop were not taken into account.
The crops were conservatively assumed to have a constant buyout cost and that the production yield would remain unchanged over time. The resulting monetary revenue represents the maximum potential agricultural production losses for any given season. The values for buyout costs and production yields were adopted from the UNDP study and modified for inflation, as presented in Table 3. Table 3. Agricultural production and maximum potential losses due to flooding, expressed in terms of annual revenue.

Culture
Yield The vulnerability of the agricultural resources and crop producing areas was defined as their direct exposure to flooding, regardless of the amount of time spent under water. Damages to crops such as apples and corn are highly dependable on the depth of flooding [34], and have been adjusted accordingly to conform to field observations and high water marks. Potatoes and unbaled hay located within the flood boundaries were considered a total loss, regardless of the corresponding depth of flooding.

Population Affected by Flooding and Temporarily Displaced
Population may be affected by flooding in various ways, from having their domiciles directly damaged to being evacuated and displaced due to impeding environmental, structural, infrastructural, or health-related hazards. In this study, only the inhabitants whose residences were directly affected by flooding are considered. The population potentially exposed to flooding within the study area is estimated at 9155 [21], with 3800 residential structures in the same domain and an average number of 2.41 inhabitants per structure.

Agricultural Production
Agricultural production and its vulnerability are analyzed through three principal crops productions (corn, potatoes, apples) and unbaled hay [21]. The seasonality of the flood event and resulting various levels of vulnerability of each crop were not taken into account.
The crops were conservatively assumed to have a constant buyout cost and that the production yield would remain unchanged over time. The resulting monetary revenue represents the maximum potential agricultural production losses for any given season. The values for buyout costs and production yields were adopted from the UNDP study and modified for inflation, as presented in Table 3. Table 3. Agricultural production and maximum potential losses due to flooding, expressed in terms of annual revenue. The vulnerability of the agricultural resources and crop producing areas was defined as their direct exposure to flooding, regardless of the amount of time spent under water. Damages to crops such as apples and corn are highly dependable on the depth of flooding [34], and have been adjusted accordingly to conform to field observations and high water marks. Potatoes and unbaled hay located within the flood boundaries were considered a total loss, regardless of the corresponding depth of flooding.

Population Affected by Flooding and Temporarily Displaced
Population may be affected by flooding in various ways, from having their domiciles directly damaged to being evacuated and displaced due to impeding environmental, structural, infrastructural, or health-related hazards. In this study, only the inhabitants whose residences were directly affected by flooding are considered. The population potentially exposed to flooding within the study area is estimated at 9155 [21], with 3800 residential structures in the same domain and an average number of 2.41 inhabitants per structure.
The vulnerability of population can be defined through displacement time, i.e., time for residents to be relocated until imminent danger is gone or until the necessary repairs are completed. The displacement duration varies with the depth of interior flooding [35] and is presented in Table 4. The cost of relocation usually includes both the one-time evacuation costs and the expenses related to temporary housing. This study considered only the temporary housing expenses estimated at approximately 10 € per person per day [36].

Valuation of the Flood Mitigation Scenarios
The effectiveness of a particular flood mitigation scenario S can be measured by the reduction of flood risk in comparison to the baseline scenario. In other words, the benefits gained by implementing specific set of flood mitigation measures are quantified in terms of avoided damages, i.e., as the difference between flood damage D 0 estimated without these measures (baseline scenario) and the damage D S with these measures (flood mitigation scenario S).
The expected annual damage EAD is the common indicator of the flood risk and can be obtained from the probability distribution of the damages [37]. In practice, it is calculated as the sum-product of the damages and their annual probabilities of exceedance: where M is the number of probabilities of flood hazard p i for which the damages D i are evaluated. The expected annual benefit from the mitigation scenario S can then be computed based on all M considered return periods T i = 1/p i , and taking into account L loss categories: where EAB S denotes expected annual benefit for flood mitigation scenario S, D 0,ij are damages estimated for the baseline scenario (existing flood protection level in the watershed), and D S,ij are damages under flood mitigation scenario S. Subscript j indicates loss category, whereas subscript i is related to flood return period. Four different loss categories are considered in this study (L = 4). For the long-term monetary analysis, it is useful to compute the total benefit as the net present value NPV S from implementing specific flood mitigation scenario S by taking into account discount rate d and planning horizon N, as follows: The long-term benefits are calculated for two planning horizons, namely 30 and 50 years, assuming the discount rate of 7%, as in the UNDP study [21].

Results
The results of the conducted flood loss assessments are complex and are examined from viewpoints of reduction of flood hazard and reduction of flood losses, i.e., the benefits from implementing various flood mitigation measures. To evaluate flood loss assessment in this research, the flood losses in flood event in May 2014 are compared to corresponding officially made estimates.

Flood Hazard
The flood hazard maps under the baseline and three alternative flood mitigation scenarios in the Tamnava watershed result from hydrodynamic modelling, as described in Section 2.3, for each considered return period (2, 10, 20, 50, 100, 200, and 1000 years). Table 5 presents peak flows at the locations of the three stream gauges under all the four flood mitigation scenarios, and for all considered return periods. The hazard map of the 100-year flood under existing flood protection system is shown in Figure 6. As described in Section 2.1, the existing levees already provide some flood protection level to certain parts of the Tamnava watershed, primarily to those in urban setting (e.g., sectors 7 , 11 and 11 '). These parts of the watershed have lower flood hazard level under current conditions than the other sectors in the watershed. Flood hazard is particularly pronounced in the most downstream parts of the watershed (e.g., sectors 9 and 12, some parts of sectors 8 and 13; see Figure 6).   Figure 7. Although the existing levees provide protection from the 25-and 50-year floods (Section 2.1), the detention basins proposed within the green scenario provide sufficient reduction of the 100-year flood, so that it can be conveyed without overtopping in these sectors. These results are obtained also for the 200-year flood, while the 1000-year flood (not shown here) causes levee overtopping in sectors 9 and 12. Reduced flood hazard is also noticed in sector 13′, where construction of the bypass canal is planned (see Figure 3). Comparison of the flood hazard maps for different return periods shows that the grey measures have effects in their immediate proximity, while the green measures (i.e., detention basins) reduce peak flows and pertinent flood levels along the downstream river sections, thus having farther-reaching effects.   Figure 7. Although the existing levees provide protection from the 25-and 50-year floods (Section 2.1), the detention basins proposed within the green scenario provide sufficient reduction of the 100-year flood, so that it can be conveyed without overtopping in these sectors. These results are obtained also for the 200-year flood, while the 1000-year flood (not shown here) causes levee overtopping in sectors 9 and 12. Reduced flood hazard is also noticed in sector 13 , where construction of the bypass canal is planned (see Figure 3). Comparison of the flood hazard maps for different return periods shows that the grey measures have effects in their immediate proximity, while the green measures (i.e., detention basins) reduce peak flows and pertinent flood levels along the downstream river sections, thus having farther-reaching effects.   Figure 7. Although the existing levees provide protection from the 25-and 50-year floods (Section 2.1), the detention basins proposed within the green scenario provide sufficient reduction of the 100-year flood, so that it can be conveyed without overtopping in these sectors. These results are obtained also for the 200-year flood, while the 1000-year flood (not shown here) causes levee overtopping in sectors 9 and 12. Reduced flood hazard is also noticed in sector 13′, where construction of the bypass canal is planned (see Figure 3). Comparison of the flood hazard maps for different return periods shows that the grey measures have effects in their immediate proximity, while the green measures (i.e., detention basins) reduce peak flows and pertinent flood levels along the downstream river sections, thus having farther-reaching effects. The flood hazard maps enable identification of the affected assets and the depth of flooding water. Figure 8 shows water depth in each building in downtown Ub due to 200-year flood event under baseline and green scenarios and illustrates how the fine-scale asset data and a fine-resolution hazard map facilitate the micro-scale assessment of damages.

Damage Assessment for the 2014 Flood Event
The losses from the May 2014 flood are estimated using the damage assessment approach described in Section 2.4 and compared to the corresponding official post-event damage assessments given in the UNDP study [21]. The comparison was possible for two asset categories: combined damages for building structures and content, and agricultural losses (population displacement was not considered in [21]). Difference between the estimated damages to the residential structures in two assessments is only 2.5% (the estimate of €37,897,822 in this study largely corresponds to the official post-disaster estimate of €38,869,942). The two estimates of damages to crops agree fully (€13,113,716 in this research vs. official estimate of €13,110,157).
These comparisons suggest that flood losses estimated in this paper can be considered valid and can provide reasonably realistic assessments of benefits gained by implementing various flood mitigation measures.

Benefits from Different Flood Mitigation Scenarios
Losses due to floods of different return periods and under different flood mitigation scenarios are computed for each river sector in the Tamnava watershed by summing the losses over the four considered loss categories (building structures and their content, agriculture, and population displacement). The losses for all sectors are presented in Figure 9 in function of the flood return period. The losses under the grey scenario are lesser than for the baseline scenario up to the 100-year return period, which is the design criterion for the elevated levees in the grey scenario. The losses for the green scenario are smaller than the losses for the baseline and grey scenarios for all return periods. Finally, the losses for the green-grey scenario show just a marginal improvement in comparison to the green scenario. Similar graphs of damage distribution are plotted for each sector (not shown in the paper), and although they show somewhat different patterns depending on the relative improvement of protection for a particular sector, they generally show that the green scenario is superior to the grey scenario, while the green-grey scenario is just marginally better than the green scenario.

Damage Assessment for the 2014 Flood Event
The losses from the May 2014 flood are estimated using the damage assessment approach described in Section 2.4 and compared to the corresponding official post-event damage assessments given in the UNDP study [21]. The comparison was possible for two asset categories: combined damages for building structures and content, and agricultural losses (population displacement was not considered in [21]). Difference between the estimated damages to the residential structures in two assessments is only 2.5% (the estimate of €37,897,822 in this study largely corresponds to the official post-disaster estimate of €38,869,942). The two estimates of damages to crops agree fully (€13,113,716 in this research vs. official estimate of €13,110,157).
These comparisons suggest that flood losses estimated in this paper can be considered valid and can provide reasonably realistic assessments of benefits gained by implementing various flood mitigation measures.

Benefits from Different Flood Mitigation Scenarios
Losses due to floods of different return periods and under different flood mitigation scenarios are computed for each river sector in the Tamnava watershed by summing the losses over the four considered loss categories (building structures and their content, agriculture, and population displacement). The losses for all sectors are presented in Figure 9 in function of the flood return period. The losses under the grey scenario are lesser than for the baseline scenario up to the 100-year return period, which is the design criterion for the elevated levees in the grey scenario. The losses for the green scenario are smaller than the losses for the baseline and grey scenarios for all return periods. Finally, the losses for the green-grey scenario show just a marginal improvement in comparison to the green scenario. Similar graphs of damage distribution are plotted for each sector (not shown in the paper), and although they show somewhat different patterns depending on the relative improvement of protection for a particular sector, they generally show that the green scenario is superior to the grey scenario, while the green-grey scenario is just marginally better than the green scenario. The expected annual damages (EAD) for the four scenarios are presented in Table 6, together with their breakdown per loss categories considered. Share of the loss categories in EAD varies across the scenarios. For example, the largest share in EAD have damages to agriculture under all flood mitigation scenarios, except for the existing conditions. Combined losses due to damages to building structures and content are slightly higher than the losses in agricultural production under the baseline scenario. These results can be explained by the fact that Tamnava is predominately rural, agricultural watershed, with low percentage of urbanized and industrial areas. Losses due to displacement of the affected population represent the smallest part of losses due to floods, with up to 5% of total EADs, and remain relatively constant under all scenarios.  Figure 10 shows EAD for individual river sectors under all flood mitigation scenarios. The greatest losses under the baseline scenario are in sectors 8,9,11, and 13′, corroborating high flood hazard in these sectors (see Figure 6). In addition, these sectors are characterized by large areas (e.g., sector 8) or by presence of urbanized or industrial zones (e.g., sector 13′, see Table 1). Sectors 8, 11, The expected annual damages (EAD) for the four scenarios are presented in Table 6, together with their breakdown per loss categories considered. Share of the loss categories in EAD varies across the scenarios. For example, the largest share in EAD have damages to agriculture under all flood mitigation scenarios, except for the existing conditions. Combined losses due to damages to building structures and content are slightly higher than the losses in agricultural production under the baseline scenario. These results can be explained by the fact that Tamnava is predominately rural, agricultural watershed, with low percentage of urbanized and industrial areas. Losses due to displacement of the affected population represent the smallest part of losses due to floods, with up to 5% of total EADs, and remain relatively constant under all scenarios.  Figure 10 shows EAD for individual river sectors under all flood mitigation scenarios. The greatest losses under the baseline scenario are in sectors 8,9,11, and 13 , corroborating high flood hazard in these sectors (see Figure 6). In addition, these sectors are characterized by large areas (e.g., sector 8) or by presence of urbanized or industrial zones (e.g., sector 13 , see Table 1). Sectors 8, 11, and 13 are not covered by grey measures, so EAD for the grey scenario is the same as for the existing conditions. Figure 10 also shows that measures of the green and green-grey scenarios result in lower damages than under grey measures, especially in the most downstream sectors 9 and 12. Under the grey scenario, the heightened levees in these sectors provide protection from the 100-year flood, so the EAD is the result of levee overtopping by more extreme events (200-and 1000-year floods). On the other hand, the detention basins in the green and green-grey scenarios efficiently reduce damages due to such extreme flood events in these sectors, and consequently result in lower EAD than under the grey scenario.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 18 of 25 and 13 are not covered by grey measures, so EAD for the grey scenario is the same as for the existing conditions. Figure 10 also shows that measures of the green and green-grey scenarios result in lower damages than under grey measures, especially in the most downstream sectors 9 and 12. Under the grey scenario, the heightened levees in these sectors provide protection from the 100-year flood, so the EAD is the result of levee overtopping by more extreme events (200-and 1000-year floods). On the other hand, the detention basins in the green and green-grey scenarios efficiently reduce damages due to such extreme flood events in these sectors, and consequently result in lower EAD than under the grey scenario. As can be seen in Figure 9 and Figure 10, the greatest benefits, i.e., reduction of losses relative to the baseline scenario, are obtained under the green and green-grey scenarios. Figure 11 shows the benefits for different return periods and indicates that the grey scenario has lower benefits than the remaining two scenarios for floods of 100-year return period and smaller. For more extreme floods, there are no benefits under the grey scenario, what is in accordance with designed heightening of the levees to the 100-year protection level. As can be seen in Figures 9 and 10, the greatest benefits, i.e., reduction of losses relative to the baseline scenario, are obtained under the green and green-grey scenarios. Figure 11 shows the benefits for different return periods and indicates that the grey scenario has lower benefits than the remaining two scenarios for floods of 100-year return period and smaller. For more extreme floods, there are no benefits under the grey scenario, what is in accordance with designed heightening of the levees to the 100-year protection level.  Table 7 presents the expected annual benefits (EAB, Section 2.5) for the considered scenarios. This table also suggests that the greatest benefits are obtained under the green and green-grey scenarios, which include detention basins that generally have greater effects on reduction of flood hazard and losses than the levees. Hence, the benefits from implementing measures within the grey scenario are considerably smaller, and reflect primarily in reduction of losses to building structures and contents. This can be explained by the fact that the levees in the grey scenario are intended for protection of populated areas in which the buildings represent the most valuable assets. Since the levees make local impact on flood hazard and loss mitigation, reduction in the remaining, more remote, agricultural parts of the watershed are substantially smaller compared to the other two scenarios. The costs of population temporary displacement are considerably lower than the losses in other categories considered, and this pattern is also exhibited in the benefits. Difference in the benefits between the green and green-grey scenarios is minor. This is because the detention basins reduce flood peaks efficiently in the downstream sections so that there is no overtopping of the existing systems of levees. In other words, the levee heightening included in the green-grey scenario has minor effect in combination with implementation of the detention basins.   Table 7 presents the expected annual benefits (EAB, Section 2.5) for the considered scenarios. This table also suggests that the greatest benefits are obtained under the green and green-grey scenarios, which include detention basins that generally have greater effects on reduction of flood hazard and losses than the levees. Hence, the benefits from implementing measures within the grey scenario are considerably smaller, and reflect primarily in reduction of losses to building structures and contents. This can be explained by the fact that the levees in the grey scenario are intended for protection of populated areas in which the buildings represent the most valuable assets. Since the levees make local impact on flood hazard and loss mitigation, reduction in the remaining, more remote, agricultural parts of the watershed are substantially smaller compared to the other two scenarios. The costs of population temporary displacement are considerably lower than the losses in other categories considered, and this pattern is also exhibited in the benefits. Difference in the benefits between the green and green-grey scenarios is minor. This is because the detention basins reduce flood peaks efficiently in the downstream sections so that there is no overtopping of the existing systems of levees. In other words, the levee heightening included in the green-grey scenario has minor effect in combination with implementation of the detention basins. The benefits achieved by implementing green, grey, and green-grey measures, computed relative to the baseline scenario are also presented for each sector in Table 8 in terms of benefits over the 50-year planning horizon (NPV), and illustrated in Figure 12 in terms of EAB. Four river sectors (7,8,11,13) are not protected by the levees, nor would be under the grey scenario, so there are no benefits under this scenario in these sections. There are also no benefits for sector 11", which is already protected up to the 100-year flood and would not be additionally protected in the grey scenario. The benefits under the green and green-grey scenarios are the same in all sections, except for section 13 , where the bypass canal is proposed (see Figure 3). As already discussed, levee heightening is not so beneficial if detention basins are implemented in the watershed. The benefits achieved by implementing green, grey, and green-grey measures, computed relative to the baseline scenario are also presented for each sector in Table 8 in terms of benefits over the 50-year planning horizon (NPV), and illustrated in Figure 12 in terms of EAB. Four river sectors (7,8,11,13) are not protected by the levees, nor would be under the grey scenario, so there are no benefits under this scenario in these sections. There are also no benefits for sector 11", which is already protected up to the 100-year flood and would not be additionally protected in the grey scenario. The benefits under the green and green-grey scenarios are the same in all sections, except for section 13′, where the bypass canal is proposed (see Figure 3). As already discussed, levee heightening is not so beneficial if detention basins are implemented in the watershed.

Discussion
The presented methodology for finding an optimal flood mitigation setup for small rural watersheds is based on a comprehensive economic evaluation of flood losses under different flood mitigation scenarios. This was a challenging task because the results and experiences from other studies, especially on flood damage estimation, cannot be readily transferred and applied to other

Discussion
The presented methodology for finding an optimal flood mitigation setup for small rural watersheds is based on a comprehensive economic evaluation of flood losses under different flood mitigation scenarios. This was a challenging task because the results and experiences from other studies, especially on flood damage estimation, cannot be readily transferred and applied to other regions [38]. Most of the existing research on flood damages is focused on urban areas, and less so on rural agricultural areas. For example, methodologies for estimating damages in buildings are much more developed than for other assets (e.g., [39]). Considering that floods cause the greatest losses to crops globally compared to other natural hazards [40], a need emerges for extensive research on flood damage and risk assessment in rural, agricultural watersheds.
Flood hazard assessment under current conditions in the Tamnava watershed shows that the existing levee system already provides some level of flood protection, primarily in the populated areas where the highest-value assets are located. The proposed mitigation measures under green and green-grey scenarios provide the greatest reduction in flood hazard and, consequently, in flood losses in the agricultural areas that presently have no flood protection, such as in sectors 7, 8, and 13. Although the share of agricultural damages in total losses depends on local conditions in the watershed [34] and prevalence of higher-value crops [41], the results presented here clearly indicate that the agricultural damages are not marginal and have a critical impact on the selection of flood protection development scenarios in predominately agricultural watershed.
Estimation of flood losses in this research is accompanied by inherent uncertainties in both flood hazard and damage assessments, as also demonstrated in many previous studies [18,38,42]. Specifically, flood hazard assessment under grey and green-grey scenario is estimated assuming that grey systems (levees and the bypass canal) fulfill their design role, and that there is no flooding of the interior areas. To accurately estimate the flood hazard, a more robust, potentially 2D hydrodynamic model capable of simulating backwater or some interior area flooding could be utilized. Although previous research shows that enhanced hydrodynamic modelling of flood hazard does not substantially reduce uncertainties in damage assessments [43], there is also evidence that flood damages occur outside of the delineated floodplain boundaries corresponding to certain flood return periods [44]. This is important for agricultural watersheds considering that the damages of this type are generally related to large areas and could be affected significantly by the flooding extent.
Uncertainties in the damage assessment can be generated by multiple causes, including vulnerability categories taken into account, spatial scale of the damage model, depth-damage relationships, and estimated asset values (e.g., [38]). In terms of the vulnerability assessment, this research can be characterized as the "lower boundary" analysis because it includes a limited number of key categories. For example, taking other crop cultures into account could increase agricultural losses [41]. Some of the additional categories that could be considered include transportation and energy infrastructure damages and loss of functionality, traffic delays, loss of critical facilities (water and wastewater systems, hospitals, communications, emergency service), lost productivity, business displacement, etc. [17]. More vulnerability categories could be taken into consideration if the adequate data were available, resulting in a more comprehensive flood loss estimates. Alternatively, using rough estimates or generalized (aggregated) data would lead to uncertain damage assessment, as shown by [45]. This is the reason for applying the "lower boundary" assessment approach in this study. If a particular scenario can be proven to be effective with the limited number of vulnerability categories, then it would be even more advantageous if the benefits from additional categories were included. Furthermore, including categories such as ecosystem services offered by green measures would also contribute to higher effectiveness of green mitigation scenarios. Despite the limited choice of vulnerability categories, the proposed flood mitigation scenarios for the Tamnava watershed are compared in relative terms, allowing their unbiased ranking.
In this paper, a comprehensive micro-scale flood damage assessment is applied, with locally adjusted depth-damage relationships and locally estimated asset values. Previous research has shown that the micro-scale approach yields realistic damage estimates [46], and has also recommended that the asset values are adjusted to the regional economic situation and property characteristics [38].
Validation of the damage models is not a usual practice in their development [42,47]. In this paper, the damage assessment approach in this paper is validated successfully for the case of the flood event in May 2014, leading to a conclusion that the approach yields plausible damage estimates. Therefore, it is expected that the uncertainties in flood hazard and limited number of vulnerability categories do not reflect the overall ranking of the scenarios.
In this study, the effects of the detention basins on peak flow reduction provide a significant advantage for the green and green-grey scenarios. This is expected because the detention basins with considerable storage volume have a principal role in reducing peak flows. The same effects in terms of flood risk reduction would be possible with a grey detention basin of similar storage volume. However, the advantage of a green detention basin should be corroborated by its ecosystem services, the value of which should be included in the assessment of benefits of a flood mitigation scenario [13]. In our case study, the green-grey scenario provides slightly greater loss reduction than the green one. This difference may not suffice to provide its advantage over the pure green scenario if the construction cost of additional grey measures is high. Therefore, for identification of an optimal flood protection development scenario, a comprehensive cost-benefit analysis needs to be performed taking into account environmental, social, and other benefits beside the economic ones [13].
The comparison of the four flood mitigation scenarios in this paper suggests that maintaining the existing flood protection system and complementing it with larger scale green measures intended for reducing flood flows could be an optimal flood protection setup for predominately agricultural watersheds. Further research is needed to perform a sensitivity analysis of the selected flood mitigation scenario to all input data, including terrain information, hydrodynamic models, and various loss categories. The findings of this research are to be tested in different watersheds and for various flood mitigation scenarios that would include a wide variety of measures with special emphasis on green, nature-based solutions. Finally, flood protection based on such solutions is important for wider and sustainable development of the rural areas [48].

Conclusions
This paper presents a comprehensive, micro-scale flood damage assessment under various scenarios of flood mitigation in a predominately rural, agricultural watershed. The objective of this research is identification of optimal flood protection setup that would result in maximum benefits, i.e., reduction of losses due to floods in such watersheds. Considering a relatively small body of research on this topic up to date, especially on valuation of damages to agricultural assets [41], and considering impacts that flood mitigation in rural watersheds can have on downstream urban areas [34], this paper presents an effort to fill this gap and provide a contribution to a better understanding of flood risk assessment and management in general.
It is demonstrated that the detention basins, as an example of green flood mitigation measures that reduce peak flood flows in downstream river sections, can significantly decrease flood hazard in rural, agricultural watersheds. In the case of the Tamnava watershed, this peak flow reduction is sufficient to enable the existing system of protection levees to operate without failure and accept design floods without overtopping. On the other hand, heightening of levees, as an example of grey flood protection measures, increases flood protection locally, but with the potential of transferring some of the flood hazard further downstream. Unlike urban watersheds, in which high-priced assets are localized in relatively small areas, rural watersheds are characterized by larger areas and long river banks that would require long and excessively expensive levee systems. Additionally, since the crops are the most valuable assets in rural watersheds, flood losses to agriculture represent a great portion of total losses in such areas, as demonstrated by a fine-scale detailed damage assessment in this paper. Therefore, with the assets distributed across large areas in rural watersheds with significantly lower asset values per unit area than in urbanized watersheds, grey flood protection measures such as the levees cannot be as effective in reducing flood risk as the green measures reducing peak flows in rural watershed, as shown in this research.
This paper confirmed that hybrid-type flood mitigation setups, i.e., combining large-scale, far-reaching green mitigation solutions that reduce peak flood flows with the existing smaller scale grey flood protection measures, can be the most effective solution for rural watersheds with prevalent agricultural land use and production.