1. Introduction
In the last century, and in particular, between the 1930s and the early 2000s, most of the Italian rivers underwent remarkable morphological changes. Two types of channel adjustments were observed, i.e., riverbed lowering, which is on average in the order of 6–7 m, and channel narrowing, in some cases up to 90% [
1,
2,
3,
4,
5]. These processes had serious repercussions on the ecological conditions of the river corridor, which includes channel, floodplain, and recent terraces. In particular, a significant functionality alteration of the lateral continuity processes of liquid and solid flows was observed [
6,
7,
8].
The causes of such channel adjustments (incision and narrowing) are still debated. However, most authors invoke a dominant control by various types of anthropogenic interventions, such as land-use changes, channelizations, constructions of dams, and in-channel sediment mining [
1,
9,
10,
11,
12]. Other authors [
13,
14] argue that climate change and related modifications in flood frequency and magnitude are the main control factors of short-term channel adjustments. According to these authors, human interventions would simply modulate the role played by such natural factors, often further exacerbating their action. Recently, Magliulo et al. [
15] showed that relatively small rivers affected by negligible in-channel human disturbances underwent short-term morphological changes fully comparable to those experienced by more anthropized rivers. Thus, they conclude that the interactions between multidecadal rainfall changes and land-use changes at the basin scale, together with variations in floods magnitude and frequency, are the main controlling factors, at least for relatively small rivers.
In this framework, the assessment of the morphological and ecological status of streams is of fundamental importance for their management and protection. To this aim, the European Water Framework Directive 2000/60/EC (WFD) is the reference legislation. The directive introduced the assessment of the hydromorphological features of streams, in addition to the biological and chemical-physical ones, to define the river ecological status [
16]. The directive envisaged the achievement, for all natural water bodies, of a “Good Ecological Status”, which reflects the conditions of good biodiversity and chemical-physical status. In practice, the ultimate aim is to identify the conditions that determine the natural equilibrium of the system and evaluate the human disturbances, which often negatively impact the equilibrium state [
17].
In 2008, the Italian National Institute for Environmental Protection and Research (ISPRA) promoted a research program with the aim of developing a comprehensive methodology for the hydromorphological analysis of the Italian streams. This methodology, called IDRAIM (stream hydromorphological assessment, analysis, and monitoring system) [
17,
18,
19,
20], pursued an integrated analysis of both the morphological quality and the risks associated with the channel dynamics, aimed at implementing both the framework directive on water and the EU Floods Directive [
21]. The Italian Ministry of the Environment also issued the monitoring and classification standards of the Water Framework Directive [
22], introducing the Morphological Quality Index (IQM), whose assessment is part of IDRAIM methodology, as a protocol for the evaluation of the morphological quality of each reach of the investigated rivers [
19].
It is widely accepted that a given classification system works better as a guide to stream management when it is developed for specific physiographic regions. The IDRAIM framework takes into account pre-existing geomorphological approaches developed in other countries but also the specific Italian context in terms of channel adjustments and human pressures. Furthermore, a basin-level perspective characterizes the method, and the analysis explicitly takes into account the temporal component in terms of reconstruction of the channel morpho-evolutionary trajectory as a basis for the interpretation of current river conditions and future trends [
19,
20]. Moreover, the method is based on an integration of GIS analysis, remote sensing, and field survey [
23,
24].
Several pre-existing studies focused on the hydromorphological quality assessment of streams [
16,
25,
26,
27,
28,
29,
30], however, in peninsular Southern Italy, such studies are lacking. This is an important gap due to the proven impact of anthropogenic activities on the recent geomorphological changes of the river channels in this area [
4,
31,
32]. The almost total absence of such kind of studies in this area can be partly explained by taking into account that, until the 2010s, an almost total lack of knowledge existed about the short-term morphological changes of rivers in Southern Italy. Such knowledge, as stated above, is of fundamental importance for the morphological quality assessment of a given river. Furthermore, the morphological quality assessment represents a first but fundamental step in the correct evaluation of the flood risks. However, the hydromorphological analysis should support and not replace classical methods of flood risk assessment based on previous inundation ranges and digital elevation models. In fact, hydromorphological analysis, such as that carried out using the IDRAIM method, considers, in an integrated way, the processes (such as riverbed lowering, lateral mobility, etc.) that control the risks associated with river morphological dynamics and induced flooding. The assessment of such processes, which are both of linear (i.e., at the reach scale) and areal type (i.e., by defining zones of different morphological river dynamics), should be integrated with the results provided by the methodologies traditionally used for the analysis of flood risk (i.e., hydraulic, topographic, and geotechnical analyses).
Due to the particular climatic and topographic conditions, the rivers of Southern Italy are often subject to catastrophic floods [
33], with deep impacts on the weak local economies. Moreover, for this reason, the application of morphological quality assessment methodologies to the rivers of Southern Italy appears to be desirable.
This study aims to provide the first contribution to fill this gap by means of the hydromorphological characterization of one of the major rivers of Southern Italy, the Calore River (
Figure 1), which suffered a catastrophic flood on 14–15 October 2015 [
34]. The analysis was carried out by using the IDRAIM methodology, which allowed, among others, the identification of the reaches with the most critical issues.
Study Area
The Calore River basin is located in Southern Italy, within the axis of the Campanian sector of the Southern Apennines fold-and-thrust belt (
Figure 1), which is the result of compressional tectonics mainly active during the Tertiary [
35].
The basin has an extension of 3050 km
2 and is located between 41°30′ N and 40°46′ N latitude and 14°27′ E and 15°20′ E longitude. The highest altitude (i.e., Accellica Mount, Massif of the Picentini Mountains) is 1810 m a.s.l., while the lowest is 55 m a.s.l. (i.e., at the confluence into the Volturno River) [
5].
The climate in the Calore River basin is temperate-warm, with mean annual rainfall of ∼1150 mm, mainly concentrated during autumn and winter [
36]. The rainiest month is November, while August is the driest. Mean annual temperature is 15.1 °C.
The Calore River flows into morphostructural depressions infilled with alluvial, volcanic, and slope deposits [
37]. In some cases (e.g., lower Calore River valley), Quaternary alluvial deposits form very thick sedimentary successions [
38].
The Calore River basin can be subdivided into three main physiographic units (
Figure 1c): (a) The Apennine Mountain range, located along the basin border, (b) the low-altitude hills, which occupy a large portion of the inner part of the basin, and (c) the alluvial valleys.
The Apennine Mountain range is shaped into limestone and dolostone, Jurassic to Cretaceous in age. Steep slopes, locally deeply dissected by tributary streams, characterize this physiographic unit. At the top of these slopes, karstified remnants of ancient erosional landscapes (paleosurfaces or Paleosuperfici Auct.) are present.
The hilly relieves are shaped into terrigenous deposits, mainly Miocene to Pliocene-aged. The hills are gently sloping and gently rolling [
5]. Frequent and clear evidence of severe water erosion and mass movements affect these relieves. Such erosional processes control the contribution of solid material to the river, which affects its morphological variations, whose assessment is fundamental for the MQI calculation that we carried out in this study. Landslides of various types are also extremely common in the Calore River basin [
39].
The alluvial valleys are shaped into Quaternary sediments and are mainly of morphostructural origin [
38,
40,
41]. Thick successions of alluvial and slope deposits, often forming coalescent and/or telescopically arranged alluvial fans, connect the slopes to the alluvial plains. The latter are characterized by different orders of river terraces. Locally, structural terraces shaped into pyroclastic deposits (i.e., Campanian Ignimbrite) and travertines are present.
The land-use in the Calore River basin [
42] is dominated by agricultural lands (64% of the total area) and, subordinately, by forests (23%). In the last 60 years, the studied basin experienced an increase in forested areas at the expense of agricultural lands. Most of the land-use changes in the study area occurred between the 1960s and the 1990s, followed by a period of substantial stability. This likely led to a decrease in soil erosion and sediment supply to the rivers between 1960 and 1990 [
42], which is confirmed by recent studies dealing with short-term river channel adjustments [
4,
5,
43], followed by a phase of overall stability since the 1990s [
42].
The Calore River, which flows through the basin, is the main left tributary of the Volturno River and one of the major rivers of the Campanian Apennines. It starts from Accellica Mount and flows down for some 110 km before joining the Volturno River. The mean annual flow discharge, measured just upstream from the confluence into the Volturno River, is ~30 m
3/s [
31]. The Calore River valley pattern is strongly controlled by tectonics, as evidenced by the four successive segments that form at 90° angles to each other [
37].
Of the entire Calore River, only the middle-lower course was investigated for the calculation of the MQI (
Figure 1c). The investigated reach is comprised between the points of coordinates 41°11′03″ N–8°27′45″ E and 41°00′03″ N–8°56′51″ E. The reach located upstream from the latter point was excluded from the analysis for reasons that will be better clarified in the “Materials and Methods” section.
2. Materials and Methods
The present research required: (i) The characterization of the fluvial system, (ii) the assessment of past morphological evolution, and (iii) the evaluation of the current conditions of the Calore River. We used the QGIS 3.24 software (OSGeo, Beaverton, OR, USA) to obtain and process all the data.
A preliminary step in the application of the methodology consisted in the manual digitization of all the environmental features that were necessary for the MQI calculation. In particular, we manually digitized: (i) The active channel of the Calore River, (ii) the centreline (i.e., the line that connects the points equally distanced from the riverbanks within the active channel), (iii) the valley axis (which was necessary for the sinuosity calculation, as will be better explained below), (iv) the fluvial bars, whose spatial distribution is necessary to assess the channel morphology), and, (v) the anthropogenic structures (e.g., bridges, check-dams, mining sites, walls).
For the digitization of the above-listed elements, we used, as background in GIS environment: (i) 1:5000-scaled digital color orthophotos, produced by the Campania Region in 2014 and georeferenced in the UTM-WGS84 coordinate system, (ii) a 10 m-cell size Digital Terrain Model (DTM), freely available and downloadable from the INGV-SURFLT website [
44], (iii) 1:25,000-scaled topographic maps from 1955, produced by the Italian Military Geographical Institute (IGMI), in *.ecw format, georeferenced in the UTM-WGS84 coordinate system. The latter were necessary to assess the variations in channel width and morphology that occurred between 1955 and 2014. The data were then validated in field.
2.1. Characterization of Fluvial System
2.1.1. Physiographic Units, River Segments (RS), and Confinement Conditions Definition
As a first step, we subdivided the Calore River basin into three physiographic units, i.e., the Apennine Mountain range (AMR), low-altitude hills (LAH), and alluvial valleys (AV) (
Figure 1), whose features were described in the “Study area” section. According to Rinaldi et al. [
20], we defined the physiographic units based on both elevation and geological data. More precisely, the 600 m a.s.l. contour line (as derived from the 10 m-cell DTM mentioned above) marked the boundary between AMR and LAH physiographic units, while the upper lithological limit of the Quaternary-aged alluvial deposits marked the boundary between AV and LAH physiographic units. We digitized the latter boundary from five Geological Maps at 1:100,000 scale, produced by Servizio Geologico d’Italia, i.e., (i) Sheet N. 173—“Benevento” [
45], (ii) Sheet N. 174—“Ariano Irpino” [
46], (iii) Sheet N. 162—“Campobasso” [
47], (iv) Sheet N. 185—“Salerno” [
48], and (v) Sheet N. 186 “S. Angelo de’ Lombardi” [
49]. These maps were previously scanned, imported into the QGIS 3.24 software, and georeferenced in the UTM-WGS84 coordinate system.
The intersections between the boundaries of the physiographic units and the Calore River allowed a first subdivision of the studied river into River Segments (RS). Such subdivision was further refined based on river confinement conditions. To this aim, two parameters were calculated, i.e., the Confinement Degree (CD) and the Confinement Index (CI).
To calculate the
CD (which expresses the lateral confinement of the river), we overlaid the riverbanks (that we previously digitized at the 1:5000 scale) on the Geological Maps (1:100,000-scaled) listed above. Thus, we split each riverbank into “confined” segments (i.e., bordered by slopes, landslides, alluvial fans, Pleistocene-aged river terraces) and “unconfined segments” (i.e., bordered by alluvial plain shaped into Holocene-aged alluvial deposits). We refined such subdivision by overlaying the previously split riverbanks on the Digital Terrain Model to identify more precisely the inner boundary of the alluvial plain. Then, the GIS software automatically calculated the length of the “confined” riverbanks, as well as the total riverbank length. We calculated the
CD of each RS by means of the following formula:
We considered the following possible cases, based on the
CD values [
23]:
Confined channel: More than 90% of the riverbanks are bordered by slopes, landslide bodies, alluvial fans, or inactive river terraces;
Partly confined channel: A part of the riverbanks ranging between 10% and 90% is bordered by slopes, landslide bodies, alluvial fans, or inactive river terraces;
Unconfined channel: Less than 10% of the riverbanks is bordered by slopes, landslide bodies, alluvial fans, or inactive river terraces.
To calculate the
CI (which expresses the longitudinal river confinement), we first digitized the alluvial plain, based on the spatial distribution of the Holocene-aged alluvial deposits. We then refined the digitized boundary by overlaying it on the DTM to take into account also the local topography. For each RS, the GIS software automatically calculated both the channel area and the alluvial plain area. Then, the “Fluvial Corridor” QGis tool automatically traced a centerline for both the alluvial plain and the river channel. The length of both centerlines was automatically calculated by the software. By dividing the alluvial plain area by the length of the alluvial plain centerline, we calculated the alluvial plain mean width. Similarly, we calculated the river channel mean width [
50]. This allowed us to calculate the
CI of each RS as follows:
According to Rinaldi et al. [
20], we used the following classes of
CI:
Because the Calore River is a single-thread channel [
31], k = 5 [
20].
Finally, we defined the overall confinement of each RS based on both
CD and
CI, according to
Table 1 [
20]:
2.1.2. Channel Morphology Assessment
To define the channel morphology, both the sinuosity index [
51] and morphological pattern were determined.
The sinuosity index calculation requested splitting the entire centerline of the Calore River into shorter segments, i.e., the previously identified RSs. Then, the QGis software automatically calculated the length of each RS centerline. The straight segment joining the extremes of each RS represented the valley axis. Similar to the centerline, the length of the valley axis for each RS was automatically calculated by the software. Finally, we calculated the sinuosity index (
SI) for each RS as follows:
To complete the characterization of the channel morphology, we also took into account the distribution of the fluvial bars within the active channel. The fluvial bars were previously digitized on the orthophotos of 2014. This allowed us to define the channel morphological pattern according to Rinaldi [
52]. By combining the SI values and the channel morphological pattern, we defined the channel morphology class for each RS according to Rinaldi et al. [
20] (
Table 2).
Unfortunately, it was not possible to assess the channel morphology for the entire Calore River. In fact, in sectors located upstream from the point of coordinates 41°00′03″ N and 8°56′51″ E, the fluvial bars were not correctly identifiable on the orthophotos due to both the negligible channel width and bar extension and, especially because they were partly or totally hidden by a very dense arboreal overhanging riparian vegetation. Thus, it was not possible to assess both the channel morphology and, subsequently, the MQI of these sectors (
Figure 1).
2.1.3. Final Definition of the River Reaches (RR)
In this step, the river reaches previously defined on the basis of confinement (both lateral and longitudinal), sinuosity, and channel morphological pattern were, where necessary, further subdivided on the basis of: (i) Presence of important hydrological discontinuities (e.g., confluences with the main tributaries), (ii) excessive length (i.e., >5 km), (iii) presence of sediment extraction sites, (iv) further variations in channel morphology. The final river reaches (hereinafter, RRs) identified in this step represented the basic elementary unit functional to the following analyses.
2.2. Assessment of the Past Morphological Evolution
In this phase, the previously obtained data were analyzed in the framework of a morpho-evolutionary framework. This was aimed at deriving an assessment of the current river morphological conditions also based on the knowledge of past conditions. In particular, we analyzed the adjustments in: (i) Channel pattern, (ii) channel width, (iii) riverbed level.
2.2.1. Assessment of Adjustments in Channel Pattern
For this parameter, in the GIS environment, we made a comparison between the morphology of the different RRs from 1955 and 2014.
According to Rinaldi et al. [
20], we considered the following possible cases:
Unchanged morphological configuration;
Limited or moderate variation in morphological configuration, i.e., between contiguous types;
Intense morphological configuration variation, i.e., between non-contiguous types.
Contiguous types refer to similar morphological configurations (e.g., a riverbed that passes from meandering to sinuous), while non-contiguous types include variations between very different morphologies (e.g., from multiple channels to single channels).
We have, therefore, classified the changes in the morphological configuration based on whether or not the channel switched to a different morphology than the initial one in the time interval under consideration (i.e., 1955–2014;
Table 3).
2.2.2. Assessment of Adjustments in Channel Width
To assess the adjustments in channel width, we first divided the Calore River channel from 1955 into the same RRs identified for the channel in 2014. The QGis software automatically calculated the area of each RR from 1955. Similarly, we also split the centerline of the Calore River from 1955 according to the RRs identified in 2014. The centerline length of each reach was automatically calculated by the software. We then calculated, for each RR from 1955, the channel mean width (
RRCMW1955) as follows [
50]:
We then compared
RRCMW1955 and
RRCMW2014 (i.e., the active channel mean width of each RR in 2014, which we previously calculated to determine the river confinement). We calculated the variation in width, expressed in meters, between
RRCMW1955 and
RRCMW2014 (∆
CMW1955–2014) as follows:
Finally, we expressed the observed variation in active channel mean width in the considered period as a percent of the
CMW in 1955:
The calculated ∆
CMW1955–2014 were classified according to Rinaldi et al. [
20], i.e.,
Channel unchanged or affected by limited variations: ∆CMW1955–2014 ≤ 15% of RRCMW1955;
Moderate narrowing or widening: 15% < ∆CMW1955–2014 < 35% of RRCMW1955;
Intense narrowing or widening: ∆CMW1955–2014 > 35% of RRCMW1955.
2.2.3. Assessment of Adjustments in Bed-Level
We derived data about bed-level changes from the literature. In particular, Magliulo et al. [
31] demonstrated that, between 1955 and 1998, the Calore River underwent a riverbed lowering ranging between 3 and 5 m. Such an amount allows classifying (according to Rinaldi et al. [
20]) the riverbed lowering experienced by the Calore River as “intense” (i.e., ranging between 3 and 6 m). Magliulo et al. [
5] report that no significant changes affected the Calore River and its main tributaries in the 2000s. Thus, the amount of the incisions reported by Magliulo et al. [
31] for the period 1955–1998 is very likely valid also for the entire 1955–2014 timespan.
2.3. Assessment of the Current Conditions of the Calore River
Morphological Quality Index (MQI) Assessment
According to Rinaldi et al. [
20], we assessed three different hydromorphological features, i.e.,
Geomorphic functionality, which is based on the observations of fluvial landforms and processes.
Artificiality, which expresses presence, frequency, and spatial continuity of anthropogenic structures and/or activities, which could affect the considered river morphological parameters;
River channel changes.
We then calculated the Morphological Alteration Index (MAI) as follows [
20]:
where S
tot is the sum of the scores and S
max is the maximum possible deviation for the given stream typology (i.e., the sum of the class C scores, reported in
Table 4, for all the questions applicable to the study case). The MAI is equal to 1 in a completely unaltered river, while it is equal to 0 in a completely altered river.
The final step was the assessment of the Calore River morphological quality by means of the calculation, for each RR, of the Morphological Quality Index (MQI), as follows:
Based on the MQI values for each RR, we defined the Morphological Quality Classes according to Rinaldi et al. [
20]. These classes are reported in
Table 5.
We finally assessed the overall morphological state of the Calore River by calculating the Mean Morphological Quality Index (
MMQI) by means of the following formula [
20]:
where:
MQIRR-i is the Morphological Quality Index of the RR “I”;
LRR-i is the length of the RR “I”.
4. Discussion
The observation of the Indicators of Functionality (FI;
Table 11), Artificiality (AI;
Table 12), and Channel Adjustments (CAI;
Table 13) clearly shows that the Calore River is characterized by a negligible overall artificiality (
Figure 8b), severe channel adjustments experienced in the time-span 1955–2014 (
Figure 9), and moderate functionality (
Figure 7).
The negligible overall artificiality is demonstrated by the fact that, according to Rinaldi et al. [
20], the highest possible value (Smax) for AI in a given RR is 72. In the case of Calore River, the highest value of the AI is 23/72, which characterizes only one RR out of thirty (i.e., RR 4.3.2;
Table 12). In the remaining RRs, the AI is lower than 5/72 in 11 reaches (i.e., 37% of the total number of RR) and lower than 10/72 in 24 RRs out-of-thirty (80%).
In contrast, the Calore River was affected by severe channel adjustments. In fact, according to Rinaldi et al. [
20], the CAI Smax is 24. In the case of the Calore River, two RRs display a very high value of such indicator (i.e., 20/24, in both RR 4.2.1 and RR 6.3.3). More generally, 57% of the total number of RRs display a high value (i.e., >17/24) of the CAI. Still more significant is the fact that 97% of the considered RRs display a CAI > 10/24.
As regards the Functionality, the Smax of the associate indicator (FI) is 54 [
20]. In the case of the Calore River, the FI highest value is relatively low, i.e., 16/54, which was obtained in two RRs out of thirty, only (i.e., RR 4.1.2 and 6.1.1). Similarly, a very low value of FI (i.e., 5/54) was obtained for one RR only. Thus, most part of the RRs display moderate-low values of FI. More precisely, 63% display FI values ranging from 6/54 to 10/54. Such percent increases to 90% if we consider intermediate FI values ranging between 6/54 and 15/54.
Based on what was stated above, the morphological quality of the Calore River seems to be negligibly affected by the type and intensity of channel adjustments experienced between 1955 and 2014, while a deeper influence is probably exercised by the degree of artificiality of the river. In fact, the morphological quality of the Calore River is generally good (as highlighted by the generally high MQI values;
Table 14), notwithstanding the intense or very intense channel adjustments experienced by the Calore River, the latter highlighted by the high values of the CAI (
Table 13). In contrast, the good morphological quality (expressed by the high values of the MQI;
Table 14) of the Calore River is well reflected in the generally low values of the AI (
Table 12), which indicates a reduced artificiality of the river. As regards Functionality, it seems to display an overall moderate influence on the morphological quality of the Calore River, compared to Artificiality. However, the river morphological quality seems to be quite sensitive to variations in Functionality, as proved by the fact that two RRs with relatively lower (i.e., “moderate”) morphological quality (i.e., 4.1.2 and 6.1.1;
Table 14) out of four are characterized by the highest values of FI (
Table 11), which indicates a relatively more compromised geomorphic functionality.
Among the features controlling the morphological quality of the Calore River, check-dams play a relevant role in altering both the sediment discharge in the reach (see A4 in
Table 12) and the longitudinal continuity of sediment and wood solids (see F1 in
Table 11). Obviously, check dams are particularly frequent along the RRs that crosses the towns of Benevento (i.e., RR 4.3.2) and Ponte (RR 6.1.1). This concurs to explain why the morphological quality of such RRs is the lowest. The presence of walls and artificial levees (A6 and A7 in
Table 12) concurs to reduce the river morphological quality in these RRs too.
Further considerations can also be made about the sediment extraction sites (A10 in
Table 12). In fact, in-channel sediment removal is considered one of the main driving forces of channel incision and narrowing [
1,
10]. In the case of the Calore River, the relationship between in-channel sediment mining and channel adjustments (i.e., narrowing and incision) seems to be proved by the fact that the RR with the highest value of CAI (i.e., RR 4.1.2;
Table 13) is limited both upstream and downstream by the two largest sediment extraction sites.
In any case, it should be stressed that the assessment of the morphological quality of the Calore River provided very encouraging results, as the morphological quality resulted “good” (
sensu [
20]) in most parts of the studied river. The morphological quality slightly decreases to “moderate” in four RRs out of thirty only, i.e., in both the urban RRs and in those RRs with sediment extraction sites. However, generally speaking, in the urban RRs, the morphological quality is often poor or very poor, due to a large number of features (mainly anthropogenic) that increase artificiality and negatively affect functionality. Thus, a “moderate” morphological quality in an urban RR, as in the case of the Calore River, can be considered quite satisfying. Furthermore, the improvement of the morphological quality in urban RRs would require the removal of human infrastructures, such as check dams and walls. Such removal is very difficult to imagine also because such infrastructures are aimed at contrasting the negative or catastrophic effects of floods and riverbed degradation.
As regards the effects of sediment removal, it should be reminded that, in Italy, the sediment removal from the riverbed has been unlawful since the 1990s. The presence of sediment extraction sites remains unquestionably an element that increases the artificiality and, accordingly, decreases the river morphological quality. However, avoiding the sediment removal from the riverbed reduces the impacts on the channel changes, namely on channel incision and narrowing, which, in turn, negatively affect the river morphological quality. Obviously, careful monitoring of sediment removal activities is strongly recommended to avoid unlawful sediment removal from the riverbed that could compromise the river morphological quality. Similarly, any anthropogenic activity that could decrease functionality or increase artificiality (
Table 9) should be avoided, if possible, or, alternatively, carefully projected or evaluated, with the aim to preserve the good current morphological quality of the Calore River.
The results obtained in this study are consistent with those provided by other studies that determined the MQI by using the IDRAIM method on Italian rivers [
25,
26,
27,
28,
30]. In fact, all these studies clearly indicate that Artificiality is the indicator that mainly controls the morphological quality of a river, while channel adjustments play a slightly more marginal role. For example, Rigon et al. [
27] found a good or moderate MQI in 86% of the Cordevole River (Northern Italy) and a poor MQI in the remaining 14%, which characterizes the river reaches in which artificial elements are more widespread. Licciardello et al. [
30] obtained less encouraging results for the Dittaino River, as 67% of the river showed high and moderate quality, while the remaining 33% showed poor morphological quality in reaches characterized by the presence of elements of artificiality, such as dams and check-dams.
As a general consideration, it should be stressed that the IDRAIM method unquestionably presents advantages and disadvantages. As regards advantages, the method at issue unquestionably takes into account the consideration and understanding of the geomorphological processes responsible for river functioning. Such understanding can be used for both classification and to support analyses of any intervention and impact, and the design of mitigation measures. Furthermore, the IDRAIM method was developed to fully comply with the WFD requirements. The main disadvantages [
30] unquestionably consist of the strong degree of subjectivity in the scoring of the indexes, which is affected by the degree of experience and training of the operators. Furthermore, MQI assessment needs to define numerous indicators, which sometimes are not suitable for the hydrological conditions of the rivers located in Southern Italy. This is probably because most of the streams used in the testing phase of the procedure are localized in Central and Northern Italy. In fact, most rivers of Southern Italy, especially the smallest ones, are characterized by irregular water flow and harsh hydrological fluctuations that are not frequent in rivers located in central-northern Italy. However, this is not the case with the Calore River. Thus, the results we obtained in this study can be considered fairly reliable.