Next Article in Journal
Decoding Chambal River Shoreline Transformations: A Comprehensive Analysis Using Remote Sensing, GIS, and DSAS
Previous Article in Journal
The Impacts of the Hydrological Regime on the Soil Aggregate Size Distribution and Stability in the Riparian Zone of the Three Gorges Reservoir, China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Radium Isotopes and Hydrochemical Signatures of Surface Water-Groundwater Interaction in the Salt-Wedge Razdolnaya River Estuary (Sea of Japan) in the Ice-Covered Period

Il’ichev Pacific Oceanological Institute, Far Eastern Branch, Russian Academy of Sciences, Vladivostok 690041, Russia
*
Authors to whom correspondence should be addressed.
Water 2023, 15(9), 1792; https://doi.org/10.3390/w15091792
Submission received: 11 March 2023 / Revised: 20 April 2023 / Accepted: 1 May 2023 / Published: 7 May 2023
(This article belongs to the Section Oceans and Coastal Zones)

Abstract

:
The interaction of surface water and groundwater is important in the ecology of coastal basins, affecting hydrological conditions, oxygen regime, carbon, and nutrient exchange. This study demonstrates a dynamic connection between the salt-wedge region and its underlying aquifer in the eutrophic estuary. In winter, this estuary is covered with ice, and the river flow is at its lowest; that is why the specific response to groundwater discharge is best marked in this season. Groundwater admixture was detected in the salt-wedge region by highly active radium isotopes: 223Ra—4.80 ± 0.42 dpm 100 L−1, 224Ra—55.37 ± 1.1 dpm 100 L−1, and 228Ra—189.71 ± 4.66 dpm 100 L−1. The temperature of groundwater and river water was about +4 °C and 0 °C, respectively; that of seawater was −1.6 °C, and temperature increased up to +2.3 °C in the surface water–groundwater interaction region. Groundwater admixture is accompanied by a lower level of oxygen concentration of 52 μmol/kg; at that time, the maximum oxygen concentration in the salt-wedge region was 567 μmol/kg. In waters with a high activity of radium isotopes, there was a maximum partial pressure of CO2—4454 μatm at the range 100–150 μatm in the salt-wedge region and also observed extremum of NH 4 + , NO 2 , and dissolved phosphorus. The surface water–groundwater interaction through anoxic sediment can form localized anaerobic areas despite the general oxygen supersaturation of eutrophic estuary waters and also cause local recycling of nutrients from bottom sediments.

1. Introduction

Surface water–groundwater (SW-GW) interaction significantly influences the biogeochemistry of rivers [1], lakes [2], and coastal salt marshes [3], and the ecology along sea coasts around the world [4]. In coastal basins, SW-GW interaction has an effect on estuary ecosystems [5,6], including the formation of water with low dissolved oxygen (DO) in salt-wedge estuaries [7]. SW-GW interaction is often accompanied by temperature anomalies [8], and thus it can affect the aquatic biota [9,10].
Traditionally known SGD (submarine groundwater discharge), which implies SW-GW interaction in the subterranean estuary [11,12], may enhance benthic flux nutrients and dissolved organic carbon (DOC) [13,14]. Anoxic groundwater inputs account for oxygen deficit in bottom waters, and radium isotopes indicate that hypoxia is caused by offshore groundwater discharge [15,16,17]. The discharge of fresh groundwater in the global water balance is only 4.5 ± 3.2% from the discharge of river water to the ocean [18], but recirculated groundwater contributes 90% of the total water flow SGD [19].
Radium is an alkaline earth element that is adsorbed in river water, but exists as dissolved Ra2+ in seawater due to ion exchanges [20]. Desorption of Ra from river suspension occurs at an early stage of the mixing zone of river and sea waters [21]. Therefore, the dissolved isotopes 224Ra, 223Ra, and 228Ra, with a half-life of 3.6 days, 11.5 days, and 5.7 years, respectively, show a non-conservative dependence on water salinity in estuaries in the presence of river suspension [22]. The second source of dissolved radium for the estuary is SGD. Sea waters influenced by the Ra source in the diurnal and synoptic time ranges are enriched only in short-lived 224Ra and 223Ra isotopes, while the accumulation of long-lived 228Ra isotopes takes years [22]. This makes it possible to establish the time scales of water dynamics and biogeochemical processes in coastal marine areas under the influence of SW-GW interaction and river runoff.
SW-GW interaction has been studied in lagoon estuaries [23,24,25,26] and in wetlands [27]. However, the study of this effect using radioisotopes in ice-covered basins is limited to only a few publications on the Arctic region [28,29]. The ice-covered estuarine ecosystems with a regime known as salt-wedge intrusion [30] have not been studied in terms of the impact of SW-GW interaction. In winter, the estuaries of the Japanese islands are subject to high water conditions [31], while the estuaries of the Russian coast are in low water and freeze-up conditions [32].
The Razdolnaya River Estuary is very eutrophic, and near-bottom hypoxia is observed during both low and high-water periods in summer [33,34]. However, during the freeze-up period, phytoplankton blooms occur with a chlorophyll “a” concentration of more than 100 µg/L, and then the concentration of DO is increased, and the partial pressure of CO2 (pCO2) is decreased [35]. The purpose of this article is to assess the impact produced by SW-GW interaction on the hydrological and hydrochemical conditions of the Razdolnaya River Estuary during the winter low water period covered by ice, based on radionuclides 223Ra, 224Ra, and 228Ra. Thus, the novelty of this study is reflected in the following aspects: we study the SW-GW interaction in the area of the salt-wedge of the estuary with a long period of low water and ice cover; this estuary is highly eutrophic and phytoplankton blooms and superoxia/hypoxia occur near the bottom; we combine geochemical tracers such as three Ra isotopes and classical hydrological profiling performed with discrete synaptic variability during the coldest period of the year and the greatest ice thickness in the area when the groundwater signature is clear.

2. Materials and Methods

2.1. Study Area

The transboundary Razdolnaya River (China–Primorsky Territory of the Russian Federation) flows into the northern part of Amursky Bay (Peter the Great Bay, Sea of Japan) (Figure 1). The estuary measures about 50 km in length and is located within the boggy Razdolnenskaya depression and the northern part of Amursky Bay (Figure 1). The Razdolnaya River catchment area is 16,800 km2. The average discharge of the river for the last 11 years has been 103.5 m3/s (http://gmvo.skniivh.ru/ (accessed on 18 May 2022)). The water regime of Razdolnaya R. is characterized by steady low winter discharge of 7.2 and 5.7 m3/s in January and February, respectively, and an absolute minimum discharge of 1.5 m3/s in February. Discharge increases after the winter low water, usually in early March, but in the harsh winters from mid-March. Spring floods can be observed in May. The peaks of spring floods are about ten times higher than the average annual discharge of the river. The absolute discharge peaks exceed 3000 m3/s during the summer and autumn floods in some years. In the course of freeze-up period from late November to early April, a salt wedge tends to penetrate the Razdolnaya River Estuary at a distance of up to 28–29 km from the river mouth bar [32,35]. During this period, the water salinity is more than 34 PSU in Amursky Bay and up to 26 PSU above the bar [32,35]. The average spring tide in Peter the Great Bay ranges from 15 to 20 cm, which makes it possible to classify the Razdolnaya River Estuary as a micro-tidal estuary with strong water stratification. The tides in the estuary of the Razdolnaya River are observed at a distance more them 20 km from the mouth river bar and were comparable with tides in Amursky Bay [36]. The currents are formed by tides in the estuary of the Razdolnaya River in the winter season. Change in the direction of the currents was observed at the surface and at the bottom layer, depending on the tidal phase. During the flood tide, the current at the upper and bottom water layer at a velocity of up to 20 cm/s was oriented to the river; during the ebb-tide and low water level, direction of the current was opposite at a velocity of up to 20 cm/s [36]. At the same time, the daily variability of the penetration range of the saltwater wedge does not occur [36].

2.2. Field Work, Hydrological Surveys, and Water Sampling

We carry out a series of five surveys at 15 monitoring stations (St.) in Razdolnaya R. estuary (Figure 1) at about a week’s intervals from 22 January 2022 to 23 February 2022 [35]. For each day, we estimated the state of the tide using data from the portal (http://esimo.oceanography.ru/tides/ (accessed on 18 May 2022)). To reduce the possible effect of diurnal variability in relation to tides, all data were taken in periods around high tide times. We probed the water with an RBR maestro multi-channel logger (RBR Ltd., Ottawa, ON, Canada) at an 8 Hz sampling rate. The following properties were logged: pressure, temperature, electrical conductivity, PAR, Chl-a fluorescence, chromophoric OM fluorescence, DO phosphorescence, and turbidity. The conversion from DO concentration to percent saturation was performed using the RBR software. On 22–23 February 2022, water samples were taken from the surface (under the ice) and bottom (0.2 m from the bottom) horizons using a 5-L Niskin bottle at 15 stations. Water samples were taken to measure radionuclides dissolved in water (223Ra, 224Ra, and 228Ra), pH, total alkalinity (TA), nutrients (silica, phosphates, nitrates, nitrites, and ammonium, Ptotal and Ntotal), salinity (S), and DOC.
We also took the water temperature in three water wells located in the estuary valley. Prior to measurements, we pumped out water from these wells to measure characteristics of the inflowing groundwater. The value of 4 °C was typical of groundwater temperatures.
We made an experiment on the desorption of radium isotopes from river suspension by adding salt to a sample of river water at St. 1 to salinity of 21.5 psu. The resulting sample was filtered one day after sampling and included in the total series of measurements of radium isotopes.
Four sediment cores were taken with a gravity tube at stations 2, 4, 5, and 10 for visual description and estimates of the activity of radium isotopes in sediments [37].

2.3. Radium Isotopes Measurements

Water samples with a volume of 44 L were filtered through 1 mm polypropylene cartridges and Mn fiber at a flow rate of maximum 1 L min−1 to obtain a Ra extraction efficiency of at least 97% [38]. Radioisotopes 224Ra and 223Ra were measured using a Radium Delayed Coincidence Counter (RaDeCC, West Columbia, SC, USA) [39]. To make correction for the supported 224Ra, we conducted the second set of measurements in 2 to 6 weeks so that the initial excessive activity of 224Ra (ex 224Ra) could reach secular equilibrium with 228Th, which was also absorbed on Mn fiber [39]. The activity of 228Ra was measured after 5–6 months, with correction for the decay of 228Th, which had initially been sorbed onto Mn fiber from original samples [38]. To calibrate RaDeCC systems for 224Ra, 223Ra, and 228Ra, measurements using 232Th standards as described by [40] were used.

2.4. Analysis of Hydrochemical Characteristics

The NH 4 + concentration was determined using the indophenols method. NO 3 , NO 2 , dissolved silicates (DSi), and dissolved inorganic phosphorus (DIP) were measured using standard colorimetric methods. Optical density was determined on a KFK-3KM photocolorimeter (Ruspribor, Saint Petersburg, Russia). Details of the methods used for the nutrient analyses are given in Grasshoff et al. [41]. The sum of the NH 4 + , NO 3 , and NO 2 concentrations is the dissolved inorganic nitrogen (DIN). The detection limit was 0.01 μmol/L for the phosphate and DIN, and 0.02 μmol/L for DSi. The total phosphorus (Ptotal) and total nitrogen (Ntotal) were determined by persulfate combustion with ultraviolet irradiation on a SKALAR SAN++ analyzer (Skalar, Breda, The Netherlands). The detection limit was 0.01 μmol/L for the Ptotal and Ntotal.
Water samples for determining DOC were prefiltered through a glass filter with pore size of 0.6 μm. The analysis used automated analyzer Shimadzu TOC-VCPN (Shimadzu, Kyoto, Japan). Salinity was also measured by salinometer Guildline Autosal 8400B (Guildline Instruments, Smiths Falls, ON, Canada) with accuracy of 0.002.
The pCO2 was calculated from the measured pH and TA using a commonly known procedure [42]. A potentiometric method was applied to determine pH on InoLab pH/ION/Cond 750 m (Laborkomplekt, Moscow, Russia). pH was measured at 10 °C using a cell without a liquid junction [43]. The precision of pH measurements was about ±0.004 pH units. TA analysis was carried out through direct colorimetric titration on Dosimat-665 (Metrohm, Herisau, Switzerland) with hydrochloric acid in an open cell according to Bruevich’s method [43,44]. TA measurements were performed with a precision of ±3 µmol/kg.
The software used for statistical analyses was MS Excel 2019. The spatial distribution maps were developed using the program Surfer 9 (Golden Software).

3. Results

3.1. Hydrological Conditions

The observations were made during a steady state salt-wedge penetration to about 29 km from the river mouth bar (Figure 2) [35]. Further, salt-wedge penetration into the estuary was limited to the sandy channel between St. 1 and St. 2, with a depth of less than 0.5 m. As our observations showed, it was frozen to the bottom. There was a general decrease in the salinity of the bottom layer, which was most pronounced for St. 10, where we observed a decrease from 30 to 25 PSU (Figure 2). Freshening of the bottom layer was accompanied by a gradient decrease in the halocline layer, which can be distinguished by the salinity jump from 8 to 25 PSU (Figure 2).
The temperature of sea and river waters was about −1.6 °C and 0 °C, respectively (Figure 3). The near-bottom layer temperature increased up to +2.3 °C at St. 2 during the observation period at the salt-wedge region [35].
The ice thickness on the transect as a whole increased during the observation period (Table 1) [35]. However, in the area of the bar on St. 12, the thickness of the ice remained almost unchanged or even decreased in early February. In this case, the water exchange rate between the estuary channel and the receiving basin decreases. The river discharge during the ice formation period is minimized, and the estuary, in this case, can be compared to a closed lake system to a certain extent.

3.2. Radium Isotopes Activity and Distribution of Hydrochemical Characteristics

As seen in Table 2, the activity of short-lived isotopes 224Ra and 223Ra increased in the sample of river water to which we added salt; however, the activity of the long-lived isotope 228Ra was unchanged. The highest total activity rate of radium isotopes 224Ra, 223Ra, and 228Ra is observed at St. 2, and an increased activity rate is observed all across the topographic depression with its salt wedge at St. 2–10. Therefore, the highest activity of all three isotopes of radium on St. 2 is mainly associated with the influence of SW-GW interaction and, secondarily, with ion-exchange processes with an increase in salinity. The names of water types are presented based on the activity of isotopes of radium in endmember water (Table 2). Ratios close to equilibrium, 224Ra/223Ra = 21.7 and 224Ra/228Ra = 0.67, were obtained for St. 4, which indicates the admixture of recent groundwater in this region. The highest activity of 228Ra was found at St. 2, with 224Ra/228Ra = 0.2. That is, in this place, 228Ra accumulated without the addition of 224Ra. Generally, we obtained the irregular distribution of the three isotopes and their ratios, probably associated with the irregular inflow of groundwater in different places and the transport of these waters by currents up and down the channel depending on the tides (Table 2).
A minimum oxygen saturation of about 21%, maximum pCO2 of 4500 µatm, and maximum radium isotope activity exist in the estuary’s near-bottom waters (Figure 4). It was a local anomaly in the St. 2 area. At the same time, water is oversaturated with oxygen, from 110 to 150%, in seawater and at stations in the region of the salt wedge (Figure 4) due to the phytoplankton bloom outbreak [35].
Figure 5 shows the non-conservative relationship of nutrients to salinity. Ntotal and DIN markers at St. 2 and in several samples from the surface layer are located above the mixing line, which indicates their addition. The remaining markers are predominantly located below the mixing line, indicating the intake of Ntotal in the mixing zone of river and sea waters. For Ptotal and DIP, St. 2 is located below the mixing line and other markers.
The maximum concentrations of NH 4 + (73.7 µmol/L) and NO 2 (2.1 µmol/L) were observed in waters with the highest activity of 228Ra at St. 2 (Figure 6). Other nutrients and DOC also increase in St. 2 (Figure 6).

4. Discussion

4.1. Identification and Control Factorsof SW-GW Interaction

The combination of high radium isotope activity rate (Table 2) and temperature anomaly (Figure 3) shows that the water of the salt wedge is formed as a result of the admixture of groundwater.
However, we accept the following assumptions. It has been reported that saline waters interact with sediments, which are a source of radium isotopes, through various processes such as erosion, diffusion, and bio irrigation for the coastal embayment [25]. Other studies confirm that the contribution of 228Ra due to erosion remains the main one for the shelf [45], including our results obtained for the summer flood in the Razdolnaya River Estuary [46]. The role of SW-GW interaction in the export of short-lived radium isotopes in the Razdolnaya River Estuary increases in some of the deepest areas, for example, the 224Ra maximum at St. 4 (Table 2). However, the contribution of pore exchange to the export of 224Ra for ichthyofauna exists in the Razdolnaya River Estuary [46]. The listed factors suggest careful use of 223Ra, 224Ra, and 228Ra as groundwater tracers if there are no additional evident indicators, such as water temperature in our case.
It is known that relatively short-period seawater recirculation in the upper aquifer is accompanied by enrichment with short-lived radium isotopes 224Ra and 223Ra. Admixing salt water in this aquifer will bring the long-lived isotope 228Ra, but long-term processes are essential for the significant accumulation of 228Ra [12,25,38,47]. The composition of stable isotopes δ18O and δD in the SW-GW interaction zone is subject to the seawater/river water ratio [48]. It is considered that the main reason for the increase in all three isotopes 223Ra, 224Ra, and 228Ra is that water recirculated into the upper aquifer during the winter period. The upper, highly permeable aquifer comprises pebbles and sand, and its thickness is 20 m [49]. This is the alluvial–marine complex that extends 30 km from the Razdolnaya River mouth bar, and seawaters penetrate into it in winter along the estuary. The low-water period and salt-wedge regime begin from September–October [33], and then the seawater penetrates for many kilometers from the bar, as seen in Figure 2. Salinization of the upper aquifer in the estuary area for 5 months from September–October to February leads to the accumulation of dissolved long-lived 288Ra isotopes in it. The subsequent recirculation of these waters will lead to an observed increase in all three isotopes with an increase in water temperature in winter (Figure 3). Salinity in the bottom layer of estuary water decreases during the low water period, since there is inflow in the estuary by fresh groundwater from the land side, but there is no inflow of cold salt water (Figure 2 and Figure 3).
As a result of SW-GW interaction, we observed an increase in all three isotopes 223Ra, 224Ra, and 228Ra but in different intensities at various stations. We believe that the SW-GW interaction in the Estuary of the Razdolnaya River is in the areas with sand sediments, similar to other permeable sand systems [50]. Exactly at St. 2 and 4, we found sand and silty sand in the core, in contrast to other areas of the Razdolnaya River Estuary, where gray silt can be found. In addition, the sandbars along the estuary located between our stations throughout the salt wedge (Figure 2) are also completely composed of sand. We believe that SW-GW interaction with the underlying aquifer occurs in the local deepest section lines of the estuary and probably also on sandbars between our stations.
The circulation of groundwater has a compound system in each local region, while the recharge of the upper aquifer is usually determined by the intensity of rainfall everywhere [51]. The favorable conditions for salt-wedge water recirculation in the upper aquifer exist during the winter discharge low period [52] because the coastal aquifer along Razdolnaya River Estuary is the least loaded with meteoric water during this period [49]. It is known that tides generate complex groundwater fluctuations in aquifers [52]. In turn, semidiurnal variations in recirculated groundwater flux also have a period of tidal sea level fluctuations [12,47,53]. In the salt-wedge region of the Razdolnaya River Estuary, semidiurnal tides are observed [36]. If the salt groundwater meets with fresh land waters in the upper aquifer, density gradients in the transitional zone cause convective circulation, and groundwater discharge occurs [54,55].
Thus, the main drivers for SW-GW interaction increasing in the salt-wedge region of Razdolnaya River Estuary during the freeze-up period are as follows: increased proportion of salt water in the upper aquifer along of estuary since autumn, slow water dynamics as a result of ice-covered and low river discharge, increased duration of the hydraulic head during the period of high tide and low tide, and the high hydraulic conductivity of the underlying aquifer.

4.2. River Water, Seawater, and Groundwater Fractions in the Bottom Water Layer

The river water, seawater, and groundwater contributions can be quantified by applying a mass balance calculation [56] based on 228Ra and salinity data. It is assumed that each sample is a mixture of marine water (fSW), river runoff (fRW), and SG (fSG). The balance is governed by the following equations:
fSW + fRW + fSG = 1
fSW × SSW + fRW × SRW + fSG × SSG = Smeas
fSW × RaSW + fRW × RaRW + fSG × RaSG = Rameas
where fSW, fRW, and fSG are the fractions of marine water, river runoff, and SG in a water parcel, and SSW, SRW, SSG, RaSW, RaRW, and RaSG are the corresponding salinity and 228Ra activity. Smeas and Rameas are the measured values of salinity and 228Ra of the water samples. A special selection of endmember values is required for the Razdolnaya River Estuary in the ice-covered period (see Table 2). The 228Ra activity of sediments is taken 270 dpm 100 L−1 in St. 2. We used this isotope also because the increase in its activity was not associated with suspended matter, as shown by the experiment on adding salt to river water (see Table 2).
Figure 7 shows the admixture of groundwater in the salt-wedge region, the fractions of which increase from the sea with decreased salinity (Figure 2). The fraction of groundwater prevails in St. 2 (Figure 7), where the water temperature was also the highest (Figure 3).

4.3. Potential Estuarine Ecosystem Response to SW-GW Interaction

We suggest that the groundwater seepage through the anoxic sediments leads to a decrease in the DO of groundwater, similar to processes shown and discussed in other papers [4,15,16,17,37]. As a result, in the course of SW-GW interaction between the upper aquifer and salt-wedge water, we see DO saturation rates in the near-bottom layer reducing toward the estuary top, and a simultaneous increase in water temperatures and radium isotopes activity (Table 2, Figure 3 and Figure 4). In the warm season, the salt wedge in the Razdolnaya River Estuary has a DO close to zero due to the intensification of production-destruction processes as well as to the high turbidity of water and, accordingly, photosynthetically active radiation close to zero of the near-bottom water layer [33]. However, suspended matter concentrations are relatively low in winter, and the photic zone normally extends to the bottom everywhere, with nutrient concentrations being very high [32]. That is why photosynthesis prevails over organic matter destruction, and DO oversaturation is observed in the Razdolnaya River Estuary in winter [35]. Not excepting that the key factor of DO-specific conditions for the Razdolnaya R. estuary is the balance of organic matter production–destruction, results (Figure 4) are indicative of SW-GW interaction influence on the formation of DO-specific and high pCO2 conditions of the salt-wedge region.
This occurs when DO is depleted in sediments, and the distraction of organic matter occurs in the following reactions: denitrification, reduction of irons (Fe) and manganese oxides, reduction of sulfates, or methanogenesis [57]. Estuary sediments Razdolnaya River Estuary are enriched with iron as a result of their flux with river water [58]. When the salt water penetration into sediments, the destruction of organic matter occurs with the reduction of sulfates and the accumulation of NH 4 + [59]. The destruction of organic matter in the sediment leads to an increase in the concentration of DIP in the pore water.
In the winter period, the dissolved nitrogen and phosphorus concentrations of the Razdolnaya R. are up to 500 and 9 µmol/L, respectively, because the flux of nutrients from municipal wastewater exists [60]. Eutrophication of river waters is accompanied by the accumulation of nutrients and organic matter in the sediments of the estuary. Therefore, the concentration of NH 4 + in the pore water of sediments reaches 2088 µmol/L, and the pCO2 value is 11,739 µatm [61]. DIP in sediments is usually sorbed, but the redox potential strongly influences the dissolution of DIP in sediments [62]. In sediments from eutrophic basins, the redox boundary can be located within a few millimeters of the upper sediment layer, and DIP is released into the pore water below this boundary [63]. Thus, Fe compounds can be a source or sink of DIP in bottom sediments depending on the oxygen content there. This can explain the local DIP extremum with the extremums of NH 4 + and NO 2 in the SW-GW interaction region at St. 2 (Figure 6 and Figure 7).

5. Conclusions

The SW-GW interaction is estimated with a geotracer-based model, including a series of hydrological profilings conducted to estimate the synaptic variability of a complex of hydrological and hydrochemical characteristics.
This study highlights the fundamental role of the SW-GW interaction in the formation of the temperature and hydrochemical regime of an ice-covered estuary when the estuary water exchange is significantly limited, and mixing due to wind is excluded. SW-GW interaction through anoxic sediment can form localized anaerobic areas despite the general oxygen supersaturation of eutrophic estuary waters, and the following important conclusions were drawn:
  • SW-GW interaction has a significant effect on the ecological state of the estuary because it is accompanied by a response in the concentration of DO, CO2, and nutrients, thus forming local recycling from sediments;
  • The SW-GW interaction can be an additional source of organic matter from sediments;
  • SW-GW interaction increases the temperature of the salt-wedge region in the ice-covered estuary.
We also obtained various ratios of 224Ra/223Ra and 224Ra/228Ra throughout the estuary. High activity of 224Ra and 223Ra in some cases was associated with desorption from suspension. The high activity of isotope 228Ra, together with temperature anomalies, most clearly indicated the admixture of groundwater in the salt wedge area. The results indicate a variety of SW-GW interactions during the freezing period, probably due to the tidal dynamics of waters and depth variability throughout the estuary.

Author Contributions

P.S.: writing—original draft, conceptualization, methodology, formal analysis, investigation, and funding acquisition. P.T. (Pavel Tishchenko): writing—review and editing, investigation, conceptualization, supervision, funding acquisition. A.C.: writing—review and editing, investigation, formal analysis, funding acquisition, and methodology. G.P.: writing—review and editing, investigation, and formal analysis. Y.B.: methodology, formal analysis, and investigation. A.L.: methodology, formal analysis, and investigation. P.T. (Petr Tishchenko): writing—review and editing, and investigation. E.S.: methodology, formal analysis, and investigation. M.S.: methodology, formal analysis, and investigation. All authors have read and agreed to the published version of the manuscript.

Funding

This work was carried out with the financial support of the Russian Science Foundation (grant no. 21-77-00028 and 19-17-00058) and programs of the POI FEB RAS (121021500052-9, AAAA-A20-120011090005-7).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Any datasets not included in this published work are available upon reasonable request from the corresponding author; however, all data produced or analyzed during this investigation are contained in it.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ledford, S.H.; Briggs, M.; Glas, R.; Margaret, A.; Zimmer, N.A. Connecting diverse disciplines to improve understanding of surface water-groundwater interactions. J. Hydrol. X 2022, 17, 100141. [Google Scholar] [CrossRef]
  2. Mackie, C.; Levison, J.; Binns, A.; O’Halloran, I. Groundwater-surface water interactions and agricultural nutrient transport in a Great Lakes clay plain system. J. Great Lakes Res. 2021, 47, 145–159. [Google Scholar] [CrossRef]
  3. Xin, P.; Wilson, A.; Shen, C.; Ge, Z.; Moffett, K.B.; Santos, I.R.; Chen, X.; Xu, X.; Yau, Y.Y.Y.; Moore, W.; et al. Surface water and groundwater interactions in salt marshes and their impact on plant ecology and coastal biogeochemistry. Rev. Geophys. 2022, 60, e2021RG000740. [Google Scholar] [CrossRef]
  4. Moore, W.S.; Joye, S.B. Saltwater intrusion and submarine groundwater discharge: Acceleration of biogeochemical reactions in changing coastal aquifers. Front. Earth Sci. 2021, 9, 600710. [Google Scholar] [CrossRef]
  5. Rocha, C.; Robinson, C.E.; Santos, I.R.; Waska, H.; Michael, H.A.; Bokuniewicz, H.J. A place for subterranean estuaries in the coastal zone. Front. Earth Sci. 2021, 250, 107167. [Google Scholar] [CrossRef]
  6. Garcia-Orellana, J.; Rodellas, V.; Tamborski, J.; Diego-Feliu, M.; van Beek, P.; Weinstein, Y.; Charette, M.; Alorda-Kleinglass, A.; Michael, H.A.; Stieglitz, T.; et al. Radium isotopes as submarine groundwater discharge (SGD) tracers: Review and recommendations. Earth-Sci. Rev. 2021, 14, 103681. [Google Scholar] [CrossRef]
  7. Santos, I.R.; Cook, P.L.M.; Rogers, L.; De Weys, J.; Eyre, B.D. The “salt wedge pump”: Convection-driven pore-water exchange as a source of dissolved organic and inorganic carbon and nitrogen to an estuary. Limnol. Oceanogr. 2012, 57, 1415–1426. [Google Scholar] [CrossRef]
  8. Anderson, M.P. Heat as a ground water tracer. Ground Water 2005, 43, 951–968. [Google Scholar] [CrossRef]
  9. Morin, A.; Lamoureux, W.; Busnarda, J. Empirical Models Predicting Primary Productivity from Chlorophyll a and Water Temperature for Stream Periphyton and Lake and Ocean Phytoplankton. J. N. Am. Benthol. Soc. 1999, 18, 299–307. [Google Scholar] [CrossRef]
  10. Shuter, B.J.; Finstad, A.G.; Helland, I.P.; Zweimüller, I.; Hölker, F. The role of winter phenology in shaping the ecology of freshwater fish and their sensitivities to climate change. Aquat. Sci. 2012, 74, 637–657. [Google Scholar] [CrossRef]
  11. Moore, W.S. The subterranean estuary: A reaction zone of ground water and sea water. Mar. Chem. 1999, 65, 111–125. [Google Scholar] [CrossRef]
  12. Burnett, W.C.; Aggarwal, P.K.; Aureli, A.; Bokuniewicz, H.; Cable, J.E.; Charette, M.A.; Kontar, E.; Krupa, S.; Kulkarni, K.M.; Loveless, A.; et al. Quantifying submarine groundwater discharge in the coastal zone via multiple methods. Sci. Total Environ. 2006, 367, 498–543. [Google Scholar] [CrossRef] [PubMed]
  13. Kim, T.-H.; Waska, H.; Kwon, E.; Suryaputra, I.G.N.; Kim, G. Production, degradation, and flux of dissolved organic matter in the subterranean estuary of a large tidal flat. Mar. Chem. 2012, 142–144, 1–10. [Google Scholar] [CrossRef]
  14. Santos, I.R.; Bryan, K.B.; Pilditch, C.A.; Tait, D.R. Influence of porewater exchange on nutrient dynamics in two New Zealand estuarine intertidal flats. Mar. Chem. 2014, 167, 57–70. [Google Scholar] [CrossRef]
  15. McCoy, C.; Viso, R.; Peterson, R.N.; Libes, S.; Lewis, B.; Ledoux, J.; Voulgaris, G.; Smith, E.; Sanger, D. Radon as an indicator of limited cross-shelf mixing of submarine groundwater discharge along an open ocean beach in the South Atlantic Bight during observed hypoxia. Cont. Shelf Res. 2011, 31, 1306–1317. [Google Scholar] [CrossRef]
  16. Peterson, R.N.; Moore, W.S.; Chappel, S.L.; Viso, F.R.; Libes, S.M.; Peterson, L.E. A new perspective on coastal hypoxia: The role of saline groundwater. Mar. Chem. 2016, 179, 1–11. [Google Scholar] [CrossRef]
  17. Guo, X.; Xu, B.; Burnett, W.C.; Wei, Q.; Nan, H.; Zhao, S.h.; Charette, M.A.; Lian, E.; Chen, G.; Yu, Z. Does submarine groundwater discharge contribute to summer hypoxia in the Changjiang (Yangtze) River Estuary? Sci. Total Environ. 2020, 719, 137450. [Google Scholar] [CrossRef]
  18. Abbott, B.W.; Bishop, K.; Zarnetske, J.P.; Minaudo, C.; Chapin III, F.S.; Krause, S.; Hannah, D.M.; Conner, L.; Ellison, D.; Godsey, S.E.; et al. Human domination of the global water cycle absent from depictions and perceptions. Nat. Geosci. 2019, 12, 533–540. [Google Scholar] [CrossRef]
  19. Kwon, E.Y.; Kim, G.; Primeau, F.; Moore, W.S.; Cho, H.-M.; DeVries, T.; Sarmiento, J.L.; Charette, M.A.; Cho, Y.-K. Global estimate of submarine groundwater discharge based on an observationally constrained radium isotope model. Geophys. Res. Lett. 2014, 41, 8438–8444. [Google Scholar] [CrossRef]
  20. Webster, I.T.; Hancock, G.J.; Murray, A.S. Modelling the effect of salinity on radium desorption from sediments. Geochim. Cosmochim. Acta 1995, 59, 2469–2476. [Google Scholar] [CrossRef]
  21. Hancock, G.J.; Murray, A.S. Source and distribution of dissolved radium in the Bega River estuary, Southeastern Australia. Earth Planet. Sci. Lett. 1996, 138, 145–155. [Google Scholar] [CrossRef]
  22. Moore, W.S. Using the Ra quartet for evaluating groundwater input and water exchange in salt marshes. Geochim. Cosmochim. Acta 1996, 60, 4645–4652. [Google Scholar] [CrossRef]
  23. Young, M.B.; Gonneea, M.E.; Fong, D.A.; Moore, W.S.; Herrera-Silveira, J.; Paytan, A. Characterizing sources of groundwater to a tropical coastal lagoon in a karstic area using radium isotopes and water chemistry. Mar. Chem. 2008, 109, 377–394. [Google Scholar] [CrossRef]
  24. Wang, X.; Du, J. Submarine groundwater discharge into typical tropical lagoons: A case study in eastern Hainan Island, China. Geochem. Geophys. Geosystems 2016, 17, 4366–4382. [Google Scholar] [CrossRef]
  25. Rodellas, V.; Garcia-Orellana, J.; Trezzi, G.; Masqué, P.; Stieglitz, T.C.; Bokuniewicz, H.; Cochran, J.K.; Berdalet, E. Using the radium quartet to quantify submarine groundwater discharge and porewater exchange. Geochim. Cosmochim. Acta 2017, 196, 58–73. [Google Scholar] [CrossRef]
  26. Tamborski, J.; van Beek, P.; Rodellas, V.; Monnin, C.; Bergsma, E.; Stieglitz, T.; Heilbrun, C.; Cochran, J.K.; Charbonnier, C.; Anschutz, P.; et al. Temporal variability of lagoon–sea water exchange and seawater circulation through a Mediterranean barrier beach. Limnol. Oceanogr. 2019, 64, 2059–2080. [Google Scholar] [CrossRef]
  27. Carol, E.; Pilar Alvarez, M.D.; Santucci, L.; Candanedo, I.; Arcia, M. Origin and dynamics of surface water-groundwater flows that sustain the Matusagaratí Wetland, Panamá. Aquat. Sci. 2022, 84, 16. [Google Scholar] [CrossRef]
  28. Lecher, A. Groundwater discharge in the Arctic: A review of studies and implications for biogeochemistry. Hydrology 2017, 441, 41. [Google Scholar] [CrossRef]
  29. Dabrowski, J.S.; Charette, M.A.; Mann, P.J.; Ludwig, S.M.; Natali, S.M.; Holmes, R.M.; Schade, J.D.; Powell, M.; Henderson, P.B. Using radon to quantify groundwater discharge and methane fluxes to a shallow, tundra lake on the Yukon-Kuskokwim Delta, Alaska. Biogeochemistry 2020, 148, 69–89. [Google Scholar] [CrossRef]
  30. Valle-Levinson, A. Contemporary Issues in Estuarine Physics; Cambridge University Press: Cambridge, UK, 2010; p. 315. [Google Scholar] [CrossRef]
  31. Funahashi, T.; Kasai, A.; Ueno, M.; Yamashita, Y. Effects of short time variation in the river discharge on the salt wedge intrusion in the Yura estuary, a micro tidal estuary, Japan. J. Water Resour. Prot. 2013, 5, 343–348. [Google Scholar] [CrossRef]
  32. Zvalinsky, V.I.; Maryash, A.A.; Shvetsova, M.G.; Sagalaev, S.G.; Tishchenko, P.Y.; Stonik, I.V.; Begun, A.A. Production and hydrochemical characteristics of ice, under-ice water and sediments in the Razdolnaya River estuary (Amursky Bay, Sea of Japan). Russ. J. Mar. Biol. 2010, 36, 270–281. [Google Scholar] [CrossRef]
  33. Tishchenko, P.Y.; Semkin, P.J.; Tishchenko, P.P.; Zvalinsky, V.I.; Barabanshchikov, Y.u.A.; Mikhailik, T.A.; Sagalaev, S.G.; Shvetsova, M.G.; Shkirnikova, E.V.; Shulkin, V.M. Hypoxia of bottom waters of the Razdolnaya River Estuary. In Doklady Earth Sciences; Pleiades Publishing: New York, NY, USA, 2017; Volume 476, pp. 1207–1211. [Google Scholar] [CrossRef]
  34. Tishchenko, P.Y.; Tishchenko, P.P.; Lobanov, V.B.; Mikhaylik, T.A.; Sergeev, A.F.; Semkin, P.Y.; Shvetsova, M.G. Impact of the transboundary Razdolnaya and Tumannaya Rivers on deoxygenation of the Peter the Great Bay (Sea of Japan). Estuar. Coast. Shelf Sci. 2020, 239, 106731. [Google Scholar] [CrossRef]
  35. Semkin, P.; Tishchenko, P.; Pavlova, G.; Barabanshchikov, Y.; Tishchenko, P.; Shvetsova, M.; Shkirnikova, E.; Fedorets, Y. O2 and CO2 Responses of the synaptic period to under-ice phytoplankton bloom in the eutrophic Razdolnaya River Estuary of Amur Bay, the Sea of Japan. J. Mar. Sci. Eng. 2022, 10, 1798. [Google Scholar] [CrossRef]
  36. Semkin, P.Y.; Tishchenko, P.Y.; Charkin, A.N.; Pavlova, G.Y.; Anisimova, E.V.; Barabanshchikov, Y.A.; Mikhailik, T.A.; Tibenko, E.Y.; Tishchenko, P.P.; Shvetsova, M.G.; et al. Isotopic, hydrological and hydrochemical indicators of submarine groundwater discharge in the estuary of Razdolnaya River (Amursky Bay, Sea of Japan) in the ice covered period. Geoecol. Eng. Geol. Hydrogeol. Geocryol. 2021, 3, 29–43. (In Russian) [Google Scholar] [CrossRef]
  37. Tamborski, J.J.; Cochran, J.K.; Bokuniewicz, H.J. Submarine groundwater discharge driven nitrogen fluxes to Long Island Sound, NY: Terrestrial vs. marine sources. Geochim. Cosmochim. Acta. 2017, 218, 40–57. [Google Scholar] [CrossRef]
  38. Moore, W.S. Fifteen years experience in measuring 224Ra and 223Ra by delayed coincidence counting. Mar. Chem. 2008, 109, 188–197. [Google Scholar] [CrossRef]
  39. Moore, W.S.; Arnold, R. Measurement of 223Ra and 224Ra in coastal waters using delayed coincidence counter. J. Geophys. Res. Ocean. 1996, 101, 1321–1329. [Google Scholar] [CrossRef]
  40. Moore, W.S.; Cai, P. Calibration of RaDeCC systems for 223Ra measurements. Mar. Chem. 2013, 156, 130–137. [Google Scholar] [CrossRef]
  41. Grasshoff, K.; Ehrhard, M.; Kremling, K. Methods of Seawater Analysis; Verlag Chemie: Weinheim, Germany, 1983; p. 419. Available online: https://scholar.google.com/scholar_lookup?title=Methods+of+Seawater+Analysis&author=Grasshoff,+K.&author=Ehrhard,+M.&author=Kremling,+K.&publication_year=1983 (accessed on 18 May 2022).
  42. Dickson, A.G.; Sabine, C.L.; Christian, J.R. (Eds.) Guide to Best Practices for Ocean CO2 Measurements; PICES Special Publication 3; PICES: Sidney, BC, Canada, 2007; 191p, Available online: http://hdl.handle.net/11329/249 (accessed on 31 October 2022).
  43. Tishchenko, P.; Zhang, J.; Pavlova, G.; Tishchenko, P.; Sagalaev, S.; Shvetsova, M. Revisiting the Carbonate Chemistry of the Sea of Japan (East Sea): From Water Column to Sediment. J. Mar. Sci. Eng. 2022, 10, 438. [Google Scholar] [CrossRef]
  44. Pavlova, G.Y.; Tishchenko, P.Y.; Volkova, T.I.; Dickson, A.; Wallmann, K. Intercalibration of Bruevich’s method to determine the total alkalinity in seawater. Oceanology 2008, 48, 438–443. [Google Scholar] [CrossRef]
  45. van der Loeff, M.R.; Kipp, L.; Charette, M.A.; Moore, W.S.; Black, E.; Stimac, I.; Charkin, A.; Bauch, D.; Valk, O.; Karcher, M.; et al. Radium isotopes across the Arctic Ocean show time scales of water mass ventilation and increasing shelf inputs. J. Geophys. Res. Oceans 2018, 123, 4853–4873. [Google Scholar] [CrossRef]
  46. Semkin, P.Y.; Tishchenko, P.Y.; Tishchenko, P.P.; Pavlova, G.Y.; Anisimova, E.V.; Barabanshchikov, Y.A.; Mikhailik, T.A.; Charakin, A.N.; Shvetsova, M.G. Radium Isotopes and Nutrients in Razdolnaya River Estuary (Amur Bay, the Sea of Japan) in the Period of Summer Flood. Water Res. 2022, 49, 429–439. [Google Scholar] [CrossRef]
  47. Taniguchi, M.; Burnett, W.C.; Cable, J.E.; Turner, J.V. Investigation of submarine groundwater discharge. Hydrol. Process. 2002, 16, 2115–2129. [Google Scholar] [CrossRef]
  48. Semkin, P.Y.; Tishchenko, P.Y.; Charkin, A.N.; Pavlova, G.Y.; Tishchenko, P.P.; Anisimova, E.V.; Barabanshchikov, Y.A.; Leusov, A.E.; Mikhailik, T.A.; Tibenko, E.Y.; et al. Discharge of salt groundwater in the Estuary of Razdolnaya River (Amursky Bay) in February 2020. Water Res. 2021, 48, 254–259. [Google Scholar] [CrossRef]
  49. Chelnokov, G.A.; Kharitonova, N.A.; Zykin, N.N.; Vereshchagina, O.F. Genesis of the mineral groundwaters of the Razdol’nenskii occurrence in Primorye. Russ. J. Pac. Geol. 2008, 2, 521–527. [Google Scholar] [CrossRef]
  50. Santos, I.R.; Burnett, W.C.; Dittmar, T.; Suryaputra, I.G.N.A.; Chanton, J. Tidal pumping drives nutrient and dissolved organic matter dynamics in a Gulf of Mexico subterranean estuary. Geochim. Cosmochim. Acta 2009, 73, 1325–1339. [Google Scholar] [CrossRef]
  51. Zhao, Y.; Zhang, Y.-K.; Yang, Y.; Li, F.; Xiao, S. Groundwater Circulation in the Xianshui River Fault Region: A Hydrogeochemical Study. Water 2020, 12, 3310. [Google Scholar] [CrossRef]
  52. Vallejos, A.; Sola, F.; Pulido-Bosch, A. Processes Influencing Groundwater Level and the Freshwater-Saltwater Interface in a Coastal Aquifer. Water Res. Manag. 2015, 29, 679–697. [Google Scholar] [CrossRef]
  53. Smith, A.J.; Turner, J.V. Density-dependent surface water-groundwater interaction and nutrient discharge in the Swan-Canning Estuary. Hydrol. Process. 2001, 15, 2595–2616. [Google Scholar] [CrossRef]
  54. Kobayashi, S.; Sugimoto, R.; Honda, H.; Miyata, Y.; Tahara, D.; Tominaga, O.; Shoji, J.; Yamada, M.; Nakada, S.; Taniguchi, M. High-resolution mapping and time-series measurements of Rn-222 concentrations and biogeochemical properties related to submarine groundwater discharge along the coast of Obama Bay, a semi-enclosed sea in Japan. Prog. Earth Planet. Sci. 2017, 4, 6. [Google Scholar] [CrossRef]
  55. Eissa, M.A.; Thomas, J.M.; Pohll, G.; Shouakar-Stash, O.; Hershey, R.L.; Dawoud, M. Groundwater recharge and salinization in the arid coastal plain aquifer of the Wadi Watir delta, Sinai, Egypt. Appl. Geochem. 2016, 71, 48–62. [Google Scholar] [CrossRef]
  56. Luo, X.; Jiao, J.J. Submarine groundwater discharge and nutrient loadings in Tolo Harbor, Hong Kong using multiple geotracer-based models, and their implications of red tide outbreaks. Water Res. 2016, 102, 11–31. [Google Scholar] [CrossRef]
  57. Regnier, P.; Dale, A.W.; Arndt, S.; LaRowe, D.E.; Mogollón, J.; Van Cappellen, P. Quantitative analysis of anaerobic oxidation of methane (AOM) in marine sediments: A modeling perspective. Earth-Sci. Rev. 2011, 106, 105–130. [Google Scholar] [CrossRef]
  58. Shulkin, V.M.; Tishchenko, P.Y.; Semkin, P.Y.; Shvetsova, M.G. Influence of river discharge and phytoplankton on the distribution of nutrients and trace metals in Razdolnaya River estuary, Russia. Estuar. Coast. Shelf Sci. 2018, 211, 166–176. [Google Scholar] [CrossRef]
  59. Wang, Y.; Jiao, J.J. Origin of groundwater salinity and hydrogeochemical processes in the confined Quaternary aquifer of the Pearl River Delta, China. J. Hydrol. 2012, 438, 112–124. [Google Scholar] [CrossRef]
  60. Mikhailik, T.A.; Tishchenko, P.Y.; Koltunov, A.M.; Tishchenko, P.P.; Shvetsova, M.G. The effect of Razdol’naya River on the environmental state of Amur Bay (the Sea of Japan). Water Res. 2011, 38, 512–521. [Google Scholar] [CrossRef]
  61. Pavlova, G.Y.; Tishchenko, P.Y.; Khodorenko, N.D.; Shvetsova, M.G.; Sagalaev, S.G. Major ion composition and carbonate equilibrium in the sediment pore water of the Razdol’naya River estuary of Amur Bay, the Sea of Japan. Russ. J. Pac. Geol. 2012, 6, 251–262. [Google Scholar] [CrossRef]
  62. Ruttenberg, K.C.; Sulak, D.J. Sorption and desorption of dissolved organic phosphorus onto iron (oxyhydr)oxides in seawater. Geochim. Cosmochim. Acta 2011, 75, 4095–4112. [Google Scholar] [CrossRef]
  63. Mort, H.P.; Slomp, C.P.; Gustafsson, B.G.; Andersen, T.J. Phosphorus recycling and burial in Baltic Sea sediments with contrasting redox conditions. Geochim. Cosmochim. Acta 2010, 74, 1350–1362. [Google Scholar] [CrossRef]
Figure 1. Map of study area: (a) Sea of Japan, (b) The Peter the Great Bay, and (c) position of monitoring stations at Razdolnaya R. estuary.
Figure 1. Map of study area: (a) Sea of Japan, (b) The Peter the Great Bay, and (c) position of monitoring stations at Razdolnaya R. estuary.
Water 15 01792 g001
Figure 2. Salinity data on transect obtained from 22 January 2022 to 23 February 2022. The dates shown in the figure indicate the day of profiling. The position of monitoring stations is shown by dotted red lines in first graph. Station site locations are indicated in Figure 1. Positive values on the X-axis—distance from the river mouth bar (St. 12) to upstream; negative values—to downstream.
Figure 2. Salinity data on transect obtained from 22 January 2022 to 23 February 2022. The dates shown in the figure indicate the day of profiling. The position of monitoring stations is shown by dotted red lines in first graph. Station site locations are indicated in Figure 1. Positive values on the X-axis—distance from the river mouth bar (St. 12) to upstream; negative values—to downstream.
Water 15 01792 g002
Figure 3. Water temperature on transect obtained from 22 January 2022 to 23 February 2022. The dates shown in the figure indicate the day of profiling. The position of monitoring stations is shown by dotted red lines in first graph. Station site locations are indicated in Figure 1. Positive values on the X-axis—distance from the river mouth bar (St. 12) to upstream; negative values—to downstream.
Figure 3. Water temperature on transect obtained from 22 January 2022 to 23 February 2022. The dates shown in the figure indicate the day of profiling. The position of monitoring stations is shown by dotted red lines in first graph. Station site locations are indicated in Figure 1. Positive values on the X-axis—distance from the river mouth bar (St. 12) to upstream; negative values—to downstream.
Water 15 01792 g003
Figure 4. Dependence between DO (% saturation) and pCO2 on the activity of 223Ra, 224Ra, and 228Ra. 23 February 2022.
Figure 4. Dependence between DO (% saturation) and pCO2 on the activity of 223Ra, 224Ra, and 228Ra. 23 February 2022.
Water 15 01792 g004
Figure 5. Dependence of nutrient concentrations on salinity. Blue lines connect markers for river and sea waters. 23 February 2022.
Figure 5. Dependence of nutrient concentrations on salinity. Blue lines connect markers for river and sea waters. 23 February 2022.
Water 15 01792 g005
Figure 6. Dependence between nutrients and DOC on the activity of 228Ra. 23 February 2022.
Figure 6. Dependence between nutrients and DOC on the activity of 228Ra. 23 February 2022.
Water 15 01792 g006
Figure 7. Distribution of water fractions in the bottom layer (RW—river water, SW—sea water, GW—groundwater). Station site locations are indicated in Figure 1. 23 February 2022.
Figure 7. Distribution of water fractions in the bottom layer (RW—river water, SW—sea water, GW—groundwater). Station site locations are indicated in Figure 1. 23 February 2022.
Water 15 01792 g007
Table 1. Ice thickness (cm).
Table 1. Ice thickness (cm).
Station Numbers
Dates123456789101112131415
22 January7570656065645656605650405348
28 January797366646565585661585142565145
4 February837568656767565762575035505549
11 February817670697070565562585040646054
23 February907966676969555561595241706660
Table 2. Depth (m) of sampling, salinity (S) (PSU), and activities radium isotopes (dpm 100 L−1). RW: river water, RW + Salt: river water sample with added salt, SWR: salt-wedge region, SW: seawater, and HTBW: high-turbidity brackish water (mixture of RW and SW). Date of sampling: 23 February 2022.
Table 2. Depth (m) of sampling, salinity (S) (PSU), and activities radium isotopes (dpm 100 L−1). RW: river water, RW + Salt: river water sample with added salt, SWR: salt-wedge region, SW: seawater, and HTBW: high-turbidity brackish water (mixture of RW and SW). Date of sampling: 23 February 2022.
St. №/Water MassDepth of SamplingSex223Raex224Ra224Ra/223Ra228Ra224Ra/228Ra
1/RW8.30.140.05 ± 0.0097.89 ± 0.19157.811.95 ± 0.820.66
1/RW + Salt8.321.54.38 ± 0.3837.27 ± 1.048.57.84 ± 0.564.76
2/SWR515.234.80 ± 0.4238.74 ± 1.098.1189.71 ± 4.660.2
4/SWR 6.821.092.55 ± 0.2355.37 ± 1.121.782.29 ± 2.10.67
5/SWR 6.324.191.36 ± 0.1827.01 ± 0.8219.979.62 ± 3.260.34
6/SWR7.824.560.70 ± 0.1130.00 ± 0.8742.949.81 ± 1.480.6
7/SWR7.324.990.79 ± 0.1322.31 ± 0.8728.281.95 ± 4.770.27
8/SWR6.725.030.33 ± 0.0720.76 ± 0.8462.961.57 ± 3.320.34
9/SWR6.825.130.79 ± 0.0615.12 ± 0.2419.168.70 ± 3.290.22
10/SWR422.881.39 ± 0.0619.64 ± 0.2914.155.44 ± 1.970.35
11/HTBW1.319.250.64 ± 0.0416.57 ± 1.3925.95.92 ± 0.422.8
12/HTBW0.519.10.41 ± 0.0310.91 ± 0.6726.621.11 ± 1.500.52
13/SW3.932.520.15 ± 0.0152.37 ± 1.6315.811.58 ± 0.920.2
14/SW8340.13 ± 0.0042.13 ± 0.0316.420.44 ± 2.120.1
15/SW131.150.08 ± 0.0151.72 ± 0.08721.523.40 ± 1.630.07
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Semkin, P.; Tishchenko, P.; Charkin, A.; Pavlova, G.; Barabanshchikov, Y.; Leusov, A.; Tishchenko, P.; Shkirnikova, E.; Shvetsova, M. Radium Isotopes and Hydrochemical Signatures of Surface Water-Groundwater Interaction in the Salt-Wedge Razdolnaya River Estuary (Sea of Japan) in the Ice-Covered Period. Water 2023, 15, 1792. https://doi.org/10.3390/w15091792

AMA Style

Semkin P, Tishchenko P, Charkin A, Pavlova G, Barabanshchikov Y, Leusov A, Tishchenko P, Shkirnikova E, Shvetsova M. Radium Isotopes and Hydrochemical Signatures of Surface Water-Groundwater Interaction in the Salt-Wedge Razdolnaya River Estuary (Sea of Japan) in the Ice-Covered Period. Water. 2023; 15(9):1792. https://doi.org/10.3390/w15091792

Chicago/Turabian Style

Semkin, Pavel, Pavel Tishchenko, Alexander Charkin, Galina Pavlova, Yuri Barabanshchikov, Andrey Leusov, Petr Tishchenko, Elena Shkirnikova, and Maria Shvetsova. 2023. "Radium Isotopes and Hydrochemical Signatures of Surface Water-Groundwater Interaction in the Salt-Wedge Razdolnaya River Estuary (Sea of Japan) in the Ice-Covered Period" Water 15, no. 9: 1792. https://doi.org/10.3390/w15091792

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