Temporal and Spatial Variability of Sediment Transport in a Mountain River : A Preliminary Investigation of the Caldone River , Italy

Sediment transport is a key evolution process of rivers and water basins. This process can pose flood hazards to nearby areas. The Eulerian and Lagrangian methods are usually employed to describe sediment transport in mountain rivers. The application of different methods was proposed by scientists to analyze specific aspects of solid transport, however a complete understanding still alludes us. After a brief review of the most common methods, the coupling of three different methods is proposed and tested in order to study sediment dynamics, and its spatial and temporal variability, in mountain rivers. Tracers, painted bed patches, and digital elevation model (DEM) comparisons are used to characterize sediment transport at both a micro-scale short-term and a macro-scale long-term level on a test reach on Caldone River, Italy. Information about travel distance, critical diameters, active width, and morphological evolution was sought. We focused on how water discharge is changing the relationships between different measurement techniques. High discharge events force the channel to behave in a unique way, while low discharge events generate more intrinsic variability. Only measurement technique coupling can overcome this issue. Results are encouraging and show the potential of a mixed Eulerian-Lagrangian approach.


Introduction
Sediment transport is a key evolution process of rivers and watersheds.In addition to floods, sediments represent a hazard that can cause major damage to infrastructure and settlements near riparian zones [1][2][3][4].A deeper knowledge of the sediment sources [5][6][7][8][9][10] is needed to, forecast and control the hazard, understand the underlying physics, and reduce the associated risk.
This condition is common to both mountain and low-gradient streams, but they are different both in sediment and water flow features, due to high variability in sediment size and sediment sources [11].Thus, the application of classic methods, such as sediment traps or turbidity sampling, to estimate sediment transport in slow current, low gradient rivers, provide reliable results.However, these techniques cannot be extended to steep mountainous rivers due to the variation of hydraulic processes in time and space.Seasonal discharge changes are higher than in lowland streams, as is the intensity of sediment transport and debris velocity.The complex condition of the river channel, influenced by boulders, debris, wood, and bedrock [12], causes highly variable bed-load transport rates [13].
Bed-load transport is the mobilization of the largest diameter classes, which cause the most relevant effects in terms of aggradation and impact on structures.
Few studies in the scientific literature have quantitatively measured bed-load transport in mountain basins [23,24].To measure bed-load transport, travel distance, critical diameter, pathways, source, and deposition areas are usually quantified.Methods to measure sediment transport are usually divided into the Lagrangian approach and Eulerian approach, and are mostly derived from laboratory testing [25][26][27][28].The Lagrangian approach refers to following the trajectories of single clasts, and the Eulerian approach refers to characterizing the sediment transport properties at a specific area.Since laboratories have a controlled environment and variables can be held constant (only changing one at a time, while natural rivers exhibit large variability), the exportability of these methods from the laboratory to the field is not straightforward.
Lagrangian methods are mainly focused on tracing pebbles; irrespective of the chosen technology, they can be found and identified in the river by in situ surveys, in order to record their position at any time.The most common identification methods include passive and active radio tags, iron or magnetic pebbles, painted pebbles, and radioactive sources [29][30][31].Tracking pebbles is important to understanding sediment mobility during average and intense flow events [32,33], since this approach is effective at most flow rates.The literature presents many examples, with some authors focusing on the relationship between peak discharge and/or duration [34], particularly in mountainous regions.Different objectives can be addressed through the use of tracers, such as transport rates, transport distances, and more probable pathways [35][36][37].Event-based surveys also provide data on mobilization thresholds, sediment sorting, depth and width of the active layer (defined as the thickness of the river bed that is moving), and source and deposition areas [11].The population of measured displacements, linked to their causative process, can be used in advanced statistical interpretation of the clasts' transport behavior.Some authors [38] noted that tracers placed in the thalweg move 15 to 30 times more than tracers entrapped in bars.This suggests that the river bed displays different velocities, with bar sediments that start moving only when adequate shear stress is reached.
Nevertheless, the use of radio tags presents some technical limitations.Passive tracers cannot be inserted into grains that have major axis <60-70 mm, while active tags need a bigger host pebble due to battery dimension.In this way, tracers represent only the largest fraction of the granulometric curve.Despite this, the collected results could be used for fine tuning or validating the numerical simulation and modeling, in comparing measured and expected transport in terms of distances and mobilized diameters [39].To guarantee a pure Lagrangian approach, tracers are normally sought with a mobile antenna.Nevertheless, examples of using fixed antennas exist and make this approach closer to a Eulerian approach [40,41].Stationary antennas provide continuous monitoring of the study reach and, thus, a larger amount of reliable data with a better understanding of the incipient motion conditions.
The Eulerian approach focuses on measuring all sediments stopping or leaving a controlled volume, while the motion of a single particle is neglected.Different techniques are used for this, including sediment traps, piezoelectric bed-load impact sensors (PBISs), digital elevation model (DEM) comparison, section evolution monitoring, scour chains, painted bed patches, and sedimentation pool bathymetric surveys.Sediment traps and PBISs have the same objective of measuring all sediments crossing a certain point of the bed, independently from their dimensions and, thus, operate at a very small, localized scale.Obviously, there is a huge difference in cost and complexity, but both present some disadvantages: Sediment traps only provide a measure of transport after an event and can be filled during major floods, while PBISs give real-time, high-resolution measurements of passing sediment but work just for >20 mm clasts and requires ad hoc complicated calibration that could worsen the results on short and intense events [42].Some authors enlarged the scale, to measure the relevant sections in a river several times in order to estimate river-bed macro-variations due to sediment transport processes [43].In an effort to improve the results, the DEM of the riverbed is also acquired as an alternative method to evaluate morphologic evolution [8].
All of these methods provide surface transport information, but not of the depth of the active layer.To investigate this aspect, the use of scour chains is commonly employed.Scour chains are metal link chains which are vertically buried in the river bed.After an event, the depth of the active layer can be evaluated as the portion of the chain that has moved.
Some authors [44] propose a method to evaluate the incipient movement of clasts in a riverbed based on the application of paint on patches of the bed.A pre-event and post-event comparison of the granulometric curve determines the clasts' entrainment by the flow; critical diameter determination is then quite straight forward.The critical diameter is commonly defined as the largest grain size moved in an event.This technique overcomes the impossibility to apply the equations developed for low gradient channels, like the Shields formula [45], to steep mountain rivers, due to the higher flow rate variability, and rough bed with coarse sediments.
Different temporal and spatial domains can be evaluated when studying sediment dynamics.One could record data during an event or after each event, and this will be defined as a short-term scale.However, it is possible to enlarge the time span to a medium or long-term scale, one hydrological year, or more.In the same way, the investigated spatial domain could be extremely small.For painted patches or sediment traps, we will speak of the micro-scale, but it can also be enlarged to a portion of the reach, comparing river bed DEM to support studies of morphological evolution trough time.The natural disposition of a river to change its morphology is often neglected, and the observation of evolutionary tendencies is not correctly interpreted or taken into account, particularly in mountain streams [46].The literature analysis highlighted how each method has the best performance in measuring specific processes over a range of time and space domains.Due to the variability in mountain rivers, the integration of different methods is mandatory to describe all the processes happening and their variation in space and time.Banks and bars behave differently from the thalweg, and long-term processes are slightly more predictable than short-term ones.Few authors proposed attempts to integrate approaches instead of focusing on a single process or method [47].
This paper presents preliminary work to test the possibility of coupling and integrating different methods to record sediment dynamics in a mountain river, summing short and long-term, micro and macro-scale observations.We propose an integration of the Lagrangian and Eulerian methods to create a new dataset to improve our understanding of the problem.The Lagrangian method consisted of tracers and the Eulerian method deployed painted patches, while a morphological analysis was used to verify changes, caused by torrent dynamics, observed during the monitoring period.

Materials
The main goal of this research is to collect data about the sediment dynamics of a mountain river, combining different analytical techniques at a micro and macro scale.Micro-scale surveys allow one to describe internal variability and behaviors of different parts of the same stream reach.
We chose the Caldone River basin, a pre-alpine basin in Italy, flowing across the city of Lecco, 50 km north of Milan, Italy.We wanted to fully characterize transport by coupling processes, with the temporal resolution being at the single-event scale.Three approaches were combined to investigate the sediment dynamics: A Lagrangian tracer-based approach, an Eulerian approach based on painted bed patches, and DEM comparison.In Figure 1, a flowchart of the work is presented.The Lagrangian approach employs RFID tracers which were searched for after any event.This is a different method than what was done by many authors [11,32,33,38,48], who used months or years as time intervals between surveys.This technique measures micro-scale short-term dynamics of pebbles in thalweg, while an event-based survey allowed for deeper understanding, since pebbles' movements are linked directly to the causal event, while the displacement database has a higher resolution than yearly time scales.A longer time span forces researchers to account for the sum of events and also the mean values of the controlling variables.The Eulerian approach combines the use of painted bed patches with the creation and comparison of DEM bed models.The aim of the former is to determine incipient motion conditions, hence it is classified as a short-term micro scale analysis [45].DEM comparison allows for the determination of whether long-term macro-scale trends of morphological evolution exist in the river [43].
[45].DEM comparison allows for the determination of whether long-term macro-scale trends of morphological evolution exist in the river [43].

Case Study
The Caldone River hydrographic basin is 24 km 2 , with elevation ranging from 197m to 2118 m above sea level (a.s.l.) at the top of Grigna Meridionale (Figure 2).In the upper part of the basin, above 400 m a.s.l., the Caldone River drains Triassic and dolomitic limestones of the Angolo, Prezzo, and Esino formations.In the middle part, the river cuts through red arenites, pelites, and yellow dolomites of the San Giovanni Bianco formation and Pleistocene glacial deposits, lacustrine and alluvial fan deposits, and debris flows [49].In the lower part of the basin, the river flows through an urbanized floodplain.The average precipitation over the city of Lecco (population 50,000) is about 1400 mm/year, concentrated in late spring and autumn.During summer, storms are common and can be intense, while winter precipitation is less intense, with snow accumulation being present only in the higher part of the basin.In its last kilometer before the outlet into the Lario Lake, the Caldone River flows within a culvert that passes below the town center (Figure 3b).The combination of short runoff time, high slope (from 1.5% to 5.5%), intense sediment transport, and a high water discharge within a densely urban area makes the Caldone River a potential hazard for the city of Lecco.Monitoring devices installed in the basin are provided in Table 1.

Case Study
The Caldone River hydrographic basin is 24 km 2 , with elevation ranging from 197m to 2118 m above sea level (a.s.l.) at the top of Grigna Meridionale (Figure 2).In the upper part of the basin, above 400 m a.s.l., the Caldone River drains Triassic and dolomitic limestones of the Angolo, Prezzo, and Esino formations.In the middle part, the river cuts through red arenites, pelites, and yellow dolomites of the San Giovanni Bianco formation and Pleistocene glacial deposits, lacustrine and alluvial fan deposits, and debris flows [49].In the lower part of the basin, the river flows through an urbanized floodplain.The average precipitation over the city of Lecco (population 50,000) is about 1400 mm/year, concentrated in late spring and autumn.During summer, storms are common and can be intense, while winter precipitation is less intense, with snow accumulation being present only in the higher part of the basin.In its last kilometer before the outlet into the Lario Lake, the Caldone River flows within a culvert that passes below the town center (Figure 3b).The combination of short runoff time, high slope (from 1.5% to 5.5%), intense sediment transport, and a high water discharge within a densely urban area makes the Caldone River a potential hazard for the city of Lecco.Monitoring devices installed in the basin are provided in Table 1.Moreover, two sediment retention basins are present in the last 5 km of the river, with volumes of around 3000 m 3 each.These pools get filled with sediment in about three years.Unfortunately, precise filling data is not available, but the volume of deposited sediments could be determined from pictures and dredging operation reports.Data is available from 1999 to 2018.Sedimentation pools protect the culvert from excessive solid transport, which can result in channel blockage, due to a reduced cross-section of the culvert entrance (Figure 3b).
However, some crucial data was lacking to characterize the basin and the river: Granulometry and water discharge.Granulometry is useful for understanding which fraction of sediment could move under certain flow conditions.Granulometry was determined at 16 points along the river stretch.A photographical technique developed by ETH Zurich [50] was applied, along with sediment size analysis.The results are reported in Table 2.

Subcatchments of the Caldone River
The watercourse consists of two main stretches, upstream and downstream of the confluence with the Grigna River (Figure 2).The upstream segment is the steepest (about 18%) and flows through narrow canyons and valleys with low debris coverage.The river is completely natural and well-confined by fluvial and glacial terraces.The lower part has milder slopes (about 4%) but locally crosses plains, where it forms bends and meanders due to the gentler slope, similar to reach 10 (Figure 2).This segment crosses the city, partially underground, and is rich in structures, levees, and check dams.From this preliminary classification, a further division into 14 sections was possible.The criteria used were adopted from [51], based on the geometrical features of the river, such as the degree of river confinement, river length, width, slope, and river bed shape.From this analysis, three stretches of the river were chosen to represent the entire course: Stretches 6, 7, and 11 (Figure 2).A choice was made to select reaches with different features and which were logistically easier to monitor.The main features are reported in Table 2.
Stretch 6 is about 1 km long, with a mean slope of 5.5%, low confinement, and low presence of artificial structures.This segment represents the most natural part of the river.The segment ends in a big sedimentation pool of about 4500 m 3 (Figure 3a).
Stretch 7 is well-confined and 375 m long, with a mean slope of 3.7%, and is just downstream of the sedimentation pool.This segment contains check dams and concrete levees.It was chosen to represent an intermediate situation between natural and artificial river bed.
Stretch 11 is 800 m long and constrained artificially.The bed is not erodible in most of the reach and its last 600 m are completely underground.The slope is limited to 1.5% and artificial walls confine the river.The analyzed segment has a natural bed which is plain and without bars or sediment accumulations.This stretch was chosen to investigate sediment dynamics within highly urbanized areas.The measured granulometric distribution showed higher coarse fraction, even though the slope is lower.This is due to the upstream presence of a sedimentation pond, which retains most of the debris coming from upstream.The finer fraction has probably been transported to the outlet in previous floods, since there is no fine material coming from upstream.In reach 11, the painted bed patches technique could not be applied, since the bed is completely covered by a water stream, even during low flow periods.

Water Discharge Measurement
Particle image velocimetry (PIV) measures the surface velocity of water using a short video, and calculates the mean velocity of each cell [52].Videos were acquired with GoPro Here 4 cameras at 30 frames per second.The images were processed using self-developed software [53].Twenty-seven measurement cells are drawn in the sections, each cell being 20 cm wide.The surface velocity is correlated to the stream velocity, thanks to two contemporary propeller measures in the river at different water stages.Estimated error in velocity measurements is about 15%, which is an acceptable value compared to those documented in earlier uses of PIV [52,53].This approach enables the construction of a stage-discharge relationship to measure discharge at any moment, based on the recorded water height.It was performed at two different points of the Caldone River, where stage measurements are available [54].All discharge events presented in this paper were calculated using the PIV technique.It was extremely important to link any displacement to its own causing event flow, since we based our observation on a short-term event scale approach.The two selected spots have regular and non-turbulent flow features, to maximize the expected PIV results.

Surveyed Events
We surveyed 18 rain events from July 2016 to November 2017, and data about tracers of event 7 and 11 were lost due to technical issues.The main characteristics of these events are reported in Table 2.The rainfall regimes varied considerably between 2016 and 2017.The rainfall in 2016 was close to average, with a rainfall deficit experienced in winter, while 2017 was a dry year with some large storms in summer, but overall there was a water deficit of about 200 mm with respect to the average of 1400 mm/year.

Lagrangian Method
The applied Lagrangian approach consists of following each pebble along the river, recording each of its movements after any significant rain event.In order to retrieve and identify each pebble in the river, radio-frequency identification (RFID) technology was adopted.Each pebble is equipped with a transponder that, if stimulated by an electromagnetic field generated by an antenna, transmits back its identification number to the antenna.A special script stores the ID, along with a GPS position and retrieval time in a database.It is, thus, possible to keep track of the particle displacements over time [11,31].The pebbles are painted yellow for ease of location and reading their RFID (Figure 4).The technical details of the equipment and pebble preparation are discussed elsewhere [2,54].
The recovery ratio of pebbles and tracers decreases over time.Tracers get lost over time.Sediment burial or transportation beyond the study reach are possible causes.In this research, a total of 504 RFID-equipped pebbles were dispersed in the river at the beginning of the three selected reaches presented in Table 2. Table 3 presents the main features of each surveyed event and the tracers' recovery ratio.Namely, 327 pebbles were placed in reach 6, 118 in reach 7, and 69 in reach 11.This uneven distribution is linked to an expected recovery ratio of each stretch, hence more pebbles have been inserted in more active reaches to compensate for a lower expected recovery ratio.Deployment spots in reach 6 are shown in Figure 4.The RFID-equipped pebbles were selected to represent the D90 of each stretch of the river, and were collected directly in the river to characterize the bed transport that affects the largest fraction of the granulometric curve.A survey conducted after each event allowed for building a database in which the displacement of each pebble is recorded with the features of the event which caused the displacement (i.e., cumulative rainfall depth, rainfall duration, and peak discharge).Since single-event displacements were often the same order of magnitude as GPS accuracy, each pebble position was also recorded on a photogrammetric map to decrease uncertainty.A total of 1948 single displacements were recorded during the research period.Table 3 reports mean travel distances for each event and the recovery ratio (present tracers/found tracers).Significant increases in recovery ratio are linked to new tracer deployment campaigns.A weak proportional relationship can be found between peak discharge, travel distance, and recovery ratio, as one would expect (Figure 5).All the displacements recorded in  A survey conducted after each event allowed for building a database in which the displacement of each pebble is recorded with the features of the event which caused the displacement (i.e., cumulative rainfall depth, rainfall duration, and peak discharge).Since single-event displacements were often the same order of magnitude as GPS accuracy, each pebble position was also recorded on a photogrammetric map to decrease uncertainty.A total of 1948 single displacements were recorded during the research period.Table 3 reports mean travel distances for each event and the recovery ratio (present tracers/found tracers).Significant increases in recovery ratio are linked to new tracer deployment campaigns.A weak proportional relationship can be found between peak discharge, travel distance, and recovery ratio, as one would expect (Figure 5).All the displacements recorded in 2016 in reach 6 were inserted in a geographical information system (GIS) layered over a georeferenced 3D model of the reaches.Pathways are clearly indicated event-by-event and correspond, in most cases, to the thalweg (Figure 6).2016 in reach 6 were inserted in a geographical information system (GIS) layered over a georeferenced 3D model of the reaches.Pathways are clearly indicated event-by-event and correspond, in most cases, to the thalweg (Figure 6).

Eulerian Method
This method relies on a semi-Eulerian approach, based on the work of Scheinegross et al. [44].In determining the area of investigation at each reach, square grids (50 cm × 50 cm) were selected in the dry streambed; specifically, focusing on banks and bars that are not normally flooded to control their activity in case of intense events.Paint was then sprayed on these areas, and pictures were taken 2016 in reach 6 were inserted in a geographical information system (GIS) layered over a georeferenced 3D model of the reaches.Pathways are clearly indicated event-by-event and correspond, in most cases, to the thalweg (Figure 6).

Eulerian Method
This method relies on a semi-Eulerian approach, based on the work of Scheinegross et al. [44].In determining the area of investigation at each reach, square grids (50 cm × 50 cm) were selected in the dry streambed; specifically, focusing on banks and bars that are not normally flooded to control their activity in case of intense events.Paint was then sprayed on these areas, and pictures were taken

Eulerian Method
This method relies on a semi-Eulerian approach, based on the work of Scheinegross et al. [44].In determining the area of investigation at each reach, square grids (50 cm × 50 cm) were selected in the dry streambed; specifically, focusing on banks and bars that are not normally flooded to control their activity in case of intense events.Paint was then sprayed on these areas, and pictures were taken before and after the occurrence of each major rainfall event.The comparison and the analysis of the images allows one to highlight the colored sediments mobilized due to the flood, and those that remained immobilized, as shown in Figure 7. before and after the occurrence of each major rainfall event.The comparison and the analysis of the images allows one to highlight the colored sediments mobilized due to the flood, and those that remained immobilized, as shown in Figure 7. Image analysis was performed using the Basegrain [50] software, which assessed the number and the dimensions (a and b-axis) of the colored pebbles in the selected area and, thus, the granulometry of the area, according to Equation (1) [50], and the number of clasts for each granulometric class.

√ •
(1 where Dn is the nominal diameter, and a and b are, respectively, the maximum and minimum axis of the pebble.In this way, initial granulometry of the area is determined.The same spot is photographed after the flood for comparison purposes.When examining post-flood pictures, some of the painted sediment within the selected area could possibly be removed.The same approach is used to determine the granulometry of remaining sediments.Comparing the latter with the pre-flood granulometry, it can be easily understood which class has lost all or most of the tracers and, therefore, evaluate the critical diameter.Eighteen test areas were set up for the study in different times (Table 4).Sixteen were set up in reach 6, because it was the most suitable due to its almost completely natural bed, while two were set up in reach 7.
This method can account for the remaining painted stones, but does not provide information on entrained pebbles.Disappeared painted pebbles could, in principle, be either entrained or covered by aggradation.In this case, the incipient motion condition could present a bias from the right value.
The incipient motion condition, as defined by Shields [45], is used to describe sediment motion.Incipient motion is defined as the point at which stabilizing forces of grains are equaled by flow thrust.Shield's equation governing the process focuses on the critical shear stress calculation, which depends on sediment and flow features.
We compared measured values with numerical modeling to test the reliability of the method and the accuracy of the numerical estimation of the incipient motion condition.HEC-RAS [56] software was used to model the hydraulics of the river and to estimate the critical diameter of the pebbles.A 1D steady flow model has been built for this purpose, covering reaches 6 and 7.The model is 1 km long and contains just 11 sections, which were measured and quite uniformly distributed.To improve model stability, additional sections were interpolated at a fixed 50 m step.The model was constrained from the valley side, where a high check dam offers a border condition (Figure 3c).Measured discharge values for each event were input into the model.We determined a critical diameter recording flow rate of each event and used it as input for these models.Steady flow analysis with different Manning's n friction coefficient values were performed to understand the sensitivity Image analysis was performed using the Basegrain [50] software, which assessed the number and the dimensions (a and b-axis) of the colored pebbles in the selected area and, thus, the granulometry of the area, according to Equation (1) [50], and the number of clasts for each granulometric class.
where D n is the nominal diameter, and a and b are, respectively, the maximum and minimum axis of the pebble.In this way, initial granulometry of the area is determined.The same spot is photographed after the flood for comparison purposes.When examining post-flood pictures, some of the painted sediment within the selected area could possibly be removed.The same approach is used to determine the granulometry of remaining sediments.Comparing the latter with the pre-flood granulometry, it can be easily understood which class has lost all or most of the tracers and, therefore, evaluate the critical diameter.Eighteen test areas were set up for the study in different times (Table 4).Sixteen were set up in reach 6, because it was the most suitable due to its almost completely natural bed, while two were set up in reach 7.
This method can account for the remaining painted stones, but does not provide information on entrained pebbles.Disappeared painted pebbles could, in principle, be either entrained or covered by aggradation.In this case, the incipient motion condition could present a bias from the right value.
The incipient motion condition, as defined by Shields [45], is used to describe sediment motion.Incipient motion is defined as the point at which stabilizing forces of grains are equaled by flow thrust.Shield's equation governing the process focuses on the critical shear stress calculation, which depends on sediment and flow features.
We compared measured values with numerical modeling to test the reliability of the method and the accuracy of the numerical estimation of the incipient motion condition.HEC-RAS [56] software was used to model the hydraulics of the river and to estimate the critical diameter of the pebbles.A 1D steady flow model has been built for this purpose, covering reaches 6 and 7.The model is 1 km long and contains just 11 sections, which were measured and quite uniformly distributed.To improve model stability, additional sections were interpolated at a fixed 50 m step.The model was constrained from the valley side, where a high check dam offers a border condition (Figure 3c).Measured discharge values for each event were input into the model.We determined a critical diameter recording flow rate of each event and used it as input for these models.Steady flow analysis with different Manning's n friction coefficient values were performed to understand the sensitivity of the model to this parameter.In normal simulations, the Manning's n resistance coefficient in such irregular rivers can only be assigned by an educated guess, and any validation requires measurements of the physical parameters of the stream.Values obtained in this study were evaluated against established standards [57,58].Table 4 provides a comparison of measured and estimated critical diameter class for different flow events and locations.The methods applied have focused on short-term, micro-scale descriptions of sediment dynamics.Our objective is to characterize the river in different time and spatial domains, however, therefore in order to enlarge the focus to macro-scale, long-term processes, an alternative technique is sought.The technique is tested on reach 6, which is characterized by boulder-cascade morphology.It shows aggregated structures, with water flowing over and around boulders and step-pools which usually have lengths similar to the width of the channel.
The procedure adopted in this study consists of the comparison between two different photogrammetric models taken in two different periods: Spring, 2016 and autumn, 2017.These periods are more than 18 months apart, with several intense-discharge events experienced between them.
AGIsoft Photoscan 1.4.2(Agisoft LLC, St. Petersburg, Russia) was used to create dense point cloud models of the river stretch.Dense point clouds are groups of points within a 3D reference system which describe the surface of an object.Each epoch has a different point cloud composed of about 1,500,000 points.The computed accuracy of point clouds is <5 cm, based on measured ground control points (GCP).Once the point clouds have been created, they are processed using CloudCompare 2.6 open source software.This permits the overlapping of dense point clouds in order to calculate relative differences between them and highlight deposition or erosion spots.
The comparison highlighted the presence of four active deposition zones where material accumulated during the year.Figure 8

Discussion
A variety of methods were successful in recording different aspects of sediment dynamics in the Caldone River.However, no single method adequately explains the complexity of sediment transport in mountain streams.A multi-time and space domain approach is necessary for a full description of such complex processes.
In this work, meso-scale, short-term Lagrangian methods are applied to a study case, along with short-term and long-term, micro and macro-Eulerian approaches.Mixing different methods allows for a complete record of active processes.These can be different at different temporal scales.Sediment movement processes can be different within the same river reach when short-term analyses are performed.The database compiled in this research is based on one survey for each rain event.However, a much larger information bank is needed from event scale surveys to allow for characterization of the activity status of the river directly, and not through hypotheses or long time span mean values related to movements and morphological changes due to diverse events.
When comparing data from Eulerian and Lagrangian short-term approaches (for critical diameter measurements, for example), we found significant disagreement between the two methods.In event eight, during a long duration, low magnitude rainfall event, the critical grain size diameter in reach six was measured by the Eulerian approach as 45-64 mm. Figure 9 displays the number of clasts before and after event eight as a percentage of the initial number of clasts.Pre-event measurements, depicted in blue, are obviously 100%; post event measurements, in orange, are variable relative to classes.In event six, for a short, intense summer storm, the critical diameter was estimated as Dc > 90.5 mm.This is evidently due to different discharge rates and water stages (i.e., 0.67 m 3 /s in event eight, and 5.88 m 3 /s in event six), proving that the painted patches are able to measure the critical diameter for a certain location during an event.The same kind of behavior is also exhibited in other events where the critical diameter is measured as a function of discharge.Minor variations are present, which are caused by differences in the positioning of painted patches.These variations are in accordance with numerical modeling.A correlation between discharge and mean travelled distance of tracers can be found.In general, higher flow rates cause larger displacement.This direct proportionality is non-existent when dealing with Dc that depends on the flow condition in the specific measured spot.The accumulation areas are three banks and one bar, and accumulated material was mostly fine.The remaining part of the clouds show no significant variation between the two surveys.A visual inspection of clouds which were demonstrated as the main structures showed that they remained unchanged, as expected, since there were no major events.However, some cyclic erosional or depositional events are common and involve the same areas, but the reach is global in an equilibrium condition.

Discussion
A variety of methods were successful in recording different aspects of sediment dynamics in the Caldone River.However, no single method adequately explains the complexity of sediment transport in mountain streams.A multi-time and space domain approach is necessary for a full description of such complex processes.
In this work, meso-scale, short-term Lagrangian methods are applied to a study case, along with short-term and long-term, micro and macro-Eulerian approaches.Mixing different methods allows for a complete record of active processes.These can be different at different temporal scales.Sediment movement processes can be different within the same river reach when short-term analyses are performed.The database compiled in this research is based on one survey for each rain event.However, a much larger information bank is needed from event scale surveys to allow for characterization of the activity status of the river directly, and not through hypotheses or long time span mean values related to movements and morphological changes due to diverse events.
When comparing data from Eulerian and Lagrangian short-term approaches (for critical diameter measurements, for example), we found significant disagreement between the two methods.In event eight, during a long duration, low magnitude rainfall event, the critical grain size diameter in reach six was measured by the Eulerian approach as 45-64 mm. Figure 9 displays the number of clasts before and after event eight as a percentage of the initial number of clasts.Pre-event measurements, depicted in blue, are obviously 100%; post event measurements, in orange, are variable relative to classes.In event six, for a short, intense summer storm, the critical diameter was estimated as D c > 90.5 mm.This is evidently due to different discharge rates and water stages (i.e., 0.67 m 3 /s in event eight, and 5.88 m 3 /s in event six), proving that the painted patches are able to measure the critical diameter for a certain location during an event.The same kind of behavior is also exhibited in other events where the critical diameter is measured as a function of discharge.Minor variations are present, which are caused by differences in the positioning of painted patches.These variations are in accordance with numerical modeling.A correlation between discharge and mean travelled distance of tracers can be found.In general, higher flow rates cause larger displacement.This direct proportionality is non-existent when dealing with D c that depends on the flow condition in the specific measured spot.In comparing these findings with the Lagrangian approach, some differences emerge.In event eight, almost half of the tracers with D > 64 mm moved, even if they were bigger in terms of critical diameter.Mean displacement was 3.4 m.Conversely, during event six, most of the tracers showed a significant displacement, with 16.6 m as the mean displacement.A check of the diameters of moved and still tracers highlighted a critical diameter between 90 and 100 mm.In this case, the Eulerian and Lagrangian approaches agree.This behavior is not completely confirmed by other events.Little tracers tend to move more than bigger ones, as expected, but a sharp difference is not visible.Reduction in tracer population results in the sample to losing significance over time.
The difference in behavior between low and high discharge events can be explained through an analysis of the activity state of the channel.As suggested by Liébault et al. [38], bars and banks act as storage areas for sediments, turning their behavior from accumulation to source, depending on flow rate.Thus, during moderate discharge events, water wanders around bars and remains in the thalweg.In this situation, the only active part is the thalweg, and tracers are better than painted patches at describing sediment dynamics.When water stage increases, the whole channel is active and the two methods described in this paper can offer a superior interpretation of sediment dynamics.
The definition of activity status can also be challenging in the time domain.Erosion and deposition patterns in space and time are interesting, and their study can lead to a better definition of the virtual velocity concept.Nowadays, virtual velocity is defined as the ratio between grain displacement and the total time elapsed between two successive surveys [59].We focused on reach six, the most active and interesting reach, and found a slow seasonal variation of the river through a DEM comparison technique.However, in the Caldone River case, no general trend was found.This means that the river is in equilibrium during the researched period, with some areas clearly subjected to alternate episodes of erosion and deposition.

Conclusions
The final aim of much research on sediment dynamics in a mountain river is to find a value for solid discharge at the basin outlet.As stated by Haschenburger and Church [59], solid discharge is In comparing these findings with the Lagrangian approach, some differences emerge.In event eight, almost half of the tracers with D > 64 mm moved, even if they were bigger in terms of critical diameter.Mean displacement was 3.4 m.Conversely, during event six, most of the tracers showed a significant displacement, with 16.6 m as the mean displacement.A check of the diameters of moved and still tracers highlighted a critical diameter between 90 and 100 mm.In this case, the Eulerian and Lagrangian approaches agree.This behavior is not completely confirmed by other events.Little tracers tend to move more than bigger ones, as expected, but a sharp difference is not visible.Reduction in tracer population results in the sample to losing significance over time.
The difference in behavior between low and high discharge events can be explained through an analysis of the activity state of the channel.As suggested by Liébault et al. [38], bars and banks act as storage areas for sediments, turning their behavior from accumulation to source, depending on flow rate.Thus, during moderate discharge events, water wanders around bars and remains in the thalweg.In this situation, the only active part is the thalweg, and tracers are better than painted patches at describing sediment dynamics.When water stage increases, the whole channel is active and the two methods described in this paper can offer a superior interpretation of sediment dynamics.
The definition of activity status can also be challenging in the time domain.Erosion and deposition patterns in space and time are interesting, and their study can lead to a better definition of the virtual velocity concept.Nowadays, virtual velocity is defined as the ratio between grain displacement and the total time elapsed between two successive surveys [59].We focused on reach six, the most active and interesting reach, and found a slow seasonal variation of the river through a DEM comparison technique.However, in the Caldone River case, no general trend was found.This means that the river is in equilibrium during the researched period, with some areas clearly subjected to alternate episodes of erosion and deposition.

Conclusions
The final aim of much research on sediment dynamics in a mountain river is to find a value for solid discharge at the basin outlet.As stated by Haschenburger and Church [59], solid discharge is evaluated as the product of virtual velocity, active width, and active layer depth.This simplification implicitly assumes a uniform state of activity for the whole channel, in time and space.However, this assumption is far from reality in complex mountain rivers.Thus, techniques able to distinguish differences inside the same reach are needed.Painted patches are useful to evaluate the activity status of bars and banks and, thus, the active width of the stream.The patches are positioned on the less active areas (Figure 4d), but are sensitive even to shallow water stages flowing over them.They cannot give information about the travelled distances, but data about this comes from tracers.The measurement of active layer thickness needs ad hoc tools and surveys.The use of scour chains, wisely distributed in the channel, could be used to investigate active layer depth.The use of a more detailed survey technique leads to evidence of how different activity states can coexist in the same stretch and how they can mix depending on the flow rate.The description of this variability and differences in the activity status can lead to a deeper understanding of sediment dynamics, and a better overall description than virtual velocity estimated from long-term, macro-scale data.

Figure 1 .
Figure 1.Flow chart of the present work.Main features of applied methods are highlighted along with the results obtained.

Figure 1 .
Figure 1.Flow chart of the present work.Main features of applied methods are highlighted along with the results obtained.

Figure 2 .
Figure 2. Caldone River watershed and main interest points.Tracers were added in segments 6, 7, and 11.A rain gauge is present in segment 11. Level gauges are present in segments 7 and 11.Part of segment 11 flows in a culvert.Red crosses stand for the two sedimentation pools, while red ticks represent granulometric analysis spots.

Figure 3 .
Figure 3. Key spots in the Caldone River: (a) Sedimentation pond in reach 6; (b) culvert entrance in reach 11; (c) level gauge in reach 7, notice the check dam used as hydraulic model constrain.

Figure 2 .Figure 2 .
Figure 2. Caldone River watershed and main interest points.Tracers were added in segments 6, 7, and 11.A rain gauge is present in segment 11. Level gauges are present in segments 7 and 11.Part of segment 11 flows in a culvert.Red crosses stand for the two sedimentation pools, while red ticks represent granulometric analysis spots.

Figure 3 .
Figure 3. Key spots in the Caldone River: (a) Sedimentation pond in reach 6; (b) culvert entrance in reach 11; (c) level gauge in reach 7, notice the check dam used as hydraulic model constrain.Figure 3. Key spots in the Caldone River: (a) Sedimentation pond in reach 6; (b) culvert entrance in reach 11; (c) level gauge in reach 7, notice the check dam used as hydraulic model constrain.

Figure 3 .
Figure 3. Key spots in the Caldone River: (a) Sedimentation pond in reach 6; (b) culvert entrance in reach 11; (c) level gauge in reach 7, notice the check dam used as hydraulic model constrain.Figure 3. Key spots in the Caldone River: (a) Sedimentation pond in reach 6; (b) culvert entrance in reach 11; (c) level gauge in reach 7, notice the check dam used as hydraulic model constrain.

Figure 4 .
Figure 4. Deployment of the tracers and their position in the stream, reach 6 (modified after [55]).(ac) Show three deployment spots in reality and on the ortophoto; (d) Is the ortophoto of reach 6 where yellow dots are tracer starting points, numbered from 1 to 4, red dots are painted bed patches, and red squares are the most active zones found with digital elevation model (DEM) comparison.

Figure 4 .
Figure 4. Deployment of the tracers and their position in the stream, reach 6 (modified after [55]).(a-c) Show three deployment spots in reality and on the ortophoto; (d) Is the ortophoto of reach 6 where yellow dots are tracer starting points, numbered from 1 to 4, red dots are painted bed patches, and red squares are the most active zones found with digital elevation model (DEM) comparison.

Figure 5 .
Figure 5. Covered distance (blue crosses) and recovery ratio (orange squares) plotted versus measured peak discharge for each event.It is possible to notice a weak direct relationship between distance and discharge and inverse between recovery ratio and discharge.Recovery ratio is also affected by new tracer insertion.

Figure 6 .
Figure 6.Positions of tracers in 2016 events for reach 6, where each color represent a date.Notice how the pebbles move downstream in time and remain in the thalweg of the river.

Figure 5 .
Figure 5. Covered distance (blue crosses) and recovery ratio (orange squares) plotted versus measured peak discharge for each event.It is possible to notice a weak direct relationship between distance and discharge and inverse between recovery ratio and discharge.Recovery ratio is also affected by new tracer insertion.

Figure 5 .
Figure 5. Covered distance (blue crosses) and recovery ratio (orange squares) plotted versus measured peak discharge for each event.It is possible to notice a weak direct relationship between distance and discharge and inverse between recovery ratio and discharge.Recovery ratio is also affected by new tracer insertion.

Figure 6 .
Figure 6.Positions of tracers in 2016 events for reach 6, where each color represent a date.Notice how the pebbles move downstream in time and remain in the thalweg of the river.

Figure 6 .
Figure 6.Positions of tracers in 2016 events for reach 6, where each color represent a date.Notice how the pebbles move downstream in time and remain in the thalweg of the river.
reports a visual comparison (a and b) of the same spot in reaches 1 and 2. Figure 8c displays a point-to-point cloud comparison of another spot, showing clear accumulation.The image depicts a portion of the dense point cloud of autumn, 2017.The points have been colored based on their distance from the spring, 2016 point cloud.Blue means no difference and the points are close, while green to orange means a growing distance from 0 m to 0.2 m.The blue shadow is background.The accumulation is clearly shown by the green area in the picture.

Figure 8 .
Figure 8. Upstream deposition area.(a) Picture taken in spring, 2016 and in (b) in fall, 2017; (c) Cloudto-cloud comparison results; color range highlights the accumulation of debris (in green).

Figure 8 .
Figure 8. Upstream deposition area.(a) Picture taken in spring, 2016 and in (b) in fall, 2017; (c) Cloud-to-cloud comparison results; color range highlights the accumulation of debris (in green).

Figure 9 .
Figure 9. Remaining clasts, as a percentage of the original sample, divided per granulometric class after event eight.The blue bar is the proportion of clasts before event eight (100%), the orange bar is for the amount after (variable percentage).

Figure 9 .
Figure 9. Remaining clasts, as a percentage of the original sample, divided per granulometric class after event eight.The blue bar is the proportion of clasts before event eight (100%), the orange bar is for the amount after (variable percentage).

Table 1 .
Main instruments available in the Caldone River.Elevation expressed as Above sea level (a.s.l.).

Table 2 .
Main features of the three studied reaches.

Table 3 .
Main features of surveyed events.Peak discharge is measured at an outlet.

Table 3 .
Main features of surveyed events.Peak discharge is measured at an outlet.

Table 4 .
Critical diameter (D c ) comparison table.Flowrate was measured in the surveyed reach.