Spatial synchrony in coral reef fish populations and the influence of climate

Can juvenile supply predict future abundance of large-bodied reef fishes?

The extent to which juvenile abundance can predict future populations of lethrinids at Ningaloo Reef was assessed using size frequency data collected over 13 consecutive years. Annual abundance of juvenile lethrinids (<5 cm TL) was highest in northern Ningaloo during La Niña years, when seawater is warmer and oceanic currents stronger. Juvenile lethrinid abundance explained 35% of the variance in 1-2 year-old Lethrinus nebulosus abundance the following year, a steeper relationship in the north suggesting greater survival of juveniles. Juvenile lethrinid abundance was also positively correlated to abundance of 1-2 year-old L. atkinsoni in the southern region of Ningaloo. Abundance of juvenile lethrinids were however poor predictors of L. nebulosus and L. atkinsoni older than 2 years of age. Post settlement processes likely weaken the link between juvenile supply and abundance of lethrinids >2 years old making it difficult to accurately quantify the overall size of future lethrinid populations.

Dispersal synchronizes giant kelp forests

Spatial synchrony is the tendency for population fluctuations to be correlated among different locations. This phenomenon is a ubiquitous feature of population dynamics and is important for ecosystem stability, but several aspects of synchrony remain unresolved. In particular, the extent to which any particular mechanism, such as dispersal, contributes to observed synchrony in natural populations has been difficult to determine. To address this gap, we leveraged recent methodological improvements to determine how dispersal structures synchrony in giant kelp (Macrocystis pyrifera), a global marine foundation species that has served as a useful system for understanding synchrony. We quantified population synchrony and fecundity with satellite imagery across 11 years and 880 km of coastline in southern California, USA, and estimated propagule dispersal probabilities using a high-resolution ocean circulation model. Using matrix regression models that control for the influence of geographic distance, resources (seawater nitrate), and disturbance (destructive waves), we discovered that dispersal was an important driver of synchrony. Our findings were robust to assumptions about propagule mortality during dispersal and consistent between two metrics of dispersal: (1) the individual probability of dispersal and (2) estimates of demographic connectivity that incorporate fecundity (the number of propagules dispersing). We also found that dispersal and environmental conditions resulted in geographic clusters with distinct patterns of synchrony. This study is among the few to statistically associate synchrony with dispersal in a natural population and the first to do so in a marine organism. The synchronizing effects of dispersal and environmental conditions on foundation species, such as giant kelp, likely have cascading effects on the spatial stability of biodiversity and ecosystem function.

The El Niño Southern Oscillation drives multidirectional inter-reef larval connectivity in the Great Barrier Reef

The El Niño Southern Oscillation (ENSO) is the strongest source of interannual global climate variability, and extreme ENSO events are projected to increase in frequency under climate change. Interannual variability in the Coral Sea circulation has been associated with ENSO, although uncertainty remains regarding ENSO's influence on hydrodynamics and larval dispersal in the adjacent Great Barrier Reef (GBR). We investigated larval connectivity during ENSO events from 2010 to 2017 throughout the GBR, based on biophysical modelling of a widespread predatory reef fish, Lutjanus carponotatus. Our results indicate a well-connected system over the study period with high interannual variability in inter-reef connectivity associated with ENSO. Larval connectivity patterns were highly correlated to variations in the Southern Oscillation Index (SOI). During El Niño conditions and periods of weak SOI, larval dispersal patterns were predominantly poleward in the central and southern regions, reversing to a predominant equatorward flow during very strong SOI and extreme La Niña conditions. These ENSO-linked connectivity patterns were associated with positive connectivity anomalies among reefs. Our findings identify ENSO as an important source of variation in larval dispersal and connectivity patterns in the GBR, which can influence the stability of population dynamics and patterns of biodiversity in the region.

The long and the short of it: Mechanisms of synchronous and compensatory dynamics across temporal scales

Synchronous dynamics (fluctuations that occur in unison) are universal phenomena with widespread implications for ecological stability. Synchronous dynamics can amplify the destabilizing effect of environmental variability on ecosystem functions such as productivity, whereas the inverse, compensatory dynamics, can stabilize function. Here we combine simulation and empirical analyses to elucidate mechanisms that underlie patterns of synchronous versus compensatory dynamics. In both simulated and empirical communities, we show that synchronous and compensatory dynamics are not mutually exclusive but instead can vary by timescale. Our simulations identify multiple mechanisms that can generate timescale‐specific patterns, including different environmental drivers, diverse life histories, dispersal, and non‐stationary dynamics. We find that traditional metrics for quantifying synchronous dynamics are often biased toward long‐term drivers and may miss the importance of short‐term drivers. Our findings indicate key mechanisms to consider when assessing synchronous versus compensatory dynamics and our approach provides a pathway for disentangling these dynamics in natural systems.

Global hotspots of coherent marine fishery catches

Although different fisheries can be tightly linked to each other by human and ecosystem processes, they are often managed independently. Synchronous fluctuations among fish populations or fishery catches can destabilize ecosystems and economies, respectively, but the degree of synchrony around the world remains unclear. We analyzed 1,092 marine fisheries catch time series over 60 yr to test for the presence of coherence, a form of synchrony that allows for phase-lagged relationships. We found that nearly every fishery was coherent with at least one other fishery catch time series globally and that coherence was strongest in the northeast Atlantic, western central Pacific, and eastern Indian Ocean. Analysis of fish biomass and fishing mortality time series from these hotspots revealed that coherence in biomass or fishing mortality were both possible, though biomass coherence was more common. Most of these relationships were synchronous with no time lags, and across catches in all regions, synchrony was a better predictor of regional catch portfolio effects than catch diversity. Regions with higher synchrony had lower stability in aggregate fishery catches, which can have negative consequences for food security and economic wealth.

Elevated temperature and CO2 have positive effects on the growth and survival of larval Australasian snapper

Rising water temperature and increased uptake of CO₂ by the ocean are predicted to have widespread impacts on marine species. However, the effects are likely to vary, depending on a species' sensitivity and the geographical location of the population. Here, we investigated the potential effects of elevated temperature and pCO₂ on larval growth and survival in a New Zealand population of the Australasian snapper, Chrysophyrs auratus. Eggs and larvae were reared in a fully cross-factored experiment (18 °C and 22 °C/pCO₂ 440 and 1040 μatm) to 16 days post hatch (dph). Morphologies at 1 dph and 16 dph were significantly affected by temperature, but not CO₂. At 1dph, larvae at 22 °C were longer (7%) and had larger muscle depth at vent (14%), but had reduced yolk (65%) and oil globule size (16%). Reduced yolk reserves in recently hatched larvae suggests higher metabolic demands in warmer water. At 16 dph, larvae at elevated temperature were longer (12%) and muscle depth at vent was larger (64%). Conversely, survival was primarily affected by CO₂ rather than temperature. Survivorship at 1 dph and 16 dph was 24% and 54% higher, respectively, under elevated CO₂ compared with ambient conditions. Elevated temperature increased survival (24%) at 1 dph, but not at 16 dph. These results suggest that projected climate change scenarios may have an overall positive effect on early life history growth and survival in this population of C. auratus. This could benefit recruitment success, but needs to be weighed against negative effects of elevated CO₂ on metabolic rates and swimming performance observed in other studies on the same population.

Synchrony is more than its top-down and climatic parts: interacting Moran effects on phytoplankton in British seas

Large-scale spatial synchrony is ubiquitous in ecology. We examined 56 years of data representing chlorophyll density in 26 areas in British seas monitored by the Continuous Plankton Recorder survey. We used wavelet methods to disaggregate synchronous fluctuations by timescale and determine that drivers of synchrony include both biotic and abiotic variables. We tested these drivers for statistical significance by comparison with spatially synchronous surrogate data. Identification of causes of synchrony is distinct from, and goes beyond, determining drivers of local population dynamics. We generated timescale-specific models, accounting for 61% of long-timescale (> 4yrs) synchrony in a chlorophyll density index, but only 3% of observed short-timescale (< 4yrs) synchrony. Thus synchrony and its causes are timescale-specific. The dominant source of long-timescale chlorophyll synchrony was closely related to sea surface temperature, through a climatic Moran effect, though likely via complex oceanographic mechanisms. The top-down action of Calanus finmarchicus predation enhances this environmental synchronising mechanism and interacts with it non-additively to produce more long-timescale synchrony than top-down and climatic drivers would produce independently. Our principal result is therefore a demonstration of interaction effects between Moran drivers of synchrony, a new mechanism for synchrony that may influence many ecosystems at large spatial scales.

The size of the annual bloom in phytoplankton can vary similarly from year to year in different parts of the same oceanic region, a phenomenon called spatial synchrony. The growth of phytoplankton near the ocean surface is the foundation of marine food webs, which include numerous commercially exploited species. And spatial synchrony in phytoplankton abundance time series can have consequences for the total production of marine ecosystems. Therefore we studied the spatial synchrony of fluctuations in green phytoplankton abundance in 26 areas in seas around the British Isles. Variation and synchrony can occur differently on long and short timescales. We used a novel wavelet-based approach to examine long- and short-timescale fluctuations separately, and we thereby show that slow synchronous fluctuations in phytoplankton can be explained by the effects of slow synchronous fluctuations in sea surface temperature and related oceanographic phenomena, and by the effects of synchronous fluctuations in a zooplankton predator. Crucially, these drivers reinforce one another in a super-additive way, the interaction constituting a new mechanism of synchrony. Future changes in the climate or changes in predation are likely to influence phytoplankton synchrony via this mechanism and hence may influence the aggregate productivity of British seas.

Taylor's power law captures the effects of environmental variability on community structure: An example from fishes in the North Sea

Taylor's power law (TPL) describes the relationship between the mean and variance in abundance of populations, with the power law exponent considered a measure of aggregation. However, the usefulness of TPL exponents as an ecological metric has been questioned, largely due to its apparent ubiquity in various complex systems. The aim of this study was to test whether TPL exponents vary systematically with potential drivers of animal aggregation in time and space and therefore capture useful ecological information of the system of interest. We derived community TPL exponents from a long-term, standardised and spatially dense data series of abundance and body size data for a strongly size-structured fish community in the North Sea. We then compared TPL exponents between regions of contrasting environmental characteristics. We find that, in general, TPL exponents vary more than expected under random conditions in the North Sea for size-based populations compared to communities considered by species. Further, size-based temporal TPL exponents are systematically higher (implying more temporally aggregated distributions) along hydrographic boundaries. Time series of size-based spatial TPL exponents also differ between hydrographically distinct basins. These findings support the notion that TPL exponents contain ecological information, capturing community spatio-temporal dynamics as influenced by external drivers.

Climatic forcing and larval dispersal capabilities shape the replenishment of fishes and their habitat-forming biota on a tropical coral reef

Fluctuations in marine populations often relate to the supply of recruits by oceanic currents. Variation in these currents is typically driven by large-scale changes in climate, in particular ENSO (El Nino Southern Oscillation). The dependence on large-scale climatic changes may, however, be modified by early life history traits of marine taxa. Based on eight years of annual surveys, along 150 km of coastline, we examined how ENSO influenced abundance of juvenile fish, coral spat, and canopy-forming macroalgae. We then investigated what traits make populations of some fish families more reliant on the ENSO relationship than others. Abundance of juvenile fish and coral recruits was generally positively correlated with the Southern Oscillation Index (SOI), higher densities recorded during La Niña years, when the ENSO-influenced Leeuwin Current is stronger and sea surface temperature higher. The relationship is typically positive and stronger among fish families with shorter pelagic larval durations and stronger swimming abilities. The relationship is also stronger at sites on the coral back reef, although the strongest of all relationships were among the lethrinids (r = .9), siganids (r = .9), and mullids (r = .8), which recruit to macroalgal meadows in the lagoon. ENSO effects on habitat seem to moderate SOI-juvenile abundance relationship. Macroalgal canopies are higher during La Niña years, providing more favorable habitat for juvenile fish and strengthening the SOI effect on juvenile abundance. Conversely, loss of coral following a La Niña-related heat wave may have compromised postsettlement survival of coral dependent species, weakening the influence of SOI on their abundance. This assessment of ENSO effects on tropical fish and habitat-forming biota and how it is mediated by functional ecology improves our ability to predict and manage changes in the replenishment of marine populations.

Environmental triggers for primary outbreaks of crown-of-thorns starfish on the Great Barrier Reef, Australia

In this paper, we postulate a unique environmental triggering sequence for primary outbreaks of crown-of-thorns starfish (COTS, Acanthaster planci) on the central Great Barrier Reef (GBR, Australia). Notably, we extend the previous terrestrial runoff hypothesis, viz. nutrient-enriched terrestrial runoff → elevated phytoplankton 'bloom' concentrations → enhanced COTS larval survival, to include the additional importance of strong larvae retention around reefs or within reef groups (clusters) that share enhanced phytoplankton concentrations. For the central GBR, this scenario is shown to occur when El Niño-Southern Oscillation (ENSO) linked hydrodynamic conditions cause the 'regional' larval connectivity network to fragment into smaller 'local' reef clusters due to low ocean current velocities. As inter-annual variations in hydrodynamic circulation patterns are not amenable to direct management intervention, the ability to reduce the future frequency of COTS outbreaks on the central GBR is shown to be contingent on reducing terrestrial bioavailable nutrient loads ~20-40%.

Phase resetting for a network of oscillators via phase response curve approach

The problem of phase regulation for a population of oscillating systems is considered. The proposed control strategy is based on a phase response curve (PRC) model of an oscillator (the first-order reduced model obtained for linearized system and inputs with infinitesimal amplitude). It is proven that the control provides phase resetting for the original nonlinear system. Next, the problem of phase resetting for a network of oscillators is considered when applying a common control input. Performance of the obtained solutions is demonstrated via computer simulation for three different models of circadian/neural oscillators.

Habitat specialization in tropical continental shelf demersal fish assemblages

The implications of shallow water impacts such as fishing and climate change on fish assemblages are generally considered in isolation from the distribution and abundance of these fish assemblages in adjacent deeper waters. We investigate the abundance and length of demersal fish assemblages across a section of tropical continental shelf at Ningaloo Reef, Western Australia, to identify fish and fish habitat relationships across steep gradients in depth and in different benthic habitat types. The assemblage composition of demersal fish were assessed from baited remote underwater stereo-video samples (n = 304) collected from 16 depth and habitat combinations. Samples were collected across a depth range poorly represented in the literature from the fringing reef lagoon (1–10 m depth), down the fore reef slope to the reef base (10–30 m depth) then across the adjacent continental shelf (30–110 m depth). Multivariate analyses showed that there were distinctive fish assemblages and different sized fish were associated with each habitat/depth category. Species richness, MaxN and diversity declined with depth, while average length and trophic level increased. The assemblage structure, diversity, size and trophic structure of demersal fishes changes from shallow inshore habitats to deeper water habitats. More habitat specialists (unique species per habitat/depth category) were associated with the reef slope and reef base than other habitats, but offshore sponge-dominated habitats and inshore coral-dominated reef also supported unique species. This suggests that marine protected areas in shallow coral-dominated reef habitats may not adequately protect those species whose depth distribution extends beyond shallow habitats, or other significant elements of demersal fish biodiversity. The ontogenetic habitat partitioning which is characteristic of many species, suggests that to maintain entire species life histories it is necessary to protect corridors of connected habitats through which fish can migrate.

Effects of ENSO-linked climate and vegetation on population dynamics of sympatric rodent species in semiarid grasslands of Inner Mongolia, China

El Niño Southern Oscillation (ENSO) linked climate has been known to be associated with several rodent species, but its effects on rodent community at both spatial and temporal scales are not well studied. In this study, we investigated the possible causal chain relating ENSO, precipitation, temperature, and vegetation index (normalized difference vegetation index, NDVI) to rodent abundance for 14 sympatric rodent species in 21 counties of semiarid grasslands in Inner Mongolia, China, from 1982 to 2006. We found that both precipitation and temperature showed a generally direct positive effect on rodent abundance in many species in the current year, but indirect effects that operate through NDVI in the current or following year could have a reverse effect on abundance. We described one ENSO-linked precipitation bottom-up chain and three ENSO-linked temperature bottom-up chains. These observed bottom-up links reveal that in El Niño years, or 1 year after La Niña years, or 2 years after El Niño years, ENSO-driven climate or vegetation factors tend to increase population abundances of many sympatric rodent species in this region. We also found time-lag effects and the life-history strategy (i.e., functional groups of hibernating behavior, activity rhythm, or food habits) also contribute to the observed complicated effects of SOI on precipitation, temperature, NDVI, and ultimately rodent abundance.

Regression analysis of spatial data

Many of the most interesting questions ecologists ask lead to analyses of spatial data. Yet, perhaps confused by the large number of statistical models and fitting methods available, many ecologists seem to believe this is best left to specialists. Here, we describe the issues that need consideration when analysing spatial data and illustrate these using simulation studies. Our comparative analysis involves using methods including generalized least squares, spatial filters, wavelet revised models, conditional autoregressive models and generalized additive mixed models to estimate regression coefficients from synthetic but realistic data sets, including some which violate standard regression assumptions. We assess the performance of each method using two measures and using statistical error rates for model selection. Methods that performed well included generalized least squares family of models and a Bayesian implementation of the conditional auto-regressive model. Ordinary least squares also performed adequately in the absence of model selection, but had poorly controlled Type I error rates and so did not show the improvements in performance under model selection when using the above methods. Removing large-scale spatial trends in the response led to poor performance. These are empirical results; hence extrapolation of these findings to other situations should be performed cautiously. Nevertheless, our simulation-based approach provides much stronger evidence for comparative analysis than assessments based on single or small numbers of data sets, and should be considered a necessary foundation for statements of this type in future.

Crucial knowledge gaps in current understanding of climate change impacts on coral reef fishes

Expert opinion was canvassed to identify crucial knowledge gaps in current understanding of climate change impacts on coral reef fishes. Scientists that had published three or more papers on the effects of climate and environmental factors on reef fishes were invited to submit five questions that, if addressed, would improve our understanding of climate change effects on coral reef fishes. Thirty-three scientists provided 155 questions, and 32 scientists scored these questions in terms of: (i) identifying a knowledge gap, (ii) achievability, (iii) applicability to a broad spectrum of species and reef habitats, and (iv) priority. Forty-two per cent of the questions related to habitat associations and community dynamics of fish, reflecting the established effects and immediate concern relating to climate-induced coral loss and habitat degradation. However, there were also questions on fish demographics, physiology, behaviour and management, all of which could be potentially affected by climate change. Irrespective of their individual expertise and background, scientists scored questions from different topics similarly, suggesting limited bias and recognition of a need for greater interdisciplinary and collaborative research. Presented here are the 53 highest-scoring unique questions. These questions should act as a guide for future research, providing a basis for better assessment and management of climate change impacts on coral reefs and associated fish communities.

Synchrony and stability of food webs in metacommunities

Synchrony has fundamental but conflicting implications for the persistence and stability of food webs at local and regional scales. In a constant environment, compensatory dynamics between species can maintain food web stability, but factors that synchronize population fluctuations within and between communities are expected to be destabilizing. We studied the dynamics of a food web in a metacommunity to determine how environmental variability and dispersal affect stability by altering compensatory dynamics and average species abundance. When dispersal rate is high, weak correlated environmental fluctuations promote food web stability by reducing the amplitude of compensatory dynamics. However, when dispersal rate is low, weak environmental fluctuations reduce food web stability by inducing intraspecific synchrony across communities. Irrespective of dispersal rate, strong environmental fluctuations disrupt compensatory dynamics and decrease stability by inducing intermittent correlated fluctuations between consumers in local food webs, which reduce both total consumer abundance and predator abundance. Strong correlated environmental fluctuations lead to (i) spatially asynchronous and highly correlated local consumer dynamics when dispersal is low and (ii) spatially synchronous but intermediate local consumer correlation when dispersal is high. By controlling intraspecific synchrony, dispersal mediates the capacity of strong environmental fluctuations to disrupt compensatory dynamics at both local and metacommunity scales.