Impacts of Birds vs. Invertebrate Predators on Rocky Intertidal Community Structure

Most studies of species interactions in rocky intertidal communities focus on invertebrate predators and herbivores interacting with sessile invertebrates and macrophytes. However, shorebirds are usually a conspicuous presence on rocky shores and eat sessile and mobile invertebrate prey, often including invertebrate predators and herbivores. Inspired by classic studies of strong bird predation effects in rocky intertidal habitats in Washington state (USA) and South Africa, I tested the effects of bird and invertebrate (sea stars, whelks) predation at multiple sites, wave exposures, and zones on the central Oregon coast from spring 1996 to fall 1997. To gain insight into the effects of birds relative to the effects of invertebrate predators, I used a crossed design, with bird exclusions (present and absent) and invertebrate predator removal (present and reduced). Compared to Washington state and South Africa, birds had little effect on the abundance of sessile or mobile prey in wave-exposed mid, wave-exposed low, and wave-protected mid zones at 2-4 sites. I suggest that differences between Oregon results and those in Washington and South Africa were driven by differences in bird abundance associated with whether the study site had resident colonies of shorebirds (primarily gulls, crows, and oystercatchers). That is, offshore islands often have resident breeding colonies such as in the Washington and South African studies, while sites in this study were all on the mainland where gulls were mostly transient visitors, while resident oystercatchers were usually limited to one or two pairs per site. Comparison with other marine and terrestrial experimental tests suggests that top-down effects of birds often vary in strength, and thus, future investigations should seek to understand the factors that underlie this variation.

Upwelling effects across different levels of biological organization: Integrating biochemical constituents, physiological performance and muscle transcriptomic response in the intertidal FISH Girella laevifrons (kyphosidae)

The physical-chemical variability of coastal upwelling creates a mosaic of environmental conditions that affect different levels of biological organization. Understanding the mechanisms that organisms use to cope with this variability is critical for addressing the challenges that climate change imposes on coastal ecosystems. This study integrates information on transcriptomic traits, metabolic performance, and the quantity of organic biomolecules in the intertidal fish Girella laevifrons from four locations with varying upwelling intensities. The results show that fish from locations with stronger upwelling intensity have higher levels of glucose, lipids, and proteins in their muscle tissue, in addition to better physiological performance compared to fish from sites with weaker upwelling intensity. Transcriptomic analyses revealed that genes associated with multicellular development and oxygen metabolism are more highly expressed in sites with stronger upwelling intensity, whereas genes related to protein ubiquitination are more expressed in sites with weaker upwelling intensity. In response to the mosaic of upwelling intensities (SAM-SST), and in-situ temperature, nutrients and oxygen variation observed in field, fish showed differential responses, suggesting local adaptations process that maximize ecological success in these areas with different physical-chemical conditions. Future studies should consider the integration of molecular tools to better understand the responses of organisms to environmental variability as upwelling intensities. This will help elucidate the complex interactions between environmental factors and biological responses, providing insights into how marine organisms might adapt to changing conditions. Understanding these mechanisms is essential for predicting the impacts of climate change on coastal ecosystems and for developing effective conservation and management strategies. The integration of transcriptomic data with metabolic and physiological performance measures offers a comprehensive approach to studying the adaptive responses of marine organisms to their dynamic environments considering the future responses in face to predict global change.

Exploring the relationship between upwelling intensity and socio-ecological attributes of marine exploitation areas for benthic resources (MEABRs), along the southern Humboldt Current system

The Eastern Boundary Upwelling Systems (EBUS) sustains some of the most productive marine systems on Earth. Within each of these systems, the upwelling process exhibits spatial and temporal variation resulting in marked differences in upwelling intensity and seasonality along extensive coastlines. The study of this variation is well needed, given the magnitude of the services provided by upwelling, and the impending impacts of global warming on EBUS. The critical link between the physical variability associated with upwelling intensity and its consequences on socio-ecological variables remain severely understudied. This study aimed to address such a gap by exploring the influence of coastal upwelling intensity on socio-ecological attributes of co-management units named Marine Exploitation Areas for Benthic Resources (MEABRs), along one the most productive ecosystems of the world: The southern Humboldt Current system. We evaluated the non-linear influence of upwelling on 1) the harvest of economically important resources, 2) the number of exploited species, 3) number and gender distribution of fishers involved, and 4) fishery activities. Our data indicated that on the one hand the annual harvest of commercial species, and all the exploited resources combined, were consistently higher in MEABRs associated with intermediate to high upwelling intensities. On the other hand, the harvest of kelp, the number of species harvested, and the number of fishers per MEABR increased towards low upwelling intensities, showing signals of fishery diversification. Interestingly, representation of female fishers increased towards high upwelling intensities, suggesting that multiple factors account for the variation in this, and other socio-ecological variables examined. Our study provides first-hand information about harvest levels and the allocation of fishery activities and gender distribution when MEABRs associated with different upwelling intensities are compared. Such information will assist in the identification of ecological and social vulnerabilities in a global warming scenario.

The importance of upwelling conditions as drivers of feeding behavior and thermal tolerance in a prominent intertidal fish

Upwelling, as a large oceanographic phenomenon, increases coastal productivity and influences all levels of biological complexity. Despite decades of research on it, much remains to be understood about the impact of upwelling on the feeding behavior and thermal tolerance of important groups such as fish. Hence, our aim was to investigate how upwelling conditions modify the feeding behavior and thermal tolerance of a prominent intertidal fish, Girella laevifrons. We collected purple mussels (Perumytilus purpuratus) from upwelling (U) and downwelling sites (DU) in central Chile, and used them as prey in feeding trials and measuring the concentration of organic matter and proteins in their tissues. We assessed fish consumption rates and growth in fish collected from the same U and DU sites, feeding on either U or DU mussels. Lastly, we assessed the thermal tolerance of U and DU fish fed with the aforementioned U vs DU mussels. We found that U mussels held higher concentrations of organic matter and proteins compared to their DU counterparts. U mussels were also selected and consumed in larger amounts than DU mussels, although the origin of the fish also influenced consumption rates. Thermal tolerance assays revealed that U fish exhibited higher maximum performance (Max.pf) and critical thermal maxima (Ctmax) and lower sensitivity to temperature changes (as measured by Q10), compared to DU fish. Altogether, these results point to a strong influence of upwelling on the quality of organisms' tissues, indirectly altering key aspects of fish feeding behavior and thermal tolerance. These findings also contribute to understanding the physiological adjustments organisms make in productive upwelling systems, and how they may change in the future with ongoing climate events.

A global synthesis of predation on bivalves

Predation is a dominant structuring force in ecological communities. In aquatic environments, predation on bivalves has long been an important focal interaction for ecological study because bivalves have central roles as ecosystem engineers, basal components of food webs, and commercial commodities. Studies of bivalves are common, not only because of bivalves' central roles, but also due to the relative ease of studying predatory effects on this taxonomic group. To understand patterns in the interactions of bivalves and their predators we synthesised data from 52 years of peer-reviewed studies on bivalve predation. Using a systematic search, we compiled 1334 studies from 75 countries, comprising 61 bivalve families (N = 2259), dominated by Mytilidae (29% of bivalves), Veneridae (14%), Ostreidae (8%), Unionidae (7%), and Dreissenidae and Tellinidae (6% each). A total of 2036 predators were studied, with crustaceans the most studied predator group (34% of predators), followed by fishes (24%), molluscs (17%), echinoderms (10%) and birds (6%). The majority of studies (86%) were conducted in marine systems, in part driven by the high commercial value of marine bivalves. Studies in freshwater ecosystems were dominated by non-native bivalves and non-native predator species, which probably reflects the important role of biological invasions affecting freshwater biodiversity. In fact, while 81% of the studied marine bivalve species were native, only 50% of the freshwater species were native to the system. In terms of approach, most studies used predation trials, visual analysis of digested contents and exclusion experiments to assess the effects of predation. These studies reflect that many factors influence bivalve predation depending on the species studied, including (i) species traits (e.g. behaviour, morphology, defence mechanisms), (ii) other biotic interactions (e.g. presence of competitors, parasites or diseases), and (iii) environmental context (e.g. temperature, current velocity, beach exposure, habitat complexity). There is a lack of research on the effects of bivalve predation at the population and community and ecosystem levels (only 7% and 0.5% of studies respectively examined impacts at these levels). At the population level, the available studies demonstrate that predation can decrease bivalve density through consumption or the reduction of recruitment. At the community and ecosystem level, predation can trigger effects that cascade through trophic levels or effects that alter the ecological functions bivalves perform. Given the conservation and commercial importance of many bivalve species, studies of predation should be pursued in the context of global change, particularly climate change, acidification and biological invasions.

Resistance of rocky intertidal communities to oceanic climate fluctuations

A powerful way to predict how ecological communities will respond to future climate change is to test how they have responded to the climate of the past. We used climate oscillations including the Pacific Decadal Oscillation (PDO), North Pacific Gyre Oscillation, and El Niño Southern Oscillation (ENSO) and variation in upwelling, air temperature, and sea temperatures to test the sensitivity of nearshore rocky intertidal communities to climate variability. Prior research shows that multiple ecological processes of key taxa (growth, recruitment, and physiology) were sensitive to environmental variation during this time frame. We also investigated the effect of the concurrent sea star wasting disease outbreak in 2013-2014. We surveyed nearly 150 taxa from 11 rocky intertidal sites in Oregon and northern California annually for up to 14-years (2006-2020) to test if community structure (i.e., the abundance of functional groups) and diversity were sensitive to past environmental variation. We found little to no evidence that these communities were sensitive to annual variation in any of the environmental measures, and that each metric was associated with < 8.6% of yearly variation in community structure. Only the years elapsed since the outbreak of sea star wasting disease had a substantial effect on community structure, but in the mid-zone only where spatially dominant mussels are a main prey of the keystone predator sea star, Pisaster ochraceus. We conclude that the established sensitivity of multiple ecological processes to annual fluctuations in climate has not yet scaled up to influence community structure. Hence, the rocky intertidal system along this coastline appears resistant to the range of oceanic climate fluctuations that occurred during the study. However, given ongoing intensification of climate change and increasing frequencies of extreme events, future responses to climate change seem likely.

The influence of upwelling on key bivalves from the Humboldt and Iberian current systems

Eastern Boundary Upwelling Systems (EBUS) deliver cold, nutrient-rich waters, influencing coastal biota from the molecular to the ecosystem level. Although local upwelling (U) and downwelling (DU) conditions are often known, their influence on body attributes of relevant species has not been systematically compared within and between EBUS (i.e., below and above regional scales). Hence, we compared the physical-chemical characteristics of U and DU sites in the Humboldt Current system (Chile) and the Iberian Current system (Portugal). We then assessed the influence of U and DU upon eight body attributes in purple mussels (Perumytilus purpuratus) and Mediterranean mussels (Mytilus galloprovincialis), from the Humboldt and Iberian systems, respectively. We hypothesized that bivalves from U sites display better fitness, as measured by body attributes, regardless of their origin (EBUS). As expected, waters from U sites in both systems showed lower temperatures and pH, and higher nitrite concentrations. We also found that mussels from U sites showed better fitness than those in DU sites in 12 out of 16 direct U vs DU comparisons. Shell length, shell volume, organic content of soft-tissues, and mechanical properties of the shell averaged consistently higher in mussels from U sites in both Current systems. In addition, total weight, soft-tissue weight, shell weight and shell thickness were all higher in the U site at the Humboldt system but had less consistent differences at the Iberian system. Altogether, most results supported our working hypothesis and indicate that U conditions support better fitted mussels. The few attributes that did not exhibit the expected U vs DU differences in the Iberian system suggest that local and species-specific differences also play a role on the attributes of these species. These results may also serve as a reference point for further studies addressing the influence of upwelling in these productive, critically important systems.

Upwelling enhances seaweed nutrient quality, altering feeding behavior and growth rates in an intertidal sea urchin, Loxechinus albus

Upwelling systems deliver nutrient-rich water into coastal ecosystems, influencing primary productivity and potentially altering seaweed-herbivore interactions. Upwelling bottom-up effects on distinct trophic levels are well-known. However, their influence on seaweed biomolecules and on algae-herbivore interactions and growth are less known. The aim of this study was threefold: i) to compare physical-chemical characteristics and nutrient levels in the water of upwelling (U) and downwelling (DU) zones, ii) to quantify their influence on the content of protein and carbohydrates in seaweed tissues of representative U and DU locations, and iii) to experimentally assess their effect on the feeding behavior and growth of a prominent intertidal herbivore, the sea urchin Loxechinus albus. Waters from U zones showed lower temperatures and pH, and higher phosphate concentrations than those from downwelling zones. Similarly, the tissue of seaweeds from a U location had significantly more proteins and carbohydrates than those from a DU location. The origin location of the sea urchins had a significant influence on consumption and growth rates: in general, those coming from a site with U conditions consumed and grew more than those coming from DU conditions. The quality of the algae was a significant factor on consumption rates, although in the case of preference trials, this factor interacted with sea urchin origin location. Our results show that the availability and quality of the food in upwelling zones has an influence on herbivore-seaweed direct interactions. However, these interactions and the growth of the sea urchins were also related to the coastal site and conditions from which the sea urchins came from. These results are relevant considering the expected impact of climate change on the world's oceans, and the importance of U zones as thermal (cold water) refuges for marine ectotherms.

Geographical variation of multiplex ecological networks in marine intertidal communities

Understanding the drivers of geographical variation in species distributions, and the resulting community structure, constitutes one of the grandest challenges in ecology. Geographical patterns of species richness and composition have been relatively well studied. Less is known about how the entire set of trophic and non-trophic ecological interactions, and the complex networks that they create by gluing species together in complex communities, change across geographical extents. Here, we compiled data of species composition and three types of ecological interactions occurring between species in rocky intertidal communities across a large spatial extent (~970 km of shoreline) of central Chile, and analyzed the geographical variability in these multiplex networks (i.e., comprising several interaction types) of ecological interactions. We calculated nine network summary statistics common across interaction types, and additional network attributes specific to each of the different types of interactions. We then investigated potential environmental drivers of this multivariate network organization. These included variation in sea surface temperature and coastal upwelling, the main drivers of productivity in nearshore waters. Our results suggest that structural properties of multiplex ecological networks are affected by local species richness and modulated by factors influencing productivity and environmental predictability. Our results show that non-trophic negative interactions are more sensitive to spatially structured temporal environmental variation than feeding relationships, with non-trophic positive interactions being the least labile to it. We also show that environmental effects are partly mediated through changes in species richness and partly through direct influences on species interactions, probably associated to changes in environmental predictability and to bottom-up nutrient availability. Our findings highlight the need for a comprehensive picture of ecological interactions and their geographical variability if we are to predict potential effects of environmental changes on ecological communities.

Biogeography of ocean acidification: Differential field performance of transplanted mussels to upwelling-driven variation in carbonate chemistry

Ocean acidification (OA) represents a serious challenge to marine ecosystems. Laboratory studies addressing OA indicate broadly negative effects for marine organisms, particularly those relying on calcification processes. Growing evidence also suggests OA combined with other environmental stressors may be even more deleterious. Scaling these laboratory studies to ecological performance in the field, where environmental heterogeneity may mediate responses, is a critical next step toward understanding OA impacts on natural communities. We leveraged an upwelling-driven pH mosaic along the California Current System to deconstruct the relative influences of pH, ocean temperature, and food availability on seasonal growth, condition and shell thickness of the ecologically dominant intertidal mussel Mytilus californianus. In 2011 and 2012, ecological performance of adult mussels from local and commonly sourced populations was measured at 8 rocky intertidal sites between central Oregon and southern California. Sites coincided with a large-scale network of intertidal pH sensors, allowing comparisons among pH and other environmental stressors. Adult California mussel growth and size varied latitudinally among sites and inter-annually, and mean shell thickness index and shell weight growth were reduced with low pH. Surprisingly, shell length growth and the ratio of tissue to shell weight were enhanced, not diminished as expected, by low pH. In contrast, and as expected, shell weight growth and shell thickness were both diminished by low pH, consistent with the idea that OA exposure can compromise shell-dependent defenses against predators or wave forces. We also found that adult mussel shell weight growth and relative tissue mass were negatively associated with increased pH variability. Including local pH conditions with previously documented influences of ocean temperature, food availability, aerial exposure, and origin site enhanced the explanatory power of models describing observed performance differences. Responses of local mussel populations differed from those of a common source population suggesting mussel performance partially depended on genetic or persistent phenotypic differences. In light of prior research showing deleterious effects of low pH on larval mussels, our results suggest a life history transition leading to greater resilience in at least some performance metrics to ocean acidification by adult California mussels. Our data also demonstrate “hot” (more extreme) and “cold” (less extreme) spots in both mussel responses and environmental conditions, a pattern that may enable mitigation approaches in response to future changes in climate.

Complementarity in spatial subsidies of carbon associated with resource partitioning along multiple niche axes

The co-occurrence of several species, all of which share similar resource requirements, remains a paradox in ecology. Here, I evaluated resource use along multiple environmental axes to understand the potential for niche partitioning and complementarity in a guild of suspension-feeding rocky shore invertebrates. I focused on the mussels Perna canaliculus, Mytilus galloprovincialis, Aulacomya maoriana, and Xenostrobus pulex, all of which coexist along the coastline of the South Island of New Zealand. I quantified the mussels' distributions, both vertically on the shore and within the three-dimensional mussel bed matrix, and used carbon (δ¹³C) and nitrogen (δ¹⁵N) stable isotope ratios to compare species' diets. Mussels exhibited niche partitioning along all resource axes. Given the mussels' use of different spatial and food resources, I evaluated the potential for complementarity with respect to mussels' roles as mediators of an important spatial subsidy, carbon inputs from the nearshore ocean into rocky-shore ecosystems. In these systems, mussels are basal species, capturing and consuming particulate organic matter in the ocean and making it available for local consumption within the benthic community. I found that mussel diversity matters; even the most productive species-Perna canaliculus-only contributed about half of the mussel-mediated carbon that accumulated over a year. Multiple co-occurring species are, therefore, likely to provide more carbon at the base of the intertidal food web than any one species living independently, and complementarity in resource use along multiple environmental axes is likely an important mechanism underlying this relationship.

Historical comparisons of body size are sensitive to data availability and ecological context

Historical comparisons of body size often lack pertinent details, including information on the sampling protocol and relevant ecological covariates that influence body size. Moreover, historical estimates of body size that rely on museum specimens may be biased towards larger size classes because of collector preferences, and thus size thresholds have been used to focus attention on maximum body size. We tested the consequences of sampling design, ecological covariates, and size thresholds on inferences of body-size change using field-contextualized historical records, rather than museum specimens. In 2014-2015, we revisited historical (1947-1963) size-frequency distributions of three gastropods (Tegula funebralis, Lottia digitalis/L. austrodigitalis, Littorina keenae) in the context of population density and tidal height. In general, gastropods declined in size. However, our inferences regarding body-size decline were tempered when the variation between sampling units was taken into consideration, resulting in greater uncertainty around the estimate of proportional change in body size. Gastropod size was correlated with population density and tidal height, and these relationships varied over time. Finally, the magnitude and direction of body-size change varied with the amount of data available for analysis, demonstrating that the use of size thresholds can lead to incomplete conclusions.

Climate shapes population variation in dogwhelk predation on foundational mussels

Trait variation among populations is important for shaping ecological dynamics. In marine intertidal systems, seawater temperature, low tide emersion temperature, and pH can drive variation in traits and affect species interactions. In western North America, Nucella dogwhelks are intertidal drilling predators of the habitat-forming mussel Mytilus californianus. Nucella exhibit local adaptation, but it is not known to what extent environmental factors and genetic structure contribute to variation in prey selectivity among populations. We surveyed drilled mussels at sites across Oregon and California, USA, and used multiple regression and Mantel tests to test the effects of abiotic factors and Nucella neutral genetic relatedness on the size of mussels drilled across sites. Our results show that Nucella at sites characterized by higher and less variable temperature and pH drilled larger mussels. Warmer temperatures appear to induce faster handling time, and more stable pH conditions may prolong opportunities for active foraging by reducing exposure to repeated stressful conditions. In contrast, there was no significant effect of genetic relatedness on prey size selectivity. Our results emphasize the role of climate in shaping marine predator selectivity on a foundation species. As coastal climates change, predator traits will respond to localized environmental conditions, changing ecological interactions.

Species co-occurrence networks: Can they reveal trophic and non-trophic interactions in ecological communities?

Co-occurrence methods are increasingly utilized in ecology to infer networks of species interactions where detailed knowledge based on empirical studies is difficult to obtain. Their use is particularly common, but not restricted to, microbial networks constructed from metagenomic analyses. In this study, we test the efficacy of this procedure by comparing an inferred network constructed using spatially intensive co-occurrence data from the rocky intertidal zone in central Chile to a well-resolved, empirically based, species interaction network from the same region. We evaluated the overlap in the information provided by each network and the extent to which there is a bias for co-occurrence data to better detect known trophic or non-trophic, positive or negative interactions. We found a poor correspondence between the co-occurrence network and the known species interactions with overall sensitivity (probability of true link detection) equal to 0.469, and specificity (true non-interaction) equal to 0.527. The ability to detect interactions varied with interaction type. Positive non-trophic interactions such as commensalism and facilitation were detected at the highest rates. These results demonstrate that co-occurrence networks do not represent classical ecological networks in which interactions are defined by direct observations or experimental manipulations. Co-occurrence networks provide information about the joint spatial effects of environmental conditions, recruitment, and, to some extent, biotic interactions, and among the latter, they tend to better detect niche-expanding positive non-trophic interactions. Detection of links (sensitivity or specificity) was not higher for well-known intertidal keystone species than for the rest of consumers in the community. Thus, as observed in previous empirical and theoretical studies, patterns of interactions in co-occurrence networks must be interpreted with caution, especially when extending interaction-based ecological theory to interpret network variability and stability. Co-occurrence networks may be particularly valuable for analysis of community dynamics that blends interactions and environment, rather than pairwise interactions alone.

Upwelling-derived oceanographic conditions impact growth performance and growth-related gene expression in intertidal fish

Growth is one of the main biological processes in aquatic organisms that is affected by environmental fluctuations such as upwelling (characterized by food-rich waters). In fish, growth is directly related with skeletal muscle increase; which represents the largest tissue of body mass. However, the effects of upwelling on growth, at the physiological and molecular level, are unknown. This study used Girella laevifrons (one of the most abundant intertidal fish in Eastern South Pacific) as a biological model, considering animals from upwelling (U) and non-upwelling (NU) areas. Here, we evaluated the effect of nutritional composition and food availability on growth performance and expression of key growth-related genes (insulin-kike growth factor 1 (igf1) and myosin heavy-chain (myhc)) and atrophy-related genes (muscle ring-finger 1 (murf1), F-box only protein 32 (atrogin-1) and BCL2/adenovirus E1B 19kDa-interacting protein 3 (bnip3)). We reported that, among zones, U fish displayed higher growth performance in response to nutritional composition, specifically between protein- and fiber-rich diets (~1g). We also found in NU fish that atrophy-related genes were upregulated with fiber-rich diet and during fasting (~2-fold at minimum respect U). In conclusion, our results suggest that the growth potential of upwelling fish may be a consequence of differential muscle gene expression. Our data provide a preliminary approach contributing on how upwelling influence fish growth at the physiological and molecular levels. Future studies are required to gain further knowledge about molecular differences between U and NU animals, as well as the possible applications of this knowledge in the aquaculture industry.

Stoichiometric Mismatch between Consumers and Resources Mediates the Growth of Rocky Intertidal Suspension Feeders

The concept of ecological stoichiometry-the balancing of elemental ratios in ecological interactions-has transformed our thinking about processes in natural systems. Here, this perspective is applied to rocky shore ecosystems to explore the consequences of variation in internal nutrient ratios across two trophic levels. Specifically, I measured the internal concentrations of carbon (C) and nitrogen (N) in mussels (Mytilus spp.) and particulate organic matter (POM) to evaluate the effects of stoichiometric mismatch-the difference in the carbon-to-nitrogen ratio (C:N) between a consumer and its resources-on mussel growth at sites on the coasts of Oregon, USA, and the South Island of New Zealand. As POM quality (i.e., Chl a, a proxy for phytoplankton availability in the POM) increased, C:N of the POM declined, but C:N of mussels increased. This resulted in a greater mismatch in C:N between mussels and their food source at low Chl a. Mussel growth across sites was positively associated with Chl a, particulate organic carbon (POC), and particulate organic nitrogen (PON) but negatively associated with stoichiometric mismatch. Overall, as the elemental ratios of consumers became more different from those of their resources, growth declined, likely due to the energetic cost associated with processing lower quality food. Furthermore, the effect of food quantity on growth depended on stoichiometric mismatch. In New Zealand, where mismatch was high-i.e., consumer C:N differed substantially from resource C:N-consumer growth was strongly affected by resource quantity (Chl a or POC). However, in Oregon, where mismatch was low, the relationship between resource quantity and growth was considerably weaker. This interaction between resource quantity and mismatch was not apparent for PON, which is consistent with variation in PON underlying variation in POM C:N and highlights the role of N in limiting growth. Previous research has neglected the importance of ecological stoichiometry as a mediator of consumer-resource interactions in rocky intertidal communities. I show that resource quality and quantity interact to determine consumer growth, highlighting the utility of ecological stoichiometry in understanding spatial subsidies in benthic marine systems.

Incorporating Context Dependency of Species Interactions in Species Distribution Models

SYNOPSIS

Species distribution models typically use correlative approaches that characterize the species-environment relationship using occurrence or abundance data for a single species. However, species distributions are determined by both abiotic conditions and biotic interactions with other species in the community. Therefore, climate change is expected to impact species through direct effects on their physiology and indirect effects propagated through their resources, predators, competitors, or mutualists. Furthermore, the sign and strength of species interactions can change according to abiotic conditions, resulting in context-dependent species interactions that may change across space or with climate change. Here, we incorporated the context dependency of species interactions into a dynamic species distribution model. We developed a multi-species model that uses a time-series of observational survey data to evaluate how abiotic conditions and species interactions affect the dynamics of three rocky intertidal species. The model further distinguishes between the direct effects of abiotic conditions on abundance and the indirect effects propagated through interactions with other species. We apply the model to keystone predation by the sea star Pisaster ochraceus on the mussel Mytilus californianus and the barnacle Balanus glandula in the rocky intertidal zone of the Pacific coast, USA. Our method indicated that biotic interactions between P. ochraceus and B. glandula affected B. glandula dynamics across >1000 km of coastline. Consistent with patterns from keystone predation, the growth rate of B. glandula varied according to the abundance of P. ochraceus in the previous year. The data and the model did not indicate that the strength of keystone predation by P. ochraceus varied with a mean annual upwelling index. Balanus glandula cover increased following years with high phytoplankton abundance measured as mean annual chlorophyll-a. M. californianus exhibited the same pattern to a lesser degree, although this pattern was not significant. This work bridges the disciplines of biogeography and community ecology to develop tools to better understand the direct and indirect effects of abiotic conditions on ecological communities.

Linking macroecology and community ecology: refining predictions of species distributions using biotic interaction networks

Macroecological models for predicting species distributions usually only include abiotic environmental conditions as explanatory variables, despite knowledge from community ecology that all species are linked to other species through biotic interactions. This disconnect is largely due to the different spatial scales considered by the two sub-disciplines: macroecologists study patterns at large extents and coarse resolutions, while community ecologists focus on small extents and fine resolutions. A general framework for including biotic interactions in macroecological models would help bridge this divide, as it would allow for rigorous testing of the role that biotic interactions play in determining species ranges. Here, we present an approach that combines species distribution models with Bayesian networks, which enables the direct and indirect effects of biotic interactions to be modelled as propagating conditional dependencies among species' presences. We show that including biotic interactions in distribution models for species from a California grassland community results in better range predictions across the western USA. This new approach will be important for improving estimates of species distributions and their dynamics under environmental change.

Keystone predation and molecules of keystone significance

Keystone species structure ecological communities and are major determinants of biodiversity. A synthesis of research on keystone species is nonetheless missing a critical component - the sensory mechanisms for behavioral interactions that determine population- and community-wide attributes. Here, we establish the chemosensory basis for keystone predation by sea stars (Pisaster ochraceus) on mussels. This consumer-resource interaction is prototypic of top-down driven trophic cascades. Each mussel species (Mytilus californianus and M. galloprovincialis) secretes a glycoprotein orthologue (29.6 and 28.1 kDa, respectively) that acts, singularly, to evoke the sea star predatory response. The orthologues (named "KEYSTONEin") are localized in the epidermis, extrapallial fluid, and organic shell coating (periostracum) of live, intact mussels. Thus, KEYSTONEin contacts chemosensory receptors on tube feet as sea stars crawl over rocky surfaces in search of prey. The complete nucleotide sequences reveal that KEYSTONEin shares 87% (M. californianus) or 98% (M. galloprovincialis) homology with a calcium-binding protein in the shell matrix of a closely related congener, M. edulis. All three molecules cluster tightly within the Complement Component 1 Domain Containing (C1qDC) protein family; each exhibits a large globular domain, low complexity region(s), coiled coil, and at least four of five histidine-aspartic acid tandem motifs. Collective results support the hypothesis that KEYSTONEin evolved ancestrally in immunological, and later, in biomineralization roles. More recently, the substance has become exploited by sea stars as a contact cue for prey recognition. As the first identified compound to evoke keystone predation, KEYSTONEin provides valuable sensory information, promotes biodiversity, and shapes community structure and function. Without this molecule, there would be no predation by sea stars on mussels.

Surfzone hydrodynamics as a key determinant of spatial variation in rocky intertidal communities

Larvae of intertidal species develop at sea and must return to adult habitats to replenish populations. Similarly, nutrients, detritus and plankton provide important subsidies spurring growth and reproduction of macroalgae and filter-feeding invertebrates that form the foundation of intertidal communities. Together, these factors determine the density and intensity of interactions among community members. We hypothesized that spatial variation in surfzone hydrodynamics affects the delivery of plankton subsidies. We compared entire zooplankton communities inside and outside the surf zone daily while monitoring physical conditions for one month each at two shores with different surfzone characteristics. Opposite cross-shore distributions of larvae and other zooplankters occurred at the two sites: zooplankton was much more abundant inside the mildly sloping dissipative surf zone (DSZ) with rip currents and was more abundant outside the steep reflective surf zone (RSZ). Biophysical numerical simulations demonstrated that zooplankters were concentrated in rip channels of the DSZ and were mostly unable to enter the RSZ, indicating the hydrodynamic processes behind the observed spatial variation of zooplankters in the surf zone. Differences in the concentration of larvae and other zooplankters between the inner shelf and surf zone may be an underappreciated, key determinant of spatial variation in inshore communities.

Interacting environmental mosaics drive geographic variation in mussel performance and predation vulnerability

Although theory suggests geographic variation in species' performance is determined by multiple niche parameters, little consideration has been given to the spatial structure of interacting stressors that may shape local and regional vulnerability to global change. Here, we use spatially explicit mosaics of carbonate chemistry, food availability and temperature spanning 1280 km of coastline to test whether persistent, overlapping environmental mosaics mediate the growth and predation vulnerability of a critical foundation species, the mussel Mytilus californianus. We find growth was highest and predation vulnerability was lowest in dynamic environments with frequent exposure to low pH seawater and consistent food. In contrast, growth was lowest and predation vulnerability highest when exposure to low pH seawater was decoupled from high food availability, or in exceptionally warm locations. These results illustrate how interactions among multiple drivers can cause unexpected, yet persistent geographic mosaics of species performance, interactions and vulnerability to environmental change.

Sea Star Wasting Disease in the Keystone Predator Pisaster ochraceus in Oregon: Insights into Differential Population Impacts, Recovery, Predation Rate, and Temperature Effects from Long-Term Research

Sea star wasting disease (SSWD) first appeared in Oregon in April 2014, and by June had spread to most of the coast. Although delayed compared to areas to the north and south, SSWD was initially most intense in north and central Oregon and spread southward. Up to 90% of individuals showed signs of disease from June-August 2014. In rocky intertidal habitats, populations of the dominant sea star Pisaster ochraceus were rapidly depleted, with magnitudes of decline in density among sites ranging from -2x to -9x (59 to 84%) and of biomass from -2.6x to -15.8x (60 to 90%) by September 2014. The frequency of symptomatic individuals declined over winter and persisted at a low rate through the spring and summer 2015 (~5-15%, at most sites) and into fall 2015. Disease expression included six symptoms: initially with twisting arms, then deflation and/or lesions, lost arms, losing grip on substrate, and final disintegration. SSWD was disproportionally higher in orange individuals, and higher in tidepools. Although historically P. ochraceus recruitment has been low, from fall 2014 to spring 2015 an unprecedented surge of sea star recruitment occurred at all sites, ranging from ~7x to 300x greater than in 2014. The loss of adult and juvenile individuals in 2014 led to a dramatic decline in predation rate on mussels compared to the previous two decades. A proximate cause of wasting was likely the "Sea Star associated Densovirus" (SSaDV), but the ultimate factors triggering the epidemic, if any, remain unclear. Although warm temperature has been proposed as a possible trigger, SSWD in Oregon populations increased with cool temperatures. Since P. ochraceus is a keystone predator that can strongly influence the biodiversity and community structure of the intertidal community, major community-level responses to the disease are expected. However, predicting the specific impacts and time course of change across west coast meta-communities is difficult, suggesting the need for detailed coast-wide investigation of the effects of this outbreak.

Patterns and Variation in Benthic Biodiversity in a Large Marine Ecosystem

While there is a persistent inverse relationship between latitude and species diversity across many taxa and ecosystems, deviations from this norm offer an opportunity to understand the conditions that contribute to large-scale diversity patterns. Marine systems, in particular, provide such an opportunity, as marine diversity does not always follow a strict latitudinal gradient, perhaps because several hypothesized drivers of the latitudinal diversity gradient are uncorrelated in marine systems. We used a large scale public monitoring dataset collected over an eight year period to examine benthic marine faunal biodiversity patterns for the continental shelf (55-183 m depth) and slope habitats (184-1280 m depth) off the US West Coast (47°20'N-32°40'N). We specifically asked whether marine biodiversity followed a strict latitudinal gradient, and if these latitudinal patterns varied across depth, in different benthic substrates, and over ecological time scales. Further, we subdivided our study area into three smaller regions to test whether coast-wide patterns of biodiversity held at regional scales, where local oceanographic processes tend to influence community structure and function. Overall, we found complex patterns of biodiversity on both the coast-wide and regional scales that differed by taxonomic group. Importantly, marine biodiversity was not always highest at low latitudes. We found that latitude, depth, substrate, and year were all important descriptors of fish and invertebrate diversity. Invertebrate richness and taxonomic diversity were highest at high latitudes and in deeper waters. Fish richness also increased with latitude, but exhibited a hump-shaped relationship with depth, increasing with depth up to the continental shelf break, ~200 m depth, and then decreasing in deeper waters. We found relationships between fish taxonomic and functional diversity and latitude, depth, substrate, and time at the regional scale, but not at the coast-wide scale, suggesting that coast-wide patterns can obscure important correlates at smaller scales. Our study provides insight into complex diversity patterns of the deep water soft substrate benthic ecosystems off the US West Coast.

Behaviorally Mediated Larval Transport in Upwelling Systems

Highly advective upwelling systems along the western margins of continents are widely believed to transport larvae far offshore in surface currents resulting in larval wastage, limited recruitment, and increased population connectivity. However, suites of larval behaviors effectively mediate interspecific differences in the extent of cross-shelf migrations between nearshore adult habitats and offshore larval habitats. Interspecific differences in behavior determining whether larvae complete development in estuaries or migrate to the continental shelf are evident in large estuaries, but they sometimes may be disrupted by turbulent tidal flow or the absence of a low-salinity cue in shallow, low-flow estuaries, which are widespread in upwelling systems. Larvae of most species on the continental shelf complete development in the coastal boundary layer of reduced flow, whereas other species migrate to the mid- or outer shelf depending on how much time is spent in surface currents. These migrations are maintained across latitudinal differences in the strength and persistence of upwelling, in upwelling jets at headlands, over upwelling-relaxation cycles, and among years of varying upwelling intensity. Incorporating larval behaviors into numerical models demonstrates that larvae recruit closer to home and in higher numbers than when larvae disperse passively or remain in surface currents.

Survival and arm abscission are linked to regional heterothermy in an intertidal sea star

Body temperature is a more pertinent variable to physiological stress than ambient air temperature. Modeling and empirical studies on the impacts of climate change on ectotherms usually assume that body temperature within organisms is uniform. However, many ectotherms show significant within-body temperature heterogeneity. The relationship between regional heterothermy and the response of ectotherms to sublethal and lethal conditions remains underexplored. We quantified within-body thermal heterogeneity in an intertidal sea star (Pisaster ochraceus) during aerial exposure at low tide to examine the lethal and sublethal effects of temperatures of different body regions. In manipulative experiments, we measured the temperature of the arms and central disc, as well as survival and arm abscission under extreme aerial conditions. Survival was related strongly to central disc temperature. Arms were generally warmer than the central disc in individuals that survived aerial heating, but we found the reverse in those that died. When the central disc reached sublethal temperatures of 31-35°C, arms reached temperatures of 33-39°C, inducing arm abscission. The absolute temperature of individual arms was a poor predictor of arm abscission, but the arms lost were consistently the hottest at the within-individual scale. Therefore, the vital region of this sea star may remain below the lethal threshold under extreme conditions, possibly through water movement from the arms to the central disc and/or evaporative cooling, but at the cost of increased risk of arm abscission. Initiation of arm abscission seems to reflect a whole-organism response while death occurs as a result of stress acting directly on central disc tissues.

Can upwelling signals be detected in intertidal fishes of different trophic levels?

For intertidal fishes belonging to three species, the herbivore Scartichthys viridis (Blenniidae), the omnivore Girella laevifrons (Kyphosidae) and the carnivore Graus nigra (Kyphosidae), mass and body size relationships were higher in individuals from an upwelling zone compared with those from a non-upwelling zone. RNA:DNA were higher in the herbivores and omnivores from the upwelling zone. Higher biomass and RNA:DNA in the upwelling intertidal fishes may be a consequence of an increased exposure to higher nutrient availability, suggesting that increased physiological conditioning in vertebrates from upwelling areas can be detected and measured using intertidal fishes of different trophic levels.

Steady as he goes: at-sea movement of adult male Australian sea lions in a dynamic marine environment

The southern coastline of Australia forms part of the worlds' only northern boundary current system. The Bonney Upwelling occurs every austral summer along the south-eastern South Australian coastline, a region that hosts over 80% of the worlds population of an endangered endemic otariid, the Australian sea lion. We present the first data on the movement characteristics and foraging behaviour of adult male Australian sea lions across their South Australian range. Synthesizing telemetric, oceanographic and isotopic datasets collected from seven individuals enabled us to characterise individual foraging behaviour over an approximate two year time period. Data suggested seasonal variability in stable carbon and nitrogen isotopes that could not be otherwise explained by changes in animal movement patterns. Similarly, animals did not change their foraging patterns despite fine-scale spatial and temporal variability of the upwelling event. Individual males tended to return to the same colony at which they were tagged and utilized the same at-sea regions for foraging irrespective of oceanographic conditions or time of year. Our study contrasts current general assumptions that male otariid life history strategies should result in greater dispersal, with adult male Australian sea lions displaying central place foraging behaviour similar to males of other otariid species in the region.

Biotic interaction strength and the intensity of selection

Although the ecological and evolutionary impacts of species interactions have been the foci of much research, the relationship between the strength of species interactions and the intensity of selection has been investigated only rarely. I develop a simple model demonstrating how the opportunity for selection varies with interaction strength, and then use the relationship between the maximum value of the selection differential and the opportunity for selection (Arnold & Wade 1984) to evaluate how selection differentials vary in relation to species interaction strength. This model predicts an initial deceleration and then an accelerating increase in the intensity of selection with increasing strength of antagonistic interactions and with decreasing strength of mutualistic interactions. Empirical data from several studies provide support for this model. These results further support an evolutionary mechanism for some striking patterns of evolutionary diversification including the latitudinal species gradient, and should be relevant to studies of eco-evolutionary dynamics.

Community ecology and the evolution of molecules of keystone significance

Molecules of keystone significance are vital in structuring ecological communities. Select bioactive compounds can cause disproportionately large effects by connecting such seemingly disparate processes as microbial loop dynamics and apex predation. Here, we develop a general theory and propose mechanisms that could lead to the evolution of keystone molecules. Introduced into a respective community by one, or only a few, autotrophic or microbial species, these compounds often originate as chemical defenses. When co-opted by resistant consumer species, however, they are used either in chemical defense against higher-order predators or as chemosensory cues that elicit courtship and mating, alarm, and predatory search. Requisite to these multifunctional properties, biosynthetic capacity evolves along with mechanisms for resistance and/or toxin storage in primary producers. Subsequently, consumers acquire resistances or tolerances, and the toxins are transferred through food webs via trophic interactions. In consumers, mechanisms eventually evolve for recognizing toxins as feeding cues and, ultimately, as signals or pheromones in chemical communication within or between species. One, or a few, active compounds can thus mediate a vast array of physiological traits, expressed differentially across many species in a given community. Through convergent evolution, molecules of keystone significance provide critical information to phylogenetically diverse species, initiate major trophic cascades, and structure communities within terrestrial, freshwater, coastal-ocean and open-ocean habitats.

Gridded meteorological data as a resource for mechanistic macroecology in coastal environments

Gridded weather data were evaluated as sources of forcing variables for biophysical models of intertidal animal body temperature with model results obtained using local weather station data serving as the baseline of comparison. The objective of the study was to determine which gridded data are sufficient to capture observed patterns of thermal stress. Three coastal sites in western North America were included in this analysis: Boiler Bay, Oregon; Bodega Bay, California; and Pacific Grove, California. The gridded data with the highest spatial resolution, the 32-km North American Regional Reanalysis (NARR) and the 38-km Climate Forecasting System Reanalysis (CFSR), predicted daily maximum intertidal animal temperature most similarly to the local weather Station data. Time step size was important for variables that change rapidly throughout the day, such as solar radiation. There were site-based differences in the ability of the model to predict daily maximum intertidal animal temperature, with the gridded data predictions being the closest to local weather station predictions in Boiler Bay, Oregon. In a review of gridded data used as part of ecological studies, there was broad use of the data across subject areas and ecosystems so the recent improvements in the spatial (from 2 degrees to 32 km) and temporal scales (from 6 hours to 1 hour) of gridded data will further add to the applicability within the ecological community particularly for mechanistic studies.

Community regulation: the relative importance of recruitment and predation intensity of an intertidal community dominant in a seascape context

Predicting the strength and context-dependency of species interactions across multiple scales is a core area in ecology. This is especially challenging in the marine environment, where populations of most predators and prey are generally open, because of their pelagic larval phase, and recruitment of both is highly variable. In this study we use a comparative-experimental approach on small and large spatial scales to test the relationship between predation intensity and prey recruitment and their relative importance in shaping populations of a dominant rocky intertidal space occupier, mussels, in the context of seascape (availability of nearby subtidal reef habitat). Predation intensity on transplanted mussels was tested inside and outside cages and recruitment was measured with standard larval settlement collectors. We found that on intertidal rocky benches with contiguous subtidal reefs in New Zealand, mussel larval recruitment is usually low but predation on recruits by subtidal consumers (fish, crabs) is intense during high tide. On nearby intertidal rocky benches with adjacent sandy subtidal habitats, larval recruitment is usually greater but subtidal predators are typically rare and predation is weaker. Multiple regression analysis showed that predation intensity accounts for most of the variability in the abundance of adult mussels compared to recruitment. This seascape-dependent, predation-recruitment relationship could scale up to explain regional community variability. We argue that community ecology models should include seascape context-dependency and its effects on recruitment and species interactions for better predictions of coastal community dynamics and structure.