Increasing variance in North Pacific climate relates to unprecedented ecosystem variability off California

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.

Anchovy boom and bust linked to trophic shifts in larval diet

Although massive biomass fluctuations of coastal-pelagic fishes are an iconic example of the impacts of climate variability on marine ecosystems, the mechanisms governing these dynamics are often elusive. We construct a 45-year record of nitrogen stable isotopes measured in larvae of Northern Anchovy (Engraulis mordax) in the California Current Ecosystem to assess patterns in food chain length. Larval trophic efficiency associated with a shortened food chain increased larval survival and produced boom periods of high adult biomass. In contrast, when larval food chain length increased, and energy transfer efficiency decreased, the population crashed. We propose the Trophic Efficiency in Early Life (TEEL) hypothesis, which states that larval fishes must consume prey that confer sufficient energy for survival, to help explain natural boom-bust dynamics of coastal pelagic fishes. Our findings illustrate a potential for trophic indicators to generally inform larval survival and adult population dynamics of coastal-pelagic fishes.

Climate Change Impacts on Eastern Boundary Upwelling Systems

The world's eastern boundary upwelling systems (EBUSs) contribute disproportionately to global ocean productivity and provide critical ecosystem services to human society. The impact of climate change on EBUSs and the ecosystems they support is thus a subject of considerable interest. Here, we review hypotheses of climate-driven change in the physics, biogeochemistry, and ecology of EBUSs; describe observed changes over recent decades; and present projected changes over the twenty-first century. Similarities in historical and projected change among EBUSs include a trend toward upwelling intensification in poleward regions, mitigatedwarming in near-coastal regions where upwelling intensifies, and enhanced water-column stratification and a shoaling mixed layer. However, there remains significant uncertainty in how EBUSs will evolve with climate change, particularly in how the sometimes competing changes in upwelling intensity, source-water chemistry, and stratification will affect productivity and ecosystem structure. We summarize the commonalities and differences in historical and projected change in EBUSs and conclude with an assessment of key remaining uncertainties and questions. Future studies will need to address these questions to better understand, project, and adapt to climate-driven changes in EBUSs.

Diverse integrated ecosystem approach overcomes pandemic-related fisheries monitoring challenges

The COVID-19 pandemic caused unprecedented cancellations of fisheries and ecosystem-assessment surveys, resulting in a recession of observations needed for management and conservation globally. This unavoidable reduction of survey data poses challenges for informing biodiversity and ecosystem functioning, developing future stock assessments of harvested species, and providing strategic advice for ecosystem-based management. We present a diversified framework involving integration of monitoring data with empirical models and simulations to inform ecosystem status within the California Current Large Marine Ecosystem. We augment trawl observations collected from a limited fisheries survey with survey effort reduction simulations, use of seabird diets as indicators of fish abundance, and krill species distribution modeling trained on past observations. This diversified approach allows for evaluation of ecosystem status during data-poor situations, especially during the COVID-19 era. The challenges to ecosystem monitoring imposed by the pandemic may be overcome by preparing for unexpected effort reduction, linking disparate ecosystem indicators, and applying new species modeling techniques.

The Covid-19 pandemic has disrupted ecosystem and biodiversity monitoring programs, including marine fisheries surveys. Here the authors combine multiple modelling approaches and data to overcome lost observational effort off the coasts of California in a diversified integrated ecosystem approach.

Synchrony erodes spatial portfolios of an anadromous fish and alters availability for resource users

Environmental forces can create spatially synchronous dynamics among nearby populations. However, increased climate variability, driven by anthropogenic climate change, will likely enhance synchrony among spatially disparate populations. Population synchrony may lead to greater fluctuations in abundance, but the consequences of population synchrony across multiple scales of biological organization, including impacts to putative competitors, dependent predators or human communities, are rarely considered in this context. Chinook salmon Oncorhynchus tshawytscha stocks distribute across the Northeast Pacific, creating spatially variable portfolios that support large ocean fisheries and marine mammal predators, such as killer whales Orcinus orca. We rely on a multi-population model that simulates Chinook salmon ocean distribution and abundance to understand spatial portfolios, or variability in abundance within and among ocean distribution regions, of Chinook salmon stocks across 17 ocean regions from Southeast Alaska to California. We found the expected positive correlation between the number of stocks in an ocean region and spatial portfolio strength; however, increased demographic synchrony eroded Chinook salmon spatial portfolios in the ocean. Moreover, we observed decreased resource availability within ocean fishery management jurisdictions but not within killer whale summer habitat. We found a strong portfolio effect across both Southern Resident and Northern Resident killer whale habitats that was relatively unaffected by increased demographic synchrony, likely a result of the large spatial area included in these habitats. However, within the areas of smaller fishing management jurisdictions we found a weakening of Chinook salmon portfolios and increased but inconsistent likelihood of low abundance years as demographic synchrony increased. We suggest that management and conservation actions that reduce spatial synchrony can enhance short-term ecosystem resilience by promoting the stabilizing effect multiple stocks have on aggregate Chinook salmon populations and overall resource availability.

Trophic ecology of yellowtail rockfish (Sebastes flavidus) during a marine heat wave off central California, USA

The yellowtail rockfish, Sebastes flavidus, is a widespread and abundant mesopredator in the California Current Large Marine Ecosystem. We utilized stomach content and stable isotope analyses to investigate the trophic ecology of this species at three sites off central California just before (August–October 2013) and during (August and October 2014) a marine heat wave. Sebastes flavidus largely consumed pelagic prey (zooplankton and micronekton). Diets were dominated by tunicates (salps and pyrosomes), pelagic crustaceans (euphausiids, hyperid amphipods, larval decapods), and fishes, with the relative contribution of these prey taxa varying spatially (sample location, longitude, depth) and temporally (year, month), based on complementary multivariate analyses. Prey-specific indices demonstrated that individual S. flavidus diet composition typically was dominated by one of these prey groups, and that prey switching occurred based on the relative availability of prey and their energetic importance. Stable isotope analysis of δ¹⁵N indicated that the S. flavidus populations sampled in 2014 had been feeding at an elevated trophic position and more variable prey spectrum relative to 2013, probably as a consequence of greater piscivory and the incorporation of temporal changes in diet composition. Because its opportunistic feeding behavior reflects the dynamism and heterogeneity of the pelagic forage preyscape, S. flavidus may be an important ecosystem indicator species. For example, the novel incorporation of pyrosomes as a large portion of the diet of S. flavidus during 2013–2014 directly related to the massive increase in pyrosome abundance in the California Current during the 2014 marine heat wave.

Krill Hotspot Formation and Phenology in the California Current Ecosystem

In the California Current Ecosystem, krill represent a key link between primary production and higher trophic level species owing to their central position in the food web and tendency to form dense aggregations. However, the strongly advective circulation associated with coastal upwelling may decouple the timing, occurrence, and persistence of krill hotspots from phytoplankton biomass and nutrient sources. Results from a coupled physical-biological model provide insights into fundamental mechanisms controlling the phenology of krill hotspots in the California Current Ecosystem, and their sensitivity to alongshore changes in coastal upwelling intensity. The simulation indicates that dynamics controlling krill hotspot formation, intensity, and persistence on seasonal and interannual timescales are strongly heterogeneous and related to alongshore variations in upwelling-favorable winds, primary production, and ocean currents. Furthermore, regions promoting persistent krill hotspot formation coincide with increased observed abundance of top predators, indicating that the model resolves important ecosystem complexity and function.

Hypoxia Tolerance of 10 Euphausiid Species in Relation to Vertical Temperature and Oxygen Gradients

Oxygen Minimum Zones prevail in most of the world's oceans and are particularly extensive in Eastern Boundary Upwelling Ecosystems such as the Humboldt and the Benguela upwelling systems. In these regions, euphausiids are an important trophic link between primary producers and higher trophic levels. The species are known as pronounced diel vertical migrators, thus facing different levels of oxygen and temperature within a 24 h cycle. Declining oxygen levels may lead to vertically constrained habitats in euphausiids, which consequently will affect several trophic levels in the food web of the respective ecosystem. By using the regulation index (RI), the present study aimed at investigating the hypoxia tolerances of different euphausiid species from Atlantic, Pacific as well as from Polar regions. RI was calculated from 141 data sets and used to differentiate between respiration strategies using median and quartile (Q) values: low degree of oxyregulation (0.25 < RI median < 0.5); high degree of oxyregulation (0.5 < RI median < 1; Q1 > 0.25 or Q3 > 0.75); and metabolic suppression (RI median, Q1 and Q3 < 0). RI values of the polar (Euphausia superba, Thysanoessa inermis) and sub-tropical (Euphausia hanseni, Nyctiphanes capensis, and Nematoscelis megalops) species indicate a high degree of oxyregulation, whereas almost perfect oxyconformity (RI median ≈ 0; Q1 < 0 and Q3 > 0) was identified for the neritic temperate species Thysanoessa spinifera and the tropical species Euphausia lamelligera. RI values of Euphausia distinguenda and the Humboldt species Euphausia mucronata qualified these as metabolic suppressors. RI showed a significant impact of temperature on the respiration strategy of E. hanseni from oxyregulation to metabolic suppression. The species' estimated hypoxia tolerances and the degree of oxyconformity vs. oxyregulation were linked to diel vertical migration behavior and the temperature experienced during migration. The results highlight that the euphausiid species investigated have evolved various strategies to deal with different levels of oxygen, ranging from species showing a high degree of oxyconformity to strong oxyregulation. Neritic species may be more affected by hypoxia, as these are often short-distance-migrators and only adapted to a narrow range of environmental conditions.

Habitat compression and ecosystem shifts as potential links between marine heatwave and record whale entanglements

Climate change and increased variability and intensity of climate events, in combination with recovering protected species populations and highly capitalized fisheries, are posing new challenges for fisheries management. We examine socio-ecological features of the unprecedented 2014–2016 northeast Pacific marine heatwave to understand the potential causes for record numbers of whale entanglements in the central California Current crab fishery. We observed habitat compression of coastal upwelling, changes in availability of forage species (krill and anchovy), and shoreward distribution shift of foraging whales. We propose that these ecosystem changes, combined with recovering whale populations, contributed to the exacerbation of entanglements throughout the marine heatwave. In 2016, domoic acid contamination prompted an unprecedented delay in the opening of California’s Dungeness crab fishery that inadvertently intensified the spatial overlap between whales and crab fishery gear. We present a retroactive assessment of entanglements to demonstrate that cooperation of fishers, resource managers, and scientists could mitigate future entanglement risk by developing climate-ready fisheries approaches, while supporting thriving fishing communities.

Climate-driven extreme events may have strong local impacts on marine organisms and fisheries. Here the authors report increased whale entanglements in the northeast Pacific following a marine heatwave, and propose compression of coastal upwelling habitat as the potential driver.

Integrating patterns of thermal tolerance and phenotypic plasticity with population genetics to improve understanding of vulnerability to warming in a widespread copepod

Differences in population vulnerability to warming are defined by spatial patterns in thermal adaptation. These patterns may be driven by natural selection over spatial environmental gradients, but can also be shaped by gene flow, especially in marine taxa with high dispersal potential. Understanding and predicting organismal responses to warming requires disentangling the opposing effects of selection and gene flow. We begin by documenting genetic divergence of thermal tolerance and developmental phenotypic plasticity. Ten populations of the widespread copepod Acartia tonsa were collected from sites across a large thermal gradient, ranging from the Florida Keys to Northern New Brunswick, Canada (spanning over 20° latitude). Thermal performance curves (TPCs) from common garden experiments revealed local adaptation at the sampling range extremes, with thermal tolerance increasing at low latitudes and decreasing at high latitudes. The opposite pattern was observed in phenotypic plasticity, which was strongest at high latitudes. No relationship was observed between phenotypic plasticity and environmental variables. Instead, the results are consistent with the hypothesis of a trade-off between thermal tolerance and the strength of phenotypic plasticity. Over a large portion of the sampled range, however, we observed a remarkable lack of differentiation of TPCs. To examine whether this lack of divergence is the result of selection for a generalist performance curve or constraint by gene flow, we analyzed cytochrome oxidase I mtDNA sequences, which revealed four distinct genetic clades, abundant genetic diversity, and widely distributed haplotypes. Strong divergence in thermal performance within genetic clades, however, suggests that the pace of thermal adaptation can be relatively rapid. The combined insight from the laboratory physiological experiments and genetic data indicate that gene flow constrains differentiation of TPCs. This balance between gene flow and selection has implications for patterns of vulnerability to warming. Taking both genetic differentiation and phenotypic plasticity into account, our results suggest that local adaptation does not increase vulnerability to warming, and that low-latitude populations in general may be more vulnerable to predicted temperature change over the next century.

Weakened portfolio effects constrain management effectiveness for population aggregates

Population diversity can reduce temporal variability in aggregate population abundances in a process known as the portfolio effect. Portfolio effects may weaken, however, due to greater synchrony among component populations. While weakened portfolio effects have been previously documented, the consequences of reduced stability on meeting conservation goals for population aggregates that are harvested (e.g., stock aggregates in fisheries) are rarely quantified. Here, we demonstrate how changes in variability within components, synchrony among components, and population productivity interact to influence the probability of achieving an array of management objectives for Fraser River sockeye salmon: a stock aggregate of high economic, ecological, and cultural value. We first present evidence that component variability and synchrony have increased over the last two decades, consistent with a weakening portfolio effect. We then parameterize a stochastic, closed-loop model that simulates the population dynamics of each stock, the fishery that harvests the stock aggregate, and the management framework used to establish mixed-stock exploitation rates. We find that while median aggregate abundance and catch through time were relatively insensitive to greater aggregate variability, catch stability and performance metrics associated with achieving management targets generally declined as component variability and synchrony increased. A notable exception we observed is that harvest control means that scale exploitation rates based on aggregate abundance may be more effective as synchrony increases. Reductions in productivity led to broad declines in performance, but also moderated the impacts of component variability and synchrony on the proportion of component stocks above management targets and catch stability. Our results suggest that even precautionary management strategies that account for declines in productivity may underestimate risk, particularly to socioeconomic objectives, if they fail to consider changes in aggregate variability. Adequately incorporating changes in portfolio effect strength may be particularly relevant when developing recovery strategies that are robust to climate change, which is likely to increase synchrony and component variability.

Predator-mediated effects of severe drought associated with poor reproductive success of a seabird in a cross-ecosystem cascade

Despite the profound impacts of drought on terrestrial productivity in coastal arid ecosystems, only a few studies have addressed how drought can influence ecological cascades across ecosystem boundaries. In this study, we examine the consequences of rainfall pulses and drought that subsequently impact the breeding success of a threatened nocturnal seabird, the Scripps's Murrelet (Synthliboramphus scrippsi). On an island off the coast of southern California, the main cause of reduced nest success for one of their largest breeding colonies is egg predation by an endemic deer mouse (Peromyscus maniculatus elusus). Mice on the island have an opportunistic diet of primarily terrestrial sources, but drastic declines in terrestrial productivity from drought might be expected to increase their reliance on marine resources, including murrelet eggs. We compiled data on terrestrial and marine productivity between 1983 and 2013 to determine how conditions in these ecosystems influence murrelet nest success. We found that the severity of drought had the strongest negative impact on murrelet nest success. We calculated that the reduction in fecundity during drought years due to increased egg predation by mice was substantial enough to produce a declining population growth rate. Nest success was much higher under normal or high rainfall conditions, depending on whether oceanic conditions were favorable to murrelets. Therefore, the more frequent and severe drought that is projected for this region could lead to an increased risk of murrelet population decline on this island. Our study highlights the need for understanding how species interactions will change through the effects of increasing drought and altered rainfall regimes under global change.

Intertidal clams exhibit population synchrony across spatial and temporal scales

Long-term datasets can be particularly useful for parsing out factors influencing populations, yet few studies have utilized continuous datasets to quantify population dynamics in bivalve molluscs. We used dynamic factor analysis on a clam biomass dataset spanning 28 yr and five distinct regions in the southern Salish Sea to determine (1) if native intertidal clam populations exhibit synchrony and (2) what environmental covariates may be correlated with these population trends. Once covariates were accounted for, the model with the most data support included three predominant trends to describe multidecadal change in clam biomass. Intraspecific synchrony was highest among Saxidomus gigantea and Leukoma staminea populations, with no clear evidence of covariance in Clinocardium nuttallii. Specifically, we quantified a pronounced decadal decline in L. staminea and an increase in S. gigantea biomass on most beaches. No beaches showed synchrony in trends across all three species, indicating that species-specific trends (regardless of location) were more common than beach-specific trends (regardless of species). Seven environmental covariates were evaluated in their capacity to explain variability in annual mean biomass. Of these, the North Pacific Gyre Oscillation lagged 4 yr prior to the observation year was most supported by the data in the best fitting model, implying that 4 yr old clam biomass is partially determined by oceanographic processes affecting larval clams. Although results suggest large-scale density-independent factors play a role in venerid clam population dynamics, it is also likely local factors account for variability not explained by our model.

Non-stationary climate-salmon relationships in the Gulf of Alaska

Studies of climate effects on ecology often account for non-stationarity in individual physical and biological variables, but rarely allow for non-stationary relationships among variables. Here, we show that non-stationary relationships among physical and biological variables are central to understanding climate effects on salmon (Onchorynchus spp.) in the Gulf of Alaska during 1965-2012. The relative importance of two leading patterns in North Pacific climate, the Pacific Decadal Oscillation (PDO) and North Pacific Gyre Oscillation (NPGO), changed around 1988/1989 as reflected by changing correlations with leading axes of sea surface temperature variability. Simultaneously, relationships between the PDO and Gulf of Alaska environmental variables weakened, and long-standing temperature-salmon and PDO-salmon covariance declined to zero. We propose a mechanistic explanation for changing climate-salmon relationships in terms of non-stationary atmosphere-ocean interactions coinciding with changing PDO-NPGO relative importance. We also show that regression models assuming stationary climate-salmon relationships are inappropriate over the multidecadal time scale we consider. Relaxing assumptions of stationary relationships markedly improved modelling of climate effects on salmon catches and productivity. Attempts to understand the implications of changing climate patterns in other ecosystems might also be aided by the application of models that allow associations among environmental and biological variables to change over time.

Rising synchrony controls western North American ecosystems

Along the western margin of North America, the winter expression of the North Pacific High (NPH) strongly influences interannual variability in coastal upwelling, storm track position, precipitation, and river discharge. Coherence among these factors induces covariance among physical and biological processes across adjacent marine and terrestrial ecosystems. Here, we show that over the past century the degree and spatial extent of this covariance (synchrony) has substantially increased, and is coincident with rising variance in the winter NPH. Furthermore, centuries-long blue oak (Quercus douglasii) growth chronologies sensitive to the winter NPH provide robust evidence that modern levels of synchrony are among the highest observed in the context of the last 250 years. These trends may ultimately be linked to changing impacts of the El Niño Southern Oscillation on midlatitude ecosystems of North America. Such a rise in synchrony may destabilize ecosystems, expose populations to higher risks of extinction, and is thus a concern given the broad biological relevance of winter climate to biological systems.

Pelagic marine refugia and climatically sensitive areas in an eastern boundary current upwelling system

Refugia are areas relatively buffered from contemporary climate change that enable the persistence of valued physical, ecological, or sociocultural resources. Spatially identifying refugia is important for conservation and applied management. Yet the concept of refugia has not been broadly extended to marine ecosystems. Here, we analyze data from a unique and long-term (1999-2015) standardized survey of pelagic marine and anadromous species off Oregon and Washington in the northern California Current to identify such refugia. We use quantitative approaches to assess locations with high species richness and community persistence relative to local and basin-scale environmental fluctuations. We have identified a potential climate change refugial zone along the continental shelf of Washington State in the Northeastern Pacific Ocean, characterized by a species-rich community with low interannual temporal community change. This region contrasts with adjacent areas to the south and offshore that have lower species richness, and higher temporal species community change. Also, using spatially variant generalized additive mixed models, we identify areas with species compositions that are more influenced by basin-scale climatic fluctuations than others. We propose that upwelling regions with retentive topographic features, such as wide continental shelves, can function as marine refugia for pelagic fauna, whereas offshore locations are potentially more climatically sensitive and experience high temporal change in species composition. Further identification of these marine refugia using in situ data for pelagic biodiversity and climatically sensitive areas can help guide management in the face of inevitable climatically driven change.

Climate mediates the success of migration strategies in a marine predator

Individual behavioural specialisation has far-reaching effects on fitness and population persistence. Theory predicts that unconditional site fidelity, that is fidelity to a site independent of past outcome, provides a fitness advantage in unpredictable environments. However, the benefits of alternative site fidelity strategies driving intraspecific variation remain poorly understood and have not been evaluated in different environmental contexts. We show that contrary to expectation, strong and weak site fidelity strategies in migratory northern elephant seals performed similarly over 10 years, but the success of each strategy varied interannually and was strongly mediated by climate conditions. Strong fidelity facilitated stable energetic rewards and low risk, while weak fidelity facilitated high rewards and high risk. Weak fidelity outperformed strong fidelity in anomalous climate conditions, suggesting that the evolutionary benefits of site fidelity may be upended by increasing environmental variability. We highlight how individual behavioural specialisation may modulate the adaptive capacity of species to climate change.

Persistence of trophic hotspots and relation to human impacts within an upwelling marine ecosystem

Human impacts (e.g., fishing, pollution, and shipping) on pelagic ecosystems are increasing, causing concerns about stresses on marine food webs. Maintaining predator-prey relationships through protection of pelagic hotspots is crucial for conservation and management of living marine resources. Biotic components of pelagic, plankton-based, ecosystems exhibit high variability in abundance in time and space (i.e., extreme patchiness), requiring investigation of persistence of abundance across trophic levels to resolve trophic hotspots. Using a 26-yr record of indicators for primary production, secondary (zooplankton and larval fish), and tertiary (seabirds) consumers, we show distributions of trophic hotspots in the southern California Current Ecosystem result from interactions between a strong upwelling center and a productive retention zone with enhanced nutrients, which concentrate prey and predators across multiple trophic levels. Trophic hotspots also overlap with human impacts, including fisheries extraction of coastal pelagic and groundfish species, as well as intense commercial shipping traffic. Spatial overlap of trophic hotspots with fisheries and shipping increases vulnerability of the ecosystem to localized depletion of forage fish, ship strikes on marine mammals, and pollution. This study represents a critical step toward resolving pelagic areas of high conservation interest for planktonic ecosystems and may serve as a model for other ocean regions where ecosystem-based management and marine spatial planning of pelagic ecosystems is warranted.

Modeling Nonresident Seabird Foraging Distributions to Inform Ocean Zoning in Central California

Seabird aggregations at sea have been shown to be associated with concentrations of prey. Previous research identified Central California as a highly used foraging area for seabirds, with locally breeding seabirds foraging close to their colonies on Southeast Farallon Island. Herein, we focus on nonresident (i.e. non-locally breeding) seabird species off of Central California. We hypothesized that high-use foraging areas for nonresident seabirds would be influenced by oceanographic and bathymetric factors and that spatial and temporal distributions would be similar within planktivorous and generalist foraging guilds but would differ between them. With data collected by the Applied California Current Ecosystem Studies (ACCESS) partnership during cruises between April and October from 2004–2013, we developed generalized linear models to identify high-use foraging areas for each of six nonresident seabird species. The four generalist species are Phoebastria nigripes (black-footed albatross), Ardenna griseus (sooty shearwater), Ardenna creatopus (pink-footed shearwater), and Fulmarus glacialis (northern fulmar). The two planktivorous species are Phalaropus lobatus (red-necked phalarope) and Phalaropus fulicarius (red phalarope). Sea surface temperature was significant for generalist species and sea surface salinity was important for planktivorous species. The distance to the 200-m isobath was significant in five of six models, Pacific Decadal Oscillation with a 3-month lag in four models, and sea surface fluorescence, the distance to Cordell Bank, and depth in three models. We did not find statistically significant differences between distributions of individual seabird species within a foraging guild or between guilds, with the exception of the sooty shearwater. Model results for a multi-use seabird foraging area highlighted the continental shelf break, particularly within the vicinity of Cordell Bank, as the highest use areas as did Marxan prioritization. Our research methods can be implemented elsewhere to identify critical habitat that needs protection as human development pressures continue to expand to the ocean.

Climate extremes are associated with invertebrate taxonomic and functional composition in mountain lakes

Climate change is expected to increase climate variability and the occurrence of extreme climatic events, with potentially devastating effects on aquatic ecosystems. However, little is known about the role of climate extremes in structuring aquatic communities or the interplay between climate and local abiotic and biotic factors. Here, we examine the relative influence of climate and local abiotic and biotic conditions on biodiversity and community structure in lake invertebrates. We sampled aquatic invertebrates and measured environmental variables in 19 lakes throughout California, USA, to test hypotheses of the relationship between climate, local biotic and environmental conditions, and the taxonomic and functional structure of aquatic invertebrate communities. We found that, while local biotic and abiotic factors such as habitat availability and conductivity were the most consistent predictors of alpha diversity, extreme climate conditions such as maximum summer temperature and dry-season precipitation were most often associated with multivariate taxonomic and functional composition. Specifically, sites with high maximum temperatures and low dry-season precipitation housed communities containing high abundances of large predatory taxa. Furthermore, both climate dissimilarity and abiotic dissimilarity determined taxonomic turnover among sites (beta diversity). These findings suggest that while local-scale environmental variables may predict alpha diversity, climatic variability is important to consider when projecting broad-scale aquatic community responses to the extreme temperature and precipitation events that are expected for much of the world during the next century.

Humpback whale diets respond to variance in ocean climate and ecosystem conditions in the California Current

Large, migratory predators are often cited as sentinel species for ecosystem processes and climate-related changes, but their utility as indicators is dependent upon an understanding of their response to environmental variability. Documentation of the links between climate variability, ecosystem change and predator dynamics is absent for most top predators. Identifying species that may be useful indicators and elucidating these mechanistic links provides insight into current ecological dynamics and may inform predictions of future ecosystem responses to climatic change. We examine humpback whale response to environmental variability through stable isotope analysis of diet over a dynamic 20-year period (1993-2012) in the California Current System (CCS). Humpback whale diets captured two major shifts in oceanographic and ecological conditions in the CCS. Isotopic signatures reflect a diet dominated by krill during periods characterized by positive phases of the North Pacific Gyre Oscillation (NPGO), cool sea surface temperature (SST), strong upwelling and high krill biomass. In contrast, humpback whale diets are dominated by schooling fish when the NPGO is negative, SST is warmer, seasonal upwelling is delayed and anchovy and sardine populations display increased biomass and range expansion. These findings demonstrate that humpback whales trophically respond to ecosystem shifts, and as a result, their foraging behavior is a synoptic indicator of oceanographic and ecological conditions across the CCS. Multi-decadal examination of these sentinel species thus provides insight into biological consequences of interannual climate fluctuations, fundamental to advancing ecosystem predictions related to global climate change.

Estuarine fish communities respond to climate variability over both river and ocean basins

Estuaries are dynamic environments at the land-sea interface that are strongly affected by interannual climate variability. Ocean-atmosphere processes propagate into estuaries from the sea, and atmospheric processes over land propagate into estuaries from watersheds. We examined the effects of these two separate climate-driven processes on pelagic and demersal fish community structure along the salinity gradient in the San Francisco Estuary, California, USA. A 33-year data set (1980-2012) on pelagic and demersal fishes spanning the freshwater to marine regions of the estuary suggested the existence of five estuarine salinity fish guilds: limnetic (salinity = 0-1), oligohaline (salinity = 1-12), mesohaline (salinity = 6-19), polyhaline (salinity = 19-28), and euhaline (salinity = 29-32). Climatic effects propagating from the adjacent Pacific Ocean, indexed by the North Pacific Gyre Oscillation (NPGO), affected demersal and pelagic fish community structure in the euhaline and polyhaline guilds. Climatic effects propagating over land, indexed as freshwater outflow from the watershed (OUT), affected demersal and pelagic fish community structure in the oligohaline, mesohaline, polyhaline, and euhaline guilds. The effects of OUT propagated further down the estuary salinity gradient than the effects of NPGO that propagated up the estuary salinity gradient, exemplifying the role of variable freshwater outflow as an important driver of biotic communities in river-dominated estuaries. These results illustrate how unique sources of climate variability interact to drive biotic communities and, therefore, that climate change is likely to be an important driver in shaping the future trajectory of biotic communities in estuaries and other transitional habitats.

Slowing down of North Pacific climate variability and its implications for abrupt ecosystem change

Marine ecosystems are sensitive to stochastic environmental variability, with higher-amplitude, lower-frequency--i.e., "redder"--variability posing a greater threat of triggering large ecosystem changes. Here we show that fluctuations in the Pacific Decadal Oscillation (PDO) index have slowed down markedly over the observational record (1900-present), as indicated by a robust increase in autocorrelation. This "reddening" of the spectrum of climate variability is also found in regionally averaged North Pacific sea surface temperatures (SSTs), and can be at least partly explained by observed deepening of the ocean mixed layer. The progressive reddening of North Pacific climate variability has important implications for marine ecosystems. Ecosystem variables that respond linearly to climate forcing will have become prone to much larger variations over the observational record, whereas ecosystem variables that respond nonlinearly to climate forcing will have become prone to more frequent "regime shifts." Thus, slowing down of North Pacific climate variability can help explain the large magnitude and potentially the quick succession of well-known abrupt changes in North Pacific ecosystems in 1977 and 1989. When looking ahead, despite model limitations in simulating mixed layer depth (MLD) in the North Pacific, global warming is robustly expected to decrease MLD. This could potentially reverse the observed trend of slowing down of North Pacific climate variability and its effects on marine ecosystems.

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.

Changing central Pacific El Niños reduce stability of North American salmon survival rates

Pacific salmon are a dominant component of the northeast Pacific ecosystem. Their status is of concern because salmon abundance is highly variable--including protected stocks, a recently closed fishery, and actively managed fisheries that provide substantial ecosystem services. Variable ocean conditions, such as the Pacific Decadal Oscillation (PDO), have influenced these fisheries, while diminished diversity of freshwater habitats have increased variability via the portfolio effect. We address the question of how recent changes in ocean conditions will affect populations of two salmon species. Since the 1980s, El Niño Southern Oscillation (ENSO) events have been more frequently associated with central tropical Pacific warming (CPW) rather than the canonical eastern Pacific warming ENSO (EPW). CPW is linked to the North Pacific Gyre Oscillation (NPGO), whereas EPW is linked to the PDO, different indicators of northeast Pacific Ocean ecosystem productivity. Here we show that both coho and Chinook salmon survival rates along western North America indicate that the NPGO, rather than the PDO, explains salmon survival since the 1980s. The observed increase in NPGO variance in recent decades was accompanied by an increase in coherence of local survival rates of these two species, increasing salmon variability via the portfolio effect. Such increases in coherence among salmon stocks are usually attributed to controllable freshwater influences such as hatcheries and habitat degradation, but the unknown mechanism underlying the ocean climate effect identified here is not directly subject to management actions.

Shifting Effects of Ocean Conditions on Survival and Breeding Probability of a Long-Lived Seabird

With a rapidly changing climate, there is an increasing need to predict how species will respond to changes in the physical environment. One approach is to use historic data to estimate the past influence of environmental variation on important demographic parameters and then use these relationships to project the abundance of a population or species under future climate scenarios. However, as novel climate conditions emerge, novel species responses may also appear. In some systems, environmental conditions beyond the range of those observed during the course of most long-term ecological studies are already evident. Yet little attention has been given to how these novel conditions may be influencing previously established environment–species relationships. Here, we model the relationships between ocean conditions and the demography of a long-lived seabird, Brandt’s cormorant (Phalacrocorax penicillatusI), in central California and show that these relationships have changed in recent years. Beginning in 2007/2008, the response of Brandt’s cormorant, an upper trophic level predator, to ocean conditions shifted, resulting in lower than predicted survival and breeding probability. Survival was generally less variable than breeding probability and was initially best predicted by the basin-scale forcing of the El Niño Southern Oscillation rather than local ocean conditions. The shifting response of Brandt’s cormorant to ocean conditions may be just a proximate indication of altered dynamics in the food web and that important forage fish are not responding to the physical ocean environment as expected. These changing relationships have important implications for our ability to project the effects of future climate change for species and communities.

Warm oceanographic anomalies and fishing pressure drive seabird nesting north

Parallel studies of nesting colonies in Mexico and the United States show that Elegant Terns (Thalasseus elegans) have expanded from the Gulf of California Midriff Island Region into Southern California, but the expansion fluctuates from year to year. A strong inverse relationship between nesting pairs in three Southern California nesting areas [San Diego saltworks, Bolsa Chica Ecological Reserve, and Los Angeles Harbor (1991 to 2014)] and Isla Rasa in the Midriff (1980 to 2014) shows that terns migrate northward when confronting warm oceanographic anomalies (>1.0°C), which may decrease fish availability and hamper nesting success. Migration pulses are triggered by sea surface temperature anomalies localized in the Midriff and, secondarily, by reductions in the sardine population as a result of intensive fishing. This behavior is new; before year 2000, the terns stayed in the Midriff even when oceanographic conditions were adverse. Our results show that terns are responding dynamically to rapidly changing oceanographic conditions and fish availability by migrating 600 km northwest in search of more productive waters.

Assessing trade-offs to inform ecosystem-based fisheries management of forage fish

Twenty-first century conservation is centered on negotiating trade-offs between the diverse needs of people and the needs of the other species constituting coupled human-natural ecosystems. Marine forage fishes, such as sardines, anchovies, and herring, are a nexus for such trade-offs because they are both central nodes in marine food webs and targeted by fisheries. An important example is Pacific herring, Clupea pallisii in the Northeast Pacific. Herring populations are subject to two distinct fisheries: one that harvests adults and one that harvests spawned eggs. We develop stochastic, age-structured models to assess the interaction between fisheries, herring populations, and the persistence of predators reliant on herring populations. We show that egg- and adult-fishing have asymmetric effects on herring population dynamics - herring stocks can withstand higher levels of egg harvest before becoming depleted. Second, ecosystem thresholds proposed to ensure the persistence of herring predators do not necessarily pose more stringent constraints on fisheries than conventional, fishery driven harvest guidelines. Our approach provides a general template to evaluate ecosystem trade-offs between stage-specific harvest practices in relation to environmental variability, the risk of fishery closures, and the risk of exceeding ecosystem thresholds intended to ensure conservation goals are met.

Climate change. Six centuries of variability and extremes in a coupled marine-terrestrial ecosystem

Reported trends in the mean and variability of coastal upwelling in eastern boundary currents have raised concerns about the future of these highly productive and biodiverse marine ecosystems. However, the instrumental records on which these estimates are based are insufficiently long to determine whether such trends exceed preindustrial limits. In the California Current, a 576-year reconstruction of climate variables associated with winter upwelling indicates that variability increased over the latter 20th century to levels equaled only twice during the past 600 years. This modern trend in variance may be unique, because it appears to be driven by an unprecedented succession of extreme, downwelling-favorable, winter climate conditions that profoundly reduce productivity for marine predators of commercial and conservation interest.

Spatio-temporal dynamics of ocean conditions and forage taxa reveal regional structuring of seabird–prey relationships

Studies of predator–prey demographic responses and the physical drivers of such relationships are rare, yet essential for predicting future changes in the structure and dynamics of marine ecosystems. Here, we hypothesize that predator–prey relationships vary spatially in association with underlying physical ocean conditions, leading to observable changes in demographic rates, such as reproduction. To test this hypothesis, we quantified spatio-temporal variability in hydrographic conditions, krill, and forage fish to model predator (seabird) demographic responses over 18 years (1990–2007). We used principal component analysis and spatial correlation maps to assess coherence among ocean conditions, krill, and forage fish, and generalized additive models to quantify interannual variability in seabird breeding success relative to prey abundance. The first principal component of four hydrographic measurements yielded an index that partitioned “warm/weak upwelling” and “cool/strong upwelling” years. Partitioning of krill and forage fish time series among shelf and oceanic regions yielded spatially explicit indicators of prey availability. Krill abundance within the oceanic region was remarkably consistent between years, whereas krill over the shelf showed marked interannual fluctuations in relation to ocean conditions. Anchovy abundance varied on the shelf, and was greater in years of strong stratification, weak upwelling and warmer temperatures. Spatio-temporal variability of juvenile forage fish co-varied strongly with each other and with krill, but was weakly correlated with hydrographic conditions. Demographic responses between seabirds and prey availability revealed spatially variable associations indicative of the dynamic nature of “predator–habitat” relationships. Quantification of spatially explicit demographic responses, and their variability through time, demonstrate the possibility of delineating specific critical areas where the implementation of protective measures could maintain functions and productivity of central place foraging predators.