Trophic amplification of climate warming

Ecosystems can alternate suddenly between contrasting persistent states due to internal processes or external drivers. It is important to understand the mechanisms by which these shifts occur, especially in exploited ecosystems. There have been several abrupt marine ecosystem shifts attributed either to fishing, recent climate change or a combination of these two drivers. We show that temperature has been an important driver of the trophodynamics of the North Sea, a heavily fished marine ecosystem, for nearly 50 years and that a recent pronounced change in temperature established a new ecosystem dynamic regime through a series of internal mechanisms. Using an end-to-end ecosystem approach that included primary producers, primary, secondary and tertiary consumers, and detritivores, we found that temperature modified the relationships among species through nonlinearities in the ecosystem involving ecological thresholds and trophic amplifications. Trophic amplification provides an alternative mechanism to positive feedback to drive an ecosystem towards a new dynamic regime, which in this case favours jellyfish in the plankton and decapods and detritivores in the benthos. Although overfishing is often held responsible for marine ecosystem degeneration, temperature can clearly bring about similar effects. Our results are relevant to ecosystem-based fisheries management (EBFM), seen as the way forward to manage exploited marine ecosystems.

Influence of climate change and trophic coupling across four trophic levels in the Celtic Sea

Climate change has had profound effects upon marine ecosystems, impacting across all trophic levels from plankton to apex predators. Determining the impacts of climate change on marine ecosystems requires understanding the direct effects on all trophic levels as well as indirect effects mediated by trophic coupling. The aim of this study was to investigate the effects of climate change on the pelagic food web in the Celtic Sea, a productive shelf region in the Northeast Atlantic. Using long-term data, we examined possible direct and indirect ‘bottom-up’ climate effects across four trophic levels: phytoplankton, zooplankton, mid-trophic level fish and seabirds. During the period 1986–2007, although there was no temporal trend in the North Atlantic Oscillation index (NAO), the decadal mean Sea Surface Temperature (SST) in the Celtic Sea increased by 0.66±0.02°C. Despite this, there was only a weak signal of climate change in the Celtic Sea food web. Changes in plankton community structure were found, however this was not related to SST or NAO. A negative relationship occurred between herring abundance (0- and 1-group) and spring SST (0-group: p = 0.02, slope = −0.305±0.125; 1-group: p = 0.04, slope = −0.410±0.193). Seabird demographics showed complex species–specific responses. There was evidence of direct effects of spring NAO (on black-legged kittiwake population growth rate: p = 0.03, slope = 0.0314±0.014) as well as indirect bottom-up effects of lagged spring SST (on razorbill breeding success: p = 0.01, slope = −0.144±0.05). Negative relationships between breeding success and population growth rate of razorbills and common guillemots may be explained by interactions between mid-trophic level fish. Our findings show that the impacts of climate change on the Celtic Sea ecosystem is not as marked as in nearby regions (e.g. the North Sea), emphasizing the need for more research at regional scales.

Biomass changes and trophic amplification of plankton in a warmer ocean

Ocean warming can modify the ecophysiology and distribution of marine organisms, and relationships between species, with nonlinear interactions between ecosystem components potentially resulting in trophic amplification. Trophic amplification (or attenuation) describe the propagation of a hydroclimatic signal up the food web, causing magnification (or depression) of biomass values along one or more trophic pathways. We have employed 3-D coupled physical-biogeochemical models to explore ecosystem responses to climate change with a focus on trophic amplification. The response of phytoplankton and zooplankton to global climate-change projections, carried out with the IPSL Earth System Model by the end of the century, is analysed at global and regional basis, including European seas (NE Atlantic, Barents Sea, Baltic Sea, Black Sea, Bay of Biscay, Adriatic Sea, Aegean Sea) and the Eastern Boundary Upwelling System (Benguela). Results indicate that globally and in Atlantic Margin and North Sea, increased ocean stratification causes primary production and zooplankton biomass to decrease in response to a warming climate, whilst in the Barents, Baltic and Black Seas, primary production and zooplankton biomass increase. Projected warming characterized by an increase in sea surface temperature of 2.29 ± 0.05 °C leads to a reduction in zooplankton and phytoplankton biomasses of 11% and 6%, respectively. This suggests negative amplification of climate driven modifications of trophic level biomass through bottom-up control, leading to a reduced capacity of oceans to regulate climate through the biological carbon pump. Simulations suggest negative amplification is the dominant response across 47% of the ocean surface and prevails in the tropical oceans; whilst positive trophic amplification prevails in the Arctic and Antarctic oceans. Trophic attenuation is projected in temperate seas. Uncertainties in ocean plankton projections, associated to the use of single global and regional models, imply the need for caution when extending these considerations into higher trophic levels.

Marine biodiversity and ecosystem function relationships:the potential for practical monitoring applications

There is an increasing demand for environmental assessments of the marine environment to include ecosystem function. However, existing schemes are predominantly based on taxonomic (i.e. structural) measures of biodiversity. Biodiversity and Ecosystem Function (BEF) relationships are suggested to provide a mechanism for converting taxonomic information into surrogates of ecosystem function. This review assesses the evidence for marine BEF relationships and their potential to be used in practical monitoring applications (i.e. operationalized). Five key requirements were identified for the practical application of BEF relationships: (1) a complete understanding of strength, direction and prevalence of marine BEF relationships, (2) an understanding of which biological components are influential within specific BEF relationships, (3) the biodiversity of the selected biological components can be measured easily, (4) the ecological mechanisms that are the most important for generating marine BEF relationships, i.e. identity effects or complementarity, are known and (5) the proportion of the overall functional variance is explained by biodiversity, and hence BEF relationships, has been established. Numerous positive and some negative BEF relationships were found within the literature, although many reproduced poorly the natural species richness, trophic structures or multiple functions of real ecosystems (requirement 1). Null relationships were also reported. The consistency of the positive and negative relationships was often low that compromised the ability to generalize BEF relationships and confident application of BEF within marine monitoring. Equally, some biological components and functions have received little or no investigation. Expert judgement was used to attribute biological components using spatial extent, presence and functional rate criteria (requirement 2). This approach highlighted the main biological components contributing the most to specific ecosystem functions, and that many of the particularly influential components were found to have received the least amount of research attention. The need for biodiversity to be measureable (requirement 3) is possible for most biological components although difficult within the functionally important microbes. Identity effects underpinned most marine BEF relationships (requirement 4). As such, processes that translated structural biodiversity measures into functional diversity were found to generate better BEF relationships. The analysis of the contribution made by biodiversity, over abiotic influences, to the total expression of a particular ecosystem function was rarely measured or considered (requirement 5). Hence it is not possible to determine the overall importance of BEF relationships within the total ecosystem functioning observed. In the few studies where abiotic factors had been considered, it was clear that these modified BEF relationships and have their own direct influence on functional rate. Based on the five requirements, the information required for immediate ‘operationalization’ of BEF relationships within marine functional monitoring is lacking. However, the concept of BEF inclusion within practical monitoring applications, supported by ecological modelling, shows promise for providing surrogate indicators of functioning.

Spatio-temporal Bayesian network models with latent variables for revealing trophic dynamics and functional networks in fisheries ecology

Ecosystems consist of complex dynamic interactions among species and the environment, the understanding of which has implications for predicting the environmental response to changes in climate and biodiversity. However, with the recent adoption of more explorative tools, like Bayesian networks, in predictive ecology, few assumptions can be made about the data and complex, spatially varying interactions can be recovered from collected field data. In this study, we compare Bayesian network modelling approaches accounting for latent effects to reveal species dynamics for 7 geographically and temporally varied areas within the North Sea. We also apply structure learning techniques to identify functional relationships such as prey–predator between trophic groups of species that vary across space and time. We examine if the use of a general hidden variable can reflect overall changes in the trophic dynamics of each spatial system and whether the inclusion of a specific hidden variable can model unmeasured group of species. The general hidden variable appears to capture changes in the variance of different groups of species biomass. Models that include both general and specific hidden variables resulted in identifying similarity with the underlying food web dynamics and modelling spatial unmeasured effect. We predict the biomass of the trophic groups and find that predictive accuracy varies with the models' features and across the different spatial areas thus proposing a model that allows for spatial autocorrelation and two hidden variables. Our proposed model was able to produce novel insights on this ecosystem's dynamics and ecological interactions mainly because we account for the heterogeneous nature of the driving factors within each area and their changes over time. Our findings demonstrate that accounting for additional sources of variation, by combining structure learning from data and experts' knowledge in the model architecture, has the potential for gaining deeper insights into the structure and stability of ecosystems. Finally, we were able to discover meaningful functional networks that were spatially and temporally differentiated with the particular mechanisms varying from trophic associations through interactions with climate and commercial fisheries.

Body mass-abundance relationships are robust to cascading effects in marine food webs

Body mass has been shown to scale negatively with abundance in a wide range of habitats and ecosystems. It is believed that this relationship has important consequences for the distribution and maintenance of energy in natural communities. Some studies have shown that the relationship between body mass and abundance may be robust to major food web perturbations, fuelling the belief that natural processes may preserve the slope of this relationship and the associated cycling of energy and nutrients. Here, we use data from a long-term experimental food web manipulation to examine this issue in a semi-natural environment. Similar communities were developed in large experimental mesocosms over a six month period. Some of the mesocosms were then subjected to species removals, based on the mean strength of their trophic interactions in the communities. In treatments where the strongest interactors were removed, a community-level trophic cascade occurred. The biomass density of invertebrates increased dramatically in these communities, which led to a suppression of primary production. In spite of these widespread changes in ecosystem functioning, the slope of the relationship between body mass and abundance remained unchanged. This was the case whether average species body mass and abundance or individual organism size spectra were considered. An examination of changes in species composition before and after the experimental manipulations revealed an important mechanism for maintaining the body mass-abundance relationship. The manipulated communities all had a higher species turnover than the intact communities, with the highest turnover in communities that experienced cascading effects. As some species increased in body mass and abundance, new species filled the available size-abundance niches that were created. This maintained the overall body mass-abundance relationship and provided a stabilising structure to these experimental communities.

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.

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

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

Species Richness-Environment Relationships of European Arthropods at Two Spatial Grains: Habitats and Countries

We study how species richness of arthropods relates to theories concerning net primary productivity, ambient energy, water-energy dynamics and spatial environmental heterogeneity. We use two datasets of arthropod richness with similar spatial extents (Scandinavia to Mediterranean), but contrasting spatial grain (local habitat and country). Samples of ground-dwelling spiders, beetles, bugs and ants were collected from 32 paired habitats at 16 locations across Europe. Species richness of these taxonomic groups was also determined for 25 European countries based on the Fauna Europaea database. We tested effects of net primary productivity (NPP), annual mean temperature (T), annual rainfall (R) and potential evapotranspiration of the coldest month (PETmin) on species richness and turnover. Spatial environmental heterogeneity within countries was considered by including the ranges of NPP, T, R and PETmin. At the local habitat grain, relationships between species richness and environmental variables differed strongly between taxa and trophic groups. However, species turnover across locations was strongly correlated with differences in T. At the country grain, species richness was significantly correlated with environmental variables from all four theories. In particular, species richness within countries increased strongly with spatial heterogeneity in T. The importance of spatial heterogeneity in T for both species turnover across locations and for species richness within countries suggests that the temperature niche is an important determinant of arthropod diversity. We suggest that, unless climatic heterogeneity is constant across sampling units, coarse-grained studies should always account for environmental heterogeneity as a predictor of arthropod species richness, just as studies with variable area of sampling units routinely consider area.

Indirect and direct controls of macroinvertebrates and small fish by abiotic factors and trophic interactions in the Florida Everglades

1. The roles of nutrients, disturbance and predation in regulating consumer densities have long been of interest, but their indirect effects have rarely been quantified in wetland ecosystems. The Florida Everglades contains gradients of hydrological disturbance (marsh drying) and nutrient enrichment (phosphorus), often correlated with densities of macroinvertebrate infauna (macroinvertebrates inhabiting periphyton), small fish and larger invertebrates, such as snails, grass shrimp, insects and crayfish. However, most causal relationships have yet to be quantified. 2.  We sampled periphyton (content and community structure) and consumer (small omnivores, carnivores and herbivores, and infaunal macroinvertebrates inhabiting periphyton) density at 28 sites spanning a range of hydrological and nutrient conditions and compared our data to seven a priori structural equation models. 3.  The best model included bottom-up and top-down effects among trophic groups and supported top-down control of infauna by omnivores and predators that cascaded to periphyton biomass. The next best model included bottom-up paths only and allowed direct effects of periphyton on omnivore density. Both models suggested a positive relationship between small herbivores and small omnivores, indicating that predation was unable to limit herbivore numbers. Total effects of time following flooding were negative for all three consumer groups even when both preferred models suggested positive direct effects for some groups. Total effects of nutrient levels (phosphorus) were positive for consumers and generally larger than those of hydrological disturbance and were mediated by changes in periphyton content. 4.  Our findings provide quantitative support for indirect effects of nutrient enrichment on consumers, and the importance of both algal community structure and periphyton biomass to Everglades food webs. Evidence for top-down control of infauna by omnivores was noted, though without substantially greater support than a competing bottom-up-only model.

Spatial Relationships among Hydroacoustic, Hydrographic and Top Predator Patterns: Cetacean Distributions in the Mid-Atlantic Bight

Effective conservation and management of top predators requires a comprehensive understanding of their distributions and of the underlying biological and physical processes that affect these distributions. The Mid-Atlantic Bight shelf break system is a dynamic and productive region where at least 32 species of cetaceans have been recorded through various systematic and opportunistic marine mammal surveys from the 1970s through 2012. My dissertation characterizes the spatial distribution and habitat of cetaceans in the Mid-Atlantic Bight shelf break system by utilizing marine mammal line-transect survey data, synoptic multi-frequency active acoustic data, and fine-scale hydrographic data collected during the 2011 summer Atlantic Marine Assessment Program for Protected Species (AMAPPS) survey. Although studies describing cetacean habitat and distributions have been previously conducted in the Mid-Atlantic Bight, my research specifically focuses on the shelf break region to elucidate both the physical and biological processes that influence cetacean distribution patterns within this cetacean hotspot.

In Chapter One I review biologically important areas for cetaceans in the Atlantic waters of the United States. I describe the study area, the shelf break region of the Mid-Atlantic Bight, in terms of the general oceanography, productivity and biodiversity. According to recent habitat-based cetacean density models, the shelf break region is an area of high cetacean abundance and density, yet little research is directed at understanding the mechanisms that establish this region as a cetacean hotspot.

In Chapter Two I present the basic physical principles of sound in water and describe the methodology used to categorize opportunistically collected multi-frequency active acoustic data using frequency responses techniques. Frequency response classification methods are usually employed in conjunction with net-tow data, but the logistics of the 2011 AMAPPS survey did not allow for appropriate net-tow data to be collected. Biologically meaningful information can be extracted from acoustic scattering regions by comparing the frequency response curves of acoustic regions to theoretical curves of known scattering models. Using the five frequencies on the EK60 system (18, 38, 70, 120, and 200 kHz), three categories of scatterers were defined: fish-like (with swim bladder), nekton-like (e.g., euphausiids), and plankton-like (e.g., copepods). I also employed a multi-frequency acoustic categorization method using three frequencies (18, 38, and 120 kHz) that has been used in the Gulf of Maine and Georges Bank which is based the presence or absence of volume backscatter above a threshold. This method is more objective than the comparison of frequency response curves because it uses an established backscatter value for the threshold. By removing all data below the threshold, only strong scattering information is retained.

In Chapter Three I analyze the distribution of the categorized acoustic regions of interest during the daytime cross shelf transects. Over all transects, plankton-like acoustic regions of interest were detected most frequently, followed by fish-like acoustic regions and then nekton-like acoustic regions. Plankton-like detections were the only significantly different acoustic detections per kilometer, although nekton-like detections were only slightly not significant. Using the threshold categorization method by Jech and Michaels (2006) provides a more conservative and discrete detection of acoustic scatterers and allows me to retrieve backscatter values along transects in areas that have been categorized. This provides continuous data values that can be integrated at discrete spatial increments for wavelet analysis. Wavelet analysis indicates significant spatial scales of interest for fish-like and nekton-like acoustic backscatter range from one to four kilometers and vary among transects.

In Chapter Four I analyze the fine scale distribution of cetaceans in the shelf break system of the Mid-Atlantic Bight using corrected sightings per trackline region, classification trees, multidimensional scaling, and random forest analysis. I describe habitat for common dolphins, Risso’s dolphins and sperm whales. From the distribution of cetacean sightings, patterns of habitat start to emerge: within the shelf break region of the Mid-Atlantic Bight, common dolphins were sighted more prevalently over the shelf while sperm whales were more frequently found in the deep waters offshore and Risso’s dolphins were most prevalent at the shelf break. Multidimensional scaling presents clear environmental separation among common dolphins and Risso’s dolphins and sperm whales. The sperm whale random forest habitat model had the lowest misclassification error (0.30) and the Risso’s dolphin random forest habitat model had the greatest misclassification error (0.37). Shallow water depth (less than 148 meters) was the primary variable selected in the classification model for common dolphin habitat. Distance to surface density fronts and surface temperature fronts were the primary variables selected in the classification models to describe Risso’s dolphin habitat and sperm whale habitat respectively. When mapped back into geographic space, these three cetacean species occupy different fine-scale habitats within the dynamic Mid-Atlantic Bight shelf break system.

In Chapter Five I present a summary of the previous chapters and present potential analytical steps to address ecological questions pertaining the dynamic shelf break region. Taken together, the results of my dissertation demonstrate the use of opportunistically collected data in ecosystem studies; emphasize the need to incorporate middle trophic level data and oceanographic features into cetacean habitat models; and emphasize the importance of developing more mechanistic understanding of dynamic ecosystems.