Next Article in Journal
Compliance with Standard Precautions among Clinical Nurses: Validity and Reliability of the Italian Version of the Compliance with Standard Precautions Scale (CSPS-It)
Previous Article in Journal
The Relationship between Nkx2.1 and DNA Oxidative Damage Repair in Nickel Smelting Workers: Jinchang Cohort Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Saltwater Intrusion Function and Preliminary Application in the Yangtze River Estuary, China

1
Department of Hyraulic Engineering, Tsinghua University, Beijing 100084, China
2
China Institute of Water Resources and Hydropower Research, Beijing 100038, China
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2019, 16(1), 118; https://doi.org/10.3390/ijerph16010118
Submission received: 26 November 2018 / Revised: 23 December 2018 / Accepted: 27 December 2018 / Published: 4 January 2019
(This article belongs to the Section Environmental Science and Engineering)

Abstract

:
More attention has been paid to saltwater-intrusion-related problems in recent years. In this research study, a saltwater intrusion function in the Yangtze River Estuary (YRE) was constructed based on the theory of the interactions between energy accumulation and impedance. A MIKE21 model was used to simulate the hydrodynamics of the YRE. Then, through the analysis of the relationships between the river discharge conditions, tidal ranges, and saltwater intrusion, it was determined that, under certain river discharge conditions, the tidal ranges and salinity levels at the stations in the southern branch (SB) of the YRE conformed to S-shaped curve characteristics. Also, the tidal ranges and salinity excessive area rate (SEAR) displayed similar characteristics. Furthermore, the river discharge conditions were also found to match the S-curve characteristics between the two aforementioned relationship features. Therefore, the saltwater intrusion function of the YRE was constructed based on the previously mentioned development rules. Also, the applied quantification methods were elaborated, and the values of the parameters were determined. As a result, the critical river discharge (more than 10,000 m3/s) was obtained, which could withstand large-scale saltwater intrusions. When the river discharge was greater than 30,000 m3/s, the area was considered to be basically without salt water intrusions, and the estuarine ecology was in an optimal state. The saltwater intrusion losses from 2005 to 2015 are also calculated. These findings have important reference value for water dispatching of the YRE in the dry season.

1. Introduction

The estuarine area of the Yangtze River has become an area of rapid economic development and progress due to the abundance of freshwater and ecological resources and its advantageous geographical position. Meanwhile, as human activities are growing, the use of water resources continues to increase. Moreover, the amounts of discharge from many rivers have been observed to have decreased [1]. In particular, the duration of the low flows during the dry seasons has increased, with some rivers even drying up. As a result, the long-term balance of conventional salt intrusion has been damaged, and abnormal phenomena continue to increase [2]. It has been observed that there are increased diffusion regions, a higher probability of occurrence, longer durations, and earlier occurrence times. Therefore, the saltwater intrusions have become increasingly serious, with negative effects on the production and lives of the population in the area, and pose threats to the fragile ecological environment of the Yangtze River Estuary (YRE) [3]. In this study, the goal was to enhance the research regarding saltwater intrusion mechanisms and control factors in order to achieve results of theoretical significance and practical need [4]. In recent years, due to the construction of upstream river hydropower stations and reservoirs, as well as the rebuilding of the estuaries, increases in water levels along the river, and global climate change, the river runoff has displayed a definite decreasing trend, and salt water intrusions have increased [5]. Previous studies were mainly based on numerical simulation and influencing factors, and the research methods have included numerical simulation methods, theoretical analyses, and statistical models [6,7,8]. Also, mathematical modelling has been carried out using one-dimensional, two-dimensional, and three-dimensional models of the salinity spatial distribution, along with physical process simulations of the influencing factors [9,10,11]. Theoretical analyses have been used in the physical model testing methods, and experimental data have been used for the theoretical analyses in order to reveal the mechanism of the intrusions [12]. The statistical models were based on the combined effects of the runoff and tidal currents on the salinity in the estuaries, and multiple regression models for the salinity and runoff and tidal current changes were subsequently established [13,14].
In recent years, the frequency and intensity of saltwater intrusion in the YRE have continuously increased. Previous related studies have included numerical simulations of the changes in salinity in the YRE [15,16,17]. Hui et al. studied the impact of saltwater intrusion on water resources [18]. These have included the “Study of Seawater Intrusion and Utilization of Freshwater Resources in the YRE” [13] and the “Study of the Intrusion of Saline Water at Extremely Low Water Runoff Levels in the YRE during Extreme Drought Years” [19]. It has been established that the tidal and river flows are the most important factors affecting the YRE. Chen Jiyu put forward that the critical flow of saltwater intrusion seriously affecting freshwater resources is 1000 m3/s in the southern branch (SB) of the YRE [20]. Tong used mathematical models to reveal the potential salt water intrusion on the surface and bottom of the Qingcaosha Reservoir estuary [21]. The abovementioned studies often performed numerical simulations in order to analyze the processes of saltwater intrusions. Therefore, there continues to be a lack of research on the combination of numerical simulations and physical mechanisms of the saltwater intrusions in the YRE.
In this study, a saltwater intrusion function of the Yangtze River estuary was constructed based on the theory of the interactions between energy accumulation and impedance. The two main hydrological factors of the saltwater intrusion, tidal range and discharge, were then analyzed. The diffusion processes of the saltwater intrusions that were formed by the aforementioned two factors were derived. Then, the functional relationships of the diffusion of the saltwater intrusion were established, which confirmed the S-shaped curve of the diffusion process. It was determined that the river discharge during the dry season in the Yangtze River estuary is not less than 10,000 m3/s. The index has important reference value for water dispatching in the middle and upper reaches of the Yangtze River.

2. Materials and Methods

2.1. Study Site

The Yangtze River is the largest river in Asia. Under the influence of flow and current, it grows the river regime of the YRE [22], as shown in Figure 1. Located in the upper reaches of the Yangtze River, the Three Gorges Dam (TGD) is the largest dam in the world [23]. The impoundment of the Three Gorges Project caused the saltwater intrusion cycle in the YRE to advance, and the intrusion intensity increased. Due to the regulation and storage of the TGD, the river discharge increased in winter, and the saltwater intrusion was slowed to some extent.
Due to the influence of the YRE, the salinity of the northern branch (NB) has been high for a long time [24]. The tide in the NB may rise to Chongming Island and enter the SB [25]. Along with the decrease in the flow into the sea, the high-salinity saltwater of NB is poured into the SB (Figure 1) [26]. The lower reaches of the SB are affected by tidal currents, resulting in saltwater intrusion. In recent years, with the increasing water intake and related water conservation construction projects, the frequency and intensity of saltwater intrusion into the YRE have increased, which poses a threat to drinking water and causes economic losses [19,20].

2.2. General Description

2.2.1. Hydrological Factors

The average annual flow value of the YRE is 9.24 × 1011 m3 [16]. The average annual inflow rate is 27,856 m3/s, and the flow data from 1957 to 2015 of Datong Hydrological Station are adopted (Figure 2). The average runoff during the dry season is 16,542 m3/s.
The TGD is China’s largest hydropower project. It began to store water on 1 June 2003. Its construction has changed the distribution of the flow into the sea [27,28]. It can be seen from Figure 2 that, after the completion of the TGD, the flood season flow of the Datong Hydrological Station decreased slightly in the flood season, with an average discharge of 39,858 m3/s, and the average runoff in the dry season increased, with an average discharge of 16,867 m3/s.
The YRE is considered to be a meso-tidal estuary. It has a half-day tidal cycle, with two climaxes and low ebbs in a lunar day (approximately 24 h and 50 min) [29]. This study selected the Takahashi Station as a representative tidal station, and subtracted the highest tidal level and lowest tidal level adjacent to each other on a daily basis to obtain the two tidal range differences. The tidal range was taken as the tidal intensity of the current day. The daily tidal range of the year changed with time. Gaoqiao Station is located in the middle of the SB in the YRE. The data regarding the tide levels was determined to be complete, and objectively reflected the changes in the tidal levels in the SB of the YRE.
The variations in the daily tidal ranges over time at the tidal stations indicated that the tidal ranges had significant half-month periodic variations, along with changes during the year [30]. The tidal range of Gaoqiao Station was determined to be between 1 and 4 m (As shown in Figure 3).

2.2.2. The Mechanism

In a general sense, the origin and development of any event are the accumulation of energy due to the growth of quantity. To a certain extent, it must be released. At the same time, in a general sense, the impact of resistance is always present during the development of an event. The accumulative process does not present a uniform linear growth law, but instead exhibits a nonlinear growth law under the influence of the impedance effect, and experiences three stages: initial, development, and end. The origin, development, and ending of the event are actually the result of the action of two opposing factors: one is the driving force that drives the accumulation and release of energy, and the other is the strength of resistance that hinders the development of the two events [31].
When there is no resistance (an ideal environment), the saltwater intrusion intensity increases into a J-shaped curve with an increase in the grade of tide. However, this ideal situation is rarely found in nature. Generally, in an imperfect environment where environmental resistance is present (river discharge), the saltwater intrusion damage is offset by resistance to a certain extent. Therefore, the saltwater intrusion intensity can be characterized as an S-shaped curve (Figure 4) with a slow increase in the early stages, an accelerated increase in the intermediate stages, and stability in the final stages.

2.3. Function Curve Construction

2.3.1. Function

The saltwater intrusion process was observed to conform to an “S” pattern of change, and was often subtle and sudden. Also, the mutation process was found to be obvious. It was determined to be similar to the loss change processes during water disasters (Figure 5) [31]. Therefore, the formula for constructing a saltwater intrusion process in this study was as shown in Formula (1).
f ( X ) = S 2 ( e 2 α ( X X D ) 1 e 2 α ( X X D ) + 1 + 1 )
where X represents the driving or impedance factor during the saltwater intrusion process; S denotes the maximum limit value during the saltwater intrusion process (selected saltwater intrusion area); the parameter α is the vulnerability coefficient; and X B represents the finish point in the saltwater intrusion process. It was determined in the saltwater intrusion process function that the higher the vulnerability coefficient α was, the greater the salinity exceeded the standard area under the same conditions; and X D is the central symmetry point of the curve. Also, the same X A and X B denote the critical values of the salinity.

2.3.2. Parametric Mathematics

In accordance with the function curve of the saltwater intrusion process shown in Figure 5 and Formula (1), the main parameters that affected the entire curve change process were the vulnerability coefficient α , the maximum saltwater intrusion area K , and the critical salinities X A , X D , and X B .
Then, according to the nature of the curve, when X = X D , the curve slope obtains the maximum value, or k m a x = α S . Therefore, the function slope had a correlation with the coefficient α, and the higher the system vulnerability was, the faster the loss change rate would be, and vice versa. At this point, from the symmetry of the S curve and the third derivative of the function at the inflection point, the critical salinities S A ,   S B ,   and   S C and the vulnerability coefficient α could be calculated as follows:
X A = X D 0.66 α ;   X B = X D + 0.66 α ;
X D = X A + X B 2 ;   α = 1.32 X B X A .
From the parameter calculation formulas, it can be seen that any two parameters that are known in the four parameters S A , S B , S C , and α could be used to solve the other two parameters. However, it should be pointed out that, although X D can be theoretically solved by X A and X B , the actual situation is more complicated. Therefore, this method only represents a simplified treatment of the practical problems.

2.4. Model Creation

MIKE21 was used to create the hydrological model. Details of MIKE21 and the model’s creation are given in Xu et al. [32].

2.5. Model Calibration and Validation

The hourly tidal level and salinity at eight stations were collected to verify the simulation (positions shown in Figure 1). The simulated water levels of the tidal stations (Baimao (Bm), Santiaogang (Stg), Yanglin (Yl), and Hengsha (Hs)) are in excellent agreement with the measured values (Figure 6).
The salinity was measured as chlorination, and the Knudsen salinity formula [33] was used for the conversion. The verification results (Z3, Z6, Y3, and Y7) are shown in Figure 7.
The model was quantitatively assessed using the skill model [34]. The model is:
s k i l l = 1 i = 1 N | M X | 2 i = 1 N ( | M X | + | M X | ) 2
where M is the computed values; X is the measured values; X is the mean value of X; and N is the number of in situ observed data. A skill value of 1 indicates that the simulation effect is exactly the same, a skill value between 0.65 and 1 means excellent, a skill value in the range of 0.5–0.65 means very good, a skill value in the range of 0.2–0.5 means good, and a skill value less than 0.2 means bad.
A time series of hourly tidal level and salinity at eight stations (positions shown in Figure 1) spanning the period from 1 March 2002 to 10 March 2002 was collected to validate the simulations. The simulation results at Bm, Stg, Yl, Hs, Z3, Z6, Y3, and Y7 are in excellent agreement with the observations (Figure 6 and Figure 7). The skill values are 0.93, 0.91, 0.94, 0.91, 0.70, 0.82, 0.78, and 0.77 at the eight stations (Table 1).
By verification, all stations’ skill values are greater than 0.65, which is excellent. Among them, four tidal station fitting values are above 0.9 (Bm, Stg, Yl, and Hs), and the salinity station skill values range from 0.70 to 0.82 (Table 1). This is because the salinity simulation depends largely on the flow rate change. The direction and intensity of seawater flow affect the diffusion of salinity [35]. Some of the upstream salinity stations have low salinity and are very sensitive to water flow. There will be some deviation during the simulation.

3. Results

3.1. Analysis of Salinity Distribution in the South Branch

Because the Yangtze River estuary reservoirs are distributed in the south branch, and the north branch salinity has been in the state of exceeding the standard for a long time, this paper mainly studies the change of salinity in the south branch.

3.1.1. Salinity Changes in Different River Discharges

The river discharges of 10,000 m3/s, 15,000 m3/s, 20,000 m3/s, 25,000 m3/s, 30,000 m3/s, and 35,000 m3/s were modelled. We did not change the other conditions. The study focused on salinity changes at eight stations (S1, S2, S3, S4, S5, S6, S7, S8), which are shown in Figure 1.
Table 2 and Figure 8 show the average salinity values of the SB under different flow conditions. It can be seen that, with the decrease of the river discharge levels, the salinity distribution of the SB in the YRE displayed a “high-low-high” trend from the top to the bottom. The “high” point shown in the upper segment was located at the S2 point near the north-south branch junction; the “low” point was located in the middle segment S3–S5. As the salinity continued to rise, the “high” point in the lower segment was located near point S1 in the mouth of the river, and was mainly influenced by the offshore tidal flows.
Figure 8 shows that there was a negative correlation between the flow into the sea and the salinity levels measured at the stations. With increased but low flow rates, the salinity at each station decreased a little. In this study, prior to the flow rates in the sea reaching 10,000 m3/s, the salinity levels at each site were measured. It was found that, when the degree exceeded the standard, a full saltwater intrusion had occurred in the southern branch. Furthermore, when the flow rate into the sea was between 10,000 and 30,000 m3/s, the salinity values of the stations decreased rapidly. When the river discharge reached more than 30,000 m3/s, it was found that all of the measuring points (with the exception of point S2) were low. Therefore, it was confirmed that no salt water intrusions occurred when the salinity critical value was below 0.45.

3.1.2. Salinity Changes in Different Tidal Ranges

The average annual dry season flow into the sea of Datong Station was confirmed to be 16,542 m3/s. A model of the flow boundaries was used to simulate six different tidal ranges (Table 3). The average salinity values of eight stations at different tidal ranges were then examined.
Table 3 and Figure 9 show the average salinity distribution at the southern branch of the YRE under the different tidal ranges. Figure 1 shows the changes in average salinity value of different stations, along with changes in tidal ranges and salinity values at the different stations. As can be seen from the chart, the salinity value of stations also increased when the tidal ranges increased. It can also be seen in the figure that, along the southern branch (from top to bottom), the stations generally showed rising trends. This was due to the fact that the measuring point S2 was affected by the salt water intrusions of the northern branch, resulting in a higher salinity than the nearby stations. Furthermore, the salinity value of the SB1 station in the upper part of the southern branch was lower than 0.45‰ in all cases, and the salinity value of the S1 station was higher than the other stations.
Figure 9 shows that there is a positive correlation between tidal range and salinity at every station. With the increase of tidal range, the salinity of station slowly increases, and the salinity intrusion intensity is very weak before the tidal range reaches 1.75 m. At the tidal range of 1.75–3.5 m, the salinity at the stations increased rapidly, and the salinity invading rate accelerated. When the tidal range reached above 3.5 m, the rising rate of salinity slowed down and stabilized, and the stations continued to be at high salinity.

3.2. Analysis of Salinity Excessive Area Rate in the South Branch

To further quantitatively study the saltwater intrusion process, the typical river discharge (8000 m3/s, 10,000 m3/s, 15,000 m3/s, 20,000 m3/s, 25,000 m3/s, 30,000 m3/s, and 35,000 m3/s), tidal range (1.5 m, 2 m, 2.5 m, 3 m, and 3.5 m), and the relationship with the excessive Saltwater intrusion rate in the south branch were analyzed. The total area of the south branch of the YRE S (1901.86 km2). According to the model calculation, under the condition of a different runoff into the sea, and the excess area of the south branch S s , the ratio of the exceeding area is:
φ = S s S .
As can be intuitively seen in Table 4 and Figure 10, the SEAR in the SB has increased when the tidal range increases under the same flow conditions. Moreover, at the same tidal range, with the increased flow levels, it was observed that the SEAR in the southern branch had also increased. When the flow rates into the sea were less than 10,000 m3/s, the saltwater began to completely invade the area. However, when the runoff from the sea reached 30,000 m3/s, basically no saltwater intrusion was observed.
Figure 11 details the relationship between the tidal ranges and the SEAR at the different flow rates. It can be seen that the tidal ranges were positively correlated with the SEAR. Under the same flow conditions, the salinity exceedance rate was found to have increased with the increases in the tidal ranges. Also, at the same tidal range, the greater the flow rate into the sea, the less the SEAR of the southern branch. It was observed that when the tidal range was less than 1.75 m, the SEAR was lower than 0.1. However, as the tidal ranges increased, the SEAR rapidly increased when the tidal range was between 1.75 and 3.5 m. When the tidal ranges were higher than 3.5 m, the SEAR tended to become stable, and the salinity exceeding the standard area had reached its maximum level.
Figure 12 shows the relationships between the tidal range under the different river discharge conditions and the SEAR in the southern branch. It can be seen in the figure that the flow into the sea was negatively correlated with the SEAR. Under the same tidal range conditions, as the flow increased, the SEAR increased. When the flow rate was less than 10,000 m3/s, the salinity exceedance rate was observed to become stable at the highest value. As the flow rate increased, it could be seen that the SEAR had rapidly dropped when the flow rate was between 10,000 and 30,000 m3/s. Then, when the flow rate into the sea was higher than 30,000 m3/s, the over-standard salinity tended to become stable, and the SEAR had reached the minimum.

3.3. Function Curve Determination

In this research study, through the coupling of the function change, quantum energy accumulation, and impedance action, the saltwater intrusion function in the YRE was constructed using both measured and simulated data from the YRE. Under the action of the flow series x 1 , x 2 , , x n and the tidal range (t), the function between the saltwater intrusion with the tidal range (t) is shown in Equation (6). Also, the function of the maximum area of saltwater intrusion with the river discharge is shown in Equations (7) and (8).
α ( x i ) = f ( t , x i ) = A ( x i ) 2 × ( e 2 α 1 × ( t t D ) 1 e 2 α 1 × ( t t D ) + 1 + 1 )
F ( t , x i ) = A max 4 × ( e 2 α 1 × ( t t D ) 1 e 2 α 1 × ( t t D ) + 1 + 1 ) × ( e 2 α 2 × ( x ^ i D x ^ i ) 1 e 2 α 2 × ( x ^ i D x ^ i ) + 1 + 1 )
A max = max ( A ( x i ) )
In the equations, x i represents the river discharge; t is the tidal range; A ( x i ) is the maximum area of the salinity intrusion under the x i condition; α 1 is the sensitivity coefficient of coverage area to the tidal range during the salt tide intrusion; t D is the central symmetry point of the curve; α 2 denotes the sensitivity coefficient of the maximum intrusion area’s influencing factors during the salt tide intrusion response process; x ^ i represents the dimensionless value of the inflow into the sea (including the measured river discharge and the average annual river discharge); and x ^ i D is a function of the dimensionless sea inflow and the salt tide intrusion response function curve’s turning point.

3.3.1. Function Determination for the SEAR and Tidal Ranges

In this study, by using multiple parameter calibrations and curve fittings, the sensitivity coefficient α 1 = 2 was successfully determined, and the relevant parameters are shown in Table 5. The curve of the function of the SEAR with the change of the tidal ranges under the different river discharge conditions is shown in Figure 13.
It can be seen from Figure 13 that, as a group, the curve of the saltwater intrusion area with the changes in the tidal ranges was within the same curve family. Its parameters had also followed certain rules. According to Table 5, the relationship between the transition tidal range t D and river discharge is as follows: when the river discharge increases, t D increases. Indirectly, the larger the discharge is, the greater the tide that causes saltwater intrusion is. The side shows that river discharge acts as an impedance in the saltwater intrusion process.

3.3.2. Function Determination for the SEAR and River Discharge Conditions

The accurate determination of the relationship between the runoff from the sea and the maximum area of the saltwater intrusion was the precondition for determining the boundary parameter A ( x i ) in the saltwater intrusion response function. In this study, through the use of multiple parameters, the parameters were obtained as follows: x ^ i D = 0.72, α 2 = 5.0, and A max = 1. Therefore, the function was as follows:
A ( x i ) = 0.5 × ( e 10.0 × ( 0.72 x ^ i ) 1 e 10.0 × ( 0.72 x ^ i ) + 1 + 1 ) .
As can be seen from Figure 14, the function fitting effect is good, and there is an S-type variation pattern between the maximum rate of salinity exceeding the standard and the flow rate into the sea.

3.3.3. Function of the Saltwater Intrusion Process in the YRE

In this study, the relationships between the intrusion areas with the salt tide and the changes in the tidal ranges, and the relationships between the maximum area of the salt tide intrusion and the flow rate into the sea, were used to obtain the theoretical function of the saltwater intrusion, along with the diffusion response process in the estuary, as shown in Equation (8). From the calculation and solution of the measured data detailed in the previous two sections, the parameters that could be determined in the theoretical function were as follows: x ^ i D = 0.72, α 1 = 2.0, α 2 = 5.0, and A max = 1. Therefore, the saltwater intrusion diffusion response function in the YRE could be further transformed into Equation (10).
F ( t , x i ) = A max 4 × ( e 2 α 1 × ( t t D ) 1 e 2 α 1 × ( t t D ) + 1 + 1 ) × ( e 2 α 2 × ( x ^ i D x ^ i ) 1 e 2 α 2 × ( x ^ i D x ^ i ) + 1 + 1 )      = 0.25 × ( e 4 × ( t t D ) 1 e 4 × ( t t D ) + 1 + 1 ) × ( e 10.0 × ( 0.72 x ^ i ) 1 e 10.0 × ( 0.72 x ^ i ) + 1 + 1 )
Among them, t D changes with the change of runoff into the sea, and, according to the relevant parameter values determined in Table 5, the relationship between the dimensionless flow into the sea and the t D can be obtained (see Figure 15). The relationship between the two is brought into Equation (10), and Formula (11) is obtained as follows:
F ( t , x i ) = 0.25 × ( e 4 × ( t ( 1.0809 x ^ i + 2.3352 ) ) 1 e 4 × ( t ( 1.0809 x ^ i + 2.3352 ) ) + 1 + 1 ) × ( e 10.0 × ( 0.72 x ^ i ) 1 e 10.0 × ( 0.72 x ^ i ) + 1 + 1 ) .
Figure 16 details the simulated scatter distribution and function fitting curved surface. It can be seen in the figure that the tidal ranges played key roles in promoting the SEAR, and the river discharge conditions played key roles in resisting the SEAR. The former (tidal ranges) were the driving forces of the saltwater intrusions. The latter (river discharge conditions) were the group resistance to the saltwater intrusions. It was determined that, under any of the flow conditions into the sea, as the tidal range increased, the SEAR consistently changed from a slow to a rapid expansion, and then tended to ease once it approached the maximum area. The three factors (tidal range, river discharge, and SEAR) were found to be related to each other and had formed an S-shaped curved surface that conformed to the principles of the interactions between the amount of accumulated energy and the impedance. The SEAR was observed to have gradually increased with the tidal range above approximately 1.75 m (in which the runoff was less than 10,000 m3/s). Also, after entering the sea with a runoff greater than 30,000 m3/s, it was observed to be difficult for the tidal ranges that exceeded 3.75 m to form an excess area.

3.4. Saltwater Intrusion Damage Assessment

According to some studies in recent years, the value of ecosystem services per unit area in the YRE is between 117.91 and 189.44 million RMB/m2 [25,26]. In this paper, the average value of 1.5 million RMB/km2 is taken for analyzing the loss. Using Formula (9), the maximum intrusion area of the SB in the YRE is calculated by the minimum sea runoff from November to April of the 2005–2015 dry season (from November to April of the following year), and the maximum loss of saltwater intrusion in the YRE is calculated. Since the NB accounts for about 24% of the YRE’s area, and the salinity of the NB’s gate is equivalent to that of the outer sea during the dry season, the salinity of the upper section of Qinglong Port is also higher than the set salinity standard value, so the salt intrusion area of the entire Yangtze River estuary is south. The sum of the intrusion area of the saltwater and the area of the NB is used to evaluate the saltwater intrusion loss in the YRE. The maximum loss in the dry season from 2005 to 2015 is shown in Table 6.

4. Discussion

This research studied the numerical simulation results of saltwater intrusion in the YRE, and studied the saltwater intrusion and diffusion response from the hydrological factors that dominate the saltwater intrusion. In future research, many problems should be explored:
The study of the saltwater intrusion and diffusion process considers the main hydrological factors in the process of saltwater intrusion: river discharges and tides. However, the environment in the estuary area is complex. Under the influence of natural and human factors, the river regime changes, and the estuary’s water conservation construction project, the change of wind stress may have a certain impact on the original saltwater intrusion process, which affects the accuracy of the diffusion response function. The above influencing factors should be gradually introduced in the future. On the basis of studying the influence of the above factors on the salt tide intrusion process, the parameters reflecting the above factors should be introduced to correct the saltwater intrusion diffusion response function.
The paper takes the YRE as an example for empirical analysis and function construction that is based on the hydrological response mechanism of saltwater intrusion. In China, in addition to the YRE, areas such as the Pearl River estuary and the Qiantang estuary are also facing serious saltwater intrusion, as well as in some estuaries of other countries. In future research, we will further expand the empirical analysis into other estuary areas. The applicability of the saltwater intrusion function in other regions should be further studied.
The calculation results of the numerical simulation can be used as the basic data of the saltwater intrusion process to ensure the richness of the research data. In order to improve the accuracy of the saltwater intrusion analysis, it is also necessary to enrich the hydrological data in the estuary area and increase the accuracy of the numerical simulation in order to more accurately evaluate the saltwater intrusion and diffusion process.

5. Conclusions

In this study, based on previous research results regarding the quantum energy accumulation and impedance action theory, the saltwater intrusion processes that occur in the YRE under the interactions of the tide rates and runoff into the sea were examined. It was concluded from the research findings that the tidal water rates were the driving forces for the saltwater intrusions. The impedance of the intrusions was determined to be consistent with the theory. Also, an S-curve function, which conformed to the theory of energy accumulation and impedance, was successfully constructed. The parameter settings and quantification methods were described in detail in the current study. A MIKE21 model was used to simulate the hydrodynamic processes and salinity diffusion in the YRE. First, a linear analysis was performed on eight stations that covered the southern branch of the YRE. The calculation results of the salinity levels from the upstream to the downstream of the southern branch of the YRE were presented. The characteristic of a “high-low-high” change trend was determined to be due to the influences of the southern branch’s high brine concentration offshore, as well as the salinity at the SB2 station near the northern branch being significantly increased. These findings were found to be in line with the research results of previous related studies. A linear analysis of the salinity exceeding standard area (SEAR) of the southern branch was carried out. The exceeding rates of the flow and salinity, and the exceeding rate of the tidal differences and salinity, were found to be in line with the change rule of the previously mentioned S curve.
Also, through the determination of the brine process function, two critical sea inlet flows were determined in this study as follows: When the seawater flows were greater than 10,000 m3/s, which could withstand large-scale saltwater intrusions; and when the seawater flows were greater than 30,000 m3/s, where no saltwater was observed to intrude, and the estuarine ecology was in an optimal state. The saltwater intrusion losses from 2005 to 2015 were also calculated. The results of this study will potentially provide important reference values for future water scheduling during the dry seasons in the upper and middle reaches of the YRE. This important information will further the comprehensive management of the estuary in the future.

Author Contributions

Conceptualization, Z.X., J.M., and Y.H.; Data curation, Y.H.; Writing (original draft), Z.X.; Writing (review & editing), J.M. and Y.H.

Funding

This research was funded by: 1. The National Key R&D Program of China (Grant No. 2016YFC0503502); 2. The National Natural Science Foundation Program of China (Grant No. 41861124006).

Acknowledgments

The authors would like to express their appreciation to the editor and the anonymous reviewers, who provided very helpful comments to improve the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chen, M.; Hu, Y.; Ma, J.; Yan, L.; Zhou, F.; Wang, Y. Diffusion response function and preliminary application of saltwater intrusion in the yangtze river estuary. Adv. Water Sci. 2017, 28, 203–212. [Google Scholar]
  2. Zhang, E.; Savenije, H.H.G.; Wu, H.; Kong, Y.; Zhu, J. Analytical solution for salt intrusion in the Yangtze estuary, China. Estuar. Coast. Shelf Sci. 2011, 91, 492–501. [Google Scholar] [CrossRef]
  3. Bhuiyan, M.J.A.N.; Dutta, D. Assessing impacts of sea level rise on river salinity in the gorai river network, Bangladesh. Estuar. Coast. Shelf Sci. 2012, 96, 219–227. [Google Scholar] [CrossRef]
  4. Yuan, Y.R.; Liang, D.; Rui, H.X. The numerical simulation and analysis of three-dimensional seawater intrusion and protection projects in porous media. Sci. China 2009, 52, 92–107. [Google Scholar] [CrossRef]
  5. Chen, S.L.; Zhang, E.F.; Guo-Chuan, G.U.; Ping, L.I. Salt intrusion in the south passage of Yangtze river Estuary during extraordinary dry year. Mar. Sci. Bull. 2009, 28, 29–36. [Google Scholar]
  6. Uncles, R.J.; Stephens, J.A. The effects of wind, runoff and tides on salinity in a strongly tidal sub-estuary. Estuar. Coasts 2011, 34, 758–774. [Google Scholar] [CrossRef]
  7. Kuang, C.; Chen, W.; Gu, J.; Su, T.C.; Song, H.; Ma, Y.; Dong, Z. River discharge contribution to sea-level rise in the yangtze river estuary, China. Cont. Shelf Res. 2017, 134, 63–75. [Google Scholar] [CrossRef]
  8. Nijs, M.A.J.D.; Pietrzak, J.D. Saltwater intrusion and ETM dynamics in a tidally-energetic stratified estuary. Ocean Model. 2012, 49–50, 60–85. [Google Scholar] [CrossRef]
  9. Kuijper, K.; Rijn, L.C.V. Analytical and numerical analysis of tides and salinities in estuaries; part II: Salinity distributions in prismatic and convergent tidal channels. Ocean Dyn. 2011, 61, 1743–1765. [Google Scholar] [CrossRef]
  10. Ralston, D.K.; Geyer, W.R.; Lerczak, J.A. Subtidal salinity and velocity in the hudson river estuary: Observations and modeling. J. Phys. Oceanogr. 2008, 38, 753. [Google Scholar] [CrossRef]
  11. Chen, W.; Chen, K.; Kuang, C.; Zhu, D.Z.; He, L.; Mao, X.; Liang, H.; Song, H. Influence of sea level rise on saline water intrusion in the yangtze river estuary, China. Appl. Ocean Res. 2016, 54, 12–25. [Google Scholar] [CrossRef]
  12. Lu, C.; Yuan, L.; Gao, S.; Chen, R.; Su, B. Experimental study on the relationship between tide strength and salt intrusion length. Adv. Water Sci. 2013, 24, 251–257. [Google Scholar]
  13. Tang, J.H.; Xu, J.-Y.; Zhao, S.W.; Liu, W.Y. Research on saltwater intrusion of the south branch of the Changjiang estuary based on measured data. Resour. Environ. Yangtze Basin 2011, 20, 677–684. [Google Scholar]
  14. Zheng, X.Q.; Yun, Y.U.; Pan, L.Z. Statistical prediction study of runoff and tide on Changjiang river saltwater intrusion. Mar. Forecasts 2014, 31, 18–23. [Google Scholar]
  15. Wu, H.; Zhu, J.; Chen, B.; Chen, Y. Quantitative relationship of runoff and tide to saltwater spilling over from the North Branch In the Changjiang estuary: A numerical study. Estuar. Coast. Shelf Sci. 2006, 69, 125–132. [Google Scholar] [CrossRef]
  16. Xue, P.; Chen, C.; Ding, P.; Beardsley, R.C.; Lin, H.; Ge, J.; Kong, Y. Saltwater intrusion into the Changjiang river: A model-guided mechanism study. J. Geophys. Res. Oceans 2009, 114. [Google Scholar] [CrossRef]
  17. Li, L.; Zhu, J.; Wu, H.; Guo, Z. Lateral saltwater intrusion in the north channel of the changjiang estuary. Estuar. Coasts 2014, 37, 36–55. [Google Scholar] [CrossRef]
  18. Hui, W.U.; Rong, Z.J. Analysis of the transport mechanism of the saltwater spilling over from the north branch in the changjiang estuary in china. Acta Oceanol. Sin. 2007, 29, 17–25. [Google Scholar]
  19. Zhu, J.R.; Hui, W.U. Numerical simulation of the longest continuous days unsuitable for water intake in the Dongfengxisha reservoir of the Changjiang estuary. J. East China Norm. Univ. 2013, 28, 2807. [Google Scholar]
  20. Chen, J.; Xu, H. Impacts of the Yangtze River Three—Gorge Hydro—Engineering Works on the Yangtze Estuaky. Pesour. Enuiron. Yangtza Valley 1995, 3. Available online: http://en.cnki.com.cn/Article_en/CJFDTOTAL-CJLY503.010.htm (accessed on 27 November 2018).
  21. Tong, C.; Zheng, J.; Zhang, C.; Guilbaud, C. Salinity Response to the Runoff from Yangtze River Basin at Qingcaosha Reservoir Area in Yangtze Estuary. In Proceedings of the ASME 2010 International Conference on Ocean, Offshore and Arctic Engineering, Shanghai, China, 6–11 June 2010; Ocean, Offshore and Arctic Engineering Division: New York, NY, USA, 2010; pp. 247–255. [Google Scholar]
  22. Yu, Y.; Wang, P.; Wang, C.; Wang, X.; Hu, B. Assessment of the multi-objective reservoir operation for maintaining the turbidity maximum zone in the Yangtze river estuary. Int. J. Environ. Res. Public Health 2018, 15, 2118. [Google Scholar] [CrossRef] [PubMed]
  23. Mei, X.; Dai, Z.; Gelder, P.H.A.J.M.; Gao, J. Linking three gorges dam and downstream hydrological regimes along the Yangtze river, China. Earth Space Sci. 2015, 2, 94–106. [Google Scholar] [CrossRef]
  24. Hou, C.; Zhu, J. Study on the quantitative relationship of the location of the tidal current limit and the river discharge in the Changjiang estuary. J. East China Norm. Univ. 2013, 28, 18–26. [Google Scholar]
  25. Sun, Z.; Xin, Y.; Xie, C.; Qi, L.I. An empirical predictive model for saltwater intrusion in the south branch influenced by tidal flow from the north branch in the yangtze river estuary. Adv. Water Sci. 2017, 28, 213–222. [Google Scholar]
  26. Chen, W.; Tang, J.; Lan, X.; Lu, L.I.; Zhao, H. Impact of salinity intrusion of north branch on water quality of intake water areas of dongfengxisha reservoir. Water Resour. Protect. 2013, 1, 011. [Google Scholar]
  27. Dai, Z.; Liu, J.T. Impacts of large dams on downstream fluvial sedimentation: An example of the three gorges dam (TGD) on the changjiang (Yangtze river). J. Hydrol. 2013, 480, 10–18. [Google Scholar] [CrossRef]
  28. Dai, Z.; Liu, J.T.; Xiang, Y. Human interference in the water discharge of the Changjiang (Yangtze river), China. Int. Assoc. Sci. Hydrol. Bull. 2015, 60, 1770–1782. [Google Scholar] [CrossRef]
  29. Luo, X.F.; Chen, Z.C. Numerical simulation of salinity in yangtze river estuary. Hydro-Sci. Eng. 2004, 2, 005. [Google Scholar]
  30. Hui, Y.L.; Rong, Z.J.; Xian, Z.S. 3-D numerical simulation of tide and tidal current fields in the Changjiang estuary, Hangzhou bay and their Adjacent sea. J. East China Norm. Univ. 2001, 3, 74–84. [Google Scholar]
  31. Chen, M.; Ma, J.; Hu, Y.; Zhou, F.; Li, J.; Yan, L. Is the s-shaped curve a general law? An application to evaluate the damage resulting from water-induced disasters. Nat. Hazards 2015, 78, 497–515. [Google Scholar] [CrossRef]
  32. Xu, Z.; Ma, J.; Wang, H.; Hu, Y.; Yang, G.; Deng, W. River discharge and saltwater intrusion level study of yangtze river estuary, China. Water 2018, 10, 683. [Google Scholar] [CrossRef]
  33. Knudsen Formula; Springer: Berlin/Heidelberg, Germany, 2014; p. 770. Available online: https://link.springer.com/referenceworkentry/10.1007/978-3-642-41714-6_110286. (accessed on 27 November 2018).
  34. Willmott, C. On the validation of models. Phys. Geogr. 1981, 2, 184–194. [Google Scholar] [CrossRef]
  35. Yang, T.J.; Wang, Y.G.; Huang, H.M.; Yuan, C.G. Research on salinity intrusion in the yangtze river estuary by numerical simulation. J. Waterw. Harb. 2013. Available online: http://en.cnki.com.cn/Article_en/CJFDTOTAL-SDGK201306003.htm (accessed on 27 November 2018).
Figure 1. The Yangtze River Estuary.
Figure 1. The Yangtze River Estuary.
Ijerph 16 00118 g001
Figure 2. The average river discharge distribution of Datong station. TGD, Three Gorges Dam.
Figure 2. The average river discharge distribution of Datong station. TGD, Three Gorges Dam.
Ijerph 16 00118 g002
Figure 3. The tidal range with time in Gaoqiao station.
Figure 3. The tidal range with time in Gaoqiao station.
Ijerph 16 00118 g003
Figure 4. The salt water intrusion process curve.
Figure 4. The salt water intrusion process curve.
Ijerph 16 00118 g004
Figure 5. A conceptual model of the saltwater intrusion function.
Figure 5. A conceptual model of the saltwater intrusion function.
Ijerph 16 00118 g005
Figure 6. Verification of the tidal levels.
Figure 6. Verification of the tidal levels.
Ijerph 16 00118 g006
Figure 7. Verification of the salinity levels.
Figure 7. Verification of the salinity levels.
Ijerph 16 00118 g007
Figure 8. The average salinity values of the stations under different flow conditions.
Figure 8. The average salinity values of the stations under different flow conditions.
Ijerph 16 00118 g008
Figure 9. The average salinity values of stations at different tidal ranges.
Figure 9. The average salinity values of stations at different tidal ranges.
Ijerph 16 00118 g009
Figure 10. The salinity excessive area of the South Branch in different tidal range and river discharge conditions. Note: (1) Q = 10,000 m3/s, T = 1.5 m; (2) Q = 10,000 m3/s, T = 2.5 m; (3) Q = 10,000 m3/s, T = 3.5 m; (4) Q = 15,000 m3/s, T = 1.5 m; (5) Q = 15,000 m3/s, T = 2.5 m; (6) Q = 15,000 m3/s, T = 3.5 m; (7) Q = 20,000 m3/s, T = 1.5 m; (8) Q = 20,000 m3/s, T = 2.5 m; (9) Q = 20,000 m3/s, T = 3.5 m; (10) Q = 30,000 m3/s, T = 1.5 m; (11) Q = 30,000 m3/s, T = 2.5 m; (12) Q = 30,000 m3/s, T = 3.5 m.
Figure 10. The salinity excessive area of the South Branch in different tidal range and river discharge conditions. Note: (1) Q = 10,000 m3/s, T = 1.5 m; (2) Q = 10,000 m3/s, T = 2.5 m; (3) Q = 10,000 m3/s, T = 3.5 m; (4) Q = 15,000 m3/s, T = 1.5 m; (5) Q = 15,000 m3/s, T = 2.5 m; (6) Q = 15,000 m3/s, T = 3.5 m; (7) Q = 20,000 m3/s, T = 1.5 m; (8) Q = 20,000 m3/s, T = 2.5 m; (9) Q = 20,000 m3/s, T = 3.5 m; (10) Q = 30,000 m3/s, T = 1.5 m; (11) Q = 30,000 m3/s, T = 2.5 m; (12) Q = 30,000 m3/s, T = 3.5 m.
Ijerph 16 00118 g010
Figure 11. The salinity excessive area of the South Branch in different tidal range and river discharge conditions.
Figure 11. The salinity excessive area of the South Branch in different tidal range and river discharge conditions.
Ijerph 16 00118 g011
Figure 12. The relationship between the river discharge and the SEAR in different tidal range conditions.
Figure 12. The relationship between the river discharge and the SEAR in different tidal range conditions.
Ijerph 16 00118 g012
Figure 13. The function curve between the tidal range and the SEAR in different river discharge conditions.
Figure 13. The function curve between the tidal range and the SEAR in different river discharge conditions.
Ijerph 16 00118 g013
Figure 14. The function curve between river discharge and maximum excessive area rate.
Figure 14. The function curve between river discharge and maximum excessive area rate.
Ijerph 16 00118 g014
Figure 15. The relationship between dimensionless river discharge and transition tidal range ( t D ).
Figure 15. The relationship between dimensionless river discharge and transition tidal range ( t D ).
Ijerph 16 00118 g015
Figure 16. The scatter diagram and saltwater intrusion function surface fitting.
Figure 16. The scatter diagram and saltwater intrusion function surface fitting.
Ijerph 16 00118 g016
Table 1. Model–observation data comparison statistics for tidal and salinity levels.
Table 1. Model–observation data comparison statistics for tidal and salinity levels.
QuantityStationsSkillCondition
Tidal currentBm0.93Excellent
Stg0.91Excellent
Yl0.94Excellent
Hs0.91Excellent
SalinityZ30.70Excellent
Z60.82Excellent
Y30.78Excellent
Y70.77Excellent
Table 2. The average salinity values of the southern branch (SB) under different flow conditions.
Table 2. The average salinity values of the southern branch (SB) under different flow conditions.
StationAverage Salinity Values with the Discharge of 16,542 m3/s (‰)Average Salinity Values with Different Flow Conditions.
10,000 m3/s (‰)15,000 m3/s (‰)20,000 m3/s (‰)25,000 m3/s (‰)30,000 m3/s (‰)35,000 m3/s (‰)
S10.050.100.050.040.030.030.02
S21.462.281.701.200.800.550.40
S30.340.760.410.260.190.150.12
S40.731.050.830.610.470.360.26
S50.691.220.800.570.440.330.24
S60.711.270.820.580.420.350.25
S70.801.670.940.650.490.350.28
S81.502.971.751.000.700.450.38
Table 3. The average salinity values of SB at different tidal ranges.
Table 3. The average salinity values of SB at different tidal ranges.
StationAverage Salinity Value at Different Tidal Range (‰)
<1.5 m (‰)1.5–2 m (‰)2–2.5 m (‰)2.5–3 m (‰)3–3.5 m (‰)>3.5 m (‰)
S10.0700.1010.1980.3430.4120.445
S20.3200.7101.3372.5242.8103.270
S30.0980.1700.475 0.7780.8681.072
S40.1810.2610.7911.3101.5401.630
S50.2300.2790.9491.3201.7401.810
S60.2310.3211.1401.5601.8002.100
S70.3100.5321.3102.2302.5702.780
S80.4120.8431.6662.5003.1203.330
Table 4. The South Branch’s salinity excessive area rate (SEAR) in different river discharge and tidal range conditions.
Table 4. The South Branch’s salinity excessive area rate (SEAR) in different river discharge and tidal range conditions.
Tidal Range
(m)
1.522.533.54
River Discharge
(m3/s)
80000.020.110.320.780.960.99
10,0000.010.060.280.740.940.98
15,00000.020.140.480.770.85
20,00000.010.040.180.390.48
25,000000.010.030.080.13
30,00000000.010.02
35,000000000
Table 5. The saltwater intrusion area and tidal changes relationship parameters.
Table 5. The saltwater intrusion area and tidal changes relationship parameters.
Q (m3/s)800010,00015,00020,00025,00030,00035,000
A i max (%)0.990.980.850.480.130.020
t D (m)2.652.752.953.053.23.6-
Table 6. The maximum loss by saltwater intrusion in the Yangtze River Estuary (YRE).
Table 6. The maximum loss by saltwater intrusion in the Yangtze River Estuary (YRE).
YearMinimum River Discharge
(m3/s)
Maximum Excess Area
(km2)
Maximum Loss
(100 Million RMB)
200511,6002062.4015.84
200696502130.8716.37
200710,0002122.2916.30
200810,3002113.8516.23
200910,9002093.4416.08
201011,2002081.1915.98
201113,8001892.0014.53
201211,2002081.1915.98
201312,2652021.4315.52
201411,2712076.5615.94
201512,3772014.7315.47

Share and Cite

MDPI and ACS Style

Xu, Z.; Ma, J.; Hu, Y. Saltwater Intrusion Function and Preliminary Application in the Yangtze River Estuary, China. Int. J. Environ. Res. Public Health 2019, 16, 118. https://doi.org/10.3390/ijerph16010118

AMA Style

Xu Z, Ma J, Hu Y. Saltwater Intrusion Function and Preliminary Application in the Yangtze River Estuary, China. International Journal of Environmental Research and Public Health. 2019; 16(1):118. https://doi.org/10.3390/ijerph16010118

Chicago/Turabian Style

Xu, Zhi, Jing Ma, and Yajie Hu. 2019. "Saltwater Intrusion Function and Preliminary Application in the Yangtze River Estuary, China" International Journal of Environmental Research and Public Health 16, no. 1: 118. https://doi.org/10.3390/ijerph16010118

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

Article Metrics

Back to TopTop