Simulating estuarine nitrous oxide production by means of a dynamic model

We describe a dynamic model developed from a commercially available modeling package (ECoS-III) to simulate estuarine dissolved inorganic nitrogen (DIN) dynamics, and consequent N(2)O production and atmospheric flux on the timescale of tidal cycles. Simulated model state variables were NH(4)(+), NO(3)(-) and N(2)O concentrations, and salinity. Model outputs were evaluated through comparison with summer field data for the Tyne estuary, UK. The model adequately reproduced the observed axial profiles of NH(4)(+), NO(3)(-) and N(2)O concentrations. Nitrification was shown to be the dominant N(2)O source and estimates of the ratios nitrification to DIN load and N(2)O emission to DIN load are considerably lower than the corresponding values adopted in global scale models of estuarine N(2)O emissions based on DIN transformations. Hence our results are consistent with the requirement imposed by atmospheric N(2)O growth rate constraints that the amount of atmospheric N(2)O arising from agriculturally related sources, including estuarine transformations of N, be revised downward.     

Spatial patterning of water quality in Biscayne Bay, Florida as a function of land use and water management

An objective classification analysis was performed on a water quality data set from 25 sites collected monthly during 1994-2003. The water quality parameters measured included: TN, TON, DIN, NH4+, NO3-, NO2-, TP, SRP, TN:TP ratio, TOC, DO, CHL A, turbidity, salinity and temperature. Based on this spatial analysis, Biscayne Bay was divided into five zones having similar water quality characteristics. A robust nutrient gradient, driven mostly by dissolved inorganic nitrogen, from alongshore to offshore in the main Bay, was a large determinant in the spatial clustering. Two of these zones (Alongshore and Inshore) were heavily influenced by freshwater input from four canals which drain the South Dade agricultural area, Black Point Landfill, and sewage treatment plant. The North Bay zone, with high turbidity, phytoplankton biomass, total phosphorus, and low DO, was affected by runoff from five canals, the Munisport Landfill, and the urban landscape. The South Bay zone, an embayment surrounded by mangrove wetlands with little urban development, was high in dissolved organic constituents but low in inorganic nutrients. The Main Bay was the area most influenced by water exchange with the Atlantic Ocean and showed the lowest nutrient concentrations. The water quality in Biscayne Bay is therefore highly dependent of the land use and influence from the watershed.     

A method to determine which nutrient is limiting for plant growth in estuarine waters--at any salinity

A method, utilising overlaid graphs for nutrients vs salinity, was developed in order to determine which nutrient is limiting for plant growth in estuarine waters-at any salinity. Dissolved inorganic nitrogen (DIN=NO(3)(-)+NO(2)(-)+NH(4)(+)) and o-phosphate (PO(4)(-)) are the main forms of N and P that are readily bio-available for plant growth in waters and these have a Redfield atomic ratio of N:P=16:1 (i.e. aquatic plants absorb N and P in the average ratio of 16 atoms of N to 1 atom of P). Graphs are prepared for (i) DIN vs salinity and (ii) o-phosphate vs salinity with the vertical scales for DIN and o-phosphate set at a ratio of N:P=16:1; when these graphs are overlaid on each other then the lowermost trendline denotes the limiting nutrient for plant/algal growth-at any salinity. The graphs also indicate the extent by which one or other of the nutrients is limiting--at any salinity. Furthermore, if there is a transition from P to N limitation somewhere along the salinity gradient, then this occurs at the salinity where the trendlines intersect. The concept was applied to three estuaries in the southeast of Ireland and the results show that, in all of these circumstances, P is the limiting nutrient throughout--except for the higher salinities (i.e. salinities 30 per thousand), where either (i) N and P may become equally limiting at salinity approximately 35 per thousand or (ii) N may become limiting at salinity 30 per thousand. Overlaid nutrients vs salinity graphs were also used to demonstrate that, in the estuaries in southeast Ireland, carbon (as dissolved inorganic carbon, DIC=CO(2)+H(2)CO(3)+HCO(3)(-)+CO(3)(2-)) is not the limiting nutrient--at any salinity.     

Simulation of submarine ground water discharge to a marine estuary: Biscayne Bay, Florida

Variable density ground water flow models are rarely used to estimate submarine ground water discharge because of limitations in computer speed, data availability, and availability of a simulation tool that can minimize numerical dispersion. This paper presents an application of the SEAWAT code, which is a combined version of MODFLOW and MT3D, to estimate rates of submarine ground water discharge to a coastal marine estuary. Discharge rates were estimated for Biscayne Bay, Florida, for the period from January 1989 to September 1998 using a three-dimensional, variable density ground water flow and transport model. Hydrologic stresses in the 10-layer model include recharge, evapotranspiration, ground water withdrawals from municipal wellfields, interactions with surface water (canals in urban areas and wetlands in the Everglades), boundary fluxes, and submarine ground water discharge to Biscayne Bay. The model was calibrated by matching ground water levels in monitoring wells, baseflow to canals, and the position of the 1995 salt water intrusion line. Results suggest that fresh submarine ground water discharge to Biscayne Bay may have exceeded surface water discharge during the 1989, 1990, and 1991 dry seasons, but the average discharge for the entire simulation period was only approximately 10% of the surface water discharge to the bay. Results from the model also suggest that tidal canals intercept fresh ground water that might otherwise have discharged directly to Biscayne Bay. This application demonstrates that regional scale variable density models are potentially useful tools for estimating rates of submarine ground water discharge.     

The characteristics of nutrients and eutrophication in the Pearl River estuary,South China

In the spring of 1998, 24-h time series and synchronization of vertical profiles of NO(3)-N, NO(2)-N, NH(3)-N, PO(4)-P, chlorophyll a, suspended substance, salinity, temperature and other chemical parameters were taken at 10 stations in the Pearl River estuary in order to analyze the status and characteristics of nutrients and eutrophication. The results indicated that dissolved inorganic nitrogen (DIN) mainly came from the four river channels in the main estuary, and NO(3)-N was the main form of DIN in most area. The concentration of DIN was general above 0.30 mg l(-1) in the estuary, and more than 0.50 mgl(-1) in most part. Phosphate from four river channels was not the main sources, but land-based sources from the area near Shenzhen Bay or along the estuary were obvious, and other land-based sources outside the estuary brought by coastal current and flood tide current were also the main contributions. The concentration of phosphate was generally about 0.015 mg l(-1) except the area near Shenzhen Bay. The ratio of N:P was generally high, and it was higher in the north than in the south. The highest ratio was higher than 300, and the lowest one was over 30. The concentration of chlorophyll a was about 0.8-7.8 mg m(-3), and turbidity and phosphate may be the main two limiting factors for algal bloom in the estuary. The concentration of nutrients decreased slightly in the past decade, but still stayed at a high level. The nutrients mainly came from domestic sewage, industrial wastewater, agriculture fertilizer and marine culture in the Pearl River estuary.     

A multivariate assessment of the coral ecosystem health of two embayments on the lee of the island of Hawai'i

Four coral-dominated coastal sites within two embayments (Kealakekua Bay and Honokōhau Bay) on the lee of the island of Hawai'i were studied to assess evidence of anthropogenic impacts in these relatively pristine locales. Nutrient-loading parameters were analyzed in relation to benthic composition data. Statistically, there were significant positive relationships between nitrate+nitrite, silicate, and ammonium with the abundance of macroalgae, coralline algae, and dead coral, and between delta(15)N and dead coral abundance. The north outside site of Kealakekua Bay and the south outside site of Honokōhau Bay appear to be most impacted by nutrient-loading factors in each bay, respectively. Comparisons with past nutrient data indicate that nutrient inputs have increased to the two bays, and that early impacts of these increased loadings are evident. It is predicted that at current nutrient-loading rates, the north sites of Kealakekua Bay and the south sites of Honokōhau Bay will exhibit evidence of further degradation in future years.     

Methane and nitrous oxide fluxes in the polluted Adyar River and estuary, SE India

We measured dissolved N(2)O, CH(4), O(2), NH(4)(+), NO(3)(-) and NO(2)(-) on 7 transects along the polluted Adyar River-estuary, SE India and estimated N(2)O and CH(4) emissions using a gas exchange relation and a floating chamber. High NO(2)(-) implied some nitrification of a large anthropogenic NH(4)(+) pool. In the lower catchment CH(4) was maximal (6.3+/-4.3 x 10(4)nM), exceeding the ebullition threshold, whereas strong undersaturation of N(2)O and O(2) implied intense denitrification. Emissions fluxes for the whole Adyar system approximately 2.5 x 10(8) g CH(4)yr(-1) and approximately 2.4 x 10(6)gN(2)O yr(-1) estimated with a gas exchange relation and approximately 2 x 10(9) g CH(4)yr(-1) derived with a floating chamber illustrate the importance of CH(4) ebullition. An equivalent CO(2) flux approximately 1-10 x 10(10)gy r(-1) derived using global warming potentials is equivalent to total Chennai motor vehicle CO(2) emissions in one month. Studies such as this may inform more effective waste management and future compliance with international emissions agreements.     

Occurrence of IRGAROL 1051 in coastal waters from Biscayne Bay,Florida, USA

Surface water samples from marinas, commercial ports and open bay areas collected from Biscayne Bay and the Miami River, Florida, USA, were analyzed for the occurrence of IRGAROL 1051 by GC/MS. The anifouling boosting herbicide was found in 80% (46/57) of the samples collected between March 1999 and September 2000. Concentrations within the bay range between non-detected (<1 ppt) and 61 ppt (ng/L) and were generally low compared with levels reported in European or Japanese waters. Aside from the elevated concentrations observed along the Miami River South Fork (61 ppt), the highest concentrations observed in the bay corresponded to marinas with high density of pleasure craft and restricted water circulation. In contrast, occurrence of IRGAROL 1051 along the commercial port or the cruise line terminal was generally lower (<1-2.2 ppt). Concentrations around Coconut Grove Marina were consistently higher (5-12 ppt) than the rest of the bay waters during the whole period of time surveyed.     

太湖梅梁湾藻类及相关环境因子逐步回归统计和蓝藻水华的初步预测

以太湖梅梁湾1992－1999年的连续监测资料为基础,运用多元逐步回归统计方法,选择水温等15项环境理化因素与藻类叶绿素a、藻类总生物量和微囊藻生物量等3项生物因素进行逐步回归分析,找出与生物因素显著相关的环境因子,建立多元逐步回归方程,预测梅梁湾藻类生物量的变化情况,初步进行了梅梁湾蓝藻水华的预测预报,结果显示,水温和总磷为梅梁湾藻类总生物量的显著相关因子,水温、硝态氮和总氮为微囊藻一物量的显著相关因子。     

The distribution and inventory of total and methyl mercury in Kuwait Bay

Sediment (surface and core) sampling was conducted during the period from August 96-December 98 from intertidal and offshore areas in Kuwait Bay. The highest (and most variable) concentrations of total mercury (T-Hg; 36500+/-34930 ng g(-1)) are encountered around the previous industrial outfall where sediments are disturbed by shipping activities. The concentrations of T-Hg are lower in the Shuwaikh Port area (650+/-210 ng g(-1)) and continue to decrease towards the northern coastline of Kuwait Bay (average concentrations in the wider Bay region are 50+/-30 ng g(-1)). These values are still above background concentrations of 15-20 ng g(-1). Calculation of the total inventory of mercury in all the surface sediments of the studied area indicates that approximately 22.5 ton is present which is similar to the estimated industrial discharges of approximately 20 ton suggesting that the contamination is largely confined to the Bay and that releases to the wider Gulf region are small. The distributions of MeHg are similar to those of T-Hg and represent ranges between 0.23% and 0.5% of the T-Hg indicating that surface sediments within Kuwait Bay contain approximately 80 kg of this chemical.     

Biotransformation of Organogenic Substances in the Aniva Bay Waters: Assessment with the Use of Mathematical Modeling

Mathematical modeling is used to study the conditions of natural biotransformation of biogenic element compounds (containing C, N, P, Si) in five areas within Aniva Bay. The input data for the model are evaluated with the use of GIS “Sakhalin Shelf” and the available reference literature on the bay oceanography. Water masses that transferred through the boundaries between the areas were evaluated based on the values of water level, river runoff, atmospheric precipitation, evaporation, and the volume of water in the bay, all of which vary as a result of water heating. The model was used to evaluate annual variations in the concentrations of organic and mineral fractions of these elements and the biomasses of planktonic organisms for the five areas in the bay. Modeling results also allowed the assessment of biogenic substance input with river runoff into Lososei Bay. Based on the internal fluxes of substances, the monthly and annual aquatic animal production was evaluated.     

Improved water quality in response to pollution control measures at Masan Bay, Korea

The total pollution load management system (TPLMS) was first applied in 2007 to the highly developed Masan Bay watershed, Korea. To evaluate the effect of TPLMS on water quality improvement, we analyzed the water qualities in rivers and bay during 2005-2010, targeting chemical oxygen demand (COD), suspended sediment (SS), total nitrogen (TN), and total phosphorus (TP) loads. Land-based pollutant loading all decreased during this period, with a significant reduction in COD and SS loads (p<0.01). The COD reduction in seawater, following the TPLMS implementation, was also significant (p<0.01). Time-lagged responses in COD and Chl-a supported an estimated seawater residence time of ~1 month. Land-based nutrient loads were also significantly reduced for TN (p<0.01) and TP (p<0.05), however, significant reductions were not observed in the bay, indicating potential alternative nutrient inputs from non-point sources into the bay system.     

Silicones in Chesapeake Bay sediments

Sediment samples from the upper reaches to the mouth of the Chesapeake Bay were analysed for extractable organic silicon (silicone) content. Organic silicon in the sediments ranged from less than 0.2 to some 36 ppm, w/w. In general, silicone tended to accumulate in high depositional areas: the northern bay, and along the shores to the south. While present, anthropogenic inputs in terms of defined source were not clear. It appears that processes in the bay act to rapidly remove silicone from the water column, probably on a function of the flocculation/precipitation of seston, especially in the northern area.     

Prediction of dissolved inorganic nitrogen (DIN) concentrations in deep,seasonally stratified lakes based on rates of DIN input and N removal processes

Mean dissolved inorganic nitrogen concentrations ([DIN]) in deep, seasonally stratified lakes with comparable DIN inputs can differ by up to a factor of 3 depending on hydraulic and morphometric properties and/or different trophic states of the lakes. In such lakes, net N sedimentation rates were estimated with two independent methods (sediment core analysis and input-output mass balances). They were higher in eutrophic lakes (Mean: 5.1; SD: ± 1.6 g m^–2 yr^–1; n = 13) than in oligotrophic lakes (1.6 ± 1.0 g m^–2 yr^–1; n = 3), but independent of [DIN]. Gaseous N loss rates to the atmosphere, as calculated from combined N- and P-mass balances from selected lakes, ranged from 0.9 to 37.4 g m^–2 yr^–1 (n = 10) and were positively correlated with [DIN]. Reduction of NO ₃ ^- to N₂ is assumed to be the main cause for gaseous N losses. A simple one-box mass balance model for [DIN], based on DIN input and rates and kinetics of N removal processes (net sedimentation and gaseous N loss) is proposed, and validated with a data base on [DIN] and DIN input in 19 deep, seasonally stratified lakes of central Europe. The model illustrated that the amount of water loading per unit surface area of a lake (called water discharge height q) is the critical parameter determining mean lake [DIN] relative to mean input [DIN]. Lakes with a q > 50 m yr^–1 have average [DIN] similar to the [DIN] of the inflows regardless of their trophic states, because input and outflow exceed lake-internal N removal processes. A high primary production favors DIN removal in lakes with q < 50 m yr^–1. It is concluded that measures to decrease primary production, e.g. by means of P removal programs, lead to an increase of [DIN] in lakes.     

Reservoirs are hotspots of nitrogen cycling in peatland catchments

This study presents input–output budgets of total dissolved nitrogen (TDN), dissolved organic N (DON) and dissolved inorganic N (DIN) for a reservoir in a peatland catchment in the south Pennines (UK). This site receives high levels of atmospheric inorganic N deposition, in the range of 26 kg N ha^?1 yr^?1. The results show that the reservoir retains ~21 to 31% of the annual TDN input (8806 ± 741 kg N). Approximately 39 to 55% of DON (3782 ± 653 kg N) and 6 to 13% of DIN (5024 ± 349 kg N) were retained/processed. A long water retention time (104 days), average annual pH of 6.5, high concentrations of DIN in the reservoir water and a deep water column suggest that denitrification is potentially a key mechanism of N retention/removal. The results also demonstrate that DON is potentially photodegraded and utilized within the reservoir, particularly during the summer season when 58 to 80% of DON input (682 ± 241 kg N) was retained, and a net export of DIN (~34 kg N) was observed. The findings therefore suggest that DON may play a more crucial role in the biogeochemistry of peat‐dominated acid sensitive upland freshwater systems than previously thought. Reservoirs, impoundments and large lakes in peatland catchments may be important sites in mediating downstream N transport and speciation. Copyright © 2016 John Wiley & Sons, Ltd.     

Estimating background and threshold nitrate concentrations using probability graphs

Because of the ubiquitous nature of anthropogenic nitrate (NO3(-)) in many parts of the world, determining background concentrations of NO3(-) in shallow ground water from natural sources is probably impossible in most environments. Present-day background must now include diffuse sources of NO3(-) such as disruption of soils and oxidation of organic matter, and atmospheric inputs from products of combustion and evaporation of ammonia from fertilizer and livestock waste. Anomalies can be defined as NO3(-) derived from nitrogen (N) inputs to the environment from anthropogenic activities, including synthetic fertilizers, livestock waste, and septic effluent. Cumulative probability graphs were used to identify threshold concentrations separating background and anomalous NO(3)-N concentrations and to assist in the determination of sources of N contamination for 232 spring water samples and 200 well water samples from karst aquifers. Thresholds were 0.4, 2.5, and 6.7 mg/L for spring water samples, and 0.1, 2.1, and 17 mg/L for well water samples. The 0.4 and 0.1 mg/L values are assumed to represent thresholds for present-day precipitation. Thresholds at 2.5 and 2.1 mg/L are interpreted to represent present-day background concentrations of NO(3)-N. The population of spring water samples with concentrations between 2.5 and 6.7 mg/L represents an amalgam of all sources of NO3(-) in the ground water basins that feed each spring; concentrations > 6.7 mg/L were typically samples collected soon after springtime application of synthetic fertilizer. The 17 mg/L threshold (adjusted to 15 mg/L) for well water samples is interpreted as the level above which livestock wastes dominate the N sources.     

Nitrogen loading to Pleasant Bay, Cape Cod: application of models and stable isotopes to detect incipient nutrient enrichment of estuaries

To test and refine methods to detect nutrient enrichment and resulting eutrophication, we applied the Waquoit Bay nitrogen loading model (NLM) and Estuarine loading model (ELM) to estuaries of Pleasant Bay that receive increasing but low N loads (25-199 kg N ha(-1) yr(-1)) from land. Contributions of wastewater to these estuaries increased from 7% to 63% as N loads increased, and modeled estimates of dissolved inorganic nitrogen in the water were within approximately 27% of measured values. N isotopic signatures in suspended and benthic organic matter and in tissue of quahogs increased as wastewater contributions to N loads increased, with clams approximately 4 per thousand heavier than organic matter, indicating that even at these low N loads, N from land-derived sources moved detectably up the food web. These results extend the application of NLM and ELM to detect incipient levels of N enrichment and demonstrate that these models can be used in conjunction with isotope measurements as the basis for food web analyses in a system exposed to relatively lower N loads than previously studied.     

Inorganic nitrogen control in wastewater treatment ponds from a fish farm (Orbetello, Italy): denitrification versus Ulva uptake

The aim of this study was to quantify the N removal efficiency of an Ulva-based phytotreatment system receiving wastewaters from a land-based fish farm (Orbetello, Italy), to identify the main biogeochemical pathways involved and to provide basic guidelines for treatment implementation and management. Fluxes of O2 and nutrients in bare and in Ulva colonised sediments were assessed by light/dark core incubations; denitrification by the isotope pairing technique and Ulva growth by in situ incubation of macroalgal disks in cages. O2 and nutrient budgets were estimated as sum of individual processes and further verified by 24-h investigations of overall inlet and outlet loads. Ulva uptake (up to 7.8 mmol Nm(-2) h(-1)) represented a net sink for water column and regenerated NH4+ whilst N removal via denitrification (10-170 micromol Nm(-2) h(-1)) accounted for a small percentage of inorganic nitrogen load (<5%). Laboratory experiments demonstrated a high potential for denitrification (over 800 microM Nm(-2) h(-1)) indicating that N loss could be enhanced. The control of Ulva standing stocks by optimised harvesting of surplus biomass may represent an effective strategy to maximise DIN removal and could result in the assimilation of approximately 50% of produced inorganic nitrogen.     

Modelling the relative impact of rivers (Scheldt/Rhine/Seine) and Western Channel waters on the nutrient and diatoms/Phaeocystis distributions in Belgian waters (Southern North Sea)

The coastal areas of the Southern North Sea (SNS) experience eutrophication problems resulting from freshwater nitrogen (N) and phosphorus (P) inputs from rivers. In particular, massive blooms of Phaeocystis colonies occur in Belgian waters. In this region, water masses result from the mixing of Western Channel (WCH) waters transported through the Straits of Dover with nutrient-rich freshwater from the Scheldt, the Rhine and Meuse, the Seine, the Thames and other smaller rivers. However, the relative contribution of the WCH and each river to the inorganic nutrient pool and the impact on the phytoplankton community structure (diatoms and Phaeocystis) are not known. In order to effectively manage the eutrophication problems, it is necessary to know: (i) the relative contribution of the WCH and of each river impacting the region and (ii) the relative effect of a N and/or P nutrient reduction on the Phaeocystis blooms. To answer these questions, sensitivity tests (1% nutrient reduction) and nutrient reduction scenarios (50% nutrient reduction) have been performed with a three-dimensional (3D) coupled physical–biogeochemical model (MIRO&CO-3D).MIRO&CO-3D results from the coupling of the COHERENS 3D hydrodynamic model with the ecological model MIRO. The model has been set up for the region between 48.5°N, 4°W and 52.5°N, 4.5°E and run to simulate the annual cycle of carbon, inorganic and organic nutrients, phytoplankton (diatoms and Phaeocystis), bacteria and zooplankton (microzooplankton and copepods) in the SNS under realistic forcing (meteorology and river inputs) for the period 1991–2003. The relative contribution of the WCH waters and of the different rivers on the inorganic nutrient pool available for phytoplankton (diatoms and Phaeocystis) growth is assessed by decreasing by 1% the nutrient (dissolved inorganic nitrogen, DIN and inorganic phosphate, PO₄) inputs from the WCH and from, respectively, the Scheldt (and smaller Belgian rivers), the Rhine/Meuse and the Seine (and smaller French rivers) [sensitivity tests]. The relative role of N and P reduction on the diatoms/Phaeocystis distribution is further explored by simulations with 50% reduction of the total (inorganic and organic) N and total P river inputs [nutrient reduction scenarios]. These scenarios allow assessing the impact of the expected 50% reduction of river nutrient inputs resulting from the implementation of nutrient reduction policy.Results of the sensitivity tests suggest that the impact of a 1% reduction of river nutrient inputs on surface nutrients (DIN and PO₄) over the Belgian Exclusive Economic Zone (EEZ) area is similar for the Seine and the Scheldt, which are in turn greater than for the Rhine. However, a hypothetical 1% reduction of nutrient input from the WCH boundary would have a higher impact than for the Scheldt. The impact of nutrient reduction is higher for DIN than for PO₄ whatever the river (contrary to the WCH). DIN is more sensitive to riverine nutrient reduction because the rivers are over enriched in DIN compared to PO₄. The sensitivity tests suggest also that a PO₄ river input reduction would result in a N:P increase and a DIN river input reduction would result in a N:P decrease but that a combined (PO₄ and DIN) input reduction would reduce the N:P ratio at sea.From 50% nutrient reduction scenarios, model results suggest that a total P reduction would induce a significant decrease of diatoms and a small (coast) to negligible (offshore) decrease of Phaeocystis biomass. On the contrary, a total N reduction would induce a significant decrease of Phaeocystis biomass and a moderate increase of diatoms. When N and P river input reductions are combined, the model predicts a significant decrease of Phaeocystis biomass in Belgian waters and a significant decrease of diatom biomass in the coastal waters and a small increase offshore. A future management plan aiming at Phaeocystis reduction should thus prioritise N reduction.     

Quantitative aspects of inorganic nutrient fluxes in the Gazi Bay (Kenya): implications for coastal ecosystems

Fluxes of dissolved inorganic nutrients: NH4+, NO2-, NO3-, PO4(3-) and Si(OH)4 from nearshore sediments of Gazi Bay were measured in situ within mangrove, seagrass and coral reef biotopes using benthic flux bell-jar chambers of cross-sectional area 0.066 m2 and volume 0.0132 m3. The objectives were: (1) to determine the influence of benthic fluxes, fluvial discharge and seasonal variations on the nutrient budget in the Bay waters; (2) to determine the effect of tidal and spatial variations on nutrient loads in the water column and (3) to establish the relative importance of the nutrient sources with regard to total community production of the Bay. The directly measured fluxes ranged from -270 to +148 micromol NH4+-N/m2/h; -60 to +63 micromol NO2(-)-N/m2/h; -79 to +41 micromol NO3(-)-N/m2/h; -79 to +75 micromol PO4(3-)-P/m2/h and +30 to +350 micromol Si(OH)4-Si/m2/h for and respectively. It was established that benthic fluxes are the major sources of dissolved inorganic NH4+, NO2- and Si(OH)4 while fluvial sources are important for NO3- and PO4(3-) into Gazi Bay waters. Seasonal variations had an appreciable effect on the PO4(3-) fluxes, N:Si ratio, river nutrient discharge, plankton productivity and important environmental factors such as salinity and temperature. Tidal and spatial variations had no significant effect on nutrient concentrations and net fluxes within the water column. The results imply that benthic fluxes are largely responsible for the nutrient dynamics of the nearshore coastal ecosystems especially where direct terrestrial inputs do not contribute significantly to the nutrient budget.