Applications of Two-Dimensional Spatial Routing Procedure for Estimating Dispersion Coefficients in Open Channel Flows

In this study, the performance of two routing procedures were evaluated to estimate the two-dimensional dispersion coefficients. The two-dimensional Stream-Tube Routing Procedure (2D ST-RP) has been widely used to obtain the dispersion coefficients from measured concentration-time curves under the frozen cloud assumption. Meanwhile, the Spatial Routing Procedure (2D S-RP) employs the spatial distributions of concentration to estimate the dispersion coefficients. The performance of the two routing methods were evaluated in aspect of the validity of the frozen cloud assumption and the applicability in the non-Fickian mixing. From the estimation results of dispersion coefficients, the results by the 2D ST-RP included errors due to skewed concentration-time curves which were created by violating the frozen cloud assumption. On the other hand, the 2D S-RP provides accurate dispersion coefficients in the same condition. The estimated results of dispersion coefficients in the meandering channel show that both the 2D ST-RP and the 2D S-RP contained errors due to the non-Fickian mixing properties of the test case. Even with the discrepancies, the 2D S-RP presented more appropriate spatial variabilities along the meander cycle than the results by the 2D ST-RP.


Introduction
Pollutants in natural streams are transported by the temporal and spatial variations of shear flows. The shear flow dispersion is dominantly governed by the vertical deviations of longitudinal and transverse velocities from the depth-averaged velocity because the scale of shear flow is much larger than the turbulent fluctuations. These mixing phenomena in shear flow can be modeled using the two-dimensional (2D) advection-dispersion equation, which was derived from the depth-averaging process of three-dimensional advection diffusion [1]. The 2D advection-dispersion equation follows the Fick's law, which defines solute mass flux as the product of the concentration gradient and dispersion coefficient. The dispersion coefficient in the Fickian model has the important role of controlling not only the numerical dispersion but also the spatial dispersion of pollutants.
The calculation method of the dispersion coefficient can be classified into velocitybased and concentration-based methods. In the velocity-based method, dispersion coefficients are calculated using vertical profiles of longitudinal and transverse velocities [1]. Shen et al. [2] and Kim [3] calculated the longitudinal dispersion coefficients from the velocity measurements using an Acoustic Doppler Current Profiler (ADCP). Furthermore, Lee and Seo [4] analyzed the spatial properties of solute mixing in meandering channels from the calculation results of dispersion coefficients using velocity profiles. As the velocity-based method includes the effects of spatially varying velocity distributions, the method has advantages in reproducing shear dispersion processes using numerical models in meandering channels [5]. However, the velocity-based method tends to underestimate dispersion coefficients compared to the concentration-based methods because the method neglects the effects of turbulent diffusion and storage zones [6,7]. Furthermore, it is difficult to obtain sufficient data from velocity structures because it requires significant time and costs. Thus, in natural streams, the bulk shear dispersion coefficients have been obtained from the concentration-based methods.
The concentration-based methods have been widely used to estimate the dispersion coefficients due to the aforementioned limitations of the velocity-based method. Previous tracer studies, which were conducted in laboratory channels [8,9] and natural streams [10,11], have been reported to evaluate shear dispersion in various flow conditions using temporal variations of concentration. From the concentration-time curves, dispersion coefficients have been estimated using the moment method or the routing procedure [12]. The moment method has been popularly used to obtain the transverse dispersion coefficient from the lateral variation of concentration-time curves [13][14][15][16]. However, the moment method is inappropriate to obtain the longitudinal dispersion coefficient for twodimensional mixing problems. Thus, the two-dimensional stream-tube routing procedure was suggested to obtain the longitudinal and the transverse dispersion coefficients simultaneously [17]. These concentration-based methods were developed under the frozen cloud assumption, which considers time invariant shear dispersion while a tracer cloud is crossing a measuring cross-section [18]. The frozen cloud assumption is essential for estimating the dispersion coefficients using concentration-time curves measured at a fixed point; however, the estimated dispersion coefficients might include errors [19]. The limitations of the field observations, in which probes are installed at a fixed point to measure concentration-time curves, are being overcome by the research of image analysis that converts a tracer cloud taken from images by unmanned aerial vehicles into a concentration field [20][21][22]. It could be possible to obtain the dispersion coefficients without the frozen cloud assumption if spatial distributions of concentration are available to use.
In this study, performance of the two-dimensional Spatial Routing Procedure (2D S-RP) and the two-dimensional Stream-Tube Routing Procedure (2D ST-RP) was evaluated by comparisons of the estimated dispersion coefficients. As aforementioned, the 2D ST-RP employs the concentration-time curves under the assumption of frozen cloud. The 2D S-RP uses the spatial distributions of concentration to estimate the dispersion coefficients and the frozen cloud assumption is not required. From these two routing procedures, applicability was evaluated against two test cases: (1) idealized solute mixing problem; (2) 2D tracer test results conducted in a meandering channel. From the comparisons of estimated dispersion coefficients, the validity of the frozen cloud assumption was evaluated. In addition, the applicability of two routing procedures was discussed for tracer mixing showing non-Fickian mixing properties.

Velocity-Based Method
Following the Fick's law, the solute transport in rivers was analyzed using the depthaveraged advection-dispersion equation under the assumption of large width to depth ratio. The depth-averaged advection-dispersion Equation is given as where C is the depth-averaged concentration; h is the water depth; t is the time; i u is the depth-averaged velocity components in i-direction;  ij is the turbulent diffusion coefficient;    j j j u u u ; j u are the time-averaged velocity components in j-direction;    c c C ; c is the concentration over flow depth. The second term in the right-hand side of Equation (1) was defined by Fischer et al. [1], who extended the one-dimensional study of Taylor [23] to the 2D mixing problems as where ij D is the dispersion coefficients in tensor form. Equation (2), which is called the Fickian dispersion model, has been widely used to simulate the 2D pollutant mixing [24][25][26]. The dispersion coefficients shown in Equation (2) can be calculated using the vertical profiles of velocity according to the theoretical derivation as where  v is the vertical turbulent diffusion coefficient. In open channel flows,  v follows the parabolic distribution as Dispersion coefficients in tensor form are efficient to simulate 2D solute mixing problems in the meandering channel [4]. However, Equation (3) is difficult to adopt in 2D models due to lack of vertical profiles of velocity. Thus, Alavian [27] suggested the coordinate transform method to define the dispersion coefficients into the tensor format. In the study by Alavian, dispersion tensor was derived using local velocity components in a fixed Cartesian coordinate and dispersion coefficients defined in the flow direction. Each component of a dispersion tensor is where L D and T D are the longitudinal and transverse dispersion coefficients, respectively; u and v are the depth-averaged velocity component in x-and y-direction, respectively. When the flow direction coincides with x-or y-direction, dispersion by xy D and yx D is neglected and solute mixing is proceeded by xx D and yy D .

Concentration-Based Method
Routing Procedure Using Concentration-Time Curves Methods for obtaining dispersion coefficients from the tracer test results can be classified into moment-based methods and routing procedures [17]. The moment-based methods have been used to estimate the transverse dispersion coefficient using continuous and transient tracer tests. The estimated dispersion coefficients using the moment method may include errors when the tracer clouds show the asymmetric distributions because values of the second moment do not provide proper relations following the Fick's law [28]. Furthermore, the moment method is not appropriate in obtaining the longitudinal dispersion coefficient. Thus, since the study by Fischer [29], several studies [17,30,31] have proposed routing procedures to compensate for the disadvantages of the moment method.
The routing procedure was first suggested by Fischer [29] for calculating the onedimensional longitudinal dispersion coefficient. Following the study by Fischer, the routing procedure can be expanded to the two-dimensional case as follows [30]: where   1 , , C s n t is the measured upstream concentration-time curves;   2 , , C s n t is the predicted downstream concentration-time curves; 1 t and 2 t are the average passage time of the tracer cloud; U is the cross-sectional average velocity. Using the twodimensional routing procedure, the longitudinal and transverse dispersion coefficients can be obtained simultaneously using the upstream and downstream concentration-time curves. However, Equation (8) is inappropriate to use for a river that includes complex channel geometries. Advancing from Equation (8), Baek and Seo [17] proposed the 2D ST-RP for adopting the routing procedure in irregularly varied channel width. The 2D ST-RP was derived by adopting the stream-tube concept as follows: ,  is the normalized shape factor defined as    ; q is the cumulative discharge;  is the normalized cumulative discharge. From Equation (9), L D and T D can be obtained using the multiple regression algorithm. In this study, the non-linear regression algorithm of MATLAB software (Mathworks, Natick, USA), which finds the dispersion coefficients to predict best-fitted concentration curves, was adopted to estimate the dispersion coefficients.

Routing Procedure Using Spatial Distributions of Concentration
Usually, tracer tests have been conducted by measuring concentrations of tracer material at a fixed point, and the measured data provide the temporal variations of concentration. Therefore, the routing procedure using the concentration-time data requires the frozen cloud assumption in which the concentration curve does not experience change of dispersion during passage through the fixed sensors. However, the routing method under the frozen cloud assumption needs to be validated to determine whether the assumption causes errors to the estimated dispersion coefficients. Thus, in this study, the routing procedure using the spatial concentration distributions was adopted for estimating the dispersion coefficients. In the case of using the routing method using the spatial concentration distributions measured at a specific moment, it is possible to calculate the dispersion coefficients without incorporating the frozen cloud assumption.
From the analytic solution of the two-dimensional advection-dispersion equation, the two-dimensional spatial routing procedure (2D S-RP) can be derived as follows: is the predicted downstream concentration field. L D and T D in Equation (10) were estimated as the non-linear regression algorithm, which is the same algorithm used for obtaining the dispersion coefficients using the 2D ST-RP method. Figure 1 shows the conceptual diagram of the two routing procedures, which are the 2D ST-RP and the 2D S-RP. The main difficulty of the 2D S-RP is that it is difficult to obtain the spatial concentration distributions in a field study. Therefore, in this study, 2D simulation results using the PDM-2D were employed for application of the 2D S-RP.

Model Descriptions
Park and Seo [32] developed the 2D Particle Dispersion Model (PDM-2D) to simulate two-dimensional pollutant mixing in rivers without input of dispersion coefficients based on the study by Fischer [29]. In the study by Fischer, the solute mixing was analyzed using step-by-step computations according to the physical interpretations of the shear dispersion instead of solving the advection-dispersion equation using numerical schemes. The PDM-2D calculates the 2D pollutant mixing using the two sequential calculation steps: (1) horizontal translation and (2) vertical mixing. Figure 2 describes the aforementioned computation procedures of the PDM-2D. As described in Figure 2, in the horizontal translation step of PDM-2D, the individual particle position was calculated as follows: where k i x is the k th particle position in i-direction;  k i x is the particle position after shear advection; sa t is the time after shear advection; z is the vertical position of each particle; t is the time step; R is the random number following the Gaussian distribution;  h is the horizontal turbulent diffusion coefficient in the assumption of isotropic turbulence. The second term of the right-hand side in Equation (11) indicates the deterministic transport by vertically varied velocity profiles, and the third term is the transport of random fluctuations by the turbulent diffusion. The solute column is stretched on the horizontal planes by the vertical profiles of longitudinal ( s u ) and transverse velocities ( n u ), as shown in Figure 2. The vertical velocity profiles were generated using theoretical and empirical formulas, which can be reproduced using the depth-averaged flow fields. In the longitudinal direction, the log-profile suggested by Rozovskii [33] was employed to generate the vertical profile as where h C is the Chezy coefficient. The log-profile was granted with several experimental and theoretical supports as the longitudinal velocity profile [34]. For the shear flow of transverse direction, the secondary flow of the first kind needs to be reproduced in meandering channels. Thus, the linear profile reported by Odgaard [35] was adopted to generate the secondary currents as follows: The linear profile shown in Equation (13) was developed under the assumptions of steady, subcritical and isotropic turbulence, and the profile is applicable in channels with constant channel width and large c r . The depth-averaged flow fields ( s u , n u ) for generating the vertical velocity profiles were computed using the depth-averaged hydrodynamic model, which is the Hydro Dynamic Model-2D (HDM-2D) using the shallow water equations [36]. The HDM-2D has been verified for several study cases, and the results were compared with the theoretical and experimental studies [37][38][39].
In the vertical mixing step, transported particles on each horizontal plane migrate on the upper or lower plane to complete the 2D mixing. The PDM-2D calculates the vertical position of particles using the vertical mixing scheme, suggested under the assumption of negligible small vertical velocity in natural streams as where L is the number of horizontal layers;  is an arbitrary integer from 1 to L. In each computation step, a fraction of particles was evenly distributed by the full column mixing scheme, and remaining particles were randomly transported by the random mixing scheme. Through these two sequential computations, a step of shear dispersion is completed, and computations of pollutant mixing can be conducted without input of the dispersion coefficients [32]. The particle distributions can be converted to concentration fields for further analysis as where m is mass of a single particle;   , i n x t is the number of particles within water column; A is the unit area includes   , i n x t particles.

Idealized Concentration Fields to Evaluate the Frozen Cloud Assumption
For validation of the frozen cloud assumption, the dispersion coefficients for an idealized solute mixing problem were estimated using both the 2D ST-RP and the 2D S-RP. The test was conducted in the straight channel of the rectangular cross-section with 2 m width. Under the unsteady flow condition, instantaneously injected solute mixing can be calculated using the analytic solution as where   0 0 , x y is the injection point of the solute; M is the inputted mass. The unsteady flow condition shown in Equation (17) was considered to evaluate the routing procedures, while the frozen cloud assumption was not maintained. In Equation (16) [41]. The spatial distribution of concentration reproduced using Equation (16) was plotted in Figure 3, in which the concentration fields show the symmetric distributions following the Fickian dispersion. The cross-section for application of the 2D ST-RP shown in Figure  3 was determined, where the peak concentration is located where the solute clouds were drawn. By assuming that nineteen sensors are installed along uniform intervals across each cross-section, the concentration time curves were extracted to observe the temporal variations of concentration under the unsteady flow condition. The concentration-time curves at the center of each cross-section were plotted in Figure 4. According to the unsteady flow condition, the concentration curves show skewed distributions, which indicate the mixing is not stably maintained while a solute cloud passes through a crosssection. From these theoretical concentration data, performance of the 2D S-RP and the 2D ST-RP was evaluated against these solute mixing problems that violate the frozen cloud assumption.

Generation of Concentration Fields in M2 Channel Using the PDM-2D Simulation Results
Two routing procedures were applied to estimate dispersion coefficients of the twodimensional tracer test results conducted in the laboratory meandering channel (M2 channel) by Seo and Park [9]. The tracer test was conducted using 100,000 ppm ( ) of the brine and ethanol mixture as a tracer material. The outline of M2 channel was depicted in Figure 5 where the channel has rectangular cross-sections with 1 m width and two curved parts, and the experimental conditions were listed in Table 1. The tracer material was instantaneously introduced at 0.5 m upstream from the first measurement section (S1), and the concentration-time curves were measured using six conductivity sensors installed across each channel cross-section. The 2D ST-RP can be applied using the measured concentration-time curves. However, in the case of the 2D S-RP, spatial distributions of concentration are required, and the data can be obtained from aerial photos of tracer cloud or simulation results. In this study, the concentration fields were obtained from the simulation results by the PDM-2D, which has been validated by comparing to the tracer studies [32].  Following the experimental conditions shown in Table 1, the depth-averaged velocity distribution in the M2 channel was computed using the HDM-2D, as shown in Figure 5a. Based on the HDM-2D simulation results, the particle tracking simulations were conducted using the PDM-2D, and the particle distributions (Figure 5b) were converted to the concentration fields using Equation (15) (Figure 5c). For the simulation results, independence of the time step was checked for Δt = 0.125, 0.25, and 0.5 s, as shown in Figure 6. There was a difference according to the randomness of the particle motion, but the overall shape of concentration-time curves presented similar. As Δt decreased from 0.5 s to 0.25 and 0.125 s, the R-squared value of the simulation results showed 0.91, 0.88, and 0.89 in S6, and also 0.99, 0.98, and 0.96 in S9, respectively. Thus, 0.5 s of time step was used for this study, as listed in Table 1. to the measurements of S6 and S9, which mean the shape of concentration curves were reasonable to represent the tracer mixing. Thus, the simulation results by the PDM-2D adequately reproduced pollutant mixing without input of the dispersion coefficients, and the simulation results were adopted to apply the 2D S-RP.

Evaluations of the Frozen Cloud Assumption
The frozen cloud assumption, which is a fundamental theory of the routing procedures using temporal concentration data, was evaluated from the idealized solute mixing problem defined in Section 2.3.1. From the concentration data shown in Figures 3  and 4, longitudinal and transverse dispersion coefficients were estimated using both the 2D ST-RP and the 2D S-RP. Figure 8 shows the comparison results between the predicted concentration fields using the 2D S-RP and the analytic solution by Equation (16). The results shown in Figure 8 were predicted from the concentration field at t = 3 s as the upstream concentration data. The results show that the predicted concentration fields tend to underestimate the high concentration area, but the entire shape of the solute cloud was adequately estimated. The R-squared values of the predicted results denoted in Figure 8 indicate that the 2D S-RP appropriately estimated the downstream concentration fields.  Figure 9 shows contours of the concentration-time curves which were calculated using the 2D ST-RP. The predictions were performed using the concentration-time data of the section S1 shown in Figure 4 as the upstream concentration data. The contour plots of C-t-y curves show asymmetric distributions by the unsteady velocity conditions. Even if the R-squared values showed quite accurate results, the predicted concentration curves have discrepancies in the low concentration part where the concentration curves create long tails, as shown in Figure 4. From these results, it can be said that the 2D S-RP is more appropriate to estimate the concentration fields than the 2D ST-RP for a mixing problem where the frozen cloud assumption is not required.

Estimations of the Dispersion Coefficients in Meandering Channel
The applicability of the 2D ST-RP and the 2D S-RP were compared by estimating dispersion coefficients of the M2 channel described in Section 2.3.2. The 2D S-RP used the simulation results shown in Figure 5 due to absent measured spatial distributions of concentration. For application of the 2D S-RP, two concentration fields located upstream and downstream against each cross-section were selected from the simulation results. The best-fitted downstream concentration fields could be obtained when the optimal dispersion coefficients were found. Figure 10 shows the comparisons of the predicted and downstream concentration fields. The predicted concentration field adequately reproduced the spatial distribution of concentration, but inconsistency occurred in the high concentration area. The R-squared values for the comparison results were in the range of 0.98-0.99. The results by the 2D ST-RP are presented in Figure 11, in which the concentration curves were plotted in space-time domain. From the comparisons, the predicted concentration-time curves showed relatively inaccurate results, with a 0.85-0.86 range of R-squared values. The dispersion coefficients obtained using the 2D S-RP were compared with the results using the 2D ST-RP. Table 4 shows the comparisons of the dimensionless dispersion coefficients in each cross-section. Both results show that the longitudinal and the transverse dispersion increase in the cross-over (S6 and S7) and the curved (S4 and S8) parts of the channel, respectively. The overall variability of the dispersion coefficients along the entire channel were found to be greater for the 2D ST-RP than the 2D S-RP.
In natural streams, the dispersion coefficient includes shear flow and storage effects due to channel irregularities [7]. However, the results shown in Table 4 can be considered to only include the effect of shear flows because the M2 channel has uniform channel width and flow depth. Thus, the results using both routing procedures could be compared with the dispersion tensors of Equation (3), which explain the shear dispersion due to the velocity structures. For the calculations, the vertical profiles of x-and y-direction velocities measured by Seo and Park [9] were employed. To convert the dispersion coefficients shown in Table 4 to the tensor format, the coordinate transform method shown in Equations (5)-(7) was used.
The comparisons of the dispersion tensors were plotted in Figure 12, in which the dispersion tensors by Equation (3)  The local errors were deduced to occur due to the discrepancies of the simulation results, which were necessary for application to the 2D S-RP, because the predicted concentration fields shown in Figure 12

Validity of the Frozen Cloud Assumption
The 2D ST-RP has been widely used to estimate the dispersion coefficients that originated from the temporal variations of concentration curves. However, due to the frozen cloud assumption, the 2D ST-RP is unsuitable for application in a mixing problem where the dispersion changes while a solute cloud crosses a measurement point. The limitations of the 2D ST-RP may lead to the inclusion of errors in estimated dispersion coefficients, whereas the 2D S-RP does not require incorporation of the assumption because the method uses the spatial distributions of concentration. Estimated dispersion coefficients using the 2D S-RP contain a history of the shear dispersion that occurred between the upstream and downstream concentration clouds. However, the 2D S-RP is limited in certain applications due to the difficulty in obtaining the spatial distribution of concentration.
To verify the validity of the frozen cloud assumption, the 2D ST-RP and 2D S-RP were applied for the ideal mixing problem. In the case of the 2D ST-RP, which requires a frozen cloud assumption, the concentration-time curve showed a relatively larger error in the continuously changing flow characteristics due to unsteady flow. The effectiveness of the frozen cloud assumption has been evaluated in previous studies [18,42,43]. Chatwin [42] mentioned that the frozen cloud assumption is not necessary when a pollutant cloud passes through a probe fast enough. In support of the explanation, Valentine and Wood [43] showed that the variance of the temporal concentration curve can be overestimated compared to the variance calculated from spatial data when the velocity of the pollutant (c) cloud is slower than the actual flow velocity. Singh and Beck [18] also found errors in results of the routing procedure by skewed concentration curves and suggested a modified one-dimensional routing procedure. On the other hand, the 2D S-RP calculates the dispersion coefficient using the spatial concentration fields captured instantly, so if the non-Fickian mixing phenomenon due to the river topography did not appear, the dispersion coefficient could be calculated with very high accuracy.
Since the flow characteristics of most natural rivers are unsteady, the advection and dispersion of the polluted cloud due to the flow velocity change while the polluted cloud passes through a channel cross-section. Previous studies have suggested that the frozen cloud assumption is maintained, assuming that the flow remains constant while passing through the polluted cloud. However, as shown in the results of this study, it can be seen that in the region where the unsteady flow characteristics appear, the S-RP using the spatial distribution of concentration shows a more suitable result.

Applicability to Non-Fickian Mixing
The dispersion coefficients for the 2D tracer test conducted in the M2 channel were compared using the two routing procedures and the theoretical method using vertical profiles of velocity. The 2D ST-RP used the concentration-time curves measured by the previous tracer study, and the 2D S-RP employed the spatial concentration distributions computed by the PDM-2D simulation results. The dispersion coefficients by the theoretical approach were considered as exact values because the additional dispersion due to channel irregularities can be neglected in the M2 channel. Each component of dispersion tensor showed periodic patterns along the channel, and the results by the 2D S-RP adequately reproduced the spatial variations of the dispersion in the meandering channel. These spatial variabilities in dispersion coefficients due to the growth and decay of secondary flows in the meander cycle have also been reported in a study by Boxall and Guymer [8,44]. However, the spatial changes by the 2D ST-RP showed discrepancies, which were deduced by the skewed concentration-time curves violating the frozen cloud assumption.
The non-Fickian mixing due to the storage zone of the river forms a long tail in the concentration-time curve, which causes difficulties with applying the mixing theory following the Fick's law. The storage zone causes an imbalance between vertical turbulent diffusion and shear advection and leads to asymmetric concentration curves [45]. This process was described in detail in the study by Denton [46]. When a pollutant cloud passes through a probe installed at a middle of water depth, concentration of the pollutant cloud located near the water surface is measured first due to fast flow velocity in the upper layer. After that, the probe measures concentration of the pollutant cloud near the channel bed, which goes to the upper layer by vertical mixing and relatively slow flow velocity. Therefore, it shows a difference from the shape of the polluted cloud observed at a specific time, and the frozen cloud assumption becomes invalid. The tracer test results of the M2 channel corresponds to the initial period that causes the imbalance of vertical diffusion and shear advection and shows a spatiotemporal asymmetric concentration distribution even though the storage zone is not included [30]. Since the 2D ST-RP and the 2D S-RP were derived from the Fickian dispersion model, both methods showed differences from the velocity-based dispersion tensor calculation results, as shown in Figure 12. However, even though the limitations originated from the Fick's law, the 2D S-RP is more advantageous for estimating the dispersion coefficient when the solute mixing shows non-Fickian characteristics.

Conclusions
In this study, two routing procedures using temporal concentration curves and spatial concentration distribution were evaluated against two solute mixing problems. Performance of the two routing methods was analyzed in two aspects: (1) validity of the frozen cloud assumption and (2) applicability to the non-Fickian solute mixing problem. The analysis results are summarized below.
• For a solute mixing problem under the unsteady flow condition, the 2D S-RP provides quite accurate estimation results of dispersion coefficients when the mixing shows the Fickian dispersion.

•
The temporal concentration curves present the non-Fickian mixing due to the unsteady flow condition even though the solute cloud shows the Fickian dispersion. Thus, the dispersion coefficients by the 2D ST-RP using the temporal data contain errors due to violation of the frozen cloud assumption.

•
Both the dispersion coefficients calculated by the 2D ST-RP and the 2D S-RP showed errors against to the results using the velocity-based method for the solute mixing in the initial period where the non-Fickian mixing occurs. These discrepancies showed that the two routing methods were derived based on the Fickian dispersion model.

•
The results obtained using the 2D S-RP in the meandering channel more accurately exhibited the spatial variability along the meander cycle of dispersion coefficients compared to the 2D ST-RP.
From the aforementioned analysis results, the 2D S-RP is able to provide more reasonable results of dispersion coefficients in both the Fickian and the non-Fickian mixing problems than the 2D ST-RP if spatial distributions of concentration are available. These conclusions are limited to the laboratory-scale application results since verifications have not been conducted on the natural river scale which has more complicated mixing properties. In future studies after the concentration measurement technique through the image analysis technique currently being studied has been established, it is expected that the applicability of the 2D S-RP can be evaluated on the field scale.

Data Availability Statement:
The data presented in this study will be shared by the authors if requested.