Expanding the Sediment Transport Tracking Possibilities in a River Basin through the Development of a Digital Platform—DNS/SWAT

: Simulation of stochastic and variable sediment transport processes within models still poses a big challenge, especially in mountainous areas. Since sediment transport, including erosion and deposition, remains an unceasing problem in many areas, sediment modeling is perceived as a possible solution. This article combines a review of the selected sediment models with a presentation of the effects of several years of research using the DNS digital platform in the Western Carpathians. The review focuses on the main advantages and gaps in selected modeling tools with particular emphasis on one of the most popular: SWAT. The description of the digital platform—DNS is an example of how to answer these gaps by combining subsequent models, methods, and databases using their best features. To accentuate the beneﬁts of such an approach, the effects of combining subsequent models (AdH/PTM) and methods (ﬁngerprinting) on a common digital DNS space are presented, on the example of the Raba River (basin). In this way, both unique possibilities of estimating the amount of contamination carried with sediment particles and their sources, as well as sequencing of sedimentation in the reservoir, taking into account its subsequent zones, were obtained.


Introduction
Sediment transport is one of the most important processes for the mass movement of soil particles within a basin area. The driving force behind this process is primarily surface runoff in the land phase and flows in the channel network in the riverbed phase [1,2]. In practice, this concept covers both transport and the erosion that initiates it [3], and the deposition that concludes it [4,5]. This process is of great importance due to its social, ecological, and economic effects [6]. It is also highly sensitive to changes in land use, climate, river regulation, operation of hydrotechnical devices, and dam reservoirs [4,5]. These are some of the main reasons why sediment transport is becoming more and more non-linear in both time and space. Their consequence is a decrease in the depth, fertility and productivity of soils, a decrease in the river's capacity, a shortened life of dam reservoirs, and thus, an increase in the risk of flooding [7]. Not surprisingly, in recent decades, this process has become an important topic on the agenda of local, national, and European policy makers, and is of scientific interest [8][9][10][11][12]. Sediment particle tracking is now one of the priority problems for water management, and environmental models were quickly identified as a possible solution [5,[13][14][15][16]. This task was seemingly simple, but it quickly became clear that river basins are dynamic and complex systems, and the processes responsible for the transport of sediments in each of the phases are stochastic and variable [17]. In mountainous areas, the problem becomes even more complicated. Observational periods are usually short [17][18][19], and the system of channels and rivers is often cut off by antierosion structures and dam reservoirs, which retain almost all, especially coarse-grained The aim of this article is to present the current state of knowledge in the field of interdisciplinary studies of sediment transport in the basin area, with a particular emphasis on the mountain areas of the Western Carpathians. The manuscript presents the possibilities of using the developed digital platform-DNS (Discharge-Nutrient-Sea) with a SWAT module, to track soil particles from the source to the deposition in the pilot basin of the Raba River located in the Polish part of the Western Carpathians. This basin, due to its dual character and presence of a dammed reservoir, highlighted both the strengths and weaknesses of SWAT. This allowed to evaluate the possibilities and effectiveness of combining models, methods, and databases as further modules of the digital platform, expanding the possibilities and reducing limitations in tracing sediments both in the land and riverbed phase.

Digital Platform-DNS with SWAT Module
The digital platform-DNS was developed at the Polish Institute of Meteorology and Water Management-National Research Institute (IMGW-PIB). The decisive factor in its development was the ability to use the best features of models, methods, and databases concentrated in one digital space ( Figure 1). Initially, the digital platform-DNS was used to simulate the transport of nutrients in basins with limited hydrological and qualitative monitoring data. This was made possible by combining the SWAT model and e.g., MIKE 11 [63], as well as IMWM-PIB monitoring databases and data-processing modules for supplementing missing measurement data, which was described in detail in [64][65][66][67]. SWAT has quickly become the most important module of the digital platform-DNS. It is a physical, semi-distributed continuous-time model that runs on a daily, monthly, or yearly time step. Its basic feature is the possibility of continuous simulation of hydrological processes (evapotranspiration, surface runoff, percolation, return flow, groundwater flow, transmission losses in the canal, retention in ponds and reservoirs, channel routing, and field drainage). However, SWAT can also simulate the transport of soil particles, both in the land and riverbed phases [68,69]. This feature quickly proved to be useful when trying to track sediment particles in the mountain basin in the Polish part of the Western Carpathians [70][71][72][73]. The key to a reliable simulation of sediment transport in SWAT is the appropriate quantity and quality of parameters and input data. Their task is to best represent the basin area with its natural and anthropogenic processes in the digital model space. The advantage and disadvantage of this tool is that it requires both the knowledge of a large number of parameters and the collection of appropriate input data [38,74].
SWAT was developed in the USA and many of its default parameters must first be adapted to the local conditions of the analyzed basin and be kept within a realistic range of uncertainty. In turn, the input data should describe the basin's topography, soil characteristics and distribution, land cover, and meteorological conditions in the basin area at an appropriate level of detail. The basic input data required to build this tool are: − A digital elevation model (DEM) for information on the watercourse stream network and length, drainage pattern of the watershed, channel width within the watershed, slope, and reach length; − A map of hydrographical divisions in order to divide the basin area into sub-basins; − land use maps and agrotechnical data for determining land cover classes; − A soil map-detailed data on soil types to identify soil classes and their physical properties such as available water content, soil texture, soil bulk density, hydraulic conductivity, and organic matter; − Meteorological data for precipitation, temperature, humidity, wind speed and direction, and solar radiation.

Land Phase
This phase is the main source of sediment in a majority of basins, which is directly related to rainfall and runoff energy [75,76], which have a decisive influence on the phenomena of detachment and initiation of soil particle transport. The two most important processes that occur in this phase are surface runoff and erosion (detailed in Chapters 2.1.1 and 2.1.2), with soil type and land use being the most important responders. The SWAT module takes these responses into account by dividing the basin into areas with unique land use, soil, and slopes, creating so-called hydrological response units (HRU). For each HRU, the individual hydrological components of the land phase, such as rainfall, surface runoff, and evapotranspiration, are then simulated. The basis here is the water balance Equation (1), which was described in detail in the theoretical SWAT documentation [77]. Based on the daily components of the hydrological cycle the daily water budget in each HRU is calculated [37,[78][79][80][81].
where: SWt is the final soil water content, SW0 is the initial soil water content, Rday is the amount of precipitation, Qsurf is amount of surface runoff, Ea is the amount of evapotranspiration, Wseep is the amount of water entering the vadose zone from the soil profile, and Qgw is the amount of return flow (all in mm H2O). The SWAT module has been successfully used many times around the world to simulate sediment yield transport in the land phase of the basin area, as confirmed by numerous publications (e.g., [82][83][84]). Nevertheless, there are some limitations to both the data and the module itself. These include the common problem with the unavailability of data on sediment yield efficiency [85]. In practice, this causes situations where the time step of sediment calibration and validation is different than for hydrology, for which monitoring data are usually more readily available even for smaller basins. This reduces the accuracy of the model in the simulation of the land-phase sediment yield. The problem is also the generalization of the SWAT model, the effect of which is, for example, allowing all eroded

Land Phase
This phase is the main source of sediment in a majority of basins, which is directly related to rainfall and runoff energy [75,76], which have a decisive influence on the phenomena of detachment and initiation of soil particle transport. The two most important processes that occur in this phase are surface runoff and erosion (detailed in Chapters 2.1.1 and 2.1.2), with soil type and land use being the most important responders. The SWAT module takes these responses into account by dividing the basin into areas with unique land use, soil, and slopes, creating so-called hydrological response units (HRU). For each HRU, the individual hydrological components of the land phase, such as rainfall, surface runoff, and evapotranspiration, are then simulated. The basis here is the water balance Equation (1), which was described in detail in the theoretical SWAT documentation [77]. Based on the daily components of the hydrological cycle the daily water budget in each HRU is calculated [37,[78][79][80][81].
where: SW t is the final soil water content, SW 0 is the initial soil water content, R day is the amount of precipitation, Q surf is amount of surface runoff, E a is the amount of evapotranspiration, W seep is the amount of water entering the vadose zone from the soil profile, and Q gw is the amount of return flow (all in mm H 2 O). The SWAT module has been successfully used many times around the world to simulate sediment yield transport in the land phase of the basin area, as confirmed by numerous publications (e.g., [82][83][84]). Nevertheless, there are some limitations to both the data and the module itself. These include the common problem with the unavailability of data on sediment yield efficiency [85]. In practice, this causes situations where the time step of sediment calibration and validation is different than for hydrology, for which monitoring data are usually more readily available even for smaller basins. This reduces the accuracy of the model in the simulation of the land-phase sediment yield. The problem is also the generalization of the SWAT model, the effect of which is, for example, allowing all eroded soil to reach the riverbed with surface runoff (without taking into account sedimentation still in the land phase), as well as simulation of sediment-bound pollutants, determined by the primary concentration of nutrients in the soil, soil loss, and the soil coefficient enrichment, which is often not applicable in many basins [86]. Another limitation of the model in this phase is the use of numerous empirical and quasi-physical equations developed on the basis of climatic conditions in the USA [87], including simplified sediment routing algorithms for erosion simulation [84,88,89].

Surface Runoff
Accurate determination of surface runoff is critical in estimating sediment transport in the land phase of the basin [90,91]. The SWAT module helps to estimate surface runoff, taking into account the related transport of organic and inorganic pollutants [92]. In this tool, it is possible to choose one of two methods of estimating this phenomenon: the Green and Ampt infiltration method [93,94], and the Soil Conservation Service Curve Number (SCS-CN) procedure [95]. The physical method of Green and Ampt assumes a sharp wetting front during infiltration with saturated soil above the unsaturated zone [96]. Nevertheless, the popularity of this method is much lower due to the required high time resolution of the input data on precipitation (time step less than 24 h) [97,98]. The SCS-CN empirical model (2 and 3), which is mainly based on soil properties, land use, and hydrological conditions, is widely used for surface runoff estimation with daily precipitation as input.
where: Q surf is the daily surface runoff (mm), Rday is the rainfall depth for the day (mm), and S is the retention parameter (mm). The retention parameter S and the prediction of lateral flow by SWAT model are defined in Equation (3): where: S is the drainable volume of soil water per unit area of saturated thickness (mm/day), and CN is the curve number. However, confidence in SCS-CN results should be limited [99,100] because this method causes problems with underestimating peak flows in relation to monitoring data (e.g., [84,101]). One of the causes is that the constant initial abstraction is used as an intrinsic parameter, which may result in inconsistency of the soil moisture accounting. In addition, the duration and intensity of rainfall are not considered; instead, the average daily rainfall is used as a SWAT input. In fact, high-intensity and even short-duration rainfall could generate more sediment than actualized in the model based on daily rainfall [86]. Therefore, after selecting this method, it is important to perform, in addition to the statistical evaluation of the model, a visual evaluation of the results, allowing a more intuitive way to characterize the model's ability to recreate the characteristics of the flow regime for a given basin area, especially floods and low flows. The few available studies [97] generally suggest that the Green and Ampt method works better for stream flow and avoids the limitations of the SCS method. The requirements for high-resolution precipitation data can be met by precipitation atlases containing accurate data from many years with a step for even every 10 min, which has already been prepared and described for the all territory of Poland by [102,103].

Erosion
Erosion is a key process that initiates the transport of sediment particles down the basin. At the same time, it is one of the most serious problems of environmental degradation, adversely affecting many natural and man-made ecosystems [104]. In agricultural basins, Appl. Sci. 2022, 12, 3848 6 of 20 soil erosion not only reduces the nutrient-rich topsoil in place, but also degrades water quality by moving sediment off site. An example of a compromise between empirical and physical algorithms in the SWAT module is the use of a modified Universal Soil Loss Equation (MUSLE) to predict soil erosion at the HRU level resulting in soil particle transport. MUSLE (4) is a modified version of USLE [105] based on the runoff characteristics as the best single indicator for predicting sediment performance at the basin outlet, and factors influencing soil erosion. By using the energy of surface runoff, not rainfall, MUSLE is suitable for use on a daily basis to estimate the sediment efficiency, and takes into account sediment deposition in HRU [106]. MUSLE eliminates the need for delivery ratios and allows the equation to be applied to individual storm events, and therefore improves prediction of sediment loading [84,[107][108][109].
where: SY = HRU sediment yield (t/day); Q = daily runoff volume (mm); qp = runoff peak discharge (m 3 /s); A = HRU area (ha); C, P, K, and LS are dimensionless factors accounting for HRU crop cover, soil protection, soil erodibility, and topography as defined in the original Universal Soil Loss Equation (USLE); and F CRFG is a dimensionless factor to account for coarse fragment cover (stoniness).
Although the effectiveness of MUSLE in sediment transport analysis has been confirmed [109,110], it should be remembered that the quality and reliability of sediment simulations are not only directly related to runoff simulations, but also show high sensitivity to rainfall characteristics. The kinetic energy of the raindrops plays an important role in soil detachment, but also causes runoff disturbance. Ignoring them by the SWAT module, in the case of some basins, may lead to a deterioration in the adjustment of simulated data to real conditions [111].

Riverbed Phase
This phase controls the movement of water and suspended sediment generation through the channel network of the watershed to the outlet [81,112,113]. The water is routed through the main canal network using the Muskingum routing method using Manning's Equation (5), and the suspended sediment transported with the water is subject to deposition and degradation processes (6-9) [114]. These two processes are computed with a simplified version of the Bagnold stream power equation, where the maximum amount of sediment transported from each sub-basin is a function of the peak channel velocity. If the concentration of sediment in the basin area, at the beginning of the time step, exceeds the maximum concentration of sediment that can be transferred by the water, additional sediment is deposited; if not, degradation occurs, which is a function of the channel erosion and the channel vegetation cover [115].
The average flow velocity, v ch (m/s), is calculated using Manning's equation: where: sl p ch is the riverbed phase of each sub-basin slope (m/m), and n is the Manning coefficient.
The maximum concentration of sediment that can be transported through in the riverbed phase of each sub-basin is closely related to the peak channel velocity at the daily time step. The peak channel velocity (m/s) and the maximum sediment concentration are calculated as: v ch,pk = PRF * v ch (6) conc sed,ch,mx = C sp * v spexp ch,pk where: PRF is the peak rate adjustment factor, v ch is the mean flow velocity, C sp is a coefficient, and spexp is an exponent. The conc sed,ch,mx (t/m 3 ) is compared to the actual concentration of sediment in the riverbed phase of each sub-basin (t/m 3 ). If conc sed,ch,i > conc sed,ch,mx , deposition dominates in the riverbed phase of each sub-basin, sediment routing and the sediment deposition sed dep is calculated as: However, if conc sed,ch,mx > conc sed,ch,i , degradation or erosion dominates in the riverbed phase of each sub-basin, sediment routing and the sediment degradation sed deg is calculated as: where: K ch is the erodibility factor in the riverbed phase of each sub-basin, and C ch is the cover factor in the riverbed phase of each sub-basin. However, the Bagnold equation, the default in the SWAT module, is considered impractical to apply across the entire streaming network [116]. The SWAT module also offers three alternative sediment tracking methods to the Bagnold method (Kodoatie, Molinas and Wu, and Yang Sand and Gravel). All these equations share the same in-stream sediment routing (sediment transport capacity is routed in main reaches/channels), but calculate the maximum sediment transport capacity (maximum concentration of transportable sediment) differently [117], and allow for the analysis of individual sediment fractions. The Yang equation was developed to simulate large particles such as sand and gravel. The Molinas and Wu equation and the Kodoatie equation are better suited for simulating small particles forming suspended sediments. In all four equations, bed erosion exceeded bank erosion [118][119][120]. Despite this choice of equations for sediment tracking in the riverbed phase, studies show that SWAT does not have a solid representation of the sediment processes in the stream [121]. The mechanisms associated with the transport of cohesive and non-cohesive sediment fractions in stream systems are different. The mentioned empirical and physical methods implemented in SWAT for the erosion, deposition, and transport of sediment in a stream system are universally applied across all sediment fractions. Consequently, a careful and critical approach to the obtained results must be exercised, especially in the case of the basin area in which the riverbed has been identified as the main source of sediment [44,122].

Using the Digital Platform-DNS/SWAT
In order to track suspended sediment transport in the land and riverbed phase of the pilot basin, the most important SWAT module was used on the digital platform-DNS, which was supported by data processing (LOAD ESTimator), and database modules (meteorological and land use). The use of the United States Geological Survey (USGS) FORTRAN-LOAD ESTimator program [123] made it possible to take into account the sediment concentration in, i.e., extreme flow situations, using the available interpolated sediment observations to fill the time gaps. LOADEST provides an estimate of the sediment load as a function of the observed daily river flow. This function is determined by performing a linear regression analysis between the available sediment concentration observations and the corresponding river flow data. In turn, the database module (meteorological and land use) made it possible to take into account the data on the current and future average monthly temperatures, and precipitation [124], as well as the results of the FORECOM project for the entire area of the Polish part of the Western Carpathians [125], taking into account the forecast changes in land use. This allowed for the creation of variant scenarios and the study of the effects of changing one, or several parameters over time.
The Raba basin area, which is over 131 km long, is part of the Polish Western Carpathians. Its sources are located in the Gorce mountain range at an altitude of 780 m a.s.l., and the estuary to the Vistula River at a height of 180 m a.s.l. [126,127]. The choice of this basin for research with the use of the digital platform-DNS/SWAT was related, i.e., by its dual nature and location ( Figure 2). The mountainous nature of its upper part manifests itself in a quick reaction to rainfall and high dynamics of surface runoff. On the other hand, in the lower part, as the height of the slopes decrease, the amount of precipitation and the dynamics of surface runoff also decrease [128]. It is also an area dominated by agriculture, which further contributes to soil loss in the area. The Raba riverbed phase is complex. The main bed runs through the bottom of an alluvial valley and is characterized by a relatively steep slope, a significant inflow of sediments from the slopes, and frequent floods [128][129][130]. Moreover, features of the basin area, such as its location in an area particularly exposed to water erosion, the proximity to one of the largest Polish urban agglomerations-Krakow, and the presence of a dam reservoir halfway along the river, make both quantitative and qualitative tracking of sediments a challenge. The Raba basin area, which is over 131 km long, is part of the Polish Western Carpathians. Its sources are located in the Gorce mountain range at an altitude of 780 m a.s.l., and the estuary to the Vistula River at a height of 180 m a.s.l. [126,127]. The choice of this basin for research with the use of the digital platform-DNS/SWAT was related, i.e., by its dual nature and location ( Figure 2). The mountainous nature of its upper part manifests itself in a quick reaction to rainfall and high dynamics of surface runoff. On the other hand, in the lower part, as the height of the slopes decrease, the amount of precipitation and the dynamics of surface runoff also decrease [128]. It is also an area dominated by agriculture, which further contributes to soil loss in the area. The Raba riverbed phase is complex. The main bed runs through the bottom of an alluvial valley and is characterized by a relatively steep slope, a significant inflow of sediments from the slopes, and frequent floods [128][129][130]. Moreover, features of the basin area, such as its location in an area particularly exposed to water erosion, the proximity to one of the largest Polish urban agglomerations-Krakow, and the presence of a dam reservoir halfway along the river, make both quantitative and qualitative tracking of sediments a challenge.  The model calibration (Myślenice and Proszówki) and validation process (Stradomka) and the preceding sensitivity analysis were carried out using the SWAT-CUP program [131], and the SUFI-2 algorithm [39,132]. The SUFI-2 algorithm tries to capture as many optimal simulations as possible that are within the 95% prediction certainty. The algorithm does not obtain a unique value for the parameters, but an interval that includes all the uncertainties of the processes in the basin. The sensitivity analysis allows us to assess the reaction of the model to the change of input parameters, allowing us to identify the most important parameters of the model for the analyzed basin. The sensitivity of parameters is measured using t-stat values, where the values are more sensitive to larger absolute values. The pvalues are used to determine the significance of the sensitivity when the parameter becomes significant if the p-values are close to zero [132].
Most of the parameters selected by the sensitivity analysis are the same for both the upper and lower parts of the Raba basin, but their significance for sediment transport simulation differs (Table 1). In mountainous areas, surface runoff (SURLAG), followed by soil erosion (USLE_K) and the value of saturated conductivity (SOL_K), which is a parameter controlling surface runoff, have the greatest impact on sediment transport simulations. In addition, specific only to the upper Raba are parameters such as the soil evaporation compensation factor (ESCO), and available soil water capacity (SOL_AWC) responsible for the so-called surface response. In contrast, in the sub-mountainous area, groundwater processes (GW_DELAY), and the ratio of soil loss within a specific support practice (e.g., contour tillage, strip tillage) to the corresponding up-and-down cultivation loss (USLE_P) have the greatest impact. The surface runoff (SURLAG) came in third. Characteristic only for this part of the basin area are the parameters responsible for the socalled subsurface response (GW_DELAY, ALPHA_BF), and reservoir parameters (RES_SED, RES_RR, and RES_NSED), which shows the important influence of a dam reservoir on the hydrology of the basin and, consequently, on sediment transport. Model calibration for the flow and sediment was carried out in the calculation profiles of Myślenice (upper part) and Proszówki (lower part), and validation for the Stradomka river, which is the right-bank tributary of the Raba. The low frequency of sediment concentration monitoring (12 times a year) was the reason why LOAD ESTimator was used to develop a regression model to reliably estimate loads. Four statistical measures were used to assess the correctness of this process: coefficient of determination (R 2 ) [133], the efficiency of the Nash-Sutcliffe (NSE) [134], the percentage of bias (PBIAS) [15], and Kling's efficiency-Gupta (KGE) [134]. Using the ranges of values for these statistical measures (Table S1), the degree of matching of the simulation results to the observed data was assessed (Table S2). On this basis, the model was considered calibrated and validated at a level that allows the use of this tool to track sediment particles in the basin area. The entire procedure is described in detail in [135].

Sediment Yield-Land Phase Simulation
Transport of sediment in the land phase of the basin in the SWAT module consists of simulating the mass of the sediment leaving individual area units (HRU) using the MUSLE (2) equation per time and per unit area, known as sediment yield (SYLD) [135]. HRUs allow the spatial heterogeneity of land use, soil, and slope to be taken into account, increasing the accuracy of SYLD forecasting, which is of particular importance in distinctly dual character basins such as the Raba. The basin area is often considered the most important explanatory factor for SYLD, due to the reduction in the rate of topsoil erosion on smaller slopes and the increased likelihood of sediment deposition as the basin area increases. In addition, especially in mountainous areas, SYLD is influenced by factors such as topography, land cover, reservoir effects, and climatic conditions. However, their relative importance in explaining the spatial variation of SYLD is still not fully understood, as it is specific to each basin area [136]. Analyses of regional SYLD patterns mainly concern larger river systems [137,138], and this further limits information on smaller basins that are often more important for sediment transport. The Western Carpathians, including the Raba River basin, can be called "hot spots" where, as a result of erosion, the transport of significant amounts of SYLD is initiated. Nevertheless, sediment yield monitoring in this area is very limited and was carried out for a limited time, only for individual tributaries of the Dobczyce dam reservoir [139]. Only the use of mathematical models has significantly increased the knowledge of the transport of sediment particles in the land phase throughout the Carpathian region [56][57][58]140,141]. The digital platform-DNS/SWAT made it possible to estimate the aggregated values of SYLD, being a function of the runoff processes, for each of the 36 Raba sub-basins, taking into account their temporal variability. Thus, spring was identified as the period of maximum average sediment yield values reaching 0.92 t/ha and 0.57 t/ha, and winter with minimum values of 0.28 t/ha and 0.21 t/ha for the upper and lower parts of the basin, respectively [138,142]. These studies can therefore both confirm the role of modeling as a reliable means of obtaining detailed information on the rate and amount of soil loss in areas such as the Carpathian region, and help understand the effect of basin factors such as slopes, land use, and soil type on the spatial variability of SYLD.

Suspended Sediment Load-Riverbed Phase
The analysis of the quantity and time variability of sediment load in the riverbed phase was a continuation of the tracking of soil particle transport in the land phase [71]. The SWAT module enables simulation of the sediment mass per time unit (sediment load), equal to the product of the flow rate and the concentration of sediment particles contained within it, on the calculation profile of the sub-basins (Figure 2). Contrary to SYLD, the concentration of sediment and the flow rate are among the basic parameters monitored by the SEM and the IMGW-PIB. This allowed the use of SWAT to estimate the sediment load flowing from both the upper (Dobczyce calculation profile) and the lower part of the Raba basin (Proszówki calculation profile) (Figure 2). The alternative Kodoatie method was selected for the calculations, which meant mainly tracking the fractions contained in the suspended sediment (<63 µm). Selected modifications to the riverbed, such as antierosion structures [143] and related disturbances in sediment transport, were also taken into account.
Ultimately, the application of the SWAT module made it possible to estimate the average load of suspended sediment flowing into the dam reservoir from the upper Raba area at 553 t/m [71]. It is worth noting that on the calculation profile of Proszówki, closing the lower part of the Raba basin (Figure 2), the suspended sediment load was estimated at an average level of 897 t/m. This is despite the dam reservoir that traps most sediment particles flowing into it from the upper Raba. However, this difference is justified and results both from the type of soil in the lower Raba, susceptible to erosion, and the dominant role of agriculture in this area, which, despite the lower slopes, favors soil loss.
The visual analysis of the results showed the negative effects of the SCS-CN method, forced by the daily time step of rainfall data for the Raba. As a result, despite the at least satisfactory results of the calibration and validation of the model (in Table S2), there are visible underestimations of the simulated values of the monitoring results, especially in the periods of peak flows (Figure 3), which also affects the matching of sediment loads.
phase was a continuation of the tracking of soil particle transport in the land phase [71]. The SWAT module enables simulation of the sediment mass per time unit (sediment load), equal to the product of the flow rate and the concentration of sediment particles contained within it, on the calculation profile of the sub-basins (Figure 2). Contrary to SYLD, the concentration of sediment and the flow rate are among the basic parameters monitored by the SEM and the IMGW-PIB. This allowed the use of SWAT to estimate the sediment load flowing from both the upper (Dobczyce calculation profile) and the lower part of the Raba basin (Proszówki calculation profile) (Figure 2). The alternative Kodoatie method was selected for the calculations, which meant mainly tracking the fractions contained in the suspended sediment (<63 μm). Selected modifications to the riverbed, such as antierosion structures [143] and related disturbances in sediment transport, were also taken into account.
Ultimately, the application of the SWAT module made it possible to estimate the average load of suspended sediment flowing into the dam reservoir from the upper Raba area at 553 t/m [71]. It is worth noting that on the calculation profile of Proszówki, closing the lower part of the Raba basin (Figure 2), the suspended sediment load was estimated at an average level of 897 t/m. This is despite the dam reservoir that traps most sediment particles flowing into it from the upper Raba. However, this difference is justified and results both from the type of soil in the lower Raba, susceptible to erosion, and the dominant role of agriculture in this area, which, despite the lower slopes, favors soil loss.
The visual analysis of the results showed the negative effects of the SCS-CN method, forced by the daily time step of rainfall data for the Raba. As a result, despite the at least satisfactory results of the calibration and validation of the model (in Table S2), there are visible underestimations of the simulated values of the monitoring results, especially in the periods of peak flows (Figure 3), which also affects the matching of sediment loads.

Digital Platform-DNS/SWAT-New Modules
The next stage of the research was an attempt to answer questions about the places of deposition of individual sediment fractions in selected zones of the dam reservoir, and the pollution transported with sediment particles and their sources. Obtaining answers to these questions, however, exceeded the capabilities of a single SWAT module. The open nature of the digital platform-DNS allowed to use its possibilities of combining models (AdH-Adaptive Hydraulic Model)/PTM-Particle Tracking Model), and methods (fingerprinting) as further modules supplementing the SWAT limitations. In the area of the Western Carpathians, the capabilities of these tools were rarely used or not used at all. Limitations also applied to monitoring data. Therefore, the database module was also expanded with the results of soil measurements-s in the basin area, and continuous measurements of the suspended sediment made by the BS-2 suspended sediment sampler permanently installed on the Wolnica River (Figure 2), which is a tributary of the dam reservoir [72].

AdH/PTM
The SWAT module made it possible to simulate suspended sediment load transport in the analyzed Raba River basin, maintaining the continuity of this process between the upper and lower part of the basin, by taking into account the dam reservoir. So, it was possible to estimate the suspended sediment load both before and after the reservoir.
However, this tool could not show the sequencing of sedimentation in the tank itself, divided into zones where this process occurs with varying intensity depending on the period and size of sediment particles [41]. Therefore, for this purpose, the AdH/PTM model was used, which was previously built and verified by the method of comparative analysis with real data [144,145]. AdH is a two-dimensional hydrodynamic modeling tool that uses the finite element method to solve two-dimensional equations of momentum conservation for water in the Euler system. On this basis, the particles' behavior over time (entrainment, advection, diffusion, settling, deposition, burying, etc.) was simulated using the PTM module. At each time step, PTM performs calculations to determine the local characteristics of the environment and the behavior of each representative particle tracked. Their combination with the SWAT module on the DNS platform allowed for precise tracking of particles of individual suspended sediment fractions from their source, through transport, and ultimately their sedimentation in each of the four zones of the dam reservoir with a monthly time step [146]. The input data for the AdH/PTP module were the SWAT simulation results for the monthly mean values of the suspended sediments obtained in the previous study (Section 3.2). The functionality of the Kodoatie method was used to determine the size of the load of individual suspended sediment fractions (mineral: CLAY-0-0.004 mm; SILT-0.004-0.062 mm, and mineral/organic: SMAG-0.03 mm) flowing to the Dobczyce dam reservoir (Myślenice calculation profile), amounting to 1.11 t/m, 96.7 t/m, and 455 t/m, respectively. These particle masses of the individual fractions were characterized by essentially different deposition sites in the reservoir. In the case of SMAG, the highest percentage of particles belonging to this fraction was deposited in the summer months (June, July, and August) in the backwater zone and the Myślenice basin (from 26 to 47%, and 52 to even 70%, respectively). In the case of the SILT fraction, most of these particles were deposited in the first two zones; however, seasonality had a large impact on this process (river-from 42 to 71%, and backwater-from 27 to 48%). In the case of CLAY, a similar seasonal pattern was observed, but most of these particles were deposited in the backwater zone and the Myślenice basin (23-47% and 30-40%, respectively). However, most important were the results concerning the particles of this fraction reaching the dam itself, and further to the lower part of the basin (even from 46 to 60%). This poses a direct threat to, e.g., drinking water, the intake of which is located in the last zone of the dam reservoir. The described results, together with their seasonal changes, can be analyzed in detail for each fraction in the developed graphical database [147].

Fingerprinting
The effectiveness of the SWAT module in simulating the transport of sediment particles, especially in the land phase of the basin, has been repeatedly confirmed. However, the problem was the absorption capacity of the finer sediment particles, which are responsible for the transport of pollutants such as heavy metals [148]. SWAT itself has limited ability to simulate such pollutants, as well as the ability to use them to determine the contribution of individual sediment sources. The fingerprinting method [72] has already been successfully combined with SWAT, but without using all the possibilities offered by these tools [30,149,150]. The effectiveness of such integration of tools can be appreciated, especially in areas under the pressure of large emitters of pollutants [151]. This is one of the reasons why the Wolnica basin (a fragment of the Raba basin) was selected for study ( Figure 2), located only 23 km from one of the largest urbanized areas in Poland. The Krakow agglomeration covers an area of approximately 3231 km 2 and is inhabited by approximately 1400 million people [152]; its concentration of particulate matter is one of the highest in Europe [153]. In the land phase of the Wolnica basin area, soil sampling sites were designated, and a suspended sediment sampler BS-2 ( Figure 4) was installed in the riverbed phase at the estuary [72].
reasons why the Wolnica basin (a fragment of the Raba basin) was selected for study (Figure 2), located only 23 km from one of the largest urbanized areas in Poland. The Krakow agglomeration covers an area of approximately 3231 km² and is inhabited by approximately 1400 million people [152]; its concentration of particulate matter is one of the highest in Europe [153]. In the land phase of the Wolnica basin area, soil sampling sites were designated, and a suspended sediment sampler BS-2 ( Figure 4) was installed in the riverbed phase at the estuary [72]. The results of these measurements broadened existing knowledge about this area and became part of the data module, which then fed the fingerprinting module. It was based on inputs from four sources of suspended sediments: grassland, forest, arable land, and residential areas, occurring in the Wolnica river basin. The research led to the selection of Pb, Zn, Cd, Cu, Ni, and Hg as the optimal fingerprints to solve the mass balance equation for suspended sediment and soil samples [72]. The temporal pattern of sediment yields for specific land-use types followed the SYLD distribution for the entire basin, with maximum values of 1.52 t/ha in May for grassland. The land-use sediment yields mainly reflected differences between surface areas for the particular land-use types, displaying the highest values of 139 t/m and 156 t/m for arable and residential land, respectively, and the lowest (2.3 t/m) for forest. Overall, adding the fingerprinting module to the digital platform-DNS/SWAT made it possible to demonstrate the dependence of SYLD distribution on meteorological and topographic conditions. Only in the case of forest areas can low SYLD values be expected. The fingerprint module shows that forests significantly change The results of these measurements broadened existing knowledge about this area and became part of the data module, which then fed the fingerprinting module. It was based on inputs from four sources of suspended sediments: grassland, forest, arable land, and residential areas, occurring in the Wolnica river basin. The research led to the selection of Pb, Zn, Cd, Cu, Ni, and Hg as the optimal fingerprints to solve the mass balance equation for suspended sediment and soil samples [72]. The temporal pattern of sediment yields for specific land-use types followed the SYLD distribution for the entire basin, with maximum values of 1.52 t/ha in May for grassland. The land-use sediment yields mainly reflected differences between surface areas for the particular land-use types, displaying the highest values of 139 t/m and 156 t/m for arable and residential land, respectively, and the lowest (2.3 t/m) for forest. Overall, adding the fingerprinting module to the digital platform-DNS/SWAT made it possible to demonstrate the dependence of SYLD distribution on meteorological and topographic conditions. Only in the case of forest areas can low SYLD values be expected. The fingerprint module shows that forests significantly change the signals of pollution sent by the system. Their location is important as they buffer the terrain-dependent flow of pollutants directly into the river.

Conclusions
This review explores the available mathematical tools for tracking sediment transport in river basins, especially considering problems in mountainous areas such as the Western Carpathians. Particular attention was paid to the limitations of using single models such as SWAT, and the possibility of eliminating them through the still-rare practices of combining subsequent models, methods, and databases/data processing in order to use their best elements. This article allows one to follow the effects of several years of research using the digital platform-DNS and its most important module-SWAT. This tool was designed to simulate, in a digital space, a model of many natural and anthropogenic processes in river basins responsible for the erosion, transport, and deposition of sediment particles. However, its effectiveness may be different for the land and riverbed phase of the basin area, which may have negative consequences in the form of limited confidence in the obtained results. SWAT also has limitations in tracking sediment deposition in dam reservoirs as well as particle-bound contaminant transfer. The use of the digital platform-DNS allowed us to broaden the spectrum of model research, first by combining the SWAT module with the fingerprinting method, and a database with a suspended sediment sampler. The result was a useful tool for tracking spatial and temporal changes in sediment yield and contaminants transported on soil particles. Both meteorological and topographic conditions were indicated as the main factors influencing the sediment yield distribution, as well as forest areas, which were the only ones distinguished by low SYLD values. In the next stage, the combined SWAT and AdH/PTM modules showed the temporal and spatial sequence of sedimentation in each zone of the dam reservoir, confirming its favorable layout and location that effectively protects the drinking water intake. The threats related to the growing amount of CLAY fraction, susceptible to pollutant absorption, reaching the last zone of the tank were also revealed. The open nature of the digital platform-DNS, and thus the possibility of using the available tools as subsequent modules provide the basis for the elimination of further limitations related to the tracking of sediment transport. Research results already available indicate that the highly uncertain simulation of sediments in the riverbed phase in the SWAT module could be successfully replaced by using a more robust physics-based approach from models such as HSPF [44]. Another problem is urbanized areas, the rapid development of which is difficult to represent in models due to the limited access of data constituting the basis for constructing variant scenarios. Meanwhile, a combination of tools such as SWAT and, for example, CLUE-S (Land Use Conversion and its Impact on a Small Regional Coverage) [43] may reduce these limitations. This will allow the effects of future land-use change scenarios to be simulated based on, inter alia, historical trends. Finally, there is a problem with database modules, especially those related to precipitation, whose current temporal resolution is insufficient and requires the use of less precise calculation methods. The solution may be the so-called precipitation atlases with high-resolution data added to the digital platform-DNS, which will significantly improve the quality of simulations in both phases of the basin. Despite the fact that even the most advanced models are not ideal, they seem to be the best tools for tracking sediment transport in river catchments today. Further work leading to an increase in the frequency of monitoring and the improvement of models and the possibility of combining them should therefore be one of the priorities of water management in the coming years.