Investigation and comparison of one-dimensional (1-D) analytical models prediction for salt intrusion condition in two selected estuaries

Abstract

In this paper, two analytical models used worldwide to assess salinity variation in alluvial estuaries are applied to the Ashtamudi estuary, a Ramsar site, southwest coast of India; and Bouregreg estuary, in northwest part of Morocco. The estuaries’ bathymetry is described by an exponential function. Both models are quite similar and use a predictive equation for the dispersion in the estuary mouth (D₀). The major difference between the two models is the use of the constant value of K?=?0.5 for the Van der Burgh coefficient (K) and the introduction of the correction factor ζ, which is a function of damping (δ) and shape (γ). The performance of these two models was evaluated by comparing their results with field measurements; this revealed that both analytical models apply well to both the estimation of salinity distribution and the prediction of salt intrusion in the Ashtamudi and Bouregreg estuaries (Ashtamudi: RMSE = 0.60–1.22 ppt; Bouregreg: RMSE = 0.92–2.71 ppt). One model agrees more with the field measurements of salinity distribution along the estuaries axis; the second underestimate and overestimate some values of salinity distribution along the estuaries. Possibly, the constant value of K?=?0.5 for the Van der Burgh coefficient (K) has applicability limits for the estuaries under tidal conditions. The specifying of the parameterization may be a field of research.     

Prediction of salinity intrusion in the sheltered estuary of Terengganu River in Malaysia using 1-D empirical intrusion model

Generally one dimensional(1-D) empirical salinity intrusion model is limited to natural alluvial estuary. However,this study attempts to investigate its ability to model a sheltered alluvial estuary of the Terengganu River in Malaysia. The constructed breakwater at the mouth of the river shelters the estuary from direct influence of the open sea. The salinity density along the estuary was collected during the wet and dry seasons for scenarios before and after the constructed breakwater. Moreover, the freshwater discharges, tidal elevations and bathymetry data were also measured as model inputs. A good fit was demonstrated between simulated and observed variables,namely salinity distribution and intrusion length for both scenarios. Thus, the results show that 1-D empirical salinity model can be utilized for sheltered estuarine condition at the Terengganu Estuary, but with an appropriate determination of an initial point. Furthermore, it was observed that the salinity intrusion in the study area is largely dependent on the freshwater discharge rather than tidal elevation fluctuations. The scale of the salinity intrusion length in the study area is proportional to the river discharge of the –1/2 power. It was appeared that the two lines of the 1-D empirical salinity model and discharge power based equation fitted well to each other, with the average predicted minimum freshwater discharge of 150 m3/s is going to be required to maintain acceptable salinity levels during high water slack(HWS) near the water intake station, which is located at 10.63 km from river mouth.     

Using salt intrusion measurements to determine the freshwater discharge distribution over the branches of a multi-channel estuary: The Mekong Delta case

The fresh water discharge is an important parameter for modelling salt intrusion in an estuary. In alluvial converging estuaries during periods of low flow, when salinity is highest, the river discharge is generally small compared to the tidal flow. This makes the determination of the fresh water discharge a challenging task. Even if discharge observations are available during a full tidal cycle, the fresh water discharge is seldom much larger than the measurement error in the tidal discharge. Observations further upstream, outside the tidal region, do not always reflect the actual flow in the saline area due to withdrawals or additional drainage. Discharge computation is even more difficult in a complex system such as the Mekong Delta, which is a multi-channel estuary consisting of many branches, over which the freshwater discharge distribution cannot be measured directly. This paper presents a new approach to determine the freshwater discharge distribution over the branches of the Mekong Delta by means of an analytical salt intrusion model, based on measurements made during the dry season of 2005 and 2006. It appears that the analytical model agrees well with observations and with a hydraulic model. This paper demonstrates that with relatively simple and appropriate salinity measurements and making use of the analytical salt intrusion model, it is possible to obtain an accurate discharge distribution over the branches of a complex estuary system. This makes the analytical model a powerful tool to analyze the water resources in tidal regions.     

Modelling salinity and circulation for the Columbia River Estuary

Simulations of the time and depth-dependent salinity and current fields of the Columbia River Estuary have been performed using a multi-channel, laterally averaged estuary model. The study simulated two periods. The first, in October 1980, with low riverflow of about 4,000m³s⁻¹, which showed marked changes in the salinity intrusion processes between neap and spring tides; and second, in spring 1981, with high riverflow varying between 7,000 and 15,000m³s⁻¹, which showed the rapid response of the salinity intrusion to changes in riverflow and that vertical mixing did not change character with increasing tidal energy because of the maintenance of stratification by freshwater flow. An extreme low flow simulation (riverflow of 2,000m³s⁻¹) showed a more partially mixed character of the estuary channels with tidal dispersion of salt across the Taylor Sands from the North Channel to the upper reaches of the Navigation Channel. Asymmetries in the non-linear tidal mean flows, in the flood and ebb circulations, and salinity intrusion characteristics between the two major channels were observed at all riverflows. The model confirms Jay and Smith's (1990) analysis of the circulation processes in that tidal advection of salt by the vertically sheared tidal currents is the dominant mechanism by which the salinity intrusion is maintained against large freshwater flows. An accurate finite-difference method, which minimized numerical dispersion, was used for the advection terms and was an important component in reasonably simulating the October neap-spring differences in the salinity intrusion. The simulations compare favorably with elevation, current and salinity time series observations taken during October 1980 and spring 1981.     

Modelling Salt Intrusion and Nitrate Concentrations in the Ythan Estuary

A one-dimensional salt intrusion model is used to investigate the hydrography of the Ythan estuary, a small shallow macrotidal estuary in the north-east of Scotland. The model simulates the longitudinal distributions of water level, salinity and total oxidized nitrogen (TON) in the estuary. The model employs upstream differencing and the Smolarkiewicz anti-diffusion scheme to avoid the numerical difficulties typically encountered when modelling strong tidal flows using centred differences. The physical mechanisms driving the simulations are the tide at the entrance to the estuary and freshwater discharge at the head. The model was calibrated against measurements of water level made at three locations in the estuary, salinity observations made at a central platform and axial salinity distributions. At both spring and neap tides, the full range of salinity observed at the central platform was simulated. However, at the midway stage between springs and neaps, the simulated peak salinity was less than that observed. This was probably due to the sensitivity of the model to the digitisation of the estuarine bathymetry.The model successfully simulated salinity distributions for periods of high and low river flow, and was used to illustrate how TON concentrations fluctuated in response to variations in river flow. The potential implications of variations in the bathymetry of the estuary on salinity and nutrient distributions were predicted to be slight. However, the four fold increase in riverine TON concentrations that has occurred over the past 30 years was shown to increase TON distributions along the entire length of the estuary. The calculated estuary flushing time was strongly dependent on river flow and varied between 11–60 h.     

A modelling study of residence time in a macro-tidal estuary

This paper outlines a numerical modelling study to predict the average residence time of a conservative tracer in a macro-tidal estuary, namely the Mersey Estuary, UK. An integrated hydrodynamic-dispersion model was used to predict the average residence time in the estuary for various tidal level and freshwater discharge conditions. The numerical model was verified against six sets of field measured hydrodynamic data, with the model-predicted water elevations and salinity levels generally agreeing well with the field measurements. The numerical model results show that in the Mersey Estuary both the tidal level and river discharge affect significantly the predicted average residence time. The value of the average residence time is also shown to be closely linked to the intensity of the residual tidal current. This is due to the fact that a large proportion of the Upper and Inner Estuary dries out during low tides, thus a significant amount of the tracer material is transported through the deep channels. An increase in the freshwater discharge causes a considerable increase in the intensity of the residual current along the main channels and thus a reduction in the average residence time. The predicted overall tracer residence time for the whole estuary is relatively short for a relatively large estuary, ranging from less than 1 day to 4 days for various tidal level and freshwater flow combinations. When the tidal range and freshwater discharge are both small, then the local tracer residence time in the upper part of the estuary can be significantly longer than the values predicted for the middle and lower reaches of the estuary.     

Salt-wedge propagation in a Mediterranean micro-tidal river mouth

The dynamics of a seasonally formed salt-wedge propagating along the micro-tidal channel of Strymon River estuary, Northern Greece, and its consequences on river water quality, are thoroughly studied through intensive sampling campaigns. The wedge is developed at the downstream river part, under the summer limited freshwater discharge conditions (Q < 30 m³/s). The geometric features of the wedge (length and thickness) appeared directly related to Strymon River discharge. A maximum intrusion length of 4.7 km along Strymon River estuary was observed under minimum river discharge of almost 6 m³/s. Relations produced from in situ data illustrate that limited river flow expands the wedge horizontally, reducing its vertical dimension, while higher flows lead to increased wedge thickness. Estuarine flushing time ranges between 0.2 and 1.5 days, exponentially dependent on Strymon River discharge. Wedge velocities depicted tidal asymmetry between tidal phases, with consistent inward motion, even under the ebb tidal stage. Strong vertical stratification prevails throughout the tidal cycle, proving the limited vertical mixing between the two layers, although higher interfacial stresses are produced in ebb. Bottom topography plays an interesting role in wedge propagation, as the presence of an underwater sill either prevents saline intrusion during flood or isolates the front of the wedge from its core at the ebb. Ecological consequences of salt-wedge propagation in Strymon River estuary are the frequent evidence of bottom hypoxic conditions and the increased TSS levels, leading to the occurrence of a turbidity maximum at the tip of the salt-wedge. Higher BOD and ammonium levels were mostly observed at the river end, associated to point and non-point pollution sources. Nitrates and silicates were found associated with freshwater fluxes, while ammonia levels were related to saline intrusions. The reduced phosphorus freshwater fluxes, resulting from phosphorus uptake at the upstream reservoir (Kerkini Lake) and the increased bottom turbidity induced by the salt-wedge seem responsible for the limited chlorophyll-a levels along Strymon River estuary.     

Investigation of saltwater intrusion and salinity stratification in winter of 2007/2008 in the Zhujiang River Estuary in China

Saltwater intrusion is a serious environmental problem in the Zhujiang River Estuary(ZRE),which threatens the water supply of fifteen million people.The hydrological observations as well as meteorological and tidal forcing in the winter of 2007/2008 were analyzed to examine the saltwater intrusion in the ZRE.The observational results suggest that the maximum vertical difference of salinity can reach 10 in the Humen Channel during neap tide,but is very small in the Hengmen Channel.The vertically averaged salinity from time series stations during spring tide is higher than that during neap tide.A three-dimensional finite difference model was developed based on the environmental fluid dynamic code(EFDC) to study the mechanism of saltwater intrusion and salinity stratification in the ZRE.By analyzing the salt transport and the temporal variation of saltwater intrusion,the authors found that the net salt transport due to the estuarine circulation during neap tide was more than that during spring tide.This caused salt to advance more into the estuary during neap tide.However,saltwater intrusion was stronger during spring tide than that during neap tide because the spring-neap variation in salt transport was small relative to the total length of the saltwater intrusion.The physical mechanism causing this saltwater intrusion was investigated by a series of sensitivity experiments,in order to examine saltwater intrusion in response to river discharge and winds.The freshwater source was a dominant influencing factor to the saltwater intrusion and controlled salinity structure,vertical stratification and length of the saltwater intrusion.The prevailing northeast monsoon during winter could increase the saltwater intrusion in the ZRE.Though the southwest wind was unfavorable to saltwater intrusion during spring tide,it could increase stratification and saltwater intrusion during neap tide.     

钱塘江河口盐度数值模拟

强潮河口盐水入侵对饮用水源地危害极大。基于平面二维水动力盐度模型, 对典型强潮河口—钱塘江的水动力及盐水入侵过程进行了数值模拟研究。结果表明枯水径流时盐度变化与潮位过程曲线类似, 潮差对盐度大小影响显著, 径流量的增加将逐渐减小其相似程度。当流量增加到一定程度后, 继续增加的一定径流量所产生的抑咸效果减弱, 水资源有效利用率降低, 此时允许水源地盐度超标并改从蓄淡避咸水库取水可有效节约水资源。盐度平面分布显示, 盐水入侵在强潮河口弯道处受涨潮流主流线影响明显, 靠近主流线一岸的盐度大于对岸, 单从盐水入侵角度考虑, 强潮河口弯道段的取水口应设置在远离涨潮流主流线一岸。钱塘江河口盐度数值模拟对于研究减轻盐水入侵对水源地危害的措施具有指导意义。     

The principal factors contributing to the flux of salt in a narrow, partially stratified estuary

Observations of the velocity and salinity structure of the Tees estuary were made at eight stations along the estuary axis between Victoria Bridge and the sea during the summer of 1975. The measurements were made on ten separate tidal periods covering neap and spring tides.The data were collected over a period of relatively low freshwater flows and the residual current was found to have a strong dependence on the Stokes drift. At the upstream stations, the residuals were more than an order of magnitude greater than the currents anticipated from the freshwater discharge. Although the mean stratification decreased as the tidal range increased, the vertical circulation was stronger on spring tides than on neaps. Vertical variations in the amplitude and phase of the tidal current results in a current which strengthens the vertical circulation. However, this effect only made a relatively small contribution to the observed vertical circulation.The relative contribution of the individual salt flux terms to the net upstream transport of salt varies along the estuary. As the estuary narrows, the contribution by the oscillatory terms dominates that from the shear in the steady state flow. Of these oscillatory terms, the correlation of velocity and salinity fluctuations plays a key rôle in the salt transport. The depth mean values make a greater contribution than deviations from the depth mean and the flux due to phase variations over depth is smaller than either of these. Since the Stokes drift is compensated by a down-stream steady state flow, it does not contribute to the tidal mean transport of salt.At the seaward end of the estuary, the salt fluxes due to the steady state vertical shear and the convariance of the tidal fluctuations act in a complementary way to counter the seaward transport of salt by the freshwater flow. With the possible exceptions of the wide or narrow reaches of the Tees, the longitudinal fluxes of salt due to transverse variations in velocity, salinity and depth and turbulent fluctuations are of secondary importance as contributors to the estuary salt budget.On both neap and spring tides, the computed total salt transports at the Newport and Victoria bridges did not match the values required for a salt balance with the corresponding freshwater flows. These fluxes were probably the cause of the observed downstream displacement of the tidal mean salinity distribution between neap and spring tides.     

One-Dimensional Hydraulic Analysis of the Effect of Sea Level Rise on Salinity Intrusion in the Sebou Estuary,Morocco

Global climate change has resulted in a gradual sea-level rise. Sea-level rise can cause saline water to migrate upstream in estuaries and rivers, thereby threatening freshwater habitat and drinking water supplies. On the other hand, sea-level rise, resulting from thermal expansion of ocean waters and increased melting of glaciers and ice caps, is one of the most apparent and widespread consequences of climate change. This phenomenon has been taken into account in all the Assessment Reports published by the Intergovernmental Panel on Climate Change (IPCC). In this paper, salinity intrusion and intrusion length due to possible sea-level rise in the Sebou estuary (Morocco) was investigated. A one-dimensional hydrodynamic-salinity transport model was used for the simulation of the salinity intrusion and associated water quality, with observed field data being used for model calibration and validation. Additionally, the model validation process showed that the model results fit the observed data fairly well. A coupled gas-cycle/climate model was used to generate the climate change scenarios in the studied area that showed sea-level rises varying from 0.3 to 0.9 m for 2100. The models were then combined to assess the impact of future sea-level rise on the salinity distribution and intrusion length in the Sebou estuary. The response of salt intrusion length to changes in important dimensional parameters are presented, showing that the salinity intrusion length is inversely correlated with the river discharge, i.e., a high river discharge results in a reduced salt intrusion and vice versa, and directly with the sea-level rise. Additionally, the magnitude and frequency of the salinity standard violations at the two pump stations were predicted for 2100, showing that the salinity violations under climate change effects can increase to ～45–48% of the times at these locations. Finally, the main objective of this simulation method is to accelerate and facilitate of systems' behavior learning in the current and future situation.     

Seawater Intrusion in Semi-Closed Convergent Estuaries (Case Study of Moroccan Atlantic Estuaries): Application of Salinity Analytical Models

Salinity is an important parameter influencing the water quality of estuaries, and can constitute a serious problem to society due to the need for freshwater for industry and agriculture. Therefore, the determination of salt intrusion length in estuaries is a challenge for managers as well as scientists in this field. The managers tend to use simple and reliable tools for salinity variation. Although 2-D and 3-D numerical models are common tools for the prediction of salinity intrusion now, analytical models of salinity variation are much more efficient, and also require minimal sets of river data. In this paper, two analytical solutions, Brockway and Savenije used worldwide to assess longitudinal salinity variation in alluvial estuaries, are applied to the Moroccan Atlantic semi-closed estuaries, i.e., Sebou and Loukkos. The solutions are derived from salt convection-dispersion equations, with different assumptions for the dispersion coefficient. The estuaries' bathymetry is described by an exponential function. The performance of these two solutions was evaluated by comparing their results with field-measured salinity data. The Brockway model's salinity predictions performs well to observations especially in downstream reaches of the two estuaries (Sebou: R² = 0.95, root mean square error [RMSE] = 1.50‰, normalized root mean square error [NRMSE] = 3.45‰; Loukkos: R² = 0.95, RMSE = 1.13‰, NRMSE = 3.01‰), while the Savenije model outperformed the Brockway's model and is better to predict salt intrusion length and salinity variation along the two estuaries (Sebou: R² = 0.97, RMSE = 1.15‰, NRMSE = 2.85‰; Loukkos: R² = 0.98, RMSE = 0.95‰, NRMSE = 1.94‰). This revealed that both analytical solutions apply well to the estimation of salinity variation and the prediction of salt intrusion in these two estuaries.     

长江口北支进入南支净盐通量的观测与计

根据 2 0 0 1年 4月 1 0— 1 3日长江口大潮期 5个潮周期 3条测量船的同步连续观测资料 ,计算了长江口北支进入南支的净盐通量为 5 4 5× 1 0 6 t,这一结果为预测长江口南支及青草沙水源地的咸潮入侵强度和开发利用长江口淡水资源提供了重要数据。     

Modeling residual circulation and stratification in Oujiang River estuary

A 3D,time-dependent,baroclinic,hydrodynamic and salinity model was implemented and applied to the Oujiang River estuarine system in the East China Sea.The model was driven by the forcing of tidal elevations along the open boundaries and freshwater inflows from the Oujiang River.The bottom friction coefficient and vertical eddy viscosity were adjusted to complete model calibration and verification in simulations.It is demonstrated that the model is capable of reproducing observed temporal variability in the water surface elevation and longitudinal velocity,presenting skill coefficient higher than 0.82.This model was then used to investigate the influence of freshwater discharge on residual current and salinity intrusion under different freshwater inflow conditions in the Oujiang River estuary.The model results reveal that the river channel presents a two-layer structure with flood currents near the bottom and ebb currents at the top layer in the region of seawater influenced on north shore under high river flow condition.The river discharge is a major factor affecting the salinity stratification in the estuarine system.The water exchange is mainly driven by the tidal forcing at the estuary mouth,except under high river flow conditions when the freshwater extends its influence from the river’s head to its mouth.     

椒江河口高混浊水混合过程分析

根据１９９１年洪季的实测资料分析了高度浑浊的椒江河口的混合过程,并探讨了水动力学和沉积动力学因素对河口混合的重要作用,调查研究表明,椒江河口最大浑浊带下的高浑浊水－浮泥层厚达１ｍ,高浑浊水－浮泥层与上覆水之间是泥跃层,泥跃层与高混浊水－浮泥层对水体稳定的作用比同期观测到的盐跃层大１７倍以上,当高浑浊水－浮泥层被侵蚀时,在高浑浊水－浮泥层中的低盐水体又增加了水体的垂向混合能力。     

Morphology, hydrography and sediment dynamics in a mangrove estuary: The Konkoure Estuary, Guinea

The Konkouré Estuary in the Republic of Guinea is a poorly understood atypical mangrove system. Sediment dynamics in tropical estuaries are controlled by a combination of processes including river discharge, morphology, salinity, erosion and deposition processes, the settling of mud, physico-chemical processes and mangrove swamps. Here we present a consistent set of data aimed at characterising the estuary and thus, increasing our understanding of tropical systems, as well as studying the impact of human intervention in the region. Water elevations, current measurements, salinity, suspended sediment concentrations, bathymetry and sediment cover are presented following a 3 year survey of the Konkouré Estuary. Here we provide conclusive evidence that the Lower Konkouré is a shallow, funnel shaped, mesotidal, mangrove-fringed, tide dominated estuary, well mixed during low river discharge. The estuary becomes stratified during high river flows and spring tides whereas a salt wedge appears during neap tides. The Konkouré Estuary has been described as hypersynchronous, and has three terminal outlets, two of which are landward-directed, attesting to a tidal pumping effect, while the third one is seaward-directed, and is controlled by the mangrove. The suspended matter is transported by the tidal effect within the middle estuary and is therefore trapped in the Turbidity Maximum zone (TMZ). The location of the TMZ is river-controlled and is correlated with residual currents but not with salinity front. A dam, constructed 130 km upstream, impacts on the hydrodynamics, and reduces the salinity intrusion by about 25%. It causes an increased low river discharge whereas its efficiency over high river flows is unclear.     

A note on salt intrusion in funnel-shaped estuaries: Application to the Incomati estuary,Mozambique

Salt intrusion in estuaries is important for ecological reasons as well as water extraction purposes. The distance salt intrudes upstream depends on a number of factors, including river discharge, tidal and wind mixing and gravitational circulation. In this paper, an analytical solution is presented for the salt intrusion in a well mixed, funnel-shaped estuary whose cross sectional area decreases exponentially (with decay coefficient β) with distance, x, inland, and in which longitudinal mixing is constant along the length of the estuary. The solution predicts that a graph of the logarithm of salinity against exp (βx) should be a straight line, with slope proportional to the mixing coefficient K_x. The solution is tested against observations from 15 surveys over a four-year period in the Incomati estuary. Good straight line fits, as predicted, are observed on all surveys, with a mean R² = 0.97. The average value of K_x for all surveys is 38 m² s⁻¹. The solution is used to make predictions about the minimum river flow required to prevent salt intruding to an extent where it causes a detrimental effect on water extraction. The minimum recommended river flow required to prevent this is 35 m³ s⁻¹. In recent years, flow has fallen below this level for several months each year.     

长江口整治工程对盐水入侵影响研究

根据实测资料分析了长江口的盐水入侵问题。采用调和常数得到外海控制潮位,用流量控制上游边界,建立了长江口、杭州湾及邻近海域正交曲线坐标系下的二维潮流和盐度数学模型。模型验证了长江口洪、枯季时大、中、小潮的潮位、流速、流向和盐度,较好地模拟了口外顺时针旋转流和口内往复流的特征,反映了外海盐水入侵和北支盐水倒灌的运移特性。在此基础上对长江口综合整治规划方案进行了研究,讨论了整治工程对减轻长江口盐水入侵的作用。     

海南岛南渡江河口的盐水入侵

南渡江是海南岛最大的河流,河口长度较短(约25 km),口外濒临琼州海峡。南渡江的盐水入侵近年来呈加剧趋势,急需对其过程与机理进行研究。本研究运用FVCOM模型,采用2009年枯季的实测数据对模型进行了验证。根据模型计算结果分析了枯季和洪季南渡江河口的流速和盐度的时空分布及盐水入侵的变化,探讨了河口不同位置驱动盐分向陆输运的机制。研究结果表明:枯季小潮时的河口环流强度大于大潮时,而洪季河口环流在大潮时更加发育,洪枯季都表现出大潮期的盐度分层(表底层盐度差)大于小潮期;口门附近,潮汐振荡输运在总的向陆盐分输运通量中占主导,而向上游方向,稳定剪切输运则表现得更重要。主河道内损失的盐分主要通过海甸溪的盐分输入进行补充。     

One-Dimensional Unsteady Analytical Solution of Salinity Intrusion in Estuaries

Based on the one-dimensional salinity transport equation with constant diffusion coefficient, and separated water flow velocity into runoff and tidal current with the single-frequency in an idealized estuary, the simplest unsteady analytical so- lution of salinity intrusion is deduced and the estimation formula of diffusion coefficient is obtained in this paper. The unsteady solution indicates that salinity process in estuaries results from the interaction of runoff and tidal current, and its amplitude is in direct proportion to the product of the velocity of runoff water and the amplitude of tidal flow velocity and in inverse proportion to the diffusion coefficient and the tidal angular frequency, and its phase lag tidal flow with 7/2 which reveals the basic features of the maximum salinity appearing after flood slack and the minimum salinity appearing before ebb slack under the effect of runoff （the advance or lag time is relative to the magnitude of runoff and tidal flow）. According to the measured flow velocity and salinity data, the salinity diffusion coefficient could be estimated. Finally, with the field data of observing sites on the deepwater navigation channel of the Yangtze Estuary, the diffusion coefficient is calculated and a comparative analysis of simulated and measured of salinity process is made. The results show that the solution can comprehensively reflects the basic characteristics and processes of salinity intrusion under the interaction of runoff and tidal flow in estuaries. The solution is not only suitable for theoretical research, but also convenient for estimating reasonable physical parameters and giving the initial condition in the salinity intrusion numerical simulation.