ORIGINAL RESEARCH article

Front. Mar. Sci., 20 October 2023

Sec. Marine Fisheries, Aquaculture and Living Resources

Volume 10 - 2023 | https://doi.org/10.3389/fmars.2023.1270603

Southern California marine protected areas promote bolder fish populations

  • Division of Natural Sciences, Pepperdine University, Malibu, CA, United States

Abstract

Marine protected areas (MPAs) aim to protect habitats and ecosystems to promote the diversity and health of marine populations. To evaluate the health of fish populations within and outside of MPAs in Southern California, we used flight initiation distance (FID). FID is the distance at which an individual will flee from a perceived predator and is a direct measure of boldness. Lower FIDs are indicative of bolder populations. Lower FID values indicate that fish will have the opportunity to dedicate less energy to fleeing from predators which could then be used for increased foraging. Data was collected from eight locations along the Southern California coast, four in protected areas and four in non-protected areas. FIDs were measured with a sonar-based range finder. All fish species that recreational fishermen could catch were targeted, with a focus on six core species spearfishermen catch commonly. Results showed that all fish sampled in MPAs have significantly lower FIDs than those in non-protected areas. An ANOVA showed no significant difference between different locations of the same protected status. The six core species had lower average FIDs in protected areas, three species had significant differences between protected and non-protected populations, and there was no significant difference in FIDs between species. These results show that Southern California MPAs are promoting healthier fish communities with respect to boldness behavior. This data and methods can be implemented in future projects to expand the range of metrics measured to inform management tools for California’s MPAs and MPAs at large.

1 Introduction

Long-term overfishing has led to the collapse of fisheries and the degradation of marine ecosystems around the world (). To protect ocean habitats and the communities they support, many countries have established marine protected areas (MPAs). MPAs limit or prevent take of marine life within their boundaries in order to allow populations to recover to healthy levels. Fish populations in California have followed the same global trends, as catch normalized by fishing pressure peaked before World War II and has declined since then (). Healthy fish populations are important for robust coastal ecosystems and the California economy. In California alone commercial and recreational fishing generated over 24 billion dollars and employed over 140,000 people in 2016 (). Fish populations decline led to California’s passage of the Marine Life Protection Act (MLPA) in 1999 which sought to create a network of MPAs to increase coherence and effectiveness in protecting the state’s marine ecosystems, and marine natural heritage, as well as to improve recreational, educational and study opportunities provided by marine ecosystems (; Saarman and Carr, 2013; ). Over the next 13 years, MPAs were created and implemented along the California coast by different agencies and with varying restrictions. In 2012 the California MPA Network (124 MPAs) was completed.

With 2022 marking the tenth year of the full implementation of the California MPA Network, there has been a decadal management review evaluating the effectiveness of the network in meeting the MLPA goals (Wertz et al., 2022; The First 10 Years: Measuring the Success of California’s Underwater Parks, 2023). Fish within the boundaries of MPAs are larger and more abundant than those outside of protected areas (). Most protected areas also had greater biomass of target species than reference sites outside of MPAs (). MPAs in Southern California have had the largest positive impact on populations within their boundaries, with greater biomass increases seen in these MPAs than in ones farther north (Sullivan-Stack et al., 2022). Reasons for this are that the Southern California MPAs were established earlier, meaning there has been a longer period to see change (Sullivan-Stack et al., 2022), and that more temperate ecosystems are more responsive to change than cooler areas (Tanaka et al., 2021). There has been a lot of work done looking at the abundance and size of marine life in protected areas, but the effects of protection status on behavior have been limited to one study focused only on barred sand bass (Paralabrax nebulifer) around Catalina Island ().

In order to evaluate the effectiveness of MPAs along the Southern California coast in promoting healthier communities of fish targeted by fishermen, we used flight initiation distance (FID). Flight initiation distance is a commonly used metric in behavioral ecology that measures the balance of current and future fitness of potential prey in response to predation events (Ydenberg and Dill, 1986). FID represents the optimization of current and future fitness; it is a metric of how an animal can best survive in the present and be able to produce and care for offspring in the future (). It has been extensively used to evaluate the effects of various anthropogenic impacts on species from reptiles to mammals (; ; McGowan et al., 2014; ). These studies have shown that increased urbanization leads to acclimatization by birds, reptiles, and mammals, showing shorter FIDs in locations where there is more human presence and that these changes can vary by species, not just phyla (; McGowan et al., 2014; ). It has also been shown that birds will modify their behavior and have longer FIDs in response to consistent threats in an area ().

An animal’s response to threats affects many different levels of animal behavior, from social structures to foraging habits (; ). These fear responses vary with the landscape, varying levels of predation and the availability of hiding places change behavioral responses, leading to the concept of a landscape of fear (), where prey’s perception of predation risk is different in different areas. These fear responses come largely from the threat of predation, so predator presence is especially important in constructing the landscape of fear. Importantly, even non-consumptive interactions between predator and prey can drastically alter prey behavior as prey will exhibit strong flight responses and dedicate less energy to foraging even with just a threat of predation (). Even without actual predation, the effect of potential predator interactions can affect community composition across trophic levels (). Not only is decreasing population sizes a problematic result of increased hunting, but changing the behavior of prey by manipulating this landscape of fear can have similarly negative effects. It has also been found that in some cases human interactions can elicit greater behavioral responses than the presence and absence of natural predators. found that the effect size of human disturbance was greater than changes from habitat type or natural predators for elk in the Rocky Mountains. This means that the effects of human predation could impact species’ behavior and by extension community interactions and composition much more than would occur from predation only by natural predators. Protected areas in both terrestrial and aquatic environments aim to limit the effects of human predation to reduce the significant changes that can result and to allow for healthy communities across trophic levels to thrive ().

Fishing, especially spearfishing, has been shown to increase fear responses of fish (; ). However, many different factors can influence fish flight responses. One of the ways that fish try to combat threats of predation is through schooling (), but studies around the effect of schooling report mixed results, with some showing that schooling fish allow a closer approach (Stankowich and Blumstein, 2005), while others say there is no effect (; Samia et al., 2019). Fish FID responses have also varied with the presence of hunting gear such as a speargun (; Tran et al., 2016). The diver’s starting distance from the fish (i.e., distance between the diver and the fish when the fish was first spotted), (Stamoulis et al., 2019), habitat complexity (Nunes et al., 2015), and fish body length () have also been shown to affect how different species respond to potential human predation. FID has been used to measure the effects of fishing pressures on fish behavior, especially comparing protected and non-protected areas (). In studies conducted along the coasts of Barbados, Brazil, and Rapa Nui fish in protected areas had lower FIDs than those in areas where fishing was allowed (; ; Petit et al., 2022). FIDs have begun to be used as a metric of the effectiveness of MPAs as they show behavioral characteristics associated with higher fitness. They can be beneficial in several different ways. First, they can be used as a metric of the effectiveness of protection. Part of the purpose of protection of fishing areas can be to change fish behavior to allow fish to be easier to catch in non-protected areas (). FID can be used to see if fish are becoming less wary, and as a result, easier to catch in fished areas around MPAs.

Additionally, changes in behavior can indicate the success of protection in promoting healthier communities of fished species. An increase in predator presence can prompt an increase in prey fish wariness (Milinski, 1986; ). As prey wariness increases, these fish dedicate less energy to foraging (), which can reduce growth rates (). As one of the metrics of the success of an MPA is the increase of fish size in its bounds (Turnbull et al., 2021), MPAs with less wary fish due to decreased fishing pressure should end up having larger fish as they can dedicate more energy to foraging. It has been shown that fish with traits correlated to increased foraging behavior, such as greater boldness and activity rates, are more likely to be caught by fishermen (). It has also been observed that populations outside of protected areas begin to evolutionarily self-select traits that lead to slower growth rates, while those inside MPAs exhibit behaviors leading to faster growth (Sørdalen et al., 2022). Larger fish, in turn, lead to more and greater viability of offspring (Qvarnström and Price, 2001; ; Pew Charitable Trusts, 2016), allowing populations to recover to pre-fishing levels. When bolder individuals are allowed to grow without the threat of human predation, they increase the health of the population through greater size and reproduction rates.

Another way that FID has been used is to measure the efficacy of MPAs in preventing marine take (; Nunes et al., 2018). Once a baseline FID is established for the species in an area, FID, along with other quantitative data, can be used over time to evaluate whether the restrictions are being followed. If in some protected areas FIDs are more consistent with fish in non-protected areas, that can indicate that the regulations in place are not being followed. As FIDs can be influenced by a host of factors, from habitat complexity (Nunes et al., 2015) to fish size () it is not enough to show the effectiveness of protection on its own. However, when used with other metrics of population health, such as biomass and biodiversity, it can inform decisions on the management and enforcement of regulations in protected areas.

Most studies using fish FID were conducted around the Pacific Islands, with some along the eastern coast of North and South America (Nunes et al., 2018). No studies have used FID of fish along the west coast of North or South America (Nunes et al., 2018). In California, unreliable visibility and surf conditions are likely reasons for this. Higher concentrations of nutrients in the water can lead to increased algae presence (; Smith, 2003; ), which contributes to lower visibility in the water. To our knowledge, this study is the first to use FID as a metric to evaluate the effectiveness of MPAs in fostering healthy fish populations in such a setting. FID has been used to evaluate the effectiveness of MPAs in Brazil and other countries (; ; Petit et al., 2022). It represents a tool to increase the scope of management strategies for protected areas in California. If well implemented, FID can be paired with more commonly collected data, such as fish biomass, size, counts, or landings to evaluate the effectiveness of MPAs. Using FID will allow managers to evaluate not only the physical changes in abundance and sizes but also the changes in fish behavior.

This study aims to evaluate the effectiveness of MPAs in Southern California in promoting healthy fish communities and to test the efficacy of FID as a metric of fish community health in Southern California waters. We hypothesized that fish populations would be healthier within the boundaries of MPAs, resulting in smaller FIDs for fish inside protected areas. To test this hypothesis FID was measured for fish species that are usually caught by spearfishermen within and outside of MPAs in Southern California.

2 Methods

2.1 Locations

The data was collected at eight different locations, four protected areas and four non-protected areas (Figure 1). The protected sites are all in State Marine Reserves (SMRs) or State Marine Conservation Areas (SMCAs). SMRs prevent the take of any living, geological, or cultural marine resource within their boundaries, while SMCAs may allow restricted take within their boundaries. The two SMCAs we sampled are Point Dume SMCA and Laguna Beach SMCA, which allow for the take of pelagic fish and incidental take of organisms during dredging, respectively. Big Dume Beach and El Matador Beach are at opposite ends of a continuous stretch of protected coast spanning two MPAs. The Point Dume SMR protects 19.50 km2 of coastal water around Point Dume. Point Dume SMCA extends over six kilometers west from the western boundary of the SMR and encloses 41.23 km2. The Point Dume SMCA allows the take of only pelagic fish, none of which were sampled in our study. Big Dume Beach is in Point Dume SMR and El Matador beach is in Point Dume SMCA. Woods Cove and Tablerock Beach are in a section of Orange County Coast protected by four separate protected areas. Woods Cove is in the Laguna Beach SMR, which covers 7.09 kilometers of coast and encloses 17.40 km2 of ocean. As an SMR, there is no take allowed within the boundaries. Tablerock is in the Laguna Beach SMCA and covers a stretch of ocean on the southern side of the SMR. It encloses 8.00 km2, and only incidental catch related to coastal management is allowed. All non-protected sites are along the Malibu coastline. Latigo is the farthest west at Latigo Point. Puerco is 4.5 km east of Latigo alongshore and off of Malibu Bluffs Park along Malibu Road. Surfrider is 2 km east of Puerco near the Malibu Lagoon, and Big Rock is 7 km east of Surfrider at the end of Big Rock Road. None of these areas have any protected status and spearfishing is common at all of the non-protected sites.

Figure 1

The reef fish we targeted have high site fidelity. These fish only leave the natural reefs they are found on originally to populate artificial reefs introduced in the area (). Since all of our sites were natural sites, we do not expect the movement of these fish between sites. In addition to this, many of our core species exhibit high site fidelity, and any movement that is seen between sites occurs at much shorter distances (i.e., 1 km) than those between any of our sites (; Valle, 1989; ) (shortest proximity between our sites is 2 km). All sampling was conducted in October and November of 2021 and between May and July of 2022 (Table 1).

Table 1

DateSiteNumber of ApproachesStatus
2021-10-29Big Dume Beach13Protected
2021-10-29Surfrider10Non-protected
2021-11-02Big Dume Beach12Protected
2021-11-12Surfrider12Non-protected
2021-11-23Woods Cove18Protected
2021-11-23Big Rock11Non-protected
2021-11-24Tablerock Beach10Protected
2022-05-13Big Dume Beach12Protected
2022-05-25Latigo9Non-protected
2022-06-20Surfrider23Non-protected
2022-06-21Big Dume Beach21Protected
2022-06-22Puerco19Non-protected
2022-07-08El Matador15Protected

Date (Year-Month-Day) and approach number of each sampling session for each site.

2.2 Fish species

Any fish encountered larger than 10 cm, indicating adult specimens, were targeted for FID measurements. The dataset includes 23 fish species (Supplemental Table 1). The only species not targeted was Hypsypops rubicundus, commonly known as the Garibaldi, because it is protected across the state and unlike the other fish observed is highly territorial. All of the species were included for all statistical tests expect the tests looking at the individual six core species (section 2.4). The following six core species were tested individually for differences between protected and non-protected populations (section 2.4): black perch (Embiotoca jacksoni), kelp bass (Paralabrax clathratus), opaleye (Girella nigricans), sargo (Anisotremus davidsonii), silver surfperch (Hyperprosopon ellipticum), and zebraperch (Kyphosus azureus). These core species were subsetted to have at least 10 total approaches with a minimum of three approaches in both protected and non-protected areas. Topsmelt (Atherinops affinis) also meets these requirements, but spearfishermen do not commonly catch it, so it was not included in our core species. The six core species are all taken by spear and line fishermen (Porzio, 2014; Stenstrom, 2018); therefore, they are good indicators of the effectiveness of fish population protection in MPAs. All of the core species are common along Southern California shores (Miller and Lea, 1972).

2.3 Sampling procedures

The data is a compilation of individual fish approaches (n=185) collected by snorkeling and freediving. Only one FID data point was collected per dive. Dive depths varied between two and eight meters. The diver swam along the bottom until they saw a target fish. The diver would swim straight at the fish at a constant pace of about one m/sec. FID was taken when the fish began to swim away from the diver faster than it was being approached, increasing the distance between it and the diver. FID measures the closest distance that the diver was able to get to the fish. Distance between the diver and the fish was taken with a RICANK Depth Finder, which is usually used to measure the distance to an object (i.e., seafloor or a fish) when ice fishing. The depth finder reports distances from 0.8 to 90 meters down to the tenth of a meter accuracy. In some cases, the distance directly to the fish could not be obtained, so the distance to the nearest rocks behind the fish was measured instead. Sometimes, the fish was too close to the diver for the range finder to get a reading, so a visual estimate was taken. This only occurred on some distances under 0.8 meters. All sampling sites were rocky bottom reefs with kelp beds; hence the distance to the nearest hard substrate behind the fish was comparable to the distance to the fish itself. In order to prevent sampling of the same fish multiple times, divers would swim unidirectionally east or west along the shore at least 10 meters away from the previous sampling location before sampling another fish of the same species. FID, species name or description, and location were recorded by the diver for all approaches. For 99 of the 185 total data points (i.e., 54%), schooling behavior was also recorded and defined as four or more fish of the same species swimming as a group and recorded by the diver. The FID of schooling fish was determined for the closest fish in the school that the diver was able to approach. Species were identified on the dive by the diver. Fish that could not be identified in the field were identified using A Field Guide to Coastal Fishes From Alaska to California () right after the dive based on the fish’s morphological characteristics. Fish that could not be identified in the field or with the field guide afterward were not included. There were only three individuals that could not be identified out of 188 total approaches (i.e., 1.5%), leaving the final data set at 185 individuals.

2.4 Statistical analysis

Our first two-way ANOVA was conducted with all FID data from all approaches (n=185) in order to see how a community responds to protected status (Table 2; Figure 2A). We conducted a two-way ANOVA using species and protected status as factors in order to account for differences in community composition in protected and non-protected areas.

Table 2

Test TypeData UsedGroups ComparedObjective
Two-way ANOVARaw FID of all approachesSpecies by protected statusSee if there are trends between protected and non-protected areas
Raw FIDs from protected sitesSpecies and site of the four protected sites: Big Dume Beach, El Matador, Tablerock Beach, and Woods CoveEvaluate whether communities of fished species had different FIDs in different protected areas
Raw FIDs from non-protected sitesSpecies and site of the four non-protected sites: Big Rock, Latigo, Puerco, and SurfriderEvaluate whether communities of fished species had different FIDs in different non-protected areas
Three-way ANOVARaw FID of all approaches for which schooling data was recordedSpecies by protected status by schooling behaviorEvaluate whether schooling behavior changes FIDs
Student’s T-testRaw FIDs of species testedEmbiotoca jacksoni: MPA vs. non-MPACompare individual species inside and outside MPA
Anisotremus davidsonii: MPA vs. non-MPACompare individual species inside and outside MPA
Hyperprosopon ellipticum: MPA vs. non-MPACompare individual species inside and outside MPA
Kyphosus azureus: MPA vs. non-MPACompare individual species inside and outside MPA
Wilcoxon Rank-sum testRaw FID of species testedParalabrax clathratus: MPA vs. non-MPACompare individual species inside and outside MPA
Girella nigricans: MPA vs. non-MPACompare individual species inside and outside MPA

All of the statistical tests conducted with the data used per analysis, the groups compared in each analysis, and the objective of the statistical test.

Figure 2

Next, we looked at the subset of the data (n=99) which includes schooling data to test whether differences in FID could be associated with schooling behavior. We used a three-way ANOVA with species, protected status, and schooling behavior as factors in order to determine whether coastal fish in Southern California vary in FID with different schooling behaviors (Table 2; Figure 2B).

We then used two two-way ANOVAs to evaluate whether fish at different sites with the same protected status (e.g, Big Dume Beach and El Matador - both protected areas) had different FIDs. We conducted two two-way ANOVAs on populations of fish within and outside of protected areas, respectively. The first two-way ANOVA used the variables of species and site on a subset of the overall data which only included approaches in MPAs (n=101, Table 2; Figure 3). The second two-way ANOVA used the same variables on a subset of the overall dataset which only included approaches in non-protected areas (n=84, Table 2; Figure 3). These tests evaluated whether populations in different sites of the same protected status had different behaviors.

Figure 3

We performed two Wilcoxon rank-sum tests and four Student’s T-tests on a subset of the data set that only included six of the species sampled. We chose these six species based on three criteria: (1) these six species had at least three approaches in each protected status, (2) these species had at least ten approaches total, and (3) all six species are commonly taken by spearfishermen in Southern California. The six species included in this test were black perch (Embiotoca jacksoni), kelp bass (Paralabrax clathratus), opaleye (Girella nigricans), sargo (Anisotremus davidsonii), silver surfperch (Hyperprosopon ellipticum), and zebra perch (Kyphosus azureus). With this group, we performed individual Student’s T-tests or Wilcoxon rank-sum tests with each species, based on whether the data passed Normality assumptions. These tests compared the FIDs of individuals approached inside MPAs to those approached outside of MPAs. For Embiotoca jacksoni, Anisotremus davidsonii, Hyperprosopon ellipticum, and Kyphosus azureus we performed Student’s T-tests because their populations both inside and outside MPAs were Normally distributed. We performed Wilcoxon rank-sum tests for Paralabrax clathratus and Girella nigricans because the data for these species was not Normally distributed. For all analyses, p-values less than 0.05 were considered statistically significant. All statistical tests described in this section are detailed in Table 2. All data analysis was done in Rstudio (Version 2023.06.2 561).

3 Results

There were 185 total approaches and FIDs recorded, 101 in protected areas and 84 in non-protected areas. There were at least nine approaches for each site, with Big Dume having the most approaches of a protected site (n=60) and Surfrider having the most data of a non-protected area (n=45) (Table 3). The two-way ANOVA conducted by species and protected status found that there was a significant difference with regards to species and protected status (Two-way ANOVA: species: df=23, F=1.88, p=0.0141 protected: df=1, F=47.5, p=1.39x10-10, Figure 2A). The three-way ANOVA looking at the effect of schooling found no relationship between schooling behavior and FID, with or without the effect of protected status considered (Three-way ANOVA: species: df=12, F=1.81, p=0.0642 schooling: df=1, F=0.399, p=0.530 protected: df=1, F=29.1, p=9.03x10-7, Figure 2B). The two two-way ANOVAs evaluating differences in FID between sites of the same protected status found that there were no differences within the four protected sites or the four non-protected sites (Two-way ANOVA in MPA sites: species: df=20, F=1.50, p=0.0821 site: df=3, F=0.571, p=0.636; Two-way ANOVA in non-MPA sites: species: df=12, F=1.55, p=0.133 site: df=3, F=0.742, p=0.531, Figure 3).

Table 3

SiteMPA StatusAverage FID ± SD (m)Sample Size
Big RockNon-protected2.1 ± 0.9311
LatigoNon-protected1.7 ± 0.659
PuercoNon-protected1.9 ± 0.9219
SurfriderNon-protected1.9 ± 1.245
Big DumeProtected1.1 ± 0.7258
El MatadorProtected0.97 ± 0.5315
TablerockProtected0.45 ± 0.3410
Woods CoveProtected0.96 ± 0.6818

Average flight initiation distance (FID) ± standard deviation (SD) in meters (m) for the eight sampling sites.

For the six core species, Embiotoca jacksoni, Anisotremus davidsonii, and Hyperprosopon ellipticum had under ten approaches in protected areas, with seven, seven, and eight approaches, respectively (Table 4). Only Embiotoca jacksoni and Anisotremus davidsonii had under ten approaches outside of MPAs, with six and three respectively, while Hyperprosopon ellipticum had ten approaches in non-protected areas (Table 4). Paralabrax clathratus had the most approaches, with 21 in MPAs and 19 in non-protected areas (Table 4). Girella nigricans had 12 in protected areas and 21 outside of protected areas, while Kyphosus azureus had 13 in MPAs and ten in non-protected areas (Table 4). All species had lower average FIDs in protected areas. When tests were performed between individuals in protected and non-protected areas for each species, there were only significant differences between Paralabrax clathratus (Pc), Girella nigricans (Gn), and Kyphosus azureus (Ka) (Students T-test: Ka t=3.35, df=21, p=0.0031; Wilcoxon Rank-Sum Test: Gn W=204, p=0.0035 Pc W=357, p=2.04x10-5, Figure 4) (Table 5). There were no significant differences between populations in protected and non-protected areas for Embiotoca jacksoni (Ej), Anisotremus davidsonii (Ad), and Hyperprosopon ellipticum (He) (Student’s T-test: Ej t=0.941, df=11, p=0.367; Ad t=0.724, df=8, p=0.490; He t=1.04, df=17, p=0.313, Figure 4) (Table 5).

Table 4

SpeciesProtected FID (m)Sample sizeNon-protected FID (m)Sample size
Embiotoca jacksoni1.2 ± 0.5671.5 ± 0.456
Paralabrax clathratus0.72 ± 0.69212.3 ± 1.219
Girella nigricans1.2 ± 0.85122.3 ± 1.121
Anisotremus davidsonii1.5 ± 0.7471.9 ± 0.853
Hyperprosopon ellipcticum1.1 ± 0.6981.3 ± 0.7710
Kyphosus azureus0.98 ± 0.56132.1 ± 0.9910

Average fid (m) with standard deviations for the six core species in both protected and non-protected areas.

Table 5

TestGroups ComparedSummary Statistics
Student’s T-testtdfp
Anisotremus davidsonii protected vs non-protected0.72480.490
Embiotoca jacksoni protected vs non-protected0.941110.367
Kyphosus azureus protected vs non-protected3.3521*0.00306
Hyperprosopon ellipticum protected vs non-protected1.04170.313
Wilcoxon Rank-Sum TestWp
Paralabrax clathratus protected vs non-protected357*2.04x10-5
Girella nigricans protected vs non-protected205*0.0034
ANOVAdfSum Sq.Mean Sq.fp
MPA vs non-MPASpecies2329.21.271.87*0.0141
Protected status132.332.347.5*1.39x10-1-
Schooling vs non-schoolingSpecies1211.40.9491.810.0642
Protected Status10.210.210.3990.528
Schooling115.315.329.1*9.03x10-7
Protected sitesSpecies2013.70.6811.600.0821
Site30.7320.2440.5710.636
Non-protected sitesSpecies1220.211.681.550.133
Site32.410.8050.7420.531

Names, descriptions, and summary statistics for all statistical analyses run in this study.

Protected, non-protected, and status refer to the protection status of an area (marine protected area vs. non-protected area). The ANOVA tests also gave results for the interactions between two variables, but we did not include them in this table as there were no significant interactions. * indicates significance.

Figure 4

4 Discussion

The purpose of this project was to evaluate the effectiveness of Southern California MPAs in promoting healthier fish communities. Using FID as a metric for health, we hypothesized that fish in MPAs would have lower FIDs than those in non-protected areas, indicating higher levels of fish boldness. Across the four protected areas studied, fish had lower FIDs than fish in the four non-protected sites to which they were compared.

There was no significant difference in FID between schooling fish (n=27) and solitary fish (n=72). These results go against a study that found that schooling decreased FIDs in reef fish (Stankowich and Blumstein, 2005), but are consistent with results found by other studies which showed schooling did not affect the FIDs of fish (; Samia et al., 2019). Adult fish are less likely to school in the first place, as schooling can lead to more competition for resources, which larger fish tend to avoid (). It could be that individuals engaged in what we defined as schooling behavior (i.e., four or more individuals in a group), were not functioning as a group but just happened to be found in the same place.

Fish in protected areas had shorter FIDs on average than those in non-protected areas (Figure 2A). This is consistent with our hypothesis that fishing pressures in non-protected areas would lead to increased wariness in fish populations and increases in fish FIDs outside of MPAs. The increase in human-based predator presence can lead to increased responsiveness from potential prey (Milinski, 1986). This increased vigilance of fish found outside MPAs can lead prey to change their foraging behavior to dedicate less energy to foraging and more to finding safety (). Fish in protected areas are not forced to deal with extra predatory pressures from humans, suggesting they can dedicate more energy to foraging and will not be forced to flee from potential predators at such a large distance. This means that fish in protected areas may have increased fitness. Decreased foraging behavior can lead to slower growth rates of individuals and populations in the long term (), which can be correlated to lower community health (). These results are consistent with other metrics of fish community health inside of MPAs, such as increases in abundance and biomass (; ). And the increases in bold behavior of fish inside of MPAs should lead to increases in the other metrics of fish community health (; Sørdalen et al., 2022). As such, FID represents a new and useful tool for evaluating the effectiveness of the MPA network in Southern California. The evaluation of fish behavior can add another dimension to monitoring methods used on the California MPA network, which historically have used physical metrics to assess effectiveness. This work is significant since behavioral studies have been constrained to one species in one location ().

When comparing FIDs of the six core species in protected versus non-protected areas, only Paralabrax Clathratus, Kyphosus azureus, and Girella nigricans had significant differences in FID between populations (Table 5). All three species are considered important game fish and are targeted by both line and spear fishermen along the Southern California coast (Sommer et al., 1995; Stenstrom, 2018; ; ; ). All three of the other species are taken by fishermen, but Hyperprosopon ellipticum and Embiotica jacksoni tend to be smaller than the other fish in the core group and are bottom feeders, meaning that their flesh is considered to be lower quality than the other fish in this grouping (; ). As a result, Hyperprosopon ellipticum and Embiotica jacksoni could be targeted less in non-protected areas, meaning that they would experience less predatory pressure from the lack of protection. Anisotremus davidsonii also had no significant differences between populations in protected and non-protected areas (Table 5). Anisotremus davidsonii is consistently targeted by spearfishermen, as it tends to be larger than Hyperprosopon ellipticum and Embiotica jacksoni and has meat that is valued more highly (). The lack of differences between protected and non-protected populations could be due to smaller sampling sizes, as there were only seven approaches in protected areas and three in non-protected areas, which were the smallest sample sizes of the core species (Table 4).

Fish behavior was consistent across different sites for both protected and non-protected areas (Figure 3). There were no significant differences between sites of the same protected status; however, Tablerock had an average FID that was less than half of the FIDs of the other protected areas (Table 3). Access to Tablerock Beach is difficult and the location is primarily known by Laguna Beach locals for being secluded and calm (). This could be the reason why the fish there seem to be less responsive to people.

There was no difference between the two sites that were partially protected and the two completely protected sites (Figure 3). There is debate about the effectiveness of partially protected marine reserves around the world (Zupan et al., 2018). While they are established with the same goal of promoting healthier marine communities, their effectiveness in accomplishing this goal has been mixed (Sciberras et al., 2013; Turnbull et al., 2021). There is some evidence that they provide increased benefit to populations when compared to completely open areas (; Sciberras et al., 2013; Sciberras et al., 2015), but other data has suggested that they are no better than completely unregulated sites (; Turnbull et al., 2021). There are many types of partially protected areas, and this variation can be correlated to the variation in effectiveness (Zupan et al., 2018). Higher levels of protection resulted in greater biomass and abundance in partially protected areas (Zupan et al., 2018). For behavioral responses, there was no difference between fully and partially protected sites within our study. However, it is important to note that the partially protected sites in this study, Tablerock Beach and El Matador, both heavily limit the species that can be captured within their boundaries, and none of the species sampled in this study were allowed to be taken by fishermen in those areas (). That could mean these function as completely protected areas for the species studied. However, a study on the response of non-target fish to protection status found that even species that spear fishermen do not target respond with shorter FIDs in MPAs versus non-protected areas (Tran et al., 2016). This same study also found that a non-target species responded differently in the presence of a speargun, with longer FIDs when approached with a spear than without one (Tran et al., 2016). It is also important to evaluate the effects of spearfishing and protection on non-target species. This work could contribute to understanding the effectiveness and drawbacks of partially protected areas.

Expanding this study to more MPAs along the California coast would be beneficial, especially to the northern part of the state. Since MPAs in that area have shown slower recovery rates than those in Southern California (Sullivan-Stack et al., 2022), understanding how MPAs affect fish behavior in their boundaries could lend useful information to site managers. It would also be useful to see how protection status changes the behavior of fish over time, not only in comparison to unprotected sites. Since the MPAs in Northern California were implemented more recently (Sullivan-Stack et al., 2022), current studies could act as a baseline, and future studies incorporating FID would give insights on changes in behavior over time, further broadening our understanding of MPA effectiveness. Models have shown that it could take over a decade for changes in biomass and abundance to become noticeable (), but there is no timeline for how long behavioral changes could take. FID should also be used to evaluate the effects of MPAs on the areas around them. MPAs have been shown to increase the amount of fish around them (i.e., spillover effect), which allows for the increased catch of fishermen (McClanahan and Mangi, 2000; Ziegler et al., 2022). There is also evidence that behavioral changes are observed in areas adjacent to MPAs, with a study showing that for a certain distance outside of an MPA, fish continued to have lower FIDs than fish that weren’t near protected areas (). However, this study was done in a tropical setting, and this type of work should also be done in temperate settings to see whether the pattern holds true in these types of ecosystems. By evaluating fish wariness in fish spilling over MPAs, managers will have more metrics to ensure connectivity between MPAs in the MPA network along the California coastline. Finally, this study shows the feasibility of using fish FIDs on the west coast of the Americas, which has not been done before (Nunes et al., 2018). FID is a useful tool for evaluating the effects of fishing on fish behavior and can be used by site managers to evaluate and help create effective marine reserve networks.

It could also be beneficial to see how acclimatization to human interactions changes FIDs of fish in non-fished areas. Terrestrial species have been shown to adapt to human interactions and have lower FIDs in areas with high human interaction (; McGowan et al., 2014; ). These species are not being hunted by humans in urban areas, meaning that there is no need for increased FIDs due to the threat of human predation. It would be helpful to see if fish in MPAs with high human traffic behave the same way. While some studies have looked at how the use of an area affects fish FIDs (), there are still gaps in the literature such as how this can change with degrees of use and in different environments. Because our study focused on human predation as the changing factor of FID, we were not able to speak to the effects of general human interactions on fish FID.

Statements

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

Ethical approval was not required for the study involving animals in accordance with the local legislation and institutional requirements because fish were observed in their natural environment and were not harmed by data collector’s presence.

Author contributions

LH: Conceptualization, Data curation, Formal Analysis, Investigation, Methodology, Project administration, Validation, Visualization, Writing – original draft. FL: Data curation, Funding acquisition, Project administration, Supervision, Writing – review & editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. Funding was received from the NSF grants: DBI 1950350, DBI 2226092, and DEB 2312723. This study was also funded by the Pepperdine University’s Seaver Dean’s Office and Division of Natural Sciences.

Acknowledgments

We would also like to thank all those who helped with data collection including Leon Yang, Conor Kramer, Bryan Kaleta, Julian Jacobs, Sydney Soquet, Caleb Nieleback, Max Sprute, and Haley Hoidal. We would finally like to thank Dr. Helen Holmlund, Dr. Javier Monzon, Mharie Brigham, and Julian Jacobs for reviewing the manuscript and providing edits.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmars.2023.1270603/full#supplementary-material

Supplementary Table 1

The number of approaches per species for all locations.

References

Summary

Keywords

MPA (marine protected area), flight initiation distance (FID), California, fish, health, community, behavior, fitness

Citation

Himes L and La Valle FF (2023) Southern California marine protected areas promote bolder fish populations. Front. Mar. Sci. 10:1270603. doi: 10.3389/fmars.2023.1270603

Received

01 August 2023

Accepted

06 October 2023

Published

20 October 2023

Volume

10 - 2023

Edited by

Carlo Pipitone, CNR - IAS Institute of Anthropic Impacts and Sustainability in Marine Environment, Italy

Reviewed by

Benjamin Toscano, Trinity College, United States; José Lino Vieira De Oliveira Costa, University of Lisbon, Portugal

Updates

Copyright

*Correspondence: Florybeth Flores La Valle,

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics