Marine protected areas promote stability of reef fish communities under climate warming.

Protection from direct human impacts can safeguard marine life, yet ocean warming crosses marine protected area boundaries. Here, we test whether protection offers resilience to marine heatwaves from local to network scales. We examine 71,269 timeseries of population abundances for 2269 reef fish species surveyed in 357 protected versus 747 open sites worldwide. We quantify the stability of reef fish abundance from populations to metacommunities, considering responses of species and functional diversity including thermal affinity of different trophic groups. Overall, protection mitigates adverse effects of marine heatwaves on fish abundance, community stability, asynchronous fluctuations and functional richness. We find that local stability is positively related to distance from centers of high human density only in protected areas. We provide evidence that networks of protected areas have persistent reef fish communities in warming oceans by maintaining large populations and promoting stability at different levels of biological organization.

Grazer behaviour can regulate large-scale patterning of community states.

Ecosystem patterning can arise from environmental heterogeneity, biological feedbacks that produce multiple persistent ecological states, or their interaction. One source of feedbacks is density-dependent changes in behaviour that regulate species interactions. By fitting state-space models to large-scale (~500 km) surveys on temperate rocky reefs, we find that behavioural feedbacks best explain why kelp and urchin barrens form either reef-wide patches or local mosaics. Best-supported models in California include feedbacks where starvation intensifies grazing across entire reefs create reef-scale, alternatively stable kelp- and urchin-dominated states (32% of reefs). Best-fitting models in New Zealand include the feedback of urchins avoiding dense kelp stands that can increase abrasion and predation risk, which drives a transition from shallower urchin-dominated to deeper kelp-dominated zones, with patchiness at 3-8 m depths with intermediate wave stress. Connecting locally studied processes with region-wide data, we highlight how behaviour can explain community patterning and why some systems exhibit community-wide alternative stable states.

eDNA metabarcoding as a biomonitoring tool for marine protected areas.

Monitoring of marine protected areas (MPAs) is critical for marine ecosystem management, yet current protocols rely on SCUBA-based visual surveys that are costly and time consuming, limiting their scope and effectiveness. Environmental DNA (eDNA) metabarcoding is a promising alternative for marine ecosystem monitoring, but more direct comparisons to visual surveys are needed to understand the strengths and limitations of each approach. This study compares fish communities inside and outside the Scorpion State Marine Reserve off Santa Cruz Island, CA using eDNA metabarcoding and underwater visual census surveys. Results from eDNA captured 76% (19/25) of fish species and 95% (19/20) of fish genera observed during pairwise underwater visual census. Species missed by eDNA were due to the inability of MiFish 12S barcodes to differentiate species of rockfishes (Sebastes, n = 4) or low site occupancy rates of crevice-dwelling Lythrypnus gobies. However, eDNA detected an additional 23 fish species not recorded in paired visual surveys, but previously reported from prior visual surveys, highlighting the sensitivity of eDNA. Significant variation in eDNA signatures by location (50 m) and site (~1000 m) demonstrates the sensitivity of eDNA to address key questions such as community composition inside and outside MPAs. Results demonstrate the utility of eDNA metabarcoding for monitoring marine ecosystems, providing an important complementary tool to visual methods.

Improving the ability of a BACI design to detect impacts within a kelp-forest community.

Distinguishing between human impacts and natural variation in abundance remains difficult because most species exhibit complex patterns of variation in space and time. When ecological monitoring data are available, a before-after-control-impact (BACI) analysis can control natural spatial and temporal variation to better identify an impact and estimate its magnitude. However, populations with limited distributions and confounding spatial-temporal dynamics can violate core assumptions of BACI-type designs. In this study, we assessed how such properties affect the potential to identify impacts. Specifically, we quantified the conditions under which BACI analyses correctly (or incorrectly) identified simulated anthropogenic impacts in a spatially and temporally replicated data set of fish, macroalgal, and invertebrate species found on nearshore subtidal reefs in southern California, USA. We found BACI failed to assess very localized impacts, and had low power but high precision when assessing region-wide impacts. Power was highest for severe impacts of moderate spatial scale, and impacts were most easily detected in species with stable, widely distributed populations. Serial autocorrelation in the data greatly inflated false impact detection rates, and could be partly controlled for statistically, while spatial synchrony in dynamics had no consistent effect on power or false detection rates. Unfortunately, species that offer high power to detect real impacts were also more likely to detect impacts where none had occurred. However, considering power and false detection rates together can identify promising indicator species, and collectively analyzing data for similar species improved the net ability to assess impacts. These insights set expectations for the sizes and severities of impacts that BACI analyses can detect in real systems, point to the importance of serial autocorrelation (but not of spatial synchrony), and indicate how to choose the species, and groups of species, that can best identify impacts.

ranacapa: An R package and Shiny web app to explore environmental DNA data with exploratory statistics and interactive visualizations.

Environmental DNA (eDNA) metabarcoding is becoming a core tool in ecology and conservation biology, and is being used in a growing number of education, biodiversity monitoring, and public outreach programs in which professional research scientists engage community partners in primary research. Results from eDNA analyses can engage and educate natural resource managers, students, community scientists, and naturalists, but without significant training in bioinformatics, it can be difficult for this diverse audience to interact with eDNA results. Here we present the R package ranacapa, at the core of which is a Shiny web app that helps perform exploratory biodiversity analyses and visualizations of eDNA results. The app requires a taxonomy-by-sample matrix and a simple metadata file with descriptive information about each sample. The app enables users to explore the data with interactive figures and presents results from simple community ecology analyses. We demonstrate the value of ranacapa to two groups of community partners engaging with eDNA metabarcoding results.

Size, growth, and density data for shallow-water sea urchins from Mexico to the Aleutian Islands, Alaska, 1956-2016.

Size, growth, and density have been studied for North American Pacific coast sea urchins Strongylocentrotus purpuratus, S. droebachiensis, S. polyacanthus, Mesocentrotus (Strongylocentrotus) franciscanus, Lytechinus pictus, Centrostephanus coronatus, and Arbacia stellata by various workers at diverse sites and for varying lengths of time from 1956 to present. Numerous peer-reviewed publications have used some of these data but some data have appeared only in graduate theses or the gray literature. There also are data that have never appeared outside original data sheets. Motivation for studies has included fisheries management and environmental monitoring of sewer and power plant outfalls as well as changes associated with disease epidemics. Studies also have focused on kelp restoration, community effects of sea otters, basic sea urchin biology, and monitoring. The data sets presented here are a historical record of size, density, and growth for a common group of marine invertebrates in intertidal and nearshore environments that can be used to test hypotheses concerning future changes associated with fisheries practices, shifts of predator distributions, climate and ecosystem changes, and ocean acidification along the Pacific Coast of North America and islands of the north Pacific. No copyright restrictions apply. Please credit this paper when using the data.

Global patterns of kelp forest change over the past half-century.

Kelp forests (Order Laminariales) form key biogenic habitats in coastal regions of temperate and Arctic seas worldwide, providing ecosystem services valued in the range of billions of dollars annually. Although local evidence suggests that kelp forests are increasingly threatened by a variety of stressors, no comprehensive global analysis of change in kelp abundances currently exists. Here, we build and analyze a global database of kelp time series spanning the past half-century to assess regional and global trends in kelp abundances. We detected a high degree of geographic variation in trends, with regional variability in the direction and magnitude of change far exceeding a small global average decline (instantaneous rate of change = -0.018 y-1). Our analysis identified declines in 38% of ecoregions for which there are data (-0.015 to -0.18 y-1), increases in 27% of ecoregions (0.015 to 0.11 y-1), and no detectable change in 35% of ecoregions. These spatially variable trajectories reflected regional differences in the drivers of change, uncertainty in some regions owing to poor spatial and temporal data coverage, and the dynamic nature of kelp populations. We conclude that although global drivers could be affecting kelp forests at multiple scales, local stressors and regional variation in the effects of these drivers dominate kelp dynamics, in contrast to many other marine and terrestrial foundation species.

Global conservation outcomes depend on marine protected areas with five key features.

In line with global targets agreed under the Convention on Biological Diversity, the number of marine protected areas (MPAs) is increasing rapidly, yet socio-economic benefits generated by MPAs remain difficult to predict and under debate. MPAs often fail to reach their full potential as a consequence of factors such as illegal harvesting, regulations that legally allow detrimental harvesting, or emigration of animals outside boundaries because of continuous habitat or inadequate size of reserve. Here we show that the conservation benefits of 87 MPAs investigated worldwide increase exponentially with the accumulation of five key features: no take, well enforced, old (>10 years), large (>100 km(2)), and isolated by deep water or sand. Using effective MPAs with four or five key features as an unfished standard, comparisons of underwater survey data from effective MPAs with predictions based on survey data from fished coasts indicate that total fish biomass has declined about two-thirds from historical baselines as a result of fishing. Effective MPAs also had twice as many large (>250 mm total length) fish species per transect, five times more large fish biomass, and fourteen times more shark biomass than fished areas. Most (59%) of the MPAs studied had only one or two key features and were not ecologically distinguishable from fished sites. Our results show that global conservation targets based on area alone will not optimize protection of marine biodiversity. More emphasis is needed on better MPA design, durable management and compliance to ensure that MPAs achieve their desired conservation value.