Next Article in Journal
Development of an IoT System for the Generation of a Database of Residential Water End-Use Consumption Time Series
Previous Article in Journal
The New Set Up of Local Performance Indices into WaterNetGen and Application to Santarém’s Network
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Proceeding Paper

Comparison between Calculation and Measurement of Total Sediment Load: Application to Nestos River †

by
Loukas Avgeris
1,*,
Konstantinos Kaffas
2 and
Vlassios Hrissanthou
1
1
Department of Civil Engineering, Democritus University of Thrace, 67100 Xanthi, Greece
2
Faculty of Science and Technology, Free University of Bozen-Bolzano, 39100 Bolzano, Italy
*
Author to whom correspondence should be addressed.
Presented at the 4th EWaS International Conference: Valuing the Water, Carbon, Ecological Footprints of Human Activities, Online, 24–27 June 2020.
Published: 13 August 2020

Abstract

:
Measurements of stream discharge, bed load transport rate and suspended sediment concentration in the Nestos River (northeastern Greece) were conducted by the Section of Hydraulic Engineering, of the Civil Engineering Department, Democritus University of Thrace. In addition to those measurements, the total sediment concentration was calculated by means of the formulas of Yang. The comparison between the calculated and measured total sediment concentration was achieved by means of several statistical criteria and the results were deemed satisfactory.

1. Introduction

Stream sediment transport still stands out as a challenging problem for hydraulic engineers. Despite the leaps that have been made in the last century in the understanding and modeling of river load transport, the problem remains largely unintelligible and insoluble, due to the complexity of the physical processes that describe it.
Sediment transport affects riverine systems, either directly or indirectly, through erosion and deposition. Excessive depositions, as a result of soil erosion, affect the cross sections, increase the risk of flooding and can lead to a deterioration of water quality. This problem can be exacerbated in the case of agricultural basins where sediments may be carriers of infectious particles due to the use of fertilizers and pesticides [1].
Total sediment transport in streams is classified into bed load transport and suspended load transport on the basis of two different motion patterns. The sum of bed load and suspended load is equal to the total load [2].
In the recent decades, the Section of Hydraulic Engineering of the Civil Engineering Department, Democritus University of Thrace (Greece), has conducted bed load transport rate measurements and suspended sediment concentration measurements at the outlet of the Nestos River basin [3], Kosynthos River basin [4,5], and Kimmeria Torrent basin [5].
This study aims to redetermine the coefficients of Yang’s (1973) and Yang’s (1979) sediment transport formulas using multiple regression. In the past, nonlinear regression equations between bed load transport rate and stream flow rate, as well as between suspended load transport rate and stream flow rate, were established for the outlet of the Nestos River basin [3]. In the present study, 111 pairs of measured stream flow rate, measured bed load transport rate or measured suspended load transport rate in the Nestos River were used. The sum of measured bed load transport rate and measured suspended load transport rate provides the measured total load transport rate from which the total sediment concentration can be estimated. Apart from those measurements, total sediment concentrations were calculated by means of Yang’s formulas [6,7]. This made it possible to compare calculated to site-measured total sediment concentrations.

2. Study Area

The Nestos River springs from the Rila Mountains in southwestern Bulgaria and is one of the main watercourses of Eastern Macedonia and Thrace (Greece). Its Greek part covers approximately 130 km and the mountainous part of the Nestos River basin extends to an area of 840 km2. The basin is covered by forest (48%), bush (20%), cultivated land (24%), urban area (2%) and areas of no significant vegetation (6%), and has an altitude between 38 m and 1747 m. The basin is divided into 20 sub-basins with coverage areas between 13 km2 and 80 km2 and the mean land slope is between 23% and 58%. The mean slope of the main streams of the sub-basins ranges between 2.5% and 20%, whereas the mean slope of the Nestos River is 0.35%.

3. Stream Flow Rate and Sediment Transport Rate Measurements

All measurements were conducted at a location between the outlet of the Nestos River basin (Toxotes) and the river’s delta [8,9]. The average width of the cross sections of all measurements is 26.7 m.
The stream flow rate measurements were conducted using the following procedure: the site cross section was divided into sub-sections and the average stream flow velocity was measured at the middle of each sub-section, at 40%, approximately, of the flow depth from the bed, using a Valeport open channel flow meter. The stream flow rate of the entire cross section was taken as the sum of the individual sub-sections stream flow rates.
The bed load transport rate measurements were conducted in the middle of each cross section using a Helley–Smith bed load sampler. In order to determine the bed load transport rate, the trapped bed load sample is dried out and the mass is divided by the trap width and the measurement time duration [8,9].
The suspended sediment concentration was determined by obtaining a sample of water at the middle of the section and subsequently filtrating the sample through retention paper filters to obtain the net weight of the suspended load [8,9].

4. Calculation of Total Sediment Concentration

“Unit stream power”, as defined by Yang [6], is the amount of dynamic energy consumed by gravitational flow per unit of time and per unit weight of water, and is expressed by the product of the flow rate and the energy slope:
d Y d t = d x d t d Y d x = u s = u n i t stream power
where Y is the elevation above a datum, equal to the potential energy per unit weight of water; x is the horizontal distance; and t is the time.

4.1. Yang (1973)

Yang’s [6] formula for the total sediment transport in a river is given by:
log c F = 5.435- 0.286log w D 50 v - 0.457log u * w + ( 1.799- 0.409 log w D 50 v - 0.314log u * w ) log ( u s w - u cr s w )
where cF is the total sediment concentration in parts per million (ppm) by weight; w is the terminal fall velocity of sediment particles (m/s); D50 is the median particle diameter (m); ν is the kinematic viscosity of water (m2/s); s is the energy slope; u is the mean flow velocity (m/s); ucr is the critical mean flow velocity (m/s); and u* is the shear velocity (m/s).
If the following auxiliary variables x 1 , x 2 , x 3 , x 4 and x 5 are considered:
x 1 = log ( w D 50 / ν ) x 2 = log ( u * / w ) x 3 = log ( u s / w - ucr s / w ) x 4 = log ( u s / w - ucr s / w ) log ( w D 50 / ν ) x 5 = log ( u s / w - ucr s / w ) log ( u * / w )
then Yang’s formula can be written as follows:
log c F = 5.435- 0.286 x 1 - 0.457 x 2 + 1.799 x 3 - 0.409 x 4 - 0.314 x 5
White et al. [10] calculated the terminal fall velocity of the particles using the following equations:
w = F ρ g D
F = ( 2 3 + 36 D * 3 ) + 36 D * 3
D * = ρ g ν 2 3 D ch
ρ = ρ F - ρ w ρ w
where F is the correction factor for suspended load; D is the grain diameter (m); D* is the Bonnefille number; Dch is the characteristic grain diameter (m); ρF is the density of sediment particles (kg/m3); and ρW is the density of water (kg/m3). The kinematic viscosity v (m2/s) of water is given by the equation:
ν = 1.78× 1 0 - 6 1 + 0.0337T + 0.00022 T 2
where T (°C) is the temperature of the water.

4.2. Yang (1979)

In 1979, Yang concluded that the critical unit stream power term in Equation (2) can be neglected without causing much error when the measured sediment concentration is greater than 20 ppm by weight [7]. The simplified unit stream power equation was derived as:
log c F = 5.165 - 0.153 log w D 50 v - 0.297 log u * w + ( 1.780 - 0.360 log w D 50 v - 0.480 log u * w ) log ( u s w )
Similarly, if the following auxiliary variables   x 1 ,   x 2 ,   x 3 ,   x 4 and   x 5 are considered:
  x 1 = log ( wD 50 / ν )   x 2 = log ( u * / w )   x 3 = log ( us / w )   x 4 = log ( us / w ) log ( wD 50 / ν )   x 5 = log ( us / w ) log ( u * / w )
Equation (10) can be written as a linear multiple regression equation:
logc F = 5 . 165 - 0 . 153   x   1 - 0 . 297   x   2 + 1 . 780   x   3 - 0 . 360   x   4 - 0 . 480   x   5

5. Development of Yang’s Equations on the Basis of Nestos River Data

On the basis of the Nestos River data, the arithmetic coefficients of the original formulas of Yang [6,7], (Equations (2) and (10)), are modified, respectively, as follows:
logc F = 2 . 595 - 0 . 560 log wD 50 v - 6 . 649 log u * w - ( 1 . 380 - 0 . 534 log wD 50 v + 2 . 315 log u * w ) log ( us w - u cr s w )
logc F = 3 . 301 - 0 . 697 log wD 50 v - 14 . 367 log u * w - ( 1 . 214 - 0 . 537 log wD 50 v + 7 . 301 log u * w ) log ( us w )
In concrete terms, the new arithmetic coefficients of Equations (13) and (14) were determined by means of the conventional least square-based regression.
The measured stream flow rate (m3/s), the measured total sediment concentration (ppm), as well as the corresponding calculated total sediment concentration (ppm), by means of Equations (13) and (14), are provided in Table 1.

6. Comparison between Calculated and Measured Total Sediment Concentration

The comparison between calculated and measured total sediment concentration is made on the basis of the following statistical criteria [11]. At this point, it should be noted that the total sediment concentration was calculated by means of both the original and the modified Yang’s formulas.

6.1. Root Mean Square Error (RMSE)

RMSE = i = 1 n ( y i - y ^ i ) 2 n
where y i is the measured total sediment concentration; y ^ i is the calculated total sediment concentration and n the number of data. The RMSE ranges between 0 and +∞. The lower the RMSE, the better the correlation between measured and calculated values.

6.2. Mean Relative Error (MRE) (%)

MRE = i = 1 n ( y i - y ^ i ) n 100
Mean Relative Error (MRE) provides the relative size of the error. It is an index of how good an approximation between the predicted and measured value is, in relation to the magnitude of the physical quantity’s value.

6.3. Nash –Sutcliffe Efficiency (NSE) [12]

NSE = 1   -   i = 1 n ( y i - y ^ i ) 2 i = 1 n ( y i - y ¯ ) 2
where y ¯ is the average value of y i . NSE indicates how well the plot of observed versus simulated data fits the line of agreement (1:1 line). Nash –Sutcliffe efficiency ranges from −∞ to 1, with 1 being the optimal value.

6.4. Linear Correlation Coefficient r

r = i = 1 n ( y i - y ¯ ) ( y ^ i - y ^ ¯ ) i = 1 n ( y i - y ¯ ) 2 i = 1 n ( y ^ i - y ^ ¯ ) 2
where y ^ ¯ is the average value of y ^ i . The coefficient r expresses the degree of mutual linear dependence between the variables y i and y ^ i , and ranges between −1 and +1. The values r = ±1 represent the ideal occasion, when the marks representing the pairs of values y i and y ^ i depicted on an orthogonal coordinate system, lie on the regression line, with positive or negative slope, respectively.

6.5. Determination Coefficient R2

The determination coefficient R2 yields the percentage of change of the calculated values, which can be explained by the linear relationship between calculated and measured values. It ranges between 0 and 1. A value of 0 states that there is no correlation, whereas the value of 1 states that the variance of the calculated values equals the variance of the measured values [11].

6.6. Discrepancy Ratio

The discrepancy ratio represents the percentage of the calculated total sediment concentration values lying between pre-determined margins of the corresponding measured total sediment concentration values. As far as the present study is concerned, the discrepancy ratio represents the percentage of the calculated total sediment concentration values that lies between the double and the half of the corresponding measured total sediment concentration values.
The values of the above-mentioned statistical criteria are displayed in Table 2 and Table 3.
The values of the RMSE and NSE, on the basis of the calibrated formulas, can be considered fairly satisfactory. Additionally, the degree of linear dependence between calculated and measured total sediment concentration is acceptable. As expected, the values of NSE and R2, on the basis of the calibrated formulas, are identical and obviously non-negative.
The plot of Figure 1 represents the discrepancy ratio between measured and calculated values of total sediment concentration. At this point, it should be noted that both coordinate axes are in logarithmic scale; therefore, the equations y = x, y = 0.5x and y = 2x are graphically represented by parallel straight lines. Especially the values of the discrepancy ratio, on the basis of the calibrated formulas, are very satisfactory.

7. Discussion—Conclusions

In this paper, an attempt was made to redefine the coefficients of Yang’s formulas based on field measurements data in the Nestos River, between 2005 and 2015. A deviation between the calculated and measured total sediment concentration was observed for this specific case. For the correct application of Yang’s formulas [6,7] to the Nestos River, the calibration of the independent variables’ coefficients was deemed necessary.
As presented above, all statistical criteria of both calibrated Yang’s formulas were improved in comparison to the ones of Yang’s original formulas. More specifically, the RMSE approached zero for the calibrated equations, whilst the MRE displayed a notable decrease. Regarding the NSE, the linear correlation coefficient, r, and the determination coefficient, R2, came closer to the optimal value. Particularly, the discrepancy ratio was very near to the optimal value (100%). Overall, the results can be considered satisfactory.
It is noted that the application of Equations (13) and (14) should be bound to the Nestos River.

Author Contributions

Conceptualization, V.H.; methodology, V.H. and K.K.; software, L.A. and K.K.; validation, L.A. and K.K.; formal analysis, L.A., K.K. and V.H.; investigation, K.K. and L.A.; resources, V.H. and K.K.; data curation, L.A.; writing—original draft preparation, L.A.; writing—review and editing, V.H., K.K. and L.A.; visualization, L.A.; supervision, V.H.; project administration, V.H. All authors have read and agree to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kaffas, K.; Hrissanthou, V.; Sevastas, S. Modeling Hydromorphological Processes in a Mountainous Basin Using a Composite Mathematical Model and ArcSWAT. Catena 2018, 162, 108–129. [Google Scholar] [CrossRef]
  2. Hrissanthou, V.; Tsakiris, G. Sediment Transport in Water Resources: I. Engineering Hydrology; Tsakiris, G., Ed.; Symmetria: Athens, Greece, 1995; pp. 537–577. (In Greek) [Google Scholar]
  3. Angelis, I.; Metallinos, A.; Hrissanthou, V. Regression analysis between sediment transport rates and stream discharge for the Nestos River, Greece. Glob. NEST J. 2012, 14, 362–370. [Google Scholar]
  4. Kaffas, K.; Hrissanthou, V. Estimate of continuous sediment graphs in a basin, using a composite mathematical model. Environ. Process. 2015, 2, 361–378. [Google Scholar] [CrossRef]
  5. Metallinos, A.; Hrissanthou, V. Regression relationships between sediment yield and hydraulic and rainfall characteristics for two basins in northeastern Greece. In Proceedings of the 6th International Symposium on Environmental Hydraulics, Athens, Greece, 1–5 July 2010; Christodoulou, G., Stamou, A., Eds.; Volume 2, pp. 899–904. [Google Scholar]
  6. Yang, C.T. Incipient motion and sediment transport. J. Hydraul. Div. ASCE 1973, 99, 1679–1704. [Google Scholar] [CrossRef]
  7. Yang, C.T. Unit stream power equations for total load. J. Hydrol. 1979, 40, 123–138. [Google Scholar] [CrossRef]
  8. Konstantinopoulou-Pesiou, C.; Sfiris, D. Measurements of Stream Discharge and Sediment Discharge of Nestos River. Diploma Thesis, Department of Civil Engineering, Democritus University of Thrace, Xanthi, Greece, 2013. [Google Scholar]
  9. Mpenekos, K.; Kassotakis, E. Measurements of Stream Discharge and Sediment Discharge in Nestos River during May and June 2014. Diploma Thesis, Department of Civil Engineering, Democritus University of Thrace, Xanthi, Greece, 2015. [Google Scholar]
  10. White, W.R.; Milli, H.; Crabbe, A.D. Sediment Transport: An Appraisal of Available Methods. Volume 2: Performance of Theoretical Methods When Applied to Flume and Field Data; INT 119; Hydraulics Research Stations: Wallingford, UK, 1973. [Google Scholar]
  11. Krause, P.; Boyle, D.P.; Bäse, F. Comparison of different efficiency criteria for hydrological model assessment. Adv. Geosci. 2005, 5, 89–97. [Google Scholar] [CrossRef]
  12. Nash, J.E.; Sutcliffe, J.V. River flow forecasting through conceptual models, Part I—A discussion of principles. J. Hydrol. 1970, 10, 282–290. [Google Scholar] [CrossRef]
Figure 1. Discrepancy ratio plot between measured and calculated values of total sediment concentration by means of: (a) calibrated Yang’s formula (1973) and (b) calibrated Yang’s formula (1979).
Figure 1. Discrepancy ratio plot between measured and calculated values of total sediment concentration by means of: (a) calibrated Yang’s formula (1973) and (b) calibrated Yang’s formula (1979).
Environsciproc 02 00019 g001aEnvironsciproc 02 00019 g001b
Table 1. Measured stream flow rate and total sediment concentration—Calculated total sediment concentration in the Nestos River.
Table 1. Measured stream flow rate and total sediment concentration—Calculated total sediment concentration in the Nestos River.
No.Stream Flow Rate (m3/s)Total Load logcF (meas.)Total Load logcF (calc.) Yang 1973Total Load logcF (calc.) Yang 1979No.Stream Flow Rate (m3/s)Total Load logcF (meas.)Total Load logcF (calc.) Yang 1973Total Load logcF (calc.) Yang 1979
114.172.40642.32722.4124313.162.48522.12432.1413
217.442.39292.20862.3058322.561.65312.18182.1209
319.502.11602.13352.1819333.952.00332.01202.0609
416.652.57692.22802.3020344.222.19902.01492.0704
518.492.31632.15672.2284353.662.37542.00992.0559
62.441.37892.06252.1273364.802.22131.91351.9778
72.731.38061.95142.0150372.062.00661.70681.7695
82.691.28521.93211.9866381.462.24931.84681.9173
92.841.27121.89491.9532391.882.10661.71641.7854
103.091.34441.92602.0052401.492.16521.83661.9067
1117.890.93001.33301.4140411.752.16391.78151.8486
1215.451.50091.45831.5803421.661.74201.83331.9043
1320.621.40761.25191.3096432.291.72001.64671.7179
1416.151.27191.43271.5588441.552.16001.83861.8937
1514.140.88081.52991.6691451.241.89191.92831.9875
1658.98−0.49770.74540.4515461.652.63121.79261.8650
1752.940.32240.88560.7877471.561.59171.84061.9286
1850.141.08320.84050.5748482.031.10851.70151.7641
1948.270.62140.91000.7553490.801.47411.94742.0033
2045.720.53410.95700.8348500.691.01021.93471.9919
2144.450.98080.95860.7920510.691.18682.01702.0654
2262.411.26580.73420.4807520.901.57531.89871.9769
2355.300.57160.81220.5924533.272.18871.96112.0192
242.622.65401.91641.9354543.702.48671.77691.8147
253.952.61322.01492.0704552.482.23562.04372.0747
264.222.50231.99732.0650562.232.61002.05332.0860
274.132.47931.97762.0225570.852.02112.38162.2171
286.202.48701.84881.9273580.642.33242.08421.9663
294.802.67151.84851.8979590.552.26242.35192.2803
303.762.53851.97072.0092602.832.09211.95002.0009
613.402.27981.93551.9988870.901.94022.02991.9359
623.292.04111.91621.9825880.881.90071.68221.6011
631.772.76042.13732.1580890.971.42801.88351.8606
641.062.78362.48872.3048900.471.75402.24802.1496
650.602.56662.47852.1120910.521.82912.30182.2275
660.392.34392.76242.3083920.292.02082.31912.1865
670.642.84552.56432.2200931.390.08032.26722.2396
682.671.60122.01612.0783941.360.29722.03201.9492
693.682.02371.86931.9422950.932.42222.72062.7267
702.451.72862.04702.0960961.152.25052.59832.6155
712.622.51352.01092.0678972.052.68522.32362.3769
722.952.58912.01452.0770980.871.93592.21072.0739
730.841.79192.05952.0582990.852.40922.71462.6828
741.761.41191.74861.81041001.062.02242.12082.0151
751.812.08261.73121.78661011.262.46992.22452.1058
761.092.28631.99242.05591020.492.57592.13052.1176
770.592.32322.32342.25691030.342.77712.39372.3719
781.062.22651.78871.84051041.341.79051.96322.0090
790.752.32691.84851.86231050.682.44272.11922.1378
801.432.28861.66531.72451060.102.98591.91982.1113
811.471.86791.71911.792710717.281.81671.58571.6423
820.642.18521.98141.98691088.242.21062.02682.1808
830.551.49342.11931.99051090.942.06631.81491.9314
840.831.91762.14372.05441109.831.95491.71341.8098
850.851.74391.98731.92281110.762.04251.84851.9434
860.511.51081.96651.9298
Table 2. Statistical criteria of Yang’s formulas—original and calibrated (1973).
Table 2. Statistical criteria of Yang’s formulas—original and calibrated (1973).
RMSEMRE (%)NSErR2Discrepancy Ratio
Original1.324−98.339−3.240−0.3380.1140.757
Calibrated0.506−31.7340.3810.6170.3810.964
Table 3. Statistical criteria of Yang’s formulas—original and calibrated (1979).
Table 3. Statistical criteria of Yang’s formulas—original and calibrated (1979).
RMSEMRE (%)NSErR2Discrepancy Ratio
Original1.390−115.039−3.673−0.4490.2010.739
Calibrated0.492−31.1850.4150.6440.4150.955
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Avgeris, L.; Kaffas, K.; Hrissanthou, V. Comparison between Calculation and Measurement of Total Sediment Load: Application to Nestos River. Environ. Sci. Proc. 2020, 2, 19. https://0-doi-org.brum.beds.ac.uk/10.3390/environsciproc2020002019

AMA Style

Avgeris L, Kaffas K, Hrissanthou V. Comparison between Calculation and Measurement of Total Sediment Load: Application to Nestos River. Environmental Sciences Proceedings. 2020; 2(1):19. https://0-doi-org.brum.beds.ac.uk/10.3390/environsciproc2020002019

Chicago/Turabian Style

Avgeris, Loukas, Konstantinos Kaffas, and Vlassios Hrissanthou. 2020. "Comparison between Calculation and Measurement of Total Sediment Load: Application to Nestos River" Environmental Sciences Proceedings 2, no. 1: 19. https://0-doi-org.brum.beds.ac.uk/10.3390/environsciproc2020002019

Article Metrics

Back to TopTop