Submarine fresh groundwater discharge into Laizhou Bay comparable to the Yellow River flux.

Near- and off-shore fresh groundwater resources become increasingly important with the social and economic development in coastal areas. Although large scale (hundreds of km) submarine groundwater discharge (SGD) to the ocean has been shown to be of the same magnitude order as river discharge, submarine fresh groundwater discharge (SFGD) with magnitude comparable to large river discharge is never reported. Here, we proposed a method coupling mass-balance models of water, salt and radium isotopes based on field data of (223)Ra, (226)Ra and salinity to estimate the SFGD, SGD. By applying the method in Laizhou Bay (a water area of ~6000 km(2)), we showed that the SFGD and SGD are 0.57 ~ 0.88 times and 7.35 ~ 8.57 times the annual Yellow River flux in August 2012, respectively. The estimate of SFGD ranges from 4.12 × 10(7) m(3)/d to 6.36 × 10(7) m(3)/d, while SGD ranges from 5.32 × 10(8) m(3)/d to 6.20 × 10(8) m(3)/d. The proportion of the Yellow River input into Laizhou Bay was less than 14% of the total in August 2012. Our method can be used to estimate SFGD in various coastal waters.

Submarine groundwater discharge (SGD) is any and all flows of water on continental margins from the seabed to the coastal ocean, regardless of fluid composition or driving force1. SGD, which is driven by both terrestrial and marine forcing components, comprises terrestrial fresh groundwater (SFGD) and re-circulated seawater (RSGD)12. As an essential part of SGD, SFGD is an important source of freshwater, nutrients, contaminants, and other chemicals to the coastal waters, and has significant impacts and implications on coastal environment and ecology34. With near half of the global population residing within 100 km of the coastline5, it is important to find more freshwater resources including near- and off-shore fresh groundwater for the relief of water scarcity in densely populated coastal megacities6. Many previous studies have found that regional scale SGD is several times greater than river discharge789. For example, Kim et al.7 estimated the magnitude of SGD into the Yellow Sea to be at most 300% of the river water input. Moore et al.8 showed that the SGD flux is probably between 0.8 and 1.6 times the river flux to the Atlantic. Moore9 reevaluated the SGD to a large section of the South Atlantic Bight and found the annual average SGD flux is three times greater than the river fluxes. However, there are few studies on large scale SFGD as compared with those on SGD, and SFGD with magnitude comparable to large river discharge is never reported. In this paper, we proposed a method for combining mass-balance models of water, salt and radium isotopes to estimate the SFGD, SGD. By applying the method in Laizhou Bay, the SFGD, SGD, and the proportion of the Yellow River input into Laizhou Bay were estimated and discussed.

Laizhou Bay, located between 37.05°N ~ 37.80°N and 118.9°E ~ 120.35°E (WGS84 reference system), is one of the three major bays in the Bohai Sea, China (Fig. 1). The natural coastline of the Bay is relatively straight and extends from the Qimu Cape to the Yellow River Estuary. Along the coastline, there are at least four different types of depositional environments: the Yellow River delta in the northwest, the alluvial plain in the southwest, the marine deposit plain in the south, and the hilly area in the east101112. The coast areas west of Hutouya (Fig. 1) are alluvial or marine deposit plains with aquifers mainly composed of permeable coarse material. The coast area east of Hutouya is hilly with a coastal plain. Groundwater occurs mainly in the Quaternary aquifers of the coastal plains. The average annual rainfall is ~640 mm and occurs mostly from June to September, while the potential evaporation is approximately 2050 mm13. The seawater depth of the Bay generally increases as the offshore distance increases, and has an average depth of ~8 m. Based on salinity measurements at surface, middle and bottom at the stations S2–S6 shown in Fig. 1, where the seawater is deepest in whole Laizhou Bay, it was found that the Bay is vertically well mixed and can be treated as a single layer in terms of the large-scale water motion.

Figure 1

Laizhou Bay and sampling stations.

The dots, triangles and pentagrams represent sampling stations for seawater (S), groundwater (G) and river water (R), respectively. Station numbers as defined in Supplementary Table S1 are marked next to each of the stations. The numbers 1, 2, and 3 in the inset indicate Laizhou Bay, Bohai Bay and Liaodong Bay, respectively. Maps were created with MAPGIS 6.7 software.

There are many rivers flowing into the Bay, including the Yellow River, Jiaolai River, Xiaoqing River and Wei River. Particularly, the Yellow River, as the second and sixth largest river in China and the world, respectively, is the largest that discharges into the Bohai Sea. However, the Yellow River discharge has significantly decreased since the 1950s due to both climate change and human activities1415. Seasonal variations of the Yellow River discharge are also significant, with the minimum occurring usually in April and May, and the maximum occurring usually in July and August (Fig. 2). The annual average discharge of the Yellow River in 2012 was 7.23 × 107 m3/d. In July and August 2012, its discharge rate was at least 50 times greater than the sum of the discharge rates of all other rivers flowing into Laizhou Bay.

Figure 2

Monthly average water discharge of the Yellow River to the Bohai Sea from 2012 to 2014.

The dashed black, blue and red lines represent annual average discharge from 1956 to 2000, in 2012, and 2013, respectively. Data from the Yellow River Water Resources Bulletin of Yellow River Conservancy Commission of the Ministry of Water Resources.

The Yellow River deposits large amounts of sediments, creating a fast-growing delta area. Its course into the Bohai Sea changed a number of times over the past several decades16. The present channel through the delta was formed artificially in 1996. The mouth of the Yellow River is located between Bohai Bay and Laizhou Bay, and the shifting of the mouth to the north may significantly reduce its direct discharge into Laizhou Bay.

As one of the three bays in the Bohai Sea, China, and with a shoreline of ~320 km and area of ~6000 km2, Laizhou Bay is an important coastal environment. As well, it is subject to a variety of environmental stresses, and hence provides an archetype of a semi-enclosed bay for which ecological functioning is a sensitive issue. Most existing studies on SGD in the Bohai Sea are restricted to the Yellow River delta1718192021, the largest estuary in the Bohai Sea. In order to estimate SFGD, SGD, and the Yellow River input into Laizhou Bay, activities of radium isotopes (223, 226Ra) and salinities were measured (Fig. 1). Natural radium isotopes are ideal tracers for effective and efficient assessment of SGD since they are conservative chemically and widely enriched in groundwater relative to surface waters782122. In August 2012, we collected eight groundwater samples, six river water samples, and 44 seawater samples for measurement of radium isotopes (223, 226Ra) and salinity (Supplementary Table S1). In May 2014, we measured the salinities at the same locations.

Results

Spatial distribution of radium isotopes and salinity

Figure 3 shows the spatial distributions of the two radium isotopes (223, 226Ra) and salinity within Laizhou Bay. The activity distributions of these two isotopes have the following common features: (1) the activities were significantly higher in the west and south than in the east and north of the Bay; and (2) the activities were very high in the estuary and near-shore areas and they generally decreased with the offshore distance (Figs. 3a and 3b). The seawater salinities in August 2012 were higher in the east than in the west of the Bay (Fig. 3c), where the activities of the radium isotopes were comparatively high. The seawater salinities in May 2014 were similar to those in August 2012 overall, but with a local low-salinity area near the middle of the eastern coast (Fig. 3d), indicating considerable SFGD since there are no river inputs in that region.

Figure 3

Contour plots.

(a), 223Ra in August 2012. (b), 226Ra in August 2012. (c), Salinity in August 2012. (d), Salinity in May 2014. Activities of radium isotopes and salinity were measured in surface water (1–2 m below the surface) of Laizhou Bay from 19 to 26 August 2012. Salinity was also measured at the same locations from 1 to 6 May 2014. Maps were created with Surfer 8.0 software.

Determination of flushing time

To explore the dynamics of coastal water, one should first estimate the flushing time T [T] for a bay, i.e., the ratio of the mass or volume of a constituent (V) to its renewal rate (Q)23, or T = V/Q. Combining the radium isotope method24 and the tidal prism method25, we obtained the following equation to estimate the flushing time T based on measurements of radium isotopes,where F(223Ra/226Ra) is the223Ra/226Ra activity ratio of the flux into the bay, I(223Ra/226Ra) is the 223Ra/226Ra activity ratio in the bay, and λ223 = 0.061d−1 is the decay constant of 223Ra. Using equation (1), we estimated T = 36.6 ± 5.3 d, which is in line with previous independent estimates172026 (see Flushing Time Model in Supplementary Information).

Assessment of SFGD and SGD

To quantify the freshwater fluxes into Laizhou Bay, we developed water-mass and salt-mass balance models that include seawater from outside the Bay, river input, SFGD, precipitation (P) and evaporation (E), all estimated over the period of T d immediately before the observation time (taken as the mid-point of the 8-d field sampling period).

Coupling the salt-mass balance model with the water-mass balance model, we derived the SFGD flux (Q) as (see Water and Salt Mass Balance Model in Supplementary Information):where r is the ratio of the Yellow River input into Laizhou Bay to the total input into the Bohai Sea (Q) during the flushing time immediately before the observation time; Q is runoff of the ith river other than the Yellow River into the Bay (Supplementary Table S2); M is the total salt mass in the Bay; and S is the salinity of the representative seawater outside of the Bay. Detailed calculations of these parameter values are described and summarized in Supplementary Table S3.

In order to quantify the fluxes of SGD into Laizhou Bay, we developed a 226Ra mass balance model within the whole bay. The inputs of 226Ra are from discharge waters (river, SGD, and seawater input from the open Bohai Sea outside of the Bay), desorption from suspended particles and diffusion from bottom sediments. Neglecting decay, gain from precipitation and loss from evaporation, the loss of 226Ra only includes mixing with the open sea. Using the steady state premise, one can ignore the variation of 226Ra storage in the Bay. In this case the SGD flux is given by (see 226Ra Mass Balance Model in Supplementary Information)where 226Rais the total input of 226Ra from the Yellow River discharge and suspended particle desorption; 226Rais the total input of 226Ra from discharges and suspended particle desorption of all rivers other than the Yellow River; 226Ra and 226Ra are the radium activity in groundwater and Bay water, respectively; and 226Ra is the background activity. Detailed calculations of these parameter values are described and summarized in Supplementary Table S4.

Figure 4a shows how the SFGD and SGD fluxes change with the proportion of the Yellow River input into Laizhou Bay (r). The annual average flux QY2012 = 7.23 × 107 m3/d of the Yellow River in 2012 was used as the reference, which is more than 20 times the fluxes of all the other rivers flowing into Laizhou Bay. One can see that both SFGD and SGD decrease as r increases. Note that SGD is much greater than SFGD since SGD includes RSGD. From Fig. 4a, since SFGD ≥ 0, we have SGD ≥ 5.13QY2012 and r ≤ 0.396. When r = 0, we obtain the upper-bound estimate of SGD as 8.57QY2012 and the corresponding value of SFGD as 0.88QY2012. In general, we have SGD = (8.57 − 8.69r)QY2012, and SFGD = (0.88 − 2.22r)QY2012.

Figure 4

Changes of SFGD and SGD with the proportion of Yellow River input into Laizhou Bay.

The SFGD and SGD are nondimensionalized by QY2012, the annual average flow rate of the Yellow River in 2012, which equals 7.23 × 107 m3/d. (a), SFGD and SGD as predicted by equations (2) and (3) using data from August 2012. The ratio of SFGD to SGD is shown on the right-axis. (b), SFGD as predicted by equation (2) using salinity data in May 2014. The blue bands bounded by the pair of dotted-lines represent the error ranges corresponding to the range of T = 36.6 ± 5.3 d, with longer flushing time corresponding to smaller SFGD or SGD. Note that the parameter r varies with time.

Seasonal variations of the Yellow River flux were significant. The minimum flux occurred in April–May and was less than one-tenth of the maximum flux in July–August (Fig. 2) during 2012 and 2013. Thus, SFGD estimation based on equation (2) using the flux of the Yellow River during April–May can effectively reduce the impact of the Yellow River. Considering that tides are the dominant forcing for the renewal of the seawater in the Bay27, it is reasonable to assume that the seasonal variation of the flushing time in Laizhou Bay is negligible. So the flushing time of T = 36.6 ± 5.3 d, which was obtained with radium data in August 2012, may also be used for April–May 2014. Figure 4b shows how the SFGD predicted by equation (2) using the salinity data observed in May 2014 changes with r, the proportion of the Yellow River input into Laizhou Bay, during April and May in 2014. As r increases from 0 to 1, the SFGD decreases from 0.79QY2012 to 0.34QY2012. As expected, the impact of the Yellow River flux on the SFGD was significantly reduced by comparison with that of August 2012 (Fig. 4a).

The precipitation in Laizhou Bay is much larger in June–August than in March–May, and so the groundwater table in the coastal area should be higher (i.e., higher hydraulic head) in August than in May owing to rainfall infiltration into ground surface. Due to the increased landward gradient in hydraulic head, it follows that SFGD (August 2012) ≥ SFGD (May 2014). Using the lower-bound estimate for SFGD in May 2014 (0.34QY2012), from Fig. 4a one can see that if SFGD in August 2012 ≥ 0.34QY2012, then r (August 2012) was less than 0.24 and SGD (August 2012) ≥6.48QY2012 (as is indicated by the thick red-dashed vertical line in Fig. 4a).

Although we do not have direct evidence to conclude that r (May 2014) ≤0.24 since r varies with time, it should be considerably less than unity. In addition, the Yellow River mouth, which is located between Bohai Bay and Laizhou Bay, is northerly-oriented, as shown in photographs taken by NASA Landsat in May 2014 or Fig. 2 of Xu et al.21. This fact may significantly limit the input of the Yellow River into Laizhou Bay. Thus, it is at least not unreasonable to assume r (May 2014) ≤0.50. From Fig. 4b one can see that r (May 2014) ≤0.50 implies SFGD (May 2014) ≥0.57QY2012. Since SFGD (August 2012) ≥SFGD (May 2014), from Fig. 4a one can see that SFGD in August 2012 was in the range (0.57 ~ 0.88)QY2012, and r (August 2012) was less than 0.14 and SGD in August 2012 was between 7.35QY2012 (as is indicated by the thin red-dashed vertical line in Fig. 4a) and 8.57QY2012.

The seawater circulation in the Bohai Sea shows apparent seasonal variations, which in turn affect the path of the Yellow River discharge28. During the summer months, the monsoonal winds blow from the south in this region, thus creating a cyclonic gyre within the Bohai Sea2628. These may limit the input of the Yellow River into Laizhou Bay and support the small estimated value of r (August 2012).

Discussion

Our tracers-orientated estimation of SGD is approximately ten times as large as SFGD. SGD is much greater than SFGD since SGD = SFGD + RSGD, and many previous studies show that tides, waves, fluid density gradient, storms, geothermal gradient, and seabed topography can result in RSGD2930313233343536. As a result, the SGD values estimated by the tracers-orientated method and mechanism-orientated method such as traditional hydrogeological approach often do not match. Li and Jiao29 reviewed the studies of tidal contribution to SGD and found that the RSGD induced by tides in the intertidal zone is at least 1 ~ 2 orders of magnitude smaller than SGD estimated by radium isotope tracers. The order of magnitude of RSGD contributed by waves is the same as that of tides3031. The density-driven RSGD is usually much less than SFGD323334. The RSGD induced by geothermal gradient and seabed topography is at least 4 orders of magnitude smaller than the SGD estimated by radium isotope tracers353637. Thus, SGD estimated by radium isotope tracers is much larger than that of the total sum of SGD induced by all the above factors/mechanisms. Solving this challenging problem needs long-term efforts involving close combination of tracers-orientated and mechanism-orientated methods.

We can compare our results with small-scale SGD studies. Taniguchi et al.19 used seepage meters to estimate SGD fluxes of 2300 m3/m/d in September 2004, 3065 m3/m/d in September 2006 for a 7-km offshore area approximately 40 km south of the Yellow River estuary. They also estimated the average SFGD of 18 and 28 m3/m/d in September 2004 and 2006, respectively. Based on their results, the average ratio of SFGD to SGD is less than 1% in both sampling periods. Our results, if using the same unit (treating SGD as a line source along the coastline and dividing the total SGD by the shoreline length of the Bay), ranges from 1662.2 m3/m/d to 1937.2 m3/m/d, which are slightly smaller than Taniguchi's estimates. The ratio of SFGD to SGD in our study, ranging from 7.5% to 9% as shown in Fig. 4a when r ≤ 0.14, however, is much higher than that of their study. Peterson et al.18 estimated SGD in the same area as Taniguchi et al.19 using radon and radium isotopes. They estimated a SGD flux of 4.5–13.9 cm/d in September 2006, most of which was recirculated seawater. Our results, if using the same unit (treating SGD as an area source and dividing the total SGD by the area of the Bay), ranges from 8.9 cm/d to 10.3 cm/d, which are within the range of Peterson's estimates for the same season. Since our tracers-orientated estimations of SGD and SFGD are based on the mass balance of radium isotopes and fresh groundwater after they entered the seawater, they are conservative in the sense that if there had no groundwater exploitations in the coastal area of Laizhou Bay1238, their values would be larger.

The Bohai Sea has three major bays: Laizhou Bay, Bohai Bay and Liaodong Bay (Fig. 1). The seawater salinity in Laizhou Bay was significantly lower than that in other two bays39. The real reason for this was not investigated in detail because it seems consistent with the freshwater discharge from the Yellow River abutting Laizhou Bay. Our investigation gives a plausible explanation for the abnormally low salinity in Laizhou Bay. We concluded that the SFGD, which accounts for 57% ~ 88% of the annual flux of the Yellow River in 2012, is a key contributor to the low salinity given that the direct input from the Yellow River to Laizhou Bay is very limited (r ≤ 14% in August 2012).

SGD has been widely recognized to be a pathway for enriching coastal waters in nutrients, carbon and metals1822. In some areas, nutrient fluxes via SGD were shown to rival those from surface waters74041. The nutrient input via SGD in the Yellow River Estuary is at least five times of that via the Yellow River20. With new understanding of our assessments of the SFGD, SGD and the Yellow River input into Laizhou Bay, the management of the Bay related to fresh groundwater resources, ecology and environment in coastal and offshore areas should be reviewed.

Although SGD has been estimated in many coastal areas all over the world based on radium isotope tracer methods2942, large scale tracers-orientated SFGD studies have not been correspondingly conducted even if SFGD is as important as, or even more important than SGD. In this study, we proposed a tracers-orientated method (using radium isotopes and salinity in seawater and coastal groundwater as tracers) to estimate SFGD. To the knowledge of the authors, this is the first time to quantify large-scale SFGD using a tracers-orientated method based on field radium and salinity measurements. Our method has potential for application considering that the seawater salinity is a common physical quantity that can be measured easily and a straightforward tracer to indicate freshwater. The proposed method can be readily applied to estimate SFGD in other coastal areas all over the world.

The main limitation of current study is most probably the steady state premise, a common approach used in all the previous tracers-orientated SGD studies by radium isotope methods78222443. Since the variation of radium storage in the whole bay approaches zero near the time when the radium mass in the whole bay reaches its maximum, and the observation period of our field work is just near such a maximum-time, the steady state assumption is approximately valid. Strict quantification of the error induced by steady state assumption, however, needs not only much more radium isotope measurement data in terms of time series, but also quantification of the seawater flow in the whole bay for a long period. This will be an interesting and challenging work in the future.

Methods

Sampling

Radium samples were collected from and adjacent to Laizhou Bay. Large volume water samples (~60 L for seawater, ~15 L for river water, respectively) were pumped and filtered through a 0.45-μm filter for Ra extraction. The coastal groundwater samples (~15 L) were taken from the nearshore zone within 100 m of the high tide mark with PushPoint samplers inserted into sediments at a depth of ~1.5 m. These samples can capture the chemical components (particularly radium isotopes) of groundwater immediately before it discharges into the sea from the aquifer. The water samples were passed slowly through Mn-fibers (~25 g) produced according to the method proposed by Moore44 for extracting radium isotopes. The flow rate was controlled not to exceed 1 L/min to ensure complete Ra adsorption on the Mn-fiber. These fibers were then washed thoroughly to remove all particles and salts and taken to the hydrogeological laboratory at The University of Hong Kong for measurements. The salinity, temperature, pH of the water samples were measured in situ using a HI9828 Model probe (HANNA).

Measurements

The long-lived radium isotope, 226Ra, was determined by a radon-in-air monitor (RAD7, Durridge Co.) as proposed by Kim et al.45. After the 223Ra measurements were completed, the Mn-fiber samples were aged for 2 - 6 weeks to allow 222Rn and its daughters to equilibrate with 226Ra. This is an indirect measurement of 226Ra based on secular equilibrium between 226Ra and 222Rn. In order to improve the test results, the determination of 226Ra in our study slightly modified the method of Kim et al.45. Before measurement we made long fiber incubation time (> 20 d) to ensure secular equilibrium between 226Ra and 222Rn and then allowed ~24 h for purging the RAD7 system prior to analysis thereby reducing the background noise. In addition, we reduced errors by using long measurement times46. The expected error of 226Ra measurements is ±7%. The short-lived radium isotope 223Ra was analyzed using a two-channel radium delayed coincidence counting system (RaDeCC)47. The expected error of 223Ra measurements is ±12%.

Models

Three models (flushing time, water and salt mass balance models) as well as the 226Ra mass balance model were used in the study to estimate the flushing time of coastal water, the freshwater discharge, the Yellow River input, and SGD, respectively. Details were shown in Supplementary Information.