Trophic dynamics in the northern Humboldt Current system: insights from stable isotopes and stomach content analyses

The northern Humboldt Current system (NHCS) off Peru is one of the most productive world marine regions. It represents less than 0.1% of the world ocean surface but presently sustains about 10% of the world fish catch, with the Peruvian anchovy or anchoveta Engraulis ringens as emblematic fish resource. Compared with other eastern boundary upwelling systems, the higher fish productivity of the NHCS cannot be explained by a corresponding higher primary productivity. On another hand, the NHCS is the region where El Niño, and climate variability in general, is most notable. Also, surface oxygenated waters overlie an intense and extremely shallow Oxygen Minimum Zone (OMZ). In this context, the main objective of this study is to better understand the trophic flows in the NHCS using both stomach content and stable isotope analyses. The study focuses on a variety of organisms from low trophic levels such as zooplankton to top predators (seabirds and fur seals). The approach combines both long-term and specific studies on emblematic species such as anchoveta, and sardine Sardinops sagax and a more inclusive analysis considering the 'global' food web in the recent years (2008 – 2012) using stable isotope analysis. Revisiting anchovy and sardine we show that whereas phytoplankton largely dominated anchoveta and sardine diets in terms of numerical abundance, the carbon content of prey items indicated that zooplankton was by far the most important dietary component. Indeed for anchovy euphausiids contributed 67.5% of dietary carbon, followed by copepods (26.3%). Selecting the largest prey, the euphausiids, provide an energetic advantage for anchoveta in its ecosystem where oxygen depletion imposes strong metabolic constrain to pelagic fish. Sardine feed on smaller zooplankton than do anchoveta, with sardine diet consisting of smaller copepods and fewer euphausiids than anchoveta diet. Hence, trophic competition between sardine and anchovy in the northern Humboldt Current system is minimized by their partitioning of the zooplankton food resource based on prey size, as has been reported in other systems. These results suggest an ecological role for pelagic fish that challenges previous understanding of their position in the foodweb (zooplanktophagous instead of phytophagous), the functioning and the trophic models of the NHCS. Finally to obtain a more comprehensive vision of the relative trophic position of NHCS main components we used stable isotope analyses. For that purpose we analyzed the δ13C and δ15N stable isotope values of thirteen taxonomic categories collected off Peru from 2008 - 2011, i.e., zooplankton, fish, squids and air-breathing top predators. The δ15N isotope signature was strongly impacted by the species, the body length and the latitude. Along the Peruvian coast, the OMZ get more intense and shallow south of ~7.5ºS impacting the baseline nitrogen stable isotopes. Employing a linear mixed-effects modelling approach taking into account the latitudinal and body length effects, we provide a new vision of the relative trophic position of key ecosystem components. Also we confirm stomach content-based results on anchoveta Engraulis ringens and highlight the potential remarkable importance of an often neglected ecosystem component, the squat lobster Pleuroncodes monodon. Indeed, our results support the hypothesis according to which this species forage to some extent on fish eggs and larvae and can thus predate on the first life stages of exploited species. However, the δ13C values of these two species suggest that anchoveta and squat lobster do not exactly share the same habitat. This would potentially reduce some direct competition and/or predation.

Pescadores y aves marinas compitiendo por el mismo recurso: estrategias de forrajeo, interacciones y consecuencias

Con el examen conjunto de datos de seguimiento satelital de la flota pesquera y de marcaje electrónico de aves marinas, se analizan las interacciones entre la actividad pesquera y el forrajeo de las aves en periodo de reproducción. Se evidencia que las aves pueden mitigar la competencia con la pesca hasta cierto punto, forrajeando más lejos o quedándose más tiempo en el mar. Sin embargo, las aves en reproducción enfrentan a la vez un alto requerimiento energético para alimentar los pichones, y viajes de forrajeo limitados en distancia y tiempo por la necesidad de atender el nido. Para optimizar el éxito reproductivo de las aves marinas se recomienda: (1) Estimar, con la ayuda de modelos ecotróficos, qué cantidad de anchoveta sería conveniente ‘reservar’ para la alimentación de las aves; (2) Establecer zonas temporalmente cerradas a la pesca, alrededor de las principales colonias en los meses de reproducción de las aves. Los rangos máximos de forrajeo observados en guanayes y piqueros sugieren un radio de 50 a 100 km alrededor de las colonias, lo cual permitiría asegurar el forrajeo de estas especies en periodo de reproducción y así favorecer la sostenibilidad de sus poblaciones.

Modeling temporal variations of Gracilaria Greville and Hypnea J.V. Lamouroux (Rhodophyta) assemblages on a midlittoral reef platform at Piedade Beach, Pernambuco State, Brazil

The diversity of algal banks composed of species out the genera Gracilaria Greville and Hypnea J.V. Lamouroux have been impacted by commercial exploitation and coastal eutrophication. The present study sought to construct dynamic models based on algal physiology to simulate seasonal variations in the biomasses of Gracilaria and Hypnea an intertidal reef at Piedade Beach in Jaboatão dos Guararapes, Pernambuco State, Brazil. Five 20 × 20 cm plots in a reef pool on a midlittoral reef platform were randomly sampled during April, June, August, October, and December/2009 and in January and March/2010. Water temperature, pH, irradiance, oxygen and salinity levels as well as the concentrations of ammonia, nitrate and phosphate were measured at the sampling site. Forcing functions were employed in the model to represent abiotic factors, and algal decay was simulated with a dispersal function. Algal growth was modeled using a logistic function and was found to be sensitive to temperature and salinity. Maximum absorption rates of ammonia and phosphate were higher in Hypnea than in Gracilaria, indicating that the former takes up nutrients more efficiently at higher concentrations. Gracilaria biomass peaked at approximately 120 g (dry weight m-2) in March/2010 and was significantly lower in August/2009; Hypnea biomasses, on the other hand, did not show any significant variations among the different months, indicating that resource competition may influence the productivity of these algae.

(Table A1) Annual biomass of plants and animals recorded on Bylot Island between 1993-2009

Determining the manner in which food webs will respond to environmental changes is difficult because the relative importance of top-down vs. bottom-up forces in controlling ecosystems is still debated. This is especially true in the Arctic tundra where, despite relatively simple food webs, it is still unclear which forces dominate in this ecosystem. Our primary goal was to assess the extent to which a tundra food web was dominated by plant-herbivore or predator--rey interactions. Based on a 17-year (1993-2009) study of terrestrial wildlife on Bylot Island, Nunavut, Canada, we developed trophic mass balance models to address this question. Snow Geese were the dominant herbivores in this ecosystem, followed by two sympatric lemming species (brown and collared lemmings). Arctic foxes, weasels, and several species of birds of prey were the dominant predators. Results of our trophic models encompassing 19 functional groups showed that <10% of the annual primary production was consumed by herbivores in most years despite the presence of a large Snow Goose colony, but that 20-100% of the annual herbivore production was consumed by predators. The impact of herbivores on vegetation has also weakened over time, probably due to an increase in primary production. The impact of predators was highest on lemmings, intermediate on passerines, and lowest on geese and shorebirds, but it varied with lemming abundance. Predation of collared lemmings exceeded production in most years and may explain why this species remained at low density. In contrast, the predation rate on brown lemmings varied with prey density and may have contributed to the high-amplitude, periodic fluctuations in the abundance of this species. Our analysis provided little evidence that herbivores are limited by primary production on Bylot Island. In contrast, we measured strong predator-prey interactions, which supports the hypothesis that this food web is primarily controlled by top-down forces. The presence of allochthonous resources subsidizing top predators and the absence of large herbivores may partly explain the predominant role of predation in this low-productivity ecosystem.

Merging Resource Availability with Isotope Mixing Models: The Role of Neutral Interaction Assumptions

Background: Bayesian mixing models have allowed for the inclusion of uncertainty and prior information in the analysis of trophic interactions using stable isotopes. Formulating prior distributions is relatively straightforward when incorporating dietary data. However, the use of data that are related, but not directly proportional, to diet (such as prey availability data) is often problematic because such information is not necessarily predictive of diet, and the information required to build a reliable prior distribution for all prey species is often unavailable. Omitting prey availability data impacts the estimation of a predator's diet and introduces the strong assumption of consumer ultrageneralism (where all prey are consumed in equal proportions), particularly when multiple prey have similar isotope values. Methodology: We develop a procedure to incorporate prey availability data into Bayesian mixing models conditional on the similarity of isotope values between two prey. If a pair of prey have similar isotope values (resulting in highly uncertain mixing model results), our model increases the weight of availability data in estimating the contribution of prey to a predator's diet. We test the utility of this method in an intertidal community against independently measured feeding rates. Conclusions: Our results indicate that our weighting procedure increases the accuracy by which consumer diets can be inferred in situations where multiple prey have similar isotope values. This suggests that the exchange of formalism for predictive power is merited, particularly when the relationship between prey availability and a predator's diet cannot be assumed for all species in a system.

Combining food web and species distribution models for improved community projections

The ability to model biodiversity patterns is of prime importance in this era of severe environmental crisis. Species assemblage along environmental gradient is subject to the interplay of biotic interactions in complement to abiotic environmental filtering. Accounting for complex biotic interactions for a wide array of species remains so far challenging. Here, we propose to use food web models that can infer the potential interaction links between species as a constraint in species distribution models. Using a plant-herbivore (butterfly) interaction dataset, we demonstrate that this combined approach is able to improve both species distribution and community forecasts. Most importantly, this combined approach is very useful in rendering models of more generalist species that have multiple potential interaction links, where gap in the literature may be recurrent. Our combined approach points a promising direction forward to model the spatial variation of entire species interaction networks. Our work has implications for studies of range shifting species and invasive species biology where it may be unknown how a given biota might interact with a potential invader or in future climate.

Climate-based empirical models show biased predictions of butterfly communities along environmental gradients

A better understanding of the factors that mould ecological community structure is required to accurately predict community composition and to anticipate threats to ecosystems due to global changes. We tested how well stacked climate-based species distribution models (S-SDMs) could predict butterfly communities in a mountain region. It has been suggested that climate is the main force driving butterfly distribution and community structure in mountain environments, and that, as a consequence, climate-based S-SDMs should yield unbiased predictions. In contrast to this expectation, at lower altitudes, climate-based S-SDMs overpredicted butterfly species richness at sites with low plant species richness and underpredicted species richness at sites with high plant species richness. According to two indices of composition accuracy, the Sorensen index and a matching coefficient considering both absences and presences, S-SDMs were more accurate in plant-rich grasslands. Butterflies display strong and often specialised trophic interactions with plants. At lower altitudes, where land use is more intense, considering climate alone without accounting for land use influences on grassland plant richness leads to erroneous predictions of butterfly presences and absences. In contrast, at higher altitudes, where climate is the main force filtering communities, there were fewer differences between observed and predicted butterfly richness. At high altitudes, even if stochastic processes decrease the accuracy of predictions of presence, climate-based S-SDMs are able to better filter out butterfly species that are unable to cope with severe climatic conditions, providing more accurate predictions of absences. Our results suggest that predictions should account for plants in disturbed habitats at lower altitudes but that stochastic processes and heterogeneity at high altitudes may limit prediction success of climate-based S-SDMs.

Diatoms as indicators of lake trophic status

Relationships between surface sediment diatom assemblages and lake trophic status were studied in 50 Canadian Precambrian Shield lakes in the Muskoka-Haliburton and southern Ontario regions. The purpose of this study was to develop mathematical regression models to infer lake trophic status from diatom assemblage data. To achieve this goal, however, additional investigations dealing with the evaluation of lake trophic status and the autecological features of key diatom species were carried out. Because a unifying index and classification for lake trophic status was not available, a new multiple index was developed in this study, by the computation of the physical, chemical and biological data from 85 south Ontario lakes. By using the new trophic parameter, the lake trophic level (TL) was determined: TL = 1.37 In[1 +(TP x Chl-a / SD)], where, TP=total phosphorus, Chl-a=chlorophyll-a and SD=Secchi depth. The boundaries between 7 lake trophic categories (Ultra-oligotrophic lakes: 0-0.24; Oligotrophic lakes: 0.241-1.8; Oligomesotrophic lakes: 1.813.0; Mesotrophic lakes: 3.01-4.20; Mesoeutrophic lakes: 4.21-5.4; Eutrophic lakes: 5.41-10 and Hyper-eutrophic lakes: above 10) were established. The new trophic parameter was more convenient for management of water quality, communication to the public and comparison with other lake trophic status indices than many of the previously published indices because the TL index attempts to Increase understanding of the characteristics of lakes and their comprehensive trophic states. It is more reasonable and clear for a unifying determination of true trophic states of lakes. Diatom specIes autecology analysis was central to this thesis. However, the autecological relationship of diatom species and lake trophic status had not previously been well documented. Based on the investigation of the diatom composition and variety of species abundance in 30 study lakes, the distribution optima of diatom species were determined. These determinations were based on a quantitative method called "weighted average" (Charles 1985). On this basis, the diatom species were classified into five trophic categories (oligotrophic, oligomesotrophic, mesotrophic, mesoeutrophic and eutrophic species groups). The resulting diatom trophic status autecological features were used in the regressIon analysis between diatom assemblages and lake trophic status. When the TL trophic level values of the 30 lakes were regressed against their fi ve corresponding diatom trophic groups, the two mathematical equations for expressing the assumed linear relationship between the diatom assemblages composition were determined by (1) uSIng a single regression technique: Trophic level of lake (TL) = 2.643 - 7.575 log (Index D) (r = 0.88 r2 = 0.77 P = 0.0001; n = 30) Where, Index D = (0% + OM% + M%)/(E% + ME% + M%); 4 (2) uSIng a' multiple regressIon technique: TL=4.285-0.076 0%- 0.055 OM% - 0.026 M% + 0.033 ME% + 0.065 E% (r=0.89, r2=0.792, P=O.OOOl, n=30) There was a significant correlation between measured and diatom inferred trophic levels both by single and multiple regressIon methods (P < 0.0001, n=20), when both models were applied to another 20 test lakes. Their correlation coefficients (r2 ) were also statistically significant (r2 >0.68, n=20). As such, the two transfer function models between diatoms and lake trophic status were validated. The two models obtained as noted above were developed using one group of lakes and then tested using an entirely different group of lakes. This study indicated that diatom assemblages are sensitive to lake trophic status. As indicators of lake trophic status, diatoms are especially useful in situations where no local trophic information is available and in studies of the paleotrophic history of lakes. Diatom autecological information was used to develop a theory assessing water quality and lake trophic status.

A comparison of diatom phosphorus transfer functions and export coefficient models as tools for reconstructing lake nutrient histories

1. We compared the baseline phosphorus (P) concentrations inferred by diatom-P transfer functions and export coefficient models at 62 lakes in Great Britain to assess whether the techniques produce similar estimates of historical nutrient status. 2. There was a strong linear relationship between the two sets of values over the whole total P (TP) gradient (2-200 mu g TP L-1). However, a systematic bias was observed with the diatom model producing the higher values in 46 lakes (of which values differed by more than 10 mu g TP L-1 in 21). The export coefficient model gave the higher values in 10 lakes (of which the values differed by more than 10 mu g TP L-1 in only 4). 3. The difference between baseline and present-day TP concentrations was calculated to compare the extent of eutrophication inferred by the two sets of model output. There was generally poor agreement between the amounts of change estimated by the two approaches. The discrepancy in both the baseline values and the degree of change inferred by the models was greatest in the shallow and more productive sites. 4. Both approaches were applied to two lakes in the English Lake District where long-term P data exist, to assess how well the models track measured P concentrations since approximately 1850. There was good agreement between the pre-enrichment TP concentrations generated by the models. The diatom model paralleled the steeper rise in maximum soluble reactive P (SRP) more closely than the gradual increase in annual mean TP in both lakes. The export coefficient model produced a closer fit to observed annual mean TP concentrations for both sites, tracking the changes in total external nutrient loading. 5. A combined approach is recommended, with the diatom model employed to reflect the nature and timing of the in-lake response to changes in nutrient loading, and the export coefficient model used to establish the origins and extent of changes in the external load and to assess potential reduction in loading under different management scenarios. 6. However, caution must be exercised when applying these models to shallow lakes where the export coefficient model TP estimate will not include internal P loading from lake sediments and where the diatom TP inferences may over-estimate TP concentrations because of the high abundance of benthic taxa, many of which are poor indicators of trophic state.

Trophic interactions in a changing world: modelling aboveground-belowground interactions

The rate and scale of human-driven changes can exert profound impacts on ecosystems, the species that make them up and the services they provide that sustain humanity. Given the speed at which these changes are occurring, one of society's major challenges is to coexist within ecosystems and to manage ecosystem services in a sustainable way. The effect of possible scenarios of global change on ecosystem services can be explored using ecosystem models. Such models should adequately represent ecosystem processes above and below the soil surface (aboveground and belowground) and the interactions between them. We explore possibilities to include such interactions into ecosystem models at scales that range from global to local. At the regional to global scale we suggest to expand the plant functional type concept (aggregating plants into groups according to their physiological attributes) to include functional types of aboveground-belowground interactions. At the scale of discrete plant communities, process-based and organism-oriented models could be combined into "hybrid approaches" that include organism-oriented mechanistic representation of a limited number of trophic interactions in an otherwise process - oriented approach. Under global change the density and activity of organisms determining the processes may change non-linearly and therefore explicit knowledge of the organisms and their responses should ideally be included. At the individual plant scale a common organism-based conceptual model of aboveground-belowground interactions has emerged. This conceptual model facilitates the formulation of research questions to guide experiments aiming to identify patterns that are common within, but differ between, ecosystem types and biomes. Such experiments inform modelling approaches at larger scales. Future ecosystem models should better include this evolving knowledge of common patterns of aboveground-belowground interactions. Improved ecosystem models are necessary toots to reduce the uncertainty in the information that assists us in the sustainable management of our environment in a changing world. (C) 2004 Elsevier GmbH. All rights reserved.

Partitioning the relative fitness effects of diet and trophic morphology in the threespine stickleback

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Exploring community assembly through an individual-based model for trophic interactions

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Trait contributions to fish community assembly emerge from trophic interactions in an individual-based model

Community ecology seeks to understand and predict the characteristics of communities that can develop under different environmental conditions, but most theory has been built on analytical models that are limited in the diversity of species traits that can be considered simultaneously. We address that limitation with an individual-based model to simulate assembly of fish communities characterized by life history and trophic interactions with multiple physiological tradeoffs as constraints on species performance. Simulation experiments were carried out to evaluate the distribution of 6 life history and 4 feeding traits along gradients of resource productivity and prey accessibility. These experiments revealed that traits differ greatly in importance for species sorting along the gradients. Body growth rate emerged as a key factor distinguishing community types and defining patterns of community stability and coexistence, followed by egg size and maximum body size. Dominance by fast-growing, relatively large, and fecund species occurred more frequently in cases where functional responses were saturated (i.e. high productivity and/or prey accessibility). Such dominance was associated with large biomass fluctuations and priority effects, which prevented richness from increasing with productivity and may have limited selection on secondary traits, such as spawning strategies and relative size at maturation. Our results illustrate that the distribution of species traits and the consequences for community dynamics are intimately linked and strictly dependent on how the benefits and costs of these traits are balanced across different conditions. © 2012 Elsevier B.V.

Estimation of chlorophyll-a concentration and the trophic state of the Barra Bonita hydroelectric reservoir using oli/landsat-8 images

Reservoirs are artificial environments built by humans, and the impacts of these environments are not completely known. Retention time and high nutrient availability in the water increases the eutrophic level. Eutrophication is directly correlated to primary productivity by phytoplankton. These organisms have an important role in the environment. However, high concentrations of determined species can lead to public health problems. Species of cyanobacteria produce toxins that in determined concentrations can cause serious diseases in the liver and nervous system, which could lead to death. Phytoplankton has photoactive pigments that can be used to identify these toxins. Thus, remote sensing data is a viable alternative for mapping these pigments, and consequently, the trophic. Chlorophyll-a (Chl-a) is present in all phytoplankton species. Therefore, the aim of this work was to evaluate the performance of images of the sensor Operational Land Imager (OLI) onboard the Landsat-8 satellite in determining Chl-a concentrations and estimating the trophic level in a tropical reservoir. Empirical models were fitted using data from two field surveys conducted in May and October 2014 (Austral Autumn and Austral Spring, respectively). Models were applied in a temporal series of OLI images from May 2013 to October 2014. The estimated Chl-a concentration was used to classify the trophic level from a trophic state index that adopted the concentration of this pigment-like parameter. The models of Chl-a concentration showed reasonable results, but their performance was likely impaired by the atmospheric correction. Consequently, the trophic level classification also did not obtain better results.

Parameterization of bio-optical models for estimating chlorophyll-a concentration in a tropical eutrophic reservoir

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)