Next Article in Journal
A Modeling Platform for Landslide Stability: A Hydrological Approach
Next Article in Special Issue
An Experimental Study on the Effects of Atomized Rain of a High Velocity Waterjet to Downstream Area in Low Ambient Pressure Environment
Previous Article in Journal
A Spatiotemporal Analysis of Nitrogen Pollution in a Coastal Region with Mangroves of the Southern Gulf of Mexico
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Application of Asymmetrical Statistical Distributions for 1D Simulation of Solute Transport in Streams

1
Department of Sanitary and Environmental Engineering, Faculty of Civil Engineering, Slovak University of Technology in Bratislava, 2412 Bratislava, Slovakia
2
Institute of Hydrology, Slovak Academy of Sciences, 2412 Bratislava, Slovakia
3
Department of Civil, Architectural and Environmental Engineering, University of Napoli Federico II, 80100 Napoli, Italy
*
Author to whom correspondence should be addressed.
Water 2019, 11(10), 2145; https://doi.org/10.3390/w11102145
Submission received: 12 September 2019 / Revised: 4 October 2019 / Accepted: 11 October 2019 / Published: 15 October 2019
(This article belongs to the Special Issue Advances in Environmental Hydraulics)

Abstract

:
Analytical solutions of the one-dimensional (1D) advection–dispersion equations, describing the substance transport in streams, are often used because of their simplicity and computational speed. Practical computations, however, clearly show the limits and the inaccuracies of this approach. These are especially visible in cases where the streams deform concentration distribution of the transported substance due to hydraulic and morphological conditions, e.g., by transient storage zones (dead zones), vegetation, and irregularities in the stream hydromorphology. In this paper, a new approach to the simulation of 1D substance transport is presented, adapted, and tested on tracer experiments available in the published research, and carried out in three small streams in Slovakia with dead zones. Evaluation of the proposed methods, based on different probability distributions, confirmed that they approximate the measured concentrations significantly better than those based upon the commonly used Gaussian distribution. Finally, an example of the application of the proposed methods to an iterative (inverse) task is presented.

1. Introduction

Solute transport in streams and rivers is strongly related to river characteristics, such as mean flow velocity, velocity distribution, secondary currents, and turbulence features. These parameters are mainly determined by the river morphology and the discharge conditions. Most natural channels are characterized by the relevant diversity of morphological conditions. In natural channels, changing river width, curvature, bed form, bed material, and vegetation are the reason for this diversity. In rivers, which are regulated by man-made constructions such as spur dikes, groins, stabilized bed, and so on, the morphological diversity is often less pronounced, and thus, flow velocities are more homogeneous.
In natural channels, some of these morphological irregularities, such as small cavities existing in sand or gravel beds, side arms and embayments, bigger obstacles, bank vegetation, and uprooted trees, can produce recirculating flows, which occur on different scales on the riverbanks and the riverbed. These irregularities act as dead zones for the current flowing in the main stream direction. In regulated rivers, groyne fields are the most important sort of dead zones. Groyne fields can cover large parts of the river, significantly affecting its flow field. Dead-water zones or dead zones can be defined as geometrical irregularities existing at the river periphery, within which the mean flow velocity in the main stream direction is approximately equal to zero (or even negative) [1,2].
Dead zones significantly modify velocity profiles in the main channel and also affect dispersive mass transport within the river by collecting and separating part of the solute from the main current [3,4]. Subsequently, the solute is slowly released and incorporated back into the main current in the stream. In addition, hyporheic exchanges between the surface waters and the adjacent/underlying groundwater may affect the transport of solutes resulting in long-time tailing behavior in the time–concentration curves [5]. As the classical advection–diffusion equation (ADE) does not consider explicitly the effect of dead zones, the concept of transient storage was introduced [6], coupling the ADE to a second equation describing the exchange between the main stream and the dead zones (see Theory).
Alternatively, the aggregated dead zones (ADZ) model was proposed [7], in which the cumulative effects of all the physical dispersive mechanisms taking place in a river stretch are related to the residence time of a single, idealized transient storage zone, the ADZ. The ADZ model is, hence, a conceptually simple, but less well-physically based approach, which was found to be robust.
Mathematical solutions for the dispersion in streams with dead zones following the transient storage approach have been extensively described [8,9,10,11,12], as well as the effects of the vegetation in streams [13,14,15]. However, the majority of these solutions are complicated due to the mathematical apparatus and their numerical applications [9,10,16,17]. Code programming could be difficult, especially in cases where the computational effort is intensive (e.g., inverse tasks). On the other hand, based on the nature of analytical solutions, the superposition principle can be applied, so the analytical solutions are widely used in engineering practice because of their ability to simulate various types of solute inputs (continuous, discontinuous, steady, unsteady, etc.).
The main objective of this paper was to provide a simple mathematical method that would allow a quick and easy one-dimensional modelling of the dispersion process in streams with the influence of dead zones. This method is based on use of various probability distributions in order to achieve better approximation of the measured concentration from field experiments. Four distributions, namely the Gaussian distribution, the Gumbel distribution, the log-normal distribution and the generalized extreme value (GEV) distribution, were considered to approximate the dispersion process from an instantaneous injection. The Gaussian distribution was considered because it is regarded as a “standard” solution, used in many scientific studies, while the Gumbel distribution was considered for its simplicity and literature references [18]. Both Gaussian and Gumbel distributions are two-parameter distributions. Assuming higher precision using the third additional parameter, we evaluated the log-normal and GEV distributions.
The results from these distributions were compared with the experimental data available in the literature, and with those collected in three streams in southern Slovakia, where a significant presence of dead zones was found (Malá Nitra, Šúrsky, and Malina streams). An application of these methods to an inverse task is finally presented.

2. Theory

The general form of the one-dimensional advection–dispersion equation, which is valid once cross-sectional mixing has been achieved, can be found in the literature, e.g., [19] in the form
C t + v x C x = D x ( 2 C x 2 ) + M s
where t is the time (s), C is the concentration of a substance (kg.m−3), Dx is the dispersion coefficient in the longitudinal direction (m2 s−1), vx is water velocity in the longitudinal direction (m.s−1), Ms is a function representing the sources of solute (kg.m−3 s−1), x is the spatial coordinate distance (m).
It is well known that downstream of the release of a tracer into a turbulent flow, two zones are established: the advective zone and the Gaussian (or equilibrium) zone [20]. The former is defined as the distance it takes a tracer particle to sample the entire flow field, and it is characterized by a skewed concentration profile due to the velocity profile in the channel. After a sufficient time, defined by a Lagrangian timescale [21], the differential longitudinal advection and the cross sectional mixing attain a balance and the skewness decays until the concentration profile is Gaussian and the variance of the tracer could increase linearly with the time [20]. However, in some cases advection dominates dispersion and the concept of frozen cloud is used by assuming that no longitudinal dispersion occurs during the time taken for the tracer to pass a sample site [20,22].
For contaminant transport in streams with transient storage, a term expressing the exchange between the transient storage (dead zones) and the main stream is usually added in (1). Furthermore, an equation for the transport inside the dead zone is needed. If we can assume the water flow and transport parameters as constant, the resulting equations can be written in the form [9,23,24]:
C t = D x ( 2 C x 2 ) v x C x m ( C C s ) + M s
n C s t = m ( C C s )
where Cs is the concentration of the substance in the storage zone (kg.m−3), m is the mass exchange coefficient between the main stream and the storage zone (s−1), and n is the ratio between the storage zone and the main stream cross-sectional area (–). If m or n becomes zero, the (2) reduces to (1).
A general form of solution for diffusion Equation (1) by [19], and eventually by [25,26] could be written as
C =   M A D x t   f ( x D x t )
where A is the cross-sectional area of the stream (m2) and f is the unknown function (“similarity solution “) [19]. The unknown function f can be determined in two different ways:
  • Based on experiment derive a curve fit to real data.
  • By analytically solving the (4).

2.1. Analytical Solution

Using the second method, we get the one-dimensional analytical solution of the advection–dispersion equation for simplified conditions and immediate solute input in the form [26]
C ( x , t ) = M 2 A π D x t   e x p ( ( x v ¯ t ) 2 4   D x   t )
where v ¯ is the average flow velocity in the stream. The simplified conditions comprise a prismatic streambed, steady and uniform flow, and conservative substance, as well as assumption of lateral and vertical homogeneous concentration of the flow (1D), so theoretically, this equation can be applied after achieving cross-sectional mixing in the river. Initial and boundary conditions are typically zero concentration for the solute in stream and no lateral inflows, but using the superposition principle and dividing the stream to sections, different (nonzero) initial and boundary conditions could also be imposed.
By comparing (4) with (5), the analytical solution of the function f could be interpreted as having a form of Gaussian normal distribution with parameters of normal distribution (e.g., standard deviation σ, etc.). In fact, this curve does not have the exact shape of a Gaussian normal distribution, rather it has an asymmetric form, caused by substitution of the standard deviation σ (which is constant in Gaussian distribution) by the time-dependent term 2 D x t .
Analytical solution of advection–dispersion equation (5) is currently considered as the standard analytical solution, but in our opinion the applicability of this solution is limited to flow without barriers, with the assumption of a symmetrical spreading of particles in a longitudinal direction. In real streams, however, this assumption is not always right; movement of particles is slowed down (delayed) by the accumulation of particles in dead zones of the streams.
For this reason, in the case of flow in real conditions with occurrence of dead zones in the stream, it is not appropriate to use the Gaussian normal distribution to approximate the function f from (4). Instead a different statistical distribution with asymmetric shape should be used.
After reviewing the literature, we selected four statistical distributions for the present study. Our selection was based mainly on the criterion of simplicity, easy mathematical formulation, and the number of distribution parameters. The first distribution is the “classical” Gaussian distribution (GaussD), according to (5). As mentioned earlier, this is considered the “standard solution”. Inspired by the work of [18] we selected two other distribution types—the Gumbel distribution (GumbelD) and the log-normal distribution (LogNormD)—modified for three parameters. The log-normal distribution was also used in previous studies and it was stated that instream advection and dispersion processes can produce only log-normal concentration distributions [12]. Finally, the generalized extreme values distribution (GEVD) was considered. Notably, the first two distributions (GaussD, GumbelD) are two-parametric, whereas the last two (LogNormD and GEVD) are three-parametric distributions. It should also be mentioned that the Gumbel distribution is a special case of the GEV distribution.
Other potential distributions (models) were the TOSS model [27] and, finally, the generalized inverse Gaussian distribution [28]. However, these are too complicated with many embedded functions (e.g., modified Bessel functions of the second kind, etc.) and complex expressions for the distribution of specific parameters and mean, median, or mode determination. Another possible distribution was the inverse Gaussian distribution [29], but this has two parameters only, so the accuracy would presumably be similar to the GumbelD approximation method. Finally, another potential approach found in the literature is that based on signal deconvolution to obtain the residence time distribution in a stream [30].

2.2. The GumbelD Approximation Method

As an appropriate approximation of the function f from the Equation (4), the formula of the Gumbel distribution is proposed. The general equation of the Gumbel distribution is [31]
p = 1 σ   e ( z + e z )
z = x μ σ
where p is the distribution probability (density), µ is the location parameter, and σ is the scale parameter.
The parameters from Equation (6), (considering Equation (4)) can be defined as follows:
σ = D x , G t
z = v ¯ . t x D x , G t
where Dx,G is the dispersion coefficient in the longitudinal direction (m2 s−1) used in the proposed model. To be dimensionally consistent, z is a dimensionless parameter and σ has the dimension of length (m). To get the results in concentration units (kg m−3) the term M/A must be included into (6). Substituting parameters from (8) and (9) into (6), the one-dimensional analytical solution has the form:
c ( x , t ) =   M A D x , G t exp   [ x v ¯ . t D x , G t exp ( x v ¯ . t D x , G t ) ]

2.3. Three-Parametric LogNormD Approximation Method

We consider, as an appropriate approximation of the function f from the (4), the expression of the three-parametric log-normal distribution. The general equation of this distribution is [32]
p = 1 σ ( x θ ) 2 π   e [ [ ln ( x θ ) μ ] 2 2 σ 2 ]       ,     f o r   x > θ p = 0   ,    f o r   x   θ  
where p is the distribution probability (density), µ is the location parameter (peak location), σ is the variance (scale) parameter, and θ is the shift (lower boundary) parameter.
The time tmax, in which the maximum (peak) concentration occurs is defined as
t m a x = e ( μ σ 2 )
As the peak time tmax is known from field measurements, Equation (12) yields
μ = ln ( t m a x t 0 ) + σ 2
It could be assumed that the parameter t0 (start time of nonzero concentrations) can be expressed as a proportion of the peak (maximal concentration) time, i.e., t0 = k tmax, where 0 < k < 1.
c ( x , t ) =   M A   ( t t 0 ) σ 2 π   e x p ( ( ln ( t t 0 ) μ ) 2 2 σ 2 )
c ( x , t ) =   M A   ( t t 0 ) D x , L N 2 π   e x p ( ( ln ( t t 0 ) μ ) 2 2 D x 2 )
c ( x , t ) =   M A   ( t k . t m a x ) D x , L N 2 π   e x p ( ( ln ( t k . t m a x ) μ ) 2 2 D x 2 )
μ = ln ( t m a x k . t m a x ) + σ 2 = ln ( t m a x   ( 1 k ) ) + σ 2

2.4. Generalized Extreme Value Distribution (GEVD) Method

Another approximation of the function f from (4) is the generalized extreme value (GEV) distribution. The general equation of the GEV distribution is [33]
p = 1 σ z ( x ) ξ + 1 e z ( x )   where   z ( x ) =   { ( 1 + ξ ( x μ σ ) ) 1 ξ         i f   ξ 0 e   x μ σ                  i f   ξ = 0
where p is the distribution probability (density), µ is the location parameter (peak location), σ is the variance (scale) parameter, and ξ is the shape parameter.
The peak time tmax is defined as
t m a x = μ + σ ( 1 + ξ ) ξ 1 ξ   for   ξ 0 .
In practical computations, we can almost always assume the case ξ ≠ 0, so then
μ = t m a x σ ( 1 + ξ ) ξ 1 ξ .
The variance (scale) parameter σ can be defined similar to previous cases (see (8))
σ =   D x , G E V t .
Then
c ( x , t ) = M A D x , G E V t   z ( t ) ξ + 1   e z ( t ) where z = ( 1 + ξ ( v ¯ . t x D x , G E V t ) ) 1 ξ .

3. Methods

3.1. Literature Data Use

The literature dataset from [34] was used to evaluate all approximation methods described above. They presented measurements from 51 field tracer experiments performed between 1959 and 1973. In the present study, we excluded the tracer experiments performed on short lengths i.e., with lengths less than 10 kilometers. It can be assumed, that in such cases the results are distorted by the influence of the initial transversal and vertical mixing. We also omitted all the experiments with noncentral tracer inputs (i.e., tracer input at the bank, input at multiple points, etc.) because of the lack of data for these particular tracer experiments.
As a result, 29 tracer experiments with overall measurements in 141 profiles (there were up to 8 measuring profiles in each tracer experiment) were selected. Statistical evaluation was performed for each profile using each of the applied approximation methods (GaussD, GumbelD, LogNormD, and GEVD). This focused on finding the best fit between the measured values and the approximation with the considered method, i.e., to find the optimal set of method parameters providing the best goodness of fit based on the minimal root mean square error (RMSE) between the measured and approximated values, i.e.,
R M S E   =   t = t 1 t = t 2 ( c m , t c a , t ) 2 n
N R M S E   =   R M S E c m a x c m i n
where RMSE is the root mean square error (kg.m−3), cm,t is the measured value (concentration) at the time t (kg.m−3), ca,t is the approximated value in the time t (kg.m−3), t1 is the measurement start time (sec) and t2 is the measurement end time (sec), NRMSE (Eq. 24) is the normalized root mean square error (–), and cmax and cmin are the maximal and minimal concentration in each particular experiment (kg.m−3).
The optimal set of parameters were analyzed and set up using the built–in Excel Solver function. The target value was the RMSE and the variable parameters were the peak time (mode), the scale parameter (to replace the term M/A in corresponding equations), and the corresponding method parameters (dispersion coefficient variance, shape parameters). The target value was the main indicator for the optimization process, looking for a minimum value for this indicator to achieve the maximum conformity of the measured and modelled data for the evaluated parameters. Limits for the peak time for the optimization process were typically set up from 0.8 up to 1.2 of the real (experimental) peak time to give some variability for the optimization process. Theoretically, there should be no limits, but by setting such limits, the convergence of the optimization procedure was secured. The scale parameter limits were set up to have a maximal difference of 0.001 (1‰) difference between the sum of the tracer in modelled and experimental values. There were just formal limits for other parameters, e.g., the dispersion coefficient should a nonzero and positive value, the lower boundary coefficient k (see (17) should be positive and less than 1 (0 < k < 1).
The considered rivers exhibited a large variability of hydraulic conditions. The length of the experiments was in the range from 2.5 up to 294 km (1.6–183 miles), while the range of discharges and velocity was from 0.48 up to 6824 m3.s−1 (17–241000 cfs—cubic feet per second) and from 0.07 up to 1.86 m.s−-1, respectively. All the data are listed in Table 1.
Examples of the graphical evaluation of the considered distributions are shown in the Figure 1, Figure 2 and Figure 3.

3.2. Field Experiments Data Use

The ability of the methods to approximate the 1D dispersion in extreme conditions of the channel with submerged vegetation in a large portion of the cross-sectional area was also tested.
For these tests, experimental data from three tracer studies in Malá Nitra, Šúrsky, and Malina streams conducted between 2012 and 2016 were used. These tracer experiments focused on the determination of longitudinal and transversal dispersion coefficients and on pollutant spreading within the primary mixing zone (prior to the transversal and vertical homogenization). In the first two streams, kitchen salt was used as tracer and the water conductivity was measured. For the experiment on the Malina stream we used a coloring agent (E133, brilliant blue) as a tracer, with subsequent spectrophotometric detection of agent concentration. In all the cases evaluated in this paper, an instantaneous input of the tracer was applied.
The discharge in all experiments was determined using the cross-section velocity method. For velocity determination the device SonTeK Flow 3D tracker was used. Water level slopes were determined using the GPS precise levelling method (with ±10 mm accuracy), and based on the water level slope, the roughness (Manning) coefficient was calculated using the Manning’s formula [35].
The Malá Nitra stream is located within the village of Veľký Kýr (N48.181799°, E18.155373). The experiments were performed in two reaches with lengths of 785 and 1340 meters. The first reach of the stream was straight and the second was slightly curved in both directions (left and right bend). The channel at the bed was about 4 m wide with a bank height of 2.5 m and a bank slope of approximately 1:2. The discharge in selected experiments used for the method tests ranged from 0.230 up to 0.235 m3.s−1. The hydraulic roughness was n = 0.035 (Manning coefficient) and the water level slope was found to be approximately 1.5‰. The stream in the considered reach is prismatic with a width of 5.5 m and a water depth in the range from 0.4 to 0.6 m. The range of the determined longitudinal dispersion coefficient was from 0.5 to 2.5 m2.s−1. The experiments were performed in high summer (July 2012), when submerged vegetation was present to a large extent.
The second tracer study was carried out on a straight reach of the Šúrsky stream, located close to the village of Svätý Jur (Slovakia, N48.232957°, E17.202934°). The field measurements were made in a 300 m to 500 m long straight reach with a relatively prismatic cross-section profile. The discharge, the channel width, the maximum depth, and the average velocity were in the range 0.38 to 0.43 m3.s−1, 4 to 5.5 m, 0.4 to 0.8 m, and 0.245 to 0.282 m.s−1, respectively. The range of the determined longitudinal dispersion coefficient was from 0.63 to 0.98 m2.s−1. Field experiments were performed in spring (March 2014) when the presence of submerged vegetation was very low.
The last field tracer experiment was performed at the Malina stream in West Slovakia in 2016. Originally, the constructed cross-section profile was significantly influenced by vegetation. The measured discharge during the experiments was 0.408 m3s−1. The water level slope, specified by levelling measurements, was about 0.45‰. The stream shape in the examined stream section can be considered prismatic, the width was around 5 m, the average depth was 0.88 m with a maximum value of about 1 m. The dispersion coefficient was 0.95 m2.s−1. The experiment was performed in high summer (August 2016), with up to approximately one meter of emerged vegetation on the banks, as well as submerged vegetation in the central part of the channel. This vegetation on the stream banks and bed created dead zones for the flow. A picture of the vegetation and hydraulic conditions in this stream is shown on Figure 4.
The main characteristics of all three streams used for tracer experiments in Slovakia are listed in Table 2.
In all the three rivers used for field experiment data, the length was smaller than the length of rivers from the literature data. Therefore, the question of the 1D well-mixed assumption should be questioned.
Vertical mixing in a river needs, typically, a length of Lmix-vert = 100*h [20], where h is the water depth. In the three Slovakian rivers, water depth h was always < 1 m (on average, 0.5, 0.6, and 0.88 m). So, Lmix-vert < 100 (from 50 to 88 m). This means in the field measurement the assumption of vertical mixing holds.
The estimation of the length for transversal mixing Lmix-transv is more difficult because no established theory exists to predict transverse mixing rate [36]. This rate is related to several boundary conditions of the mixing process (riverbank/riverbed irregularities, vegetation, meandering, bends, islands, confluences). The vegetation is expected to affect Manning coefficient [37], which is related to the friction factor, the shear velocity, and, ultimately, to the transverse mixing rate.
For scaling analysis, again, the length for transverse mixing is Lmix-transv = c1*(W²/*h), where W is channel width and c1 is a numerical constant related to the rate of the process of transversal mixing [38]. In the three Slovakian rivers, c1 could be in the order of 20–40, so the Lmix-transv should be in the order of 150–200 m, which is quite a lot shorter than the minimum distance of the field experiments (Table 2). Hence, the assumption of complete mixing in the cross-section should be valid.

4. Results

4.1. Literature Data

The evaluation of the results from the approximation methods is based on comparison of the normalized root mean square error (NRMSE), which is the standard root mean square error (RMSE) divided by the maximal difference of measured values (cmax − cmin) [39]. Since the minimal concentration in all experiments was set up to be zero, the NRMSE = RMSE/cmax. For all experiments, the NRMSE was determined to evaluate the fit between the measured and approximated values. Results are listed in Table 3.
If the NRMSE for the GaussD method will be regarded as 100%, values less than 100% demonstrate that the approximation method was more accurate than the GaussD method, whereas values over 100% demonstrate that the GaussD method was more accurate. Such relative results are shown in Table 4.
The number of values exceeding the assigned accuracy limit are listed in Table 5. An accuracy limit, e.g., 40%, means that the NRMSE of the particular method was less than 40% of the NRMSE of the GaussD method (=100%). So, for example, for the 40% accuracy limit, the GEVD method was more precise than the GaussD method in 76 of total 141 experiments (53.9%). Accuracy limit of 100% means the same precision as the GaussD method.
Further evaluation focused on the relationship (correlation) between the Gaussian (GaussD method) dispersion coefficient and the dispersion coefficient for examined approximation methods. The results are listed in Table 6.
Another subject of this study was the question of whether any dependency between each method performance and stream conditions (distance, hydraulic parameters, etc.) exists. Results of such statistical evaluation are listed in Table 7. Similarly, the correlation between the stream conditions and the values of parameters used in each approximation method were determined (Table 8).
The ranges of dispersion parameters for each approximation method are presented in Table 9.

4.2. Field Experiments

Table 10, Table 11, Table 12 and Table 13 list the results of the proposed analysis for the field experiments. The comparison was that which was already used in the previous section.
The results of the approximation methods for the Malina stream, where a large portion of submerged and emerged vegetation was observed, are shown in Figure 5 and Figure 6. It is also noted that the considered methods were able to simulate the 1D dispersion process in such conditions.

5. Discussion

5.1. Literature Data

As it can be seen from Table 3, the approximation accuracy of the examined approximation methods was generally better than that of the GaussD method. An average value of NRMSE for each approximation method from all experiments was determined. If we assume the average NRMSE of the GaussD method to be 100%, the relative NRMSE average of the GumbelD method was 53.1%, the LogNormD and GEVD methods were even better—40.6 and 38.2%, respectively (see Table 4). The extreme (min, max.) values of the NRMSE were caused mainly by the extreme accuracy (or inaccuracy) of the GaussD approximation method with the measured values, i.e., in cases where the NRMSE was lower than 2% or larger than 10%.
The standard deviation (the 68.3 percentile) was also significantly better for the LogNormD and GEVD methods—51.3 and, 46.2%, respectively vs. 67.6% for the GumbelD approximation method.
Table 3 shows that all the three approximation methods were better than the GaussD method, but, in turn, the LogNormD and the GEVD were clearly better than the GumbelD approximation method. The GEVD approximation method was better almost in the whole limit range (0%–100%) than the LogNormD method. In particular, the LogNormD and GEVD methods were less accurate than the GaussD method in only two cases (1.4%) and three cases (2.1%), respectively. This happened in cases, when the time–concentration curve was perfectly symmetric (presumably a prismatic streambed with no obstacles), so the GaussD method could be fully applied.
Further results were obtained in the correlation examination between the specific approximation methods. Table 6 shows a very high correlation between the dispersion coefficients of GaussD, GumbelD, and GEVD methods (r2 > 0.99), whereas all the correlations of the LogNormD method were quite weak (r2 was about 0.2). This means that the dispersion coefficients for the GumbelD and GEVD methods were very close to our understanding of the role of the dispersion coefficient, as identified from the standard Gauss–based analytical solution of the ADE (Equation (5)). Even the values of the dispersion coefficient between the GaussD and GEVD methods were quite similar, probably due to a linear dependency between them.
The weak correlation between the LogNormD and the GaussD method should be considered a major disadvantage of the LogNormD method, as it will create practical difficulties when using this method. This is due to the dependence of the form of the approximation curve on the parameter t0 (see (15), whereas the dispersion coefficient is very difficult to estimate due its relatively small range (average 0.413, std. dev. 0.1). On the other hand, the parameter t0 can be very easily measured in tracer experiments as it is the time of the first nonzero concentration measured in the stream.
It should be noted that the values of the dispersion coefficient for the GumbelD and the GEVD methods were very similar. In fact, the Gumbel distribution, having the shape parameter ξ equal to zero (see (18)), is a special case of the GEV distribution. The results showed that the nonzero value of the shape parameter allows better fit of the approximation curve to the measured data, both for standard and distorted hydraulic conditions in the stream (flood, dead zones, etc.).
As it can be seen in Table 7, there are no noticeable or significant dependencies between the hydraulic parameters of examined streams and the particular method accuracy expressed by the NRMSE parameter. On the other hand, some correlations can be observed between the longitudinal dispersion coefficients Dx, Dx,G, and Dx,GEV (Gauss, Gumbel, and GEV methods) on stream distance from injection point, stream velocity, and stream width (see Table 8). The parameter ξ, used in the GEV method (see (22)) has only a weak correlation to the stream hydraulic parameters, but graph analysis revealed that higher values of this parameter correspond with the long “tail” of the time-concentration curve. Hence, it is possible to assume that this parameter characterizes the extent of dead zones in the stream. Very weak (almost none) correlations were found between the dispersion parameters obtained using the LogNorm approximation method, which can be regarded as a serious disadvantage of this method (difficulty with parameters prediction).

5.2. Field Experiments

The results of the statistical test performed on field data are presented in Table 10, Table 11, Table 12 and Table 13 The results were generally very similar to those of the literature data and confirmed the above conclusions.
Higher accuracy of the LogNormD and GEVD approximation methods was observed, i.e., a better fit of the modelled data to the experimental data. In addition, the average precision of all examined methods increased almost by 100%. This statement also applies to the minimal and maximal values, as well as for the standard deviation (see Table 10 and Table 11).
The exceeding values test also documents the higher precision of proposed approximation methods using the field experiment data. Table 12 shows that no experiments were less precise than 60% and 50% of the GaussD method using the GumbelD, the LogNormD, and GEVD methods, respectively. The correlation results between the dispersion coefficients were almost identical to the cases using literature data (Table 13). The enhanced statistical analysis (similar as in Table 7, Table 8 and Table 9) was not possible for the field experiments carried in Slovakia, due to the low variety of hydraulic parameters.

6. Application of the Proposed Approximation Method to Inverse Tasks

The initial idea to explore different approximation methods for use in solute transport tasks originated within an ongoing research grant focused on localization of pollution sources on rivers (VEGA 1/0805/16, see Funding), i.e., solving the inverse task. One of the methods examined for solving such a task was the so-called “brute force method”, which is based on simulation of all possible scenarios (contaminant distribution in time, distances, etc.) and comparison of the results with measured data. The best fit between the simulated and measured data set up the position of the pollution source. Such an approach requires a large number of simulations to be performed and the results are highly dependent on the simulation precision and adaptation of local conditions.
Use of the GaussD method (5) led to unsatisfactory results (overall precision about 30–50%), so the idea was to use much more precise simulation methods, which are fast and able to describe the pollutant breakthrough curve more precisely, in conditions of real streams with transient (dead) zones. The solute (tracer) source localization experiment was performed on real data from the Malina stream experiment. Both numerical experiments were performed using our own simulation software LOPS (LOcalisation of Pollution Sources) with the time step of 10 minutes and space (longitudinal) step of 50 m using two different simulation methods: GaussD and the GumbelD methods. The longitudinal dispersion coefficients were determined by the results of the tracer experiment and the measured breakthrough curve was intentionally shifted in time by 3600 seconds. In the first case (GaussD method, (5)), the best fit between the measured and modelled data was achieved at a distance of 1750 m, whereas the exact distance was 1,435 m. The inaccuracy was quite large (22%) and mainly due to the inaccuracy of the model results (using GaussD method) compared with the field data.
The same conditions, but different approximation method—GumbelD (Equation (11)) were used in the second case. The fit between the modelled and measured data was much better than in the previous case—the difference sum was 0.229 in the first case, in the second case it was 0.024, i.e., approximately ten times less. Also, the solute source distance was determined almost exactly at 1450 m.

7. Conclusions

The application of alternative forms of the 1D analytical solution for the hydrodynamic dispersion in rivers with dead zones was explored in this study. Besides the “classical” Gauss distribution, Gumbel, log-normal, and generalized extreme value distributions were considered to evaluate their ability to approximate the dispersion process from an instantaneous injection. For this evaluation, both literature data and experimental data collected in three small streams in Slovakia were applied.
For both datasets, the results of this study show that all three approximation methods provided a higher accuracy than the Gaussian–based method. The highest accuracy was obtained using the LogNormD and the GEVD methods, having very similar precision. In particular, the LogNormD and GEVD methods were less accurate than the GaussD approximation method in only 0.7% and 2.1% of the considered cases, respectively. The GumbelD method was the least accurate when compared with the LogNormD and GEVD methods, but still more accurate than the GaussD method by around 53%.
A partial comparison of the presented methods with the dead zone model was presented in [40], where the GaussD and the GumbelD were compared with the aggregated dead zone (ADZ) model, represented by the OTIS-P model [8]. As referred in this paper, the performance of both methods (GumbelD and ADZ) was comparable, so it can be assumed that it will be similar also in the case of the LogNormD and GEVD methods. In our view, the relative ease of use of the methods used is clearly in favor of the proposed methods, for their use is sufficiently simple in a commonly available spreadsheet editor, no specialized SW is needed.
All the three examined approximation methods are suitable, especially for streams with significant influence of dead zones, which cause a strong asymmetry in the shape of the concentration distribution curves over time. These approximation methods can be shortly characterized as follows:
  • The Gumbel method is very easy to compute, but is less accurate because this method does not allow justification of the time concentration curve shape expressing the extent of dead zones along the examining stream reach;
  • The LogNormD approximation method is more precise than Gumbel and Gauss methods, but its parameters are not directly determinable and there is considerable uncertainty in their determination (there is a simple estimation of the t0 parameter, but it has difficult predictable dispersion and location parameters (µ, Dx,LN), which have no direct physical meaning);
  • The GEVD approximation method is the most precise method, the dispersion coefficient (Dx,GEV) has similar meaning and way of determination as the dispersion coefficients Dx and Dx,G. It also has relatively predictable shape coefficient ξ (higher values correspond with longer “tails” of the time concentration curves, i.e., it expresses the extent of the dead zones). The correlation between the dispersion coefficient used in this method (Dx,GEV) and the dispersion coefficient Dx, used in Equation (5) can be a subject of further research, with the aim of determination of the Dx,GEV value, as well as the shape coefficient (ξ) value.
The application of the Gaussian distribution (GaussD method) in cases of streams with large a extent of dead zones is problematic and its use leads to significant differences between measured and simulated values. More detailed and precise methods for the dispersion simulation in streams with dead zones can be found in the literature [10], but these methods are quite difficult to apply as they are based on a solution of partial differential equations and require a high number of parameters.
An added advantage of the novel approach herein proposed is its simplicity and low computational effort. This allows its application for fast and simple simulations of contaminant dispersion in rivers and streams using the discrete advection–diffusion equation (ADE) technique, as described by [41], and applying the superposition principle, which is a common feature of all analytical solutions. These features have been used successfully to solve inverse tasks, where a large number of simulations are required. As documented in the discussion, such approach has significantly increased the accuracy of the inverse task solution.

Author Contributions

Conceptualization, M.S.; data curation, Y.V.; investigation, M.S. and Y.V.; methodology, C.G.; software, M.S.; validation, C.G.; writing—original draft, M.S.; writing—review and editing, M.S., Y.V. and C.G.

Funding

This paper was prepared with the support of the Scientific Grant Agency VEGA within the scientific project Nr. VEGA 1/0805/16 “Localisation of accidental point sources of pollution in watercourses based on-line monitoring data”, project of the H2020 program SYSTEM, grant agreement nr. 787128. It was also supported by the Scientific Grant Agency APVV within the implementation of the project no. APVV-18-0205 “Management of crisis situations in water supply with respect to climate change”. Carlo Gualtieri acknowledges the financial support from the Erasmus Plus (Erasmus+) Program of the University of Napoli Federico II for the period 2017–2018.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Weitbrecht, V. Influence of Dead-Water Zones on the Dispersive Mass-Transport in Rivers, Dissertationsreihe am Institut fuer Hydromechanik der Universitat Karlsruhe ed.; Universitatsverlag Karlsruhe: Karlsruhe, Germany, 2004; Volume 2004/1. [Google Scholar]
  2. Gualtieri, C. Numerical simulation of flow patterns and mass exchange processes in dead zones. In Proceedings of the iEMSs Fourth Biennial Meeting: International Congress on Environmental Modelling and Software (iEMSs 2008), Barcelona, Spain, 1 July 2008. [Google Scholar]
  3. Valentine, E.M.; Wood, I.R. Longitudinal dispersion with dead zones. J. Hydraul. Div. ASCE 1977, 103, 975–990. [Google Scholar]
  4. Valentine, E.; Wood, I. Experiments in Longitudinal Dispersion with Dead Zones. J. Hydraul. Div. ASCE 1979, 105, 999–1016. [Google Scholar]
  5. Tonina, D. Surface water and streambed sediment interaction: The hyporheic exchange. In Fluid Mechanics of Environmental Interfaces; Gualtieri, C., Mihailovic, D., Eds.; CRC Press: Boca Raton, FL, USA, 2012; pp. 255–294. [Google Scholar]
  6. Bencala, K.; Walters, R.A. Simulation of solute transport in a mountain pool-and-riffle stream: A transient storage model. Water Resour. Res. 1983, 19, 718–724. [Google Scholar] [CrossRef]
  7. Beer, T.; Young, P.C. Longitudinal dispersion in natural streams. J. Environ. Eng. ASCE 1983, 109, 1049–1067. [Google Scholar] [CrossRef]
  8. Runkel, R.L. One-Dimensional Transport with Inflow and Storage (OTIS): A Solute Transport Model for Streams and Rivers, Water-Resources Investigations Report 98–4018 ed.; Geological Survey: Denver, CO, USA, 1998. [Google Scholar]
  9. de Smedt, F.; Brevis, W.; Debels, P. Analytical solution for solute transport resulting from instantaneous injection in streams with transient storage. J. Hydrol. 2005, 315, 25–39. [Google Scholar] [CrossRef]
  10. de Smedt, F. Analytical solutions for transport of decaying solutes in rivers with transient storage. J. Hydrol. 2006, 330, 672–680. [Google Scholar] [CrossRef]
  11. de Smedt, F. Analytical solution and analysis of solute transport in rivers affected by diffusive transfer in the hyporheic zone. J. Hydrol. 2007, 339, 29–38. [Google Scholar] [CrossRef]
  12. Deng, Z.-Q.; Jung, H.-S.; Ghimire, B. Effect of channel size on solute residence time distributions in rivers. Adv. Water Resour. 2010, 33, 1118–1127. [Google Scholar] [CrossRef]
  13. Shucksmith, J.D.; Boxall, J.B.; Guymer, I. Effects of emergent and submerged natural vegetation on longitudinal mixing in open channel flow. Water Resour. Res. 2010, 46. [Google Scholar] [CrossRef]
  14. Shucksmith, J.D.; Boxall, J.B.; Guymer, I. Determining longitudinal dispersion coefficients for submerged vegetated flow. Water Resour. Res. 2011, 47. [Google Scholar] [CrossRef]
  15. Murphy, E.; Ghisalberti, M.; Nepf, H. Model and laboratory study of dispersion in flows with submerged vegetation. Water Resour. Res. 2007, 43. [Google Scholar] [CrossRef]
  16. Davis, P.; Atkinson, T. Longitudinal dispersion in natural channels: 3. An aggregated dead zone model applied to the River Severn, UK. Hydrol. Earth Syst. Sci. 2000, 4, 373–381. [Google Scholar] [CrossRef]
  17. Czernuszenko, W.; Rowinski, P.; Sukhodolov, A. Experimental and numerical validation of the dead-zone model for longitudinal dispersion in rivers. J. Hydraul. Res. 1998, 36, 269–280. [Google Scholar] [CrossRef]
  18. González-Pinzón, R.; Haggerty, R.; Dentz, M. Scaling and predicting solute transport processes in streams. Water Resour. Res. 2013, 49, 4071–4088. [Google Scholar] [CrossRef] [Green Version]
  19. Socolofsky, S.A.; Jirka, G.H. Special Topics in Mixing and Transport Processes in the Environment, 5th ed.; Texas A&M University: College Station, TX, USA, 2005; p. 172. [Google Scholar]
  20. Rutherford, J. River Mixing; John Wiley & Sons: Chichester, UK, 1994; p. 348. [Google Scholar]
  21. Shucksmith, J.; Boxall, J.; Guymer, I. Importance of advective zone in longitudinal mixing experiments. Acta Geophys. 2007, 55, 95–103. [Google Scholar] [CrossRef]
  22. Chanson, H. Environmental Hydraulics of Open Channel Flows; Jordan, H., Ed.; Elsevier Butterworth-Heinemann Linacre House: Oxford OX2 8DP 200 Wheeler Road, Burlington, NJ, USA, 2004; p. 430. [Google Scholar]
  23. Nordin, C.; Troutman, B. Longitudinal dispersion in rivers: The persistence of skewness in observed data. Water Resour. Res. 1980, 16, 123–128. [Google Scholar] [CrossRef]
  24. Czernuszenko, W.; Rowinski, P. Properties of the dead-zone model of longitudinal dispersion in rivers. J. Hydraul. Res. 1997, 35, 491–504. [Google Scholar] [CrossRef]
  25. Martin, J.L.; McCutcheon, S.C. Hydrodynamics and Transport for Water Quality Modeling; CRC Press, Inc.: Boca Raton, FL, USA, 1998; p. 816. [Google Scholar]
  26. Fischer, H.B.; List, E.; Koh, R.; Imberger, J.; Brooks, N. Mixing in Inland and Coastal; Academic Press: New York, NY, USA, 1979. [Google Scholar]
  27. Cvetkovic, V. The tempered one-sided stable density: A universal model for hydrological transport? Environ. Res. Lett. 2011, 6, 034008. [Google Scholar] [CrossRef]
  28. Koudou, A.E.; Ley, C. Characterizations of GIG laws: A survey. Probab. Surv. 2014, 11, 161–176. [Google Scholar] [CrossRef]
  29. Folks, J.L. Inverse Gaussian Distribution. In Encyclopedia of Biostatistics; Colton, P.A.A.T., Ed.; John Wiley & Sons, Ltd.: Hoboken, NJ, USA, 2005. [Google Scholar]
  30. Gooseff, M.N.; Benson, D.A.; Briggs, M.A.; Weaver, M.; Wollheim, W.; Peterson, B.; Hopkinson, C.S. Residence time distributions in surface transient storage zones in streams: Estimation via signal deconvolution. Water Resour. Res. 2011, 47. [Google Scholar] [CrossRef] [Green Version]
  31. Loaiciga, H.A.; Leipnik, R.B. Analysis of extreme hydrologic events with Gumbel distributions: Marginal and additive cases. Stoch. Environ. Res.Risk Assess. 1999, 13, 251–259. [Google Scholar] [CrossRef]
  32. Aristizabal, R.J. Estimating the Parameters of the Three-Parameter Lognormal Distribution, FIU Electronic Theses and Dissertations ed.; Florida International University: Miami, FL, USA, 2012. [Google Scholar]
  33. Koutsoyiannis, D. Probability and Statistics for Geophysical Processes; National Technical University of Athens: Athens, Greece, 2008. [Google Scholar]
  34. Nordin, C.F.; Sabol, G. Empirical Data on Longitudinal Dispersion, Water-Resources Investigations Report, 74-20th ed.; Geological Survey: Reston, VA, USA, 1974; p. 332. [Google Scholar]
  35. Limerinos, J.T. Determination of the Manning coefficient for measured bed roughness in natural channels. Water Supply Paper 1970, 1898-B, 53. [Google Scholar]
  36. Gualtieri. RANS-based simulation of transverse turbulent mixing in a 2D geometry. Environ. Fluid Mech. 2010, 10, 137–156. [Google Scholar] [CrossRef]
  37. Shields, D.; Coulton, K.G.; Nepf, H. Representation of vegetation in two-dimensional hydrodynamic models. J. Hydraul. Eng. 2017, 143, 02517002. [Google Scholar] [CrossRef]
  38. Bouchez, J.; Lajeunesse, E.; Gaillardet, J.; France-Lanord, C.; Dutra-Maia, P.; Maurice, L. Turbulent mixing in the Amazon River: The isotopic memory of confluences. Earth Planet. Sci. Lett. 2011, 311, 448–450. [Google Scholar] [CrossRef]
  39. Ranatunga, T.; Tong, S.T.; Yang, J. An approach to measure parameter sensitivity in watershed hydrological modelling. Hydrol. Sci. J. 2017, 62, 76–92. [Google Scholar] [CrossRef]
  40. Sokáč, M.; Velísková, Y.; Gualtieri, C. An approximate method for 1-D simulation of pollution transport in streams with dead zones. J. Hydrol. Hydromech. 2018, 66, 437–447. [Google Scholar] [CrossRef]
  41. Guymer, I. A National Database of Travel Time, Dispersion and Methodologies for the Protection of River Abstractions; R&D Technical Report P346; Environment Agency: Bristol, UK, 2002. [Google Scholar]
Figure 1. Graphical evaluation of approximation methods—case 28, Red River, March 13, 1972, site 1, distance 5 miles (8.045 km), see [34].
Figure 1. Graphical evaluation of approximation methods—case 28, Red River, March 13, 1972, site 1, distance 5 miles (8.045 km), see [34].
Water 11 02145 g001
Figure 2. Graphical evaluation of approximation methods—case 15, Salt Creek, September 6, 1972, site 6, distance 37.62 miles (60.5 km), see [34].
Figure 2. Graphical evaluation of approximation methods—case 15, Salt Creek, September 6, 1972, site 6, distance 37.62 miles (60.5 km), see [34].
Water 11 02145 g002
Figure 3. Graphical evaluation of approximation methods—case 20, Comite River, November 4, 1968, site 5, distance 49 miles (78.8 km), see [34].
Figure 3. Graphical evaluation of approximation methods—case 20, Comite River, November 4, 1968, site 5, distance 49 miles (78.8 km), see [34].
Water 11 02145 g003
Figure 4. Emerged and submerged vegetation in the Malina stream, West Slovakia, near the village of Zohor (GPS 48.334830, 16.967463).
Figure 4. Emerged and submerged vegetation in the Malina stream, West Slovakia, near the village of Zohor (GPS 48.334830, 16.967463).
Water 11 02145 g004
Figure 5. Graphical evaluation of approximation methods—Malina Stream, August 3, 2016, site 6, distance 1.435 km.
Figure 5. Graphical evaluation of approximation methods—Malina Stream, August 3, 2016, site 6, distance 1.435 km.
Water 11 02145 g005
Figure 6. Graphical evaluation of approximation methods—Šúrsky Stream, March 20, 2014, experiment 12-II, distance 0.406 km.
Figure 6. Graphical evaluation of approximation methods—Šúrsky Stream, March 20, 2014, experiment 12-II, distance 0.406 km.
Water 11 02145 g006
Table 1. Main characteristics of the considered rivers/streams [34].
Table 1. Main characteristics of the considered rivers/streams [34].
Nr.Case Nr. *River/Stream NameDateMin. DistanceMax. DistanceMin. DischargeMax. DischargeMin. VelocityMax. Velocity
(km)(km)(m3.s1)(m3.s1)(m.s1)(m.s1)
11Anteniam Creek27 May 19692.5729.611.161.780.120.22
22Anteniam Creek24 March 19702.5766.695.1012.180.350.51
37Monocacy River7 June 196810.3034.2715.1518.550.370.38
49Monocacy River25 September 19687.4833.793.063.200.080.11
510Conococheague Creek6 May 19694.4219.876.826.940.350.41
612Conococheague Creek30 April 19704.4233.8729.4530.580.620.87
713Chattahoochee River20 April 197116.0176.59107.60179.530.730.97
814Chattahoochee River4 May 197110.49104.60140.17140.170.751.09
915Salt Creek River6 September 19729.2551.812.464.080.290.34
1019Bayou Anacoco River15 June 196911.4237.972.012.710.150.16
1120Comite River4 November 19686.7678.840.761.130.080.15
1221Bayou Bartholomew25 June 19713.22117.464.118.100.120.15
1322Amite River21 October 19689.98148.035.669.910.180.33
1424Tangipahoa River24 April 19698.2193.975.7818.690.220.34
1525Tangipahoa River15 September 19698.2193.973.4510.850.160.27
1626Red River7 April 19715.74193.08230.22249.470.570.63
1727Red River14 April 19718.05159.29139.60184.910.460.64
1828Red River13 March 19728.05159.29187.46251.740.430.52
1929Red River12 June 197212.07199.52107.60166.500.450.51
2030Sabine River9 September 19697.8879.81110.72130.260.400.80
2131Sabine River4 February 197217.22209.65311.49433.250.570.84
2232Sabine River16 April 197222.2085.920.481.330.090.10
2333Sabine River16 April 197217.06121.001.026.910.070.11
2434Sabine River16 April 197221.08113.436.177.080.120.15
2535Mississippi River15 September 196535.40204.343398.026796.040.620.81
2636Mississippi River11 March 196854.71294.452633.472973.270.951.10
2737Mississippi River7 August 196854.71294.456824.366824.361.361.48
2838Wind/Righorn River21 March 19719.17181.3354.9368.810.830.95
2939Wind/Righorn River29 June 19719.17181.33215.21254.851.511.86
Min. value 2.5719.870.481.130.070.10
Max. value 54.71294.456824.366824.361.511.86
* means the numbering according to [34]
Table 2. Main characteristics of the streams used for tracer experiments in Slovakia.
Table 2. Main characteristics of the streams used for tracer experiments in Slovakia.
Nr.Stream/
Channel Name
DateMin. DistanceMax. DistanceMin. DischargeMax. DischargeMin. VelocityMax. Velocity
(km)(km)(m3.s−1)(m3.s−1)(m.s−1)(m.s−1)
1Malá Nitra3–4 July 20120.7851.3400.2300.2350.1310.163
2Šúrsky stream 19–20 March 20140.3970.5060.3800.4300.2450.282
3Malina3 August 20161.4351.4350.4080.4080.0760.076
Table 3. Comparison of approximation methods’ accuracy.
Table 3. Comparison of approximation methods’ accuracy.
Statistical ParameterGaussDGumbelDLogNormDGEVD
%%%%
NRMSE—Average6.403.402.602.44
NRMSE—Min1.620.710.690.70
NRMSE—Max15.5511.006.806.42
Std. dev.2.491.781.281.15
Table 4. Relative comparison of approximation methods’ accuracy (relative to Gaussian distribution (GaussD) = 100%).
Table 4. Relative comparison of approximation methods’ accuracy (relative to Gaussian distribution (GaussD) = 100%).
Statistical Parameter GaussDGumbelDLogNormDGEVD
%%%%
Relative NRMSE—Average10053.140.638.2
Relative NRMSE—Min10040.642.542.8
Relative NRMSE—Max10070.743.741.3
Std. Dev.10067.651.346.2
Table 5. Exceeding values.
Table 5. Exceeding values.
LimitGumbelDLogNormDGEVDGumbelDLogNormDGEVD
---%%%
0%0000%0% 0%
10%0100.0%0.7%0.0%
20%215131.4%10.6%9.2%
40%38637627.0%44.7%53.9%
50%709710949.6%68.8%77.3%
60%9812312569.5%87.2%88.7%
80%12013713785.1%97.2%97.2%
100%12614013889.4%99.3%97.9%
Table 6. Correlation between the dispersion coefficients.
Table 6. Correlation between the dispersion coefficients.
MethodGaussDGumbelDLogNormDGEVD
GaussD-0.99880.21310.9983
GumbelD0.9988-0.18970.9967
LogNormD0.21310.1897-0.2225
GEVD0.99830.99670.2225-
Table 7. Correlation (rxy) between the NRMSE and hydraulic parameters of examined streams.
Table 7. Correlation (rxy) between the NRMSE and hydraulic parameters of examined streams.
NRMSE of MethodQ–DischargeX–DistanceV–VelocityW–River WidthD–Mean DepthS0–SlopeU*–Shear VelocityA–Unitless Disp. Coef.
GaussD0.2770.1730.0880.2120.049–0.280–0.1580.401
GumbelD0.2320.144–0.0650.1260.038–0.201–0.2300.378
LogNormD0.2480.272–0.0280.1820.096–0.202–0.1540.271
GEVD0.2500.309–0.0230.2410.136–0.228–0.1660.285
Note: the unitless dispersion coefficient was defined as a = D x d . u * (s−1)
Table 8. Correlation (rxy) between the approximation method parameters and hydraulic parameters of examined streams.
Table 8. Correlation (rxy) between the approximation method parameters and hydraulic parameters of examined streams.
ParameterEquationQ—DischargeX—DischargeV—VelocityW—River WidthD—Mean DepthS0—SlopeU*—Shear VelocityA—Unitless Disp. Coef.
Dx(5)0.2480.4160.6000.5930.357–0.1210.2240.269
Dx,G(10)0.2310.4160.6040.6090.368–0.1090.2400.239
Dx.LN(16)0.227–0.0960.0610.065–0.047–0.189–0.1050.348
k=t0/ttmax (17)0.1180.215–0.0160.060–0.057–0.382–0.3420.166
Dx,GEV(22)0.2320.4090.5940.5680.341–0.1210.2120.305
ξ(22)0.261–0.0190.0770.083–0.030–0.214–0.1460.371
Table 9. Ranges of dispersion parameters for each approximation method.
Table 9. Ranges of dispersion parameters for each approximation method.
Dispersion ParameterEquationAverageMin. ValueMax. ValueStd. Deviation
Dx(5)92.270.973353.66299.88
Dx,G(10)151.661.655733.00507.40
Dx.LN(16)0.4130.1001.1400.181
k = t0/ttmax (17)0.8350.1340.9750.102
Dx,GEV(22)163.481.686314.14560.83
ξ(22)0.030–0.3540.8600.182
Table 10. Comparison of approximation methods’ accuracy.
Table 10. Comparison of approximation methods’ accuracy.
Statistical ParameterGaussDGumbelDLogNormDGEVD
%%%%
NRMSE—Average9.684.813.523.01
NRMSE—Min5.321.721.391.08
NRMSE—Max15.7011.5910.7411.00
Std. Dev.3.192.692.502.68
Table 11. Relative comparison of approximation methods’ accuracy (relative to Gauss = 100%).
Table 11. Relative comparison of approximation methods’ accuracy (relative to Gauss = 100%).
Statistical ParameterGaussDGumbelDLogNormDGEVD
%%%%
Relative NRMSE—Average10049.736.431.1
Relative NRMSE—Min10032.226.220.3
Relative NRMSE—Max10073.868.470.1
Std. Dev.10084.578.483.9
Table 12. Exceeding values.
Table 12. Exceeding values.
LimitGumbelDLogNormDGEVDGumbelDLogNormDGEVD
---%%%
0%161616100.0%100.0%100.0%
10%1410587.5%62.5%31.3%
20%91156.3%6.3%6.3%
40%31118.8%6.3%6.3%
50%20012.5%0.0%0.0%
60%0000.0%0.0%0.0%
80%0000.0%0.0%0.0%
100%0000.0%0.0%0.0%
Table 13. Correlation between the dispersion coefficients.
Table 13. Correlation between the dispersion coefficients.
MethodGaussDGumbelDLogNormDGEVD
GaussD-0.99650.26750.9993
GumbelD0.9965-0.22880.9981
LogNormD0.26750.2288-0.2687
GEVD0.99930.99810.2687-

Share and Cite

MDPI and ACS Style

Sokáč, M.; Velísková, Y.; Gualtieri, C. Application of Asymmetrical Statistical Distributions for 1D Simulation of Solute Transport in Streams. Water 2019, 11, 2145. https://doi.org/10.3390/w11102145

AMA Style

Sokáč M, Velísková Y, Gualtieri C. Application of Asymmetrical Statistical Distributions for 1D Simulation of Solute Transport in Streams. Water. 2019; 11(10):2145. https://doi.org/10.3390/w11102145

Chicago/Turabian Style

Sokáč, Marek, Yvetta Velísková, and Carlo Gualtieri. 2019. "Application of Asymmetrical Statistical Distributions for 1D Simulation of Solute Transport in Streams" Water 11, no. 10: 2145. https://doi.org/10.3390/w11102145

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop