A New Trophic State Index for Lagoons

1 Earth Sciences and Hydrology Division, Marine and Earth Sciences Group, Remote Sensing Applications and Image Processing Area, Space Applications Centre, Indian Space Research Organization, Ahmedabad, Gujarat 380 015, India 2 Centre for Earth Observation Science, Department of Environment and Geography, Clayton H. Riddell Faculty of Environment, Earth, and Resources, University of Manitoba, 463 Wallace Building, Winnipeg, MB, Canada R3T 2N2


Introduction
It is much easier to convey to the public the status and nature of a lagoon through an index, which provides the scientific accord of eutrophication and character of the lagoon.The monitoring of water quality, especially in the inland water bodies, aquifers, lakes, reservoirs, and lagoons, is important as it helps with the management of the eutrophication and productivity of the water body.Carlson [1] made early attempts to define an index, termed trophic state index (TSI), which could tell about the trophic status and nature of the lakes.Based on this trophic index, further classification of eutrophication of a lake into oligotrophic, mesotrophic, and eutrophic was possible.Because of water pollution and contamination of inland and coastal water bodies, it has become important to develop methods to determine such index of eutrophication of coastal lagoons that provides better spatial and temporal coverage.
Carlson [1] presented the trophic condition on a 0-100 scale; however, the author reported three different TSI values, that is, total phosphorus (TP), chlorophyll-a (Chl-a), and Secchi disk transparency depth (SDD) [2].To eliminate this confusion, Osgood [3] suggested using differences in the indices to assess the water quality in lakes.The basis of the development of an index of trophic state depends on the characteristic features and peculiarity of eutrophication.Different from Carlson's and others' approaches, O'Boyle et al. [4] have recently suggested pH-and dissolved oxygenbased index of status of eutrophication.Chalar et al. [5] report a completely different TSI based on benthic invertebrates, which is little dependent on environment variables.The variability and dominance of specific eutrophication parameters in different lakes hold the key to using a suitable index because different lakes may have different sources of nutrients, depth, and size.So far, there is no universally applicable trophic index reported which could adequately predict the productive character of a coastal lagoon.
However, many trophic scales, water quality indices, and measures of an assessment of water quality and trophic status of lakes, estuaries, and marine coastal areas are abundant in the literature.For example, Vollenweider et al. [6] proposed a new trophic index (TRIX) based on Chl-, oxygen saturation, mineral, total nitrogen, and phosphorus for coastal marine waters.TRIX is scaled from 0 to 10, which covers a wide range of trophic conditions.Analogously, a turbidity index (TRBIX) serves as a complimentary water quality index.TRIX and analogous other indices are most useful for saline coastal marine waters but are not well suited for a brackish water lagoon.Padisák et al. [7] have developed an assemblage index, , to assess ecological status of lake types established by the Water Framework Directive (WFD). varies from 0 to 5 and provides a 5-grade qualification as required by WFD.But  is primarily a phytoplankton assemblage index that neglects water quality parameters.OSPAR, Oslo/Paris convention (for the Protection of the Marine Environment of the North-East Atlantic), has adopted COMPP (Comprehensive Procedure) for the identification of eutrophication status of OSPAR Maritime Area [8].Although, OSPAR COMPP provides a holistic procedure for trophic state assessment, it does not set any thresholds for the trophic status indicator variables (e.g., Chl-a, water clarity, etc.), and there are no differential weightings across indicators.Bricker et al. [9] have suggested an analogous integrated methodology for the Assessment of Estuarine Trophic Status (ASSETS) for eutrophication status of estuaries and coastal areas [9].It considers a whole suite of nutrients and water quality parameters for a detailed assessment but suffers from lengthy and computation demanding procedure for trophic assessment.However, it has a limitation to only estuaries and coastal waters.An index that is easy to implement, a first-hand approach for a quicker assessment of trophic status of a brackish water lagoon, is still a requirement.The most common estimates of TSI have their basis on various in situ water quality data collected in the lake [1].The spatial and temporal changes that occur in the lake constrain these estimates.Wezernak et al. [10] reported a method using Chl-a and SDD as parameters for the determination of TSI.Thiemann and Kaufmann [11] determined TSI based on Chl-a estimated using field spectrometer data for lakes in Germany.Walker Jr. [12] suggested a TSI based on depletion of dissolved oxygen from in situ measurements as a measure of status of water quality of lakes.Analogous to Carlson's [1] trophic indices, Kratzer and Brezonik [13] estimated TSI for total nitrogen based on in situ measurements for various lakes in Florida, USA.Nitrite-nitrogen (NO 2 -N) is an important element in a lagoon system because it is present in plant tissues, compost, manure, soil, food, fertilizer, mineral, water, and yellow substances.Literature on a TSI primarily based on nitrogen (nitrite or nitrate or ammonia) is rare.This paper does not use ammonia to compute TSI as it nitrifies into nitrite by bacterial oxidation [4].This paper presents a new detailed model of TSI based on NO 2 -N with intercomparison to TSI based on Chl-a, and SDD for Chilika lagoon's eutrophication status.

Study Area and Sampling Data
The study area lies between 19 ∘ 28  -19 ∘ 54  N and 85 ∘ 05  -85 ∘ 38  E in the coastal regions of the state of Odisha, India (Figure 1).Eastern Ghats hills surround the eastern and western margins of the lagoon.There are 52 small rivulets in the north that drain fresh water into the lagoon, main rivers being Daya and Bhargavi Rivers, the distributaries of Mahanadi River [14].The Bay of Bengal seawater enters the lagoon through a mouth midway along the east coast length of the lagoon.
Chilika lagoon, situated on the east coast of India in the state of Odisha, is the second largest brackish water lagoon in the world [15].Ramsar Convention designated it as a Wetland of International Importance, since 1 October 1981.The pear-shaped lagoon sprawls about 64 km in length NE-SW and about 5-18 km in width from S-N covering an area of 1165-906 km 2 during the monsoon and summer, respectively.Chilika lagoon is a shallow water body with the deepest point of ∼3 m in the southern end and the shallowest point of ∼ 60 cm in the northern end.It opens in the Bay of Bengal through an artificially opened mouth approximately 200 m wide.This new artificial mouth was opened in 2007 by cutting and dredging the sand spit, to increase the eutrophication and protect wild life habitat, which were declining due to decreasing salinity after the closure of northern opening into the Bay of Bengal [14,16].
Water sample data collection campaign occurred during 2001-2003 at 27 different locations, evenly distributed in the lagoon to cover the entire lagoon area (Figure 1).Stations 15,15A,14,13,12,16,11,10,and 17 in the southern sector were comparatively deeper.Stations in the central sector were 1, 1A, 2, 2A, 8A, 9, 9A, 20A, and 21.The northern sector stations 3, 3A, 4, 7, 7A, 8, 19, 19A, and 20 were shallow, characterized by very high suspended sediment concentration.The sample analysis, that is, filtering, weighing, and so forth assisted in retrieving Chl-a, NO 2 -N, and SDD values.The samples were filtered through GF/C Whatman filter paper for the analysis of Chl-a.NO 2 -N data were analyzed according to the method described by Strickland and Parsons [17].Generally, Chilika lagoon has lower values of nitrite except during monsoon when there are higher nitrite values resulting from release by decomposed freshwater weeds into the water column [18].Table 1 provides the range and general statistics of Chl-a, NO 2 -N, and SDD.

Materials and Methods
The use of many variables responsible for the eutrophication of lagoon may question the selection of two-three variables for determining the trophic state of the lagoon.However, the evaluation of only a few variables is the basis of the need for a TSI that can provide the trophic status of the lagoon and to avoid the large computations involving multiparameter statistical analysis to arrive at an indicator of trophic state of the lagoon.A TSI should have its basis on only a few variables that one can easily measure and analyze with optimum effort.
In the present analysis, Chl-a, NO 2 -N, and SDD have been chosen as the potential candidates for estimating a TSI which can cover a substantial range of concentrations of each of the three parameters and is applicable to any brackish water lake, lagoon, inland water body, and reservoir.One can easily measure and understand the SDD, because it is dependent on direct human vision.SDD and Chl-a are representatives of the optical scattering and absorption in the lagoon waters.The underwater light availability in the lagoon controls the general growth/health of the microorganisms, biomass, flora, and fauna.Nitrogen is one of the major nutrient components of the food chain of a lagoon, rather to say, it lies at the base of food web and provides food for all living organisms in a lagoon.

Model A. This model computes a TSI based on the
Carlson's model [1].The relationship among three parameters, that is, Chl-a, NO 2 -N, and SDD, is nonlinear and is computed as follows (no transformation of the data values into logarithmic and so forth was performed, only the nontransformed actual data values have been used): where,  sample points constitute the data for regression;  ( The model underestimates the computed values of Chl-a and NO 2 -N in the absence of upper limits of the SDD.

Model B.
This model is a slight modification of Model A.
Defining the maximum limit of SDD and computing the TSI for Chl- and NO 2 -N from the empirical equations construct this model It is safe to consider SDD max as 200 cm for Chilika lagoon.Let us consider Solving ( 1) and ( 4), we get  = 1.5475 and  = 1.5650: 0.09 ≤ Chl- ≤ 20.00, 25 ≤ SDD ≤ 120.

Results and Discussion
A large amount of data collected over three-year span during monsoon, summer, and postmonsoon provides a basis of the results in this study.Thus, the data contain large variability required for formulating an ameliorated TSI for a brackish water lagoon.

Model A.
Carlson's index [1] provides a basis for Model A taking the entire data set into consideration.The computed TSI values are underestimated by this model and produce negative numbers for Chl-a and NO 2 -N. Figure 2 shows a plot of three TSI computed for Chl-a, NO 2 -N, and SDD using Model A. The variability in Chl-a and NO 2 -N does not reflect in the computed TSI whereas SDD is open-ended.Tables 1 and 2 provide the corresponding data values.Model A is obviously not a good TSI estimator for the Chilika lagoon.

Model B.
Here, the SDD has a limit of a maximum depth of 200 cm for Chilika lagoon.The choice of SDD max relies on the observational data in the lagoon.One can elucidate the physical significance of this limitation on SDD through the penetration of light into the lagoon waters.The light availability, in turn, controls the growth of organisms, hence, the eutrophication in the lagoon.There is a requirement of further modifications and regulations to arrive at a reasonable TSI, which consider the entire range of concentrations of Chla, NO 2 -N, and SDD from the surface to the SDD max in the lagoon.This model produces substantial negative values of the computed Chl-a and NO 2 -N; thus, it is not suitable for a good TSI (Table 3).Figure 3 shows the results of Model B. The computed TSI clearly elucidates the variability in Chl-a and NO 2 -N.However, Model B also underestimates NO 2 -N beyond 21 cm of SDD.This implies that the model considers only the surface concentration of NO 2 -N up to a depth of 21 cm (Table 3).(Figure 4).For instance, Figure 5 shows the observed TSI validated/computed using the Model C from the observed data from the Chilika lagoon [20].Each data point represents a mean value taken over several observations at select locations within the lagoon during monsoon, postmonsoon, and summer periods.Note that the data ranges provided in the Table 1 account for the entire data at all stations to create statistics.The observed data cover a large seasonal variability, for example, monsoon, summer, and postmonsoon, and temporal variability from 2001 to 2003.The computed TSI values go through a validation at each station in the lagoon.The physically observed eutrophication state and the predicted TSI-based health status are in good agreement.
A limited range of Chl-a varying between 0.09 and 20.00 mg m −3 and NO 2 -N between 0.41 and 1.07 g L −1 enables computation of reasonable values of TSI that covers greater ranges of concentration of Chl-a and NO 2 -N (Figure 4).The high values of Chl-a and NO 2 -N signify higher TSI in contrast to high values of transparency signifying low TSI.The variability in Chl-a and NO 2 -N is well accounted for by the computed TSI (Table 4).
The TSI provided in the Model C has many advantages over previous trophic indices.One can increase the scale of TSI (0 to 100) if the concentrations of observed Chl-a or NO 2 -N are even higher.The used ranges of Chl-a, NO 2 -N, and SDD provide sufficient information to construct a generalized TSI, which is likely to provide better estimates of the status of eutrophication in a lagoon.These parameters fundamentally account for the availability of organisms (Chl-a), food (nitrogen and phosphorus), and light (SDD), which are typical requirements for a healthy lagoon.In addition, these three parameters provide the biological (Chla), chemical (NO 2 -N), and physical (SDD) grounds for establishing the trophic state.Each of the three TSI obtained plays a role in biological, chemical, and physical contribution to the general eutrophication in the lagoon.One should consider the necessary abatement of controlling parameters for inverting the results to achieve the required TSI according to a standard.For example, Chilika Development Authority opened a new artificial mouth to the Bay of Bengal to increase the salinity of the lagoon and to control the eutrophication in the lagoon [16].addressed in previous studies.Yang et al. [2] reported TSI values ranging between 0 and 100 for Chl-a concentrations between 17.05 and 264.54 g L −1 and SDD between 12 and 257 cm compared to the ones in this paper (Table 1).Obviously, previous studies either lack the ability to cover a large range of Chl-a or do not consider a nitrogen-based TSI.

Spatial and Temporal Variability within the Lagoon.
Having discussed the three models, it is now imperative to advert to the spatial and temporal variability of Chla, NO 2 -N, and SDD within the lagoon.This is important for examination of computed TSI for different seasons and locations in the lagoon to assist in relevant management plans.The lagoon divides into three sectors spatially, that is, northern, southern, and central sectors, and temporally the lagoon experiences three different seasons during a year, that is, monsoon (July-October), postmonsoon (November-March), and summer (April-June).The mean Chl-a concentration is highest (8.33 mg m −3 ) in the central sector during summer.The mean NO 2 -N concentration is highest (15.98 g L −1 ) in the central sector during monsoon.The SDD is generally the lowest in the northern sector in any given season; that is, the water is clearer in the southern sector in all seasons.The proposed TSI accounts for the expected concentrations of Chl-a and NO 2 -N and the SDD for any season covering all sectors of the lagoon.

Conclusions
This study offers new Chl-a, NO 2 -N, and SDD based TSI expressions for the estimation of status of eutrophication in a lagoon system analogous to Carlson's index.Nitrogen being the fourth most abundant element found in the organisms is suitable for computing TSI of the lagoon.Nitrite is preferable because it decreases the fluorescence affecting the underwater photosynthesis and thereby controls primary productivity of Chilika lagoon.Nitrite is short-lived and readily converts into nitrate that does not affect fluorescence/photosynthesis [19].NO 2 -N is preferable to TP because of its close relationship with other trophic state criteria such as Chl-a and SDD.One can compute TSI using Chl-a, NO 2 -N, and SDD; each of which alone is capable of providing an estimate of trophic status of the lagoon.The three parameters account for the biological, chemical, and physical characteristics of the lagoon.It will be possible to estimate the TSI of freshwater and brackish water lagoons and other water bodies using the new expressions taking into consideration the spatial and temporal variability in the dataset.TSI (SDD) provides first-hand information about the trophic state of the lagoon.
Depending on the requirement and data availability, one can use alternative TSI (Chl-a) and TSI (NO 2 -N) to account for the biological and chemical contributions to eutrophication status.The estimated TSI can account for Chl-a and NO 2 -N up to 322.18 mg m −3 and 61.99 g L −1 , respectively.The TSI, based on these three parameters, serves as a complimentary index to assess nitrite levels and thereby primary productivity and as a predictive tool for lagoon management and prior to conducting field programs to monitor the health of a lagoon.

Figure 2 :
Figure 2: TSI versus corresponding values of Chl-a, NO 2 -N, and SDD computed using Model A. Most Chl-a and NO 2 -N values are underestimated which fall below zero.Note that the SDD is openended.

4. 3 .
Model C. The intersected area among the three curves corresponding to Chl-a, NO 2 -N, and SDD provides the most common TSI values in a typical lagoon

4. 4 .
Comparison with Other Indices.Several different indices, having their basis on various eutrophication parameters, are available in the literature since Carlson[1] (see Introduction).

Figure 3 :Figure 4 :Figure 5 :
Figure 3: TSI versus corresponding values of Chl-a, NO 2 -N, and SDD computed using Model B. Some NO 2 -N values are underestimated which fall below zero.The SDD has a limit of 200 cm.

Table 1 :
[19]range and general statistics of in situ data used in Model A computation for TSI.Reasons for Choosing Nitrite for TSI Estimation.NO 2 -N is preferable for TSI estimation because of its greater abundance in a trophic system and its relation to the other criteria of trophic state such as Chl-a and SDD.The presence of nitrite in biomass rich water decreases fluorescence thus affecting the Chl-a and photosynthesis[19], and thereby regulating primary production of the lagoon.A Chl-a based trophic index is likely to underestimate the trophic state of the lagoon in the presence of nitrite.A nitrite-based index can overcome this problem.Nitrate is not preferable, though it is more stable than nitrite, because nitrate does not have any influence on the Chl-a; thus, it does not affect the fluorescence and photosynthesis[19].Secondly, SDD-based TSI is likely to be confused with water quality index; however, the nutrientbased index should serve as a better indicator of the trophic status in a lagoon.Finally, the computed TSI values undergo a validation at each station in the lagoon to check if the computed TSI values are consistent with the observations in the lagoon.

Table 2 :
TSI and corresponding values of Chl-a, NO 2 -N, and SDD computed using Model A. It is noticeable that this model gives negative numbers (shown as "-" dash), which implies the nonapplicability of this model for Chl-a concentrations less than 0.04 mg m −3 and NO 2 -N concentrations less than 0.71 g L −1 .

Table 3 :
TSI and corresponding values of Chl-a, NO 2 -N, and SDD computed using Model B. It is noticeable that this model gives negative numbers (shown as "-" dash), which implies the nonapplicability of this model for NO 2 -N concentrations less than 0.13 g L −1 .Model B shows slight improvement over Model A; however, the values of computed indices are underestimated in this model.

Table 4 :
TSI and corresponding values of Chl-a, NO 2 -N, and SDD computed using Model C. It is noticeable that this model covers large ranges of NO 2 -N between 0.13 and 61.99 g L −1 and Chl-a between 1.00 and 322.18 mg m −3 .Model C shows significant improvement over Model A and Model B.