4.1. Major ion chemistry
Major ion concentrations in the surface water of Devi River estuary are given in Table 1. Sodium is the most dominant cation having a concentration of 3593 ± 3105 mg/l during the studied period (Fig. 2; Table 1). Concentration of Mg ranges from 26 to 1006 mg/l, 24 to 967 mg/l and 30 to 992 mg/l during summer, monsoon and winter season respectively. Potassium is the least abundant cation, and varies between 34 and 497 mg/l during summer, 3 and 107 mg/l in monsoon, and 16 and 360 mg/l in the winter period. The dissolved Ca concentration in the estuary ranges from 12 to 584 mg/l. Chloride is the most abundant anion in the estuarine water with a concentration of 7107 ± 6504 mg/l. The PO4 concentration is found be the lowest, and varies between 0.30 and 0.89 mg/l. Thus, the major ion concentrations in the surface water are in the order of: Cl > Na > SO4 > Mg > Ca > K > NO3, DSi > PO4. The observed concentrations are well within the range obtained in other global estuaries (Cohen et al. 1999; Patra et al. 2012; Zhang et al. 2020).
Table 1
Concentration of major ions in the surface water of Devi River estuary.
|
Summer
|
Monsoon
|
Winter
|
Na
|
5927.65 ± 3063.17
(379.00–9495.00)
|
1056.34 ± 1235.65
(15.50–4063.00)
|
3796.08 ± 2561.71
(23.10–6478.00)
|
K
|
298.12 ± 159.63
(34.00–497.00)
|
37.33 ± 31.91
(3.00–107.00)
|
220.41 ± 136.69
(16.00–360.00)
|
Ca
|
316.78 ± 186.23
(32.40–584.00)
|
156.59 ± 108.87
(12.00–293.00)
|
250.98 ± 131.67
(21.80–412.70)
|
Mg
|
352.55 ± 288.39
(26.10–1006.00)
|
296.16 ± 359.86
(23.50–967.00)
|
543.79 ± 344.31
(30.40–991.70)
|
Cl
|
10964.88 ± 5747.30
(619.00–18123.00)
|
2075.14 ± 2367.06
(73.00–7240.00)
|
8281.06 ± 5635.36
(58.70–14124.00)
|
SO4
|
1829.96 ± 953.76
(74.27–2774.00)
|
358.55 ± 344.50
(49.61–1090.81)
|
949.63 ± 599.28
(39.52–1669.1)
|
NO3
|
9.61 ± 4.78
(4.71–17.38)
|
13.93 ± 4.60
(4.03–19.43)
|
9.78 ± 5.46
(4.13–17.98)
|
PO4
|
0.67 ± 0.21
(0.33–0.89)
|
0.33 ± 0.01
(0.32–0.36)
|
0.32 ± 0.01
(0.30–0.34)
|
SiO2
|
11.81 ± 5.12
(6.29–20.32)
|
18.31 ± 4.55
(9.10–21.64)
|
4.16 ± 1.10
(3.04–5.96)
|
4.2. Spatial and seasonal distribution
All major ion concentration, except NO3 and DSi, are found to increase towards the mouth of the estuary, and show significant positive correlations with salinity in all seasons (Fig. 3–5; Table S1). Similar behavior of these ions has been reported in other estuaries (Ramanathan et al. 1993; Kumar et al. 2015; Wang et al. 2009). Concentration of NO3 and DSi decreases from upper to lower estuary. This is due to the fluvial inputs of NO3 and DSi to the estuary. Domestic sewage, manure and agricultural fertilizers contribute largely to the NO3 concentration in rivers (Chen et al. 2021). However, primary producers and microbial denitrification rapidly utilize this NO3 in the ocean water (Singh and Ramesh 2011; Bartoli et al. 2021). This causes lower NO3 concentration towards the mouth of the estuary. The DSi in rivers and coastal waters are mainly derived by weathering and erosion of continental crust (Tréguer and De La Rocha 2013). Biological uptake of DSi by diatoms removes a significant part of this riverine DSi which results in a decrease in its concentration towards the lower estuary (Tréguer et al. 1995). Thus, a strong negative correlation between DSi and salinity is observed. Similar decrease in concentration of NO3 and DSi with increasing salinity have been observed in various other studies (Padmavathi and Satyanarayana 1999; Fatema et al. 2015; Sanders and Laanbroek 2018).
Seasonally, significantly higher concentration of Na, K, Ca, Cl, SO4 and PO4 have been observed during summer period in comparison to winter while the lowest concentrations are found in monsoon season (Fig. 2). Low water discharge during the dry seasons (summer and winter) leads to greater residence time. This causes the dissolved ionic species to be retained within the estuary (John et al. 2020). Also, higher temperatures in summer results in the increased rates of evaporation over precipitation (Montagna et al. 2018). Thus, the combined influence of lower freshwater input and higher evaporation rates have led to peak concentration of these ions in surface water during summer. On the other hand, increased precipitation and high freshwater discharge during monsoon season has a diluting effect on these major ions and causes a reduction in ionic concentration during monsoon (Sumner and Belaineh 2005). Contrastingly, the concentrations of NO3 and DSi are found to be relatively lower in the dry periods in comparison to the wet period. Greater particle supply via increased erosion rates mostly occur during monsoon period (Mangalaa et al. 2017). This increases the greater availability of silicate minerals for dissolution, which releases DSi into the water column (Asano et al. 2003; Sospedra et al. 2018). Significantly, higher terrestrial runoff during monsoon have led to increased nitrate concentrations (Hershey et al. 2021). Such pattern in seasonal variation of major ion chemistry has been observed in other estuaries (Sarin et al. 1985; Jordan et al. 1991; Padmavathi and Satyanarayana 1999).
4.3. Saturation indices of minerals
The saturation index (SI) is a theoretical indicator for understanding the dissolution and/or precipitation potential of a mineral in an aqueous system. It is estimated using the ionic activity potential (IAP) and solubility product (Ksp) following Eq. 1. When SI < 0, the solution is considered to be under-saturated with respect to the mineral, and leading to its dissolution. However, when SI > 0, the solution is oversaturated with respect to the particular mineral, and its precipitation is expected. An aqueous solution with SI = 0 will promote neither dissolution or precipitation as the solution is considered to be in equilibrium.
The saturation index of calcite (SICal) and aragonite (SIArag) ranges from − 1.73 to 0.98 and − 1.87 to 0.83 respectively, and having the highest SI during summer and lowest in winter in the Devi estuary (Akhtar et al. 2021). The SIDol varies between − 2.84 and 2.68 and is also found to be have the highest value in summer and lowest during monsoon. Spatially, increase in the SICal, SIArag and SIDol are observed towards the mouth of estuary during summer and winter (Fig. 6). Similar, distribution of SI along the salinity gradient has been reported in several other estuaries (Jiang et al. 2010; Wit et al. 2018; Li et al. 2020). Significant positive correlation between pH and SI values is found in the estuary (Table S1). Dependency of SI on pH of water is well reported by various authors (Abril et al. 2003; Li et al. 2020). An increase in pH towards the mouth of estuary is evident from its positive correlation with salinity (Table S1). River water is generally characterized by low pH due to the production of carbonic acid via organic matter decomposition (Wang et al. 2011; Hossain and Marshall 2014). Thus, lowering of pH in the upper estuary may have resulted in a decrease in SI values. Additionally, strong positive correlation between dissolved oxygen saturation (DO%) and SI of calcite (r = 0.69 to 0.83), aragonite (r = 0.69 to 0.83) and dolomite (r = 0.65 to 0.86) are observed in the study area (Table S1). The influence of heterotrophic processes such as organic matter degradation on DO% has been previously reported in the Devi estuary (Akhtar et al. 2021). During decomposition of organic matter, a decrease in DO% is accompanied by a reduction in pH which further leads to the lowering of SI. Hence, a coupling between carbonate saturation state in the estuarine water and the process of organic matter degradation is inferred. Similar coupled behavior has been reported in several other estuaries (Abril et al. 2003; Guo et al. 2021). Increase in pH towards the mouth of the estuary is due to the buffering effect of seawater, which causes higher SI in the water. Both SI and pH show unsystematic spatial variation in the monsoon season.
4.4. Mixing behavior
Distribution of dissolved ions in estuaries is primarily governed by the mixing of two largely different water masses i.e. fresh river- and saline sea-water (Telesh and Khlebovich 2010). Hence, concentration of these dissolved constituents tends to fall on a straight line which is referred to as the theoretical dilution line (TDL) when plotted against the salinity. However, various biogeochemical processes are responsible for the removal or addition of the ionic constituents within the estuary that causes deviation of ionic concentrations from the TDL (Liss 1976). Based on this concept, the mixing behavior of major ions in the Devi River estuary has been evaluated and subsequently their percentage of addition or removal (PA/PR) were calculated using Eq. 2–4. The PA/PR of the major ions in Devi estuary are given in Table 2.
Table 2
Percentage addition and removal (PA/PR) of dissolved ions in the estuary.
|
SiO2
|
Na
|
K
|
Ca
|
Mg
|
Cl
|
SO4
|
NO3
|
PO4
|
Summer
|
Average
|
1.94
|
-4.86
|
-10.00
|
-3.03
|
-38.45
|
-7.59
|
4.45
|
-0.03
|
-2.45
|
SD
|
13.49
|
12.06
|
28.92
|
34.05
|
32.31
|
12.51
|
9.84
|
7.23
|
6.13
|
Monsoon
|
Average
|
0.60
|
-6.26
|
-5.86
|
167.63
|
9.62
|
3.52
|
4.55
|
-3.23
|
0.11
|
SD
|
6.51
|
8.63
|
27.79
|
164.93
|
31.85
|
4.00
|
16.53
|
7.88
|
0.68
|
Winter
|
Average
|
1.88
|
0.17
|
6.09
|
30.13
|
-2.46
|
-0.91
|
7.15
|
0.20
|
-0.07
|
SD
|
10.07
|
4.74
|
16.87
|
52.03
|
28.12
|
3.46
|
15.79
|
7.89
|
0.66
|
Where, Nmix and Nm are the theoretical and measured concentration of constituents at a particular salinity while Nr and No are the concentration of ion in river- and sea-water respectively. The Smix, Sr and So are the salinity in the sample, river- and sea-water respectively whereas ff is the mixing ratio of fresh- and saline-water in the estuary.
The concentration of Na and Cl are mostly plotted along the TDL, which indicates their conservative behavior throughout the studied period (Fig. 3–5; Table 2). Behavior of Na and Cl is considered to remain conservative in both seawater and estuarine water, as they generally do not take part in biochemical processes (Borole et al. 1979; Mora et al. 2017). Apart from these, the distributions of SO4, NO3 and PO4 in the estuary are also mostly conservative in the estuary. Hence, it can be inferred that distribution of these ions are mostly hydro-dynamically controlled by the mixing of fresh- and saline-water rather than any biogeochemical process. Conservative behavior of these major ions has been reported in several estuaries worldwide (Balls 1992; Patra et al. 2012). The DSi shows enrichment in the upper reaches of the estuary and gets depleted towards the mouth of the estuary. Intense weathering due to land use and/or land cover changes has been reported to increase the DSi concentrations in rivers (Hartmann et al. 2011). This may have led to the excess DSi concentrations in the upper estuary. However, as discussed in Section 4.2, incorporation of silica in the frustules of marine diatoms present in coastal water acts as a sink for the DSi supplied by fluvial input. The presence of centric diatoms in the estuarine sediment is confirmed in the SEM microphotograph and EDX spectra (Fig. 7a). Hence, biogenic removal occurring near the mouth of estuary can explain the depletion of DSi. This contrasting behavior of DSi between the upper and lower estuary is more prominent during the dry seasons (Fig. 3, 5). This is due to less supply of DSi to the coast because of decreased water flow. However, during monsoon, increased water discharge leads to continuous supply of DSi to the mouth of estuary. This compensates the removal of DSi by diatom uptake and causes less deviation from the TDL. Such non-conservative behavior of DSi in estuaries has been widely reported (Bell 1984; Amann et al. 2014; Akshitha et al. 2021).
Behavior of Ca and Mg along the salinity gradient shows marked variation in different seasons. During summer, negative deviation of Ca (PA/PR = -3.03 ± 34.05) and Mg (PA/PR = -38.45 ± 32.31) from the dilution line is observed (Fig. 3; Table 2). This depletion is caused due to the removal of Ca and Mg from water column as a result of precipitation of carbonates such as calcite, aragonite and dolomite. Precipitation has been found to be an important ongoing process in the estuary during the summer season (Fig. 6). Depletion in the concentration of K (-10.00 ± 28.92) relative to the TDL indicates its removal from water. This is possibly due to incorporation of K in the crystal lattices of carbonates (calcite and aragonite) (White 1977; Ishikawa and Ichikuni 1984). Positive correlations of K with calcite and aragonite further supports this inference (Table S1). However, K mostly shows conservative behavior during monsoon and winter. The monsoon period is characterized by significant non-conservative behavior of Ca in the estuary though Mg shows mostly conservative behavior in this season with slight enrichment at higher salinity. The enrichment of Ca and Mg coincides with the negative saturation indices during monsoon indicates that carbonates dissolution contributes the excess Ca and Mg ions in the surface water. However, the significant enrichment of Ca during this period cannot be explained by the dissolution of carbonate minerals solely. Hence, other factors such as groundwater discharge may have contributed to the enrichment of Ca in this period. Increase in groundwater input into estuaries have been reported during the monsoon period (Rahaman and Singh 2012; Selvam et al. 2021). Lagomasino et al. (2015) also observed excess Ca in the Sian Ka’an Biosphere Reserves estuary due to the combined effects of dissolution of carbonates and influx of Ca-rich groundwater into the estuary. Slight enrichment and depletion of Ca and Mg are observed in the upper and lower estuary, respectively in the winter period. Therefore, during winter, the SI shows mostly negative values in the upper estuary, and positive values in the lower estuary (Fig. 6). Dissolution and/or precipitation of carbonates have been reported to influence the non-conservative behavior of Ca and Mg in estuaries (Ramanathan et al. 1993; Jarvie et al. 2000; Patra et al. 2012).
4.5. Nutrient stoichiometry
Nutrient availability and their stoichiometry in estuarine water plays an important role in modulating the phytoplankton growth, biomass and species composition (Kilham 1990; Domingues et al. 2011). Spatial and temporal availability of nutrient can be changed drastically in estuary due to mixing of fresh river water and saline ocean water. In natural water, the stoichiometric ratio between the nitrate and phosphate (NO3:PO4) is found to have a constant value of 16:1, which is referred to as the Redfield ratio (Redfield 1958). Brzezinski (1985) further observed that DSi also acts as an important nutrient and exhibits constant ratios with nitrate and phosphate, i.e. DSi:PO4 = 10:1 and DSi:NO3 = 1:1. The availability of NO3, PO4 and DSi in these stoichiometric ratios are considered to be optimal for the growth of phytoplankton (Figler et al. 2021). The constancy of these nutrient ratios is maintained by the continuous biogeochemical cycling of nutrients into and out of their bioavailable forms. Limitation or non-availability of a particular nutrient has profound effect on the phytoplankton growth and their assemblage. Therefore, these ratios are important for assessing the condition of phytoplankton growth. Further, any perturbation in these ratios can be used to determine the nutrient limitation in localized areas (von Oheimb et al. 2010).
The NO3:PO4 ratio ranges from 5.38 to 52.89 in summer followed by 12.27 to 59.32 during winter and 11.26 to 60.36 in monsoon. The highest DSi:PO4 ratio is observed in monsoon season (55.45 ± 15.55) while the lowest in winter (13.26 ± 4.05). Estuarine water is characterized by DSi:PO4 ratio of 23.56 ± 19.40 in the summer period. During summer and winter, high NO3:PO4 values are observed in upper parts of the estuary whereas most of the middle and lower estuary samples show values closer to that of the Redfield ratio (Fig. 8a). Input from various anthropogenic sources accounts for high NO3 content in the upper estuary. However, microbial denitrification in the lower estuary removes NO3 from the water column and brings the NO3:PO4 ratio close to the Redfield value (Lenton and Watson 2000; Howarth and Marino 2006). Another explanation for the wide variation in NO3:PO4 values between upper and lower estuary could be the greater availability of PO4 towards estuary mouth (Caraco et al. 1989; Blomqvist et al. 2004). Under the anoxic conditions prevailing in the lower estuary, PO4 bound to Fe3+ oxide is easily released into the water column due to reduction of Fe3+ to Fe2+. Further, the Fe2+ precipitates in the form of sulphides in sediments underlying the sulphate rich seawater (Hartzell and Jordan 2012). This prevents the escape of Fe2+ to the overlying oxic layers and inhibits further complexing with PO4. Framboidal pyrite has been identified in sediments near mouth of the Devi estuary which supports the above geochemical process (Fig. 7b). However, low concentrations of sulphate in the fresh water causes the Fe2+ to form particulate ferrous compounds with PO4, and thereby reducing the availability of latter (Caraco et al. 1989).
Similar spatial distribution of DSi:PO4 ratio is also observed during summer and winter period. This indicates that productivity is limited by nutrient PO4 in the upper estuary during the dry periods (Fig. 8b). An increase in chlorophyll concentration towards the mouth of the estuary during dry period in Devi estuary has also been reported by Akhtar et al. (2021). The availability of nutrients in accordance to the Redfield ratio supports proliferation of primary producers in the lower estuary. Similar PO4 limitation in the fresh water dominated low salinity zone have been reported in some of the major estuaries such as Neva estuary (Russia), Pearl river estuary (China) and Chesapeake Bay (US) (Pitkänen and Tamminen 1995; Yin et al. 2004; Hartzell and Jordan 2012). Contrastingly, NO3:PO4 higher than the Redfield ratio has been observed during monsoon period throughout the estuary. Phosphate limitation has been widely reported in fresh river water (Howarth and Marino 2006). This is attributed to increase in NO3 because of higher nitrogen fixation rates and relatively lower availability of PO4 due to complexation with Fe-oxides in freshwater systems (Blomqvist et al. 2004). The estuary remains dominated by fresh water in monsoon season due to high river discharge. This leads to the influx of PO4-limited fresh water into the estuary and increases the NO3:PO4 ratio. Thus, PO4 is stoichiometrically the limiting nutrient during the monsoon season. The significantly high values of DSi:PO4 in rainy season also indicate a PO4-limitation in the estuary. This is further supported by the lowest chlorophyll concentration (0.35 ± 0.17 µg/l) in the Devi estuary during the monsoon period (Akhtar et al. 2021).
The DSi:NO3 ratio has primary control over the growth of diatom in estuaries (Egge and Aksnes 1992; Conley et al. 1993; Choudhury and Bhadury 2015). It is found to be 1.29 ± 0.21, 1.37 ± 0.26 and 0.52 ± 0.18 in summer, monsoon and winter respectively. Thus, this ratio is observed to be approximately 1:1 during summer and monsoon that suggests nutrient availability is favorable for diatom growth. However, low DSi:NO3 ratios in upper estuary during winter represents DSi-limitation.