A challenge modern-day ecologists and resource managers face is how to separate natural variations in populations from changes caused by human activities (e.g., climate change). Long-term monitoring programs provide valuable information to assist in this endeavor. This study details the initial findings of a long-term monitoring program initiated in 2004 to monitor changes in rocky intertidal communities within Redwood National and State Parks, located in northern California. Permanent plots were established at three sites using protocols developed by the Multi-Agency Rocky Intertidal Network, a consortium that monitors rocky intertidal communities along the western coast of North America. Replicate plots were established to monitor changes in abundance of key intertidal taxa, including mussels (Mytilus californianus), barnacles (Chthamalus dalli and Balanus glandula), red alga (Endocladia muricata), and rockweeds (Pelvetiopsis limitata and Fucus gardneri). Plots were sampled annually since they were established. Results from the first 15 years of this study indicate that all taxa exhibited substantial short-term (annual) variation, with barnacles and E. muricata exhibiting the most. For barnacles, such variations were correlated with measures of recruitment. Except for P. limitata, all other target taxa experienced at least one period of large-scale major change, where abundances decreased dramatically (> 50% of long term mean) and simultaneously in most plots. However, in almost all cases abundances recovered, resulting in no apparent long-term changes. For the few instances where long-term changes were detected, it is possible this result may be an artifact of the analytical methods used to assess them. The potential implications of this finding are discussed.
Introduction
It can be a challenge to assess the effects of human activities on ecosystems, a problem of increasing interest and concern, because populations vary naturally over time (e.g., Andrewartha and Birch 1954). Whether trying to assess the effects of oil spills (e.g., Peterson et al. 2003, Barron et al. 2020, Lewis et al. 2020), overfishing (e.g., McCauley et al. 2015, Gordon et al. 2018, Sumaila and Tai 2020), invasive species (e.g., Molnar et al. 2008; Simberloff 2013; Iacarella et al. 2019a, b), or human-induced climate change (e.g., Harley et al. 2006; Hawkins et al. 2008, 2016; Bindoff et al. 2019), it is important to know the “natural”, pre-disturbance condition of an ecosystem to assess potential effects of human activity. Long-term monitoring studies can provide the information necessary to address this issue (e.g., Magnuson 1990).
Lindemayer and Likens (2010) define a long-term monitoring study as a field study in which data are repeatedly and continuously collected for a period of at least 10 years. Numerous authors have discussed the importance of such studies (e.g., Magurran et al. 2010; Hughes et al. 2017; Miner et al. 2018, 2021). They can establish a baseline set of conditions against which future changes to the community can be measured. Such data are useful not only to interpret short-term ecological experiments (Magnuson and Bowser 1990), but to inform policy makers and resource managers tasked with mitigating the effects of natural and human-caused changes (e.g., Lindenmayer et al. 2012, Schuurman et al. 2020).
Recognizing the importance and benefits of long-term monitoring, in 1999 the National Park Service (NPS) established the Inventory and Monitoring Program, a nation-wide effort to provide park managers with scientific information about both the natural resources within the parks and how they are changing over time (Fancy et al. 2009). In 2004, the Klamath Inventory and Monitoring Network (KLMN) of the NPS, a unit of this monitoring network, initiated a long-term monitoring study to document the status and temporal dynamics of common intertidal invertebrate and algal taxa found within Redwood National and State Parks (RNSP). Established in 1994, this approximately 560 sq km park, located in northern California, is visited by an estimated 700,000 people annually. Although perhaps best known for its old-growth redwood forests, RNSP also includes 56 km of coastline.
Here, we present some initial findings from this ongoing long-term monitoring study. We utilized data from the first 15 years (2004 to 2018) of this study to describe the temporal changes in abundance of a select group of common intertidal organisms. We then used several different analyses to detect both short- and long-term patterns and trends in species abundance, and then discuss potential mechanisms driving these patterns. Finally, we offer insights into ways of improving future analyses of long-term change.
Methods
This ongoing study was initiated in 2004, when permanent plots were established at three sites within RNSP (Figure 1) to monitor changes in abundance of several common and ecologically important intertidal taxa. The site at Enderts Beach consists of several rocky benches separated by boulder-filled channels, while the site at Damnation Creek is an extensive channel-cut bench surrounded by a boulder and cobble beach. Both sites experience moderate to high wave exposure, and due to accessibility issues, receive few visitors. In comparison, the site at False Klamath Cove consists of several small benches within a large boulder field. It is somewhat more protected from waves and more accessible than the other two sites.
At each site, permanent plots were established to monitor abundance of six target species: 1) California mussel (Mytilus californianus Conrad), 2) red turfweed alga (Endocladia muricata Postels & Ruprecht), 3) common brown rockweed (Fucus gardneri Linnaeus), 4) dwarf rockweed (Pelvetiopsis limitata Gardner), 5) small acorn barnacle (Chthamalus dalli Pilsbury), and 6) common acorn barnacle (Balanus glandula Darwin). The plots were established using protocols developed by the Multi-Agency Rocky Intertidal Network (MARINe), a consortium of various governmental agencies, universities, and private organizations that monitors the intertidal community at sites along the western coast of North America; a detailed description of these protocols and their rationale can be found in Ammann and Raimondi (2008). Briefly, five permanent plots, each measuring 50 × 75 cm, were established per taxon; the corners of each plot were marked using stainless steel bolts and marine epoxy. The plots were installed where the abundance of targeted taxa was high (ideally > 80% cover); since each of these taxa is dominant at a different height on the shore, the monitoring plots were spread throughout the intertidal zone. At False Klamath Cove (FKC), plots were established for all target taxa, with one set of plots used to monitor the two barnacle species. However, at Enderts Beach (END) and Damnation Creek (DMN), only four plot types were established; F. gardneri was not common enough to establish plots at Enderts Beach, nor was P. limitata at Damnation Creek.
Since 2004, the abundance (% cover) of taxa within each plot have been assessed annually in May or June of each year (note: Enderts Beach was inaccessible in 2017). Initially, this was done in the lab by identifying and counting individuals under a grid of 100 uniformly spaced points overlaid onto photographs taken of the plots using imaging software (points spaced the digital equivalent of 5 cm apart widthwise, 7.5 cm apart lengthwise). However, starting in 2007, abundances were measured directly in the field by placing a quadrat, strung with line to create a grid of 100 uniformly spaced points (spaced approximately 5 × 7.5 cm apart as noted above) over the permanent plots. Again, the individuals under each point were identified and counted. Although there can be multiple taxa under each point due to layering, those on the bottom layer(s) were usually not visible in the early photographs scored in the lab. Therefore, for the sake of consistency, only the data for the top layer were considered for this study. The exception was the mussels in the mussel plots, where their presence was scored regardless of layer. Note that during the first two years of the study, barnacles were not identified to species, but they were from 2006 onward. Thus, percent barnacle cover represents the combined abundance of both C. dalli and B. glandula (except for recruitment, see below). This paper includes the data collected from 2004 to 2018.
In 2006, five 10 × 10 cm plots, one adjacent to each barnacle plot, were established to monitor barnacle recruitment. The corners of the plots were marked with marine epoxy, and then the plots were scraped clean. A 2-cm border around the plots was also scraped to minimize edge effects. The following year (2007), all barnacle recruits in the plots were identified and counted, after which the plots were again scraped clean. This was repeated every year until 2016.
Because mussel beds often consist of multiple rather than single layers of mussels, percent cover alone often does not provide a complete measure of mussel abundance. Thus, beginning in 2010, the size distribution of mussels in the M. californianus plots was also measured annually (but beginning in 2011 at END). The size measurements were taken by placing a small 10 × 10 cm quadrat in the four corners and center of each plot and measuring the length of the mussels under two corners of each quadrat. If there was more than one layer of mussels, as many additional individuals as possible were measured under each corner without physically altering the bed. Thus, a minimum of 10 individuals were measured per plot (2 corners × 5 quadrats), but sometimes more. Beginning in 2011 the thickness of the mussel bed was also determined by measuring the distance from the primary substrate (i.e., rock) to the top of the top-most mussel, with the aid of a metal skewer, in the center of each 10 × 10 cm quadrat. All measurements were taken to the nearest 1 cm.
Temporal Analyses
Several analyses were used to investigate temporal changes in abundance of the target taxa in their designated plots during the period from 2004 to 2018. Short-term annual changes were investigated by first using correlation analysis to compare abundances in sequential years: year(x) versus year(x+1). Second, periods of major change were identified by determining when the change in abundance between sequential surveys exceeded 50% of the long-term mean.
Two different statistical analyses were used to detect longer term changes in abundance. First, two-way randomized block analysis of variance (ANOVA) was used to detect non-linear changes in abundance over short (3 years), medium (5 years), or long (7 years) time scales. For each, the 15 annual surveys were divided into chronological groups containing the appropriate number of samples (e.g., for the short time scale there were 5 groups, each containing 3 surveys). Note that because there were 15 surveys, to maintain a balanced comparison (i.e., to compare equal time periods), the year 2011 was not used for the long time scale (7 years per group) analysis. Prior to analysis, the average abundance of the target taxa was first calculated for each replicate plot within each group; this mean was then log(x + 1) transformed. Any significant differences detected by the ANOVA were further investigated using Tukey's a posteriori test.
For the second analysis, we used the unweighted linear mixed model of Piepho and Ogutu (2002) to detect any linear trends in the abundance of the target taxa; Starcevich et al. (2018b) found this method was robust for trend detection. The analyses were done using the R package TrendNPS (Starcevich et al. 2018b), and all data were log(x + 1) transformed prior to analysis. When possible, the model included a site-level slope variable, a random variable that captured spatial variation at the plot level (Starcevich et al. 2018a). However, when the model was singular or failed to converge (likely due to over parameterization), this variable was omitted. Tests for significance were conducted by comparing the Wald statistic for the trend coefficient (i.e., slope) against the t-distribution, with degrees of freedom determined by the Satterthwaite method (Satterthwaite 1946).
Results
A total of 42 taxa were counted within the plots during the study, but the vast majority (76%) of these taxa had mean cover < 2%. The three most commonineachplot,whichincludedsomecombination of the five target taxa, plus open space and the two crustose algae (Ralfsia spp. and Petrocelis [the tetrasporophyte phase of the red algae Mastocarpus spp.]) accounted for, on average, 77 to 98% of all points scored (Table 1). Open space, which in most cases was second in abundance after the target taxa in each plot, tended to be strongly, and negatively, correlated with cover of the target taxa (Table 2).
The target taxa exhibited varying degrees of annual variation in abundance (Figure 2), with the barnacle and Endocladia plots appearing to exhibit the most (Figure 3). The pattern of variation often differed among the sites. For example, barnacle abundances at both Enderts Beach and Damnation Creek were lowest in 2010, while at False Klamath Cove, they were lowest at the start (2004) and end (2018) of the 15-year study period (Figure 2). E. muricata abundances started high and then dropped precipitously at all three sites, but only at Enderts Beach and False Klamath Cove did abundances quickly recover (Figure 2). Starting in 2014, percent cover again declined at Enderts Beach before beginning to recover (Figure 2). In the mussel (Mytilus) plots, mussel abundances were dramatically lower at Enderts Beach in 2015 and 2016, which coincided with a reduction in mussel bed thickness (Figure 4); there were no changes at the other two sites (Figure 2). Although less striking, site differences were also evident in the cover of F. gardneri and P. limitata in their respective plots (Figure 2). For the barnacle plots, the observed variation in cover was correlated with annual differences in barnacle recruitment (Figure 5).
While many of the short-term (annual) changes were small (< 10%), there were several instances where the differences were large (> 40%; Figure 3). Major changes occurred with some frequency at a small scale (i.e., occurring in < 2 of the 5 replicate plots), but also occasionally at a larger scale (> 3 of the plots) (Figure 6). Undoubtedly, the frequency of such changes contributed to the fact that, with a few exceptions, there was no temporal autocorrelation among samples (Table 3), this despite the fact that the taxa were sessile and abundances were measured in fixed, not random, plots.
In the barnacle plots, a major loss in abundance was often immediately (the next year) followed by a major increase, but this was less common for the other taxa. There were several instances (e.g., 2005 at Enderts Beach and False Klamath Cove in the barnacle plots, and 2009 at Enderts Beach in the Endocladia plots) when some replicate plots exhibited major losses, while the others were exhibiting major increases (Figure 6). There were also several instances (e.g., 2005, 2010, 2012, and 2015) when major losses were recorded both across plot types within a site and across sites (Figure 6).
TABLE 1.
The top three most abundant taxa or open space found in the top layer of each plot type at the three study sites. Numbers in parentheses are mean abundance (percent cover) over the course of the study. Non-coralline crust was Petrocelis plus Ralfsia spp. NP = plots of this type not present.
For longer term changes, the results of the ANOVA analysis found significant differences in abundance at the short (3 yrs/group) and medium (5 yrs/group) time scales for many of the taxa in their designated plots (Figure 7, Table 4). In almost all cases, these differences suggest abundances were either higher or lower in the middle, rather than the beginning or end, of the study (Figure 7); E. muricata and P. limitata at Enderts Beach were the exceptions, as both had abundances lower at the end of the study than the beginning (Figure 7). However, at the longest time scale (7 yrs/group), there were only two instances when abundances differed between the first 7 years versus last 7 years of the study (Figure 7, Table 4): during the second half of the study, barnacle cover at Damnation Creek was higher, and at Enderts Beach P. limitata cover was lower. For P. limitata, this trend is consistent with the findings of the short and medium time scales, but for the barnacles, it is not. Similarly, the abundance of E. muricata at Enderts Beach, which appeared to decrease over time at the short and medium time scales, was only marginally different when examined at the longer time scale (Table 4).
Perhaps not unexpectedly, the results of the linear trend analysis most closely matched those of the long time scale (7 years per group) ANOVA analysis in that, in almost all cases, there were no significant linear temporal trends in abundance (Table 5). However, there were two important differences. First, neither of the two cases where a significant linear trend was detected (E. muricata and M. californianus at Enderts Beach) were significantly different in the ANOVA analysis. And second, neither P. limitata at Enderts Beach, nor the barnacles at Damnation Creek, which exhibited significantly different mean abundances in the two halves of the study, exhibited a significant linear trend (Table 5). In three cases (P. limitata, E. muricata, and M. californianus at Enderts Beach), there was a large-scale, major decline in abundance, albeit in different years, during the second half of the study (Figure 6). That both M. californianus and P. limitata abundances had recovered by the end of the study (assessed by comparing abundance in 2018 versus mean abundance from 2004 to 2011 in a paired t-test [n = 5], t = 1.53, P > 0.05 for M. californianus, and t = 0.44, P > 0.05 for P. limitata), but E. muricata had not (abundance in 2018 versus mean abundance from 2004 to 2011, t = 3.34, P < 0.05) contributed to these different results. Recovery of M. californianus was primarily due to an influx of small mussels (Figure 4).
TABLE 2.
Correlation coefficients (Spearman's rho) comparing percent cover of open space to the percent cover of the target taxa for the three study sites. Data were log(x +1) transformed prior to analysis. a = P < 0.1, * = P < 0.05, ** = P < 0.01, *** = P < 0.001. Study sites are END = Enderts Beach, DMN = Damnation Creek, FKC = False Klamath Cove. NP = plots of this type not present.
Discussion
During the 15-year study period, all the target taxa showed some degree of temporal variation in abundance, including several instances of very large change. For the barnacles, there is evidence that variations in abundance were, in part, determined by variations in recruitment (see below). For the other taxa, the frequent incidents of both small- and large-scale major declines in abundance suggest that disturbance played a major role in determining their abundance. In general, disturbances are unexpected events that disrupt the structure of communities. However, here we use the term to represent events that kill or damage individuals (sensu Sousa 1984a). Although extreme temperatures can be agents of disturbance (e.g., Sebens and Lewis 1985, Suchanek 1985, Robles et al. 1995, Harley 2008, Harley and Paine 2009), in the rocky intertidal system along this coast, physical disturbances are usually related to high wave action and wave-borne objects striking the shore (e.g., Dayton 1971, 1973; Harger and Landenberger 1971; Paine 1979, 1988; Paine and Levin 1981; Sousa 1979, 1984b; Shanks and Wright 1986). Because such events create open space on the rock, they both renew an important and often limited resource (i.e., open space) and reset the successional sequence. Thus, disturbances play a key role in determining the structure in rocky intertidal systems (e.g., Levin and Paine 1974; Sousa 1979, 1984b; Paine and Levin 1981; Guichard et al. 2003).
Given the observed strong, negative correlations between the abundance of the target taxa and the amount of open space in the plots, sudden drops in abundance of the target taxa most likely represent disturbance events. As such, small-scale disturbances, as evidenced by losses in one or two replicate plots, appeared to happen with some frequency, particularly in the barnacle, Endocladia, and Pelvetiopsis plots. Why they were less common in the Fucus and Mytilus plots is unknown. For mussels, perhaps the protection offered by their shells (e.g., Burnett and Belk 2018, Crane et al. 2021), along with their ability to increase attachment strength in response to increased hydrodynamic forces (e.g., Witman and Suchanek 1984, Bell and Gosline 1997, Carrington 2002), offer them heightened levels of protection from some disturbances. For F. gardneri, because it is a canopy-forming species, the measures of its abundance reflect the size of the canopy, not the space it occupies on the rock. Because this canopy is quite flexible, it may be able to avoid damage from certain types of disturbances. Alternatively, if the density of F. gardneri is high enough, or individuals big enough, neighboring individuals drape over each other during low tide. Under these conditions, the loss of individuals due to disturbances may have more of an effect on the thickness of the canopy than on the measures of percent cover.
There were also instances when disturbances occurred on a larger scale, impacting multiple plots of a given taxon, or in several cases, plots of multiple taxa within and across sites. Again, the Fucus and Mytilus plots were the least affected. Although such instances represent periods of greater loss, when more individuals spread out over larger areas were eliminated from the population, given the sampling interval, it is unknown whether these disturbances happened all at once (e.g., one winter storm), or occurred over a longer period of time (e.g., one winter season). Since previous authors have had success linking offshore measures of oceanographic conditions to disturbances onshore (Denny 1995, Blanchette et al. 2002), this is clearly a topic deserving of further research.
Curiously, one large-scale event that apparently did not have much of an impact on the target taxa was the 2014 to 2016 marine heatwave (MHW). This MHW was an extended period of anomalously high water temperatures along much of the coastline of western North America (Gentemann et al. 2017); during this period we observed mean monthly water temperatures 1 to 3 °C above normal at the study sites. This MHW had major impacts on many nearshore and coastal communities (Cavole et al. 2016, Weitzman et al. 2021), including the rocky intertidal system, where it has been linked to the onset of sea star wasting syndrome (SSWS) and the subsequent coast-wide decline in sea star numbers (Eisenlord et al. 2016, Miner et al. 2018; however, see Menge et al. 2016). For example, numbers of adult ochre sea stars (Pisaster ochraceus Brandt) declined by more than 50% at the three study sites at the start of the MHW (Ammann et al. 2017). Thus, it is somewhat surprising that, except for widespread small-scale disturbances in 2015 and a major decline in mussels at Enderts Beach, the target taxa apparently were not greatly impacted by this MHW. Most marine taxa occupy the fullest range of latitudes allowed by their thermal tolerances (Sunday et al. 2012), and the target taxa at the three study sites are not close to their southern distributional limits (Table 6). This suggests that the above-normal temperatures experienced during the MHW were not high enough to cause widespread harm. This is supported by published temperatures tolerances (Table 6), which are much greater than the maximum water temperatures (< 17 °C, Ammann et al. 2017) observed at the study sites during this MHW. Since climate change is expected to cause an increase in the frequency and intensity of MHWs over the next century (Oliver et al. 2017, Frölicher et al. 2018, Collins et al. 2019, Darmaraki et al. 2019), how the target taxa respond to future MHWs is of interest. Similarly, since intertidal organisms must also deal with air temperatures during low tides, how they respond to climate-induced changes in air temperatures is also worthy of attention.
While the MHW may have played a role in the decline in mussel abundance at Enderts Beach, it could have also facilitated their subsequent recovery. While smaller gaps in mussel beds created by disturbances commonly fill in via encroachment of the surrounding mussels (Paine and Levin 1981, Sousa 1984b), the recovery of larger gaps depends more upon recruitment (Connell and Slatyer 1977, Paine and Levin 1981, Conway-Cranos 2012). The success of this process can depend on the actions of the ochre sea star (Pisaster ochraceus), a major predator on mussels (Paine 1966, 1974; Menge et al. 1994); previous studies have documented that this predator can eliminate patches of new mussel recruits (Sousa 1984b, Robles et al. 1995). However, Menge et al.(2016) found that the rate of predation on mussels was greatly reduced during the MHW, which they attributed to the loss of P. ochraceus due to SSWS. Since the loss of this predator has been linked to the post-SSWS downward (i.e., towards the water) expansion of mussel beds along the California coast, due in part to recruitment (Moritsch 2021), it is plausible that it also contributed to the large influx of small mussels in the mussel plots at Enderts Beach beginning in 2016.
Although the recovery of mussels at Enderts Beach demonstrates the importance of recruitment following disturbances, the significant correlation between recruitment in the barnacle clearings and cover in the plots suggests that, for barnacles, recruitment plays a key role in the temporal dynamics of their populations. This finding is consistent with previous studies that concluded intertidal populations along the stretch of coastline encompassing the three study sites are recruitment limited (Connolly and Roughgarden 1998, Connolly et al. 2001, Broitman et al. 2008). For such populations, not enough recruitment occurs for open space to become limited. Consequently, post-settlement mortality is density independent, and variation in adult abundance reflects variation in recruitment (Caley et al. 1996). Additionally, although B. glandula can live up to 10 years, and C. dalli for up to 3 years (Newman and Abbott 1980), the significant correlation between recruitment and cover also suggests that the populations in the plots were comprised largely of a single generation. This further implies not only that some processes, such as winter storms or predation, remove most of the population each year, but that the long-term persistence of these populations is reliant on continual recruitment.
Long-term Changes
Despite evidence of disturbance, in most cases there were no significant long-term changes in abundance, most likely because following each disturbance was a period of “recovery”; this was true even for those exceptions that exhibited long-term declines (see below). In the barnacle plots, such periods often occurred immediately following the decline, but for the other taxa, it often took a little longer. Previous studies have found that not only are barnacles often the first to colonize following a disturbance, but that their presence can facilitate the appearance of later successional stages (Dayton 1971, Farrell et al. 1991, Van Tamelen and Stekoll 1996, Berlow 1997, Kim 1997, Menge et al. 2011). Thus, the apparent differences in the rate of recovery could reflect successional dynamics. Since reproductive seasons vary among the taxa, the timing (i.e., season) of the disturbance event could also play a role (Kim and DeWreede 1996). Alternatively, since algae can expand via vegetative growth (Sousa 1980, 1984b; Sousa et al. 1981; Conway-Cranos 2012), it could simply reflect a different process of recovery (i.e., outward expansion of existing individuals versus colonization by planktonic larvae), one that operates at a different time scale.
TABLE 3.
Correlation coefficients (Spearman's rho) comparing taxa abundance at time(t) versus time(t + 1). Data were log(x + 1) transformed prior to analysis. a = P < 0.1, * = P < 0.05, ** = P < 0.01. Study sites are END = Enderts Beach, DMN = Damnation Creek, FKC = False Klamath Cove. NP = plots of this type not present.
For those cases that demonstrated long-term declines (M. californianus, E. muricata, and P. limitata, all at Enderts Beach), the timing and frequency of the disturbances played an important role in this outcome. For example, in the case of M. californianus, there was only one instance of a major disturbance, and it occurred near the end of the study. Although abundances had recovered by the last census (2018), it was depressed for long enough to reduce the mean abundance of the second half of the study relative to the first. As for P. limitata and E. muricata, both experienced two major disturbances, one in each half of the study. P. limitata recovered from both, but because it took longer to recover after the second disturbance, abundances were lower for a longer period during the second half of the study than during the first. Although its abundance had not yet fully recovered by the end of the study, the same was basically true for E. muricata. Thus, while in all three cases the statistical tests suggest overall abundances were lower in the second half of the study, which is suggestive of long-term change, arguably this result may have less to do with the actual state of the populations and more to do with the methodology of the analyses used to detect these changes (discussed further in the next section).
There are two additional points worth noting. First, the highest abundances of E. muricata at Enderts Beach over the course of the study were observed at the start. Although this could be an artifact of the criteria used to establish the plots (i.e., they were installed where abundances were high), it could also be an example of a shift in baseline abundance (that is, the level around which abundances vary over time). In this case, after both postulated major disturbance events, abundances of E. muricata increased, but after the second, they did not return to the same level. Whether, given more time, abundances will eventually do so, or whether they will remain at a lower level, can only be determined by continued monitoring.
Second, if this change truly represents a shift in the baseline, it only took place at one site, Enderts Beach. If it were due to some large-scale phenomenon, like climate change, similar changes would be expected at the other sites. Since this was not the case, the shift likely reflects a change in local, not global, conditions.
Implications
To understand how human activities affect natural communities, there has been renewed interest in observational ecology in general (Sagarin and Pauchard 2010) and a recognition of the importance of long-term monitoring studies (e.g., Lovett et al. 2007; Magurran et al. 2010; Lindenmayer et al 2012; Hughes et al 2017; Miner et al. 2018, 2021). By collecting data both continuously (i.e., annually in the current study) and over an extended period of time (Lindenmayer and Likens 2010), such studies provide a more complete picture of the temporal dynamics of a population. Both aspects of such studies are important. Because data are collected continuously, it is possible to determine whether an observed change is more extreme than ‘normal’, whether it is temporary or permanent, or whether or how often similar changes have occurred before. Such knowledge provides important context, making it possible to identify those changes that are unexpected or extreme. Also when data are collected over long periods of time, it is possible to establish a baseline measure of abundance, a value against which any observed change can be compared. Because all populations exhibit some degree of short-term variation, such information is essential for detecting long-term changes, like those predicted for climate change (e.g., Harley et al. 2006; Hawkins et al. 2008, 2016; Bindoff et al. 2019).
TABLE 4.
F-statistics of 2-way analysis of variances comparing mean abundance among groups defined using three different intervals: 3, 5, or 7 years per group for study sites. Results of Tukey's multiple comparison tests are presented in Figure 5. df = degrees of freedom. a = P < 0.1, * = P < 0.05, ** = P < 0.01, *** = P < 0.001. Study sites are END = Enderts Beach, DMN = Damnation Creek, FKC = False Klamath Cove.
TABLE 5.
Annual percent change estimates of linear mixed models (Piepho and Ogutu 2002) used to detect linear trends of target taxa at study sites over time. Data were log(x+1) transformed prior to analysis, and results (i.e., estimates of slopes) were back-transformed. CI = confidence interval; * = slope significantly different than 0 (P < 0.1); ϕ = the random site-level slope variable was retained in the model. Study sites are END = Enderts Beach, DMN = Damnation Creek, FKC = False Klamath Cove.
Results from the first 15 years of this ongoing study reinforce the importance of both aspects of long-term monitoring studies. Consider the mussel plots at Enderts Beach, one of the two instances (the other being E. muricata at Enderts Beach) that exhibited a significant linear change in abundance over the course of the study. Because the data were collected continuously over a long period of time, we know: 1) this population experienced only one major disturbance event, the only mussel population to do so, 2) this event occurred towards the end of the 15-year time period, and 3) abundance returned to pre-disturbance levels relatively quickly. Thus, while the mixed linear model suggested a decline in abundance, a closer look at the temporal dynamics of the population, possible because this was a long-term study, offered a different interpretation. Specifically, except for the short period encompassing the disturbance/recovery event (2015 to 2018), abundances were relatively stable; abundances were no different in the short and medium time scales ANOVA analyses. Again, it was the timing of the one major disturbance event that led to the significant result. This suggests that care should be taken when using a linear model to analyze data that are either not necessarily linear or that have longer cyclical patterns (e.g., Denny and Paine 1998, Burnaford et al. 2014), which when analyzed over shorter time periods, could give alternate interpretations.
TABLE 6.
Southern distributional limits and thermal tolerances of the target taxa. NI = no information.
This study illustrates the importance of defining an appropriate baseline. Since baselines represent the condition against which future population change will be assessed, they should typify the average condition of the population, not an arbitrary point in time (e.g., the starting point of a project). An appropriate baseline could depend upon the question being addressed. For example, determining the effects of an oil spill could require a shorter temporal baseline than one to detect the effects of long-term climate change. We recommend that the minimum length of time needed to define an appropriate baseline, especially in a disturbance-prone community, include at least one major disturbance and recovery event. This may be taxa or site dependent; for example, based on this study, the baseline for M. californianus may need to be longer than one for E. muricata. Defining an appropriate baseline, especially one to detect longer term changes, should also consider periodic or cyclical events known to affect the intertidal communities; these include the El Nino-Southern Oscillation, the North Pacific Gyre Oscillation, and the Pacific Decadal Oscillation (e.g., Menge et al. 2009, 2011; Spiecker and Menge 2022), as well as the nodal tidal cycle (e.g., Denny and Paine 1998, Burnaford et al. 2014). While obtaining an appropriate baseline may encompass a length of time that can be both logistically and financially challenging, to assess the effects of long-term climate change, it may be necessary. Thus, there is a need for committed long-term monitoring studies, like the National Park Service Inventory and Monitoring Program.
Finally, the results from the current study support the idea that, when linked with other biological or physical measurements, long-term monitoring studies can help identify potentially important processes and mechanisms that determine population size. For example, combining the long-term abundance data with the data from the recruitment study revealed that the dynamics of the barnacle populations at the study sites were very dependent upon recruitment. Given the important role barnacles play in this system, from important prey item (e.g., Connell 1970, Spight 1982, Palmer 1984), to competitor (e.g., Choat 1977, Van Tamelen 1987), to facilitator (Dayton 1971, Farrell et al. 1991, Van Tamelen and Stekoll 1996, Berlow 1997, Kim 1997, Menge et al. 2011), fluctuations in barnacle abundance can have important consequences for the community. Similarly, by combining long-term abundance data with data on sea surface temperatures, Miner et al. (2021) were able to show that ocean temperatures play an important role in the stability of intertidal communities along the western coast of North America. Examples such as these demonstrate that, when combined with other relevant biological and physical data, long-term monitoring studies can be a powerful tool in helping further our understanding of what drives the temporal dynamics of natural populations and communities.
In conclusion, many resource managers, like those in the National Park Service, face the difficult task of both preserving natural resources and encouraging public use (e.g., recreation, fishing, etc.) of those same resources. In an environment of ever-increasing human influence, including climate change, there is a need for a more nuanced understanding of “baseline” when making management decisions (Schuurman et al. 2020). The types of information collected by long-term monitoring can help lead to better, more informed decisions about how to best achieve these often conflicting goals.
Acknowledgments
Funding for this effort was provided by the National Park Service through the Klamath Inventory and Monitoring Network, managed by Alice Chung-Maccoubrey. The principal investigator for this monitoring project is Dr. Peter Raimondi. We also thank the many non-author researchers, technicians, students, and volunteers who helped collect and manage data, including Allison Snyder, David Anderson, Heather Brown, Kyle Max, and Terri Williams from NPS, and Maya George, Melissa Douglas, Laura Anderson, Nathaniel Fletcher, Rani Gaddam, and Melissa Miner from UCSC. The manuscript was improved thanks to input from three reviewers. Authors have no conflicts of interest to declare.