# Dynamics of Dissolved Oxygen in Relation to Saturation and Health of an Aquatic Body: A Case for Chilka Lagoon, India

Authors: B. S. R. V. Prasad; P. D. N. Srinivasu; P. Sarada Varma; A. V. Raman; Santanu Ray
Journal: Journal of Ecosystems
 (2014)

---

## Abstract

Dissolved oxygen (DO) is essential for an aquatic ecosystem since it controls the biological productivity. Here, we propose a unidimensional dynamic model for DO by incorporating biological (photosynthesis, respiration, and mineralization), physical (atmospheric reaeration) and chemical (nitrification) processes so characteristic of shallow coastal water bodies. The analytical study of the proposed model is focussed on supersaturation and undersaturation of oxygen in the water body. The controllability of the ecosystem health has also been investigated. Model results indicate that, while undersaturation of oxygen is largely governed by nitrification and Net Ecosystem Metabolism (NEM), the supersaturation is controlled by photosynthetic activity. The model results are corroborated with observed data collected from Chilka lagoon, India. Subsequently, a biogeochemical model to study the DO variations in Chilka lagoon has been constructed. The model is properly calibrated and validated with observed data. Two independent sets of data (2004-2005 and 2005-2006) were used for model calibration and validation and Chi-square tests supported its robustness (R2=0.982 and 0.987; P<0.05). The model was used to evaluate independently the influence of individual taxa (diatoms, microphytobenthos, and cyanophyceans) on DO variations. Simulations indicate the vital role of microphytobenthos in lagoon DO dynamics and the overall wellbeing.


---

## Body

## 1. Introduction

Dissolved oxygen (DO) is an essential component which determines the water quality and trophodynamics of an aquatic system [1–3]. A fluctuation of DO near its saturation indicates relatively healthy waters [4]. In a healthy aquatic lagoon ecosystem the concentration of oxygen must always be above the chronic criterion for growth (4.8 mg L−1) developed by the U.S. EPA [5, 6]. A fall of oxygen below this chronic criterion can cause adverse effects on the health of the aquatic system. The system can suffer with hypoxia/anoxia conditions, resulting in depletion of fish stocks and other important aquatic animals [7, 8]. These losses can have significant damaging effects on ecological health, economic health, and stability of aquatic systems [8, 9].DO dynamics are complex in nature and are affected by many physical, chemical, and biological processes. The important factors that affect DO dynamics in an aquatic environment are atmosphere-water surface exchange, photosynthesis, respiration, and mineralisation [5, 10, 11]. The other processes that affect the oxygen dynamics in the water column include nitrification, sediment-water exchange. The availability of the oxygen in the water column also depends on topography of the system, the nature of the soil, and so forth ([4] and references there in).Many experimental studies are carried out to study the conditions that are responsible for hypoxia/anoxia in costal ecosystems [3, 12–14]. These studies pointed out that factors such as excess of nutrient load (eutrophic condition) and nitrification are the primary reasons for the hypoxic/anoxic conditions. To support these conclusions, many theoretical and empirical models (ranging from simple one-dimensional to three-dimensional) are constructed [10, 15–18]. All these studies concentrate on specific ecosystems (such as Cartagena Bay, Chesapeake Bay, Gulf of Mexico, Lake Zurich, and Virginia Estuaries, etc.) and to our knowledge no general theory exists which presents the influence of ecological and environmental factors on the qualitative behaviour of DO. Such a study would be useful not only to understand the strength of the external factors on the dynamics of DO but also to design control strategies to elevate the DO levels in the ecosystems.In the present paper we construct a universal model to represent the DO dynamics that takes into consideration the essential factors that influence the concentration of DO in aquatic systems. This work presents a novel approach to evaluate the conditions under which an aquatic system experiences super/under saturated conditions. This analysis also highlights the influence of Net Ecosystem Metabolism (a vital concept related to studies on aquatic ecosystems which is introduced by Odum [19]) on the saturation characteristics of the aquatic systems. The theoretical study indicates that the system will always be driven to undersaturated state whenever Net Ecosystem Metabolism (NEM) is negative. On the other hand, when NEM is positive the system can exist in any one of the saturated states depending on the strengths of other processes (such as mineralisation, nitrification, and reaeration). The model analysis also brings forward control strategies to steer the system from undersaturated state to supersaturated state. To our knowledge these types of analysis and results have not been presented in the literature so far.The main objective of the current study is to develop a universal model for DO dynamics and to bridge the theoretical results obtained from the abstract model with that of simulation modelling studies for better understanding of the ecosystem dynamics. This approach has been successfully applied to examine the causes and identify the crucial processes that are responsible for the saturation states observed in Chilka Lagoon (the largest brackish water lagoon in Asia, situated on the east-coast of India). By this study, it also has been possible to identify the key biogeochemical processes taking place in the Chilka Lagoon, based on which it is hypothesised that MPB are playing a key role in maintaining the health of the Chilka Lagoon.The paper is structured as follows. An abstract model representing the oxygen dynamics in a water column is formulated in next section. This is followed by a section that contains detailed analysis of the model which highlights conditions in terms of the system parameters that ensure the presence of oxygen in water column. In Section4 the conditions for the system to be super-/undersaturated are inferred and the results pertaining to system controllability are presented. The consequence of variations in the involved parameters on the health of the ecosystem is studied in Section 5. The theoretical findings are corroborated with the observations made for Chilka Lake in Section 6. In Section 7 appropriate forms have been attributed to the coefficient functions involved in the general model (by taking into consideration various biogeochemical events taking place in Chilka Lagoon), enabling simulation of DO trends observed in Chilka Lagoon and evaluation of vital parameters for the considered aquatic system. Finally, discussions and conclusions are presented in Section 8.


## 2. Model Description and Formulation

In this section we construct an abstract model to represent the dynamics of DO in an aquatic system by taking into consideration various important processes that influence DO in the (system) water column. These processes include atmosphere-water surface exchange, photosynthesis, nitrification, mineralisation of detritus, respiration by plant (phytoplankton, macrophytes, and microphytobenthos (MPB)) and animal (zooplankton, fish, etc.) community, and loss due to other processes. Figure1 presents a schematic diagram that describes the effect of these processes on DO in the water column.Figure 1
Schematic diagram presenting various important processes and their effects on the oxygen concentration in the water column.There are a few other significant pathways which can also influence the oxygen dynamics in aquatic systems. These include ground water incursion and entrainment of low DO from the lake sediments. These processes are hard to measure in times, in particular, the DO contribution through the ground water. Moreover, the addition of entrainment of DO from sediments increases model dimension from one to two which makes the system relatively complex and hard to track. To avoid this complexity, we have added a loss/gain term to our model and calibrated this term to fit the model results. This seems to be a convenient approach as it gives a rough estimate on how much of DO is being lost/gained due to other processes which are hard to estimate using available tools. The other important aspect which influences the DO is stratification of water column in certain periods of time. In our model we have neglected this effect due to the fact that Chilka is a shallow lagoon (mean depth< 1.5 m), which is seen through absence of any appreciable column variability in water temperature and/or salinity. Often, the strong surface winds overpower vertical gradients, if any, in the lake.IfO stands for DO concentration in the water column, then the following equation represents its associated dynamics:
(1)dOdT=±(Atmosphere-watersurfaceexchange)+Photosynthesis-(Nitrification+Respirationmmll+Mineralisation+Lossduetootherprocesses).The mathematical formulation of the various processes involved in the above equation and the description of parameters used in the model formulation are presented in Tables1 and 2, respectively.Table 1
Mathematical description of various processes involved in (1).
Process
Mathematical formulation
References
Atmosphere-water surface exchange
K
a
(
O
ss
-
O
)
[1, 4, 5, 10, 11, 15, 18]
Photosynthesis
K
2
G
P
P
[1, 5, 10, 15, 37]
Nitrification
K
1
O
A
/
(
K
Nit
+
O
)
[15, 18]
Respiration
K
3
P
+
K
4
Z
[1, 10]
Mineralisation
K
5
D
P
P
O
+
K
6
D
Z
Z
O
[37]
Loss due to other processes
K
Loss
OTable 2
Parameters of the model equation along with their units.
Parameter
Description
Units
O
Oxygen concentration in water column
mgO2L-1
Oss
Saturating oxygen concentration
mgO2L-1
K
2
Photosynthesis rate constant
(
mg
O
2
L
-
1
)
×
(
mgC
L
-
1
)
-
1
G
Phy
Specific growth rate of plant community
tim
e
-
1
P
Biomass of aquatic plant community
mgC
L
-
1
K
a
Reaeration constant
tim
e
-
1
A
Ammonium nitrogen concentration
mgN
L
-
1
K
1
Nitrification rate constant
(
mg
O
2
L
-
1
)
×
(
mgN
L
-
1
)
-
1
×
tim
e
-
1
K
Nit
Half-saturation constant for oxygen  limitation of nitrification
mg
O
2
L
-
1
K
3
Respiration rate constant for aquatic  plant community
(
mg
O
2
L
-
1
)
×
(
mgC
L
-
1
)
-
1
×
tim
e
-
1
K
4
Respiration rate constant for aquatic  animal community
(
mg
O
2
L
-
1
)
×
(
mgC
L
-
1
)
-
1
×
tim
e
-
1
Z
Biomass of aquatic animal community
mgC
L
-
1
K
5
Mineralization rate constant for aquatic  plant community
(
mgC
L
-
1
)
-
1
K
5
Mineralization rate constant for aquatic  animal community
(
mgC
L
-
1
)
-
1
D
P
Mortality rate constant for aquatic  plant community
tim
e
-
1
D
Z
Mortality rate constant for aquatic  animal community
tim
e
-
1
K
Loss
Loss rate oxygen due to other processes
tim
e
-
1Incorporating the above mathematical expressions presented in Table1 into (1), we obtain the following one-dimensional model representing the dynamics of DO in the water column:
(2)dOdT=Ka(Oss-O)+K2GPP-K1AO(KNit+O)-K3P-K4Z-K5DPPO-K6DZZO-KLossO.The above equation can be conveniently rearranged to obtain the following simpler representation:(3)dOdT=-K1AO(KNit+O)-S1O+S2-S3+KaOss,
where
(4)S1=K5DPP+K6DZZ+Ka+KLoss,(5)S2=K2GPP,(6)S3=K3P+K4Z.
Here the expression S1 represents efflux of DO due to mineralisation of dead organic matter, exchange with atmosphere, and other processes, S2 represents the influx of DO by photosynthetic activity, and S3 represents the efflux of DO due to respiration by the plant and animal communities.


## 3. Model Analysis

The general model formulated in the previous section is analysed to study the existence and nonexistence of equilibrium solution for the considered system. This study enables the identification conditions under which the system moves towards oxic/anoxic state.To minimise the complexity in the analysis we first reduce the number of parameters in the model (3) by the process of nonimpersonalisation using the following transformations:
(7)x=OKNit,t=KaT.
These transformations reduce the model (3) to the form
(8)dxdt=-ax(1+x)-bx+c,
where
(9)a=K1AKaKNit,(10)b=S1Ka,(11)c=(S2-S3)(KaKNit)+OssKNit=c1+OssKNit.From the model description presented in the previous section and the formulation presented in (9), the positive parameter a is directly proportional to the maximum rate of nitrification and ammonium concentration in the system. The cumulative effect of loss due to mineralisation process, exchange with atmosphere, and other processes is represented by b. It can be easily observed that this parameter b is always greater than unity (cf. (4), (10)). The parameter c is given to be the sum of two terms given by c1 and Oss/KNit (cf. (11)). While the latter term is always positive, the former term represents the net ecosystem metabolism (NEM) [20] and can take either sign. Here NEM represents the net effect of production (due to photosynthesis) and consumption (due to respiration by all communities).If the net gain due to photosynthetic activity and diffusion from atmosphere is greater than the loss due to respiration, then we havec to be positive else it is negative. Therefore positivity of NEM (c1) implies positivity of c and negativity of c implies negativity of NEM (c1) as the term Oss/KNit is always positive.The model analysis gives us valuable information regarding the cause and relations of various events that take place in water bodies. It also helps in deriving solutions for vital problems such as avoidance of anoxic conditions to prevent ecological and economic losses in the ecosystem. As a part of this analysis we shall first study the existence of equilibrium solutions for the considered model and investigate their qualitative behaviour. The considered model admits at most one positive equilibrium solution given by(12)x*=(-[a+b-c]+[a+b-c]2+4bc)2b.The stability of the equilibrium depends on the sign of the derivative of-(ax)/(1+x)-bx+c (cf. RHS of (8)) evaluated at the equilibrium. This equilibrium will be stable if this sign is nonpositive and unstable otherwise [21]. Further this equilibrium is asymptotically stable if the derivative sign is negative. Clearly, the existence and stability nature of the equilibrium solution depend on the involved parameters a, b, and c. We know that a is always positive, b is greater than unity, and c can take any sign. Thus for a given b>1, it is enough to study the qualitative behaviour of (8) in the first and fourth quadrants of ac-space.Figure2 presents the relationship between these parameters and the existence of equilibrium solution. Note that the curves a+b-c=0 and c=0 divide the first and fourth quadrants of ac-space into three significant regions given by
(13)RegionI:={(a,c):a+b-c<0,c>0},RegionII:={(a,c):a+b-c>0,c>0},RegionIII:={(a,c):a+b-c>0,c≤0}.Figure 2
Figure representing the relationship between the parametersa, b, and c and the existence of equilibrium solution. In Region I and II, the system admits unique equilibrium. The system admits no equilibrium in Region III. The magnitude of equilibrium DO is higher in Region I when compared to Region II. This is due to domination of the processes present in c over the processes that are involved in a and b.It is easy to observe that (8) admits a unique globally asymptotically stable equilibrium solution when the parameters belong to Regions I and II and it does not admit any equilibrium solution in Region III. We have the following result.Theorem 1.
For a givenb>1, (8) admits a unique globally asymptotically stable equilibrium if the parameters a and c belong to the first quadrant of the ac-space and it admits no equilibrium solution if they belong to the fourth quadrant including c=0.Theorem1 is a significant result in the sense that it gives the essential conditions under which the water column supports DO. This result along with the expression for the equilibrium solution (12) and Figure 2 helps us to derive conditions for improving DO levels in the water column.In Region III we havec to be negative indicating negativity of NEM for the aquatic system. In this case the system loses oxygen and gets driven to anoxic conditions as cumulative gain due to the fact that photosynthetic and atmospheric diffusion processes are dominated by cumulative loss due to respiration, nitrification, and mineralisation.In Regions I and II we havec to be positive representing the dominance of photosynthetic and atmospheric diffusion processes over the process of respiration. In this case the system supports DO (reflected by the existence of asymptotically stable equilibrium) but its magnitude depends upon strength of c relative to a+b.If the system parameters belong to Region I, then the characteristic of the aquatic system is such that the activities of photosynthesis and atmospheric diffusion dominate not only the respiration process (i.e.,c>0) but also the activities of nitrification and mineralisation (i.e., a+b<c). In this case the aquatic system retains higher levels of oxygen (owing to higher magnitude of the equilibrium).In Region II although the photosynthetic and atmospheric diffusion dominate the respiration process (i.e.,c>0), the cumulative effect of respiration, nitrification, and mineralisation dominates the cumulative effect of photosynthesis and atmospheric diffusion processes (i.e., a+b>c). Thus the system could tend to hypoxic conditions (reflected by the smaller magnitude of the equilibrium).Here, it is important to note that the level of DO continues to decrease as the value ofc decreases. The system gets driven to hypoxia followed by anoxia in the vicinity of c=0 as c crosses the a-axis. Further decrease in c from this level (c=0) makes the system free from DO. These findings are in accordance with the experimental studies carried out on hypoxia/anoxia [8, 13, 14].


## 4. Regions of Supersaturation, Undersaturation

In this section a novel approach is presented to derive the conditions under which the equilibrium state of DO will be in two possible states given by supersaturation or undersaturation. By definition, oxygen in a water body is said to be super- (under)saturated if DO concentration is greater (lesser) than the saturating oxygen in the water body. To study the conditions under which the DO will be super- (under)saturated, we first identify the ranges on the parameters that lead the system to the respective saturation states. For this we proceed as follows.For a given saturating oxygen level Oss of the water body, the equilibrium state x (under the assumption that it exists) of the system (8) satisfies
(14)ax(1+x)+bx-c1=OssKNit.Herex is nothing but O/KNit, a proportion of the DO concentration and c1 represents the NEM of the system. Treating Oss as a variable and defining y=Oss/KNit let us consider the following equations:
(15)y=ax(1+x)+bx-c1≡f(x),(16)y=x,
in the positive quadrant of the xy-space.The curve (15) is a representative of the set of all admissible equilibrium solutions of (8). Henceforth we shall label this curve (15) as admissible equilibrium curve. Note that the equilibrium of the system for a given saturating oxygen Oss is given by f-1(Oss/KNit)KNit.Observe that the line (16) divides the positive quadrant into two regions such that one of them (with x<y) represents region of undersaturation and the other (with x>y) represents region of supersaturation. Figures 3 and 4 present the curves (15), (16) in the positive quadrant of the xy-space for the cases c1>0 and c1<0, respectively. These figures give further information on the saturation characteristics of the equilibrium state and associated controllability aspects.Figure 3
Regions of supersaturation and undersaturation for the casec1>0(NEM>0). The dashed line represents the line (16), which is the boundary between the supersaturation and undersaturation regions. The solid line represents the curve (15). (xc,yc) is the intersection point of (15) and (16). Observe here that the system equilibrates at super- (under) saturated state if y<(>)yc.Figure 4
Regions of supersaturation and undersaturation for the case whenc1<0(NEM<0). The dashed line represents the line (16), which is the boundary between the supersaturation and undersaturation regions. The solid line represents the curve (15). Observe here that (15) and (16) do not admit any intersection and the admissible equilibrium curve lies completely in undersaturated region.Note that the curves (15) and (16) intersect in the interior of the positive quadrant only if c1>0 (Figure 3). On the contrary, if c1 is nonpositive, then there does not exist any such intersection in the positive quadrant and the curve (15) lies completely above (16) (cf. Figure 4). Observe that the admissible equilibrium curve (15) is a hyperbola and intersects y-axis at -c1 having a positive slope b+a/(1+x)2. This curve tends to infinity as x→∞. The slope of this curve highlights the vital role the mineralisation and nitrification processes play in improving the equilibrium state of the water body. This aspect will be discussed in the next section.In order to study the relationship between the eventual state of oxygen with the saturating oxygen level, let us consider the intersection point of the curves (15) and (16) in the positive quadrant of xy-space. The x component of this intersection point is given by the positive root (xc) of the quadratic equation:
(17)(b-1)x2+[a+(b-1)-c1]x-c1=0,
where
(18)xc=(12(b-1))×{[a+(b-1)-c1]2+4c1(b-1)-[a+(b-1)-c1]+[a+(b-1)-c1]2+4c1(b-1)}.Let(xc,yc)(=(xc,xc)) be the intersection of the curves (15) and (16) which is represented in Figure 3. Thus for positive NEM (i.e, c1>0), the admissible equilibriums curve can be divided into two significant segments given by
(19)Segment1={(x,y):y=f(x)withy<yc},Segment2={(x,y):y=f(x)withy>yc}.IfO¯ss is the saturation level of the water body with y¯=O¯ss/KNit and x¯=O¯/KNit is the corresponding equilibrium of (8), then (x¯,y¯)∈ Segment 1 (Segment 2) implies that the water body will be in super- (under)saturation state eventually.For the case of negative NEM (i.e.,c1<0) the system always equilibrates in undersaturated state as the curve (23) lies completely in the region of undersaturation (cf. Figure 4). From these observations, we conclude that while a water body may equilibrate in a supersaturated state or undersaturated state when the NEM is positive, it always equilibrates in undersaturated state if NEM is negative.


## 5. Ecosystem Health and System Controllability

In this section we study the consequences of variations in the involved parameters (a, b, and c) on the health of the aquatic system. We also find strategies to drive the equilibrium state of the system from undersaturation to supersaturation state. From the graph of (15) (Figures 3 and 4), observe that the DO equilibrium level decreases as the saturating oxygen decreases. Also from the qualitative properties of this curve (in particular slope) it is easy to see that for a given saturating oxygen concentration the corresponding equilibrium DO concentration increases with decrease in either a or b or both. This clearly indicates that the DO concentration in the system decreases with an increase in the ammonium concentration in the system (which increases a) or increase in the mineralisation activity (which increases b) or both.From the analysis presented in the previous section, it is clear that an equilibrium state of the considered system cannot be driven to supersaturation state when NEM is negative. Thus the possibility of transition from undersaturation to supersaturation exists only when NEM is positive. Observe that, for a given representative of saturating oxygen (y) if the corresponding DO equilibrium (x) satisfies f(x)>yc(f(x)<yc), then the system would be under- (super)saturated at the equilibrium.Thus given that the system is undersaturated at a given saturation levely, the way to drive the system from undersaturated state to supersaturated state is to alter the function f in such a way that the corresponding DO equilibrium state, denoted by xa, satisfies the supersaturation condition (f(xa)<yc). The monotonicity property of f(x) ensures such a possibility through a sufficient reduction in its slope (which is b+a/(1+x)2). This can be achieved by reducing the coefficients a and/or b sufficiently. It can be easily observed that reduction in b decreases the slope of f(x) faster when compared to reduction in a.For the purpose of illustration, we assume the following hypothetical values (but, these values can also be found using field experiments) for the nondimensional parametersa, b, and c of the system (8):
(20)a=8.00,b=2.00,c1=9.35.
The admissible equilibrium curve (15) with the foresaid parameter values is presented in Figures 5 and 6 with dashed lines. Taking the saturation level (y) of the system to be 7, we observe that the DO in the system equilibrates in undersaturated state with a magnitude of 4.85784 (cf. Figures 5 and 6).Figure 5
Figure illustrating the possibility of driving the equilibrium state of the system from undersaturation to supersaturation state by decreasing the ammonium concentration (i.e., parametera). The DO concentration in the water column is plotted on x-axis and the saturating oxygen concentration is plotted on y-axis. The dashed line and the dot-dash line represent the admissible equilibrium curves with (a,b,c1)=(8.00,2.00,9.35) and (0.4,2.00,9.35), respectively. Under the assumption that y=7, the system with a=8 equilibrates in undersaturated state (x=4.85784) and it equilibrates in supersaturated state (x=8) when this parameter is reduced to 0.4.Figure 6
Figure illustrating the possibility of driving the equilibrium state of the system from undersaturation to supersaturation state by decreasing the mineralisation (i.e., parameterb). The DO concentration in the water column is plotted on x-axis and the saturating oxygen concentration is plotted on y-axis. The dashed line and the dot-dash line represent the admissible equilibrium curves with (a,b,c1)=(8.00,2.00,9.35) and (8.0,1.155,9.35), respectively. Under the assumption that y=7, the system with b=2 equilibrates in undersaturated state (x=4.85784) and it equilibrates in supersaturated state (x=8) when this parameter is reduced to 1.155.Following the theory, the system can be driven to a supersaturated state (sincec1>0), by varying the parameter a or b or both. Fixing the values of b and c1 it can be easily inferred from (18) that the system will equilibrate in supersaturated state if a is made to stay below 2.6857. Similarly, fixing a and c1, the system can be made to equilibrate in supersaturated state by ensuring that b is less than 1.3357.Suppose we wish that system should equilibrate in supersaturated state withx=8.00, (18) guarantees such eventual state for the choices of parameters (a,b,c1)=(0.4,2.0,9.35) or (8.0,1.155,9.35). The simulations demonstrating the move from undersaturation state to supersaturation state for both the above said choices are presented in Figures 5 and 6. The required reductions in the parameters a and b to reach a state with x=8 clearly demonstrate the stronger effect of reducing b over a in driving the state to supersaturation.In the case, wherec1<0 (i.e., system with negative NEM), the system can be driven from anoxic condition to hypoxic and further to conditions of fairly high DO concentration by reducing a and b sufficiently. However, it would not be possible to drive the system to supersaturated state by altering a and b.Thus, irrespective of the sign of NEM, reduction in ammonium concentration(a) and mineralisation activity (b) enhances the health of the ecosystem and enriches the waters with oxygen. At this juncture, it is worth noting that, for the case c1>0 it is also possible to drive the system to supersaturation state by reducing the saturating oxygen level to a state below yc. One of the ways to reduce the saturating oxygen level in the water body is to increase the salinity of the system; other possibilities are increase in temperature and decrease in partial pressure [5]. Thus increasing the salinity of the water body is yet another solution to increase the health of the water bodies.From (9)–(11), it can be noted that the values of the parameters a, b, and c1 are dependent on the factors such as reaeration coefficient, ammonium concentration, nitrification rate constant, half-saturation constant for nitrification, plankton biomass, mortality rates of plankton groups, and losses due to other processes, and so forth. Having knowledge of all these parameter (either through field surveys or experimental studies) for a real-life system under consideration, it is possible to determine the values of a, b, and c1, from which we can derive the management strategies to drive the system to supersaturation conditions. Please refer to Section 7.3 and Figure 12 for evaluation of these coefficient values and their seasonal dependent variations for the case of Chilka Lagoon.


## 6. A Few Lessons in the Light of the Analysis

Analysis of the abstract model presented in the previous sections can be used to answer several questions pertaining to water bodies such as Chilka Lake, a brackish water lagoon on the east-coast of India. In this section the events that are responsible for the observed DO levels in Chilka Lagoon during 2004-2005 are interpreted using the analysis presented in previous sections.


### 6.1. Study Area

Chilka Lake (19°28′–19°54′N 85°06′–85°36′E) (Figure 7) is situated on the eastern seaboard of India some 350 km south of Kolkata. The lagoon is pear shaped and covers a total area of about 1000 km2 during wet season (August–October), but the volume is reduced by nearly 60% during summer months (April-May). Topographically, Chilka Lake can be divided into two major regions—the main area and the outer channel. The main area, which is about 65 km long and 20 km at its greatest breadth, constitutes the lagoon proper and can be arbitrarily divided into south, central, and north sectors. In general, the lagoon is shallow (~1.5 m) particularly in the north sector where considerable silting has taken place. Overall, the topography is not even since large variations in depth occur as a result of heavy deposition of terrigenous materials from catchments and river discharge. In the outer channel, which connects the lagoon proper to the sea (Bay of Bengal), depth variations are even marked owing to erosion and accretion processes. There are many small and large islands scattered throughout Chilka Lagoon. The largest is the Nalaban (~35 km2), which is covered by considerable marsh vegetation. During monsoon when the lagoon receives enormous quantities of freshwater from the catchments and surrounding rivers, most of these islands are inundated. However, during summer when freshwater discharge is negligible and evaporation exceeds precipitation, the lagoon recedes to a considerable distance exposing a vast expanse of wetland area. Chilka Lagoon is unique not only from the point of its shear extent but also hydrologically as it is nature's most intriguing habitat supporting a wealth of marine, brackish, and freshwater life. In the year 1981 Chilka Lake was designated as a Ramsar Site—a wetland of international importance, because of its rich biodiversity and socioeconomic value.Figure 7
Chilka Lake—the study area.


### 6.2.  Data Collection Procedures

#### 6.2.1. Depth Measurements

Time series observation on the bathymetry of the lagoon has been obtained from 36 GPS fixed locations (cf. Figure7) at frequent intervals of time (every 15-day interval) during the period of study May 3, 2004 to September 22, 2006 using the conventional depth measurement technique (weighted line and a graduated bamboo pole) described by [22].


#### 6.2.2. Salinity Measurements and Local Weather Data

Salinity was determined according to Knudsen’s method [23] at all the 36 GPS locations during the study period on a monthly basis. A complex combination of freshwater discharge, evaporation, wind condition, and tidal inflow of seawater govern the spatial and temporal changes in salinity of the lagoon. In October the entire lagoon is occupied with freshwater. Salinities between 1.8 to 4.6 PSU in the main area and between 7 to 8 PSU in the outer channel (cf. Figure 7) were noticed. In February salinities rise (5–20 PSU) due to decrease in freshwater flow and they are maximum (16.5–34.7 PSU) in May. Meteorological data (rainfall, wind speed, atmospheric temperature, and relative humidity), recorded at Gopalpur station (some50 km south of Chilka), were obtained from Indian Meteorological Department, Bhubaneswar. The annual maximum temperature (34.8°C) is recorded in May and minimum temperature (15.6°C) in December. The annual rainfall in the region is 1123 mm. Rainfall is very high (151–298 mm) during June–September and the least (0.4–4.9 mm) during January–March. The lagoon experiences Southwest and Northeast monsoons during June–September and November-December, respectively. Wind speeds are high (6.5–9m/s) during monsoon months (March–August) and low (4.8–6.4 m/s) during winter months (January-February). Wind direction is mostly N and NE during postmonsoon months (October-November) and SW during monsoon months. During the study period the mean monthly atmospheric minimum temperature was 21.51°C and maximum temperature was 30.29°C; mean total rainfall is amounted to 1514.1 mm the bulk of which 72.2% was recorded during the months July–September of every year corresponding to the southwest monsoon.


### 6.3. Model Application to Chilka

Figure8 presents the Lake-wide mean DO (along with computed DO saturation), ammonium concentration, and phytoplankton biomass for the period July 19, 2004 to July 18, 2005. The DO saturation for Chilka is computed using the formula [5]
(21)ln⁡(Osf)=-139.34411+(1.575701×105)Ta-(6.642308×107)Ta2+(1.243800×1010)Ta3-(8.62149×1011)Ta4,ln⁡(Oss)=ln⁡(Osf)-S((2.1407×103)Ta21.7674×10-2-(1.0754×101)Ta+(2.1407×103)Ta2),
where Oss is the saturating oxygen concentration in salt water at 1 atm [mg L−1], Osf is the saturating oxygen concentration in fresh water at 1 atm [mg L−1], Ta is absolute temperature [K], and S is the salinity [ppt] of the lagoon.Figure 8
Field observations of oxygen concentration in water column along with the computed saturating oxygen, ammonium concentration, nitrate concentration, phytoplankton biomass, depth of water column, and inflow from rivers for the period July 19, 2004 to July 18, 2005. Region A in the figure represents the active flow period (September 2004 to November 2004). Region B in the figure represents the lean flow period (February 2005 to April 2005).The period September–November, 2004, in Figure8 represents the active river flow period and the period February–April, 2005, represents the lean flow conditions. Here, active flow period represents the period during which the lake receives enormous amount of freshwater from the catchments.During the period of active flow, it is observed that the lake experiences high turbidity (mean 34.25 NTU) and low salinity (mean 9.44 PSU) conditions. In this period the findings of Gupta et al. [24] reported that the NEM of the Chilka Lagoon is negative, which is attributed to reduced activity of phytoplankton and MPB. According to the theory developed in Sections 3 and 4, the lake is supposed to experience undersaturation conditions during active flow period, where the NEM is negative. The observations on DO presented in Figure 8 confirm this fact. The fall in ammonium concentrations during this period indicates the activity of nitrification, which further reduces the oxygen concentration.During lean flow period, Gupta et al. [24] reported that the NEM of the Chilka is positive, which is due to high photosynthetic activity by autotrophic organisms. Again, according to the theory the lake can experience either of the saturation states during lean flow periods depending on the ambient conditions. Here undersaturated conditions prevail if the processes of nitrification and mineralisation dominate NEM. However supersaturation conditions can be observed if the photosynthetic activity by autotrophic organisms dominates the processes of nitrification and mineralisation.From Figure8 it may be seen that in the first half of the lean flow period the phytoplankton activity dominates nitrification and mineralisation processes. Thus the system experiences supersaturated conditions, validating the concepts presented in the model. In the latter half of the lean flow period, increase in ammonium concentrations may be attributed to mineralisation activity. This observation may be substantiated by the decrease in phytoplankton concentration following the bloom. Thus the oxygen concentrations during this period almost reach the level of saturation, which is also in accordance with the theory developed in Sections 3 and 4.Coming to intermittent conditions (particularly during the periods August 2004, November 15, 2004 to January 14, 2005), we notice that the system is in supersaturated state even though the phytoplankton activity is low. From the theory, supersaturated state is possible only when NEM is positive. So it would appear that, during the said periods, the low phytoplankton activity is compensated by the activities of other autotrophic organisms such as seagrasses and MPB which seem to be driving NEM to a positive state. This is an important observation from the perspective of ecosystem health. This indicates that the seagrasses, MPB, and other plant community (other than phytoplankton) also play a vital role in maintaining the health of the Chilka Lagoon.


## 6.1. Study Area

Chilka Lake (19°28′–19°54′N 85°06′–85°36′E) (Figure 7) is situated on the eastern seaboard of India some 350 km south of Kolkata. The lagoon is pear shaped and covers a total area of about 1000 km2 during wet season (August–October), but the volume is reduced by nearly 60% during summer months (April-May). Topographically, Chilka Lake can be divided into two major regions—the main area and the outer channel. The main area, which is about 65 km long and 20 km at its greatest breadth, constitutes the lagoon proper and can be arbitrarily divided into south, central, and north sectors. In general, the lagoon is shallow (~1.5 m) particularly in the north sector where considerable silting has taken place. Overall, the topography is not even since large variations in depth occur as a result of heavy deposition of terrigenous materials from catchments and river discharge. In the outer channel, which connects the lagoon proper to the sea (Bay of Bengal), depth variations are even marked owing to erosion and accretion processes. There are many small and large islands scattered throughout Chilka Lagoon. The largest is the Nalaban (~35 km2), which is covered by considerable marsh vegetation. During monsoon when the lagoon receives enormous quantities of freshwater from the catchments and surrounding rivers, most of these islands are inundated. However, during summer when freshwater discharge is negligible and evaporation exceeds precipitation, the lagoon recedes to a considerable distance exposing a vast expanse of wetland area. Chilka Lagoon is unique not only from the point of its shear extent but also hydrologically as it is nature's most intriguing habitat supporting a wealth of marine, brackish, and freshwater life. In the year 1981 Chilka Lake was designated as a Ramsar Site—a wetland of international importance, because of its rich biodiversity and socioeconomic value.Figure 7
Chilka Lake—the study area.


## 6.2.  Data Collection Procedures

### 6.2.1. Depth Measurements

Time series observation on the bathymetry of the lagoon has been obtained from 36 GPS fixed locations (cf. Figure7) at frequent intervals of time (every 15-day interval) during the period of study May 3, 2004 to September 22, 2006 using the conventional depth measurement technique (weighted line and a graduated bamboo pole) described by [22].


### 6.2.2. Salinity Measurements and Local Weather Data

Salinity was determined according to Knudsen’s method [23] at all the 36 GPS locations during the study period on a monthly basis. A complex combination of freshwater discharge, evaporation, wind condition, and tidal inflow of seawater govern the spatial and temporal changes in salinity of the lagoon. In October the entire lagoon is occupied with freshwater. Salinities between 1.8 to 4.6 PSU in the main area and between 7 to 8 PSU in the outer channel (cf. Figure 7) were noticed. In February salinities rise (5–20 PSU) due to decrease in freshwater flow and they are maximum (16.5–34.7 PSU) in May. Meteorological data (rainfall, wind speed, atmospheric temperature, and relative humidity), recorded at Gopalpur station (some50 km south of Chilka), were obtained from Indian Meteorological Department, Bhubaneswar. The annual maximum temperature (34.8°C) is recorded in May and minimum temperature (15.6°C) in December. The annual rainfall in the region is 1123 mm. Rainfall is very high (151–298 mm) during June–September and the least (0.4–4.9 mm) during January–March. The lagoon experiences Southwest and Northeast monsoons during June–September and November-December, respectively. Wind speeds are high (6.5–9m/s) during monsoon months (March–August) and low (4.8–6.4 m/s) during winter months (January-February). Wind direction is mostly N and NE during postmonsoon months (October-November) and SW during monsoon months. During the study period the mean monthly atmospheric minimum temperature was 21.51°C and maximum temperature was 30.29°C; mean total rainfall is amounted to 1514.1 mm the bulk of which 72.2% was recorded during the months July–September of every year corresponding to the southwest monsoon.


## 6.2.1. Depth Measurements

Time series observation on the bathymetry of the lagoon has been obtained from 36 GPS fixed locations (cf. Figure7) at frequent intervals of time (every 15-day interval) during the period of study May 3, 2004 to September 22, 2006 using the conventional depth measurement technique (weighted line and a graduated bamboo pole) described by [22].


## 6.2.2. Salinity Measurements and Local Weather Data

Salinity was determined according to Knudsen’s method [23] at all the 36 GPS locations during the study period on a monthly basis. A complex combination of freshwater discharge, evaporation, wind condition, and tidal inflow of seawater govern the spatial and temporal changes in salinity of the lagoon. In October the entire lagoon is occupied with freshwater. Salinities between 1.8 to 4.6 PSU in the main area and between 7 to 8 PSU in the outer channel (cf. Figure 7) were noticed. In February salinities rise (5–20 PSU) due to decrease in freshwater flow and they are maximum (16.5–34.7 PSU) in May. Meteorological data (rainfall, wind speed, atmospheric temperature, and relative humidity), recorded at Gopalpur station (some50 km south of Chilka), were obtained from Indian Meteorological Department, Bhubaneswar. The annual maximum temperature (34.8°C) is recorded in May and minimum temperature (15.6°C) in December. The annual rainfall in the region is 1123 mm. Rainfall is very high (151–298 mm) during June–September and the least (0.4–4.9 mm) during January–March. The lagoon experiences Southwest and Northeast monsoons during June–September and November-December, respectively. Wind speeds are high (6.5–9m/s) during monsoon months (March–August) and low (4.8–6.4 m/s) during winter months (January-February). Wind direction is mostly N and NE during postmonsoon months (October-November) and SW during monsoon months. During the study period the mean monthly atmospheric minimum temperature was 21.51°C and maximum temperature was 30.29°C; mean total rainfall is amounted to 1514.1 mm the bulk of which 72.2% was recorded during the months July–September of every year corresponding to the southwest monsoon.


## 6.3. Model Application to Chilka

Figure8 presents the Lake-wide mean DO (along with computed DO saturation), ammonium concentration, and phytoplankton biomass for the period July 19, 2004 to July 18, 2005. The DO saturation for Chilka is computed using the formula [5]
(21)ln⁡(Osf)=-139.34411+(1.575701×105)Ta-(6.642308×107)Ta2+(1.243800×1010)Ta3-(8.62149×1011)Ta4,ln⁡(Oss)=ln⁡(Osf)-S((2.1407×103)Ta21.7674×10-2-(1.0754×101)Ta+(2.1407×103)Ta2),
where Oss is the saturating oxygen concentration in salt water at 1 atm [mg L−1], Osf is the saturating oxygen concentration in fresh water at 1 atm [mg L−1], Ta is absolute temperature [K], and S is the salinity [ppt] of the lagoon.Figure 8
Field observations of oxygen concentration in water column along with the computed saturating oxygen, ammonium concentration, nitrate concentration, phytoplankton biomass, depth of water column, and inflow from rivers for the period July 19, 2004 to July 18, 2005. Region A in the figure represents the active flow period (September 2004 to November 2004). Region B in the figure represents the lean flow period (February 2005 to April 2005).The period September–November, 2004, in Figure8 represents the active river flow period and the period February–April, 2005, represents the lean flow conditions. Here, active flow period represents the period during which the lake receives enormous amount of freshwater from the catchments.During the period of active flow, it is observed that the lake experiences high turbidity (mean 34.25 NTU) and low salinity (mean 9.44 PSU) conditions. In this period the findings of Gupta et al. [24] reported that the NEM of the Chilka Lagoon is negative, which is attributed to reduced activity of phytoplankton and MPB. According to the theory developed in Sections 3 and 4, the lake is supposed to experience undersaturation conditions during active flow period, where the NEM is negative. The observations on DO presented in Figure 8 confirm this fact. The fall in ammonium concentrations during this period indicates the activity of nitrification, which further reduces the oxygen concentration.During lean flow period, Gupta et al. [24] reported that the NEM of the Chilka is positive, which is due to high photosynthetic activity by autotrophic organisms. Again, according to the theory the lake can experience either of the saturation states during lean flow periods depending on the ambient conditions. Here undersaturated conditions prevail if the processes of nitrification and mineralisation dominate NEM. However supersaturation conditions can be observed if the photosynthetic activity by autotrophic organisms dominates the processes of nitrification and mineralisation.From Figure8 it may be seen that in the first half of the lean flow period the phytoplankton activity dominates nitrification and mineralisation processes. Thus the system experiences supersaturated conditions, validating the concepts presented in the model. In the latter half of the lean flow period, increase in ammonium concentrations may be attributed to mineralisation activity. This observation may be substantiated by the decrease in phytoplankton concentration following the bloom. Thus the oxygen concentrations during this period almost reach the level of saturation, which is also in accordance with the theory developed in Sections 3 and 4.Coming to intermittent conditions (particularly during the periods August 2004, November 15, 2004 to January 14, 2005), we notice that the system is in supersaturated state even though the phytoplankton activity is low. From the theory, supersaturated state is possible only when NEM is positive. So it would appear that, during the said periods, the low phytoplankton activity is compensated by the activities of other autotrophic organisms such as seagrasses and MPB which seem to be driving NEM to a positive state. This is an important observation from the perspective of ecosystem health. This indicates that the seagrasses, MPB, and other plant community (other than phytoplankton) also play a vital role in maintaining the health of the Chilka Lagoon.


## 7. Development of Biogeochemical Model for DO in Lake Chilka

In this section, the general model discussed so far is applied to Chilka Lagoon by attributing appropriate forms to the involved coefficient functions. These functions are derived by taking into consideration various biogeochemical events taking place in the lagoon. Throughout this modelling exercise carbon is taken to be the currency to express plant and animal biomass.The reaeration coefficient (Ka) in the atmosphere-water surface exchange process is calculated using formula [5, 25]
(22)Ka=3.93U0H3/2+(0.728Uw1/2-0.317Uw+0.0372Uw2)H,
where U0 is the mean current velocity in m/sec, Uw is the wind speed measured 10 m above the water surface in m/sec, and H is the mean depth in meters.Saturating oxygen levels (Oss) for Chilka Lagoon are computed using formulae (21).Diatoms, MPB, and cyanophyceans are found to be dominating groups of primary producers in Chilka Lagoon during the period of study. It is well known that the diatoms and MPB utilise the available ammonia [5, 26, 27] in the water column for their growth during the photosynthetic activity and the amount of oxygen liberated in this process is proportional to their growth rates. Thus, the rate of oxygen liberation due to photosynthetic activity by diatoms and MPB can be given by [5]
(23)(32)(12)[GPhyPhy+GMPBMPB],
where GPhy and GMPB are the growth rates of diatoms and MPB, respectively. It is well known that cyanophyceans prefer nitrate and convert nitrate to ammonium, which is subsequently used for photosynthesis activity [5]. Thus, the rate of oxygen liberation due to photosynthesis activity by cyanophyceans is given by [5]
(24)[(32)(12)+(48)(14)anc]GCPhyCPhy,
where anc is nitrogen to carbon conversion factor and GCPhy is the growth rate of cyanophyceans.The nitrification process can be represented by the expression [5]
(25)(64)(14)K12θ12(WTemp-20)(OKNit+O)NH4,
where WTemp is the water temperature in °C, K12 is the nitrification rate KNit is the half-saturation constant for oxygen limitation during the nitrification activity, and NH4 is the ammonium concentration in the water column.The respiration process is the reverse process of photosynthesis and is modelled by multiplying the concentrations of phytoplankton and zooplankton by their respective respiration rates. Thus this process is mathematically represented as(26)(32)(12)[K1R1Phy+K1R2CPhy+K2RMPB+Z1RZoo],
where K1R1, K1R2, K2R, and Z1R are the respiration rate constants for diatoms, cyanophyceans, MPB, and zooplankton, respectively.The mineralisation process in aquatic system is governed by particulate organic matter mineralisation. During this process, one molecule ofO2 is consumed per micromole of particulate matter. The particulate matter depends on the mortality rates of phytoplankton and zooplankton. The mineralisation process is temperature dependent and this effect is modelled as Arrhenius type expressions [5, 17]. Thus the mineralisation process can be represented by the following expression:
(27)1(12)exp⁡KT(WTemp-20)[DPhyPhy+DMPBMPB+DCPhyCPhy+DZooZoo]O,
where KT is the mineralisation temperature correction coefficient. DPhy, DMPB, DCPhy, and DZoo are the mortality rates of diatoms, MPB, cyanophyceans, and zooplankton, respectively.Finally, the loss due to other processes is modelled as a linear function of available DO with a constant loss rate,KLoss. Since, quantification of this loss is difficult, calibration has been done to determine this rate constant.Thus the dynamics of DO in the water column of Chilka Lagoon can be represented by(28)dOdt=Ka(Oss-O)︸Atmospheric-WaterExchange-(6414)K12θ12(WTemp-20)(O(KNit+O))NH4︸Nitrification+(3212)[GPhyPhy+GMPBMPB]︸PhotosynthesisbyDiatomsandMPB+[(3212)+(4814)anc]GCPhyCPhy︸PhotosynthesisbyCyanophyceans-(3212)[K1R1Phy+K1R2CPhy+K2RMPB]︸Respirationbyplantcommunity-(3212)Z1RZoo︸RespirationbyZooplankton-(112)exp⁡KT(WTemp-20)[DPhyPhy+DMPBMPB+DCPhyCPhy+DZooZoo]O︸MineralisationbyPlantandAnimalCommunities-KLossO︸LossduetoOtherProcesses.The monthly average values of salinity, biomass of plant and animal communities, water temperature, and ammonia obtained from the periodical field surveys are interpolated with Hermite interpolation to obtain their daily estimates and are presented in Figure9. These estimates are used as inputs for the model. The description of the parameters and their units along with the assumed values are presented in the Table 3.Table 3
List of parameters.
Parameter (units)
Description
Range
Value
Source
Parameters
O (mgO2L-1)
Dissolved oxygen
Oss (mgO2L-1)
Dissolved oxygen saturation concentration
NH4 (mgNL-1)
Ammonium nitrogen concentration
Phy (mgCL-1)
Phytoplankton carbon concentration
Zoo (mgCL-1)
Zooplankton carbon concentration
WTemp (°C)
Water temperature
Rate parameters
Ka (day-1)
Reaeration rate at 20°C
Computed
K12 (day-1)
Nitrification rate constant at 20°C
0.005–0.54
0.54
[5]
K1R1 (day-1)
Diatoms respiration rate constant at 20°C
0.05–0.2
0.1
[5]
K2R (day-1)
MPB respiration rate constant at 20°C
0.05–0.2
0.1
[5]
K1R2 (day-1)
Cyanophyceans respiration rate constant at 20°C
0.05–0.2
0.147
[5]
Z1R (day-1)
Zooplankton respiration rate constant at 20°C
0.084
[15]
KLoss (Assumed) (day-1)
Loss rate due to other processes
1.86
—
Temp correction coefficients
θ12 (none)
Nitrification temperature correction coefficient
1.045–1.08
1.08
[5]
θ1R (none)
Respiration temperature correction coefficient
1.04–1.08
1.08
[5]
KT (none)
Mineralisation temperature correction coefficient
0.063
[5]
Half-saturation constants
KNIT (mgO2L-1)
Half-saturation constant of nitrification
1-2
2
[5]
Growth parameters
GPhy (day-1)
Specific growth rate constant of Diatoms
0.2–4
1.3
[5]
GMPB (day-1)
Specific growth rate constant of MPB
0.2–4
1.45
[5]
GCPhy (day-1)
Specific growth rate constant of Cyanophyceans
0.2–4
2.041
[5]
Mortality parameters
DPhy (day-1)
Mortality rate constant of Diatoms
0.005–0.2
0.1
[5]
DMPB (day-1)
Mortality rate constant of MPB
0.005–0.2
0.05
[5]
DCPhy (day-1)
Mortality rate constant of Cyanophyceans
0.005–0.2
0.145
[5]
DZoo (day-1)
Mortality rate constant of Zoo
0.01
[5]
Conversion factors
anc (mgNmg-1C)
Nitrogen to carbon ratio
0.15–0.25
0.1509
[15]Frames (a), (b), and (c) present the daily interpolated biomass of diatoms, cyanophyceans, and MPB in terms of carbon, respectively, for the study period of June 18, 2004 to September 22, 2006. Frames (d) and (e) present the concentrations of ammonium and nitrate in the water column. Frames (f), (g), and (h) present the variations in salinity, turbidity, and temperature in the lagoon.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)


### 7.1. Sensitivity Analysis

Sensitivity analysis attempts to provide a measure of the sensitivity of either parameters, or forcing functions to the state variables of greatest interest in the model [17, 28]. Sensitivity analysis is performed using the following formula:
(29)S=(∂x/x)(∂p/p),
where S is sensitivity, x is state variable (here O), and p is parameter. ∂x and ∂p are change of initial values of state variable, parameters, and forcing functions, respectively, at ±10% level and ±20% levels. Sensitivity analysis results are summarized in Table 4. Parameters, which are not possible to determine from field survey experiments are calibrated using first set data and further validated using second set data.Table 4
Sensitivity analysis of different parameters for DO (SDO) at ±10% and ±20% levels.
Parameter
Sensitivity
at ±10% level
at ±20% level
θ
12
−0.0004
−0.0011
G
Phy
+0.0358
+0.0716
G
CPhy
+0.0443
+0.0886
G
MPB
+0.0287
+0.0575
K
1
R
1
−0.0058
−0.0115
K
1
R
2
−0.0338
−0.0675
K
2
R
−0.0041
−0.0082
θ
1
R
−0.6133
−1.6854
Z
1
R
−0.0135
−0.0271
D
Phy
−0.0011
−0.0023
D
CPhy
−0.0074
−0.0148
D
MPB
−0.0004
−0.0008
D
Zoo
−0.0003
−0.0006
K
T
−0.0054
−0.0108
K
Loss
−0.0400
−0.0800


### 7.2.  Model Calibration and Validation

Model calibration is done by adjusting the selected parameters such as growth rates, loss rates in the model to obtain a best fit between the model calculations and the monthly average field data (Set number 1) collected during first year (June 18, 2004 to June 27, 2005). The monthly average field data (Set number 2) collected during second year (July 31, 2005 to 22, September 2006) is used to validate the model.


### 7.3.  Results

The simulated values of DO for Chilka are presented in Figure10. Clearly the results of simulation reflect the trends observed with the field data indicating that the chosen parameters are closer to reality.Figure 10
Figure representing the simulated values of DO along with the observed DO values of Chilka Lagoon along with the saturating oxygen concentrations for the entire period of study. The observations during the period of June 18, 2004 to June 27, 2005 (Set number 1) are used to calibrate the parameter values. The observations during the period July 31, 2005 to September 22, 2006 (Set number 2) are used to validate the proposed model.The goodness of fit for statistical significance between the field data and modelled results is performed using Chi-square test with the help of SPSS Software. The performance of the model is quantified from the RMSE value of the Chi-square test. The RMSE value for the model output is found to be 0.3028. Further, for the first and second years of data collection the Chi-square values are found to be 0.982 and 0.987, respectively, (P<0.05) and for the entire simulation period the Chi-square value is found to be 0.993 (P<0.05). The model validation results using regression analysis are presented in Figure 11.Figure 11
Figure representing the regression analysis plot with 95% confidence interval along with the regression statistics.Figure representing variations in the simulated values of the nondimensional parametersa (representative of nitrification), b (representative of net effect of mineralisation and loss due to other processes), c1 (representative of NEM), c (cumulative effect of NEM and diffusion from atmosphere), and a+b-c (cumulative effect of loss due to nitrification, mineralisation, and respiration and gain due to photosynthesis and diffusion from atmosphere, resp.) for Chilka Lagoon during entire period of study June 18, 2004 to September 22, 2006.
(a)
(b)
(c)
(d)
(e)The nondimensional parametersa, b, c1, and c for the model (8) have been computed for the period of simulation and are presented in Figure 12. From Figure 12(a) we observe that the parameter a which is representative of nitrification activity belongs to the range (7.7179×10-5, 3.1467×10-3) and the parameter b that represents the cumulative effect of mineralisation, surface exchange, and loss due to other processes ranges between 1.0815 and 2.7773 (cf. Figure 12(b)). From Figure 12(c) we observe that the parameter c1 which is representative of NEM varies between 0.20091 and 9.891 and is low during monsoon and high during recovery period (December–April) which coincides with the occurrence of massive cyanophyceans bloom. Figure 12(d) presents the variations in the parameter c which is the cumulative effect of NEM and atmospheric diffusion. The range for this parameter during the period of study is found to be (3.5685, 13.543). We also observe that the term a+b-c is negative all through the period of simulation ranging from −10.872 to −2.4444 (cf. Figure 12(e)) and it is inversely related to the oxygen present in the system (as predicted by the theory).


### 7.4. Influence of Individual Plankton Groups on the DO Dynamics

The biogeochemical model developed to simulate the DO trends in Chilka is used to assess the strength of each of the individual plankton groups in oxygenating the waters of Chilka. The expected scenarios for DO in Chilka Lagoon in the presence of only one of the plankton groups among diatoms, cyanophyceans and MPB are presented in Figure13.Simulation presenting the expected (modelled) DO trends in relation to saturating oxygen concentration if only one of the three dominant groups of algae observed in Chilka during the period of study June 18, 2004 to September 22, 2006. Frame (a) presents the situation when only Diatoms were present. Frame (b) presents the situation when only Cyanophyceans were present and the situation when only MPB were present is depicted in Frame (c).
(a)
(b)
(c)From Figure13 we clearly observe that the system would have been undersaturated all through the period if either diatoms or MPB alone were present in the system. Even in the case of cyanophyceans (Figure 9(b)) the system would be undersaturated during the simulation period except for a short duration (January 14 to April 14, 2005) during the massive bloom formation.Figure14 presents the expected DO trends if two out of the three plankton groups were present in the system. These numerical simulations clearly indicate that the supersaturation conditions observed are due to the presence of both diatoms and MPB (cf. Figure 13).Simulation presenting the expected (modelled) DO trends in relation to saturating oxygen concentration if one of the three dominant groups of algae are absent during the period of study June 18, 2004 to September 22, 2006. Frame (a) presents the situation when Diatoms were absent in the system. Frame (b) presents the situation when Cyanophyceans were absent from the system and the situation when MPB were absent from the system is depicted in Frame (c).
(a)
(b)
(c)Observe that the system cannot be driven to supersaturation in these periods in the absence of either diatoms or MPB (cf. Figures12 and 14). The pronounced supersaturation observed during January–March 2006 is due to the cyanophyceans bloom (cf. Figure 10). From Figures 10 and 13 it can be inferred that the near saturation state of DO after the massive cyanophyceans bloom is due to the presence of MPB and diatoms in the system else the system would have experienced hypoxic conditions during June 2006. Thus the health of Chilka Lagoon appears to be maintained by the presence of ammonium consuming plant community, that is, diatoms and MPB.


## 7.1. Sensitivity Analysis

Sensitivity analysis attempts to provide a measure of the sensitivity of either parameters, or forcing functions to the state variables of greatest interest in the model [17, 28]. Sensitivity analysis is performed using the following formula:
(29)S=(∂x/x)(∂p/p),
where S is sensitivity, x is state variable (here O), and p is parameter. ∂x and ∂p are change of initial values of state variable, parameters, and forcing functions, respectively, at ±10% level and ±20% levels. Sensitivity analysis results are summarized in Table 4. Parameters, which are not possible to determine from field survey experiments are calibrated using first set data and further validated using second set data.Table 4
Sensitivity analysis of different parameters for DO (SDO) at ±10% and ±20% levels.
Parameter
Sensitivity
at ±10% level
at ±20% level
θ
12
−0.0004
−0.0011
G
Phy
+0.0358
+0.0716
G
CPhy
+0.0443
+0.0886
G
MPB
+0.0287
+0.0575
K
1
R
1
−0.0058
−0.0115
K
1
R
2
−0.0338
−0.0675
K
2
R
−0.0041
−0.0082
θ
1
R
−0.6133
−1.6854
Z
1
R
−0.0135
−0.0271
D
Phy
−0.0011
−0.0023
D
CPhy
−0.0074
−0.0148
D
MPB
−0.0004
−0.0008
D
Zoo
−0.0003
−0.0006
K
T
−0.0054
−0.0108
K
Loss
−0.0400
−0.0800


## 7.2.  Model Calibration and Validation

Model calibration is done by adjusting the selected parameters such as growth rates, loss rates in the model to obtain a best fit between the model calculations and the monthly average field data (Set number 1) collected during first year (June 18, 2004 to June 27, 2005). The monthly average field data (Set number 2) collected during second year (July 31, 2005 to 22, September 2006) is used to validate the model.


## 7.3.  Results

The simulated values of DO for Chilka are presented in Figure10. Clearly the results of simulation reflect the trends observed with the field data indicating that the chosen parameters are closer to reality.Figure 10
Figure representing the simulated values of DO along with the observed DO values of Chilka Lagoon along with the saturating oxygen concentrations for the entire period of study. The observations during the period of June 18, 2004 to June 27, 2005 (Set number 1) are used to calibrate the parameter values. The observations during the period July 31, 2005 to September 22, 2006 (Set number 2) are used to validate the proposed model.The goodness of fit for statistical significance between the field data and modelled results is performed using Chi-square test with the help of SPSS Software. The performance of the model is quantified from the RMSE value of the Chi-square test. The RMSE value for the model output is found to be 0.3028. Further, for the first and second years of data collection the Chi-square values are found to be 0.982 and 0.987, respectively, (P<0.05) and for the entire simulation period the Chi-square value is found to be 0.993 (P<0.05). The model validation results using regression analysis are presented in Figure 11.Figure 11
Figure representing the regression analysis plot with 95% confidence interval along with the regression statistics.Figure representing variations in the simulated values of the nondimensional parametersa (representative of nitrification), b (representative of net effect of mineralisation and loss due to other processes), c1 (representative of NEM), c (cumulative effect of NEM and diffusion from atmosphere), and a+b-c (cumulative effect of loss due to nitrification, mineralisation, and respiration and gain due to photosynthesis and diffusion from atmosphere, resp.) for Chilka Lagoon during entire period of study June 18, 2004 to September 22, 2006.
(a)
(b)
(c)
(d)
(e)The nondimensional parametersa, b, c1, and c for the model (8) have been computed for the period of simulation and are presented in Figure 12. From Figure 12(a) we observe that the parameter a which is representative of nitrification activity belongs to the range (7.7179×10-5, 3.1467×10-3) and the parameter b that represents the cumulative effect of mineralisation, surface exchange, and loss due to other processes ranges between 1.0815 and 2.7773 (cf. Figure 12(b)). From Figure 12(c) we observe that the parameter c1 which is representative of NEM varies between 0.20091 and 9.891 and is low during monsoon and high during recovery period (December–April) which coincides with the occurrence of massive cyanophyceans bloom. Figure 12(d) presents the variations in the parameter c which is the cumulative effect of NEM and atmospheric diffusion. The range for this parameter during the period of study is found to be (3.5685, 13.543). We also observe that the term a+b-c is negative all through the period of simulation ranging from −10.872 to −2.4444 (cf. Figure 12(e)) and it is inversely related to the oxygen present in the system (as predicted by the theory).


## 7.4. Influence of Individual Plankton Groups on the DO Dynamics

The biogeochemical model developed to simulate the DO trends in Chilka is used to assess the strength of each of the individual plankton groups in oxygenating the waters of Chilka. The expected scenarios for DO in Chilka Lagoon in the presence of only one of the plankton groups among diatoms, cyanophyceans and MPB are presented in Figure13.Simulation presenting the expected (modelled) DO trends in relation to saturating oxygen concentration if only one of the three dominant groups of algae observed in Chilka during the period of study June 18, 2004 to September 22, 2006. Frame (a) presents the situation when only Diatoms were present. Frame (b) presents the situation when only Cyanophyceans were present and the situation when only MPB were present is depicted in Frame (c).
(a)
(b)
(c)From Figure13 we clearly observe that the system would have been undersaturated all through the period if either diatoms or MPB alone were present in the system. Even in the case of cyanophyceans (Figure 9(b)) the system would be undersaturated during the simulation period except for a short duration (January 14 to April 14, 2005) during the massive bloom formation.Figure14 presents the expected DO trends if two out of the three plankton groups were present in the system. These numerical simulations clearly indicate that the supersaturation conditions observed are due to the presence of both diatoms and MPB (cf. Figure 13).Simulation presenting the expected (modelled) DO trends in relation to saturating oxygen concentration if one of the three dominant groups of algae are absent during the period of study June 18, 2004 to September 22, 2006. Frame (a) presents the situation when Diatoms were absent in the system. Frame (b) presents the situation when Cyanophyceans were absent from the system and the situation when MPB were absent from the system is depicted in Frame (c).
(a)
(b)
(c)Observe that the system cannot be driven to supersaturation in these periods in the absence of either diatoms or MPB (cf. Figures12 and 14). The pronounced supersaturation observed during January–March 2006 is due to the cyanophyceans bloom (cf. Figure 10). From Figures 10 and 13 it can be inferred that the near saturation state of DO after the massive cyanophyceans bloom is due to the presence of MPB and diatoms in the system else the system would have experienced hypoxic conditions during June 2006. Thus the health of Chilka Lagoon appears to be maintained by the presence of ammonium consuming plant community, that is, diatoms and MPB.


## 8. Discussions and Conclusions

Dissolved oxygen (DO) is known to be a vital component in determining the health of aquatic systems. Much work was carried out on the events that are responsible for hypoxic/anoxic conditions in the aquatic systems [20, 29–33], characterised by differing topographical conditions, physical and chemical attributes. Excess nutrient inflows, nitrification, and denitrification events are the main causes for the observed conditions. A few key research concerns raised in this regard include “When does large-scale hypoxia occur?” “What amount of nutrient load reductions is needed to reduce the hypoxia?” [7, 34]. We show that mathematical models coupled with biology, physics, and chemistry can be used to address the above issues.In this paper a theoretical model of fairly general nature has been constructed to understand the dynamics of DO in the aquatic systems. Forcing factors for the model include atmosphere-water surface exchange, photosynthesis, respiration, nitrification, and mineralisation. In the proposed model the concentrations of plant and animal community in the water body are taken as parameters. This gives the advantage of studying the responses of the system with respect to variations in these concentrations and enhances the predictive potential of the model under various scenarios involving conditions on the concentrations of nutrient, plant, and animal community. The main interest of the study is to bridge the theoretical results obtained from the abstract model with that of simulation modelling studies for better understanding of the ecosystem dynamics. This type of approach has proved to be of great help in understanding system dynamics of Chilka Lagoon, India.The analysis of the theoretical model revealed several interesting phenomena that are responsible for decrease and increase in DO of aquatic systems. It is found that Net Ecosystem Metabolism (NEM), nitrification and/or mineralisation of detritus could play a key role in the eventual state of the system. The analysis indicates that the effect of mineralisation is stronger than the effect of nitrification in reducing the oxygen concentration in the system. For a system where in the processes of photosynthesis and surface exchange dominate the process of respiration, it is found that the saturation characteristics depend on the magnitude of nitrification and/or mineralisation rates. Supersaturation conditions prevail if the said rates are small enough. On the other hand, the system can get driven to hypoxic conditions at higher rates of nitrification and/or mineralisation. This observation is in conformity with studies done on hypoxia [13, 14].Further, the analysis reveals that an aquatic system can get driven to anoxic conditions even if respiration exceeds net contribution by photosynthesis and atmosphere-water surface exchange. In such a case the analysis indicates that the health of the system can only be improved (DO can be enhanced) by external supply of oxygen to water column as observed in the studies [35, 36]. Apart from the above findings, this study provides a few management strategies to drive the system away from hypoxic conditions towards higher DO concentrations. These strategies can be used to steer the system from undersaturation to supersaturation whenever NEM is positive. Such a transition is not possible for a system with negative NEM although the DO concentration can be enhanced.After a thorough analysis of the general model an effort is made to corroborate the theoretical findings with the field observations made for Chilka Lagoon, situated in East-Coast of India. Insights drawn from the theoretical analysis are used to explain the oscillations observed in DO for Chilka. From these investigations it is learnt that the macrophytes, submerged vegetation play a vital role in maintaining the health of the Chilka Lagoon. It is found out that during the active flow period (September–November) of 2005, the Chilka experienced undersaturation conditions due to NEM being less (in fact close to zero) and high nitrification activity. During the lean flow conditions (February–April) of 2005, Chilka experienced supersaturation conditions due to enhanced NEM (as a result of high photosynthetic activity of plant community) and reduced nitrification activity.Having learnt the importance of the plant species in Chilka Lagoon other than phytoplankton, a biogeochemical model has been developed to simulate the dynamics of DO in Chilka Lagoon. The construction of this model is done based on insights gained from the general model proposed for DO dynamics. Biogeochemical events taking place in the lagoon are taken into consideration when assigning forms to various coefficients in the model. Values for some of the involved parameters are drawn from the literature and a few have been calibrated. The results of the simulation are validated with that of the observed trends of DO in Chilka. Statistical tests such as Chi-square test for goodness of fit, regression analysis, and RMSE are conducted for establishing the ability of the proposed biogeochemical model to imitate the observed trends. From these simulations it has become possible to not only to estimate vital parameters present in the proposed model but also to capture the seasonal variations associated with these parameters. The proposed biogeochemical model has enabled evaluation of strength of each of the plant community in oxygenating water of Chilka. As a consequence, it is hypothesised that diatoms and MPB are mainly responsible for maintaining the DO levels of Chilka around saturation during the period of study.Finally, we wish to conclude that inclusion of macrophyte dynamics into the considered biogeochemical model would take the study much closer to reality. This calls for a model of higher dimension since the nutritional requirements of macrophytes are different from that of phytoplankton and MPB. Due to lack of data related to macrophytes for the case of Chilka, the simulation study is confined only to a model that involves explicit roles of phytoplankton and MPB. It is noteworthy that the considered model itself is able to replicate the system dynamics sufficiently close to reality enabling us to assess the roles of phytoplankton and MPB in the DO dynamics associated with Chilka Lagoon.



---

*Source: 526245-2014-02-20.xml*