1036 Clean – Soil, Air, Water 2012, 40 (10), 1036–1046
Zhiming Zhang1 Research Article Baoshan Cui1 Xiaoyun Fan1 Kejiang Zhang2 Wetland Network Design for Mitigation of Saltwater Hui Zhao3 Intrusion by Replenishing Freshwater in an Estuary Honggang Zhang1
In this paper, a wetland network is designed to mitigate saltwater intrusion based on 1State Key Joint Laboratory of Environmental Simulation and the relationship between river discharge and salinity in Modaomen waterway (MDMW), Pollution Control, School of the Pearl River estuary (PRE) of southern China. The designed network consists of Environment, Beijing Normal existing and expanded wetlands, river channels and their tributaries. The freshwater University, Beijing, P. R. China stored in wetlands can be diverted into river channels to adjust salinity gradients when 2Xinjiang Research Center of Water and Wastewater Treatment, Xinjiang saltwater intrusion reaches the predefined locations. In the MDMW, two exponential Deland Co., Ltd., Urumqi, P. R. China regression models are established between the freshwater discharge and the saltwater 3Institute of Mountain Hazards and intrusion length at both high water slack (HWS) and low water slack (LWS), respectively. Environment, Chinese Academy of The wetland network can effectively mitigate the saltwater intrusion when the fresh- Sciences, Chengdu, P. R. China water is diverted from wetlands into river channels according to the regressive models. The method of wetland network can save over 50% more freshwater (varying from 46.1 to 56.4% at HWS and from 53.4 to 61.8% at LWS) than the emergent water allocation schedule (EWAS) implemented in 2005 to protect against saltwater intrusion in the same area. Wetland network design facilitates water resource management in the PRE and can also be easily generalized to other estuaries.
Keywords: Discharge–salinity relationships; Ecosystem; Modeling; Storage capacity Received: January 11, 2012; revised: April 16, 2012; accepted: May 22, 2012 DOI: 10.1002/clen.201100735
1 Introduction and recharging freshwater in coastal aquifers [24, 25]. This process is usually time-consuming and negatively influences the surrounding An estuarine ecosystem, which features the land–ocean interactions ecosystems due to the pumped saltwater. With the relationships of an abundant population of flora and fauna, has the highest service between the freshwater discharge and salinity gradient, many values compared to other ecosystems [1]. Meanwhile, it is a sensitive efforts have also been employed to divert large quantities of fresh- ecosystem to the rapid industrialization and economic development water from upstream into estuarine or coastal zones to mitigate and sea level rise [2–5]. With the water shortages and sea level rise, saltwater intrusion. According to these relationships, an emergent researchers found that the salinity gradient was an important water allocation schedule (EWAS) has been implemented to mitigate factor in determining these estuarine ecosystem services [6–14]. saltwater intrusion since 2005 in Pearl River Basin, China [33, 34]. Over the past decades, salinity gradients and freshwater inflow in The frequency and scope of saltwater intrusion are serious due to the an estuary have been significantly disturbed by anthropogenic activi- drought in Pearl River Basin. Saltwater intrusion is a threat to the ties [15–19] and climate change [20]. For example, saltwater intrusion safety of water supply to the surrounding areas. Fortunately, will increase water salinity if less freshwater is discharged into authorities implement the EWAS which tried to divert the upstream estuaries arising from the dam construction projects [15, 16, 21, 22]. freshwater into the tidal rivers during dry season of 2005. There Freshwater in estuaries is also usually intercepted and used for is a significant negative correlation between the river discharge irrigation, maintenance of ecosystem integrity, and domestic and and water salinity in Pearl River estuary (PRE) [35]. The saltwater industrial purposes such that less freshwater is available for miti- intrusion will be mitigated by increased the river discharge, which gation of saltwater intrusion. can ensure water supply to the cities of Zhuhai and Macao. Although From the literature, researches about the saltwater intrusion have this schedule is effective, it is costly due to a lack of sufficient mainly focused on groundwater in coastal aquifers [23–25] and freshwater. surface water in river channels [26–32]. Previous studies indicated The freshwater discharge is significantly associated with the sa- the control effects of saltwater intrusion by pumping out saltwater linity gradient in many estuaries around the word. There are negative correlations between the saltwater intrusion length and river dis- charge [26, 27, 30, 36]. The saltwater intrusion follows a negative Correspondence: Professor B. Cui, State Key Joint Laboratory of power law to the river discharge in northern San Francisco Bay [37]. Environmental Simulation and Pollution Control, School of Environment, Becker et al. [26] compared three regression models (linear, expo- Beijing Normal University, No. 19 Xinjiekouwai Street, Beijing 100875, P. R. China nential, and power-law regressions) to determine the optimal E-mail: [email protected] relationship between the saltwater intrusion length and river dis- charge. Therefore, a simple and effective relationship between the Abbreviations: EWAS, emergent water allocation schedule; HWS, high water slack; LWS, low water slack; MDMW, Modaomen waterway; PRE, freshwater discharge and salinity gradient along the river channel Pearl River estuary should be established.
ß 2012 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim www.clean-journal.com Clean – Soil, Air, Water 2012, 40 (10), 1036–1046 Wetland Network Design for Mitigation of Saltwater Intrusion 1037
Considering the water saving and the mitigation of saltwater 228N–228400 N, 1138 E–1138400 E) is one of the main outlets of the intrusion, how to allocate the upstream freshwater in the estuarine Pearl River (Fig. 1) [28, 50]. and coastal regions is another important question. Poulter et al. [4] The average annual runoff of MDMW is approximately presented a network-based method to identify the major shoreline 8.84 1011 m3 accounting for 28.2% of the total runoff of the entry points of saltwater to the drainage river networks. These Pearl River [46, 51]. The MDMW is a typical river channel with points are useful in the managing and monitoring of saltwater. intensive runoff and neap range in PRE. The great diurnal range The network-based method is nearly natural and cost-effective is approximately 0.86 m [51]. There are six fixed gauging stations [38]. With the application of graph theory to network analysis, (Baiqing, Daao, Zhuyin, Xihekou, Zhupaisha, and Guadingjiao) along wetland network approaches are widely used in watershed storm the MDMW from upstream to downstream. These stations are far water management [39, 40] and regional drought and flood control away from the mouth of the river inland with a distance of 0, 13.9, [38]. This methodology proves to be effective and has few negative 22.9, 29.6, 44.5, and 49.1 km, respectively (Fig. 1 and Tab. 1). effects on surrounding ecosystems. Typically, there are a number of wetlands which are connected by tributaries/river channels in an estuary. These wetlands include marsh, lakes, and reservoirs that 3 Methodology form a network, together with river channels, to regulate hydro- 3.1 Data collection logical processes [38–44]. In a network, the freshwater can be diverted into river channels from wetlands. The increased river In this study, the cross-sectional area and river width which were discharge then decreases the range of saltwater intrusion along defined at the tidal averaged water level were obtained from the river [26–30]. A new equilibrium can be established between measured data at these six fixed gauging stations. The cross-sectional water and salinity in river channels [24]. area and river width of river channels were used to determine the The objectives of this paper are as follows: (i) to investigate the shape of MDMW. The daily mean discharge at Wuzhou station was relationship between the freshwater discharge and the saltwater used to estimate the upstream freshwater discharge. Hourly average intrusion length; (ii) to design a wetland network to mitigate chlorosity of the six fixed gauging stations was used to calculate the saltwater intrusion; and (iii) to compare water allocations between river water salinity at each station. These data were obtained from network method and EWAS and evaluate the performance of wet- the publication [52]. Figure 1 illustrates the location of these land networks in different scenarios. stations.
2 Study area 3.2 Wetland network design The Pearl River is the third longest river in China. PRE is a multi- The framework to design a wetland network is illustrated in Fig. 2. A channel estuary with a drainage density of 0.68–1.07 km/km2 [45–47]. quantitative relationship between freshwater discharge and salinity There is a mixed diurnal and semi-diurnal tide with two troughs was identified, and then, the locations of wetlands were determined and two peaks in a day [48, 49]. The Modaomen waterway (MDMW, based on the distribution of water areas, which included ponds,
Figure 1. (a) Map of China. (b) The PRE and its river network. (c) The MDMW G1, G2, G3, G4, G5, and G6 are gauging stations named as Guadingjiao, Zhupaisha, Xihekou, Zhuyin, Daao, and Baiqing, respectively. S1 is Lianshiwan sluices.
ß 2012 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim www.clean-journal.com 1038 Z. Zhang et al. Clean – Soil, Air, Water 2012, 40 (10), 1036–1046
Table 1. Characteristics of eight segments in the MDMW
Segment Cross-section Distance Length of Qb) (km)a) segment (km) (m3/s)
– X1 50.4 – 1040 X1X2 X2 40.2 10.2 1360 X1X3 X3 40.2 10.2 1360 X2X4 and X3X4 X4 30.6 9.6 1750 X4X5 X5 22.0 8.6 2220 X5X6 X6 16.9 5.1 2600 X6X7 X7 12.7 4.2 3010 X7X8 X8 8.9 3.8 3510 a) Distance away from the cross section Xi (i ¼ 1,2,...,8) to the Figure 3. Variations of the cross-sectional area and the river width of Gudingjiao station (G1). the MDMW with distance away from the estuary mouth (Part of data from b) It is the freshwater discharge at the Wuzhou station. When Ou et al. [57]). the practical freshwater discharge equals the values given in the fifth column of freshwater discharge obtained from Eq. (8), the saltwater intrusion will not reach the corresponding segments.
From Eq. (1), a one-dimensional steady state advection-diffusion equation [27, 58, 59] was used to predict the saltwater intrusion as discarded ditches, and aquaculture water near the river banks. The follows: wetlands were then connected to the existing river channels to form @ðASÞ @ðAvSÞ @ @S a wetland network. Isolated river channels and wetlands were con- þ ¼ AKx (2) nected by new-building river segments. The discharge and the stor- @t @x @x @x age capacity of the wetlands were then evaluated. Finally, if the where S (%) is the salinity at the upstream boundary of the river, t (s) wetland network design was finished, the performance of wetlands is time, v (m/s) is velocity, and Kx is the longitudinal salinity distri- in different scenarios was then evaluated. bution coefficient that is independent of the distance from the mouth. Under steady state [27, 59], by introducing v ¼ Q/A, Eq. (2) 3.2.1 Relationship between freshwater discharge and can be written as: saltwater intrusion length QS dS ¼ K (3) Researches [10, 27, 53–55] found that the cross-sectional area of river A x dx channels exponentially reduced with increasing distance away 3 from the river mouth in some estuaries. This regression model where Q (m /s) is the river discharge from upstream. In this study, the was identified in the MDMW by Chen et al. [56] and Ou et al. [57] freshwater discharge (Q) was collected from the Wuzhou station and was given by: where the hydrological process was not influenced by the tides [49]. The Wuzhou station is approximately 340 km away from the river mouth (see Fig. 1b). Substituting Eq. (1) into Eq. (2), a new A ¼ A0 expð x=aÞ (1) equation can obtain the following [27]:
2 where A (m ) is the cross-sectional area of the river channels and A0 is aQ lnðS=S0Þ¼ ðÞexpðx=aÞ 1 (4) the cross-sectional area at river mouth (G1). x (km) is the length of the KxA0 cross-section away from the river mouth. a (km) is the convergence where S (%) indicates the salinity at the river mouth. The saltwater length which is the length scale of the exponential function against 0 intrusion length of the PRE was set as x ¼ x , where the salinity (S) measured data. It can be seen very clearly in Fig. 3 that convergence L was 0.5%. The threshold of salinity (S ¼ 0.5 %) at the upstream length (a) can be determined by the regression model (Eq. (1)). boundary of the river was extracted from the Standards for Drinking Water Quality in China (GB5749-2006). The saltwater intru- sion length is then determined by: A K lnð0:5=S Þ x ¼ aln 1 0 x 0 (5) L aQ
where xL (km) is the saltwater intrusion length.
In Eq. (5), the negative ln(0.5/S0) indicated that the saltwater
intrusion length (xL) increased with greater cross-sectional area of
the river mouth (A0) and with the higher mixing coefficient (Kx) and decreased with the increase of the river discharge (Q). In this paper, the salinity was converted from the chlorosity using Han’s function, which indicated the salinity was proportional to the chlorosity [60]. The maximum and minimum values were used to calculate the Figure 2. The flow chart of designing a wetland network for mitigation of saltwater intrusion length at high water slack (HWS) and low water saltwater intrusion. slack (LWS), respectively [61]. This model might underestimate the
ß 2012 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim www.clean-journal.com Clean – Soil, Air, Water 2012, 40 (10), 1036–1046 Wetland Network Design for Mitigation of Saltwater Intrusion 1039 range of saltwater intrusion [27, 62, 63]. Thus, a parameter (f) was used to calibrate Eq. (5). The calibrated intrusion length is given as: xL0 ¼ fxL (6) where xL0 (km) is the calibrated saltwater intrusion length. The regression model between upstream freshwater discharge and the calibrated saltwater intrusion length was then determined by: