https://orcid.org/0000-0001-7784-2933
Figures
Abstract
Bull kelp, Nereocystis luetkeana, is an iconic kelp forest species of the Northeast Pacific that provides a wide range of ecosystem services to coastal marine species and society. In northern California, U.S.A., Nereocystis abundance declined sharply in 2014 and has yet to recover. While abiotic and biotic stressors were present prior to 2014, the population collapse highlights the need for a better understanding of how environmental conditions impact Nereocystis. In this study, we used a newly-developed, satellite-based dataset of bull kelp abundance, proxied by canopy cover over 20 years, to test the hypothesis that winter oceanographic conditions determine summer Nereocystis canopy cover. For the years before the collapse (1991 through 2013), wintertime ocean conditions, synthesized in a Multivariate Ocean Climate Indicator (MOCI), were indeed a good predictor of summer Nereocystis canopy cover (R2 = 0.40 to 0.87). We attribute this relationship to the effects of upwelling and/or temperature on nutrient availability. South of Point Arena, California, winter ocean conditions had slightly lower explanatory power than north of Point Arena, also reflective of spring upwelling-driven nutrient entrainment. Results suggest that the Nereocystis gametophytes and/or early sporophytes are sensitive to winter oceanographic conditions. Furthermore, environmental conditions in winter 2014 could have been used to predict the Nereocystis collapse in summer 2014, and for kelp north of Point Arena, a further decline in 2015. Importantly, environmental models do not predict changes in kelp after 2015, suggesting biotic factors suppressed kelp recovery, most likely extreme sea urchin herbivory. Conditions during winter, a season that is often overlooked in studies of biophysical interactions, are useful for predicting summer Nereocystis kelp forest canopy cover, and will be useful in supporting kelp restoration actions in California and perhaps elsewhere in the world.
Citation: García-Reyes M, Thompson SA, Rogers-Bennett L, Sydeman WJ (2022) Winter oceanographic conditions predict summer bull kelp canopy cover in northern California. PLoS ONE 17(5): e0267737. https://doi.org/10.1371/journal.pone.0267737
Editor: Eliza Heery, Smithsonian - University of Washington, UNITED STATES
Received: September 23, 2021; Accepted: April 13, 2022; Published: May 5, 2022
Copyright: © 2022 García-Reyes et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All data used in this study is publicly available; the details of where can it be accessed are included in the manuscript and its Supporting Information. Derived data used in this analysis can be found at: https://github.com/farallon-institute/Garcia-Reyes_etal_BullKelp.
Funding: M. García-Reyes and W.J. Sydeman acknowledges the partial financial support from The Angell Foundation and L. Rogers-Bennett from the National Science Foundation (Biological Oceanography Award 2023664).
Competing interests: The authors have declared that no competing interests exist.
Introduction
Canopy-forming kelps structure temperate and subpolar coastal ecosystems by providing habitat, food, and shelter to multitudes of fishes, invertebrates, seabirds, and marine mammals [1]. Due to their extensive spatial coverage and efficient CO2 uptake, kelp forests are also a key component in biogeochemical processes [2], which in turn provide a range of ecosystem services to people including direct exploitation, sustainable fisheries, recreational activities [3–5], and carbon sequestration [6]. In general, the highest rates of kelp growth occur in late winter and spring due to nutrient availability, and they are lowest in late summer and fall, when temperatures are higher and nutrient availability decreases [5]. Kelp are highly sensitive to their environment; often they are negatively impacted by higher temperatures [7], reduced nutrient availability [8], and strong wave and storm activity [9], as well as direct impact from human activities (see [5, 10] and references therein).
In the northeast Pacific, bull kelp (Nereocystis luetkeana) is the dominant canopy-forming species north of ~38°N, from northern California to Alaska. This annual kelp species is characterized by large floating pneumatocysts fastened by a single long stipe with its canopy residing mostly on the water’s surface. Nereocystis has historically formed large canopies along the central-northern California coast particularly between ~38 and 39.5°N and showed periodic short-term fluctuations in abundance (Fig 1) [11]. The Nereocystis annual life cycle involves the settlement of swimming zoospores produced by mature sporophytes, which develop into the over-wintering microscopic haploid gametophyte stage [12]. During winter, male and female gametes are released from the gametophytes, fertilize, and then develop into diploid zygotes, which grow into new sporophytes in the following spring [13]. Tiny sporophytes grow rapidly over the spring and can reach up to 25 m in length, forming the large adult kelp canopy visible on the surface in the summer. Sporulation occurs in late summer and fall when Nereocystis blades mature, but before storms and wave action dislodge the algae [14]. In lab experiments, Nereocystis has shown sensitivity to temperature, with an optimal growing temperature of ~11.9°C [15], and a reported thermal upper limit of 18–20°C [16, 17]. However, a recent experimental study suggested only partial reductions in blade growth rates at temperatures of 20°C [15]. Early life history stages (from germination to sporophyte development) have shown similar thermal tolerance below 18–20°C [17, 18]. There is limited understanding, however, on Nereocystis reproduction, dispersal, gametophyte sensitivity, and growth in relation to regional oceanography [19], particularly concerning the effect of seasonal oceanographic conditions on this annual canopy-forming kelp.
(A) Study region from 38–40°N on the northern California Coast, U.S.A. Green dots indicate kelp canopy data locations. Extensive sandy beaches are located north of Point Arena and Fort Bragg; no measurable kelp canopy occurs there. (B) Summer kelp canopy extent data (km2), 1991–2020, aggregated regionally for north (orange) and south (gray) of 39°N (Point Arena). Data derived from [34], freely available under CC BY 4.0 License at https://doi.org/10.6073/pasta/89b63c4b49b80fb839613e9d389d9902.
This lack of understanding of the relationship between survival of various life history stages and ocean conditions was made more apparent when the Nereocystis population of northern California collapsed in 2014 to less than 10% of its 2008 areal extent [20]. Further, this decline in kelp canopy has persisted (Fig 1), with some limited and spatially patchy recovery in 2021 (Rogers-Bennett personal observation). Multiple stressors were identified as present during the period of collapse [20] including: 1) a population decline of the sunflower sea star (Pycnopodia helianthoides), an important sea urchin predator [21], 2) grazing pressure from a rapid increase in purple sea urchin densities [20], and 3) potential thermal stress caused by an intense multi-year marine heatwave (MHW) that impacted a broad geographic region in the northeastern Pacific [22] and reached the California coast in summer 2014 [23]. Increased urchin density is a well-documented factor in kelp forest collapse around the globe [10, 24], yet in northern California, purple urchin densities in 2014 were not anomalously high (see [11]). Rogers-Bennett and Catton [20] as well as McPherson et al. [11] noted that kelp recovery has been hindered by the unprecedented increases in purple urchin densities in northern California.
The role that the 2014–2016 MHW [23] played on the Nereocystis decline is poorly understood, partly because long-term biophysical relationships between Nereocystis and the environment have not been adequately studied. Extreme temperatures are a known stressor for kelps in general [5]. However, sea surface temperatures during the 2014 MHW only sporadically reached magnitudes of 17°C around Point Arena, California (Fig 2 in [25]), which is below the documented Nereocystis thermal tolerance limit of 18–20°C [15–18]. Moreover, despite the fact that the MHW impacted the entire North American west coast [23], kelp forests only collapsed in some areas of California [21, 26, 27]. McPherson et al. [11] suggest that the high summer temperatures and low spring nutrient content in the ocean water may have been a key player in the 2014 Nereocystis collapse.
(A) Linear regression between central California MOCI (CenCal MOCI) and kelp canopy extent north of Point Arena, 1991–2013, red-dotted lines indicate the confident intervals; R2 indicates the explained variance and the color indicates the year. (B) Time series of summer kelp canopy extent north of Point Arena: measured data are in red, modeled pre-collapse data are black, and predicted data for collapsed years are shown in grey asterisks (*). Modeled and predicted kelp uses the linear regression in (A). (C) Same as (A) but for winter southern California MOCI (SoCal MOCI). (D) Same as (B) but model and predicted data used the linear regression with southern California MOCI shown in (C). Regressions details are in Table 2.
Importantly, summer kelp abundance may be impacted by oceanographic conditions during earlier seasons, thereby affecting early life history stages. Of particular interest are wintertime conditions, which in the California Current are highly variable [28] and are known to play a large role in the productivity of the pelagic ecosystem in subsequent seasons [29, 30]. For example, groundfish growth and seabird survival in the central California Current are negatively impacted by warm winters and weak late-winter/early-spring upwelling [28, 31]. Winter ocean conditions may explain variation in Nereocystis canopy cover in this region, impacting early life history stages, but to date these ideas have yet to be investigated.
In this study, we use a new 20-year annual-scale dataset of kelp canopy cover to examine the effects and predictive capacity of seasonal oceanographic conditions on Nereocystis canopy cover in northern California. In particular, we test the hypothesis that winter oceanographic conditions determine summer Nereocystis canopy cover, and that conditions in winter 2014 played a role in the dramatic decline observed later that year. To test this hypothesis, we model the relationship between fall, winter, spring, and summer ocean conditions, synthesized by the Multivariate Ocean Climate Indicator (MOCI) [32, 33] on Nereocystis canopy cover variability in summer. We also examine potential mechanistic linkages by considering how temperature and upwelling-driven nutrient availability may affect variability in kelp abundance. Finally, we suggest that biophysical models may inform kelp restoration efforts in the future.
Data and methods
The canopy data [34] is based on Landsat 5, 7, and 8 imagery providing the area (m2) covered by bull and giant kelp canopy along the coast. The spatial resolution is 30 x 30 m per pixel extending from Baja California Sur, Mexico, to Washington State, U.S.A., and it is resolved quarterly (seasonal) from March 1984 to December 2020. While kelp canopy cover measured in this way cannot be explicitly identified as Nereocystis, it is assumed that the majority of the observed canopy north of 38°N is indeed Nereocystis as it is the dominant, canopy-forming species in the region. In contrast, giant kelp (Macrocystis pyrifera) is dominant between San Francisco Bay and northern Baja California, but north of San Francisco, its presence is limited to small quantities in sheltered bays as it is more susceptible to wave forcing. Although the presence of perennial giant kelp (or other canopy-forming species) could lead to an overestimation of Nereocystis in our satellite imagery dataset, its extent is small in comparison, as illustrated by the difference in canopy cover between winter and summer/fall (S1 Fig).
We extracted and use Nereocystis data from summer (July–September), 1991–2020, for the coastal region between 38°N and 40°N, when the Nereocystis canopy has reached its annual peak in aerial extent (S1 Fig). Our study was also divided into two regions: north and south of Point Arena (39°N); Nereocystis canopy cover differs between these regions (S1 Fig). Importantly, Point Arena marks a division in oceanographic conditions due to a combination of topography and change in coastal orientation [35]. This results in the strongest upwelling in the California Current [36] downstream (south) of Point Arena, where upwelling plumes of cold and nutrient-rich water dominate ocean conditions [37, 38]. Coincidentally, a large stretch of sandy beach on the north side of Point Arena is devoid of kelp, producing a natural break in the kelp canopy data.
To estimate the total extent of the kelp canopy within a region required aggregating the coverage value of all pixels, however, we cannot distinguish between a pixel with no data due to obscuring cloud cover and a pixel with zero kelp cover. To ensure we did not count missing data as zeros we imposed thresholds of ‘available’ data to do these aggregations. First, we divided the region into bands of 0.1-degree latitude and summed the available data for each summer in each band. For each band, we calculated the number of pixels containing data each year and then the median number of pixels with data across all years in the data set. In any year, if a band had < 90% of the median number of pixels, it was assigned as a missing value. In this way, we accounted for areas with missing pixels due to data quality or cloud cover. Next, we summed the data in each region, north and south of Point Arena (referred to also as northern and southern kelp, respectively). For this we also imposed a threshold of only 10% of missing data, otherwise the summer was assigned a missing value. Years with missing data due to these gaps were 1993, 1996, 1998, 2010, and 2012 for both regions, and 1992 and 1996 for the northern region only.
To test that winter conditions influenced summer Nereocystis canopy, and also if winter is the only or most relevant season, we examined summer kelp co-variance with seasonal oceanographic conditions, lagged (fall, winter, and spring) and not-lagged (summer). We used the MOCI, which tracks the main mode of variability in ocean conditions along three regions (northern: north of 38°N, central: 34.5–38°N, and southern: south of 34.5°N) of the California coast [32, 33]. MOCI is calculated as the first principal component of sea and air surface temperature, sea level pressure, and alongshore winds from NOAA buoys over the continental shelf, the Bakun upwelling index [39], sea level from NOAA shore stations, and the climate oscillations: Multivariate El Niño Southern Oscillation (ENSO) Index (MEI)[40], Pacific Decadal Oscillation (PDO)[41], and North Pacific Gyre Oscillation index (NPGO)[42], see S2 Fig. Seasonal MOCI values used in the analyses were: fall (OND: October, November, and December; lagged for the year previous to the summer kelp canopy data (labeled fall-1)), winter (JFM: January, February, and March), spring (AMJ: April, May, and June), and summer (JAS: July, August, and September).
To gain insight into how ocean conditions influence kelp growth, we investigated two potential mechanisms: thermal stress and nutrient availability [5]. To test the effect of thermal stress, we used sea surface temperature (SST) data from NOAA buoy 46014 off Point Arena (labeled N14) and buoy 46013 south off Jenner (N13, see map in S2 Fig). Note that SST is already included in MOCI, but here we test a specific mechanism not possible by using MOCI alone. In this region, nutrients become available to the coastal areas through mixing and upwelling, and nutrient (nitrate) concentration has a tight inverse relationship with water temperature [36, 43]. We used two environmental variables to track nutrient availability for kelp: SST for overall nutrient content, and the Biologically Effective Upwelling Transport Index (BEUTI), an oceanographic index that tracks nutrient input at the coast driven by coastal upwelling [36]. We analyzed this index at 37°N, 39°N, and 41°N as it shows differences in values along the northern California coast associated with centers of upwelling. The monthly values of these indices were averaged seasonally to coincide with MOCI seasons. Time series of MOCI, SST, and BEUTI indices and data details can be found in S3–S10 Figs.
The shared period of kelp and environmental data was 1991–2020, which we divided into kelp pre-collapse (1991–2013) and post-collapse (2014–2020) periods. We investigated the relationships between summer kelp canopy cover and MOCI using univariate regression models, and then between kelp canopy and individual environmental variables using linear regressions; we assumed significance at p < 0.05 and used AIC values (Akaike Information Criterion) [44] to select the best models. We cross-validated each model and its predictability skills by running the selected model and removing each year at a time; average R2 for all “boot-strapped” models was taken as precision. Predictability errors were estimated in the same way: the model was run without one year at a time, and that year was predicted using that model. The error in prediction was calculated as absolute (measured—predicted)/measured, and expressed as a percentage. Finally, full pre-collapse models (1991–2013) were used to predict kelp values for 2014 onward. All analyses were performed in Python 3.7 using the scikit-learn package. Code is available in GitHub (https://github.com/farallon-institute/Garcia-Reyes_etal_BullKelp).
Results
Kelp canopy south of Point Arena covered ~5x more area than north of Point Arena (Figs 1 and S1), and though interannual covariability between regions was substantial (⍴ = 0.87, p < 0.001, ~50% of shared variance), we considered these regions separately. No statistically significant linear trends were observed in summer kelp canopy time series from 1991 to 2013, although it is worth noting the reduced number of data points due to gaps and the large variability in canopy extent. During the pre-collapse period, however, environmental variables showed statistically significant trends in winter consistent with increasing upwelling and decreasing water temperature (S1 Table). There were not significant trends when the years 2014–2020 were included.
Univariate linear regressions between seasonal MOCI and summer kelp canopy cover resulted in strong significant relationships for winter and spring MOCI, with R2 values ranging from 0.40 to 0.87 (Tables 1 and 2 shows regression details for selected models). Weak relationships were found with the previous year’s fall, and no significant relationships were established for the summer. The highest R2 values were found for winter southern and central California MOCI for kelp north of Point Arena (Table 1 and Fig 2) and for northern and central California MOCI for kelp south of Point Arena; surprisingly, the strongest correlation was northern kelp with southern California MOCI, but only slightly better. Spring relationships were slightly weaker than winter for kelp south of Point Arena and distinctly lower in the north, although spring R2 values were similar in magnitude for both regions. Univariate models for strong winter MOCI relationships are shown in Table 2 and illustrated in Figs 2 and 3 and S11. We did not include spring, since the correlation between winter and spring MOCI was high for the north and central California regions (S2 Table). Averaged R2 values from the cross-validation and predictability errors analysis for these models are reported in Tables 2 and 3, respectively. Cross-validated R2 values were consistent with, if not equal to, those of the model including all years (1991–2013). Average predictability error varied from 11% to 36%, and we note that there was large variability in the predictability for each year (8 to 52%). Also notable was that despite higher R2 values for models of kelp north of Point Arena, they had larger predictability errors than those for kelp south of Point Arena, especially with northern and central California MOCI.
Same as Fig 2 for regressions between kelp canopy south of Point Arena and northern California MOCI (NorCal MOCI; (A) and (B)) and central California MOCI, (CenCal MOCI; (C) and (D)).
Nutrient availability
Linear regressions of summer kelp canopy with the environmental variables reflecting potential nutrient availability showed strong relationships with winter temperature, similar to those with winter MOCI, but relationships were weaker for BEUTI (Tables 1 and 2 and Fig 4). South of Point Arena, spring relationships were significant and comparable to those of winter MOCI for BEUTI at 37°N (Table 1). In both regions, spring relationships with BEUTI were better than those for SST. In the summer, these variables did not have significant relationships with Nereocystis abundance, indicating notably that kelp canopy coverage is not related to concurrent changes in temperature or upwelling-driven nutrient entrainment. Additionally, no significant relationships were found for previous year fall environmental conditions. Univariate models for strong relationships are shown in Table 2: winter SST at buoy N14 (similar to N13) was the strongest in both regions, and for south of Point Arena, also the model for spring BEUTI 37°N.
Linear regression for the pre-collapse period 1991–2013, red-dotted lines indicate the confident intervals; R2 indicates the explained variance and the color indicates the year. (A) Summer kelp canopy extent north of Point Arena and winter SSTN14; summer kelp canopy extent south of Point Arena and (B) winter SSTN14, (C) spring BEUTI37N. Regressions details are in Table 2. SSTN14: sea surface temperature at buoy N14, BEUTI37N: BEUTI index at 37°N.
Collapse and post-collapse years
Nereocystis canopy cover declined sharply beginning in summer 2014, declined further in 2015, and has shown little recovery as of 2020. It is worth noting that 2014 was not the lowest value of canopy extent in either region before 2015; 2014 was the 4th lowest year in the north, and 6th in the south. However, before 2014, kelp north of Point Arena showed resilience after occasional years with extremely low canopy cover years. In contrast, south of Point Arena, Nereocystis abundance has historically been relatively high and more stable before 2015.
To explore the role of oceanographic conditions on the 2014 Nereocystis population collapse and subsequent years’ lack of recovery, we used the selected MOCI-kelp regression models in Table 2 to predict kelp canopy extent values for years 2014–2020 (Table 3). For kelp north of Point Arena, the model including winter southern California MOCI predicted the 2014 collapse with very low error; for southern kelp, the winter central California MOCI model predicted the decrease observed in 2014. Other models had a larger predictive error. After the 2014 collapse (from 2015 on), the models had poor predictive powers and the predictive errors became much larger. We found that for northern kelp, models underestimated kelp in 2016, and also in 2015 for southern California MOCI, as the environmental conditions in those winters were highly unfavorable due to the marine heatwave. In the south, kelp models showed a further decrease in 2015 as well, but they largely overestimated the kelp canopy coverage. In 2017, models predicted a recovery that did not materialize in the field; a further decrease of kelp canopy was observed instead. By 2020, however, winter models and measured kelp extent for both regions showed a slight recovery, although modeled values still overestimated the measured values by at least an order of magnitude. This suggests that in the post-collapse time period, something other than oceanographic conditions impeded kelp recovery.
Discussion
Nereocystis’ relationship with oceanographic conditions
We investigated lagged relationships between northern California Nereocystis canopy cover during its peak season in summer and oceanographic conditions in the preceding seasons. We focused on the period 1991–2013, prior to the long-term decline in kelp abundance that began in 2014, to understand environmental drivers of abundance. As hypothesized, winter oceanographic conditions, tracked by MOCI, predicted kelp canopy in the summer. While there was some covariance in canopy variability in regions north and south of Point Arena (39°N), we found differences in the potential oceanographic mechanisms explaining kelp canopy extent. Remarkably, north of Point Arena, winter ocean conditions alone explained up to 87% of the interannual variability in kelp canopy cover. South of Point Arena, where kelp canopy is greater, winter conditions explained 57% of the variance. We surmised these relationships are driven by upwelling and temperature effects on nutrients.
Nereocystis is an annual species, and as such, all life history stages are likewise exposed to seasonal oceanographic conditions through the year [17]. We showed that in northern California, Nereocystis is most sensitive to oceanographic changes in winter, during its gametophyte life history stage. Of the variables considered, winter sea surface temperature (SSTN14), similar to MOCI, showed the highest R2 values, lowest AIC, and was able to predict summer canopy cover later in the year. Lethal impacts of warm temperatures on the Nereocystis gametophyte stage have been previously reported by Vadas [18] and Muth et al. [17] at 18–20°C. However, during the study period 1991–2020, seasonal or daily mean winter SST did not reach lethal limits (S6–S7 and S12 Figs), and indeed the maximum daily SST value encountered in winter was only 14.9°C. This suggests that mechanistically, thermal stress of gametophytes does not explain the correlation between ocean conditions and kelp canopy cover. However, temperature is also an indicator of nutrient availability. In the central California Current, temperature has a strong negative relationship with nitrate concentration [36, 43], and ocean water becomes nitrate limited at temperatures around 13°C. The importance of nutrients to kelp is well established [5, 45]. As spring approaches, nutrients are essential for gametogenesis and sporophyte production [46]. It is therefore likely that nutrient limitation associated with higher, but not lethal, temperatures (and perhaps weaker upwelling) could have led to the initial depletion of the Nereocystis population in northern California in early 2014. This interpretation is consistent with previous work that shows winter conditions have a large impact on the benthic and pelagic ecosystems of the central California Current (e.g., [28, 29, 33]).
Nutrient-related mechanisms help to explain the surprising result that the best predictor for kelp canopy extent north of Point Arena was the southern California MOCI, and not the central or northern MOCIs. All MOCIs include similar indicators: temperature, upwelling, and regional climate indices, but each regional MOCI represents a different ‘combination’ of variables according to the region’s variability. The central and northern California MOCI show greater influence from upwelling-related variables than does the southern California MOCI, for which temperature variables show greater influence [33]. The strong relationship of northern kelp with SSTN14 also supports the idea that the temperature-driven southern California MOCI is the best predictor of canopy cover.
As with the northern area, kelp canopy south of Point Arena shows a close association with winter conditions (in this case the northern California MOCI). There were also good relationships between kelp and winter SST, and the spring nutrient entrainment index (BEUTI37N). This suggests that in addition to winter temperatures, springtime upwelling-driven nutrient availability, as reflected by BEUTI, plays an important role in determining kelp canopy cover. This is not surprising as upwelling is the dominant oceanographic process in the region south of Point Arena due to the changes in coastal orientation there that leads to very strong upwelling-favorable winds [35, 36, 38]. The explained variance for spring BEUTI is slightly lower than winter SST, suggesting that while spring nutrient input to the coastal zone by upwelling is important to kelp growth, it is not as limiting as winter temperatures and associated winter nutrient availability. This is likely due to variable, yet prevalent entrainment of nutrients in spring due to upwelling (see S10 Fig), while during winter, the large SST variability could lead to very limited nutrient conditions (S6 and S7 Figs). Finally, the lack of significant relationships between summer kelp canopy extent and summer environmental variables suggests that by summer, Nereocystis abundance as measured by canopy cover has already been set by previous oceanographic conditions.
Nereocystis canopy predictability
We established that oceanographic conditions in winter predict summer Nereocystis canopy cover from 1991–2013, with precision (predictability error) of 11 to 36% (Table 3). We then examined how ocean conditions predicted the collapse of the kelp forest in 2014, and particularly 2015, which was unprecedented in the observational record. The work presented here shows that regional ocean conditions in winter 2014, particularly temperature, were strongly unfavorable and could have been sufficient to cause the decline in kelp canopy extent observed by summer 2014, although more work is needed to elucidate the mechanism. McPherson et al. [11] also suggested that spring and summer environmental conditions have limited capability (explaining 28.3% of the variance) predicting the Nereocystis decline in 2014.
The warm conditions of the marine heatwave in winter 2015 likely contributed to the lack of kelp recovery. The low kelp canopy extent observed in 2015 north of Point Arena was accurately predicted by the winter environmental conditions. However, south of Point Arena, larger kelp canopy extent was predicted by winter conditions than was observed (Fig 3). This indicates that oceanographic conditions probably were not the only cause of continued decline in 2015 south of Point Arena. After 2015, as we have shown, ocean conditions no longer predicted kelp canopy cover (Fig 3). The models predicted a modest recovery of canopy extent beginning in 2017 as the MHW abated [23], but observed canopy actually declined further. In addition to the intense herbivory by large numbers of sea urchins, it is possible that cumulative conditions (i.e., 2+ years of unfavorable winter warming) could have had additional negative impacts on the kelp population, as these persisting poor conditions extended beyond the full Nereocystis annual life history cycle. Such an extended time of poor environmental conditions was unprecedented. It is worth noting that 2018–2020 changes in observed kelp canopy mirror the changes predicted by the models (declining in 2019 and increasing in 2020), although at a much-dampened scale due to increased herbivory. This suggests that Nereocystis is responding to a certain extent to winter oceanographic conditions. Grazing pressure due to increased sea urchin population that established while Nereocystis was decimated is the most likely explanation for the lack of kelp recovery [11, 20]. Further modeling efforts that include biotic factors are needed to help predict kelp extent in the post-collapse sea urchin barrens period. The models presented here, however, show the importance of oceanographic conditions to Nereocystis in the early life history stages.
Implications for restoration
Kelp restoration planning would benefit from knowing which regions and years are best for restoration. Spatial decisions about where to conduct kelp restoration might also be improved by fine scale temperature information to aid in selecting cold water “climate refugia” [47]. Knowing which years are predicted to be favorable for the early life history stages of kelps so that various kelp restoration measures could be enacted with a higher likelihood of success. If oceanographic conditions signal a poor kelp year then the expense and effort of restoration measures, such as sea urchin control, might best be spent in future years with environmental conditions more favorable for kelp. Likewise, kelp restoration methods such as setting out twine lines that have been seeded with kelp [48] or green gravel seeding methods [49] would be best done in years when the oceanographic conditions are favorable for young kelp. In poor years, when warm winters are expected for the California coasts, proactive steps to protect and conserve Nereocystis could be taken. There is potential for successful restoration of local Nereocystis populations using an informed application of seeding and urchin control methods in careful consideration of environmental conditions. Our work highlights the need to consider complex oceanographic factors in different regions at seasonal scales to inform restoration and conservation efforts. Looking ahead, we suggest that winter MOCI may be particularly useful for predicting good years for kelp recruitment and survival in northern California, and could also be used to further investigate the influence of ocean conditions on other kelp such as Macrocystis and canopy forming kelps worldwide.
Supporting information
S1 Fig. Average canopy extent by season by region.
Lines indicate the standard deviation around the mean canopy extent.
https://doi.org/10.1371/journal.pone.0267737.s001
(PNG)
S2 Fig. Location of regions and data included in MOCI.
MOCI (Multivariate Ocean Oscillation Indicator) is a synthesized indicator of main variability mode of oceanographic conditions in southern, central and northern California, found at http://www.faralloninstitute.org/moci.
https://doi.org/10.1371/journal.pone.0267737.s002
(PNG)
S3 Fig. Northern California MOCI.
Time series of winter, spring summer and fall values of Northern California MOCI (NorCal MOCI). Note that fall is not lagged in this and following plots as it is on the analysis with kelp canopy.
https://doi.org/10.1371/journal.pone.0267737.s003
(PNG)
S4 Fig. Central California MOCI.
Time series of winter, spring summer and fall values of Central California MOCI (CenCal MOCI).
https://doi.org/10.1371/journal.pone.0267737.s004
(PNG)
S5 Fig. Southern California MOCI.
Time series of winter, spring, summer and fall values of Southern California MOCI (SoCal MOCI).
https://doi.org/10.1371/journal.pone.0267737.s005
(PNG)
S6 Fig. Sea Surface Temperature at buoy N14.
Time series of winter, spring, summer and fall averages of SST for buoy 46014 (SSTN14), located at 39.23°N 123.97°W. Data from https://www.ndbc.noaa.gov/. Gaps in the data have been filled using reanalysis data from the NOAA’s Optimal Interpolation SST dataset and neighboring buoys.
https://doi.org/10.1371/journal.pone.0267737.s006
(PNG)
S7 Fig. Sea Surface Temperature at buoy N13.
Time series of winter, spring, summer and fall averages of SST for buoy 46013 (SSTN13), located at 38.25°N 123.30°W. Data from https://www.ndbc.noaa.gov/. Gaps in the data have been filled using reanalysis data from the NOAA’s Optimal Interpolation SST dataset and neighboring buoys.
https://doi.org/10.1371/journal.pone.0267737.s007
(PNG)
S8 Fig. BEUTI at 41°N.
Time series of winter, spring, summer and fall averages of Biologically Effective Upwelling Transport Index (BEUTI, indicator of nutrients influx to the surface layer, integrating upwelling and temperature) at 41°N. Data from: https://oceanview.pfeg.noaa.gov/products/upwelling/cutibeuti.
https://doi.org/10.1371/journal.pone.0267737.s008
(PNG)
S9 Fig. BEUTI at 39°N.
Time series of winter, spring, summer and fall averages of Biologically Effective Upwelling Transport Index (BEUTI, indicator of nutrients influx to the surface layer, integrating upwelling and temperature) at 39°N. Data from: https://oceanview.pfeg.noaa.gov/products/upwelling/cutibeuti.
https://doi.org/10.1371/journal.pone.0267737.s009
(PNG)
S10 Fig. BEUTI at 37°N.
Time series of winter, spring, summer and fall averages of Biologically Effective Upwelling Transport Index (BEUTI, indicator of nutrients influx to the surface layer, integrating upwelling and temperature) at 37°N. Data from: https://oceanview.pfeg.noaa.gov/products/upwelling/cutibeuti.
https://doi.org/10.1371/journal.pone.0267737.s010
(PNG)
S11 Fig. Linear regression between summer canopy coverage south of Point Arena and winter northern California MOCI.
(A) Linear regression and scatter plot between winter northern California MOCI (NorCal MOCI) and kelp canopy extent south of Point Arena, 1991–2013, red-dotted lines indicate the confident intervals; R2 indicates the explained variance and the color indicates the year (B) Time series of summer kelp canopy extent south of Point Arena: measured data are in red, modeled pre-collapse data are black, and predicted data for collapsed years are shown in grey asterisks (*). Modeled and predicted kelp uses the linear regression in (A).
https://doi.org/10.1371/journal.pone.0267737.s011
(PNG)
S12 Fig. Histogram of daily SST for N14.
Daily SST values during winter (January-March) for buoy N14 for the period of study 1991–2020.
https://doi.org/10.1371/journal.pone.0267737.s012
(PNG)
S1 Table. Linear trends for environmental variables in winter.
Only significant p<0.05 values shown (except BEUTI at 39°N), comparing two periods: pre-collapse (1991–2013), and entire period (1991–2020). Color indicates the sign of the trend.
https://doi.org/10.1371/journal.pone.0267737.s013
(PDF)
S2 Table. Ranked correlations of environmental variables between seasons.
1 season lag in the first three columns and 2 seasons lag in the last column. Only statistically significant correlations p<0.05 are shown.
https://doi.org/10.1371/journal.pone.0267737.s014
(PDF)
References
- 1. Teagle H, Hawkins SJ, Moore PJ, Smale DA. The role of kelp species as biogenic habitat formers in coastal marine ecosystems. J Exper Mar Biol Ecol. 2017 Jul 1;492:81–98.
- 2. Carney LT, Waaland JR, Klinger T, Ewing K. Restoration of the bull kelp Nereocystis luetkeana in nearshore rocky habitats. Mar Ecol Prog Ser. 2005 Nov 4;302:49–61.
- 3. Vásquez JA, Zuñiga S, Tala F, Piaget N, Rodríguez DC, Vega JA. Economic valuation of kelp forests in northern Chile: values of goods and services of the ecosystem. J Appl Phycol. 2014 Apr;26(2):1081–8.
- 4. Bertocci I, Araújo R, Oliveira P, Sousa‐Pinto I. Potential effects of kelp species on local fisheries. J Appl Ecol. 2015 Oct;52(5):1216–26.
- 5.
Wernberg T, Krumhansl K, Filbee-Dexter K, Pedersen MF. Status and trends for the world’s kelp forests. In World seas: An environmental evaluation 2019 Jan 1 (pp. 57–78). Academic Press.
- 6. Filbee-Dexter K, Wernberg T. Substantial blue carbon in overlooked Australian kelp forests. Sci Rep. 2020 Jul 23;10(1):1–6. pmid:31913322
- 7. Simonson EJ, Scheibling RE, Metaxas A. Kelp in hot water: I. Warming seawater temperature induces weakening and loss of kelp tissue. Mar Ecol Prog Ser. 2015 Oct 14;537:89–104.
- 8. Schmid M, Fernández PA, Gaitán-Espitia JD, Virtue P, Leal PP, Revill AT, et al. Stress due to low nitrate availability reduces the biochemical acclimation potential of the giant kelp Macrocystis pyrifera to high temperature. Algal Res. 2020 May 1;47:101895.
- 9. Cavanaugh KC, Siegel DA, Reed DC, Dennison PE. Environmental controls of giant-kelp biomass in the Santa Barbara Channel, California. Mar Ecol Prog Ser. 2011 May 16;429:1–7.
- 10. Steneck RS, Graham MH, Bourque BJ, Corbett D, Erlandson JM, Estes JA, et al. Kelp forest ecosystems: biodiversity, stability, resilience and future. Env Cons. 2002 Dec;29(4):436–59.
- 11. McPherson ML, Finger DJ, Houskeeper HF, Bell TW, Carr MH, Rogers-Bennett L, et al. Large-scale shift in the structure of a kelp forest ecosystem co-occurs with an epizootic and marine heatwave. Commun Biol. 2021 Mar 5;4(1):1–9. pmid:33398033
- 12. John BD. Alternation of generations in algae: its complexity, maintenance and evolution. Biol Rev. 1994 Aug;69(3):275–91.
- 13. Dobkowski KA, Flanagan KD, Nordstrom JR. Factors influencing recruitment and appearance of bull kelp, Nereocystis luetkeana (phylum Ochrophyta). J Phycol. 2019 Feb;55(1):236–44. pmid:30565227
- 14. Bodkin JL. Fish assemblages in Macrocystis and Nereocystis kelp forests off central California. Fish Bull. 1986 Oct 1;84(4):799–808.
- 15. Supratya VP, Coleman LJ, Martone PT. Elevated temperature affects phenotypic plasticity in the bull kelp (Nereocystis luetkeana, Phaeophyceae). J Phycol. 2020 Dec;56(6):1534–41. pmid:32666523
- 16. Lüning K, Freshwater W. Temperature tolerance of Northeast Pacific marine algae 1. J Phycol. 1988 Sep;24(3):310–5.
- 17. Muth AF, Graham MH, Lane CE, Harley CD. Recruitment tolerance to increased temperature present across multiple kelp clades. Ecology. 2019 Mar;100(3):e02594. pmid:30615200
- 18. Vadas RL. Ecological implications of culture studies on Nereocystis luetkeana. J Phycol. 1972 Jun;8(2):196–203. pmid:30615200
- 19. Castorani MC, Reed DC, Raimondi PT, Alberto F, Bell TW, Cavanaugh KC, et al. Fluctuations in population fecundity drive variation in demographic connectivity and metapopulation dynamics. Proc R Soc B: Biol Sci. 2017 Jan 25;284(1847):20162086. pmid:28123088
- 20. Rogers-Bennett L, Catton CA. Marine heat wave and multiple stressors tip bull kelp forest to sea urchin barrens. Sci Rep. 2019 Oct 21;9(1):1–9. pmid:30626917
- 21. Hamilton SL, Bell TW, Watson JR, Grorud‐Colvert KA, Menge BA. Remote sensing: generation of long‐term kelp bed data sets for evaluation of impacts of climatic variation. Ecology. 2020 Jul;101(7):e03031. pmid:32108936
- 22. Bond NA, Cronin MF, Freeland H, Mantua N. Causes and impacts of the 2014 warm anomaly in the NE Pacific. Geophys Res Lett. 2015 May 16;42(9):3414–20.
- 23. Gentemann CL, Fewings MR, García‐Reyes M. Satellite sea surface temperatures along the West Coast of the United States during the 2014–2016 northeast Pacific marine heat wave. Geophys Res Lett. 2017 Jan 16;44(1):312–9.
- 24. Filbee-Dexter K, Scheibling RE. Sea urchin barrens as alternative stable states of collapsed kelp ecosystems. Mar Ecol Prog Ser. 2014 Jan 9;495:1–25.
- 25. Sanford E, Sones JL, García-Reyes M, Goddard JH, Largier JL. Widespread shifts in the coastal biota of northern California during the 2014–2016 marine heatwaves. Sci Rep. 2019 Mar 12;9(1):1–4. pmid:30626917
- 26. Cavanaugh KC, Reed DC, Bell TW, Castorani MC, Beas-Luna R. Spatial variability in the resistance and resilience of giant kelp in southern and Baja California to a multiyear heatwave. Front Mar Sci. 2019 Jul 23;6:413.
- 27. Beas‐Luna R, Micheli F, Woodson CB, Carr M, Malone D, Torre J, et al. Geographic variation in responses of kelp forest communities of the California Current to recent climatic changes. Glob Change Biol. 2020 Nov;26(11):6457–73. pmid:32902090
- 28. Black BA, Schroeder ID, Sydeman WJ, Bograd SJ, Wells BK, Schwing FB. Winter and summer upwelling modes and their biological importance in the California Current Ecosystem. Glob Change Biol. 2011 Aug;17(8):2536–45.
- 29. Schroeder ID, Sydeman WJ, Sarkar N, Thompson SA, Bograd SJ, Schwing FB. Winter pre-conditioning of seabird phenology in the California Current. Mar Ecol Prog Ser. 2009 Oct 30;393:211–23.
- 30. Samhouri JF, Andrews KS, Fay G, Harvey CJ, Hazen EL, Hennessey SM, et al. Defining ecosystem thresholds for human activities and environmental pressures in the California Current. Ecosphere. 2017 Jun;8(6):e01860.
- 31. García-Reyes M, Sydeman WJ, Thompson SA, Black BA, Rykaczewski RR, Thayer JA, et al. Integrated assessment of wind effects on Central California’s pelagic ecosystem. Ecosystems. 2013 Aug;16(5):722–35.
- 32. Sydeman WJ, Thompson SA, García-Reyes M, Kahru M, Peterson WT, Largier JL. Multivariate ocean-climate indicators (MOCI) for the central California Current: Environmental change, 1990–2010. Prog Oceanogr. 2014 Jan 1;120:352–69.
- 33. García-Reyes M, Sydeman WJ. California multivariate ocean climate indicator (MOCI) and marine ecosystem dynamics. Ecol Ind. 2017 Jan 1;72:521–9.
- 34.
Santa Barbara Coastal LTER, Bell, T, K. Cavanaugh, D. Siegel. SBC LTER: Time series of quarterly NetCDF files of kelp biomass in the canopy from Landsat 5, 7 and 8, since 1984 (ongoing) ver 14; 2021. Environmental Data Initiative. Available, under CC BY 4.0 License, from: https://doi.org/10.6073/pasta/89b63c4b49b80fb839613e9d389d9902. (Accessed April 27, 2021).
- 35. Winant CD, Dorman CE, Friehe CA, Beardsley RC. The marine layer off northern California: An example of supercritical channel flow. J Atmos Sci. 1988 Dec 1;45(23):3588–605.
- 36. Jacox MG, Edwards CA, Hazen EL, Bograd SJ. Coastal upwelling revisited: Ekman, Bakun, and improved upwelling indices for the US West Coast. J Geophys Res. 2018 Oct;123(10):7332–50.
- 37. Largier JL, Lawrence CA, Roughan M, Kaplan DM, Dever EP, Dorman CE,. WEST: A northern California study of the role of wind-driven transport in the productivity of coastal plankton communities. Deep Sea Res Pt II. 2006 Dec 1;53(25–26):2833–49.
- 38. García‐Reyes M, Largier JL. Seasonality of coastal upwelling off central and northern California: New insights, including temporal and spatial variability. J Geophys Res. 2012 Mar;117(C3).
- 39.
Bakun A. Coastal upwelling indices, west coast of North America, 1946–71. 1973 NOAA Technical Rprt. NMFS SSRF-671.
- 40.
Wolter K, Timlin MS. Monitoring ENSO in COADS with a seasonally adjusted principal. In Proc. of the 17th Climate Diagnostics Workshop, Norman, OK, NOAA/NMC/CAC, NSSL, Oklahoma Clim. Survey, CIMMS and the School of Meteor., Univ. of Oklahoma 1993 (Vol. 52).
- 41. Mantua NJ, Hare SR. The Pacific decadal oscillation. J Oceanogr. 2002 Feb;58(1):35–44.
- 42. Di Lorenzo E, Schneider N, Cobb KM, Franks PJ, Chhak K, Miller AJ, et al. North Pacific Gyre Oscillation links ocean climate and ecosystem change. Geophysi Res Lett. 2008 Apr;35(8).
- 43. García-Reyes M, Largier JL, Sydeman WJ. Synoptic-scale upwelling indices and predictions of phyto- and zooplankton populations. Prog Oceanogr. 2014 Jan;120:177–188.
- 44. Akaike H. A new look at the statistical model identification. IEEE Trans Auto Cont. 1974 Dec;19(6):716–23.
- 45. Zimmerman RC, Kremer JN. 1984. Episodic nutrient supply to a kelp forest ecosystem in Southern California. J Mar Res 42:591–604.
- 46. Springer YP, Hays CG, Carr MH, Mackey MR. Toward ecosystem-based management of marine macroalgae—the bull kelp, Nereocystis luetkeana. Oceanogr Mar Biol Ann Rev. 2010;48:1–42.
- 47. Wilson KL, Tittensor DP, Worm B, Lotze HK. Incorporating climate change adaptation into marine protected area planning. Glob Change Biol. 2020 Jun;26(6):3251–67. pmid:32222010
- 48.
Merrill JE, Gillingham DM. Bull Kelp Cultivation Handbook. 1991. NCRI Publ. #NCRI-T-91-011.
- 49. Fredriksen S, Filbee-Dexter K, Norderhaug KM, Steen H, Bodvin T, Coleman MA, et al. 2020. Green gravel: a novel restoration tool to combat kelp forest decline. Scientific Reports 10:3983. pmid:32132550