Marine Pollution Bulletin 146 (2019) 502–508

Contents lists available at ScienceDirect

Marine Pollution Bulletin

journal homepage: www.elsevier.com/locate/marpolbul

Effects of chronic pesticide exposure on an epibenthic oyster reef community T

a,b,⁎ b b,c
Samantha M. Iliff , Rachel J. Harris , Elizabeth W. Stoner
a
Harriet L. Wilkes Honors College, Florida Atlantic University, 5353 Parkside Drive, Jupiter, FL 33458, USA
b
Loxahatchee River Environmental Control District, 2500 Jupiter Park Drive, Jupiter, FL 33458, USA
c
Department of Natural and Applied Sciences, Bentley University, 175 Forest Street, North Waltham, MA 02452, USA

ARTICLE INFO ABSTRACT

Keywords: In recent decades, oyster reefs have been deteriorating throughout North America as a result of multiple in-
Americorophium amphipod teracting anthropogenic stressors, including pesticide pollution. Here we elucidated the potential chronic effects
Carbaryl of the commonly utilized pesticide, carbaryl, on oyster reef communities in the Loxahatchee River Estuary in
Crustacean southeast Florida. Though carbaryl had a limited effect on total epifaunal community diversity, species richness
Pesticide
and evenness, the results of this experiment indicate that carbaryl significantly shifted crustacean community
Loxahatchee River Estuary
Oyster reef
composition, resulting in a substantial loss in total crustacean abundance. One crustacean in particular,
Americorophium spp. (tube building amphipod), was significantly less abundant within the carbaryl treatment,
driving the shift in crustacean community composition. Ultimately, our results signal that pesticide pollution in
estuaries will negatively impact crustaceans. Over time, this may shift benthic community composition, po-
tentially disrupting species interactions and threatening valuable economic and ecosystem services.

1. Introduction complete loss of function (Beck et al., 2011).

One stressor contributing to this decline may be pollution in the
Coastal marine and estuarine ecosystems are among the most pro- form of pesticide runoff. There is growing evidence that pesticides can
ductive ecosystems globally (Kennish, 2002; Luisetti et al., 2014). Oy- accumulate to toxic levels in coastal marine habitats and degrade
ster reefs in particular play a critical role in providing valuable eco- benthic invertebrate communities (Bay et al., 2005; Anderson et al.,
system services such as carbon sequestration, filtration of water along 2007; Ranasinghe et al., 2007; Holmes et al., 2008; Werner et al., 2010;
with improvement of water quality, and habitat provisioning for epi- Phillips et al., 2010a; Phillips et al., 2010b; Lao et al., 2010; Anderson
fauna including crustaceans, mollusks, and other oysters (Grabowski et al., 2014). However, there are few studies describing the effects of
and Peterson, 2007; Layman et al., 2014). Further, oyster reefs harbor a chronic pesticide exposure on oyster reef communities in tropical-sub-
wide diversity of small cryptic species that stimulate secondary pro- tropical estuaries. There are several mechanistic pathways by which
duction and drive the transfer of energy within oyster reefs (Edgar and pesticides are transported into estuaries. Pesticides are widely used in
Aoki, 1993; Adams, 1976; Yeager and Layman, 2011; Abeels et al., agricultural settings as well as in many homes and residential gardens
2012). Because of their proximity to developed coastlines and urban and may enter estuarine and coastal environments in stormwater runoff
centers, oyster reefs are particularly susceptible to non-point and point following heavy rain events, or through daily water use activities
source nutrient and chemical pollution from wastewater and storm- (Stoner and Arrington, 2017). Once in a water body, pesticides can
water runoff, septic tank discharges, and groundwater flow, which can exert both lethal and sub-lethal effects on aquatic organisms. For ex-
have deleterious effects on population, community, and ecosystem ample, sediment and water samples collected from a Californian estuary
dynamics (Kennish, 2002). As a result, oyster reefs have been rapidly contained highly toxic levels of organophosphate and pyrethroid pes-
declining throughout the world, largely due to anthropogenic activities ticides driving mortality in various crustaceans in acute toxicity
acting in concert with shifting environmental conditions (e.g., fresh- bioassessments (Anderson et al., 2014). Additionally, Schroeder-Spain
water discharge) (Jackson, 2008; Beck et al., 2011; Zu Ermgassen et al., et al. (2018) examined the effects of pesticides on the blue crab (Cal-
2012). An estimated 85% of oyster reefs have been lost globally, and linectes sapidus) which is both an important predator and prey species in
many North American bays have lost over 99% of oyster reefs with estuarine systems. Following acute exposure to pesticides, blue crab

⁎
Corresponding author currently at: Department of Ecology and Evolutionary Biology, University of Michigan, 1105 North University Ave., Ann Arbor, MI 48109,
USA
E-mail addresses: siliff@umich.edu (S.M. Iliff), rachel.harris@lrecd.org (R.J. Harris), estoner@bentley.edu (E.W. Stoner).

https://doi.org/10.1016/j.marpolbul.2019.06.060
Received 1 March 2019; Received in revised form 15 June 2019; Accepted 20 June 2019
Available online 10 July 2019
0025-326X/ © 2019 The Authors. Published by Elsevier Ltd. This is an open access article under the CC BY license
(http://creativecommons.org/licenses/by/4.0/).
S.M. Iliff, et al. Marine Pollution Bulletin 146 (2019) 502–508

mortality occurred at high concentrations of pesticides (100 μg/L), and plaster blocks containing carbaryl. The carbaryl blocks were made by
exposure to low concentrations (1–10 μg/L) had sub-lethal effects on mixing Plaster of Paris, 80% carbaryl wettable powder and cold water
the crabs, resulting in altered behaviors, decreased muscle function, (Whalen et al., 2013). As a procedural control, additional plaster blocks
and increased reaction time (Schroeder-Spain et al., 2018). As such, were made using the same method without the addition of carbaryl.
pesticide pollution can influence invertebrate populations directly via The blocks were allowed to dry and harden and were then removed
toxicity or indirectly through decreased foraging, or by increasing from their molds and affixed to strings (Fig. 2).
predation risk (Schroeder-Spain, 2017). Although there is a need to To test the effects of carbaryl on oyster reef community composition
understand the effects of acute pesticide toxicity on valuable species, in situ, we designed a field experiment to quantify how carbaryl ex-
few studies have examined the effects of chronic pesticide exposure on posure may influence epifaunal community composition. We tested
estuarine communities. three treatments in this experiment: carbaryl additions (using carbaryl
One common pesticide in the United States that could affect benthic dissolution blocks), a plaster-control treatment (using dissolution
estuarine organisms is carbaryl, a broad-spectrum organocarbamate blocks made without carbaryl), and a double control group with no
(Donaldson et al., 2002). Carbaryl works by inhibiting acet- treatment of dissolution blocks. This double control method was im-
ylcholinesterase in arthropods, causing death by disrupting the nervous plemented to analyze the procedural effects of plaster dissolution on
system (Dumbauld et al., 2001; Relyea, 2005). Commercially available response variables. To simulate the hard-bottom microhabitat of an
as Sevin Insecticide, carbaryl was ranked as the number one most oyster reef, a string of five, 100 cm2 travertine tiles was hung off of T-
commonly used insecticide in the U.S. home and garden market (Grube shape structures constructed with PVC pipes (Fig. 2). We elected to use
et al., 2011). Although insecticides are more commonly used in ter- non-invasive removable tiles rather than a natural oyster reef in order
restrial systems to deter insect predation on plants, carbaryl has also to easily collect and quantify mobile and highly cryptic reef species
been used in intertidal oyster farming practices (Simenstad and Fresh, without removing or damaging a section of natural reef. This method
1995). Since the 1960s oyster growers in Willapa, Bay (Washington accurately simulates an oyster reef, allowing for oyster spat recruitment
State, USA), have treated farmed oyster beds with carbaryl in alter- and settlement (Metz et al., 2015). Each tile was vertically separated by
nating cycles of short-term exposure to protect their stock from bur- a 6.35 mm spacer. The lower tile of each string hung 30 cm above the
rowing thalassinid shrimp (Feldman et al., 2000; Dumbauld et al., sediment in the shallow sub-intertidal zone, and the top-most tile of
2001). Numerous studies have evaluated the lethal and sub-lethal ef- each string was at least 30 cm below the surface at low tide. Each dis-
fects of carbaryl on marine organisms (Estes, 1986; Fulton and Key, solution block was then hung directly adjacent (~3 cm) to a string of
2001; Canty et al., 2007; Schroeder-Spain et al., 2018); however, there recruitment tiles.
is a paucity of data regarding effects of chronic carbaryl exposure in Eighteen recruitment structures were deployed 2 m apart at the field
estuarine ecosystems, especially oyster reefs. It is critical to identify the site. The 2 m structure separation was chosen based on previously es-
impacts of chronic exposure of pesticides because it is likely that pro- tablished methods to locally increase carbaryl concentrations relative to
longed exposure to pesticides could substantially degrade ecosystems controls (see Poore et al., 2009; Cook et al., 2011; and Whalen et al.,
though cascading and multiple-stressor effects. To this end, the goal of 2013). In total, thirty-six strings of recruitment tiles were hung, with
this study was to determine whether chronic carbaryl exposure would one on each end of each T structure. Strings were randomly assigned to
alter oyster reef community composition in a sub-tropical estuarine a treatment, and then assigned to a T structure. Dissolution blocks were
waterbody. In having a better understanding of the effects of carbaryl replaced every 3–5 days and the time and date were recorded. Every
exposure on oyster reefs, we are providing resource managers with a dissolution block was weighed before deployment and after retrieval to
holistic understanding needed to fully assess anthropogenic impacts on measure dissolution rates. Sessile and epibenthic organisms were al-
oyster reefs. lowed to colonize the tiles for 28 days from September 13th–October
11th 2017, which is a time frame that allows for recruitment and set-
2. Methods tlement of eastern oyster (Crassostrea virginica) spat in the Loxahatchee
River estuary (Metz et al., 2015).
2.1. Study area
2.3. Sample collection, processing, and data analysis
We conducted this experiment in the Southwest Fork of the sub-
tropical Loxahatchee River Estuary in Jupiter, Florida (Fig. 1). At this After the 28-day experiment, both sessile and vagile fauna asso-
site, salinities ranged from 0.1 to 34.7 ppt with a mean of 11.2 ppt and ciated with the recruitment tiles were collected by gently enclosing the
water temperatures ranged from 25 to 32 °C (available at: https:// entire string of tiles in a gallon bag before removing the T structures. All
loxahatcheeriver.org/river/river-keeper/), providing suitable oyster unattached organic material was rinsed off of the tiles and strained
habitat requirement (Barnes et al., 2007). The Loxahatchee River Es- through a 500 μm sieve. All strained material was preserved in 70%
tuary covers roughly 4 km2, drains a watershed of approximately ethanol and stained with a 10% Rose Bengal dye solution. All of the
700 km2 and flows into the Atlantic Ocean at the Jupiter Inlet (Ridler epifauna were then identified to the lowest taxonomic level practical,
et al., 2006). The watershed covers a variety of habitats including riv- quantified, and recorded. After the tiles had been rinsed, they were
erine wetland and cypress communities, and estuarine mangroves, oy- allowed to dry overnight. A dissecting microscope was used to analyze
ster reefs, and seagrass beds. Although the system is currently flow and enumerate all sessile fauna on the middle three tiles as outlined in
managed, the combination of altered freshwater inputs and deterior- Haven and Fritz (1985).
ating water quality (see Stoner and Arrington, 2017) have driven shifts The initial and final weights of the dissolution blocks, as well as the
in oyster bed community composition and ecosystem function (Van date and time of deployment and collection, were then used to de-
Arman et al., 2005; Roberts et al., 2008; Parker et al., 2013; Jud, 2014), termine the rate of dissolution of blocks (Thompson and Glenn, 1994)
emphasizing a need for effective resource management and restoration. and dissolution of carbaryl (80% by dry weight). The following equa-
Although water quality is regularly monitored in the Loxahatchee River tion was used where D is dissolution rate, Wi is initial weight, Wf is final
Estuary, monitoring for pesticides has been limited. weight, and θ is total time (in days) dissolving:
Wi Wf
2.2. Experimental design D=

In order to test the effects of carbaryl on an estuarine oyster reef, we Dcarbaryl = D × 0.8 × 0.10
used a method developed by Poore et al. (2009) using slow-dissolving

503
S.M. Iliff, et al. Marine Pollution Bulletin 146 (2019) 502–508

Fig. 1. Map of the study site, indicated with a bullseye, within the Loxahatchee River Estuary, Jupiter, FL.

When comparing overall mean density between treatments, each (set at 999 unique permutations) were used to identify significant dif-
tile string was treated as a discrete sample in which mean standardized ferences between treatments and similarity percentages (SIMPER) were
faunal densities were calculated across the middle three tiles for each used to verify the individual species driving community differences
tile string. One-way ANOVAs were used to test for difference in in- between treatments (Anderson, 2001). Community composition ana-
dividual species densities between treatments, and the non-parametric lysis was conducted using Primer PERMANOVA+v7.0 and all other
test, Kruskal-Wallis test, was used when unequal variance was identi- analyses using with SPSS v24 statistical software.
fied. p-Values ≤ 0.05 were considered statistically significant, while p-
values between 0.05 and 0.10 were considered marginally significant. 3. Results
Community composition dissimilarities were analyzed on a com-
munity wide basis (total fauna) and within the crustacean community. The initial and final weights of the dissolution blocks, as well as the
Multidimensional scaling (MDS) plots were used to visualize the simi- date and time of deployment and collection, were used to determine the
larity and dissimilarity of macrofaunal community compositions (based rate of dissolution of blocks as well as dissolution of carbaryl. When the
on a Bray-Curtis resemblance matrix) by carbaryl, plaster, and control carbaryl dissolution blocks were deployed, each contained approxi-
treatments. The vectors on the MDS plots (set at a Pearson's R Value of mately 15 g of carbaryl by dry weight and ended with approximately
0.75) were used to identify which species were driving the similarities 8 g of carbaryl by dry weight, with an average carbaryl dissolution rate
within and dissimilarity between treatments. Permutational Analysis of of 2.64 ± 0.3 g/day; the average dissolution of carbaryl blocks and
Variance (PERMANOVA) and post-hoc pair-wise PERMANOVA tests plaster blocks was quite similar at 32.9 g/day and 30.0 g/day,

504
S.M. Iliff, et al. Marine Pollution Bulletin 146 (2019) 502–508

Fig. 2. Illustration of site layout showing paired experimental

treatments: control (□), plaster ( ) and carbaryl ( ). Inset
image of individual T-structure with attached recruitment
tiles and carbaryl dissolution blocks.

respectfully. The carbaryl blocks dissolved approximately 9% faster

than the plaster blocks. The differences in dissolution rates in both
plaster and carbaryl ( ± 5 g day) varied between deployment periods,
likely due to changing environmental conditions throughout the 28-day
study period.
After 28 days, 27 taxa were collected from or observed on the re-
cruitment tiles. The most abundant and frequently observed in-
vertebrates were molluscs and arthropods (Table 1). The most abundant
organisms were tube-forming serpulid worms (Hydroides sp.), eastern
oyster (C. virginica), ivory barnacles (Balanus eburneus), gammaridean
amphipods (Hourstonius laguna), and tube building amphipods (Amer-
icorophium spp.). Mean taxa richness, diversity, and evenness (based on
the Shannon-Wiener Diversity Index) were lowest within the carbaryl
treatment, although not significantly so (Kruskal-Wallis test, p > 0.05;
Table 1). The overall mean density of individuals per sample was
slightly lower within the carbaryl treated samples, showing a 5% de-
crease in total fauna compared to control samples. However, within the
crustacean community, the mean density of individuals was 43% less in
the carbaryl treatment compared to the control treatment (Fig. 3).
Multidimensional scaling (MDS) plots were used to visualize the Fig. 3. Mean density of total organisms and crustaceans between treatments.
dissimilarity between and similarities within the community composi- Percentages indicate difference in organism density relative to the control
tions of the carbaryl, plaster, and control treatments where the vectors treatment. Error bars indicate 95% confidence interval.

Table 1
Mean abundance (individuals/m2) ± standard deviation of the most-abundant animals found on recruitment tiles between treatments. Species are listed by clas-
sification and order of abundance within the control treatment. * indicates marginal difference between treatments (ANOVA or Kruskal-Wallis, 0.05 < p < 0.10).
** indicates significant difference between treatments (ANOVA or Kruskal-Wallis, p < 0.05).
Scientific name Common name Control Plaster Carbaryl

Mollusca
C. virginica Eastern oyster 18,560 ± 1506 19,440 ± 4343 20,450 ± 7321
G. demissa Ribbed mussel 1450 ± 633 1060 ± 396 1425 ± 608
Neritina spp. Nerite snail 118 ± 84 152 ± 104 76 ± 39
Arthropoda
**Americorophium spp. Tube building amphipod 4392 ± 3197 2406 ± 765 998 ± 175
*B. eburneus Ivory barnacle 4960 ± 1306 8170 ± 3366 5616 ± 826
H. laguna Amphipod 2952 ± 589 3026 ± 1263 2708 ± 1082
Insecta Insect larvae 2182 ± 1122 1814 ± 2293 1084 ± 760
S. stanfordi Amphipod 1574 ± 740 1402 ± 1059 1540 ± 704
B. amphitrite Purple striped barnacle 976 ± 625 2273 ± 1518 1670 ± 215
*Isopoda Isopod 176 ± 65 190 ± 45 110 ± 49
Other
Membranipora spp. Encrusting Bryozoan 440 ± 203 403 ± 130 496 ± 215
Nereididae Nereid worms 74 ± 30 80 ± 70 68 ± 33
Hydroides sp.a Serpulid worm 55.6% 50.4% 50.6%

Notes: Sessile organisms were quantified on the tile surfaces from the middle three tiles totaling an area of 0.06 m2. The vagile organisms were rinsed off the entire
string of tiles from an area of 0.1 m2.
a
Serpulid worms were quantified using percent cover because densities often exceeded 50,000 individuals/m2.

505
S.M. Iliff, et al. Marine Pollution Bulletin 146 (2019) 502–508

insignificant (pair-wise PERMANOVA, p < 0.10; Fig. 4a). Among

crustaceans, a significant dissimilarity was observed between the car-
baryl group and both control groups (pair-wise PERMANOVA,
p < 0.05; Fig. 4b). The primary drivers of these dissimilarities were the
lack of Americorophium spp. and relative increase in abundance of ta-
naids (Sinelobus stanfordi) within the carbaryl treatment. As indicated
by the vectors within the MDS plots (Fig. 4), one crustacean in parti-
cular, Americorophium spp., drove the dissimilarity in community
compositions between treatments. Independent sample Kruskal-Wallis
test results indicated that Americorophium spp. was the only individual
species to show a significantly lower density within the carbaryl
treatment, compared to the controls (p < 0.05; Fig. 5). The mean
density of Americorophium spp. was 77% less in the carbaryl treatment
when compared to the control, driving 56% of the overall dissimilarity
in community composition between control and carbaryl treatments.
Only one other species, Grandidierella bonnieroides, an amphipod, was
marginally lower (Kruskal-Wallis tests, p = 0.06–0.10) in carbaryl
treatments.

4. Discussion

In this study, we considered chronic carbaryl exposure on oyster

reef fauna in a dynamic estuarine system. Overall, carbaryl influenced
oyster reef communities, though the effects varied by species. The
foundation species in the oyster reef community, C. virginica, was un-
affected by carbaryl at the levels tested in this study. However, carbaryl
exposure had a strong effect on crustacean community composition,
reducing crustacean abundance by 43%. This trend may largely be
explained by the deleterious effects carbaryl exerted on tube-building
amphipods. Interestingly, carbaryl had less of an effect on overall oyster
reef community composition in this experiment.
Short term exposure to pesticides can lead to a reduction in the
Fig. 4. Multidimensional scaling (MDS) plot illustrating dissimilarity in the
activity and viability of cells in filter feeders (e.g., Mytilus edulis; Canty
community composition between treatments. The distance between two points et al., 2007). Here, we did not measure a dramatic decline in the
is proportional to degree of dissimilarity (Pearson's R = 0.7, stress < 0.15). common, filter feeding organisms encountered throughout the study (C.
PERMANOVA used to measure significant difference between treatments. (a) virginica, G. demissa, B. eburneus, B. amphitrite, Hydroides sp. or Mem-
represents dissimilarity between total fauna in the community (p > 0.05). (b) branipora spp.). Future studies should more closely evaluate potential
represents dissimilarity in crustacean community (p < 0.05). effects of chronic pesticide exposure on these organisms, especially C.
virginica, which is a foundation species and facilitates the development
of oyster reef habitat for many species. However, we did measure a
significant decline in Americorophium spp., which filter-feeds on phy-
toplankton (Croteau et al., 2005). Although evaluating specific me-
chanisms driving species-level differences from carbaryl loading were
beyond the scope of this study, it is possible that the mode of filter
feeding specific to Americorophium spp. may have contributed to the
decreased numbers found in the presence of carbaryl. In addition,
previous studies have shown that carbaryl is less stable in aqueous
solutions and becomes concentrated in soils, subjecting organisms
living within sediments to higher concentrations than organisms in the
water column (Kanazawa, 1975; Estes, 1986). Since Americorophium
spp. build their soft tubes out of local sediments and organic materials,
this may explain the particularly strong impact of carbaryl on Amer-
icorophium spp. measured in this study. Whether by specific feeding
mode or sediment reworking activity, it is likely that the physiology and
ecology of Americorophium spp. are driving the differences between
treatments. For instance, Anderson et al. (2014) found that Americor-
ophium spp. in the San Francisco Estuary occurred in less polluted and
least toxic field sites and Thompson et al. (2012) lists Americorophium
spp. as sensitive taxa in tidal freshwater habitats. As investigation into
Fig. 5. Mean density of Americorophium spp. between treatments. * indicates
pesticide pollution proceeds, the specific pesticide tolerances of these
significant difference in density between control and carbaryl groups at the
α = 0.05 level (Kruskal-Wallis test, p = 0.008).
species will need to be measured. Based on the detected sensitivity of
Americorophium spp. observed in this study, and similarly reported in
Anderson et al. (2014), we recommend these species could be used as
indicate drivers of dissimilarity between treatments (Fig. 4). Among the potential bioindicators to detect pesticide pollution, primarily, in oligo-
entire community of fauna found on the tiles, the dissimilarity between mesohaline (0–18 ppt) habitats, where they occur in South Florida
the carbaryl treatment and the control treatment was marginally (Evans et al., 2010).

506
S.M. Iliff, et al. Marine Pollution Bulletin 146 (2019) 502–508

Within oyster reefs, crustaceans are among the most abundant taxa, Pesticides such as carbaryl have the potential to significantly alter
second only to bivalves (Bahr Jr, 1976; Bahr and Lanier, 1981, Wells, crustacean communities, which would have cascading impacts on oy-
1961). It has been well documented that small crustaceans are a key ster reefs, diminishing their economic and ecological value. In order to
link in the transfer of energy between primary producers and secondary better understand these processes, we must place future efforts on the
consumers and that these species compose a large proportion of the effects of combined stressors, including pesticides, on oyster reef
diets of juvenile fish and other mid-level predators (Adams, 1976; Edgar community composition and overall ecosystem function.
and Aoki, 1993; Yeager and Layman, 2011; Abeels et al., 2012). A
stable isotope study by Abeels et al. (2012) revealed that crustaceans CRediT authorship contribution statement
are included in all lower trophic levels and were consumed by every
secondary consumer studied on a sub-tropical oyster reef. An in- Samantha M. Iliff: Conceptualization, Data curation, Formal ana-
troduction of pesticides into an estuarine ecosystem could potentially lysis, Investigation, Methodology, Project administration, Resources,
impact the abundance and composition of crustacean communities. A Software, Supervision, Visualization, Writing - original draft, Writing -
shift at this trophic level could result in cascading effects throughout review & editing. Rachel J. Harris: Conceptualization, Methodology,
the ecosystem including a reduction in biomass, a loss of diversity, and Project administration, Resources, Supervision, Writing - review &
a decrease in available prey for fish and other secondary consumers. editing. Elizabeth W. Stoner: Formal analysis, Investigation,
Although there were no significant differences in species diversity Methodology, Project administration, Resources, Software, Supervision,
and richness between the treatments, there are several mechanisms that Visualization, Writing - review & editing.
could potentially explain these variable results. Based on the average
salinity (11.2 ppt) and the dissolved fraction of the dissolution blocks, Acknowledgements
we estimate an average flow velocity of 15 cm/s over the course of this
study (Thompson and Glenn, 1994), which falls within the average The authors wish to thank the Loxahatchee River Environmental
range 5–25 cm/s of flow velocities reported on oyster reefs (Reidenbach Control District, United States for funding. We also wish to thank
et al., 2013), and within the maximum current velocities (17 cm/s) Katherine Banger, Owen Silvera, Aaron Kent, Kyle Vincent, Liberty
previously measured in the South West fork of the Loxahatchee River Boyd, Susan Noel, Helen Johnson, David Porter, and Jerry Metz for field
Estuary (Patra, 2003). Although low compared to high energy systems, and data assistance; Eric Lazo-Wasem for his assistance identifying
over time, these flow rates would have caused the carbaryl to dissipate epifauna; Jon Moore and Jim Wetterer for edits to the original under-
rapidly; averaging 2.6 ( ± 0.3) g carbaryl/day. Thus, we regularly re- graduate thesis manuscript; and an anonymous review for their com-
placed the blocks (every 3–5 days) to achieve chronic exposure over the ments and suggestionson.
28-day study period. We did not measure carbaryl in ambient waters
during this study. However, with a solubility of 110,000 μg/L in es- Funding
tuarine waters (Poore et al., 2009), and a dissolution rate of 2.6 g/day,
we estimate that the average concentrations of carbaryl were approxi- This work was supported by the Loxahatchee River District.
mately 8800 μg carbaryl/L; although the concentrations would have
varied greatly throughout the tidal cycle (0–110,000 μg/L). On average Declaration of Competing Interest
the carbaryl concentrations reported here are much higher than 100 μg/
L reported in well-mixed estuarine waters (Schroeder-Spain, 2017), and None.
even on the high end of acute (24–48 h) exposure limits reported for
invertebrate species (30–7200 μg carbaryl/L; Feldman et al., 2000). References
These mean concentrations of carbaryl (8800 μg/L) seem quite high
compared to the acute toxicity studies, yet, unlike the acute toxicity Abeels, H.A., Loh, A.N., Volety, A.K., 2012. Trophic transfer and habitat use of oyster
studies, this experiment was not conducted in a closed system. Thus, Crassostrea virginica reefs in southwest Florida, identified by stable isotope analysis.
Mar. Ecol. Prog. Ser. 462, 125–142.
despite regularly replacing the carbaryl dissolution blocks it is likely Adams, S.M., 1976. Feeding ecology of eelgrass fish communities. Trans. Am. Fish. Soc.
that the variations in flow direction and velocities dissipated the ef- 105 (4), 514–519.
fectiveness of our carbaryl treatments (Whalen et al., 2013). Therefore, Anderson, M.J., 2001. A new method for non-parametric multivariate analysis of var-
iance. Austral. Ecol. 26 (1), 32–46.
it is possible that some applications of carbaryl could enter the oyster Anderson, B.S., Hunt, J.W., Phillips, B.M., Tjeerdema, R.S., 2007. Navigating the TMDL
reef at higher concentrations than those used here. In turn, this would Process: Sediment Toxicity. 02-WSM-2. Water Environment Research Foundation, pp.
magnify impacts to the oyster reef community. 194.
Anderson, B., Phillips, B., Hunt, J., Siegler, K., Voorhees, J., Smalling, K., Kuivila, K.,
Oyster reefs are one of the most ecologically and economically va- Hamilton, M., Ranasinghe, A., Tjeerdema, R., 2014. Impacts of pesticides in a Central
luable habitats within estuaries. It is not well known why oyster reefs California estuary. Environ. Monit. Assess. 186 (3), 1801–1814.
are in decline, yet it is apparent that a combination of stressors could be Armhurst, K.L., Grosby, D.G., 1991. Fate of carbaryl, 1-naphthol and atrazine in seawater.
Pac. Sci. 45, 314–320.
in effect. Here, we were able to detect the negative effects of chronic
Bahr Jr., L.M., 1976. Energetic aspects of the intertidal oyster reef community at Sapelo
pesticide exposure on crustaceans. Additional studies must focus on Island, Georgia (USA). Ecology 57 (1), 121–131.
acute versus chronic effects of pesticides on estuarine communities. Our Bahr, L.M., Lanier, W.P., 1981. The Ecology of Intertidal Oyster Reefs of the South
exposure time was limited to 28 days and it is possible that carbaryl Atlantic Coast: A Community Profile (No. 81/15). US Fish and Wildlife Service.
Barnes, T.K., Volety, A.K., Chartier, K., Mazzotti, F.J., Pearlstine, L., 2007. A habitat
could have greater impacts on estuarine community composition over a suitability index model for the eastern oyster (Crassostrea virginica), a tool for re-
larger time period. Furthermore, carbaryl could be a relatively mild storation of the Caloosahatchee Estuary, Florida. J. Shellfish Res. 26 (4), 949–960.
pesticide in marine environments with variable toxicity and a relatively Bay, S.M., Mikel, T.M., Schiff, K., Mathison, S., Hester, B., Young, D., Greenstein, D.,
2005. Southern California Bight Regional Monitoring Program: i. Sediment Toxicity.
short half-life (Armhurst and Grosby, 1991). Future work must consider Technical Report 451. Southern California Coastal Water Research Project,
the response of both individual species and the entire community to Westminster.
various degrees of pesticides to determine specific toxicities and toler- Beck, M.W., Brumbaugh, R.D., Airoldi, L., Carranza, A., Coen, L.D., Crawford, C., Defeo,
O., Edgar, G.J., Hancock, B., Kay, M.C., Lenihan, H.S., 2011. Oyster reefs at risk and
ances. Additionally, mixtures of pesticides commonly used in agri- recommendations for conservation, restoration, and management. Bioscience 61 (2),
cultural and gardening practices including pyrethroids and organo- 107–116.
phosphates could increase the toxicity to marine species, lingering in Canty, M.N., Hagger, J.A., Moore, R.T.B., Cooper, L., Galloway, T.S., 2007. Sublethal
impact of short term exposure to the organophosphate pesticide azamethiphos in the
the water and sediment, or bioaccumulating through trophic levels. marine mollusc Mytilus edulis. Mar. Pollut. Bull. 54 (4), 396–402.
Future studies should incorporate other pesticides such as herbicides Cook, K., Vanderklift, M.A., Poore, A.G., 2011. Strong effects of herbivorous amphipods
and algaecides to consider the impacts on sensitive marine habitats. on epiphyte biomass in a temperate seagrass meadow. Mar. Ecol. Prog. Ser. 442,

507
S.M. Iliff, et al. Marine Pollution Bulletin 146 (2019) 502–508

263–269. Res. 32 (3), 695–708.

Croteau, M.N., Luoma, S.N., Stewart, A.R., 2005. Trophic transfer of metals along Patra, R.R., 2003. Sediment Management in Low-energy Estuaries: Loxahatchee, Florida
freshwater food webs: evidence of cadmium biomagnification in nature. Limnol. (Doctoral dissertation). University of Florida.
Oceanogr. 50 (5), 1511–1519. Phillips, B.M., Anderson, B.A., Hunt, J.W., Siegler, K., Voorhees, J.P., McNeill, K., 2010a.
Donaldson, D., Kiely, T., Grube, A., 2002. Pesticides Industry Sales and Usage 1998 and Pyrethroid and organophosphate pesticide-associated toxicity in two coastal water-
1999 Market Estimates. US Environmental Protection Agency, Office of Prevention, sheds (California, USA). Environ. Toxicol. Chem. 31, 1595–1603.
Pesticides, and Toxic Substances, Office of Pesticide Programs, Washington, DC. Phillips, B.M., Anderson, B.S., Voorhees, J.P., Hunt, J.W., Holmes, R.W., Mekebri, A.,
Dumbauld, B.R., Brooks, K.M., Posey, M.H., 2001. Response of an estuarine benthic Connor, V., Tjeerdema, R.S., 2010b. The contribution of pyrethroid pesticides to
community to application of the pesticide carbaryl and cultivation of Pacific oysters sediment toxicity in four urban creeks in California, USA. J. Pestic. Sci. 35, 302–309.
(Crassostrea gigas) in Willapa Bay, Washington. Mar. Pollut. Bull. 42 (10), 826–844. Poore, A.G., Campbell, A.H., Steinberg, P.D., 2009. Natural densities of mesograzers fail
Edgar, G.J., Aoki, M., 1993. Resource limitation and fish predation: their importance to to limit growth of macroalgae or their epiphytes in a temperate algal bed. J. Ecol. 97
mobile epifauna associated with Japanese Sargassum. Oecologia 95 (1), 122–133. (1), 164–175.
Estes, P.S., 1986. Cardiovascular and Respiratory Responses of the Ghost Shrimp, Ranasinghe, J.A., Barnett, A.M., Schiff, K., Montagne, D.E., Brantley, C., Beegan, C.,
Callianassa Californiensis Dana, to the Pesticide Carbaryl and its Hydrolytic Product Cadien, D.B., Cash, C., Deets, G.B., Diener, D.R., Mikel, T.M., Smith, R.W., Velarde,
1-naphthol. M.S (Thesis). Oregon State University, Corvallis, Oregon (digital ar- R.G., Watts, S.D., Weisberg, S.B., 2007. Southern California Bight 2003 Regional
chives). Monitoring Program: III. Benthic Macrofauna. Technical Report 529. Southern
Evans, D.L., Strom, D.G., Mosura-Bliss, E.L., 2010. Spatial distribution of benthic mac- California Coastal Water Research Project, Costa Mesa.
roinvertebrates in the Pithlachascotee River during low-flow conditions with em- Reidenbach, M.A., Berg, P., Hume, A., Hansen, J.C.R., Whitman, E.R., 2013.
phasis on relationships with salinity. In: Technical Report. Water & Air Research, Inc. Hydrodynamics of intertidal oyster reefs: the influence of boundary layer flow pro-
(For Southwest Florida Water Management District). cesses on sediment oxygen exchange. Limnol. Oceanogr. 3, 225–239.
Feldman, K.L., Armstrong, D.A., Dumbauld, B.R., SeWitt, T.H., Doty, D.C., 2000. Oysters, Relyea, R.A., 2005. The impact of insecticides and herbicides on the biodiversity and
crabs, and burrowing shrimp: review of an environmental conflict over aquatic re- productivity of aquatic communities. Ecol. Appl. 15 (2), 618–627.
sources and pesticide in Washington State's (USA) coastal estuaries. Estuaries 23 (2), Ridler, M.S., Dent, R.C., Arrington, D.A., 2006. Effects of two hurricanes on Syringodium
141–176. filiforme, manatee grass, within the Loxahatchee River estuary, Southeast Florida.
Fulton, M.H., Key, P.B., 2001. Acetylcholinesterase inhibition in estuarine fish and in- Estuar. Coasts 29 (6), 1019–1025.
vertebrates as an indicator of organophosphorus insecticide exposure and effects. Roberts, R.E., Hedgepeth, M.Y., Alexander, T.R., 2008. Vegetational responses to salt-
Environ. Toxicol. Chem. Int. J. 20 (1), 37–45. water intrusion along the Northwest Fork of the Loxahatchee River within Jonathan
Grabowski, J.H., Peterson, C.H., 2007. Restoring oyster reefs to recover ecosystem ser- Dickinson State Park. Florida Sci. 71 (4), 383–397.
vices. Ecosyst. Eng. Plants Protists 4, 281–298. Schroeder-Spain, K.J., 2017. Dazed and Confused: Pesticides Alter Physiology, Behavior,
Grube, A., Donaldson, D., Kiely, T., Wu, L., 2011. Pesticides Industry Sales and Usage. US and Predator–Prey Interactions of Juvenile and Adult Blue Crabs (Callinectes sapidus)
EPA, Washington, DC. (Doctoral dissertation). Texas A&M University, Corpus Christi.
Haven, D.S., Fritz, L.W., 1985. Setting of the American oyster Crassostrea virginica in the Schroeder-Spain, K., Fisher, L.L., Smee, D.L., 2018. Uncoordinated: effects of sublethal
James River, Virginia, USA: temporal and spatial distribution. Mar. Biol. 86 (3), malathion and carbaryl exposures on juvenile and adult blue crabs (Callinectes sa-
271–282. pidus). J. Exp. Mar. Biol. Ecol. 504, 1–9.
Holmes, R.W., Anderson, B.S., Phillips, B.M., Hunt, J.W., Crane, D., Mekebri, A., Simenstad, C., Fresh, K., 1995. Influence of intertidal aquaculture on benthic commu-
Blondina, G., Nguyen, L., Connor, V., 2008. Statewide investigation of the role of nities in Pacific Northwest estuaries: scales of disturbance. Estuaries 18 (1), 43–70.
pyrethroid pesticides in sediment toxicity in California's urban waterways. Environ. Stoner, E.W., Arrington, D.A., 2017. Nutrient inputs from an urbanized landscape may
Sci. Technol. 42, 7003–7009. drive water quality degradation. Sustain. Water Qual. Ecol. 9, 136–150.
Jackson, J.B., 2008. Ecological extinction and evolution in the brave new ocean. Proc. Thompson, T.L., Glenn, E.P., 1994. Plaster standards to measure water motion. Limnol.
Natl. Acad. Sci. 105 (Supplement 1), 11458–11465. Oceanogr. 39, 1768–1779.
Jud, Z.R., 2014. Anthropogenic Disturbances in Estuarine Ecosystems: The Effects of Thompson, B., Weisberg, S.B., Melwani, A., Lowe, S., Ranasinghe, J.A., Cadien, D.B.,
Altered Freshwater Inflow, Introduction of Invasive Species, and Habitat Alteration in Dauer, D.M., Diaz, R.J., Fields, W., Kellogg, M., Montagne, D.E., 2012. Low levels of
the Loxahatchee River, FL (FIU Electronic Theses and Dissertations). pp. 1197. agreement among experts using best professional judgment to assess benthic condi-
Kanazawa, J., 1975. Distribution of carbaryl and 3, 5-xylyl methlcarbamate in aquatic tion in the San Francisco Estuary and Delta. Ecol. Indic. 12 (1), 167–173.
model ecosystem. J. Agric. Food Chem. 23 (4), 760–763. Van Arman, J., Graves, G.A., Fike, D., 2005. Loxahatchee watershed conceptual ecolo-
Kennish, M.J., 2002. Environmental threats and environmental future of estuaries. gical model. Wetlands 25 (4), 926–942. https://link.springer.com/article/10.1672/
Environ. Conserv. 29 (1), 78–107. 0277-5212(2005)025[0926:LWCEM]2.0.CO;2.
Lao, W., Tsukada, D., Greenstein, D.J., Bay, S.M., Maruya, K.A., 2010. Analysis, occur- Wells, H.W., 1961. The fauna of oyster beds, with special reference to the salinity factor.
rence, and toxic potential of pyrethroids, and fipronil in sediments from an urban Ecol. Monogr. 31 (3), 239–266.
estuary. Environ. Toxicol. Chem. 29, 843–851. Werner, I., Deanovic, L.A., Markiewicz, D., Khamphanh, M., Reece, C.K., Stillway, M.,
Layman, C.A., Jud, Z.R., Archer, S.K., Riera, D., 2014. Provision of ecosystem services by Reece, C., 2010. Monitoring acute and chronic water column toxicity in the Northern
human-made structures in a highly impacted estuary. Environ. Res. Lett. 9 (4), Sacramento–San Joaquin Estuary, California, USA, using the euryhaline amphipod,
044009. Hyalella azteca: 2006 to 2007. Environ. Toxicol. Chem. 29 (10), 2190–2199.
Luisetti, T., Turner, R.K., Jickells, T., Andrews, J., Elliott, M., Schaafsma, M., Watts, W., Whalen, M.A., Duffy, J.E., Grace, J.B., 2013. Temporal shifts in top-down vs. bottom-up
2014. Coastal zone ecosystem services: from science to values and decision making; a control of epiphytic algae in a seagrass ecosystem. Ecology 94 (2), 510–520.
case study. Sci. Total Environ. 493, 682–693. Yeager, L.A., Layman, C.A., 2011. Energy flow to two abundant consumers in a sub-
Metz, J.L., Stoner, E.W., Arrington, D.A., 2015. Comparison of substrates for eastern tropical oyster reef food web. Aquat. Ecol. 45 (2), 267–277.
oyster (Crassostrea virginica) spat settlement in the Loxahatchee River estuary, Zu Ermgassen, P.S., Spalding, M.D., Blake, B., Coen, L.D., Dumbauld, B., Geiger, S.,
Florida. J. Shellfish Res. 34 (3), 861–865. Grabowski, J.H., Grizzle, R., Luckenbach, M., McGraw, K., Rodney, W., 2012.
Parker, M.L., Arnold, W.S., Geiger, S.P., Gorman, P., Leone, E.H., 2013. Impacts of Historical ecology with real numbers: past and present extent and biomass of an
freshwater management activities on eastern oyster (Crassostrea virginica) density and imperiled estuarine habitat. Proc. R. Soc. B 279 (1742), 3393–3400.
recruitment: recovery and long-term stability in seven Florida estuaries. J. Shellfish

508