Seagrass ecosystem multifunctionality under the rise of a flagship marine megaherbivore.

Large grazers (megaherbivores) have a profound impact on ecosystem functioning. However, how ecosystem multifunctionality is affected by changes in megaherbivore populations remains poorly understood. Understanding the total impact on ecosystem multifunctionality requires an integrative ecosystem approach, which is especially challenging to obtain in marine systems. We assessed the effects of experimentally simulated grazing intensity scenarios on ecosystem functions and multifunctionality in a tropical Caribbean seagrass ecosystem. As a model, we selected a key marine megaherbivore, the green turtle, whose ecological role is rapidly unfolding in numerous foraging areas where populations are recovering through conservation after centuries of decline, with an increase in recorded overgrazing episodes. To quantify the effects, we employed a novel integrated index of seagrass ecosystem multifunctionality based upon multiple, well-recognized measures of seagrass ecosystem functions that reflect ecosystem services. Experiments revealed that intermediate turtle grazing resulted in the highest rates of nutrient cycling and carbon storage, while sediment stabilization, decomposition rates, epifauna richness, and fish biomass are highest in the absence of turtle grazing. In contrast, intense grazing resulted in disproportionally large effects on ecosystem functions and a collapse of multifunctionality. These results imply that (i) the return of a megaherbivore can exert strong effects on coastal ecosystem functions and multifunctionality, (ii) conservation efforts that are skewed toward megaherbivores, but ignore their key drivers like predators or habitat, will likely result in overgrazing-induced loss of multifunctionality, and (iii) the multifunctionality index shows great potential as a quantitative tool to assess ecosystem performance. Considerable and rapid alterations in megaherbivore abundance (both through extinction and conservation) cause an imbalance in ecosystem functioning and substantially alter or even compromise ecosystem services that help to negate global change effects. An integrative ecosystem approach in environmental management is urgently required to protect and enhance ecosystem multifunctionality.

Demographic changes in Pleistocene sea turtles were driven by past sea level fluctuations affecting feeding habitat availability.

Pleistocene environmental changes are generally assumed to have dramatically affected species' demography via changes in habitat availability, but this is challenging to investigate due to our limited knowledge of how Pleistocene ecosystems changed through time. Here, we tracked changes in shallow marine habitat availability resulting from Pleistocene sea level fluctuations throughout the last glacial cycle (120-14 thousand years ago; kya) and assessed correlations with past changes in genetic diversity inferred from genome-wide SNPs, obtained via ddRAD sequencing, in Caribbean hawksbill turtles, which feed in coral reefs commonly found in shallow tropical waters. We found sea level regression resulted in an average 75% reduction in shallow marine habitat availability during the last glacial cycle. Changes in shallow marine habitat availability correlated strongly with past changes in hawksbill turtle genetic diversity, which gradually declined to ~1/4th of present-day levels during the Last Glacial Maximum (LGM; 26-19 kya). Shallow marine habitat availability and genetic diversity rapidly increased after the LGM, signifying a population expansion in response to warming environmental conditions. Our results suggest a positive correlation between Pleistocene environmental changes, habitat availability and species' demography, and that demographic changes in hawksbill turtles were potentially driven by feeding habitat availability. However, we also identified challenges associated with disentangling the potential environmental drivers of past demographic changes, which highlights the need for integrative approaches. Our conclusions underline the role of habitat availability on species' demography and biodiversity, and that the consequences of ongoing habitat loss should not be underestimated.

The population genomic structure of green turtles (Chelonia mydas) suggests a warm-water corridor for tropical marine fauna between the Atlantic and Indian oceans during the last interglacial.

The occasional westward transport of warm water of the Agulhas Current, "Agulhas leakage", around southern Africa has been suggested to facilitate tropical marine connectivity between the Atlantic and Indian oceans, but the "Agulhas leakage" hypothesis does not explain the signatures of eastward gene flow observed in many tropical marine fauna. We investigated an alternative hypothesis: the establishment of a warm-water corridor during comparatively warm interglacial periods. The "warm-water corridor" hypothesis was investigated by studying the population genomic structure of Atlantic and Southwest Indian Ocean green turtles (N = 27) using 12,035 genome-wide single nucleotide polymorphisms (SNPs) obtained via ddRAD sequencing. Model-based and multivariate clustering suggested a hierarchical population structure with two main Atlantic and Southwest Indian Ocean clusters, and a Caribbean and East Atlantic sub-cluster nested within the Atlantic cluster. Coalescent-based model selection supported a model where Southwest Indian Ocean and Caribbean populations diverged from the East Atlantic population during the transition from the last interglacial period (130-115 thousand years ago; kya) to the last glacial period (115-90 kya). The onset of the last glaciation appeared to isolate Atlantic and Southwest Indian Ocean green turtles into three refugia, which subsequently came into secondary contact in the Caribbean and Southwest Indian Ocean when global temperatures increased after the Last Glacial Maximum. Our findings support the establishment of a warm-water corridor facilitating tropical marine connectivity between the Atlantic and Southwest Indian Ocean during warm interglacials.

Ecological regime shift drives declining growth rates of sea turtles throughout the West Atlantic.

Somatic growth is an integrated, individual-based response to environmental conditions, especially in ectotherms. Growth dynamics of large, mobile animals are particularly useful as bio-indicators of environmental change at regional scales. We assembled growth rate data from throughout the West Atlantic for green turtles, Chelonia mydas, which are long-lived, highly migratory, primarily herbivorous mega-consumers that may migrate over hundreds to thousands of kilometers. Our dataset, the largest ever compiled for sea turtles, has 9690 growth increments from 30 sites from Bermuda to Uruguay from 1973 to 2015. Using generalized additive mixed models, we evaluated covariates that could affect growth rates; body size, diet, and year have significant effects on growth. Growth increases in early years until 1999, then declines by 26% to 2015. The temporal (year) effect is of particular interest because two carnivorous species of sea turtles-hawksbills, Eretmochelys imbricata, and loggerheads, Caretta caretta-exhibited similar significant declines in growth rates starting in 1997 in the West Atlantic, based on previous studies. These synchronous declines in productivity among three sea turtle species across a trophic spectrum provide strong evidence that an ecological regime shift (ERS) in the Atlantic is driving growth dynamics. The ERS resulted from a synergy of the 1997/1998 El Niño Southern Oscillation (ENSO)-the strongest on record-combined with an unprecedented warming rate over the last two to three decades. Further support is provided by the strong correlations between annualized mean growth rates of green turtles and both sea surface temperatures (SST) in the West Atlantic for years of declining growth rates (r = -.94) and the Multivariate ENSO Index (MEI) for all years (r = .74). Granger-causality analysis also supports the latter finding. We discuss multiple stressors that could reinforce and prolong the effect of the ERS. This study demonstrates the importance of region-wide collaborations.

Population recovery changes population composition at a major southern Caribbean juvenile developmental habitat for the green turtle, Chelonia mydas.

Understanding the population composition and dynamics of migratory megafauna at key developmental habitats is critical for conservation and management. The present study investigated whether differential recovery of Caribbean green turtle (Chelonia mydas) rookeries influenced population composition at a major juvenile feeding ground in the southern Caribbean (Lac Bay, Bonaire, Caribbean Netherlands) using genetic and demographic analyses. Genetic divergence indicated a strong temporal shift in population composition between 2006-2007 and 2015-2016 (ϕST = 0.101, P < 0.001). Juvenile recruitment (<75.0 cm straight carapace length; SCL) from the north-western Caribbean increased from 12% to 38% while recruitment from the eastern Caribbean region decreased from 46% to 20% between 2006-2007 and 2015-2016. Furthermore, the product of the population growth rate and adult female abundance was a significant predictor for population composition in 2015-2016. Our results may reflect early warning signals of declining reproductive output at eastern Caribbean rookeries, potential displacement effects of smaller rookeries by larger rookeries, and advocate for genetic monitoring as a useful method for monitoring trends in juvenile megafauna. Furthermore, these findings underline the need for adequate conservation of juvenile developmental habitats and a deeper understanding of the interactions between megafaunal population dynamics in different habitats.