Are closure terms appropriate or necessary descriptors of zooplankton loss in nutrient–phytoplankton–zooplankton type models?

Our current knowledge of plankton ecology ascribes a large proportion of zooplankton losses to zooplankton cannibalism and carnivory, rather than via the activity of higher trophic levels beyond the plankton. However, planktonic ecosystem models, such as the widely used nutrient–phytoplankton–zooplankton (NPZ) type models, typically represent all zooplankton losses by mathematically (rather than biologically) justified closure functions. Even where it is assumed that these closure functions include zooplanktonic cannibalism and carnivory, these processes are not explicitly implemented within the grazing function of the zooplankton. Here it is argued that this representation of zooplankton losses through “closure” terms within planktonic food web models is neither appropriate nor necessary. The general consequences of implementing a simple function incorporating zooplankton cannibalism and carnivory (intra-guild predation) within a planktonic food web model are compared against models implementing different types of traditional closure functions. While the modelled biomass outputs may appear similar, the fate of annual primary production and f-ratios vary widely. There appears no justification for the continued use of traditional closure term to depict zooplankton loss processes on biological or modelling arguments. To do so can seriously misrepresent the fate of primary production and thence trophic dynamics.     

Food web structure and trophic interactions of the tropical highland lake Hayq, Ethiopia

We generated a mass-balance model to figure out the food web structure and trophic interactions of the major functional groups of the Ethiopian highland Lake Hayq. Moreover, the study lay down a baseline data for future ecosystem-based investigations and management activities. Extensive data collection has been taken place between October 2007 and May 2009. Ecotrophic efficiency (EE) of several functional groups including phytoplankton (0.8) and detritus (0.85) was high indicating the utilization of the groups within the system. However, the EE of Mesocyclops (0.03) and Thermocyclops (0.30) was very low implying these resources were rather a ‘sink’ in the trophic hierarchy. Flows based on aggregated trophic level sensu Lindeman revealed the importance of both phytoplankton and detritus to higher trophic levels. The computed average transfer efficiency of 11.5% for the first four trophic levels was within the range for highly efficient African lakes. The primary production to respiration (P/R) ratio (1.05) of Lake Hayq indicates the maturity of the ecosystem. We also modeled the food-web by excluding Tilapia and reduced phytoplankton biomass to get insight into the mass balance before Tilapia was introduced. The analysis resulted in a lower system omnivory index (SOI = 0.016) and a reduced P/R ratio (0.13) that described the lake as immature ecosystem, suggesting the introduction of Tilapia might have contributed to the maturity of the lake. Tilapia in Lake Hayq filled an ecological empty niche of pelagic planktivores, and contributed for the better transfer efficiency observed from primary production to fish yield.     

The effects of seston food quality on planktonic food web patterns

In planktonic food webs, the conversion rate of plant material to herbivore biomass is determined by a variety of factors such as seston biochemical/elemental composition, phytoplankton cell morphology, and colony architecture. Despite the overwhelming heterogeneity characterizing the plant–animal interface, plankton population models usually misrepresent the food quality constraints imposed on zooplankton growth. In this study, we reformulate the zooplankton grazing term to include seston food quality effects on zooplankton assimilation efficiency and examine its ramifications on system stability. Using different phytoplankton parameterizations with regards to growth strategies, light requirements, sinking rates, and food quality, we examined the dynamics induced in planktonic systems under varying zooplankton mortality/fish predation, light conditions, nutrient availability, and detritus food quality levels. In general, our analysis suggests that high food quality tends to stabilize the planktonic systems, whereas unforced oscillations (limit cycles) emerge with lower seston food quality. For a given phytoplankton specification and resource availability, the amplitude of the plankton oscillations is primarily modulated from zooplankton mortality and secondarily from the nutritional quality of the alternative food source (i.e., detritus). When the phytoplankton community is parameterized as a cyanobacterium-like species, conditions of high nutrient availability combined with high zooplankton mortality led to phytoplankton biomass accumulation, whereas a diatom-like parameterization resulted in relatively low phytoplankton to zooplankton biomass ratios highlighting the notion that high phytoplankton food quality allows the zooplankton community to sustain relatively high biomass and to suppress phytoplankton biomass to low levels. During nutrient and light enrichment conditions, both phytoplankton and detritus food quality determine the extent of the limit cycle region, whereas high algal food quality increases system resilience by shifting the oscillatory region towards lower light attenuation levels. Detritus food quality seems to regulate the amplitude of the dynamic oscillations following enrichment, when algal food quality is low. These results highlight the profitability of the alternative food sources for the grazer as an important predictor for the dynamic behavior of primary producer–grazer interactions in nature.     

Adaptations in a hierarchical food web of southeastern Lake Michigan

Two issues in ecological network theory are: (1) how to construct an ecological network model and (2) how do entire networks (as opposed to individual species) adapt to changing conditions? We present a novel method for constructing an ecological network model for the food web of southeastern Lake Michigan (USA) and we identify changes in key system properties that are large relative to their uncertainty as this ecological network adapts from one time point to a second time point in response to multiple perturbations. To construct our food web for southeastern Lake Michigan, we followed the list of seven recommendations outlined in Cohen et al. [Cohen, J.E., et al., 1993. Improving food webs. Ecology 74, 252–258] for improving food webs. We explored two inter-related extensions of hierarchical system theory with our food web; the first one was that subsystems react to perturbations independently in the short-term and the second one was that a system's properties change at a slower rate than its subsystems’ properties. We used Shannon's equations to provide quantitative versions of the basic food web properties: number of prey, number of predators, number of feeding links, and connectance (or density). We then compared these properties between the two time-periods by developing distributions of each property for each time period that took uncertainty about the property into account. We compared these distributions, and concluded that non-overlapping distributions indicated changes in these properties that were large relative to their uncertainty. Two subsystems were identified within our food web system structure (p < 0.001). One subsystem had more non-overlapping distributions in food web properties between Time 1 and Time 2 than the other subsystem. The overall system had all overlapping distributions in food web properties between Time 1 and Time 2. These results supported both extensions of hierarchical systems theory. Interestingly, the subsystem with more non-overlapping distributions in food web properties was the subsystem that contained primarily benthic taxa, contrary to expectations that the identified major perturbations (lower phosphorous inputs and invasive species) would more greatly affect the subsystem containing primarily pelagic taxa. Future food-web research should employ rigorous statistical analysis and incorporate uncertainty in food web properties for a better understanding of how ecological networks adapt.     

Evidence for the dominance of indirect effects in 50 trophic ecosystem networks

Indirect effects are powerful influences in ecosystems that may maintain species diversity and alter apparent relationships between species in surprising ways. Here, we applied network environ analysis to 50 empirically-based trophic ecosystem models to test the hypothesis that indirect flows dominate direct flows in ecosystem networks. Further, we used Monte Carlo based perturbations to investigate the robustness of these results to potential error in the underlying data. To explain our findings, we further investigated the importance of the microbial food web in recycling energy-matter using components of the Finn Cycling Index and analysis of environ centrality. We found that indirect flows dominate direct flows in 37/50 (74.0%) models. This increases to 31/35 (88.5%) models when we consider only models that have cycling structure and a representation of the microbial food web. The uncertainty analysis reveals that there is less error in the I/D values than the ±5% error introduced into the models, suggesting the results are robust to uncertainty. Our results show that the microbial food web mediates a substantial percentage of cycling in some systems (median = 30.2%), but its role is highly variable in these models, in agreement with the literature. Our results, combined with previous work, strongly suggest that indirect effects are dominant components of activity in ecosystems.     

A thermodynamic model of freshwater Antarctic lake ice

Antarctic lakes with simple plankton ecosystems are believed to be sensitive biological indicators of climate change. Models of the physical environment, in particular the ice layer, support understanding of how the ecosystems respond to meteorological variables. This paper describes how data from a previously reported automatic measuring probe and meteorological data from Davis station were used to develop a detailed thermodynamic model of the ice layer on Crooked Lake, one of the largest and deepest freshwater lakes in Antarctica. The general model structure is similar to a previously reported model of sea ice but with modifications specific to the Antarctic freshwater lake case informed by the data. The model inputs are atmospheric variables as well as water temperature, ice albedo and the radiation extinction coefficient for the ice. Heat and radiation fluxes at the ice–air and ice–water boundaries are calculated using equations chosen for their suitability for the Antarctic. In the case of shortwave radiation, equations were fitted to data from the automatic probe. Using the heat fluxes to establish boundary conditions, and incorporating the known thermodynamic properties of ice, the temperature profile within the ice and the resulting growth and melt of the ice can be calculated. The model uses a largely mechanistic approach, with most equations taken from established thermodynamic theories or empirical studies and only one adjustable parameter related to the sensible heat flux from the water, which is not easily calculated from the available data. It was found to accurately reproduce ice temperature and ice thickness data for the year 2003, with r² = 0.89, n = 2005. Finally, the model was simplified to run with air temperature as the only input variable and was shown to perform well—this suggests that freshwater lake ice is affected more by air temperature than any other variable, and is therefore a useful indicator of climate change in its own right.     

Calculating ecological carrying capacity of shellfish aquaculture using mass-balance modeling: Narragansett Bay, Rhode Island

Increasing growth in the aquaculture industry demands ecosystem-based techniques for management if that growth is to be ecologically sustainable and promote equity among users of the ecosystems in which it occurs. Models of carrying capacity can be used to responsibly limit the growth of aquaculture in increasingly crowded coastal areas. Narragansett Bay, Rhode Island, USA is one such crowded coastal region experiencing a rapid increase in bivalve aquaculture. An ecosystem mass-balance model was used to calculate the ecological carrying capacity of bivalve aquaculture. Cultured oyster biomass is currently at 0.47 t km⁻² and could be increased 625 times without exceeding the ecological carrying capacity of 297 t km⁻². This translates to approximately 38,950 t of harvested cultured oysters annually which is 4 times the total estimated annual harvest of finfish. This potential for growth is due to the high primary productivity and large energy throughput to detritus of this ecosystem. Shellfish aquaculture has potential for continued growth and is unlikely to become food limited due, in part, to the large detritus pool.     

A system-dynamic model on the competitive growth between Potamogeton malaianus Miq. and Spirogyra sp.

A system-dynamic model has been built to evaluate the competition between submerged macrophytes Potamogeton malaianus Miq. (PM) and filamentous green algae Spirogyra sp. (SP). The data background is based on a spring–summer and an autumn–winter experiment carried out in artificial field ponds. The experiments had the aim to acquire a knowledge base necessary to a successful restoration of submerged macrophyte vegetation in Lake Taihu, China by use of P. malaianus Miq. The model mainly focuses on variations in water volume; biomass dynamics of P. malaianus Miq., Spirogyra sp. and zoobenthos; nutrients cycling between water column, P. malaianus Miq., Spirogyra sp., zoobenthos, detritus and sediment. Sixteen state variables are included in the model: biomass of P. malaianus Miq., Spirogyra sp. and zoobenthos; nitrogen in sediments, detritus, in P. malaianus Miq., in Spirogyra sp. and zoobenthos; total dissolved nitrogen; phosphorus in sediments, detritus, in P. malaianus Miq., in Spirogyra sp. and in zoobenthos; total dissolved phosphorus, and water volume of the experiment pond. The calibration and validation of the model show a good accordance with the results of the spring–summer experiment and the autumn–winter experiment.     

Two-way coupling versus one-way forcing of plankton and fish models to predict ecosystem changes in the Benguela

‘End-to-end’ models have been adopted in an attempt to capture more of the processes that influence the ecology of marine ecosystems and to make system wide predictions of the effects of fishing and climate change. Here, we develop an end-to-end model by coupling existing models that describe the dynamics of low (ROMS–N₂P₂Z₂D₂) and high trophic levels (OSMOSE). ROMS–N₂P₂Z₂D₂ is a biogeochemical model representing phytoplankton and zooplankton seasonal dynamics forced by hydrodynamics in the Benguela upwelling ecosystem. OSMOSE is an individual-based model representing the dynamics of several species of fish, linked through opportunistic and size-based trophic interactions. The models are coupled through a two-way size-based predation process. Plankton provides prey for fish, and the effects of predation by fish on the plankton are described by a plankton mortality term that is variable in space and time. Using the end-to-end model, we compare the effects of two-way coupling versus one-way forcing of the fish model with the plankton biomass field. The fish-induced mortality on plankton is temporally variable, in part explained by seasonal changes in fish biomass. Inclusion of two-way feedback affects the seasonal dynamics of plankton groups and usually reduces the amplitude of variation in abundance (top-down effect). Forcing and coupling lead to different predicted food web structures owing to changes in the dominant food chain which is supported by plankton (bottom-up effect). Our comparisons of one-way forcing and two-way coupling show how feedbacks may affect abundance, food web structure and food web function and emphasise the need to critically examine the consequences of different model architectures when seeking to predict the effects of fishing and climate change.     

Integrative trophic network assessments of a lentic ecosystem by key ecological approaches of water chemistry, trophic guilds, and ecosystem health assessments along with an ECOPATH model

The objective of this study was to describe the trophic structure and energy flow in a lentic ecosystem in South Korea. Physicochemical water conditions were evaluated along with the reservoir ecosystem health using a multimetric IBI model. Nutrient analyses of the reservoir showed a nutrient rich and hypereutrophic system. Guild analysis revealed that tolerant and omnivorous species dominated the ecosystem. Tolerant fish, as a proportion of the number of individuals, were associated (R² > 0.90, p < 0.01) with TN and TP, the key indicators of trophic state in lentic ecosystems. The mean Reservoir Ecosystem Health Assessment (REHA) score was 19.3 during the study, which was judged as in ‘fair to poor’ condition. A trophic analysis of the reservoir estimated by the ECOPATH model shows that most activity in terms of energy flow occurred in the lower part of the trophic web, where there was intensive use of primary producers as a food source. Consequently, of the 10 consumer groups, nine fell within trophic levels <2.8. Trophic levels (TL) estimated from the weighted average of prey trophic levels varied from 1.0 for phytoplankton, macrophytes, and detritus to 3.25 for the top predator, Pseudobagrus fulvidraco. Our integrated approach to trophic network analysis may provide a key tool for determining the effects of nutrient influx on energy flow pathways in lentic ecosystems.     

A process-based model designed for filling of large data gaps in tower-based measurements of net ecosystem productivity

In this paper we present a simple hybrid gap-filling model (GFM) designed with a minimum number of parameters necessary to capture the ecological processes important for filling medium-to-large gaps in Flux data. As the model is process-based, the model has potential to be used in filling large gaps exhibiting a broad range of micro-meteorological and site conditions. The GFM performance was evaluated using “Punch hole” and extrapolation experiments based on data collected in west-central New Brunswick. These experiments indicated that the GFM is able to provide acceptable results (r² > 0.80) when >500 data points are used in model parameterization. The GFM was shown to address daytime evolution of NEP reasonably well for a wide range of weather and site conditions. An analysis of residuals indicated that for the most part no obvious trends were evident; although a slight bias was detected in NEP with soil temperature. To explore the portability of the GFM across ecosystem types, a transcontinental validation was conducted using NEP and ancillary data from seven ecosystems along a north-south transect (i.e., temperature–moisture gradient) from northern Europe (Finland) to the Middle East (Israel). The GFM was shown to explain over 75% of the variability in NEP measured at most ecosystems, which strongly suggests that the GFM maybe successfully applied to forest ecosystems outside Canada.     

Optimization of ecosystem model parameters using spatio-temporal soil moisture information

Parameters in process-based terrestrial ecosystem models are often nonlinearly related to the water flux to the atmosphere, and they also change temporally and spatially. Therefore, for estimating soil moisture, process-based terrestrial ecosystem models inevitably need to specify spatially and temporally variant model parameters. This study presents a two-stage data assimilation scheme (TSDA) to spatially and temporally optimize some key parameters of an ecosystem model which are closely related to soil moisture. At the first stage, a simplified ecosystem model, namely the Boreal Ecosystem Productivity Simulator (BEPS), is used to obtain the prior estimation of daily soil moisture. After the spatial distribution of 0–10 cm surface soil moisture is derived from remote sensing, an Ensemble Kalman Filter is used to minimize the difference between the remote sensing model results, through optimizing some model parameters spatially. At the second stage, BEPS is reinitialized using the optimized parameters to provide the updated model predictions of daily soil moisture. TSDA has been applied to an arid and semi-arid area of northwest China, and the performance of the model for estimating daily 0–10 cm soil moisture after parameter optimization was validated using field measurements. Results indicate that the TSDA developed in this study is robust and efficient in both temporal and spatial model parameter optimization. After performing the optimization, the correlation (r²) between model-predicted 0–10 cm soil moisture and field measurement increased from 0.66 to 0.75. It is demonstrated that spatial and temporal optimization of ecosystem model parameters can not only improve the model prediction of daily soil moisture but also help to understand the spatial and temporal variation of some key parameters in an ecosystem model and the corresponding ecological mechanisms controlling the variation.     

Resistance and re-organization of an ecosystem in response to biological invasion: Some hypotheses

Using a dynamic model of Lake Chozas developed by Marchi et al. (2011), we tested three hypotheses about recovery of the indigenous community and water quality after radical changes caused by introduction of an invasive allochthonous crayfish, Procambarus clarkii:

1.

Can the lake resist the pressure of an invasive species, like P. clarkii, by adaptation?

2.

Can the ecosystem recover when all the crayfish are removed and low phosphorus concentrations persist in inflow water?

3.

Does the simulated recovery of submerged vegetation occur at a total phosphorus concentration below 100 mg TP m⁻³, as estimated by Scheffer et al. (1993), Scheffer (1997), Jeppesen et al. (1998) and Zhang et al. (2003)?

We obtained the following answers:

1.

Lake Chozas can at least partly resist by adaptation. A combination of possible parameter changes could lead to a significant increase in eco-exergy.

2.

Removal of the phosphorus represented by crayfish (by harvesting) implies complete recovery of the lake and its eco-exergy, albeit not necessarily with the same organisms having the same properties.

3.

The expected hysteresis created by introduction and harvesting of crayfish is observed under the following conditions: phytoplankton dominance at total phosphorus ≥ about 200-250 mg TP m⁻³ and submerged vegetation returns at total phosphorus < 100 mg TP m⁻³.

     

Ecological role of stomatopods (mantis shrimps) and potential impacts of trawling in a marine ecosystem of the southeast coast of India

Ecosystems are balanced by nature and each component in the system has a role in the sustenance of other components. A change in one component would invariably have an effect on others. Stomatopods (mantis shrimps) are common and ecologically important predatory crustaceans in tropical marine waters. The ecological role of mantis shrimps and potential impacts of trawling in a marine ecosystem were estimated using Ecopath with Ecosim (EwE) Version 5.0 software, by constructing a mass balanced Ecopath model of Parangipettai (Porto Novo) ecosystem. Based on fisheries information from the region, 17 ecological groups were defined including stomatopods. Both primary and secondary data on biomass, P/B, Q/B and diet composition were used as basic inputs. The mass balanced model gave a total system throughput of 14,756 t km⁻² year⁻¹. The gross efficiency of 0.000942 indicated higher contribution of lower food chain groups in the fishery though the mean trophic level was 3.08. The immature and developing stage of the ecosystem was indicated by the ratio of total primary production and total respiration (1.832) and the net system production (2643.30 t km⁻² year⁻¹). Key indices (flow to detritus, net efficiency and omnivory index), split mortality rates and mixed trophic impact of different ecological groups were obtained from the model. A flow diagram was constructed to illustrate the trophic interactions, which explained the biomass flows in the ecosystem with reference to stomatopods. Two temporal simulations were made, with 10 year durations in the mass balanced Ecopath model by using ecosim routine incorporated in EwE software. The effect of decrease in biomass of stomatopods in the ecosystem was well defined, in the first run with increase in stomatopod fishing mortality, and the group showed a high positive impact on benthopelagic fish biomass increase (129%). The simulation with increase in trawling efforts resulted in the biomass decline of different ecological groups as elasmobranchs to 1%, stomatopods to 2%, crabs and lobsters to 36%, cephalopods to 63%, mackerel to 78%, and shrimps to 89%. Present study warns stomatopod discards and further increase in trawling efforts in the region and it explained the need for ecosystem based fisheries management practices for the sustainability of marine fisheries.     

Switching between herbivory and carnivory by the planktonic marine copepod Calanus pacificus

Adult females of the omnivorous copepod Calanus pacificus, collected from the plankton off La Jolla, California, USA (June, 1978), fed disproportionately on the prey in greatest relative abundance when given mixtures of diatoms (Thalassiosira fluviatilis) and copepod (C. pacificus) nauplii as food. This switch from herbivorous to carnivorous behavior may be significant in nature during the decline of phytoplankton blooms. More generally, the widespread omnivorous habit among pelagic animals suggests a responsive and flexible trophic organization which contributes to the resiliency of planktonic communities in a dynamic physical environment.Contribution No. 1223 from the Department of Oceanography, University of Washington, Seattle, Washington 98195, USA     

Trophic relationships in the recently impounded Bagré reservoir in Burkina Faso

The trophic dynamics of Bagré reservoir which has been recently impounded in Burkina Faso was based on the data collected during 1997–1998 period using the Ecopath model and software. Total fish biomass is 22.63 t km⁻² and mainly represents trophic levels (TLs) 2 and 3. The trophic food chain is relatively long and the overall transfer efficiency is quite low. Grazing foodweb based on primary producers is prominent in the reservoir ecosystem and detritus plays a less significant role. Seasonal and long-term variations in water quality have significant influences on the lower TLs clearly showing a bottom-up functioning of the ecosystem. Environmental degradations, such as siltation occurring in the lake, suggest possible risks in limiting ecosystem productivity.     

Spatial and temporal hydrodynamic and water quality modeling analysis of a large reservoir on the South Carolina (USA) coastal plain

Two-dimensional, 31-segment, 61-channel hydrodynamic and water quality models of Lake Marion (surface area 330.7 km²; volume 1548.3×10⁶ m³) were developed using the WASP5 modeling system. Field data from 1985 to 1990 were used to parameterize the models. Phytoplankton kinetic rates and constants were obtained from a related in situ study; others from modeling literature. The hydrodynamic model was calibrated to estimates of daily lake volume; the water quality model was calibrated for ammonia, nitrate, ortho-phosphate, dissolved oxygen, chlorophyll-a, biochemical oxygen demand, organic nitrogen, and organic phosphorus. Water quality calibration suggested the model characterized phytoplankton and nutrient dynamics quite well. The model was validated (Kolmogorov–Smirnov two-sample goodness-of-fit test at P<0.05) by reparameterizing the nutrient loading functions using an independent set of field data. The models identified several factors that may contribute to the spatial variability previously reported from other research in the reservoir, despite the superficial absence of complex structure. Sensitivity analysis of the phytoplankton kinetic rates suggest that study site-specific estimates were important for obtaining model fit to field data. Sediment sources of ammonia (10–60 mg m⁻² day⁻¹) and phosphate (1–6 mg m⁻² day⁻¹) were important to achieve model calibration, especially during periods of high temperatures and low dissolved oxygen. This sediment flux accounted for 78% (nitrogen) and 50% (phosphorus) of the annual load. Spatial and temporal variability in the lake, reflected in the calibrated and validated models, suggest that ecological factors that influence phytoplankton productivity and nutrient dynamics are different in various parts of the lake. The WASP5 model as implemented here does not fully accommodate the ecological variability in Lake Marion due to model constraints on the specification of rate constants. This level of spatial detail may not be appropriate for an operational reservoir model, but as a research tool the models are both versatile and useful.     

Food resources of postlarval brown shrimp (Penaeus aztecus) in a Texas salt marsh

Field and laboratory experiments were conducted to identify the sources of food in the natural diet of postlarval brown shrimp (Penaeus aztecus Ives). A series of enclosures placed in East Lagoon (29°20N; 94°45W) on Galveston Island, Texas, USA, in May 1985, were used to evaluate the individual and combined contribution of Spartina alterniflora detritus, epiphytes of S. alterniflora, plankton, and demersal fauna in terms of differences in shrimp growth and carbon assimilation (stable carbon-isotope analysis). Demersal fauna (harpacticoid copepods, amphipods, tanaids and polychaete annelids), and plankton (>0.095 mm) accounted for approximately 53 and 47% of the growth of the postlarvae (11 to 22 mm rostrum-telson length), respectively, while the autochthonous plant substrates, S. alterniflora detritus and epiphytes, contributed little. Laboratory experiments confirm that a mixed diet consisting of both animal protein and phytoplankton promotes maximum growth. Our results indicate that plankton may be an important allochthonous source of carbon contributing to the growth and development of shrimp in the salt marsh.     

Local calibration of coliforms parameters of water quality problem at Igapó I Lake, Londrina, Paraná, Brazil

This article presents results concerning the local calibration of the transport parameters (longitudinal and transversal diffusions and decay coefficient) for a two-dimensional problem of water quality at Igapó I Lake, located in Londrina, Paraná, Brazil, using fecal coliforms as an indicator of water quality. The simulation of fecal coliforms concentrations all over the water body is conducted by means of a structured discretization of the geometry of Igapó I Lake, together with the finite difference and finite element methods. By using the velocity field, modeled by the Navier-Stokes and Poisson equations, the flow of fecal coliforms is described by means of a transport model, which considers advective and diffusive processes, as well as a process of fecal coliforms decay. In the checkpoint, the longitudinal and transversal diffusion coefficients and the coliforms decay coefficient that best fitted the value of the fecal coliforms concentration were D_x = D_y = 0.001 m²/h and k = 0.5 d⁻¹ = 0.02083 h⁻¹. A qualitative and quantitative analysis of the numerical simulations conducted in function of the diffusion coefficients and of the coliforms decay parameter provided a better understanding of the local water quality at Igapó I Lake.