Nutrient budgets (C, N and P) and trophic dynamicsof a Brazilian tropical estuary: Barra das Jangadas

EARTH SCIENCES

Nutrient budgets (C, N and P) and trophic dynamicsof a Brazilian tropical estuary: Barra das Jangadas

Carlos E. D. Noriega; Moacyr Araujo

Laboratório de Oceanografia Física Estuarina e Costeira, Departamento de Oceanografia da Universidade Federal de Pernambuco, LOFEC/DOCEAN/UFPE, Avenida Arquitetura s/n, Cidade Universitária, 50740-550 Recife, PE, Brasil

^{Correspondence to}

ABSTRACT

This paper focuses on the nutrient dynamics of a tropical estuary on the northeastern Brazilian coast, studied using the LOICZ biogeochemical budgeting protocol. We describe the methodology and assumptions underlying this model. Input data (monthly for rainfall, evaporation, river discharge, and concentrations of salt, phosphorus and nitrogen) were obtained during field campaigns in the Barra das Jangadas Estuary (BJE) over a 5 years period (1999 to 2003). Mass balance results indicate large inputs of nutrients to the system. The model shows that the seasonal variation of the Net Ecosystem Metabolism (NEM) indicates that the system passes from a stage of organic matter liquid production and mineralization during the dry season (-0.5 mmoles C m^-2 d^-1) to liquid mineralization during the rainy season (-19 mmoles C m^-2 d^-1). We suggest that the system varies slightly between autotrophy and heterotrophy during the year due to the rainfall regime, human activities in the basin (density population and sugarcane plantations), and associated DIP riverine loads. High per capita loads of N and P indicate a high population density and high runoff. The application of flux balance modeling was useful to understand the nutrient dynamics of this typical small tropical estuary.

Key words: Brazil, Barra das Jangadas Estuary, seasonal variations, heterotrophy, autotrophy, carbon dioxide.

RESUMO

Este trabalho se focalizou na dinâmica de nutrientes de um estuário tropical na costa nordeste brasileira, usando o protocolo LOICZ de balanços biogeoquímicos. Nós descrevemos a metodologia e os pressupostos subjacentes a este modelo. Os dados de entrada (precipitação e evaporação mensal, vazão do rio, e as concentrações de sal, fósforo e nitrogênio) foram obtidos durante as campanhas de campo no estuário de Barra das Jangadas - Brasil durante um período de 5 anos (1999 a 2003). Os resultados indicam grandes entradas de nutrientes ao sistema. O modelo mostrou que a variação sazonal do Metabolismo do Ecosistema (NEM) indica que o sistema passa de uma fase de produção de líquido da matéria orgânica, durante a estação seca (-0,5 mmoles C m^-2 d^-1) para uma mineralização líquida durante a estação chuvosa (-19 mmol m^-2 C d^-1). Sugerimos que o sistema varia ligeiramente entre autotrófica e heterotrófica durante o ano, devido ao regime de chuvas, as atividades antrópicas na bacia (densidade populacional e as plantações de cana de açúcar), e as cargas ribeirinhas de DIP associadas. A alta carga de N e P per capita, indica uma alta densidade populacional e um alto runoff. A aplicação da modelagem de balanço de fluxos foi útil para o entendimento da dinâmica de nutrientes em um pequeno estuário tipicamente tropical.

Palavras-chave: Brasil, Barra das Jangadas Estuário, variações sazonais, heterotrofia, autotrofia, dióxido de carbono.

INTRODUCTION

Estuaries are dynamic systems, in which biological populations fluctuate according to natural cycles. Water quality also varies, particularly as seasonal and annual climatic patterns change. In these systems tracking environmental changes can be challenging, and distinguishing impacts caused by human actions from natural variations can be even more difficult (Marone et al. 2005). Under normal estuarine spatial and temporal constraints, reactive materials, such as nutrients, behave non-conservatively due to modifications by biological recycling and chemical transformations acting independently of simple physical advection and mixing (Dale and Prego 2005).

Furthermore, estuaries are areas in which anthropogenic effects, such as increased nutrient loads, have their most direct influence, and where there is a danger of adverse impacts. Most of these impacts results from a complex chain of events varying over different scales in space and time, which can be ultimately attributed to the accumulation of anthropogenic nitrogen and phosphorus in river water on its way to the ocean (Tappin 2002).

The magnitude of these fluxes is such that the transfer of organic matter from land to ocean via rivers is a key link in global carbon cycles. Due to the intense anthropogenic disturbance, estuaries are of ten considered to be net heterotrophic ecosystems and act as a source of CO₂ (Biswas et al. 2004, Mukhopadhyay et al. 2006). The increased nutrient load leads to eutrophication, enhances net ecosystem production, and shifts the system towards increased autotrophy (Gattuso et al. 1998). On the other hand, respiration of the organic carbon leads to increased heterotrophy. Additionally, light may become limiting for primary production in the upper part of estuaries (Irigoien and Castel 1997); respiration is then the dominant metabolic process, and an oxygen-depleted zone may occur, stimulating various anaerobic processes. It is well known that the chemical transformation pathways for nitrogen and phosphorus differ markedly from one another (Schlesinger 1997). In addition to being present in inorganic and organic dissolved forms, nitrogen is involved in biotic reactions and is the primary constituent of the atmosphere.

Besides direct uptake and release with respect to organic matter, the biotic processes of nitrogen fixation and denitrification actively move nitrogen between among the atmosphere (as nitrogen gas (N₂) and nitrous oxide (N₂O)) and both organic and inorganic forms of fixed nitrogen. Both nitrate (NO₃) and ammonia (NH₃) are highly soluble in water, and dissolved ammonia readily ionizes to ammonium (NH₄). Nitrate is an important byproduct of combustion, while ammonia is a highly volatile byproduct of animal waste. As a result, atmospheric transport and both wet and dry deposition are important pathways by which these materials are delivered to the landscape (Meyers et al. 2001). By contrast, phosphorus is involved in biotic reactions, primarily through the relatively simple (though still highly complex) pathways of organic production and oxidation. Phosphorus is also involved in various important mineral reactions (including both precipitation-dissolution of various forms of the mineral group apatite and adsorption- desorption reactions). In general, phosphorus is very particle-reactive and is taken up or released from the particles under changing conditions of pH, redox, and ionic strength. It has no significant gas phase.

The scatter in the loading ratio probably reflects, in large part, different chemical reaction pathways for DIN and DIP. The only real overlap in the reaction pathways for nitrogen and phosphorus involves production and oxidation of organic matter.

Because the composition ratio of nitrogen to phosphorus for most terrestrial organic matter is close to the DIN:DIP loading ratio we observed (approximately 19:1), decomposition of organic matter apparently dominates the inorganic nutrient loading, both in absolute range and loading ratio (Smith et al. 2003).

A close link is generally found between ecosystem metabolism and terrestrially derived nutrients in temperate ecosystems. It remains difficult to assess completely a function of estuarine ecosystem in response to the input of terrestrial nutrients in tropical area largely because of confounding physical and biogeochemical factors (Eyre and McKee 2002). Therefore, it is of interest to know whether a shallow coastal water body is a carbon source or sink, particularly in tropical areas Nutrient budgets can provide valuable information as to whether the system is a net exporter or importer of nutrients and can therefore determine its trophic status. Smith and Hollibaugh (1997) used the term "trophic status" to describe the net balance (net respiration or net synthesis) of organic carbon in an ecosystem. The results of these budgets and the use of stoichiometric tools provide estimates of processes such as net production/ respiration and nitrogen fixation/denitrification (Gordon et al. 1996). To assess carbon sources and sinks through process studies is not a simple task (Gordon et al. 1996). However, proposed guidelines for the Land- Ocean Interactions in the Coastal Zone (LOICZ) programme to assess non-conservative nutrient fluxes and carbon budgets for well boundary-defined coastal systems. This steady state budgeting method provides an alternative method to evaluate the biogeochemical metabolism and fate of nutrients and carbon in coastal systems when direct measurements of productivity and respiration are not available. Net nutrient fluxes in the coastal zone can be also determined from budget calculations, which is essential to evaluate the effects of riverine discharges on coastal function and carbon metabolism.

The rivers of the Northeast and East are marked by a pattern of seasonal flow typically unimodal, but differ in amplitude. As the climate states, the rivers of the Northeast are subjected to marked seasonal variability, with high intakes of pulses and floods during the wet season flows and low to negligible in the dry season (Knoppers et al. 2009).

In tropical ecosystems, mangrove-fringed estuaries play important roles in global processes, economic issues, political concerns and conservation strategies. Among numerous other processes, these tropical ecosystems affect the global carbon cycle (Lal et al. 2000).

Studies in tropical regions are of paramount importance for understanding the diversity of processes that occur at annual and seasonal scales and how these affect the biogeochemical cycles of the elements in these regions. Urbanization, industrialization, deforestation, agriculture, mining, and engineering works (i.e., dredging and damming) have changed the hydrological balance, material yields, and the water quality of estuarine systems, including those of the tropical Brazilian coast (Knoppers et al. 1999). This part of Brazil (Lat. 2ºS to 22ºS), harbors about 50 small and 7 medium-sized river estuaries subjected to either humid or semiarid climates (Ekau and Knoppers 1999), and includes the Barra das Jangadas.

The objective of the present work was to characterize and model the cultural eutrophication of the BJE to establish the mass balance of N and P throughout the year, considering detailed and complete (rainy and dry season) datasets from 1999 to 2003.

MATERIALS AND METHODS

STUDY AREA

The Barra das Jangadas Estuary (BJE) is the union of the lower course of the Pirapama and Jaboatão rivers, in the state of Pernambuco (Brazil). These basins cover semi-arid areas until they reach regions of intense urbanization along the coast, where they receive domestic and industrial effluents without previous treatment (CPRH 2003) (Noriega et al. 2009). Both rivers suffer the impact of domestic and agricultural effluents, mainly from the sugarcane agro-industry, under the form of high Biochemical Oxygen Demand (BOD), especially between November and March (dry season). During the rainy season, the higher freshwater discharge is the product of increased precipitation throughout the hydrographic basin (Araujo et al. 1999, Noriega et al. 2005a) and the controlled outflow from Pirapama Dam (Araujo et al. 2008).

The hydrographic basins of Jaboatão, Pirapama and other small rivers add up to 1000 km² of drainage area. BJE is a small estuarine area in which these rivers converge (8.7ºS-8.8ºS and 34.4ºW-34.8ºW). The estuary extends for approximately 13 km², with an average depth of 2.6 m (Branco 2002, Noriega et al. 2009) (Fig. 1).

The climate is typically tropical, hot and humid. The air temperature is 26±2.8ºC, and the mean annual precipitation and evaporation are around 1.5 and 1.2 m, respectively (Araujo et al. 1999). The rainfall regime is subdivided into two well-defined periods: the dry season (September-February), when the precipitation is exceeded by evaporation; and the rainy season (March- August), when rainfall dominates evaporation (Fig. 2).

The drainage basin includes areas originally covered by the Atlantic Rain Forest, and is presently occupied by sugar-cane and high density populated areas (1100 inhabitants km^-2) (IBGE 2000). Despite the deforestation of the margins and the large volume of industrial and domestic effluents received, the estuary itself is surrounded by relatively well-preserved and highly productive mangrove forests. Organic matter pollution by the sugar-cane agroindustry substantially increases during the harvest and milling season, which is from September to February. CPRH (2003) reported high BOD in the harvest periods of 69.6 mg L-1 (Jaboatão River) and 152 mg L-1 (Pirapama River). The polluting organic load sources are represented mainly by domestic sewage in the Jaboatão river (14.46 t BOD d^-1) and by agro-industrial activities in the Pirapama river (24.13 t BOD d^-1) (CPRH 2003). Algal blooms are now more frequent during the year and consist of several species of Cyanophyceae, mainly Microcystis aeruginosa, Oscillatoria sp and Euglena sp (Euglenophyta), suggesting some degree of permanent impact on the environment (Branco 2002).

The river runoff is strongly controlled by rainfall (Fig. 2), with an average discharge of 15 m³ s^-1 (annual average) (SECTMA 1999). The tidal regime is semidiurnal, with a mean amplitude of 1.3 m (neap tides) and 1.8 m (spring tides) (Araujo et al. 1999). The estuary is well mixed, being classified as type 1 with an absence of vertical stratification (Araujo et al. 1999, Noriega et al. 2009).

SAMPLING AND COMPILATION OF EARLY DATA

The nutrient and salinity data used in this study are monthly, annual and seasonal average concentrations calculated from 1999 to 2003 for the BJE and Jaboatão and Pirapama rivers (CPRH 2003, Branco 2002, 2006, Noriega et al. 2005a, b, 2009). Coastal concentrations were obtained from BNDO (2004). These data were used to construct an annual nutrient budget using the LOICZ approach, as proposed by Gordon et al. (1996). The construction of balances (annual, wet and dry period) was made through the average of the monthly balance sheets, following the methodology proposed by Webster et al. (2000).

Data for river runoff were obtained from SECTMA (1999). In order to obtain monthly estimates for the years 2000-2003, a Schreiber's model modified by Holland (1978) was applied using the measured monthly precipitation and air temperature of the watershed, to calculate surface runoff, by calculating the differences between precipitation and evaporation over a drainage basin valid for tropical and temperate regions (Gordon et al. 1996). Meteorological data from 1999- 2003 were obtained from INMET (INMET 2003).

WATER, SALT, AND FLUX CALCULATIONS FOR DIN AND DIP: THE LOICZ MODEL

The "Land Ocean Interactions in the Coastal Zone (LOICZ) Core Project" of the IGBP, established in 1993, is dedicated to understand the role of coastal subsystems in the functioning of the world oceans, including the role of the coastal zones and in the disturbed and undisturbed cycles of carbon, nitrogen and phosphorus (Gordon et al. 1996). The advantage of the LOICZ model is that extensive datasets are not required, so it is a suitable model for the Brazilian Northeast, where water quality data for most estuaries are extremely limited. The model is considered robust and uses a widely applicable, uniform methodology to provide information on the CNP fluxes in estuaries. Within the context of LOICZ biogeochemical modeling, the primary question to be addressed concerns the role of the coastal zone as a source or sink for carbon, nitrogen, and phosphorus (Wepener 2007).

The LOICZ biogeochemical model is based on the mass balance of water and materials (Gordon et al. 1996, Smith et al. 2005). Water and salt are assumed to not undergo significant biogeochemical transformations within the system, while nutrients behave as nonconservative compounds due to biogeochemical processing within the system. Hence, salt budgets and known water inputs and outputs are used to estimate water exchange between the system and the adjacent sea. The mass balance of essential non-conservative nutrients, namely dissolved inorganic phosphorus (DIP) and nitrogen (DIN), allows estimates to be made of rates of biological transformations and ecosystem processes, such as the net ecosystem metabolic (NEM) - i.e., the difference between primary production and community respiration - and the net nitrogen budget, which is assumed to depend on the difference between the nitrogen fixation and denitrification rates.

The water budget can be easily estimated using measurements of runoff (V_R), precipitation (V_P), groundwater (V_G), sewage or other inputs (V_O), and evaporative outflow (V_E). The compensating outflow or inflow that balances the water volume in the system is called the residual flow (V_R). The seawater volume necessary to maintain the salinity in the system (mixing flow, V_X) can be estimated using the conservative salt budget. The salt budget is calculated using the salinity difference between the system and the adjacent sea. DIP and DIN budgets are calculated from water budgets and concentration data. Deviations of budgets/concentrations (ΔDIP and ΔDIN) from predicted values are assumed to depend on non-conservative processes or internal transformations, and basically represent the net difference between nutrient sources and sinks.

In the LOICZ model, both NEM and the net nitrogen budget are calculated from ΔDIP and the molar C:N:P ratios of the reacting organic matter, generally that of the dominant primary producers, but other material (e.g., sewage) may be considered if judged to be significant. This assumes that ΔDIP depends only on biological transformations.

The molar C:N:P ratios can be considered as the link among the cycles of these elements in the production and respiration processes, and deviations from the expectations based on these ratios can be quantitatively assigned to other processes. The net nitrogen budget is the difference between ΔDIN and ΔDIN expected from ΔDIP and the C:N:P ratios, and is considered to be the difference between nitrogen fixation and denitrification (nfix-denit). Care is required in interpreting ΔDIP because it is affected by benthic fluxes and sorptive processes with suspended materials, as well as biotic processes in the water column. These effects may be especially important in shallow and turbid water bodies. Moreover, the C:N:P stoichiometry and its effects on ecosystem processes vary greatly among primary producers (Hessen et al. 2004).

Linear regression was used to observe the correlation between the model results and variables associated with these calculations. All the analyses were produced with the statistical software STATISTICA 8.0 for Windows.

RESULTS AND DISCUSSION

WATER AND SALT BALANCES

Water flow, salinity, and nutrient concentrations for the BJE, based on the monthly averages of samples taken from January 1999 to December 2003 (Table I), were successfully fed into the model. The residual water flux (V_R, Gordon et al. 1996) from this system, which is necessary to balance the freshwater outflow, was approximately 2 times greater in the rainy season than in the dry season. V_R occurs as a result of river runoff, precipitation influx and evaporation outflow. The amount of freshwater flowing (V_Q) into the estuary was estimated at 638 × 10³ and 1366 × 10³ m³, respectively (Fig. 3). Rainfall in the study area was seasonal. About 75% of rainfall occurred during the rainy season. The rain volume over the whole estuary area (13 km²) was found to be 41 and 125 × 10³ m³ d^-1, respectively, with an annual value of 83 × 10³ m³ d^-1. Evaporation from the water surface was calculated to be 1.4 times greater in the dry season than in the rainy season (Fig. 3).

The residence time of water was estimated by dividing the estuary volume by the sum of the mixing exchange flux (V_X) and the residual flux (V_R). Therefore, the time of total water exchange in the BJE ranged between 7-21 days. As expected, higher residence times were observed during the dry season (>13 days), while during the rainy months they did not exceed 9 days (Table II).

DIP AND DIN BUDGET

The calculated four major components of the material balance (V_RDIP_R, V_X (DIP_OC-DIP_SYS), V_GDIP_G, V_QDIP_Q), in the estuary indicated that BJE acts as a source for DIP (import-export = 1220 mol d^-1) (Fig. 4). Throughout the year, phosphorus (P) budgets in the BJE change sign according to the season. During the dry season, the balances are negative during four months (Table II). This implies that the BJE sequesters the difference between import and export, and acts as a sink for P during these months. This period of the year shows the smallest potential of exportation and higher residence times. The mixing outflow of DIP from this system is substantially larger than the residual inflow, and demonstrates that there must be DIP production (ΔDIP) of approximately + 57 mol d^-1 in the system. We assume that this represents decomposition of organic matter. We have observed that there is very high release of DIP, especially from the sediments associated with sugar cane wastes, so this and other organic discharges into the system are assumed to support the high non-conservative flux of DIP.

This period (September-March) is when the sugarcane harvest and milling occur. During the rainy season, the sign is positive. This implies that the BJE acts as a source for P, mainly in July, which is a characteristically high rainfall month (Table II).

Seasonal variations of river runoff (V_QDIPQ) DIP were found to be 2396 and 5782 moles d^-1 during the dry and rainy seasons, respectively. SECTMA (1999) indicated a residual organic pollution load during the sugar-cane harvest of 5000 and 2000 kg of BOD d^-1 for the Jaboatão and Pirapama rivers, respectively, which represents 57.8 and 37.8 kg of DIP d^-1, according to the coefficients of San Diego-McGlone et al. (2000). The total estimated for this study was 95.6 kg of DIP d^-1, while our estimates were 74.2 kg of DIP d^-1 for the period of the sugar-cane harvest (dry season).

According to CPRH (2003), the Jaboatão River showed concentrations of total phosphorus 1.3 times higher than the Pirapama River during the dry season (0.71 / 0.56 mg L-1) from 1999 to 2003, and 1.5 times higher during the rainy season. The Jabotão River had a P total 1.3 times higher than the Pirapama River in the dry season, and 1.5 times higher in the rainy season (5 years monitoring) (CPRH 2003).

The daily DIP load from the watershed was 4089 moles d^-1 (annual average) (Table II). This represents a riverine load in the BJE of 1492 moles of DIP per km^-2 per yr^-1 (or 46 kg of DIP per km^-2 per yr^-1) ((Load DIP_Q/Watershed area) × 365).

During the dry season, the flux per unit area of catchment was 874 moles of DIP per km^-2 per yr^-1, and during rainy season it was estimated at 2110 moles of DIP per km^-2 per yr^-1. According to Smith et al. (2003), the average concentrations of 6 mmoles m-3 correspond to a high population density (1000 people per km^-2) and high runoff (V_Q) per unit area (1 m yr^-1) in excess of 6300 moles per km^-2 per yr^-1. These authors used a regression model to describe DIN and DIP exportation by analyzing 165 systems for which DIN and DIP flux data were available (http://data.ecology. su.se/MNODE/). In the present study, the June load reached 3154 moles per km^-2 per yr^-1, a value 0.5 times lower than that estimated by Smith et al. (2003). The region presents a per capita load for the hydrographic basin of 1.4 moles of DIP per person per yr^-1 ((Load DIPQ/Population Watershed) × 365) or 0.04 kg per person per yr^-1 of DIP (annual average), a value that reflects the high population density and low runoff, according to Smith et al. (2003). The population density of the Jaboatão and Pirapama basins is about 1100 people per km^-2, with a total population of 1 100 000 hab. Bidone and Lacerda (2002) estimated a daily riverine load for estuarine or riverine areas in Northeast Brazil around 0.002 kg of P hab-1 d^-1. This value was computed taking into account 200-250 L hab-ΔD^-1 as the typical water consumption for developing countries. This is actually overestimated compared to the real water consumption of 100-150 L hab-ΔD^-1 in the Brazilian northeastern coastal area (I.M. Abreu et al., unpublished data). These figures would yield 0.01 kg of DIP per person per yr^-1, considering the San Diego-McGlone et al. (2000) transformations from P to DIP, which is 4-fold smaller than our previously calculated value. During the dry season, this value was nearer to that estimated by Bidone and Lacerda (2002) (0.02 kg of DIP per person per yr^-1). Therefore, during the rainy season, it was 6 times higher. Our estimates identify a high organic load as a result of unplanned activities in the hydrographic basin, such as disposal of domestic sewage, uncontrolled land runoff, and industrial and agro-industrial effluents.

The calculated four major components of the material balance (V_RDIN_R, V_X (DIN_OC-DIN_SYS), V_GDIN_G, V_QDIN_Q), in the estuary indicated that BJE acts as a sink for DIN (import-export = -46 046 mol d^-1) (Fig. 5).

The nitrogen budget in the BJE shows a higher inflow contribution during the rainy season, due to freshwater runoff and reduced residual water flow (Fig. 5). ΔDIN is negative in all months, suggesting that inputs are higher than outputs and indicating a net sink of inorganic nitrogen in the BJE. Seasonal variations of river runoff DIN (V_QDIN_Q) were found to be 46 061 and 106 978 moles d^-1 during the dry and rainy seasons, respectively. According to SECTMA (1999), the daily loads of BOD during the sugar-cane harvest (dry season) for the Jaboatão and Pirapama rivers are 5000 and 2000 kg of BOD d^-1, respectively, which represents 950 and 380 kg of DIN d^-1, according to the coefficients of San Diego-McGlone et al. (2000). The total daily load estimated for this study was 1092 kg of DIN d^-1, while our estimates were 687 kg of DIN d^-1 for the period of the sugar-cane harvest (dry season).

According to CPRH (2003), the Jaboatão River shows concentrations of ammonium 4 times higher than the Pirapama River (3.37/0.82 (dry season); 3.0/ 0.70 mg L-1 (rainy season)). The daily DIN load from the watershed was 78 019 moles d^-1 (annual average) (Table II), which would represent a yield of 28 477 moles of DIN per km^-2 per yr^-1 (400 kg of DIN per km^-2 per yr^-1), a value two times greater than the mean estimate (Smith et al. 2003).

During the dry season, the flux per unit area of the catchment was 17 907 moles of DIN per km^-2 per yr^-1, while during the rainy season it was estimated as 39 047 moles of DIN per km^-2 per yr^-1. The computed per capita load was 26 moles of DIN per person per yr^-1 (0.4 kg of DIN per person per yr^-1), a value 9 times greater than the scenario with high density and low runoff of Smith et al. (2003). According to Bidone (2000), the nitrogen (N) load for the regional scenario is 0.01 kg of N hab-ΔD^-1, equivalent to 0.004 kg of DIN hab-1 d^-1.

STOICHIOMETRIC CALCULATIONS OF NET SYSTEM METABOLISM

The evolution of the BJE metabolism shows a tendency towards heterotrophy (Fig. 6). The negative net ecosystem metabolism ( p-r) values indicate that the system is heterotrophic, with a net loss of organic matter from the BJE of -10 mmoles C m^-2 d^-1 (annual average). We believe that these values can vary if we consider other rates (42:12:1 SanDiego-McGlone et al. (2000) for the waste load ratio), or with mangroves dominating the net production (rate of 1000:11:1 Smith and Camacho (2000)). In either cases, if the DIP uptake primarily represents net organic metabolism, rather than sorption or precipitation of inorganic P, this system is net heterotrophic. During the dry season (September- February), we observed a slightly heterotrophy of -0.5 mmoles C m^-2 d^-1 (seasonal average). The long residence time retains materials long enough to react internally during the dry season. In addition, the water quality is enhanced by a slightly deeper euphotic depth (Noriega et al. 2005, Branco 2002), which also favors the phytoplankton community. In the rainy season, the biggest nutrient contribution to the rivers occurs (-19 mmoles C m^-2 d^-1; seasonal average). Mukhopadhyay et al. (2006) suggest that tropical estuaries with a shallow photic zone dominated by physical processes could cause the phytoplankton to not reach their maximum growth rates, which could contribute to the phytoplanktonic production of the estuary being limited. The objective is to modify the riverine flux of nutrients before it is released to the coastal water. These values demonstrate that outputs are higher than inputs at the BJE, with highlights on the mineralization of organic matter and a net source of CO₂ to the atmosphere.

The seasonal differences between heterotrophy (January to August) and autotrophy (September to December) indicate an extension of this second condition (Fig. 6).

Gordon et al. (1996) pointed out that p (primary production) and r (respiration) are within about 10% of one another. Assuming that p is known, this implies that the quantity ( p-r) = ±0.1 p. The lack of direct measures of primary productivity in the studied area was approached through the following: (i) mean annual and seasonal values of regional systems with biological characteristics similar to the phytoplanktonic biomass and species taxa (Passavante and Feitosa 2004), and (ii) studies in the literature about primary productivity for tropical systems (Berger 1989). The regional productivity is 128 mmoles C m^-2 d^-1 (annual mean), 101 mmoles C m^-2 d^-1 (rainy mean) and 155 mmoles C m^-2 d^-1 (dry mean) (Fig. 6). The literature reports values ranging from 15 to 399 mmoles C m^-2 d^-1. So, the estimates from regional averages represent an appropriate value to validate the results from the present study. The (p-r) estimate of -10 mmoles C m^-2 d^-1 (annual mean) represents 10% of the primary production, which is considered appropriate. In this way, respiration would represent -138 mmoles C m^-2 d^-1, and p/r = 0.92, which means that the system uses 8% more organic matter than it produces. During the dry season, the value of (p-r) is -0.5 mmoles C m^-2 d^-1, considering that phytoplanktonic primary production is 155 mmoles C m^-2 d^-1, p/r = 0.99. It follows that r is approximately 155.5 mmoles C m^-2 d^-1. That is, the system produces about 0.1 more organic matter than it uses. On the other hand, during the rainy season, (p-r) = -19 mmoles C m^-2 d^-1 and primary production is 101 mmoles C m^-2 d^-1, and so the value of r is 120 mmoles C m^-2 d^-1, with p/r = 0.84. The system consumes about 16% more organic matter than it produces in this second case.

A simple linear regression was used to relate pr values to the residence time (TR) (P<0.05) (Fig. 7). Longer residence times indicate that the system remains closer to 0, with a small trend towards autotrophy. On the other hand, shorter residence times show oscillations between heterotrophy and autotrophy (Fig. 7).

Rainfall often favors heterotrophic aquatic metabolism due to the increase in the contribution of terrestrial organic lixiviation (Ram et al. 2003). However, rainfall intensification also increases nitrogen and phosphorus loads in estuaries (Schindler 1978), which would benefit autotrophic metabolism, especially in urban and agricultural areas. During periods with opposite rainfall characteristics, metabolism seems to oscillate between light autotrophy and light heterotrophy.

The nitrogen fixation and denitrification are important processes in coastal systems. Again, because the major source of reacting matter is unclear, two N/P ratios are used. The decomposing material has a mean C/P of 106/1, and N/P of 16/1, which is near the value of N/P of 11/1 quoted for mangrove litter (Gordon et al. 1996). Based on this ratio of N/P, we estimated that ΔDINobs ΔDINesp ΔDIP* 11 = -5 mmoles m^-2 d^-1 (annual average). Smith and Camacho (2000) estimated that the differences between N fixation and denitrification are in general close to zero (with a dominance of denitrification), and that values above 5 moles m^-2 yr^-1 are rare. Our results in general suggest denitrification (Fig. 8).

The nitrogen fixation process is ordinarily slow in marine systems (<1 mmoles m^-2 d^-1), according to Swaney and Smith (2003), although they suggested that some coral reef, mangrove and tropical seagrass communities may exhibit rates >20 times this upper limit. As a general rule, few systems have nitrogen fixation faster than this rate. The value reported for the BJE in Februarywas low, submitting to this limit, and indicating that the adjacent mangrove forest did not accelerate this fixation in the estuary.

The apparently high denitrification during the rainy season (-7 mmoles m^-2 d^-1) indicates high benthic respiration (driven by high loads with labile organic matter such as sewage). Typical rates in benthic systems are around 0.5-2 mmoles N m^-2 d^-1. Systems with high benthic respiration may have denitrification rates >10 mmoles m^-2 d^-1 (Swaney and Smith 2003). During the dry season, denitrification is lower (-2 mmoles m^-2 d^-1) than in the rainy period. Other tropical estuaries, such as the PiauíRiver Estuary (Brazil), presented a denitrification rate of -0.13 mmoles m^-2 d^-1, while the Sergipe River Estuary (Brazil) seems to fix nitrogen at 0.1 mmoles m^-2 d^-1 (Souza 2000).

CONCLUSIONS

We used a bulk modeling approach to evaluate the nutrient budgets (C, N and P) and the trophic state of a tropical estuarine system (BJE). Results show that variations in the annual cycle of the net ecosystem metabolism from 1999-2003 depend on seasonal forces, such as basin-scale runoff and DIP loads. Results obtained through mass balance indicate large amounts of anthropogenic nutrient inputs to the system. These loads act as a source for dissolved inorganic phosphorus during the dry and rainy season. The loads of dissolved inorganic nitrogen act as sinks throughout the year. During the winter, the BJE basin exceeded the values reported for DIP and DIN (moles km^-2 yr^-1) in the literature for basins of up to 1000 km². These seasonal oscillations of heterotrophy and autotrophy show a moderate tendency to heterotrophy, indicating that the system passes to liquid production stages of organic matter when production surpasses mineralization (September-December) and liquid mineralization stages (March-August).

The linear regression between p-r and the residence time shows lower entropy in the dry season and autotrophy at lower rates than during the high residence times but, this needs to be confirmed in future studies.

Also evident is the importance of denitrification in the BJE, which establishes that the system is a net denitrificator at moderate rates, probably in association with the degradation of labile organic matter originated from sewage.

Consequently, both heterotrophy and denitrification are enhanced by the production of carbon and nitrogen during the rainy season, whereas heterotrophic systems mainly depend on the inputs or loads of organic carbon of the adjacent systems.

We considered that high-density human occupation in the basin contributes significantly to N and P emissions throughout the year. High per capita loads of N and P indicate a scenario of high population density and high runoff. However, it seems important to recognize that ignoring the uptake and release of nutrients (N and P) by abundant mineral particles in the estuary may cause errors in nutrient balances, although the evidence reported here helps us to understand the main processes driving the metabolism of poorly studied typical small low-latitude estuaries.

ACKNOWLEDGMENTS

We thank the Pernambuco State Water Resources Agency (SRH) (Secretaria de Recursos Hídricos) and Pernambuco State Environmental Agency (CPRH) (Agência Estadual de Meio Ambiente e Recursos Hídricos) for their cooperation regarding the field data used in this work. The authors would like to thank the Brazilian National Council of Scientific and Technological Development - CNPq (Conselho Nacional de Desenvolvimento Científico e Tecnológico) under the scope of the Project BIO-NE (grant 558143/2009-1). We are also grateful to Dr. Monica Costa for her manuscript review and insightful comments.

Manuscript received on November 3, 2009; accepted for publication on October 27, 2010

• ARAUJO M, COSTA MF, AURELIANO JT AND SILVA MA. 2008. Mathematical modelling of hydrodynamics and water quality in a tropical reservoir, Northeast Brazil. BJAST 12: 19-30.
• ARAUJO M, MEDEIROS C AND RIBEIRO C. 1999. Energy balance and time-scales of mixing and stratification in the Jaboatão estuary, NE-Brazil. Braz J Oceanogr 47: 145-154.
• BERGER WH (ED). 1989. Global maps of ocean productivity. In: SCHLESINGER WH (Ed), Biogeochemistry an analysis of global change, New York: Academic Press, New York, USA, p. 301-307.
• BIDONE ED. 2000. Análise econômica-ambiental aplicada à contaminação de águas fluviais de pequenas bacias costeiras do Estado do Rio de Janeiro. In: ESTEVES FA AND LACERDA LD (Eds), Ecologia de restingas e lagoas costeiras, Rio de Janeiro, Ed. UFRJ, Rio de Janeiro, Brasil, p. 371-394.
• BIDONE ED AND LACERDA LD. 2002. A preliminary approach of the link between socio-economic and natural indicators into a driver-pressure-impact-response framework case study: Guanabara Bay Basin, Rio de Janeiro, Brazil. In: LACERDA LD ET AL. (Eds), South American Basins: LOICZ global change assessment and synthesis of river catchment - coastal sea interaction and human dimensions, Texel, LOICZ Reports and Studies No. 21, Texel, The Netherlands, 212 p.
• BISWAS H, MUKHOPADHYAY SK, DE TK, SEN S AND JANA TK. 2004. Biogenic controls on the air-water carbon dioxide exchange in the Sundarban mangrove environment, northeast coast of Bay of Bengal, India. Limnol Oceanogr 49(1): 95-101.
• BNDO - BANCO NACIONAL DE DADOS OCEANOGRÁFICOS. 2004. Serviço de banco de dados oceanográficos. <http://www.mar.mil.br/dhn/chm/bndo/>
• BRANCO ES. 2002. Variação Sazonal e Espacial da Biomassa Fitoplanctônica Relacionada com Parâmetros Hidrológicos no Estuário de Barra das Jangadas (Jaboatão dos Guararapes - Pernambuco - Brasil). Trop Ocean 30: 79-96.
• BRANCO ES. 2006. Variação sazonal das algas planctônicas correlacionadas com parâmetros ambientais no estuário de Barra das Jangadas (Jaboatão dos Guararapes - PE - Brasil). Bol Téc Cient CEPENE 14: 17-23.
• CPRH - AGÊNCIA ESTADUAL DE MEIO AMBIENTE E RECURSOS HÍDRICOS. 2003. Relatório de monitoramento de bacias hidrográficas do estado de Pernambuco 1999- 2003, Recife. <http://www.cprh.pe.gov.br>
• DALE A AND PREGO R. 2005. Net autotrophy and heterotrophy in the Pontevedra Ria upwelling system (NW Iberian margin). Cienc Mar 31: 213-220.
• EKAU W AND KNOPPERS B. 1999. An introduction to the pelagic system of the north-east and east Brazilian shelf. Arch Fish Mar Res 47: 113-132.
• EYRE BD AND MCKEE LJ. 2002. Carbon, nitrogen and phosphorus budgets for a shallow subtropical coastal embayment (Moreton Bay, Australia). Limnol Oceanogr 47: 1043-1055.
• GATTUSO JP, FRANKIGNOULLE M AND WOLLAST R. 1998. Carbon and carbonate metabolism in coastal aquatic ecosystems. Annual Review of Ecol and Syst 29: 405-434.
• GORDON JR DC, BOUDREAU PR, MANN KH, ONG J-E, SILVERT WL, SMITH SV, WATTAYAKORN G, WULFF F AND YANAGI T. 1996. LOICZ Biogeochemical Modelling Guidelines. LOICZ Reports and Studies No. 5, 96 p.
• HESSEN DO, AGREN GI, ADERSON TR, ELSER JJ AND DE RUITER PC. 2004. Carbon sequestration in ecosystems: the role of stoichiometry. Ecology 85: 1179-1192.
• HOLLAND HD. 1978. The chemistry of the atmosphere and oceans. New York, Willey, 351 p.
• IBGE - INSTITUTO BRASILEIRO DE GEOGRAFIA E ESTATÍSTICA. 2000. Censo demográfico. <http://www.ibge.gov.br/cidadessat/default.
  » link
• INMET - INSTITUTO NACIONAL DE METEOROLOGIA. 2003. Relatório mensal de dados meteorológicos 1999- 2003. <http://www.inmet.gov.br>
• IRIGOIEN X AND CASTEL JC. 1997. Light limitation and distribution of Chlorophyll pigments in the highly turbid estuary: the Gironde (SW France). Estuar Coast Shelf Sci 44: 507-517.
• KNOPPERS BA, FIGUEIREDO AG AND EKAU W. 1999. The coast and shelf of east and northeast Brazil and material transport. Geo-Mar Lett 19: 171-178.
• KNOPPERS BA, SOUZA WFL, EKAU W, FIGUEIREDO AG AND SOARES-GOMES A. 2009. A Interface Terra-Mar do Brasil. In: PEREIRA RC AND SOARES-GOMES A (Eds), Biologia Marinha. 2Ş ed., Rio de Janeiro: Interciência, p. 529-553.
• LAL R, KIMBLE JM AND STEWART BA. 2000. Global climate change and tropical ecosystems, Boca Raton: CRC Press, 438 p.
• MARONE E, MACHADO E, LOPES R AND DA SILVA ET. 2005. Land-Ocean fluxes in the Paranaguá Bay estuarine system, southern Brazil. Braz J Oceanogr 53(3/4): 169- 181.
• MEYERS T, SICKLES J, DENNIS R, RUSSELL RK, GALLOWAY J AND CHURCH T. 2001. Atmospheric nitrogen deposition to coastal estuaries and their watersheds. In: VALIGURA RA, ALEXANDER RB, CASTRO MS, MEYERS TP, PAERL HW, STACEY PE AND TURNER RE (Eds), Nitrogen loading in coastalwater bodies: An atmospheric perspective, Washington: American Geophysical Union, p. 401-410.
• MUKHOPADHYAY SK, BISWAS H, DE TK AND JANA TK. 2006. Fluxes of nutrients from tropical River Hooghly at the land-ocean boundary of Sundarbans NE Coast of Bay of Bengal, India. J Marine Syst 62: 9-21.
• NORIEGA CED, MUNIZ K, ARAUJO MC, TRAVASSOS RK AND NEUMANN-LEITÃO S. 2005a. Fluxos de nutrientes inorgânicos dissolvidos em um estuário tropical - Barra das Jangadas - PE, Brasil. Trop Ocean 33: 129-139.
• NORIEGA CED, MUNIZ K, FEITOSA FA, FLORES-MONTES MJ, GREGO CK, SCHEIDT G AND SILVA HP. 2005b. Distribuição espacial da biomassa fitoplanctônica e sua relação com os sais nutrientes em um estuário tropical (Barra das Jangadas - PE - Brasil). Arq Cien Mar 38: 5-21.
• NORIEGA CED, MUNIZ K, FLORES-MONTES MJ, MACÊ- DO S, ARAUJO M, FEITOSA FA AND LACERDA S. 2009. Series temporales de variables hidrobiológicas en un estuario tropical (Brasil). Rev Biol Mar Oceanog 44(1): 93-108.
• PASSAVANTE JZ AND FEITOSA FA. 2004. Dinâmica da produtividade fitoplanctônica na zona costeira marinha. In: ESKINAZI-LEÇA E, NEUMANN-LEITAO S AND COSTA MF (Eds), Oceanografia: Um cenário tropical, Recife: Edições Bagaço, Recife, Brasil, p. 425-440.
• RAM PAS, NAIR S AND CHANDRAMOHAN D. 2003. Seasonal shift in net ecosystem production in a tropical estuary. Limnol Oceanogr 48: 1601-1607.
• SAN DIEGO-MCGLONE ML, SMITH SV AND NICOLAS VF. 2000. Stoichiometric interpretations of C:N:P ratios in organic waste materials. Mar Pollut Bull 40: 325-330.
• SCHINDLER DW. 1978. Factors regulating phytoplankton production and standing crop in worlds freshwaters. Limnol Oceanogr 23: 478-486.
• SCHLESINGER WH. 1997. Biogeochemistry: An analysis of global change. 2nd ed. San Diego: Academic Press. 430 p.
• SECTMA - SECRETARIA DE CIÊNCIA, TECNOLOGIA E MEIO AMBIENTE. 1999. Plano estadual de recursos hídricos, PERHPE. <http://www.inmet.gov.pe.br>
• SMITH SV AND CAMACHO V. 2000. Flujos de CNP en la zona costera. <http://nest.su.se/mnode/Methods/powerpoint/LOICZoverview_sp.ppt>
• SMITH SV AND HOLLIBAUGH JT. 1997. Annual cycle and interannual variability of net and gross ecosystem metabolism in a temperate climate embayment. Ecol Monogr 67: 509-533.
• SMITH SV, SWANEY DP, BUDDEMEIER RW, SCARSBROOK MR, EATHERHEAD MA, HUMBORG C, ERIKSSON H AND HANNERZ F. 2005. River nutrient loads and catchment size. Biogeochemistry 75: 83-107.
• SMITH SV ET AL. 2003. Humans, hydrology, and the distribution of inorganic nutrient loading to the ocean. Bio- Science 53: 235-245.
• SOUZA MFL. 2000. Rio Sergipe and PiauíRiver Estuaries. In: SMITH SV, DUPRA V, MARSHALL CROSSLAND JI AND CROSSLAND CJ (Eds), Estuarine systems of the South American region: carbon, nitrogen and phosphorus fluxes, Texel, LOICZ Reports and Studies 15, Texel, The Netherlands, p. 6-17.
• SWANEY DP AND SMITH SV. 2003. Guidelines for constructing nutrient budgets of coastal systems. In: CROSSLAND CJ, KREMER HH, LINDEBOOM HJ, MARSHALLCROSSLAND JI AND LE TISSIER MDA (Eds) Coastal fluxes in the anthropocene, Berlin, Springer, Berlin, Deutshland, p. 110-111.
• TAPPIN AD. 2002. An examination of the fluxes of nitrogen and phosphorus in temperate and tropical estuaries: Current estimates and uncertainties. Estuar Coast Shelf Sci 55: 885-901.
• WEBSTER IT, PARSLOW JS AND SMITH SV. 2000. Implications of spatial and temporal for biogeochemical budgets of estuaries. Estuaries 23: 341-350.
• WEPENER V. 2007. Carbon, nitrogen and phosphorus fluxes in four sub-tropical estuaries of northern KwaZulu-Natal: Case studies in the application of a mass balance approach. Water SA 33: 203-214.

Publication Dates

• Publication in this collection
  03 June 2011
• Date of issue
  June 2011