Model-Based Evaluation of Urban River Restoration: Conflicts between Sensitive Fish Species and Recreational Users

Urban rivers are socioecological systems, and restored habitats may be attractive to both sensitive species and recreationists. Understanding the potential conflicts between ecological and recreational values is a critical issue for the development of a sustainable river-management plan. Habitat models are very promising tools for the ecological evaluation of river restoration projects that are already concluded, ongoing, or even to be planned. With our paper, we make a first attempt at integrating recreational user pressure into habitat modeling. The objective of this study was to analyze whether human impact is likely to hinder the re-establishment of a target species despite the successful restoration of physical habitat structures in the case of the restoration of the Isar River in Munich (Germany) and the target fish species Chondostroma nasus L. Our analysis combined high-resolution 2D hydrodynamic modeling with mapping of recreational pressure and used an expert-based procedure for modeling habitat suitability. The results are twofold: (1) the restored river contains suitable physical habitats for population conservation but has low suitability for recruitment; (2) densely used areas match highly suitable habitats for C. nasus. In the future, the integrated modeling procedure presented here may allow ecological refuge for sensitive target species to be included in the design of restoration and may help in the development of visitor-management plans to safeguard biodiversity and recreational ecosystem services.


Introduction
Centuries of human activities have directly or indirectly degraded, damaged, transformed, or entirely destroyed aquatic ecosystems [1], threatening worldwide riverine ecosystem services [2], Sustainability 2018, 10, 1747; doi:10.3390/su10061747www.mdpi.com/journal/sustainabilitye.g., water security for fauna, flora, and humans [3].Restoration activities are essential to re-establish the functional capacities of ecosystems for providing and maintaining biological and cultural diversity [1,2,4].River restoration is a term applied to a wide range of activities concerned with "repairing waterways that can no longer perform essential ecological and social functions" [5].Project surveys in the United States, Japan, Australia, Germany, and France have shown that improving instream physical habitats through hydro-morphological changes is one of the most common goals of restoration [6][7][8][9][10][11]. Hydro-morphological river restoration intends to recreate natural forms of the river course or at least mitigates the human impact, for instance removing artificial embankments and meandering the river bed.However, despite significant hydro-morphological changes, indicator species may not recover [12][13][14][15], which suggests that stressors other than hydro-morphological degradations still affect the biota in restored river sections [10,12].Urban rivers are showpieces of a socioecological system [11], and their restoration often faces conflicts between ecological quality and recreational uses.They are particularly impoverished in habitat diversity and quality [16,17], and their dysfunctions have been described as the "urban river syndrome" [18].At the same time, urban rivers have important cultural functions [19], and their restoration (or rehabilitation) is often driven by citizens' demands for more near-natural landscapes and for a greater recreational potential [11,20].Rather than ecological restoration in the sense of a return to pre-disturbance conditions, which is unlikely in urban areas [21], urban river restoration often aims at integrating ecological demands into urban planning [20].According to the socioecological approach of the "River Culture Concept" [2], improvements in ecosystem functions and biodiversity will have a positive effect on human culture, e.g., riverine recreational activities [22].Junker and Buchecker [23] also showed that users respond very positively to morphological river restoration.With population growth and increasing demands for recreational opportunities in urban areas, an increase in the use of river beaches for leisure activities has been predicted since the 1990s [24], especially on urban rivers with improving water quality [2].However, recreational activities cause pressures on ecosystems; for example, they negatively impact wildlife, soil, and vegetation [25][26][27][28][29][30].Therefore, conflicts between ecological restoration and recreational uses are likely to exist [31,32], but impacts of recreational users have been overlooked in comparison to those of other perturbations [33][34][35].
While project assessment is an important procedure for feedback and guidance, few river restoration projects implement an evaluation procedure.Studies showed that in the United States, less than 10% of river restorations have been monitored [6].Similar observations have been made in Bavaria, Germany [36].Recent studies identified five difficulties faced in the evaluation of river-restoration projects: 1.
Requisite data are missing because pre-and post-restoration monitoring is absent for many projects or occurs over a very short period [9,36].

2.
Control sites are inappropriate for assessing the success of restoration [9,37].

3.
Biological indicators have serious limitations [38][39][40], such as a lack of identification of the causes of failures in restoration.4.
Current evaluation procedures investigate ecological responses, but social evaluation remains lacking [9,36,41].Few studies have performed parallel investigations into both social and biological aspects, e.g., water quality and project acceptance, but without the integration of both aspects into a single evaluation method, the resulting conclusions remain speculative.

5.
The time span needed for the re-establishment of a target species is often longer than the monitoring period [9,12,13].
Furthermore, the evaluation method should inform about potential conflicts, integrate social and ecological assessments, and investigate the quality of the physical habitat.
Habitat models are very promising tools and may contribute to solving this part of the problem.Habitat models are appropriate tools to investigate physical habitats and their ecological functions.Habitat simulation models are increasingly used in water management and the investigation of habitat changes caused by hydropeaking [42][43][44], weir removal [45,46], and the presence of reservoirs.The models also help in the decision-making process for choosing the most efficient restoration design [47].The Physical Habitat Simulation (PHABSIM) model was the first widely available physical habitat model and was used in the 1990s to assess restoration measures [48], but it uses univariate functions and ignores interactions among physical habitat variables, which cause major limitations [49].The Computer Aided Simulation Model for Instream Flow Requirements (CASiMiR) solves these limitations by using an approach based on fuzzy logic, sets, and rules to integrate interactions between physical habitat variables and enable the prediction of habitat selection for different fish species and life cycle stages in defined river types [49][50][51][52][53].However, habitat models such as CASiMiR are solely based on physical hydro-morphological variables that affect the ecological function of habitats, whereas human impacts by recreational uses of river banks on the availability of instream habitats for the target fauna has not yet been considered.
Fish habitats are a common indicator of the quality of and richness in aquatic habitats [9,[54][55][56][57].Modeling of physical fish habitats is a suitable tool for assessing the impact (positive or negative) of hydro-morphological changes on physical instream habitats.While impacts of hydro-morphological changes on fish habitats are well described in the literature, little information exists about the response of fish to nonfishing recreational pressure.Some studies showed that anthropogenic noise negatively affects the functions of fish habitats; for example, road traffic (around 85 dB) disturbs fish migration [58], piling sounds (>120 dB) disturb communication [59] and cause total dispersion of fish schools or a drift in their habitat from a depth of 15-20 m [60], and ship noise (around 85 dB) inhibits larval development and increases rates of larval death due to predators [61].However, we found no study that investigated the impact of recreationists along or in rivers on fish habitats.
Through a novel methodological approach, our study made a first attempt at integrating social aspects (human use of riverbanks) into habitat modeling for a sensitive fish species to detect potential conflicts between physical fish habitats and recreational uses.We investigated, in the case of the common nase (Chondrostoma nasus L.) in the restored section of the Isar River in the center of Munich, (i) whether hydro-morphological urban river restoration succeeded in providing physical habitats for different life-cycle stages of the target species, (ii) how urban recreational pressure is distributed with regard to fish habitats and how it modifies the availability of highly suitable habitats, and (iii) whether urban recreational pressure is likely to explain the absence of C. nasus from the studied restored river section.

Study Area
The Isar River drains part of the Karwendel Mountains (Northern Alps), crosses the Munich conurbation (Germany), and joins the Danube River (Figure 1).The pre-alpine river has been used since prehistoric times as a trade route from the Alps, but only using rafts.Major environmental issues began in the 1920s with the construction of 28 hydroelectric power plants.The whole river has been canalized, river water has been diverted several times, and the Sylvenstein Reservoir (1954)(1955)(1956)(1957)(1958)(1959) has been built in the upstream part of the river to mitigate flood risk.
The discharge of the Isar River in Munich (Figure 2) is minimal during the winter (11-202 m 3 /s between 1959 and 2012; www.hnd.bayern.de),although brief and slight increases in water quantity may occur, which are mostly due to rain events (annual mean maximum discharge (HQ) in winter = 450 m 3 /s).In spring, the discharge increases due to snowmelt (>63.8 m 3 /s), although substantial decreases (to as little as 8.63 m 3 /s) may happen sporadically due to late freeze and snowfall.During summer, minor flash flood events (mean maximum discharge (HQ) in summer = 395 m 3 /s) caused by summer storms in the Alps are frequent but they only occasionally reach a very high discharge (650-1050 m 3 /s).Dry periods may also cause minimal discharges during summer, although regulation by the upstream Sylvenstein Reservoir usually avoids these to ensure water supply for hydroelectric power plants and cooling water for the nuclear power plant.At the gauging site in Munich, NQ locally occurred because of important water diversion to supply hydropower plants (Figure 2).Since the 1990s, restoration projects have been carried out to improve the ecological status of the river and the esthetics of its riverscape in order to decrease flood risk and to increase recreational potential.For example, the project "New Life for the Isar" (1999-2011) (Figure 3) was carried out by the Bavarian Water Agency in collaboration with the Munich city government to restore 8 km of the Isar River crossing the city of Munich in such a manner that high ecological quality could be expected in the future [62,63].The project extended from the Großhesseloher Bridge (48°4′29.59′′N, 11°32′25.83′′E) to Museum Island (48°7′41.42′′N, 11°34′46.88′′E).The restoration focused on the eastern side of the river, while the western side remained unchanged because of topographic limits, namely, a very steep wooded slope ending at a stone wall in the river that separates the river from the housing area.One of the goals of the restoration project was the improvement of the habitats for the endemic fish species, such as the emblematic common nase (C.nasus).The project focused on hydro-morphological improvements such as the removal of the concrete embankment, the creation  Since the 1990s, restoration projects have been carried out to improve the ecological status of the river and the esthetics of its riverscape in order to decrease flood risk and to increase recreational potential.For example, the project "New Life for the Isar" (1999-2011) (Figure 3) was carried out by the Bavarian Water Agency in collaboration with the Munich city government to restore 8 km of the Isar River crossing the city of Munich in such a manner that high ecological quality could be expected in the future [62,63].The project extended from the Großhesseloher Bridge (48°4′29.59′′N, 11°32′25.83′′E) to Museum Island (48°7′41.42′′N, 11°34′46.88′′E).The restoration focused on the eastern side of the river, while the western side remained unchanged because of topographic limits, namely, a very steep wooded slope ending at a stone wall in the river that separates the river from the housing area.One of the goals of the restoration project was the improvement of the habitats for Since the 1990s, restoration projects have been carried out to improve the ecological status of the river and the esthetics of its riverscape in order to decrease flood risk and to increase recreational potential.For example, the project "New Life for the Isar" (1999-2011) (Figure 3) was carried out by the Bavarian Water Agency in collaboration with the Munich city government to restore 8 km of the Isar River crossing the city of Munich in such a manner that high ecological quality could be Sustainability 2018, 10, 1747 5 of 27 expected in the future [62,63].The project extended from the Großhesseloher Bridge (48  88 E).The restoration focused on the eastern side of the river, while the western side remained unchanged because of topographic limits, namely, a very steep wooded slope ending at a stone wall in the river that separates the river from the housing area.One of the goals of the restoration project was the improvement of the habitats for the endemic fish species, such as the emblematic common nase (C.nasus).The project focused on hydro-morphological improvements such as the removal of the concrete embankment, the creation of seminatural fishways (Figure 4) at the 400 m long Flaucher site (Figure 4), and the construction of near-natural river-bottom ramps (Figure 4).Special care was also taken to reproduce a near-natural waterscape, and one island (Willow Island) (Figure 4) was created to satisfy citizens' esthetic demands.Long-term monitoring confirmed that the good chemical and biological (macrozoobenthos) conditions required for the re-establishment of C. nasus had been achieved.However, electrofishing showed no recovery of C. nasus despite the potential for rapid recovery of the species [64][65][66].
Sustainability 2018, 10, x FOR PEER REVIEW 5 of 26 of seminatural fishways (Figure 4) at the 400 m long Flaucher site (Figure 4), and the construction of near-natural river-bottom ramps (Figure 4).Special care was also taken to reproduce a near-natural waterscape, and one island (Willow Island) (Figure 4) was created to satisfy citizens' esthetic demands.Long-term monitoring confirmed that the good chemical and biological (macrozoobenthos) conditions required for the re-establishment of C. nasus had been achieved.However, electrofishing showed no recovery of C. nasus despite the potential for rapid recovery of the species [64][65][66].

Biology of Chondrostoma nasus
Chondrostoma nasus (L.), named the common nase, is a rheophilic migratory cyprinid fish that occurs in the drainage basins of the Black Sea, southern North Sea, and Baltic Sea.It inhabits moderate-to fast-flowing, medium-sized to large rivers with rock or gravel beds and grazes on benthic algae using a characteristic horny layer on the lower lip.Owing to their preference for benthic of seminatural fishways (Figure 4) at the 400 m long Flaucher site (Figure 4), and the construction of near-natural river-bottom ramps (Figure 4).Special care was also taken to reproduce a near-natural waterscape, and one island (Willow Island) (Figure 4) was created to satisfy citizens' esthetic demands.Long-term monitoring confirmed that the good chemical and biological (macrozoobenthos) conditions required for the re-establishment of C. nasus had been achieved.However, electrofishing showed no recovery of C. nasus despite the potential for rapid recovery of the species [64][65][66].

Biology of Chondrostoma nasus
Chondrostoma nasus (L.), named the common nase, is a rheophilic migratory cyprinid fish that occurs in the drainage basins of the Black Sea, southern North Sea, and Baltic Sea.It inhabits moderate-to fast-flowing, medium-sized to large rivers with rock or gravel beds and grazes on

Biology of Chondrostoma nasus
Chondrostoma nasus (L.), named the common nase, is a rheophilic migratory cyprinid fish that occurs in the drainage basins of the Black Sea, southern North Sea, and Baltic Sea.It inhabits moderateto fast-flowing, medium-sized to large rivers with rock or gravel beds and grazes on benthic algae using a characteristic horny layer on the lower lip.Owing to their preference for benthic algae, nase often dwell in shallow, light-flooded habitats, where they can be easily seen.C. nasus is a potamodromous species that needs diverse and very closely connected habitats during its lifecycle.According to literature statements and expert descriptions, the nase population performs different types of movement during the year [64,65,[67][68][69][70][71][72][73][74], which are also confirmed for the Isar River [65].At the end of winter, adult nase often migrate more than 10 km upstream to reach shallow water to feed and recover after the migration and to wait until the environmental conditions are optimal for spawning.When the water temperature reaches 12 • C (between March and May), the adults perform a short migration to nearby well-oxygenated spawning areas in the main channel or tributaries.After spawning, they migrate to summer habitats downstream in the main channel.Before the first freeze-up of winter, the adults move to nearby wintering habitats.Around two weeks after spawning, eggs hatch and larvae drift to shallower habitats near the spawning area.Juveniles (1-3 years old) group in shallow water, e.g., slow-flowing anabranches, but gradually migrate to nearby summer habitats as adults (daily migration to seasonal migration).Accordingly, six habitats (Figure 5  types of movement during the year [64,65,[67][68][69][70][71][72][73][74], which are also confirmed for the Isar River [65].At the end of winter, adult nase often migrate more than 10 km upstream to reach shallow water to feed and recover after the migration and to wait until the environmental conditions are optimal for spawning.When the water temperature reaches 12 °C (between March and May), the adults perform a short migration to nearby well-oxygenated spawning areas in the main channel or tributaries.After spawning, they migrate to summer habitats downstream in the main channel.Before the first freezeup of winter, the adults move to nearby wintering habitats.Around two weeks after spawning, eggs hatch and larvae drift to shallower habitats near the spawning area.Juveniles (1-3 years old) group in shallow water, e.g., slow-flowing anabranches, but gradually migrate to nearby summer habitats as adults (daily migration to seasonal migration).Accordingly, six habitats (Figure 5 and Table 1) have been identified by local fish biologists and cross-validated against the abovementioned literature [64,65,67-75].Table 1.List of habitats associated with lifecycle stages of Chondrostoma nasus (L.) described in terms of their physical characteristics.C. nasus is a sensitive species that is locally threatened by morphological deterioration, e.g., damming and the destruction of spawning sites [64].It is protected by the Berne Convention on the Conservation of European Wildlife and Natural Habitats [76].According to the Red List of the International Union for Conservation of Nature and Natural Resources, C. nasus is classified as a species of "least concern" [75].Historically, all lifecycle stages were found in the Isar River in Munich, although spawning areas were mostly located in the tributaries [65].A population that indicates a good ecological status should reach 20 adults of C. nasus per 100 m river section [77].In the period  C. nasus is a sensitive species that is locally threatened by morphological deterioration, e.g., damming and the destruction of spawning sites [64].It is protected by the Berne Convention on the Conservation of European Wildlife and Natural Habitats [76].According to the Red List of the International Union for Conservation of Nature and Natural Resources, C. nasus is classified as a species of "least concern" [75].Historically, all lifecycle stages were found in the Isar River in Munich, although spawning areas were mostly located in the tributaries [65].A population that indicates a good ecological status should reach 20 adults of C. nasus per 100 m river section [77].In the period from 1995 to 2012, the species' abundance dropped from 40 to fewer than 5 individuals per 100 m in the river sections upstream and downstream from Munich [77], and tributaries of the Danube near Munich (Isar River, Inn River, Vils River, etc.) lost 41% of adult fish of reproductive age (>30 cm long, 4 to around 20 years old) [77].In 2012, no C. nasus were found in the investigated urban section of the Isar River.
C. nasus are able to recolonize restored river stretches in a very short time span [65].Studies on the Danube showed that reconnection of habitats using a nature-oriented scheme was successful, with 46% of the source species pool present in the study area after only two months [66].Reconnection between the main river channel and spawning areas also enabled high recruitment of juveniles only two years after restoration [78].

The Substrate
The substrate properties were determined by field measurements in 2013 and were assumed to be constant over time.Because the water at mean low discharge (MNQ) was shallow and clear, we used a classic survey procedure by boat based on a 5-meter grid (Figure 6) to characterize the substrate (N = 1628).Nine substrate types were visually distinguished on the basis of the grain size of the dominant component: 1.
Organic matter or detritus 2.
Rock or concrete.

The 2D Hydro-Morphological Model
To investigate the conflict between a suitable habitat for fish and intensive recreational uses, we established a 2D hydro-morphological model for 8 kilometers of the restored stretch of the Isar River, namely, from the Großhesseloher Bridge (48°4′29.59″N, 11°32′25.83″E) to Museum Island (48°7′41.42″N, 11°34′46.88″E).For generation of the grid and pre-and postprocessing of the simulated data, SMS software (Surface Modeling System, Aquaveo, Provo, Utah, USA) was employed, while for the hydraulic simulations the software Hydro_AS-2D version 3, was used, which solved the shallow-water equation using a finite-volume discretization.The digital elevation model of this area was based on cross-sectional data provided by the local water authority (Wasserwirtschaftsamt München, http://www.wwa-m.bayern.de/).In addition, local topographical measurements were conducted to consider topographical details.The topographical data were

The 2D Hydro-Morphological Model
To investigate the conflict between a suitable habitat for fish and intensive recreational uses, we established a 2D hydro-morphological model for 8 kilometers of the restored stretch of the Isar River, namely, from the Großhesseloher Bridge (48 • 4 29.59"N, 11 • 32 25.83" E) to Museum Island (48 • 7 41.42"N, 11 • 34 46.88"E).For generation of the grid and pre-and postprocessing of the simulated data, SMS software (Surface Modeling System, Aquaveo, Provo, Utah, USA) was employed, while for the hydraulic simulations the software Hydro_AS-2D version 3, was used, which solved the shallow-water equation using a finite-volume discretization.The digital elevation model of this area was based on cross-sectional data provided by the local water authority (Wasserwirtschaftsamt München, http://www.wwa-m.bayern.de/).In addition, local topographical measurements were conducted to consider topographical details.The topographical data were interpolated to an unstructured grid of 142,000 elements for the study region (Figure 7) and 6000 elements for the reference site.The average spatial resolution was in both reaches 8.0 m in flow direction and 3.0 m in lateral direction, which represented a recommended ratio between length and width of elements.The model was extensively calibrated and validated by comparing measured water levels (provided by the local water authority, Wasserwirtschaftsamt München, http://www.wwa-m.bayern.de/)and simulated water levels [79] for four discharges ranging from 65 m 3 /s to 782 m 3 /s by adapting the roughness values to verify adequate model performance.The roughness values were adapted towards best model performance for two discharges (calibration) and verified for the other two discharges (validation).The deviations between measured and simulated water levels ranged from +8.15 cm to −0.21 cm with a mean deviation of 1.68 cm, which proved to be an adequate model performance (R 2 = 0.98) for both low-and high-flow conditions.The simulated scenarios were: 1.
While scenarios 2-5 were the most representative of the common situation, the other scenarios may inform about stressful situations.Scenario 1 simulated major dryness and enabled the estimation of dewatering risk as fish stranding.Scenarios 4-9 informed about the existence of refuge for the fishes during flood events, which depend of the diversity of the floodplain form.This also informed about the recreational potential of the area.Since the frequency of perturbations are important for resilience capacity of the species, we investigated different flood events from HQ1 to HQ100.

Recreational Pressure
Only expert knowledge can provide information on the response of species to recreational pressure [80].This is commonly acquired by on-site data collection, such as user counting and While scenarios 2-5 were the most representative of the common situation, the other scenarios may inform about stressful situations.Scenario 1 simulated major dryness and enabled the estimation of dewatering risk as fish stranding.Scenarios 4-9 informed about the existence of refuge for the fishes during flood events, which depend of the diversity of the floodplain form.This also informed about the recreational potential of the area.Since the frequency of perturbations are important for resilience capacity of the species, we investigated different flood events from HQ1 to HQ100.

Recreational Pressure
Only expert knowledge can provide information on the response of species to recreational pressure [80].This is commonly acquired by on-site data collection, such as user counting and surveys of user distribution [81][82][83].Riverine recreational activities on the Isar River are water-based from May to October, e.g., boating and swimming, and land-based throughout the year, e.g., lying in the sun, walking, and cycling [84].In accordance with common practice in studies of recreational intensity [34,36,[81][82][83]85], user pressure was evaluated by counting users on a limited number of sampling days.A preliminary user survey during a sunny nonworking day in June at two sites inside the study river section showed that maximum user numbers may be found between noon and 2 p.m. in the south of the study area and between noon and 4 p.m. in the north of the study area (Figure 8).Therefore, users were counted along 8 kilometers of restored river stretch during 10 sunny days (3 during spring, 3 during summer, 3 during autumn, and 1 during winter) between noon and 3 p.m., beginning in the south of the study area.We counted all the recreational users we encountered while walking along the river.This was carried out by using a counter and by differentiating land-based and water-based activities in 10-m intervals (Figure 6), which were delimited using 10-m ropes on both sides of the river, either on the embankment (5-m broad sections) or in the water (between the embankment and the middle of the aquatic area, i.e., 12-35 m broad sections).In the water, users roamed through sections with lengths of 50-100 m by swimming and wading.User densities were classified into five ranks of recreational user density (Table 2; for the ranking procedure, see Section 3.3).User intensity maps were digitized using SMS 10 and acted as an additional parameter for the modeling procedure when testing the impact of users on the availability of fish habitats (see objectives (ii) and (iii) in Section 1).

Habitat Suitability Model
The suitability of habitats for C. nasus was modeled using CASiMiR.This software was designed to determine the suitability of habitats for target species using hydraulic and morphological characteristics.The CASiMiR procedure uses the three main parameters that determine fish habitat preferences [48,86]: water depth, flow velocity, and bed substrate type.Using them, CASiMiR calculates the habitat suitability using three indicators: weighted usable area (WUA), hydraulic habitat suitability index (HHS), and the Habitat Suitability Index (HSI).Both the WUA and the HHS represent functions that relate the habitat suitability to the flow regime.The HHS is obtained by dividing the WUA by the total wetted area, which leads to an index ranging from 0 to 1.The HHS thus eliminates the influence of the size of the wetted area and enables a direct comparison between scenarios [87]; for example, between scenarios for a single study area with different discharge levels.Highly suitable habitats were defined as those with an HSI >0.6.CASiMiR uses a multivariate fuzzy logic approach to link these abiotic attributes with the habitat requirements of fish.Therefore, the overlapping fuzzy magnitudes of the descriptive physical properties are formulated in terms of linguistic categories, i.e., "very high", "high", "medium", "low", and "very low" (Table 2).This approach has proven to be an excellent modeling technique for ecological purposes because the overlapping fuzzy sets allow researchers to deal with uncertain and imprecise information, which commonly occurs in ecological investigations [88].Because little published information on quantified habitat preferences of C. nasus on the Isar River near Munich is available, the physical limits of the categories were partly based on expert knowledge.This concerned, in particular, the behavior of C. nasus when facing manmade stressors such as walkers, swimmers, and large groups of humans on the embankment.Therefore, individual interviews were conducted with seven experts, i.e., fish biologists belonging to the Chair of Aquatic Systems Biology of the Technical University of Munich (N = 2), the Bavarian Water Agency (N = 1), and the Bavarian State Research Center for Agriculture (N = 2), as well as with two fishermen from the NGO Isarfischer.They were asked to describe the C. nasus habitats on the Isar River near Munich considering all the possible combinations of the variables (Table 2) to establish fuzzy habitat suitability sets (Table 3 and Table S1), and to define the response of the species to recreational uses.Despite the fact that we interviewed the biologists in a nondirective manner, namely without list of possible answers, no major deviations between the interviewees were observed.The greatest difference between the interviewees regarded the limits of "low" velocity, namely, 0.5-0.7 or 0.4-0.6 m/s.We considered the deviating entries as the tolerance limits of the values.To validate the resulting habitat quality rules, a reference site with an established population of C. nasus was used.The reference site was located 5 kilometers downstream of the restored reach (Figure 1).According to the last survey performed by the Bavarian fish monitoring authorities, a robust population of C. nasus remains at this site [65,77].The reference site was close to the restored  [0.6 to 0.9] Very high >0.9

Habitat Suitability Model
The suitability of habitats for C. nasus was modeled using CASiMiR.This software was designed to determine the suitability of habitats for target species using hydraulic and morphological characteristics.The CASiMiR procedure uses the three main parameters that determine fish habitat preferences [48,86]: water depth, flow velocity, and bed substrate type.Using them, CASiMiR calculates the habitat suitability using three indicators: weighted usable area (WUA), hydraulic habitat suitability index (HHS), and the Habitat Suitability Index (HSI).Both the WUA and the HHS represent functions that relate the habitat suitability to the flow regime.The HHS is obtained by dividing the WUA by the total wetted area, which leads to an index ranging from 0 to 1.The HHS thus eliminates the influence of the size of the wetted area and enables a direct comparison between scenarios [87]; for example, between scenarios for a single study area with different discharge levels.Highly suitable habitats were defined as those with an HSI >0.6.CASiMiR uses a multivariate fuzzy logic approach to link these abiotic attributes with the habitat requirements of fish.Therefore, the overlapping fuzzy magnitudes of the descriptive physical properties are formulated in terms of linguistic categories, i.e., "very high", "high", "medium", "low", and "very low" (Table 2).This approach has proven to be an excellent modeling technique for ecological purposes because the overlapping fuzzy sets allow researchers to deal with uncertain and imprecise information, which commonly occurs in ecological investigations [88].Because little published information on quantified habitat preferences of C. nasus on the Isar River near Munich is available, the physical limits of the categories were partly based on expert knowledge.This concerned, in particular, the behavior of C. nasus when facing manmade stressors such as walkers, swimmers, and large groups of humans on the embankment.Therefore, individual interviews were conducted with seven experts, i.e., fish biologists belonging to the Chair of Aquatic Systems Biology of the Technical University of Munich (N = 2), the Bavarian Water Agency (N = 1), and the Bavarian State Research Center for Agriculture (N = 2), as well as with two fishermen from the NGO Isarfischer.They were asked to describe the C. nasus habitats on the Isar River near Munich considering all the possible combinations of the variables (Table 2) to establish fuzzy habitat suitability sets (Table 3 and Table S1), and to define the response of the species to recreational uses.Despite the fact that we interviewed the biologists in a nondirective manner, namely without list of possible answers, no major deviations between the interviewees were observed.The greatest difference between the interviewees regarded the limits of "low" velocity, namely, 0.5-0.7 or 0.4-0.6 m/s.We considered the deviating entries as the tolerance limits of the values.To validate the resulting habitat quality rules, a reference site with an established population of C. nasus was used.The reference site was located 5 kilometers downstream of the restored reach (Figure 1).According to the last survey performed by the Bavarian fish monitoring authorities, a robust population of C. nasus remains at this site [65,77].The reference site was close to the restored section; it did not support recreational uses, and it has a high value for spawning activities.A 2D hydro-morphological model of a 2500-m long reference site, namely a 60,000 elements mesh, was created using the same procedure as for the 2D hydro-morphological model of the restored area (presented in Section 3.1.2).We performed the expert-based habitat modeling procedure on the reference site with the existing C. nasus population.The model outcomes showed suitable physical habitats for C. nasus (Figure 9).This result validated the habitat quality fuzzy rules and set of the model.The proportion of suitable wetted area supporting each habitat slightly varies between the reference and the restored sites.This variation can be explained by the length difference between the sites.
The particular aims of this study were to integrate the survey of recreational pressure as supplementary input data into the CASiMiR interface and to integrate the influence of recreational users on habitat suitability into the fuzzy sets and rules.Our interviews with fish biologists revealed the following: (a) C. nasus is for most of the year a very shy species that is easily scared by walkers at the riverside.Even moderate recreational pressure, namely, two or more users within the area nearest to the wetted area (10-m long river section and 5-m cross section), has a negative impact on the riverside habitat used by C. nasus (juveniles and adults during the summer); (b) C. nasus is rapidly scared away by fishermen walking in the river, and even single swimmers and inflatable boats have a high impact on habitat preference; (c) Remote users (far from the wetted area) have either no or a very low, indirect impact on C. nasus; (d) If well-suited habitats are no longer usable owing to disturbances, C. nasus may, within a narrow range of tolerance, use deeper habitats with higher flow velocities; however, this behavior may cause physiological stress; (e) The species is fearless during the annual reproduction migration and the pre-reproduction period.section; it did not support recreational uses, and it has a high value for spawning activities.A 2D hydro-morphological model of a 2500-m long reference site, namely a 60,000 elements mesh, was created using the same procedure as for the 2D hydro-morphological model of the restored area (presented in Section 3.3.2).We performed the expert-based habitat modeling procedure on the reference site with the existing C. nasus population.The model outcomes showed suitable physical habitats for C. nasus (Figure 9).This result validated the habitat quality fuzzy rules and set of the model.The proportion of suitable wetted area supporting each habitat slightly varies between the reference and the restored sites.This variation can be explained by the length difference between the sites.The particular aims of this study were to integrate the survey of recreational pressure as supplementary input data into the CASiMiR interface and to integrate the influence of recreational users on habitat suitability into the fuzzy sets and rules.Our interviews with fish biologists revealed the following: (a) C. nasus is for most of the year a very shy species that is easily scared by walkers at the riverside.Even moderate recreational pressure, namely, two or more users within the area nearest to the wetted area (10-m long river section and 5-m cross section), has a negative impact on the riverside habitat used by C. nasus (juveniles and adults during the summer); (b) C. nasus is rapidly scared away by fishermen walking in the river, and even single swimmers and inflatable boats have

Model Analysis
The model was run in two different ways: (a) exclusively using the tool for modeling physical habitat suitability applied to the 8 kilometers of studied river stretch (response to objective (i) in Section 1) and (b) adding "recreational pressure" as an additional component to identify potential conflicts between recreational uses and ecological quality (response to objectives (ii) and (iii) in Section 1).
We modeled the habitats in 108 scenarios, which comprised combinations of six habitat types (type 1 for Adults spawning, type 2 for Larvae, type 3 for Juveniles, type 4 for Adults in winter, type 5 for Adults in pre-reproduction phase, and type 6 for Adults in summer and autumn), nine discharge levels (NQ, MNQ, MQ, HQ1, HQ2, HQ5, HQ10, HQ50, and HQ100), and two modalities (with and without user pressure).The outputs of the model for each of the 108 scenarios were threefold: (a) Habitat suitability maps, which were composed of the hydro-morphological grid, in which each 5-m element had a calculated habitat suitability index (HSI).The HSI is the most common index for describing biological responses to abiotic attributes and represents the suitability of a habitat for a target species and life stage.The HSI has scalar values between 0 and 1, where 1 at a given grid element represents the most suitable habitat and 0 represents the most unsuitable habitat.(b) A table listing the surface of the wetted area for each HSI value.(c) A table listing the weighted usable area (WUA) and hydraulic habitat suitability index (HHS) for each scenario.
First, to estimate the physical habitat suitability post-restoration for each of the six habitat types and for each discharge level, we first performed descriptive analyses of the WUA and HHS values.Second, to define potential conflict areas, the geographical distribution of suitable physical habitats and discharge level was indicated on physical habitat suitability maps, and the spatial distribution of users was described.Finally, to investigate the influence of recreational pressure on physical habitat suitability, we performed Mann-Whitney-Wilcoxon and Ansari-Bradley tests to compare the medians and variances of the HSI at 200 randomly chosen elements of the grid for the 108 scenarios.Analyses were performed on all scenarios.However, some Figures and Tables presented only the results at mean annual discharge.This choice had been made to increase the readability.Considering Figure 2, mean annual discharge was considered as the common situation.Results for the other scenarios are available in the supplementary material.

Physical Habitat Suitability
Physical habitat modeling results in terms of WUA in function of discharge (Figure 10 and Table 4) showed that highly suitable habitats (HIS > 0.6) exist for adults during the winter, summer, and pre-reproduction periods and for juveniles.Few areas with suitable habitat quality for larvae and spawning adults occurred at the restored Isar River stretch (Table 4).The areas with suitable habitats were located at the site Flaucher and at near-natural river bottom ramps with honeycomb-shaped structures.On the scale of the whole modeled river section, the HHS varied between 0.01 and 0.55, according to all investigated lifecycle stages and discharge levels (Table 5).The highest rates were found at mean low discharge and mean discharge for juveniles (HHS = 0.4 and 0.3), pre-reproduction adults (HHS = 0.3 and 0.2), and adults during summer (HHS = 0.2 and 0.2).The highest HHS was reached for juveniles and pre-reproduction adults at minimum discharge-55% and 35%, respectively, of the wetted area was highly suitable.Suitable physical habitats for pre-reproduction adults and for juveniles largely decreased with an increase in discharge.Suitable habitats for larval development and spawning activities were rare at any discharge level (<10% of the wetted area).

Spatial Distribution of Suitable Habitats
At NQ, MNQ, and MQ, highly suitable physical habitats for adults in winter and summer (>80%) were mostly located in the southern two-thirds of the study area and predominantly along the western riverside.The Flaucher site had high physical habitat suitability for adult fish during summer (Figure 11f) and winter (Figure 11d), except in the seminatural fishway (more than 80% of which had an HSI of <0.1).River sections (100-200 m) upstream of weirs had low HSI values (mostly <0.1).The area surrounding Willow Island (800 m upstream and 1000 m downstream) also had low HSI values (>90% of the area had an HSI of <0.1) at NQ and MNQ, but HSI values increased in this area at higher discharges.During flood events, suitable habitats for adults were relocated to the western riverside, namely, in the flooded recreational meadows.For NQ, MNQ, and MQ, suitable physical habitats for juveniles and pre-reproduction adults were found in the whole restored section, even in backwaters upstream of weirs.The areas surrounding the Flaucher (Figure 11c,e) and Willow Island had high density of habitat patches with an HSI >0.9.During flood events, no relocation of habitats for juveniles was found.Suitable physical habitats for larvae and spawning were very limited within the whole restored river section for all investigated discharge levels, even within the Flaucher (Figure 11a,b).Surprisingly, single spots with medium HSI values (0.4-0.7) for both habitats were located on the near-natural river bottom ramps.
Sustainability 2018, 10, x FOR PEER REVIEW 13 of 26 Figure 2, mean annual discharge was considered as the common situation.Results for the other scenarios are available in the supplementary material.

Physical Habitat Suitability
Physical habitat modeling results in terms of WUA in function of discharge (Figure 10 and Table 4) showed that highly suitable habitats (HIS > 0.6) exist for adults during the winter, summer, and pre-reproduction periods and for juveniles.Few areas with suitable habitat quality for larvae and spawning adults occurred at the restored Isar River stretch (Table 4).The areas with suitable habitats were located at the site Flaucher and at near-natural river bottom ramps with honeycomb-shaped structures.On the scale of the whole modeled river section, the HHS varied between 0.01 and 0.55, according to all investigated lifecycle stages and discharge levels (Table 5).The highest rates were found at mean low discharge and mean discharge for juveniles (HHS = 0.4 and 0.3), pre-reproduction adults (HHS = 0.3 and 0.2), and adults during summer (HHS = 0.2 and 0.2).The highest HHS was reached for juveniles and pre-reproduction adults at minimum discharge-55% and 35%, respectively, of the wetted area was highly suitable.Suitable physical habitats for pre-reproduction adults and for juveniles largely decreased with an increase in discharge.Suitable habitats for larval development and spawning activities were rare at any discharge level (<10% of the wetted area).

Table 4.
Wetted area and relative percent of the wetted area with not suitable (HIS < 0.1), low suitability (HSI between 0.1 and 0.3), medium suitability (HSI between 0.3 and 0.6), and high suitability (HIS > 0.6) at mean annual discharge.even in backwaters upstream of weirs.The areas surrounding the Flaucher (Figure 11c,e) and Willow Island had a high density of habitat patches with an HSI >0.9.During flood events, no relocation of habitats for juveniles was found.Suitable physical habitats for larvae and spawning were very limited within the whole restored river section for all investigated discharge levels, even within the Flaucher (Figure 11a,b).Surprisingly, single spots with medium HSI values (0.4-0.7) for both habitats were located on the near-natural river bottom ramps.

Spatial and Temporal Distribution of Recreational Pressure
The number of users inside the study area increased at the end of April and decreased at the end of summer (Figure 12).The mean number of users during the study period was 599 (standard deviation (sd) = 211) per daily observation, among whom 23 (sd = 20) were swimmers.The maximum value, which was observed on a June day, was 774 users, including 60 swimmers.In summer, namely late July and August, the number of users dropped mainly due to major rain and flood events.
User density varied inside the study area (Figure 13) (from 0 to 0.24 users per square meter).Most of the users were found on the eastern side (>90%).Recreational user pressure was high in five areas: (1) 600 m up-and downstream from Willow Island during the summer user density reached, on average, 0.12 users per square meter over the whole area.Single spots were found with 0.2 users

Spatial and Temporal Distribution of Recreational Pressure
The number of users inside the study area increased at the end of April and decreased at the end of summer (Figure 12).The mean number of users during the study period was 599 (standard deviation (sd) = 211) per daily observation, among whom 23 (sd = 20) were swimmers.The maximum value, which was observed on a June day, was 774 users, including 60 swimmers.In summer, namely late July and August, the number of users dropped mainly due to major rain and flood events.per square meter, whereas the minimum value of 0.01 users per square meter was recorded in winter.Users mostly sat in small groups on the embankment close to the water and were present only on the eastern side of the river; (2) At the designed stairs, 0.18 (maximum) and 0.04 (minimum) users per square meter were found in summer and in winter, respectively; however, no swimmers were observed; (3) Around the Flaucher, user density in summer reached 0.24 users per square meter, and many users were found in the water (on average, N = 21); (4) On river-bottom ramps, user density locally increased to 0.12 users per square meter (maximum value); (5) From the Großhesseloher Bridge to 1000 m downstream, user density was lower (0.04 users per square meter on average), but a quarter of the users were swimmers.The western side of the river within the whole restored section was under low recreational user pressure.User density varied inside the study area (Figure 13) (from 0 to 0.24 users per square meter).Most of the users were found on the eastern side (>90%).Recreational user pressure was high in five areas: (1) 600 m up-and downstream from Willow Island during the summer user density reached, on average, 0.12 users per square meter over the whole area.Single spots were found with 0.2 users per square meter, whereas the minimum value of 0.01 users per square meter was recorded in winter.Users mostly sat in small groups on the embankment close to the water and were present only on the eastern side of the river; (2) At the designed stairs, 0.18 (maximum) and 0.04 (minimum) users per square meter were found in summer and in winter, respectively; however, no swimmers were observed; (3) Around the Flaucher, user density in summer reached 0.24 users per square meter, and many users were found in the water (on average, N = 21); (4) On river-bottom ramps, user density locally increased to 0.12 users per square meter (maximum value); ( 5) From the Großhesseloher Bridge to 1000 m downstream, user density was lower (0.04 users per square meter on average), but a quarter of the users were swimmers.The western side of the river within the whole restored section was under low recreational user pressure.many users were found in the water (on average, N = 21); (4) On river-bottom ramps, user density locally increased to 0.12 users per square meter (maximum value); ( 5) From the Großhesseloher Bridge to 1000 m downstream, user density was lower (0.04 users per square meter on average), but a quarter of the users were swimmers.The western side of the river within the whole restored section was under low recreational user pressure.2).

Influence of Recreational Users on Availability of Habitats for C. nasus
An overall comparison of the WUA values showed that, at MQ, 8% of highly suitable habitats for C. nasus may be lost owing to recreational pressure.Statistical analyses of the 200 randomly chosen points showed that the quantity of suitable habitats and the habitat quality of the restored river section for C. nasus decreased when recreational pressure was integrated into the evaluation procedure for juveniles, spawning adults, and larvae, but no significant differences were found for adults during summer, autumn, winter, and pre-reproduction periods (Figure 14).Potentially highly suitable habitats for juveniles, spawning activities, and larval development lost important parts of their surface areas due to recreational pressure (9%-15.5%,20%-25%, and 71%-76%, respectively), whereas potentially highly suitable habitats for adults during summer, autumn, winter, and pre-reproduction periods lost only a small part of their surface areas due to recreational pressure (0%-9%, 0%-4%, and 0%-3%, respectively).The largest losses of potentially suitable habitats for juveniles owing to recreational pressure were found on the low-flow secondary arm flowing on the side of Willow Island, on the eastern part of the Flaucher, and just downstream from the Großhesseloher Bridge (Figure 15).For example, the Flaucher site lost more than 60% of its highly suitable habitats for juveniles when recreational pressure was integrated into the model.All the relatively suitable habitats for spawning activities and larval development were located in areas with high user and swimmer densities.
An overall comparison of the WUA values showed that, at MQ, 8% of highly suitable habitats for C. nasus may be lost owing to recreational pressure.Statistical analyses of the 200 randomly chosen points showed that the quantity of suitable habitats and the habitat quality of the restored section for C. nasus decreased when recreational pressure was integrated into the evaluation procedure for juveniles, spawning adults, and larvae, but no significant differences were found for adults during summer, autumn, winter, and pre-reproduction periods (Figure 14).Potentially highly suitable habitats for juveniles, spawning activities, and larval development lost important parts of their surface areas due to recreational pressure (9%-15.5%,20%-25%, and 71%-76%, respectively), whereas potentially highly suitable habitats for adults during summer, autumn, winter, and prereproduction periods lost only a small part of their surface areas due to recreational pressure (0%-9%, 0%-4%, and 0%-3%, respectively).The largest losses of potentially suitable habitats for juveniles owing to recreational pressure were found on the low-flow secondary arm flowing on the eastern side of Willow Island, on the eastern part of the Flaucher, and just downstream from the Großhesseloher Bridge (Figure 15).For example, the Flaucher site lost more than 60% of its highly suitable habitats for juveniles when recreational pressure was integrated into the model.All the relatively suitable habitats for spawning activities and larval development were located in areas with high user and swimmer densities.

Physical Habitat Suitability
The modelling results show that the restored river section provides four out of six suitable physical habitats that C. nasus requires during the stages of its lifecycle, i.e., habitats for adult fish during summer, autumn, and winter, pre-reproduction periods as well as for juveniles.However, the results of the habitat suitability model also showed three major issues that may explain the nonrecovery of the species.
First, we found that the restored river section has very limited suitable habitats for spawning activities and larval development.However, the physical habitat model suggests that near-natural

Physical Habitat Suitability
The modelling results show that the restored river section provides four out of six suitable physical habitats that C. nasus requires during the stages of its lifecycle, i.e., habitats for adult fish during summer, autumn, and winter, pre-reproduction periods as well as for juveniles.However, the results of the habitat suitability model also showed three major issues that may explain the nonrecovery of the species.First, we found that the restored river section has very limited suitable habitats for spawning activities and larval development.However, the physical habitat model suggests that near-natural manmade elements, e.g., near natural river-bottom ramps with honeycomb-shaped structures, may have locally positive effects on physical habitat suitability for C. nasus, even if they did not succeed in recreating highly suitable habitats.River-bottom ramps result from the removal of linear low weirs and have honeycomb-shaped structures constructed with natural large rocks to disperse the hydraulic head over a short distance but with a gentle slope.These provide a high diversity of microhabitats with different velocities, depths, and substrata, e.g., for sensitive benthic invertebrates and fish in gravel bar rivers [89].However, although slight improvements were found, the functionality of the created habitats remains to be proven for the nase.Cyprinidae species require recently deposited clean fine-gravel substratum for spawning and larval development, and further field research is needed to test if the sediments in the ramps are sufficiently clean to allow attachment of nase eggs.Another constraint for the re-establishment of the nase population in the Isar River is that tributaries that historically served as additional spawning sites are today dammed, canalized, and partly buried.Our findings suggest that the priority of the restoration/conservation strategy should be a long-term solution based on the improvement of both the migration potential (namely, the re-establishment of longitudinal and transverse connectivity between the main channel, secondary arms, and tributaries) and the restoration of physical habitats for spawning activities and larval development.
Second, the restored Isar River did not provide sufficient slow-flowing anabranches suitable for juvenile C. nasus during high-discharge events.According to the Isar River seasonal flow regime, major flood events may happen during reproduction or shortly afterward.Larvae and young fish are particularly sensitive to drift during flood events [65,70,72,73].One anabranch has been recreated (Willow Island), but the goals of the restoration were set to maintain broad flat recreational grassland instead of meandering the riverbed and creating an undulating riverscape.As a consequence, the water velocity remains too high, and the sediment dynamics in the floodplain are too low to provide suitable habitats for recruitment.Even if the direct impact of the recreation is removed during the reproduction period and the early development stage of juveniles, the choice of the project goals, namely, to "provide flat grassland for recreational use", results in a low suitability of the habitats for these critical life stages.As a consequence, the population of C. nasus in the Isar River is ageing [65].
Natural recruitment is the most promising option for conservation [77], as previous attempts of restocking C. nasus have failed [74].Our findings suggest that the restoration/conservation strategy should include the improvement of both the migration potential (namely, the re-establishment of longitudinal and transverse connectivity between the main channel, secondary arms, and tributaries) and the restoration of physical habitats for spawning activities and larval development by permitting more sediment dynamics.We predict that if recovery occurs, the southern section will have higher densities than the northern one due to its higher habitat suitability values.

Conflicts between Wildlife and Recreational User Pressure
Our results suggest that different conflicts occur according to lifecycle-specific habitat types of C. nasus and the periodic preferences of recreational users.We found that the highest user densities partly occurred in suitable habitats for juvenile C. nasus, i.e., low-flow zones and the submerged tops of gravel bars.However, the modelled habitats for juveniles were rather widely distributed over the restored river stretch, and the overall impact of users may remain limited.User densities were also high in the rare potential spawning areas, e.g., river-bottom ramps, causing a significant decrease in the WUA.However, the spawning period (beginning of May) may occur shortly before the user density ranges from high to very high (May to October), but overlapping of the areas used by C. nasus and recreationists remains likely.Our results suggest that users have a low impact (a) on adult fish during winter, because user density is significantly lower and C. nasus inhabits deeper water than during summer and autumn [64,65,69]; (b) on adult fish during summer and autumn, because users observed in our study were mostly located on the eastern side and suitable physical habitats were on the western side; and (c) on adults during the pre-reproduction period, because fish biologists described C. nasus as less sensitive to disturbance at this time.We suggest that conflicts may be avoided if user distribution can be wisely influenced by the design of restoration and by adequate guidance of users.
Tools for managers attempting to mitigate recreation-wildlife interactions are regulations, public education, and a user-management plan.First, we generally advise against the formulation of regulations forbidding the recreational use of larger riverine areas within the city because positive public support for sustainable management and restoration projects is strongly driven by perception, communication, possibilities to participate, and the usefulness of the restoration outcomes.However, in case of great importance, e.g., last occurrence of a population, restrictions may be the only solution.Acceptance of this restriction is efficient if the abstinence of using a protected site becomes a collective activity of respect towards an acknowledged heritage.Second, informing and educating the public can be addressed by an attitude-change strategy.Urban river beaches are prime sites for creating encounter places of man and nature, which help to re-establish emotional linkages, create motivation, and change decisional values in favor of maintaining or restoring ecosystem integrity [2].A study on the Danube showed that only 40% of recreational users were aware that wildlife is disturbed by recreational activities, e.g., off-trail walking [90].Educational work to increase the awareness and knowledge of the public seems to be a sustainable tool for avoiding conflicts between users and wildlife.However, a strategy for attitudinal change may be ineffective.A study showed that only 5% of recreationists change their behavior after educational work.Furthermore, a study in the United States showed that people are becoming less supportive toward fish and restoration efforts targeting fish species.Accordingly, we consider that this strategy may be very useful to provide sites for environmental education in cities and enable attitude change even beyond the visited site.However, highly sensitive species may not be able to sustain the disturbances at highly frequented sites.The targets for restoration projects need to be well defined, explained, and chosen by the public [15].Third, managing flows of visitors applying a user-management plan may relocate users and separate sections of the river between recreationists and ecological refuges for wildlife [91].Studies have shown that user preferences for waterscapes are driven by both visual values [23,[92][93][94][95], e.g., parameters of fascination, vividness, and naturalness [93][94][95], and usability [96].In our case study, it may be relevant to create suitable habitats for wildlife on the western side of the river, which is almost inaccessible to citizens because of the topography of the area.However, user preferences are complex because every individual may perceive landscapes in a different way and may be differently affected by environmental and societal stimuli.
Interest in the prediction of user density and distribution for implementing a motivational management strategy that regulates users has grown considerably in recent years.User distribution maps or models should be integrated into habitat suitability models to create a predictive tool to help choose the best restoration design, to define refuges for species, and to investigate future scenarios (climate change, industrialization, increases in user density, etc., see next paragraph).However, this area of scientific inquiry involves a broad range of disciplines and their respective fields of knowledge.More studies on recreational user preferences and recreational impacts on wildlife remain to be carried out to provide a strong baseline for the design of user-management plans based on an integrated framework for coexistence between recreationists and biodiversity.

Benefits of Integrating User Pressure into the Physical Habitat Model
The modeling of fish habitats using CASiMiR provided significant results for evaluating the quality of the physical habitats for C. nasus in the studied river section.A novelty was the inclusion of recreational user pressure as a new parameter in the habitat-modeling procedure, which showed clear differences from the results of modeling physical habitats only.The modeling procedure helped to identify conflict areas and the lifecycle stage (juveniles) that may be most affected by users.However, the method still has some limitations.Considering the precision of the assessment of recreational pressures, we have to state that occasional observations have revealed much higher figures of riverbank-users than the results from a standardized procedure presented here.In order to reduce bias from the high fluctuation of users (driven by working time and weather conditions), user surveys should be more frequent and long term.We suggest the use of drones or aerial photographs to reduce cost.Another critical point is the evaluation of species tolerance limits, namely, tolerated user intensities and frequencies.In this study, the habitat preferences and impacts of recreationists were estimated by expert statements of fish biologists because little scientific literature was available concerning these points.Despite the fact that expert evaluations of ecological quality may be as trustworthy as assessments made by experimental field investigations [97], future research designed specifically to evaluate the impact of recreationists on wildlife would be helpful.
Our study has demonstrated the existence of habitats with suitable flow velocities, depths, and substrata for all lifecycle stages of C. nasus historically observed in this Isar River section, and it suggests that the reduction of adequate spawning and juvenile habitats by physical destruction and by user pressure (on the juvenile habitats) is responsible for the absence of recruitment.However, it does not explain fully the absence of C. nasus in the restored Isar River section.Other habitat variables such as temperature, food sources, and predators as well as habitat availability on the scale of the catchment could be taken into account to complete the picture.

Conclusions
This study delivered a model of the suitability of habitats for the indicator species Chondrostoma nasus (L.) on the restored Isar River in Munich (Germany), including, as a novel feature, the impact of recreational users as a supplementary parameter in the habitat suitability model.The findings of the study are threefold.First, the research showed that the restoration of the studied river section succeeded in providing physically suitable habitats for all lifecycle stages of the species in general.It also showed that the establishment of habitats for spawning and larval development, the creation of refuges for sensitive lifecycle stages during flood events, and the re-establishment of a wide diversity of habitats has not been achieved within the northern part of the restored river section.Second, the model showed that recreational users may cause decreases in highly suitable habitats for sensitive lifecycle stages of C. nasus in identified zones of conflict, especially during summer, when recreational user density is high.We suggest to combine two approaches: (a) providing sections of riverine beaches as encounter sites of man and nature in the sense of the River Culture Concept [2] with (b) a user-management plan directing flows of visitors out of highly sensitive sites by making them less or not accessible.. Third, the developed method enabled an integrative evaluation of habitat quality for the target species.The modeling approach, which combines hydro-morphological-biological modelling with user pressure, has been shown to be a useful tool for increasing the efficiency of ecological restoration design, defining the best location for the re-introduction of sensitive species, and supporting user-management plans in urban contexts.Further research about the behavior of fish facing recreational pressures, variables such as temperature and the food-web composition, and habitat availability on the catchment scale is needed to identify the ultimate causes for the failure to re-establish C. nasus.We argue that there is an urgent need for a shift in the design of restoration projects and their evaluation toward holistic and interdisciplinary approaches that embrace the principles of socioecological systems.

Supplementary Materials:
The following are available online at http://www.mdpi.com/2071-1050/10/6/1747/s1.Table S1: Fuzzy sets and rules for the studied habitats, namely for adults (during the winter, the summer and autumn, during pre-reproduction, and spawning), larvae, and juveniles C. nasus at the Isar river (Germany) considering the input variables (i.e.velocity, water depth, Substratum, recreational pressure).Table S2: Wetted area and relative percent of the wetted area with not suitable (HSI<0.1),low suitability (HSI between 0.1 and 0.3), medium suitability (HSI between 0.3 and 0.6) and high suitability (HSI>0.6)at MQ. Table S2 presents all the scenarios.Table S3: Weighted Usable Area (WUA), Hydraulic Habitat Suitability index (HHS), and Mean Habitat Suitability Index (mean HSI) value for each habitat and all scenarios.

26 Figure 1 .
Figure 1.Locations of Munich, of the study area, and of the reference site for the habitat suitability model.

Figure 2 .
Figure 2. Discharges of the Isar River in Munich recorded at the gauging station inside the study area (a) in 1999 and (b) between 1999 and 2006 by the Bavarian Water Agency.HQ100, HQ10, and HQ1: 100-, 10-year, and 1-year maximum discharge; MQ: mean discharge; NQ: minimum discharge.

Figure 1 . 26 Figure 1 .
Figure 1.Locations of Munich, of the study area, and of the reference site for the habitat suitability model.

Figure 2 .
Figure 2. Discharges of the Isar River in Munich recorded at the gauging station inside the study area (a) in 1999 and (b) between 1999 and 2006 by the Bavarian Water Agency.HQ100, HQ10, and HQ1: 100-, 10-year, and 1-year maximum discharge; MQ: mean discharge; NQ: minimum discharge.

Figure 2 .
Figure 2. Discharges of the Isar River in Munich recorded at the gauging station inside the study area (a) in 1999 and (b) between 1999 and 2006 by the Bavarian Water Agency.HQ100, HQ10, and HQ1: 100-, 10-year, and 1-year maximum discharge; MQ: mean discharge; NQ: minimum discharge.

Figure 3 .
Figure 3. Photograph of the restored Isar River section in Munich (Aude Zingraff-Hamed in August 2015).

Figure 4 .
Figure 4. Localization and photograph of the cited site of importance at the restored Isar River section in Munich (Photograph took by Aude Zingraff-Hamed in March 2013).

Figure 3 .
Figure 3. Photograph of the restored Isar River section in Munich (Aude Zingraff-Hamed in August 2015).

Figure 3 .
Figure 3. Photograph of the restored Isar River section in Munich (Aude Zingraff-Hamed in August 2015).

Figure 4 .
Figure 4. Localization and photograph of the cited site of importance at the restored Isar River section in Munich (Photograph took by Aude Zingraff-Hamed in March 2013).

Figure 4 .
Figure 4. Localization and photograph of the cited site of importance at the restored Isar River section in Munich (Photograph took by Aude Zingraff-Hamed in March 2013).

Table 1 .
List of habitats associated with lifecycle stages of Chondrostoma nasus (L.) described in terms of their physical characteristics.

Figure 6 .
Figure 6.Grid design used for the substratum and user surveys.

Figure 6 .
Figure 6.Grid design used for the substratum and user surveys.

Figure 7 .
Figure 7. Mapping of the (a) Velocity and of the (b) Water Depth at mean annual discharge (displaying 3.5 kilometers of the 8 kilometers of the restored river section including the Flaucher).

Figure 7 .
Figure 7. Mapping of the (a) Velocity and of the (b) Water Depth at mean annual discharge (displaying 3.5 kilometers of the 8 kilometers of the restored river section including the Flaucher).

Figure 8 .
Figure 8. User numbers on two emerged gravel bars (100-m long and 20-m wide) during a sunny nonworking day in June at two sites: (A) in the north of the study area (48 • 7 37.85" N, 11 • 34 40.15"E) and (B) in the south of the study area (48 • 4 45.07" N, 11 • 32 29.25"E).

Table 2 .
List of variables used for the fuzzy sets.m (±5 cm) to 0.2 m (±5 cm) medium 0.2 m ±5 cm) to 0.5 m (±10 cm) high 0.5 m (±10 cm) to 1.15 m (±25 cm) very high Above 1.25 m (±25 cm) Velocity very low 0 to 0.4 m/s (±0.1 m/s) low 0.5 m/s (±0.1 m/s) to 0.7 m/s (±0.1 m/s) medium 0.75 m/s (±0.1 m/s) to 0.9 m/s (±0.1 m/s) high 1 m/s (±0.15 m/s) to 1.5 m/s (±0.25 m/s) very high Above 1.75 m/s (±0.25 m/s) Substratum grain size low Organic material, detritus, silt, clay, loam, sand (<6 mm) medium Gravel from 6 mm to 12 cm high Large stones (12-20 cm) very high Boulders (>20 cm), rock User pressure (per 10-m long section) very lowNo users low One user on the riverbank and no users in/on the water medium Two or more users on the riverbank or one user in/on the water high One user in/on the water and >1 users on the riverbank very high More than two users in the water and more than two users on the riverbank

Figure 9 .
Figure 9. Habitat suitability map of the Reference site for (a) adults during summer, (b) adults during winter, (c) juveniles, (d) adults during the pre-reproduction period, (e) larvae, and (f) spawning adults at mean annual discharge.

Figure 9 .
Figure 9. Habitat suitability map of the Reference site for (a) adults during summer, (b) adults during winter, (c) juveniles, (d) adults during the pre-reproduction period, (e) larvae, and (f) spawning adults at mean annual discharge.

Figure 10 .
Figure 10.WUAs: (a) for spawning adults, (b) for larvae, (c) for juveniles, (d) for adults during winter, (e) for adults during the pre-reproduction period, and (f) for adults during summer and autumn.Flow acronyms defined in Section 3.1.2.

Figure 11 .
Figure 11.Habitat suitability map of the Flaucher site for (a) adults during summer, (b) adults during winter, (c) juveniles, (d) adults during the pre-reproduction period, (e) larvae, and (f) spawning adults at annual mean discharge.

Figure 11 .
Figure 11.Habitat suitability map of the Flaucher site for (a) adults during summer, (b) adults during winter, (c) juveniles, (d) adults during the pre-reproduction period, (e) larvae, and (f) spawning adults at annual mean discharge.

Figure 12 .
Figure 12.Plot of user numbers inside the whole study area on the 10 investigated days highlighting the number of swimmers.

Figure 12 .
Figure 12.Plot of user numbers inside the whole study area on the 10 investigated days highlighting the number of swimmers.

Figure 12 .
Figure 12.Plot of user numbers inside the whole study area on the 10 investigated days highlighting the number of swimmers.

Figure 13 .
Figure 13.User density map of the Isar River in Munich indicating on the eastern side spots with high to very high user densities separated by long stretches with very low user densities, and on the western side, mostly very low user densities (categories are as defined in Table2).

Figure 14 .
Figure 14.Weighted usable area (WUA) in square meters of the restored river stretch at annual mean discharge for the six habitats of Chondrostoma nasus L. with and without the inclusion of recreational pressure.

Figure 14 . 26 Figure 15 .
Figure 14.Weighted usable area (WUA) in square meters of the restored river stretch at annual mean discharge for the six habitats of Chondrostoma nasus L. with and without the inclusion of recreational pressure.

Figure 15 .
Figure 15.Habitat suitability map for juveniles of Chondrostoma nasus at mean annual discharge considering or ignoring recreational pressure in the case of (a) the Flaucher area, (b) Willow Island, and (c) the area near the Großhesseloher Bridge (HSI = Habitat Suitability Index).

Table 1 )
and have been identified by local fish biologists and cross-validated against the abovementioned literature [64,65,67-75].

Table 3 .
Example of a fuzzy rule set describing the habitat requirements of Chondrostoma nasus (L.) for adult fish in the Isar River during summer.(FullTableS1).
Rule 4: IF velocity "Very low" AND depth "Medium" AND substratum "Very high" AND recreational pressure "Medium" THEN HSI "Low" VL = very low; L = low; M = medium; H = high; and VH = very high.

Table 5 .
TableS2presents all the scenarios.Weighted Usable Area (WUA), Hydraulic Habitat Suitability index (HHS), and Mean Habitat Suitability Index (mean HSI) value for each habitat and scenarios MNQ (mean low discharge), MQ (mean annual discharge), and HQ (Annual mean maximum discharge), with and without users.TableS3presents all the scenarios.