Revised predictive equations for salt intrusion modelling in estuaries

Introduction Conclusions References Tables Figures


Introduction
Predictive methods to determine salinity profiles in estuaries can be very useful to water resources managers, particularly when applied to ungauged estuaries where only a minimal amount of data is available.Before any decision is made on collecting detailed field observations, it is useful to obtain a first estimate of the strength and range of the salt intrusion in the area of interest.Such estimate can be made if there are predictive equations available to compute the longitudinal salinity profile along the Introduction

Conclusions References
Tables Figures

Back Close
Full estuary.With reliable predictive equations, water managers are able to estimate how far salt water intrudes into the river system under different circumstances, and more importantly, how interventions may change this situation.
The one-dimensional salt intrusion model of Savenije (1993b) makes use of the Van der Burgh and dispersion equation to represent the longitudinal variation of the salinity.
The Van der Burgh and dispersion coefficient at the ocean boundary are obtained by calibration of the simulated salinity curve to observations.Savenije (1993b) established a predictive equation for each of these parameters, so that the longitudinal salinity distribution could be estimated when data were lacking or to monitor the impact of interventions, such as dredging or fresh water withdrawal.The predictive equations have subsequently been modified and tested by several researchers including Savenije (2005), Nguyen and Savenije (2006), Kuijper and van Rijn (2011), and Shaha and Cho (2009).
In this paper, we shall revisit the predictive equations in the light of new insights on how friction and estuary shape affect tidal mixing, makes use of new observation in 7 previously ungauged estuaries in Malaysia that were sampled through a consistent approach.As a result, we present the fully revised and more accurate predictive equations for the Van der Burgh coefficient and for the boundary value of the dispersion at a well identifiable location, based on tidal average (TA) condition.

One-dimensional analytical salt intrusion model
The analytical one-dimensional salinity model developed by Savenije (1993bSavenije ( , 2005Savenije ( , 2012)), presented below, is used to simulate the salinity profile in the estuaries studied.In a steady state situation, the partial temporal derivative in the salt balance equation is zero.Considering a constant fresh water discharge Q f [L where S = S(x) [M L −3 ] and D = D(x) [L 2 T −1 ] are the salinity and dispersion at tidal average (TA) condition.Since discharge has a negative value, the absolute value of Q f is taken in Eq. ( 1). S f [M L −3 ] represents the fresh water salinity.Making used of the Van der Burgh equation in combination with the salt balance equation, Savenije (2005Savenije ( , 2012) ) described the relation between dispersion and salinity to be: in which K [-] is defined as the Van der Burgh coefficient (shape factor).Substituting Eq. (1) into Eq.( 2), the differential equation for the tidal average longitudinal salinity distribution is expressed as: Integration of Eq. (3) leads to: The symbols S [M L −3 ] and D [L 2 T −1 ] are the steady state salinity and dispersion coefficient at location x, while S 0 [M L −3 ] and D 0 [L 2 T −1 ] are the salinity and dispersion at the estuary mouth.In alluvial estuaries, the variation of the estuaries shape over the distance upstream can be expressed as an exponential function (Savenije, 2005(Savenije, , 2012;;Nguyen and Savenije, 2006;Zhang et al., 2011) as:  5) into the Eq. ( 2) and integration gives: Here, β [-] is the dispersion reduction rate.At the salt intrusion limit (upstream) where only fresh water discharge exist, the dispersion coefficient becomes zero and x is equal to the salt intrusion length L [L].Hence, the intrusion length is expressed by: Equations ( 4) to ( 9) are the general equations used to compute the longitudinal salinity distribution.
3 Existing predictive equations

Van der Burgh's coefficient
Van der Burgh's coefficient K is also known as the "shape factor" of the salinity curve (Savenije, 1993a).Based on salinity measurements of 15 estuaries, Savenije found that K is strongly related to the geometry (the convergence length a or b and the width B [L]) and its influence is more significant at the tail of the salinity curve (upstream).Moreover, Savenije (1986Savenije ( , 1989) ) observed that every estuary had its own characteristic value of K , ranging from zero to one.Assuming that the Van der Burgh coefficient is Introduction

Conclusions References
Tables Figures

Back Close
Full −2 ] is the gravity acceleration.More than 10 years later, Savenije (2005) and Nguyen and Savenije (2006) made used of an expended database, modified the predictive equation involving more parameters: The symbols and A [L 2 ] refer to the tidal excursion, tidal range, and a boundary value for the cross-sectional area, respectively.This relation had a correlation of 0.96 and seemed very promising.However, as can be seen from the equation, the Chezy roughness C [L 0.5 T −1 ] and damping δ H [L −1 ] had to be computed from tidal dynamics analysis.

Dispersion coefficient
Dispersion is not a physical parameter; it is rather the product of averaging, representing the mixing of saline and fresh water in an estuary as a result of residual circulation induced by density gradients (gravitational circulation) and tidal movement.In salt intrusion modelling, the definition of dispersion is often unclear as it is scale dependent and not directly measurable.The role of dispersion is only meaningful if it is related to the appropriate temporal and spatial scale of mixing, which here we identify as the tidal period (time scale), tidal excursion (longitudinal mixing length), estuary width (lateral Introduction

Conclusions References
Tables Figures

Back Close
Full mixing length) and depth (vertical mixing length).A physically based description of the dispersion would allow the analytical solution of the salt intrusion profile.Dispersion due to gravitational circulation has been studied since 1957, as summarized by Fischer (1976).This type of dispersion is also known as density driven dispersion between the two main sources: sea water and fresh river water.Schultz and Simmons (1957) were some of the first to relate buoyancy to mixing in estuaries, whereby they introduced the ratio between fresh water discharge and tidal volume to represent the degree of stratification.This ratio is also known as the Canter-Cremers number N [-] as defined by Harleman and Abraham (1966).The buoyancy effect or stratification in an estuary can also be represented by the Estuarine Richardson Number N r [-] which is the ratio of potential energy of the buoyant fresh water to the kinetic energy of the tide: where ρ [M L −3 ] is the water density, ∆ρ [M L −3 ] is the density difference over the salt intrusion length, and υ [L T −1 ] is tidal velocity amplitude.The difference between N and N r lies in the densimetric Froude number F d [-] which is expressed as: Since then, researchers have tried to look for a relation between dispersion and estuarine numbers.Laboratory results of WES flume, Delft flume and Daniels (1974) indicated an agreement with the result of Fischer (1972) in computing the salt intrusion length, using shear velocity instead of mean velocity in the Estuarine Richardson number.Subsequently, the relationship between the dispersion and modified N r also gave good correlation for all the other cases (mostly flume experiments).Thatcher and Harleman (1972) suggested that the longitudinal dispersion is proportional to the salinity gradient and included this in his one dimensional analytical salt intrusion model, Introduction

Conclusions References
Tables Figures

Back Close
Full which later was used by Fischer (1972) to model the vertical salinity and velocity distribution.A disadvantage of all these methods was that they did not account for convergence (implicitly assuming an infinitely large convergence length) and that the tidal excursion, as the most important mixing length scale, was missing in the derivations.Deriving the dimensionless dispersion coefficient from scaling the steady state salt balance equation, Savenije (2005) developed the following empirical predictive relation for the longitudinal dispersion at the estuary mouth for HWS: The estuary shape was represented by the ratio of the averaged depth h [L] to the convergence length a, while the dispersion was made dimensionless by the tidal velocity amplitude and tidal excursion which was not considered in any of the earlier studies.The applicability of these predictive equations has been widely tested in many estuaries including multi-channel estuaries.Kuijper and van Rijn (2011) later modified the predictive equation including the dimensionless friction (C 2 /g).The equation was divided into two depending on the types of channel -prismatic and convergent: Convergent channel: Prismatic channel: where u [L T −1 ] is the fresh water velocity.These equations can be used to calculate dispersion locally at any location.However, the weakness is that it is no longer predictive as calibration is required to determine the p values to fit the measurement data.l 1 and l 2 are the coefficients for the convergent and prismatic channel, respectively.

Salt intrusion length
Several researchers have tried to develop a general relation for the salt intrusion length.
The development of such predictive equations was done empirically based on a reasonable amount of data.A pioneer effort was made by van der Burgh (1972), making use of prototype information from the Dutch and German estuaries.His equation for the salt intrusion length as summarized by Savenije (1992Savenije ( , 1993bSavenije ( , 2005) ) is as follows: with: and: In this equation, L TA [L] is the salt intrusion length at TA situation, F [-] is the Froude number, and P t [L 3 ] is the tidal flood volume.
A year later, Rigter (1973) established another empirical equation for the salt intrusion length at LWS using experimental data of Delft Hydraulic Laboratory and the WES flume.number F d [-].Fischer (1974) re-conducted the empirical analysis using the same batch of data from Rigter (1973) and came out with a slightly different result.
Here, it can be seen that despite using the same parameters, the powers of f D , F d and N are different.
About 20 years later, Van Os and Abraham (1990) established a similar equation as Rigter (1973) with a slightly different coefficient: All these methods were based on flume data with prismatic geometry.Savenije (1993bSavenije ( , 2005Savenije ( , 2012) ) who explicitly accounted for channel convergence and the tidal excursion, developed a predictive equation for the salt intrusion length at HWS.The reasoning was that the maximum salt intrusion length occurs during HWS, which is most important for water resources management.Based on Eq. ( 14), the equation reads: where υ 0 [L T −1 ] is the tidal velocity amplitude at the mouth.It is worth noting that Savenije follows Van der Burgh's equation, with an additional shape indicator referring to the area convergence length a.
Most of the empirical equations discussed above are based on LWS except for Van der Burgh's and Savenije's methods which are based on TA and HWS, respectively.However, they can easily be brought in agreement with each other by adding E/2 or E to L HWS , respectively.Here, we aim to develop a universal predictive equation for Introduction

Conclusions References
Tables Figures

Back Close
Full estimating the Van der Burgh and dispersion coefficient for TA condition, which can be applied in the salt intrusion model to predict the salinity profile for any estuary worldwide under different tidal and flood conditions.

Methods
In this paper, the main focus is on the mixing mechanisms which lead to longitudinal dispersion in estuaries: the tide and density driven dispersion.Key parameters are developed based on measurable parameters of geometry, tidal hydraulics and fresh water discharge.In total 89 measurements data of 30 estuaries worldwide have been used to develop the predictive equations.Measurements in 7 newly surveyed estuaries were collected from 2011 to 2013 in Malaysia (Gisen et al., 2014), whereas the remaining were compiled by revisiting existing data available in the database of Savenije (2005) and from professional reports.The locations of the estuaries studied are displayed in Fig. 1.
Adjustments have been made to the geometry (see Fig. S1 in the Supplement) and salinity analysis for some of the estuaries to ensure consistency in the input data used.
The entire dataset was split into two: reliable and less reliable data.The reliable dataset have been used to develop the predictive equations, whereas the less reliable ones have been used for verification purposes.The study was performed based on Savenije (1993bSavenije ( , 2005Savenije ( , 2012))'s method for predicting K and D 0 with some modifications.The modifications include: -All geometry and tide information used refers to the well identifiable inflection point x 1 as the boundary condition.
-Analyses were performed on TA condition instead of HWS, which is consistent with the geometry information.
-Estuary roughness and the ratio of estuary width to river width have been added in the predictive equations.Introduction

Conclusions References
Tables Figures

Back Close
Full -The parameters chosen are mostly independent and easy to observe without the need for prior calibration.
Although the predictive equations are based on the tidal average (TA) situation, one can still compute the salinity distribution for both low water slack (LWS) and high water slack (HWS) by shifting the curve over E/2 in seaward and landward direction.

Selecting the dimensionless ratios
Revising the parameters selected by Savenije (1993bSavenije ( , 2005)), we found that the latter contained some parameters that required tidal dynamics analysis while one of the ratios was not dimensionless.The followings are the dimensionless ratios selected for the revised predictive equation for the Van der Burgh coefficient: where B f [L] is the river regime width, and λ 1 = gh 1 /r s [L] is the wave length at the inflection point with r s [-] being the storage width ratio (defined as the ratio between storage width and stream width).The symbols represent the estuary width, tidal excursion, averaged estuary depth, tidal range and width convergence length at the inflection point x 1 .It is worth noting that the roughness was obtained through calibration using the tidal dynamics solution of Cai et al. ( 2012) which makes use of observed tidal damping.In the above equation, it can be seen that all parameters used have been defined at the inflection point x 1 .It is also important to note that the convergence length adopted is of the second reach not the first part of the estuary.Generally the tidal indicators E and H, are defined at the mouth.In order to obtain the tidal excursion and tidal damping at the inflection point, a projection can be made considering tidal damping as follows (Kuijper and van Rijn, 2011): Introduction

Conclusions References
Tables Figures

Back Close
Full where the damping factor δ H also follows from the tidal dynamics simulation of Cai et al. (2012).The values of H 1 and E 1 used in the dimensionless ratios represent the condition of spring tide, where υ is considered to be closed to 1 m s −1 (Bruun and Gerritsen, 1960;Pethick, 1984;Langbein, 1963).This is to ensure that K is time-independent representing a general characteristic of an estuary.As a result, E essentially reflects the tidal period as described in Eq. ( 27) (see also Table 1).
For the dispersion coefficient, eight dimensionless ratios have been selected with 18 different types of equations including the one of Savenije (1993bSavenije ( , 2005) ) as benchmark.
The dispersion coefficient is represented in dimensionless form as: where the gauging station was also considered in the analysis.Stepwise multiple regression analysis has been used to identify the best combination of the dimensionless ratios in predicting K and D 1 .The efficiency of the established equations was examined by comparing the correlation coefficient R 2 and the standard error SE.The predicted results calculated by the most suitable equations were plotted against the calibrated values to evaluate their predictive performance (see Fig. S2).

Substitution of predictive equations in the salt intrusion model
Since the predictive dispersion is computed at the inflection point x 1 , reverse calculation has to be done to obtain the dispersion at the mouth.This is necessary to enable the simulation of the longitudinal salinity distribution starting from the mouth to the salt intrusion limit.Inverse integration from x 1 to x = 0 of Eq. ( 2) yields for the dispersion at the estuary mouth: with: and: where β rev [-] is the reversed dispersion reduction rate, whereas and α 1 [L −1 ] are the cross-sectional area, dispersion coefficient and mixing number at the inflection point, respectively.It is important to note that the convergence length a 1 [L] applied in Eqs. ( 30) and ( 31) is of the first section of the estuary.The relation between dispersion and salinity is then expressed by: Figures

Back Close
Full where S 1 [M L −3 ] refers to the salinity at the inflection point.Substituting the tidally average dispersion coefficient into the general form of the salt intrusion length of Savenije (1993bSavenije ( , 2005) ) yields: with: Note that all parameters used in these equations refer to the inflection point.We obtain the salinity profile at HWS and LWS by moving the salinity curve over E/2 in the upstream and downstream direction.Similarly, the maximum salt intrusion length can be obtained by shifting the intrusion length at TA in landward direction by half of the tidal excursion at the mouth as: and the LWS intrusion length by moving the tidal excursion seaward.reliable dataset (see Table S2 and S3 in the Supplement).This distinction was made based on the following criteria.
Criteria for classifying estuaries as reliable: the estuary is generally in steady state condition; the fresh water discharge is estimated, observed or measured correctly; the estuary is alluvial and undisturbed; complete measurement data for tidal dynamics and salinity analysis are available.
Criteria for classifying estuaries as less reliable: -The estuary is not in steady state particularly during low river discharge.This depends on the ratio of the time scale of system response to the time scale of discharge reduction (see Savenije, 2012) (NSS).
-The estimation of the fresh water discharge is uncertain (UQ).
-The estuary may not be alluvial (e.g.dredged, modified or constricted by rocky banks) (NA).
-Information on tidal dynamics and salinity is lacking or unclear (IL).
The estuaries that fall under category NSS, UQ, NA, and IL are listed in Table 1.It is worth noting that only the reliable set is used in regression analysis.The less reliable ones are merely plotted for verification purpose.Introduction

Conclusions References
Tables Figures

Back Close
Full Results from the stepwise multiple regression analyses show that the best combinations of the dimensionless ratios to represent the Van der Burgh predictive equation are: or: where Eq. ( 39) is the simplified form.The correlation coefficient R 2 and the standard error SE obtained for predictive K equation is 0.72 and 0.11, respectively.If we had used the cross-sectional area convergence a 2 instead of b 2 , then the correlation would be slightly poorer.Hence the width convergence is a better indicator, which is fortunate because it is easier to determine.From the equation, we can see that the parameters that have the most influence on the Van der Burgh coefficient is the tidal period, follows by the gravity acceleration, the storage width ratio and tidal velocity amplitude which have the power of 0.97, 0.93, 0.84 and 0.71, respectively.The importance of the friction appears to be minor, which is also fortunate as C is not directly observable.The estuary to river width ratio shows higher power than the convergence length, which indicates that the width is a better shape indicator.Finally, we should realize that 0 ≤ K ≤ 1 according to Savenije (1993a).For prismatic channel where b 2 becomes infinity, K approaches zero, implying constant dispersion.Introduction

Conclusions References
Tables Figures

Back Close
Full Figure 2 shows the plot of the predicted K against the calibrated value.All the reliable data points appear to fall close to the perfect agreement line.About half the unreliable data points were outliers particularly the Gambia (30) and Tejo (27) Estuary which lie much further away from the perfect agreement line.This is not strange in the sense that the Tejo Estuary is not entirely alluvial, and its narrow and deep mouth caused by a rock outcrop formation turns it into a fjord type estuary.As for the Gambia, it is an unsteady state estuary.Nevertheless, for the rest of the outliers we believe that they would fit better if good data had been available.The results are summarized in Table 1.

Predictive equation for the dispersion coefficient D
In this study, 18 combinations of the dimensionless ratios were established by multiple regression method of which the results are displayed in Table S1 (equations) and Fig. S2 (correlations and standard error).By observing the exponent, it can be seen that the power of the Estuarine Richardson number N r varies little, indicating the clear correlation with N r compared to the other parameters.The next parameter that has a high exponent is the dimensionless roughness, of which the inclusion improves the correlation.As for the rest of the dimensionless ratios, it appears that the contribution is minimal.Hence, the best equations chosen for further analysis are: More information about the equations tested is provided in the Supplement (Table S1).The Eqs. ( 40), ( 41) and ( 42) correspond with Eqs.(R2), (R4) and (R9) in Table S1.Introduction

Conclusions References
Tables Figures

Back Close
Full It is interesting to note that the performance of the benchmark equation of Savenije (1993bSavenije ( , 2005) ) (Eq.R1) is rather poor, with R 2 and SE of 0.67 and 0.33.These significant differences may be caused by the homogenisation of the input information (e.g.geometry), and the use of selective data for calibration.With more or less equal performance, it is decided that the simplest equation with the best performance is the most attractive one.Therefore, we conclude that Eq. ( 41) is the best to predict the tidal average dispersion coefficient at x 1 .This is also theoretically the most attractive, since laboratory experiments have demonstrated that both N r and the roughness are key parameters.Nevertheless, if the Chezy roughness is unknown, then Eq. ( 40) can be applied.Equations ( 40), ( 41) and ( 42) have an R 2 of 0.84, 0.86 and 0.80 with SE of 0.14, 0.13 and 0.15, respectively.We can also conclude that although estuary shape is the key in defining K , the dispersion boundary condition D 1 appears to be determined by hydraulic parameters.Figure 3 displays the plots of the predicted D 1 and α 1 against the calibrated values for both the reliable and less reliable datasets using Eqs.( 40), ( 41) and ( 42).Here, it is shown that all the reliable data points fall nicely within the range of a factor 1.5.Some of the less reliable data points are also within or near the range except several obvious outliers such as the Delaware, Schelde, Pungue, and Tejo.This is because the Pungue is often in an unsteady state condition, while the Schelde is dredged, and the Tejo and Delaware are not completely alluvial.In addition, the doubt on the accuracy of the discharge data is also one of the factors contributing to poor results.It can be seen that all the predictive equations selected have underestimated the values of the dispersion coefficient for the outlying data points, indicating a possible underestimation of the river discharge.
Comparing the outliers in both plots, it appears that the unreliable data are distributed closer to the reference lines if the dispersion is represented in term of the mixing number.This implies that the fresh water discharge is partly to blame for the discrepancy.
The data used for the regression and results of the predicted dispersion are tabulated in Table S2 in the Supplement.

Modified predictive equation for maximum salt intrusion length L HWS
Comparison between the predicted and calibrated salt intrusion length has been done for HWS condition instead of TA.This is because the salt intrudes furthest into the river system at HWS, and the maximum intrusion is the information water managers are most interested in.Substituting the predictive dispersion Eqs. ( 40), ( 41) and ( 42) into the general form for salt intrusion length yields: 43) 44) Figure 4 shows the performance of these equations in predicting the maximum salt intrusion length.In the plots using Eqs.( 43), ( 44), and ( 45), all data points fall within the range of factor 1.5 except the Solo Estuary.The list of data and the results are summarized in Table S3 in the Supplement.It appears that the predictive equations overestimated the intrusion length in the Solo Estuary.This may be due to the nearly prismatic shape of the channel which has a very long convergence length of 226 km.

Longitudinal salinity profiles
The salinity curve can be computed by applying Eqs. ( 43), (44), and ( 45) with the different dispersions calculated by each of the predictive measure developed.Considering the substantial amount of salinity measurements available, only the salinity profiles of the 7 newly surveyed estuaries are discussed.The plots of all salinity profiles will be uploaded as electronic material at the website -http://salinityandtides.com.demonstrates the performance of the simulated longitudinal salinity distribution with and without calibration of K and D 1 .
From the salinity curve comparison, it appears that all the predictive equations do not performed very well for Kurau and Bernam estuaries.This may be caused by the uncertainty in discharge data.The Kurau and Bernam discharge calculation were based on the discharge observed in a small part of the catchments of about 12 and 20 % of the total area, respectively (Gisen and Savenije, 2014).Thus, it is possible that we may have underestimate the discharge draining into the Kurau Estuary, and overestimate the one for Bernam Estuary.It is also interesting to note that Eq. ( 45) works better in predicting the salinity distribution for some of the estuaries such as the Perak, Linggi and Endau estuary.As for most of the cases, Eq. ( 44) appears to give the best fit.The difference in the performance of these equations suggests that there is a possibility that the equations are subject to improvement if more reliable measurements are available.Thus, it is appropriate to retain the three Eqs.( 43), ( 44), and (45) for consideration.

Discussion
Before Savenije's (1993a) effort to develop predictive equations for the Van der Burgh and dispersion coefficient, these parameters could only be obtained by calibration.Without site measurements, it was impossible to make any estimate of the salinity distribution along an estuary.The predictive equations of Savenije (1993aSavenije ( , 2005) ) were able to estimate the value of K and D reasonably well in reference to the calibration data.However, after re-evaluating and re-analysing the available data, we found that the equations do not work as well for all estuaries.
In this study, we have collected an additional 32 salinity profiles from 16 new estuaries for consideration in the analysis.Moreover, the measurements were split into two datasets to make sure that only the reliable data were used for establishing the revised equations.In previous work, the data were not split.The selection process is important so that the results are not influenced by incomplete or uncertain data.Re-examining Introduction

Conclusions References
Tables Figures

Back Close
Full the available measurements from the old database ensures that all data used are accessible and consistent.The new compilation also provides a section containing important information about each measurement (see electronic additional material -salinity worksheet at http://salintyandtides.com).
Another important modification in this work is the change in the selected boundary condition.In this research, we decided to process the cross-sectional data in reference to the tidal averaged situation (TA).whereas previous methods were based on HWS and LWS, which led to inconsistencies because the geometry during low and high water can be different from tidal average situation.Moreover, in this study we fixed the location of the downstream boundary at the inflection point x 1 and not at the estuary mouth (adopted by all the earlier researchers).The reasons and advantages of moving the downstream boundary to the inflection point are: to eliminate the difficulty of determining the exact location of the estuary mouth; to reduce the effect from wind and waves; to eliminate the dilemma of which geometry parameters to use in the predictive equation.
In Savenije's (1993aSavenije's ( , 2005) ) and Kuijper and van Rijn's (2011) predictive model, the cross-sectional area convergence length applied to calculate the salt intrusion length was the weighted value obtained from an iteration process.Hence, with the change of the downstream boundary to x 1 , this process is no longer needed and the predictive measures are more consistent.
The new set of dimensionless ratios proposed in this study to establish the predictive equation for K contains mostly measurable independent parameters.The selection was made based on the existing equations, considering only the parameters that are easy to obtain.It is worth noting that the ratio (1 − δ H b) has been removed from the equation because the damping changes from spring to neap tide.Furthermore, it also decreases or decreases towards upstream and is highly influenced by fresh water discharge.The river to estuary width ratio has been added in the new equation as an Introduction

Conclusions References
Tables Figures

Back Close
Full additional geometry indicator besides the depth and convergence length.This ratio appeared to have great influence on K .
For the predictive dispersion equation, the ratio of the depth to the convergence length is no longer important, but the longitudinal length scale E and velocity amplitude υ remain important in the scaling of D 1 .The elimination of h/a allows the new equation to be applied also in prismatic channels.In the old equation, when a 2 approaching infinity, the calculation became invalid.Since Kuijper and van Rijn (2011) suggested that the friction parameter is related to the vertical mixing, g/C 2 has been included in this new equation and it indeed improved the correlation.Savenije (2005) did not consider roughness in his predictive equation for dispersion.
Although some improvements and simplicity have been introduced in this study, there are limitations in using the new equations.Until now, we have only taken into account single network estuaries.Furthermore, it has implicitly been assumed that no water is entering or leaving the tributaries in the estuary region.If these are large tributaries or large areas draining on the estuary, then these should be accounted for.From the plot of Van der Burg's coefficient, we found that the performance in predicting K is rather low.This indicates that the equation has to be used with caution.Another constraint in using the developed equations is the friction factor.The Chezy roughness is not directly measurable and can only be obtained by calibration using a tidal dynamics model.However, if this information is impossible to get, it can be neglected (the correlation only decrease to 0.70 for the predictive equation without roughness).If cross-sectional area information is lacking, then b 1 can be used to replace a 1 .For the depth estimates, one can made used of the method presented by Gisen and Savenije (2014) which links h 1 to the bankfull discharge.

Conclusions
Calibrating K and D 1 is only possible if measurement of the salinity distribution is available.In a situation where data are limited, a predictive equation is required to estimate Introduction

Conclusions References
Tables Figures

Back Close
Full The analysis based on tidal average conditions enables the entire process to be carried out consistently, whereby model and data errors can be reduced.The obtained salt intrusion can easily be converted from TA to HWS by adding half of the tidal excursion.The performance of the predictive equation for K is rather good with a R 2 value of 0.72.For the dispersion, the correlation of 0.86 seems very promising.All the reliable data points fall within a factor of 1.5 for both the predicted K and D 1 results.Some less reliable ones are also within this range.This indicates that the predictive equations developed are appropriate to be applied in getting a first estimate of the K and D 1 .Subsequently, the longitudinal salinity distribution in an estuary can be estimated.
Hence, these tools can be very helpful for water managers and engineering to make preliminary estimates on the salt intrusion in an estuary of interest and to analyse the impact of interventions.Finally, it is recommended to collect more reliable measurements to strengthen the development of the empirical relationships.New data are also required for validation purposes.1.Data used to develop the predictive equation for the Van der Burgh coefficient K .

No
Estuary Full Discussion Paper | Discussion Paper | Discussion Paper | and tidally averaged cross-sectional area A [L 3 ], the salt balance equation for tidal average condition (TA)can then be written as: Discussion Paper | Discussion Paper | Discussion Paper | where a [L] and b [L] representing the cross-sectional area and width convergence length, A 0 [L 2 ] and B 0 [L] are the cross-sectional area and width at the mouth, B [L] is the width of estuary at distance x [L] (towards upstream).Substituting the exponential relation of Eq. ( Discussion Paper | Discussion Paper | Discussion Paper | not time-dependent, Savenije (1993b) established an empirical predictive equation for K as: K = 0.16 × 10 −6 h 0h 0 [L], H 0 [L], and B 0 [L] are the depth, tidal range and width at the estuary mouth, respectively.The symbol T [T] represents the tidal period, while b [L] is the width convergence length, and g [L T Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | 20) It is important to note that Van der Burgh's coefficient K is replaced by the Darcy-Weisbach's roughness f D = 8g/C 2 [-] and F is represented by the densimetric Froude Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | N r1 [-] is the Estuarine Richardson number with υ 1 [L T −1 ] being the tidal velocity amplitude, both at the inflection point.It is important to note that the values taken for E 1 and H 1 in the dispersion analysis are based on the real-time data captured during measurements and the depth is referring to the depth at the inflection point.In general, the density different between the saline and fresh water is taken as (25/35)S 0 kg m −3 and the fresh water density as 1000 kg m −3 .The fresh water discharge data were adjusted for the 7 newly surveyed estuaries so that the runoff contribution downstream of Introduction Discussion Paper | Discussion Paper | Discussion Paper | Data Data were divided into two categories: reliable and less reliable.There are 47 measurements grouped under the reliable dataset, and 38 measurements under the less 753 Discussion Paper | Discussion Paper | Discussion Paper | equation for the Van der Burgh coefficient K Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | the desired variables.A good predictive equation should be simple (parameters can be easily measured) and efficient.The predictive equations established in this study consist of mostly measurable independent parameters.Options are suggested for the case in which data are very limited.The adjustment of the downstream boundary to the inflection point has clarified the selection of the right geometry parameters to be used and the position of the downstream boundary.
Discussion Paper | Discussion Paper | Discussion Paper | Table Fig. 5.1 Global map showing the locations of the estuaries studied.

Figure 1 .Figure 2 .
Figure 1.Global map showing the locations of the estuaries studied.

Figure 3 .
Figure 3. Performance of the predictive equations for the dispersion coefficient (left panel) and mixing number (right panel) against calibrated values.