Harmful Algae 10 (2011) 310–318
Contents lists available at ScienceDirect
Harmful Algae
journal homepage: www.elsevier.com/locate/hal
Stormwater nutrient inputs favor growth of non-native macroalgae (Rhodophyta) on O’ahu, Hawaiian Islands
Brian E. Lapointe *, Bradley J. Bedford
Center for Marine Ecosystem Health, Harbor Branch Oceanographic Institute at Florida Atlantic University, 5600 US 1 North, Ft. Pierce, FL 34946, United States
ARTICLE INFO ABSTRACT
Article history: In Hawaii, blooms of native and non-native macroalgae (limu) have become increasingly problematic in Received 27 August 2010 recent decades. Although the role of human vectors in introducing non-native macroalgae is well Received in revised form 23 November 2010 documented, the ecological role of nutrient pollution in facilitating blooms of these species is not. This Accepted 24 November 2010 study assessed the effects of stormwater discharges on the diversity, abundance, and nutrient content (C, Available online 1 December 2010 N, P and d15N) of native and non-native limu at three sites in the intertidal zone at Ewa Beach, O’ahu. The results showed that native limu species diversity and abundance decreased with proximity to a Keywords: stormwater outfall (ASWO), whereas non-native species abundance increased. Limu tissue d15N values Hawaii at all three sites were within the range reported for sewage N. 15N, %N, and N:P ratios all increased with Limu d Native proximity to the ASWO, supporting the hypothesis that stormwater was a primary source of N Nitrogen enrichment in the study area. In contrast to N, limu %P showed little change among the sites, suggesting Non-native that the generally high N:P ratios, indicative of P-limitation, resulted from high N:P ratios from the Macroalgae upland watershed. Abundance and tissue %N of the non-native rhodophyte Acanthophora spicifera Phosphorus increased with proximity to the ASWO and were strongly correlated (r2 = 0.94) compared to native Stormwater rhodophytes, indicating that stormwater N enrichment provided this invader a competitive advantage (lower C:N ratio) over native limu. These results indicate that the spread of non-native macroalgae in oligotrophic coral reef regions can be facilitated by anthropogenic nutrients in stormwater runoff, thereby threatening native species and ecosystem services. ß 2010 Elsevier B.V. All rights reserved.
1. Introduction aquaculture, has become extremely abundant, especially on Maui where large accumulations accumulate on beaches, interfering Non-native macroalgal invasions are a major driver of coastal with tourists’ use of beaches as a result of malodorous ecosystem change worldwide (UNEP, 2006; Williams and Smith, decomposition (Huisman et al., 2007). 2007). The geographically isolated Hawaiian Islands are especially Throughout the Hawaiian archipelago, there is growing concern vulnerable to biological invasions, which have resulted in about the displacement of native seaweeds, known as limu in the significant impacts on biodiversity (Staples and Cowie, 2001). Hawaiian language, by non-native species. Non-native invasive Biological invasions in the marine environment include introduc- limu compete with and displace native limu species important to tions of non-native seaweeds that date back to the 1950s when the Hawaiians for food, medicine, and religious purposes (Abbott, non-native rhodophyte Acanthophora spicifera was accidentally 1984). Russell (1992) documented how non-native A. spicifera and introduced to Pearl Harbor, O’ahu (Doty, 1961). Today, A. spicifera is H. musciformis displaced native populations of Laurencia nidifica considered the most pervasive non-native algal species in Hawaii and Hypnea cervicornis. In his seminal work on biological invasions, (Huisman et al., 2007). Since the 1950s, over 20 species of non- Elton (1958) emphasized the importance of human-mediated native macroalgae have been introduced to the Hawaiian Islands, vectors, especially physical transport. Humans have subsequently but only about five species have become established and form been recognized as the primary vector in the global epidemic of extensive blooms that alter coastal ecosystems (Russell, 1992; biotic invasions in aquatic ecosystems (Carlton and Geller, 1993). Rodgers and Cox, 1999; Smith et al., 2002). Recently, the non- Indeed, increasing evidence shows that human activities facilitate native rhodophyte Hypnea musciformis, imported from Florida for the physical spread of non-native species in Hawaiian coastal waters, despite recent progress in the prevention of non-native limu introductions (Staples and Cowie, 2001). * Corresponding author. Tel.: +1 772 242 2276; fax: +1 772 468 0757. Increased urbanization of upland watersheds is a major E-mail address: [email protected] (B.E. Lapointe). mechanism increasing nutrient pollution of coastal waters and
1568-9883/$ – see front matter ß 2010 Elsevier B.V. All rights reserved. doi:10.1016/j.hal.2010.11.004 B.E. Lapointe, B.J. Bedford / Harmful Algae 10 (2011) 310–318 311 may facilitate invasions of non-native limu in Hawaiian coastal nitrogen sources on coral reefs, particularly that from human waters. Anthropogenic nutrient pollution of coastal waters has sewage, which can also be a significant source of nitrogen been widely recognized as a common factor linking an array of enrichment in urban stormwater runoff (Wanielista and Yousef, problems, including harmful algal blooms, dead zones, seagrass 1993; Dillon and Chanton, 2008). and coral reef die-offs, declining fisheries, and marine mammal and seabird deaths (ECOHAB, 1997; NRC, 2000; Howarth et al., 2. Materials and methods 2000; MEA, 2005; HARRNESS, 2005; UNEP, 2006). In Hawaiian coastal waters, Soegiarto (1972) and Johannes (1975) first 2.1. Selection of sampling sites suggested a linkage between nutrient pollution from sewage and the expansion of the non-native rhodophyte A. spicifera in To test the hypothesis that stormwater nutrient pollution Kaneohe Bay, O’ahu. Blooms of the native invasive chlorophyte affects the relative abundance of non-native versus native limu in Dictyosphaeria cavernosa, which overgrew and killed corals in intertidal communities, three study sites were chosen along a Kaneohe Bay, O’ahu, were also linked to nutrient enrichment from gradient of exposure to stormwater discharge in the Ewa Beach sewage (Banner, 1974; Smith et al., 1981). Following sewage area on O’ahu (Fig. 1). The urbanized Ewa Beach area had a diversion from Kaneohe Bay in the late 1970s, nutrient concentra- stormwater drainage system constructed when this neighborhood tions and D. cavernosa biomass both decreased (Hunter and Evans, was developed in the 1970s. Currently, several stormwater outfalls 1995), demonstrating the ecological importance of nutrient discharge into the intertidal zone where longshore currents in the enrichment to Kaneohe Bay. Since then, however, blooms of nearshore area generally flow from the urbanized Ewa Beach area non-native limu have expanded throughout the Hawaiian Islands westwardly towards One’Ula State Beach Park. The easternmost (Russell, 1992; Rodgers and Cox, 1999; Smith et al., 2002), sampling site (Amio) in our study, located near the Amio Street especially in coastal waters adjacent to urbanized watersheds. stormwater outfall (ASWO), was chosen to be representative of Although sewage has been identified as a significant nitrogen direct stormwater impacts (Fig. 2a). A second site near Papipi Road source supporting blooms of native and non-native limu (Dailer (Papipi), 625 m west of Amio, was chosen as a site less impacted et al., 2010), studies have not addressed the importance of nutrient by the ASWO discharges to the east (Fig. 2b). The third site at Kaloi enrichment (nitrogen, N and phosphorus, P) from urban storm- Gulch (Kaloi) in One’Ula State Beach Park, west of Papipi and water runoff to the relative abundance of native and non-native 975 m from Amio, was selected as a reference site least impacted limu, such as A. spicifera, in Hawaiian coastal waters. Urban by stormwater discharges to the east (Fig. 2c). The three intertidal stormwater runoff can contain relatively high concentrations of sites were sampled during low (minus) tides, March 3–7, 2008. ammonium, nitrate, total N, soluble reactive P, and total P, which combined with the high volumes of stormwater following rain 2.2. Sampling for taxonomic composition of limu communities events, can account for considerable nutrient loads to coastal waters (Wanielista and Yousef, 1993). Four separate 30 m survey transects were established end-to- We posed the hypothesis that nutrients from stormwater runoff end in the low intertidal zone at each site using Keson fiberglass would affect the nutrition and relative abundance of native and survey tapes. The four transects were independent, had no spatial non-native limu. To test this hypothesis, we studied intertidal limu overlap, and resulted in a total surveyed length of 120 m at each communities at three locations in Ewa Beach, O’ahu. The study was site. Qualitative collections of conspicuous limu species were multi-faceted, and involved measuring the following variables in sampled along the four transects at each site; specimens were intertidal communities at the three sites: macroalgal species identified according to Abbott (1999), Abbott and Huisman (2004), presence, percent cover of abundant taxa, tissue C:N:P contents to and Huisman et al. (2007). gauge the degree of N versus P limitation (Atkinson and Smith, Limu communities along the four transects at each site were 1983; Lapointe et al., 1992), and stable nitrogen isotope ratios quantitatively sampled by photogrammetric techniques, using a (d15N) to discriminate between anthropogenic and natural Nikon Coolpix 5000 camera to obtain high resolution digital color nitrogen sources such as sewage, fertilizers, upwelling, or nitrogen images along the transects. This non-destructive method yields fixation (Lapointe, 1997; Lapointe et al., 2005; Derse et al., 2007; parallax-free sampling that generates highly reproducible quanti- Parsons et al., 2008; Dailer et al., 2010). Risk et al. (2009) concluded tative data, and is one of the most widespread and sophisticated that the measurement of stable nitrogen isotopes in macroalgae techniques for permanently recording marine algal standing stocks 2 [()TD$FIG]provides a cost-effective and objective means of quantifying (Littler and Littler, 1985). Ten digital images (0.1 m quadrat size)
Fig. 1. Satellite image showing the Ewa Beach study sites. 312[()TD$FIG] B.E. Lapointe, B.J. Bedford / Harmful Algae 10 (2011) 310–318
Fig. 2. Photographs showing: (a) Amio stormwater outfall pipe (ASWO); (b) intertidal zone at Papipi; (c) transect sampling tape at Kaloi Gulch; (d) red limu at Amio showing dark pigmentation indicative of nitrogen enrichment; (e) Sargassum echinocarpum at Papipi; and (f) mixed native limu community at Kaloi Gulch showing reduced pigmentation indicative of nitrogen limitation. were recorded, at 3 m intervals along each of the four 30 m transport to the lab. At the ASWO pipe, replicate (n = 2 per species) transects, yielding 40 photo-quadrats per site. At some locations samples of two limu species, S. echinocarpum and L. majuscula, along the transects, lack of suitable substrate prevented sampling. were also collected, to establish a stormwater N isotope source At these locations, the closest suitable substrate was sampled, signal. The freshly collected limu were processed in Dr. Robert typically within 1 m of the original location. Images were Richmond’s laboratory at the Kewalo Marine Lab, Honolulu, where analyzed using the Randomized Point Count method (Littler and the specimens were cleaned of epibionts, rinsed in deionized water Littler, 1985). Ten random points were overlayed on each quadrat and dried in a lab oven at 65 8C for 48 h. Dried limu samples were image (displayed on a high resolution LCD monitor) and the limu or ground to powder, using a mortar and pestle, and stored in plastic other biota/substrate beneath each point identified. In the case of screw-cap vials for later analysis (Lapointe et al., 2005). some limu that were difficult to identify to species (such as Limu d15N ratios were measured to aid in identification of the N Laurencia spp.) identification was made only to genus. Small source(s) in stormwater sustaining the limu communities (Heaton, microfilamentous turfs <2 cm high were scored as either ‘‘algal 1986; Dawson et al., 2002). Previous studies of d15N ratios in turf’’ (greens) or ‘‘red turf’’ depending on pigmentation. macroalgae experiencing variable inputs of stormwater runoff in Sarasota Bay, FL, reported values ranging from +1.3% to +13.3% 2.3. Sampling limu for nutrient contents and stable nitrogen (Dillon and Chanton, 2008). Analyses of replicate limu tissue isotope ratios samples (n = 4 per sample) were performed at the University of Maryland Center for Environmental Science, Horn Point Laborato- To quantitatively compare nutrient availability among the ry, Cambridge, MD. Samples were packed in tin capsules and three sites, four limu species common to all three sites were analyzed using a Sercon isotope ratio mass spectrometer. The collected from each site for tissue C:N:P contents and d15N standard used for stable nitrogen isotope analysis was d15N in air. 15 analysis. These included two phaeophytes, Padina sanctae-crucis d N values (%) were calculated as [(Rsample/Rstandard) À 1] Â 100; and Sargassum echinocarpum, and two rhodophytes, Laurencia where R = 15N/14N. majuscula and A. spicifera. Replicate (n = 2 per species), composite Limu samples were also analyzed for percent C, N, and P, at samples (10–20 thalli per replicate) were collected at each site, Nutrient Analytical Services, Chesapeake Biological Laboratory, placed in Whirl-Pak baggies, and held over ice in a cooler for University of Maryland System, Solomons, MD. C and N were B.E. Lapointe, B.J. Bedford / Harmful Algae 10 (2011) 310–318 313 measured on an Exeter Analytical, Inc., CE-440 Elemental Analyzer; Table 1 P was measured following the methodology of Asplia et al. (1976) Taxonomic inventory at each site, including site totals. on a Technicon Autoanalyzer II (Keefe et al., 2004). Kaloi Papipi Amio
Phaeophyta 2.4. Statistical analyses Asteronema breviarticulatum X Colpomenia sinuosu XX X Macroalgal benthic cover (%), C:N:P ratios, and d15N(%) value Dictyota sandvicensis XX X Dictyota acutiloba XX comparisons among the three study sites (df = 2) and between Dictyopteris plagiogramma X Phyla (df = 1) were performed in SPSS 11 for Mac (www.spss.com) Hincksia indica X using the Generalized Linear Model procedure (GLM; Type III sum- Hydroclathrus clathratus X of-squares). When values among sites were found by GLM to differ Padina sanctae-crucis XX X significantly (p 0.05), Tukey’s HSD (THSD) post hoc testing was Sargassum echinocarpum XX X Sargassum obustifolium X performed to identify the source(s) of difference(s). Correlative Sphaceleria novae-hollandiae XX comparison of A. spicifera tissue N (%) with A. spicifera benthic cover Stypopodium flabelliforme X (%) was performed using least-squares linear regression in Chlorophyta Microsoft Excel. Avrainvillea amadelpha X Bornetella sphaerica X Caulerpa microphysa XX 3. Results Caulerpa antonensis X Caulerpa taxifolia X 3.1. Taxonomic composition of intertidal limu communities Chaetomorpha antennina XX Cladophora catenata X Codium edule XX X A total of 52 limu species were collected and identified from the Codium arabicum XX three study sites. The Kaloi site, farthest from stormwater Codium reediae X discharges, had 39 spp., the highest species richness among sites; Dictyosphaeria versluysii XX X Papipi, intermediate in distance from the stormwater discharges, Halimeda opuntia XX X had 30 species; Amio, closest to the stormwater inputs, had 25 Halimeda discoidea XX Microdictyon setchellianum X species (Table 1 and Fig. 3). Neomeris vanbosseae XX Benthic cover at Kaloi was dominated by native Laurencia spp. Ulva fasciata XX (34.5 Æ 5.4%) and P. sanctae-crucis (33.8 Æ 7.1%) with smaller Rhodophyta amounts of algal turf (11.3 Æ 2.9%), S. echinocarpum (7.8 Æ 6.6%), Acanthophora specifera XX X Asparagopsis taxiformis XX X and Dictyota spp. (5.8 Æ 3.5%), and trace amounts of the non-native Botryocladia skottsbergii XX invasive rhodophyte A. spicifera (0.25 Æ 0.5%) and other taxa (Figs. 2f Gracilaria coronopifolia XX X and 4). At Papipi, S. echinocarpum (see Fig. 2e; 20.3 Æ 17.7%), Laurencia Grateloupia phuquocensis XX spp. (13 Æ 9.8%) and Pterocladiella spp. (12.5 Æ 13.9%) were the most Griffithsia heteromorpha X abundant taxa. However, the native invasive chlorophyte Ulva fasciata Hypnea spinella X Hypnea chordaceae XX (8.8 Æ 16.2%) and the non-native invasive rhodophytes A. spicifera Hypnea musciformis XX (7.5 Æ 5.1%) and H. musciformis (2.5 Æ 2.6%) also contributed signifi- Jania micrarthrodia X cantly to overall cover. At Amio, A. spicifera was the single most Jania micrarthrodia XX abundant limu species (29.8 Æ 11.5%), followed by Laurencia spp. Liagora sp. X Laurencia spp. X (25 Æ 12%), S. echinocarpum (18.8 Æ 16%), H. musciformis (8.3 Æ 1.0%) Laurencia dotyi X and other taxa (Fig. 2d). Laurencia majuscula XX X Total benthic cover did not vary significantly among sites, Laurencia mcdermidiae XX X averaging 95.9 Æ 2.4% overall (Fig. 3). However, cover of the non- Plocamium sandvicense XX native rhodophytes A. spicifera and H. musciformis was significantly Portieria hornemanni X Pterocladiella capillaceae XX higher at Amio than at Kaloi and Papipi (p 0.005, THSD; Fig. 4). Per- Pterocladiella caerlescens XX species average non-native cover (Fig. 3) also was significantly higher Spyridia filamentosa XX X at Amio (19.0 Æ 13.8%) than at Kaloi (0.1 Æ 0.4%) and Papipi Tricholgloea requienii XX (5.0 Æ 4.6%; p 0.008, THSD). In contrast, average native species Wrangelia requienni X Total taxa 39 30 25 cover was significantly lower at Amio (3.4 Æ 8.3%) than at Kaloi (8.5 Æ 13.6%; p = 0.030, THSD). These decreases in native species cover indicate a loss of native species abundance, with a correspond- spicifera) and Phaeophyta (brown limu: P. sanctae-crucis and S. ing increase in non-native invasive species abundance with proximity echinocarpum) for all variables except C:P ratios. Therefore, to the stormwater discharges (Fig. 3). This pattern was maintained in comparisons among sites were made separately for the two Phyla the summed benthic cover data (Fig. 3), which showed that total non- – Rhodophyta and Phaeophyta. Overall and individual species native benthic cover was significantly higher at Amio (36.3 Æ 13.2%) means (Æ1 SD) for all measured variables are shown in Table 2. In one than at Kaloi (0.3 Æ 0.5%) and Papipi (10.0 Æ 5.0%; p 0.004, THSD); or both Phyla, significant trends in d15N values, %N, N:P ratios, %P conversely, total native cover was significantly lower at Amio (rhodophytes only) and C:N ratios (rhodophytes only) between Kaloi (59.8 Æ 12.5%) than at Kaloi (96.3 Æ 1.5%) and Papipi (85.3 Æ 5.6%; and the AWSO were observed (Fig. 5). However, no significant p 0.004, THSD). differences among sites were observed for %C or C:P ratios, either overall or by Phylum (Fig. 5). d15N values in both rhodophytes and 3.2. Nutrient composition of intertidal limu communities phaeophytes increased significantly from Kaloi to the ASWO (p 0.007, THSD; as indicated in Fig. 5). Rhodophyte d15N averaged Comparison of pooled species among sites can become +11.5 Æ 2.0%, and ranged from +10.4 Æ 1.2% at Kaloi and +10.0 Æ confounded when tissue nutrient composition differs significantly 1.3% at Papipi to +12.7 Æ 1.6% at Amio and +14.2 Æ 0.6% at the among taxa. In this study, tissue nutrient composition differed ASWO. Phaeophyte d15N averaged +11.0 Æ 1.3% and ranged from significantly between Rhodophyta (red limu: L. majuscula and A. +9.8 Æ 1.1% at Kaloi to +11.4 Æ 1.3% at the ASWO. 314[()TD$FIG] B.E. Lapointe, B.J. Bedford / Harmful Algae 10 (2011) 310–318
Tissue C (% organic) in both rhodophytes and phaeophytes was similar among sites, as were tissue C:P ratios (Fig. 5). Rhodophyte C averaged 28 Æ 3%, ranging from 27 Æ 2% at Amio to 29 Æ 2% at Papipi; C:P averaged 1046 Æ 306% and ranged from 1218 Æ 506% at Kaloi to 717 Æ 3% at the ASWO. Phaeophyte C averaged 25 Æ 4% and ranged from 23 Æ 5% at Amio to 29 Æ 1% at the ASWO; C:P averaged 1199 Æ 337% and ranged from 1079 Æ 273% at Kaloi to 1641 Æ 18% at the AWSO. Tissue N (%) of both rhodophytes and phaeophytes increased from Kaloi to the ASWO (p 0.026; THSD), but was significantly higher (p = 0.012, GLM) and had a wider range in the rhodophytes than the phaeophytes (Fig. 5). Rhodophyte N averaged 2.0 Æ 0.6% and ranged from 1.5 Æ 0.2% at Kaloi to 3.5 Æ 0.1% at the AWSO. This increase in tissue N was most evident in the invasive, non-native A. spicifera, which increased from 1.5 Æ 0.3% at Kaloi Gulch to 3.5 Æ 0.1% at the ASWO (Table 2) and correlated strongly with benthic cover (r2 = 0.94; Fig. 6). Phaeophyte N averaged 1.2 Æ 0.4% and ranged from 0.8 Æ 0.1% at Kaloi to 1.7 Æ 0.2% at the ASWO. Tissue C:N ratios decreased from Kaloi to the ASWO, but the decrease was significant only for the rhodophytes (p 0.013), which averaged 16 Æ 5 and ranged from 22 Æ 5 at Kaloi to 10 Æ 0.2 at the AWSO (Fig. 5). Phaeophyte C:N averaged 27 Æ 10 and ranged from 35 Æ 13 at Kaloi to 20 Æ 2 at the AWSO. Tissue P (%) showed no consistent trends among sites, and only rhodophyte P differed significantly among sites (Fig. 5) with those at the ASWO having the highest P (0.10 Æ 0.00%; p 0.040, THSD). Rhodophyte P averaged 0.07 Æ 0.02% and was the lowest at Kaloi (0.06 Æ 0.01%), but not significantly lower than at Amio or Papipi. Phaeophyte P averaged 0.05 Æ 0.01% and ranged from 0.04 Æ 0.00% at the ASWO to 0.06 Æ 0.01% at Papipi. Fig. 3. (a) Mean (per species) native and non-native benthic limu cover (%) and Tissue N:P ratios also increased, in both rhodophytes (p = 0.037) numbers of taxa; (b) total benthic cover (%) of native and non-native limu at the and phaeophytes (p = 0.038), between Kaloi and the ASWO (Fig. 5). three study sites. Values are means Æ SE for percent cover (n = 4). Values with matching symbols differ significantly. Rhodophyte N:P averaged 70 Æ 19 and ranged from 54 Æ 12 at Kaloi to 88 Æ 22 at Amio. This increase in N:P ratio was again most evident in A. spicifera, which ranged from 45 Æ 6 at Kaloi to 107 Æ 1 at Amio (Table 2). Phaeophyte N:P ratios averaged 49 Æ 17 and ranged from [()TD$FIG] 32 Æ 4 at Kaloi to 82 Æ 8 at the ASWO.
4. Discussion
Multiple lines of evidence from our study support the hypothesis that nutrients from urban stormwater runoff support non-native limu invasions in the Ewa Beach area. First, as the non- native benthic cover increased with proximity to the stormwater outfall, native species diversity and abundance decreased. The two non-native rhodophytes, A. spicifera and H. musciformis, both increased in abundance with proximity to stormwater discharges at AWSO. Second, the abundance of A. spicifera, the most pervasive alien invader in Hawaii’s coastal waters (Russell, 1992; Huisman et al., 2007), was strongly correlated with tissue %N (r2 = 0.94), suggesting stormwater N enhances growth and abundance of this species in N-limited coastal waters of O’ahu (cf. Larned, 1998). Third, mean N:P ratios increased significantly in both red and brown limu from Kaloi to Amio, which further indicated increased N availability from ASWO discharges. Fourth, d15N values increased significantly between Kaloi and AWSO, pointing to stormwater discharges as the source of 15N enrichment. The highest tissue %N, %P, N:P ratio, and d15N values all occurred in limu collected from the ASWO, which provides compelling evidence that stormwater discharges were the primary source of both N and P enrichment in our study area. The significant increase in N:P ratio between the Kaloi and Amio sites indicates that the stormwater discharges have a relatively high N:P ratio. N:P ratios reported for stormwater on O’ahu are Fig. 4. Benthic cover (%) of the 10 most common taxa at each of the three study sites. typically low, ranging from 7.3 (Presley and Jamison, 2009)to16 Values are means Æ SE (n = 4). (USGS, 2010), and reflect relatively high P inputs from human and/ B.E. Lapointe, B.J. Bedford / Harmful Algae 10 (2011) 310–318 315
Table 2 Tissue C, N, and P (%); C:N, C:P, and N:P molar ratios; and d15N values (%) of limu collected at the study sites, including pooled values. Values are means Æ 1 SD.
n %C %N %P C:N C:P N:P d15N(%)
Kaloi Acanthophora spicifera 2 23.1 Æ 1.3 1.5 Æ 0.3 0.07 Æ 0.01 18 Æ 3 811 Æ 24 45 Æ 6 9.8 Æ 1.3 Laurencia majuscula 2 31.4 Æ 3.8 1.4 Æ 0.1 0.05 Æ 0.00 26 Æ 1 1625 Æ 326 63 Æ 9 10.9 Æ 0.8 Sargassum echinocarpum 2 29.0 Æ 0.7 0.8 Æ 0.1 0.06 Æ 0.02 46 Æ 8 1306 Æ 127 29 Æ 2 10.4 Æ 1.2 Padina sanctae-crucis 2 18.9 Æ 1.8 0.9 Æ 0.1 0.06 Æ 0.00 24 Æ 0 853 Æ 53 35 Æ 2 9.1 Æ 0.5 Kaloi mean 8 25.6 Æ 5.5 1.1 Æ 0.4 0.06 Æ 0.00 24 Æ 0 853 Æ 53 43 Æ 14 10.1 Æ 1.1 Papipi Acanthophora spicifera 2270Æ 0.6 2.1 Æ 0.2 0.08 Æ 0.00 15 Æ 1 906 Æ 50 62 Æ 8 10.2 Æ 1.1 Laurencia majuscula 2 31.1 Æ 0.1 2.2 Æ 0.2 0.08 Æ 0.01 17 Æ 2 1068 Æ 161 65 Æ 16 9.9 Æ 1.7 Sargassum echinocarpum 2 25.1 Æ 0.9 1.1 Æ 0.1 0.05 Æ 0.00 27 Æ 0 1218 Æ 30 45 Æ 2 12.6 Æ 1.0 Padina sanctae-crucis 2 25.0 Æ 1.8 1.7 Æ 0.1 0.07 Æ 0.01 17 Æ 0 956 Æ 41 57 Æ 2 11.2 Æ 1.1 Papipi mean 8 27.0 Æ 2.8 1.8 Æ 0.5 0.07 Æ 0.01 19 Æ 5 1037 Æ 145 57 Æ 11 11.0 Æ 1.6 Amio Acanthophora spicifera 2 25.1 Æ 1.0 2.9 Æ 0.5 0.06 Æ 0.01 10 Æ 2 1110 Æ 171 107 Æ 1 11.7 Æ 1.2 Laurencia majuscula 2 28.2 Æ 0.2 2.1 Æ 0.2 0.07 Æ 0.00 16 Æ 1 1087 Æ 58 69 Æ 1 13.8 Æ 1.1 Sargassum echinocarpum 2 27.6 Æ 0.4 1.1 Æ 0.1 0.05 Æ 0.01 31 Æ 3 1601 Æ 256 52 Æ 4 11.3 Æ 0.4 Padina sanctae-crucis 2 18.4 Æ 1.4 1.1 Æ 0.1 0.06 Æ 0.01 20 Æ 4 821 Æ 110 41 Æ 2 11.6 Æ 0.6 Amio mean 8 24.8 Æ 4.2 1.8 Æ 0.8 0.06 Æ 0.01 19 Æ 8 115 Æ 326 67 Æ 27 12.1 Æ 1.3 ASWO Laurencia majuscula 2 28.4 Æ 0.6 3.5 Æ 0.0 0.10 Æ 0.00 10 Æ 0 717 Æ 3 76 Æ 2 14.2 Æ 0.6 Sargassum echinocarpum 2 28.6 Æ 0.6 1.7 Æ 0.2 0.05 Æ 0.00 20 Æ 2 1641 Æ 18 82 Æ 8 15.1 Æ 0.4 ASWO mean 4 28.8 Æ 0.5 2.6 Æ 1.1 0.07 Æ 0.03 15 Æ 6 1179 Æ 534 79 Æ 6 14.7 Æ 0.7 Grand mean 28 26.2 Æ 4.0 1.7 Æ 0.8 0.06 Æ 0.02 21 Æ 9 1123 Æ 326 59 Æ 21 11.7 Æ 1.9 or animal waste, fertilizers, or other sources. However, high N:P Hawaiian Islands in recent decades (Smith et al. 2002; Huisman ratios ranging from 39 to 247 have been reported for human- et al., 2007). These blooms are generally located in bays (Kaneohe impacted streams in Waimana, O’ahu; the relatively low P Bay) and coastlines (Waikiki, Lahaina) impacted by urbanization concentrations in these streams probably reflected the high iron and increased stormwater and/or wastewater nutrient loads from content of Hawaiian soils, which effectively immobilizes P in human activities (Dailer et al., 2010). In southwest Florida, massive groundwater (Laws and Ferentinos, 2003). Similarly high limu N:P blooms of red drift macroalgae, including Gracilaria spp. and ratios were associated with stormwater discharges in our study, as Hypnea spp., have emerged in coastal waters that are enriched by N:P ratios averaged 43 at Kaloi, and increased to 79 at the ASWO, pulsed discharges of both N and P from the Caloosahatchee and indicating a shift towards stronger P-limitation along this gradient Peace rivers, respectively (Lapointe and Bedford, 2007). of stormwater enrichment (Atkinson and Smith, 1983; Lapointe Among all the native and non-native limu species assessed in et al., 1992). Elsewhere on O’ahu, in Kaneohe Bay, the mean N:P this study, A. spicifera responded most to increasing nutrient ratio of 30 limu species was 44 (Atkinson and Smith, 1983; Smith, enrichment from the stormwater inputs. Between the Kaloi and 1994), a value remarkably similar to the mean value of 43 at Kaloi Amio sites, tissue %N doubled (1.5–2.9%) while C:N ratio decreased in this study. Elevated limu N:P ratios at Amio (67) and the ASWO (18–10) to non-limiting levels in A. spicifera; in comparison, the (79) were consistent with those in other human-impacted, P- native L. majuscula had relatively small increases in %N (1.4–2.1%) limited, carbonate-rich marine environments. For example, high and its C:N ratio was relatively high at both sites (26–16), N:P ratios (>100) occur in groundwaters in the Florida Keys suggesting nitrogen limitation (Lapointe and Bedford, 2010). On impacted by septic tanks (Lapointe et al., 1990). the Belizean Barrier Reef, experimental N pulses significantly Nutrient availability is a major factor affecting competition increased %N and photosynthesis of A. spicifera, a native species in among limu in tropical oligotrophic settings but has been poorly this location (Lapointe et al., 1987). The %N of the Belizean A. documented for introduced invasive macroalgae (Williams and spicifera ranged from 0.48% to 0.69% dry weight, 3-fold lower Smith, 2007). Although Soegiarto (1972) and Johannes (1975) than the 1.5–2.9% range found in A. spicifera from Ewa Beach (Table noted an apparent link between sewage pollution and the 2). These N-enriched conditions in the Ewa Beach study area, expansion of A. spicifera in Kaneohe Bay four decades ago, studies especially at Amio, support increased productivity and cover of the examining the relationship(s) among nutrient source(s), tissue non-native A. spicifera at the expense of native limu (Russell, 1992). C:N:P ratios and abundance of this pervasive invader have not been The superior ability of A. spicifera to assimilate nitrogen allows this conducted. Rhodophytes, such as A. spicifera and Hypnea spp., are invader to compete favorably with native species in eutrophic favored in nutrient-enriched tropical coastal waters where they environments and explains why this species is such a pervasive can outcompete other species that dominate under more invader in Hawaiian waters experiencing land-based nutrient oligotrophic conditions (Lapointe et al., 2004). The significantly enrichment. Similarly, on coral reefs off highly populated higher %N, %P, and N:P ratios in the rhodophytes compared to the southeast Florida, the non-native Caulerpa brachypus had a phaeophytes in this study, along a gradient of nutrient enrichment, significantly lower C:N ratio than native species, enabling this illustrate the competitive advantage of fast-growing rhodophytes alga to overgrow and displace native species in these nutrient- that sequester growth-limiting nutrients when they become enriched habitats (Lapointe and Bedford, 2010). available during episodic stormwater discharge events. Our In urban areas that rely on cesspools and septic tanks for human observations are consistent with early aquaculture studies that sewage disposal, stormwater runoff can be enriched with human found highly productive rhodophytes are capable of rapid uptake waste through the process of infiltration and inflow of contami- and storage of N (D’Elia and DeBoer, 1978). This phenomenon is not nated groundwater into the stormwater collection system restricted to rhodophytes in the Ewa Beach area, but is also (Wanielista and Yousef, 1993). Concentrations of total nitrogen evidenced by blooms of non-native rhodophytes – A. spicifera, H. in stormwater on O’ahu have ranged from 40 mM(USGS, 2010)to musciformis, Gracilaria salicornia and Kappaphycus alvarezii – that 70 mM(Presley and Jamison, 2009). Although these concentrations have become common around urbanized coastlines of the are lower than that typical of sewage effluent, the high volume 316[()TD$FIG] B.E. Lapointe, B.J. Bedford / Harmful[()TD$FIG] Algae 10 (2011) 310–318
Fig. 6. Linear regression showing the relationship between Acanthophora spicifera tissue N contents (%; n = 2) and benthic cover (%; n = 4). Values are means Æ SE.
similar to wetland plant species on O’ahu, which average +8.2 Æ 1.2% (n = 11; Table 3); d15N values decrease in wetland species on other less populated Hawaiian Islands where lower levels of urbanization and d15N values occur (Bruland and MacKenzie, 2010). Laws et al. (1999) reported values of +6% for particulate nitrogen in downstream coastal waters off Ewa Beach, which they attributed to wastewater (cesspools) enrichment of groundwaters on the Ewa plain. Other studies in the Hawaiian Islands have reported variable d15N values in coastal macroalgae, depending on the type and degree of anthropogenic influence: highly elevated values of +17.8 to 50.1% off Maui (Lahaina, Kihei, Kahului) were indicative of wastewater effluent (Dailer et al., 2010) compared to lower values of À0.5% off Kauai that reflect fertilizer N as well as relatively low inputs of sewage N (Derse et al., 2007). Similar patterns have been reported globally. On coral reefs in the vicinity of sewage outfalls off highly populated southeast Florida, invasive blooms of native and non-native chlorophytes had elevated d15N values ranging from +8 to +10% (Lapointe et al., 2005). In the Buccoo Reef Complex in Tobago, West Indies, macroalgal blooms adjacent to a sewage outfall had d15N values ranging from +10 to +12% (Lapointe et al., 2010). All these studies support the recent conclusions of Risk et al. (2009) that measurement of stable nitrogen isotopes in macroalgae provides a cost-effective and objective means of quantifying sewage stress on coral reefs.
5. Summary and conclusions
Fig. 5. Tissue chemistry (d15N, %C, %N, %P, C:N, C:P, N:P) of red (Rhodophyta) and Our results showed that native limu species diversity and brown (Phaeophyta) limu at three study sites and the Amio stormwater outfall pipe abundance decreased with proximity to stormwater discharges, (ASWO). Values are means Æ SE (n = 4 for Kaloi, Papipi, and Amio; n = 2 for ASWO). whereas non-native limu increased. The highest tissue %N, %P, N:P Dumbbells (dots connected by a line) connect statistically similar values; unconnected ratio and d15N values all occurred in limu collected from the ASWO, columns differ significantly. which provides compelling evidence that cumulative impacts from episodic stormwater discharges were the primary source of associated with stormwater runoff during episodic events can nutrient enrichment in the study area. The abundance of the result in high overall N loadings to coastal waters. This mechanism non-native invasive A. spicifera was strongly correlated with tissue could explain why the d15N values that we observed in exposed %N (r2 = 0.94), suggesting that stormwater N conferred a competi- macroalgae were elevated to values typical of sewage N (Heaton,
1986). In general, limu that rely on natural N from nitrogen fixation Table 3 15 have d N values close to atmospheric N (0%; France et al., 1998) Tissue d15N values (%) of coastal wetland plants in the Hawaiian and become progressively enriched as a result of assimilation of Islands from Bruland and MacKenzie (2010). d15N values are fertilizer (0 to +3%) and wastewater (+3 to +15%) N sources means Æ S.E.; ‘‘n’’ values are the number of sites sampled per island. (Heaton, 1986; Lapointe, 1997; Costanzo et al., 2001; Lapointe Island d15N(%) n et al., 2004, 2005; Lapointe and Bedford, 2007). The significant 15 Oahu 8.2 Æ 1.2 11 increase in d N values from +10 to +15% between Kaloi and Maui 7.0 Æ 1.0 6 Amio supports the hypothesis that stormwater discharges from Molokai 4.4 Æ 0.1 2 ASWO were a primary source of d15N enrichment to limu in the Kauai 3.4 Æ 0.1 8 study area. The high background d15N value of +10% in this study Hawaii 3.0 Æ 0.2 7 Mean 5.2 Æ 0.2 15 is typical for limu growing on wastewater N (Risk et al., 2009) and B.E. Lapointe, B.J. Bedford / Harmful Algae 10 (2011) 310–318 317 tive advantage on this species compared with native species, Costanzo, S.D., O’Donohue, M.J., Dennison, W.C., Loneragan, N.R., Thomas, M., 2001. A new approach for detecting and mapping sewage impacts. Mar. Poll. Bull. 42, particularly phaeophytes. N:P ratios in limu were relatively high in 149–156. the study area, indicating P-limitation, and increased in proximity D’Elia, C.F., DeBoer, J.A., 1978. Nutritional studies of two red algae. 2. Kinetics of to ASWO discharges. Finally, d15N values increased with increasing ammonium and nitrate uptake. J. Phycol. 14, 266–272. 15 proximity to the ASWO and limu at all sites were in the range of Dailer, M.L., Knox, R.S., Smith, J.E., Napier, M., Smith, C.M., 2010. Using d N values in algal tissue to map locations and potential sources of anthropogenic sewage N, suggesting that human and/or animal waste from nutrient inputs on the island of Maui, Hawaii, USA. Mar. Poll. Bull. 60, 655– stormwater runoff influenced the entire study area. 671. Blooms of both native and non-native limu have been Dawson, T.E., Mambelli, S., Plamboeck, A.H., Templer, P.H., Tu, K.P., 2002. Stable isotopes in plant ecology. Annu. Rev. Ecol. Ecol. Syst. 33, 507–559. increasingly problematic in the Hawaiian Islands over the past Derse, E., Knee, K.L., Wankel, S.D., Kendall, C., Berg, C.J., Paytan, A., 2007. Identifying five decades (Doty, 1961; Johannes, 1975; Smith et al., 1981, 2002; sources of nitrogen to Hanalei Bay, Kauai, using the nitrogen isotope signature Russell, 1987, 1992; Rodgers and Cox, 1999; Huisman et al., 2007). of macroalgae. Environ. Sci. Technol. 41, 5217–5223. Dillon, K.S., Chanton, J.P., 2008. Nitrogen stable isotopes of macrophytes assess The results of this study indicate that the spread of non-native stormwater nitrogen inputs to an urbanized estuary. Estuaries Coasts 31, 360– rhodophytes in Ewa Beach is facilitated by nutrient pollution, in 370. this case from stormwater discharges. However, nutrient enrich- Doty, M.S., 1961. Acanthophora, a possible invader of the marine flora of Hawaii. Pac. Sci. 15, 547–552. ment via submarine groundwater discharges can also deliver ECOHAB, 1997. The ecology and oceanography of harmful algae blooms. In: anthropogenic nutrients that support invasive blooms of limu, Anderson, D.M. (Ed.), A National Research Agenda. Woods Hole Oceanographic such as on Maui where injection wells are used to dispose of Institution, Woods Hole. treated sewage (Dailer et al., 2010). Sea turtles forage on the Elton, C.S., 1958. The Ecology of Invasions by Animals and Plants. Methuen and Co., London. expanding blooms of non-native invasive limu and records of France, R., Holmquist, J., Chandler, M., Cattaneo, A., 1998. Evidence for nitrogen fibropapilloma indicate a cluster of the disease in Hawaiian coastal fixation associated with macroalgae from a seagrass–mangrove–coral reef waters impacted by human nitrogen enrichment (Van Houtan system. Mar. Ecol. Prog. Ser. 167, 297–299. HARRNESS, 2005. Harmful Algal Research and Response: A National Environmental et al., 2010). With an increase in non-native invasive limu, the Science Strategy 2005–2015. Woods Hole Oceanographic Institution, Woods corresponding decrease in native limu abundance and diversity Hole. will also diminish reef ecosystem services (cf. Ruddle and Heaton, T.H.E., 1986. Isotopic studies of nitrogen pollution in the hydrosphere and atmosphere: a review. Chem. Geol. 59, 87–102. Johannes, 1985) and threaten traditional cultural uses of native Howarth, R., Anderson, D., Cloern, J., Elfring, C., Hopkinson, C., Lapointe, B., Malone, limu species. T., Marcus, N., McGlathery, K., Sharpley, A., Walker, D., 2000. Nutrient pollution Globally, non-native macroalgal invasions are expected to of coastal rivers, bays, and seas. Ecology 7, 1–15 (issues). Huisman, J.M., Abbott, I.A., Smith, C.M., 2007. Hawaiian Reef Plants. Report # UNIHI- increase in coming years due to the anticipated increase in SEAGRANT-BA-03-02. University of Hawaii Sea Grant College Program, Hono- unintended introductions of non-native organisms as a side effect lulu. of globalization (UNEP, 2006). Land-based nutrient pollution has Hunter, C.L., Evans, C.W., 1995. Coral reefs in Kaneohe Bay, Hawaii: two centuries been identified as an important factor facilitating invasions of non- of western influence and two decades of data. Bull. Mar. Sci. 57 (1), 501– 515. native Caulerpa in southeast Florida (Lapointe et al., 2005, 2010) Jaubert, J.M., Chisholm, J.R.M., Minghelli-Roman, A., Marchioretti, M., Morrow, J.H., and the Mediterranean Sea (Chisholm et al., 1997; Jaubert et al., Ripley, H.T., 2003. Re-evaluation of the extent of Caulerpa taxifolia development 2003). Hence, the role of land-based nutrient pollution in Hawaii, in the northern Mediterranean using airborne spectrographic sensing. Mar. Ecol. Prog. Ser. 263, 75–82. and elsewhere, must be recognized if coastal ecosystems are to be Johannes, R.E., 1975. Pollution and degradation of coral reef communities. In: Wood, managed effectively in the future. E.J.F., Johannes, R.E. (Eds.), Tropical Marine Pollution. Elsevier Scientific Pub- lishing Company, Oxford, pp. 13–51. Keefe, C.W., Blodnikar, K.L., Boynton, W.R., Clark, C.A., Frank, J.M., Kaumeyer, N.L., Acknowledgements Weir, M.M., Wood, K.V., Zimmermann, C.F., 2004. Nutrient Analytical Services Laboratory Standard Operating Procedures. Special Publication Series SS-80-04- We would like to acknowledge the assistance of Dr. Phil CBL. Chesapeake Biological Laboratory, Center for Environmental and Estuarine Studies, University of Maryland System, Solomons. McGillivary, Mike Lee, Henry Chang Wo, David Kimo Frankel, the Lapointe, B.E., 1997. Nutrient thresholds for bottom-up control of macroalgal Office of Hawaiian Affairs, Dr. Robert Richmond of the University of blooms on coral reefs in Jamaica and southeast Florida. Limnol. Oceanogr. Hawaii, Manoa, and residents of Ewa Beach that provided access to 42, 1119–1131. Lapointe, B.E., Bedford, B.J., 2007. Drift rhodophyte blooms emerge in Lee County, their properties. The paper was improved by the constructive FL, USA: evidence of escalating coastal eutrophication. Harmful Algae 6, 421– comments of several anonymous reviewers. This is contribution # 437. 1817 from the Harbor Branch Oceanographic Institute at Florida Lapointe, B.E., Bedford, B.J., 2010. Ecology and nutrition of invasive Caulerpa Atlantic University, Ft. Pierce, FL.[SS] brachypus f. parvifolia blooms on coral reefs off southeast Florida, USA. Harmful Algae 9, 1–12. Lapointe, B.E., Littler, M.M., Littler, D.S., 1987. A comparison of nutrient-limited References productivity in macroalgae from a Caribbean Barrier Reef and from a mangrove ecosystem. Aquat. Bot. 28, 243–255. Abbott, I.A., 1984. Limu: An Ethnobotanical Study of Some Hawaiian Seaweeds, 4th Lapointe, B.E., O’Connell, J.D., Garrett, G.S., 1990. Nutrient couplings between on- edition. National Tropical Botanical Garden, Kauai. site sewage disposal systems, groundwaters, and nearshore surface waters of Abbott, I.A., 1999. Marine Red Algae of the Hawaiian Islands. Bishop Museum Press, the Florida Keys. Biogeochemistry 10, 289–307. Honolulu. Lapointe, B.E., Littler, M.M., Littler, D.S., 1992. Nutrient availability to marine Abbott, I.A., Huisman, J.M., 2004. Marine Green and Brown Algae of the Hawaiian macroalgae in siliciclastic versus carbonate-rich coastal waters. Estuaries 15, Islands. Bishop Museum Press, Honolulu. 75–82. Asplia, I., Agemian, H., Chau, A.S.Y., 1976. A semi-automated method for the Lapointe, B.E., Barile, P.J., Yentsch, C.S., Littler, M.M., Littler, D.S., Kakuk, B., 2004. The determination of inorganic, organic, and total phosphate in sediments. Analyst relative importance of nutrient enrichment and herbivory on macroalgal com- 101, 187–197. munities near Norman’s Pond Cay, Exumas Cays, Bahamas: a ‘‘natural’’ enrich- Atkinson, M., Smith, S.V., 1983. C:N:P ratios of benthic marine plants. Limnol. ment experiment. J. Exp. Mar. Biol. Ecol. 298 (2), 275–301. Oceanogr. 28 (3), 568–774. Lapointe, B.E., Barile, P.J., Littler, M.M., Littler, D.S., 2005. Macroalgal blooms on Banner, A.H., 1974. Kaneohe Bay, Hawaii: urban pollution and a coral reef ecosys- southeast Florida coral reefs II. Cross-shelf discrimination of nitrogen sources tem. In: Proc. 2nd Int. Symp. Coral Reefs, vol. 2. Great Barrier Reef Committee, indicates widespread assimilation of sewage nitrogen. Harmful Algae 4, 1106– Brisbane. 1122. Bruland, G.L., MacKenzie, R.A., 2010. Nitrogen source tracking with d15N content of Lapointe, B.E., Langton, R., Bedford, B.J., Potts, A.C., Day, O., Hu, C., 2010. Land-based coastal wetland plants in Hawaii. J. Environ. Qual. 39, 409–419. nutrient enrichment of the Buccoo Reef Complex and fringing coral reefs of Carlton, J.T., Geller, J.B., 1993. Ecological roulette: the global transport of non- Tobago, West Indies. Mar. Poll. Bull. 60, 334–343. indigenous marine organisms. Science 261, 78–82. Larned, S.T., 1998. Nitrogen versus phosphorus limited growth and sources of Chisholm, J.R.M., Fernex, F.E., Mathieu, D., Jaubert, J.M., 1997. Wastewater dis- nutrients for coral reef macroalgae. Mar. Biol. 132, 409–421. charge, seagrass decline and algal proliferation on the Cote d’ Azur. Mar. Poll. Laws, E.A., Ferentinos, L., 2003. Human impacts on fluxes of nutrients and sediment Bull. 34, 78–84. in Waimanalo stream, O’ahu, Hawaiian Islands. Pac. Sci. 57, 119–140. 318 B.E. Lapointe, B.J. Bedford / Harmful Algae 10 (2011) 310–318
Laws, E., Ziemann, D., Schulman, D., 1999. Coastal water quality in Hawaii: the Russell, D.J., 1992. The ecological invasion of Hawaiian reefs by two marine red importance of buffer zones and dilution. Mar. Environ. Res. 48, 1–21. algae Acanthophora spicifera (Vahl) Boerg. and Hypnea musciformis (Wulfen) J. Littler, M.M., Littler, D.S., 1985. Nondestructive sampling. In: Littler, M.M., Littler, Ag., and their association with two native species, Laurencia nidifica J. Ag. and D.S. (Eds.), Handbook of Phycological Methods, Ecological Field Methods: Hypnea cervicornis. J. Agric. ICES Mar. Sci. Symp. 194, 110–125. Macroalgae. Cambridge University Press, Cambridge, pp. 161–175. Smith, A.A.C., 1994. The effect of nutrient loading on C:N:P ratios of marine (MEA) Millenium Ecosystem Assessment, 2005. Ecosystems and Human Well- macroalgae in Kanehoe Bay, Hawaii. MS Thesis, University of Hawaii, Honolulu. Being: Biodiversity Synthesis. World Resources Institute, Washington, DC. Smith, S.V., Kimmerer, W.J., Laws, E.A., Brocknand, R.E., Walsh, T.W., 1981. Kaneohe (NRC) National Research Council, 2000. Clean Coastal Waters: Understanding and bay sewage diversion experiment: perspectives on ecosystem responses to Reducing the Effects of Nutrient Pollution. Ocean Studies Board, Water Science nutritional pertubation. Pac. Sci. 35 (4), 279–402. and Technology Board, National Academy Press, Washington, DC. Smith, J.E., Hunter, C.L., Smith, C.M., 2002. Distribution and reproductive character- Parsons, M.L., Walsh, W.J., Settlemier, C.J., White, D.J., Ballauer, J.M., Ayotte, P.M., istics of nonindigenous and invasive marine algae in the Hawaiian Islands. Pac. Osada, K.M., Carman, B., 2008. A multivariate assessment of the coral ecosystem Sci. 56, 299–315. health of two embayments on the lee side of the island of Hawaii. Mar. Poll. Bull. Soegiarto, A., 1972. The role of benthic algae in the carbonate budget of the modern 56, 1138–1149. reef complex, Kaneohe Bay. Thesis, University of Hawaii, Honolulu. Presley, T.K., Jamison, M.T.J., 2009. Rainfall, discharge, and water quality data during Staples, G.W., Cowie, R.H., 2001. Hawaii’s Invasive Species. Bishop Museum Press, stormwater monitoring, July 1, 2008, to June 30, 2009, Halawa stream drainage Honolulu. basin and the H-1 storm drain, Oahu, Hawaii. United States Geological Survey, (UNEP) United Nations Environmental Programme, 2006. Marine and Coastal Open-File Report, 2009-1162. . Ecosystems and Human Well-being: A Synthesis Report Based on the Findings Risk, M.J., Lapointe, B.E., Sherwood, O.A., Bedford, B.J., 2009. The use of d15Nin of the Millenium Ecosystem Assessment. . assessing sewage stress on coral reefs. Mar. Poll. Bull. 58, 793–802. United States Geological Survey, 2010. http://nwis.waterdata.usgs.gov/hi/nwis/ Rodgers, S.K., Cox, E.F., 1999. Rate of spread of introduced rhodophytes Kappa- qwdata/?site_no=212353157533001&agency_cd=USGS. phycus ahwarrzii, Kappaphycus striatum,andGracilaria salicornia and their Van Houtan, K., Hargrove, S.K., Balzas, G.H., 2010. Land use, macroalgae, and tumor- current distributions in Kaneohe Bay, Oahu, Hawaii. Pac. Sci. 53, 232–241. forming disease in marine turtles. PLOS One 5 (9), e12900. Ruddle, K., Johannes, R.E., 1985. The traditional knowledge and management of Wanielista, M.P., Yousef, Y.A., 1993. Stormwater Management. John Wiley and Sons, coastal systems in Asia and the Pacific. UNESCO, Jakarta. New York. Russell, D.J., 1987. Introductions and establishment of alien marine algae. Bull. Mar. Williams, S.L., Smith, J.E., 2007. A global review of the distribution, taxonomy, and Sci. 42, 641–642. impacts of introduced seaweeds. Annu. Rev. Ecol. Evol. Syst. 38, 327–359. See discussions, stats, and author profiles for this publication at: https://www.researchgate.net/publication/313480082
Sources and dispersal of land-based runoff from small Hawaiian drainages to a coral reef: Insights from geochemical signatures
Article · February 2017 DOI: 10.1016/j.ecss.2017.02.013
CITATIONS READS 0 60
2 authors, including:
Curt D Storlazzi United States Geological Survey 126 PUBLICATIONS 2,046 CITATIONS
SEE PROFILE
Some of the authors of this publication are also working on these related projects:
Coral Reefs View project
Quantifying the long-term anthropogenic and natural effects of sediment resuspension along the coast, shelf and slope View project
All content following this page was uploaded by Curt D Storlazzi on 01 December 2017.
The user has requested enhancement of the downloaded file. Estuarine, Coastal and Shelf Science 188 (2017) 69e80
Contents lists available at ScienceDirect
Estuarine, Coastal and Shelf Science
journal homepage: www.elsevier.com/locate/ecss
Sources and dispersal of land-based runoff from small Hawaiian drainages to a coral reef: Insights from geochemical signatures
* Renee K. Takesue , Curt D. Storlazzi
U.S. Geological Survey, Pacific Coastal and Marine Science Center, 2885 Mission Street, Santa Cruz, CA 95060, USA article info abstract
Article history: Land-based sediment and contaminant runoff is a major threat to coral reefs, and runoff reduction efforts Received 12 September 2016 would benefit from knowledge of specific runoff sources. Geochemical signatures of small drainage Received in revised form basins were determined in the fine fraction of soil and sediment, then used in the nearshore region of a 3 February 2017 coral reef-fringed urban embayment on southeast Oahu, Hawaii, to describe sources and dispersal of Accepted 6 February 2017 land-based runoff. The sedimentary rare earth element ratio (La/Yb) showed a clear distinction between Available online 8 February 2017 N the two main rock types in the overall contributing area, tholeiitic and alkalic olivine basalt. Based on this geochemical signature it was apparent that the majority of terrigenous sediment on the reef flat origi- Keywords: fi Coral reef nated from geologically old tholeiitic drainages. Sediment from one of ve tholeiitic drainages had a fl Runoff distinct geochemical signature, and sediment with this signature was dispersed on the reef at 2 km Sediment provenance west and 150 m offshore of the contributing basin. Sediment and the anthropogenic metals Cd, Pb, and Rare earth elements Zn were entrained in runoff from the most heavily urbanized region of the watershed. Although Anthropogenic metals anthropogenic Cd and Zn had localized distributions close to shore, anthropogenic Pb was found asso- USA ciated with fine sediment on the westernmost part of the reef flat and 400 m offshore, illustrating how Hawaii trade-wind-driven sediment transport can increase the scale of runoff impacts to nearshore commu- Oahu nities. Our findings show that sediment geochemical signatures can provide insights about the source Maunalua Bay and dispersal of land-based runoff in shallow coastal environments. The application of such knowledge to watershed management and habitat remediation efforts can aid in the protection and restoration of runoff-impacted coastal ecosystems worldwide. Published by Elsevier Ltd.
1. Introduction smothers live corals and exposes them to contaminants, and it is associated with lower live coral cover, lower species diversity, and Coastal areas have long been desirable locations for human degraded fisheries (e.g., Fabricius, 2005; Knowlton, 2001; Rogers, settlement and economic activity (Mee, 2012). Human use of 1990). Runoff prevention and control measures in watersheds up- coastal areas and watersheds can, however, exert a heavy toll on stream of coral reefs are recognized globally as means to prevent ecosystems by altering natural processes and habitats and over- further degradation and facilitate recovery of runoff-impacted coral using resources (Jackson, 2008; Lotze et al., 2006). For example, reefs (e.g., Bartley et al., 2014; Hughes et al., 2010; Richmond et al., increased erosion and runoff of sediment, nutrients, and contami- 2007). In addition, the coupling of land-based runoff management nants can degrade coastal water quality and lead to the loss of vital with an understanding of local hydrodynamic controls on near- ecosystems including coral reefs (Fabricius, 2005; Pandolfi et al., shore sediment and contaminant transport is an important 2005), seagrasses (Short and Wyllie-Echeverria, 1996; Waycott component of comprehensive and effective remediation strategies et al., 2009), and kelp forests (Jackson et al., 2001). Land-based for runoff-impacted reefs (Done, 1995; Hunter and Evans, 1995; sediment and contaminant runoff is harmful to coral reefs in Restrepo et al., 2016). Recent studies have demonstrated the use many ways: it inhibits photosynthesis and larval recruitment, it of sediment-geochemical tracers in identifying sources of land- derived sediment to the coastal zone (Araújo et al., 2002; Prego et al., 2009, 2012; Roussiez et al., 2013; Smith et al., 2008). The goals of this study were to identify geochemical signatures of * Corresponding author. E-mail addresses: [email protected] (R.K. Takesue), [email protected] terrigenous sediment and trace metal runoff to a coral reef-fringed (C.D. Storlazzi). http://dx.doi.org/10.1016/j.ecss.2017.02.013 0272-7714/Published by Elsevier Ltd. 70 R.K. Takesue, C.D. Storlazzi / Estuarine, Coastal and Shelf Science 188 (2017) 69e80 urbanized embayment from several short, steep drainages; to use Hawaii Kai Marina complex (Coles et al., 2002). Urban stream these signatures to identify sediment sources and infer nearshore sediment and roadside soil in southeast Oahu have been found to transport; and to describe the distribution of anthropogenic metals contain elevated anthropogenic trace metals (De Carlo et al., 2005; in urbanized basins and the reef flat. The trace metals cadmium Sutherland and Tolosa, 2000). (Cd), copper (Cu), lead (Pb), and zinc (Zn) are enriched by anthro- Maunalua Bay lies in the lee of Koolau Ridge relative to the di- pogenic activities, particularly the operation of motor vehicles rection of the northeast trade winds. Rainfall in the watershed is (Alloway, 1995). Insights from geochemical signatures about sour- higher in winter than in summer due to a higher frequency of ces of land-based sediment and contaminants and their nearshore southerly (Kona) storms and other low pressure disturbances dispersal can help runoff management efforts target priority (Giambelluca et al., 2013; Oki and Brasher, 2003). At high elevations contributing areas and guide nearshore habitat remediation of on Koolau Ridge where many streams have their headwaters, mean sediment-impacted coastal ecosystems. annual rainfall is approximately 1500e2000 mm, whereas rainfall on the urbanized coastal plain is about half that amount (Giambelluca et al., 2013; Oki and Brasher, 2003). Sediment 2. Site description retention structures were built in valleys above residential areas to control runoff and stream channels have been straightened and 2.1. Environmental setting hardened to varying extents. Storm runoff is flashy in nature (Tomlinson and De Carlo, 2003) and exacerbated by the high degree Maunalua Bay is an urbanized embayment on the southeast of impervious surface in urban areas (Wolanski et al., 2009). shore of the Island of Oahu, Hawaii, U.S.A. (Fig. 1). Hawaiian wa- The Maunalua Bay reef flat ranges from 0.2 to 1.0 km wide and is tersheds are generally composed of several small, steep valleys approximately 1 m deep and 10 km long. It is subject to water with streams that enter a larger body of water. Maunalua Bay is one quality impairments due to elevated nutrients and chlorophyll, such body, receiving runoff from 10 small drainage basins with a non-native algae, low live coral cover, and diminished fish and total contributing area of 57 km2. A major highway parallels about seagrass communities (Coles et al., 2002). half of the shoreline (Fig. 1) and dense residential development Currents in Maunalua Bay are driven by winds and tides (Presto occupies the coastal plain, valley floors, and some ridges. Kuapa et al., 2012; Storlazzi et al., 2010). At the surface to a depth of 1 m, Pond, a 2 km2 shallow lagoon in the east part of the watershed, was the prevailing trade winds drive westward transport (Presto et al., breached permanently, its marshland filled, and developed into the
Fig. 1. Shaded relief map of southeast Oahu showing Maunalua Bay (white line marks the shoreline) and its watershed (dashed line). Drainage basins are identified with numbers and abbreviations described in Section 3.1. Streams in basins 1e5 are shown by name. The sub-basin 3-WPE-W (Wiliwilinui) is denoted by a ‘w’ inside a gray triangle. Black triangles show USGS gaging stations (GS) on Wailupe and Waiakeakua Streams. ‘terr’, terrestrial; ‘ns’, nearshore. Black lines on the reef flat show reef transects. The 5 m isobath is shown for reference. Inset shows the location of Maunalua Bay in the Main Hawaiian Islands. R.K. Takesue, C.D. Storlazzi / Estuarine, Coastal and Shelf Science 188 (2017) 69e80 71
2012), and the corresponding transport of suspended sediment from Hawaii Loa Ridge to address the local perception that con- over the reef flat is to the west (Storlazzi et al., 2010). Wind-driven struction activities made it a source of runoff to the reef. The lower currents are weakest in the eastern part of the bay in the lee of Koko reaches of Wailupe, Niu, and Kuli'ou'ou Streams were sampled from Head and strongest in the middle of the bay offshore of Wailupe highway bridges, and marina bottom sediment from a small boat, ® Stream (Storlazzi et al., 2010). using a hand-deployed Petite Ponar benthic sampler. The upper 0e5 cm of seabed sediment on the reef flat was collected by hand at 2.2. Geology low tide in acid-cleaned polypropylene jars along shore-normal transects offshore of Black Point, Wai'alae Nui Stream, Wiliwilinui The uplands around Maunalua Bay consist of layered tholeiitic Stream, Wailupe Stream, Hawaii Loa Ridge, Kuli'ou'ou Stream, and basalt lava flows from the Makapuu stage of Koolau Volcano Portlock at distances of 0, 50, 100, 250, and 400 m, where possible. (Haskins and Garcia, 2004), which erupted 2e3Ma(Jackson et al., Sampling locations and the numbers of samples are shown in Fig. 1 1999; Stearns and Vaksvik, 1935). Lava beds dip 4e10 to the south and Table 1. toward Maunalua Bay (Wentworth and Winchell, 1947). After about The 10 small basins that drain to Maunalua Bay and Hawaii Loa a million years of quiescence, subsidence, and erosion of Koolau Ridge are for convenience denoted with numbers increasing from Volcano (Gramlich et al., 1971), alkalic basalts erupted from rift west to east and an abbreviation: Wai'alae Nui (1-WAN), Wai'alae zones and vents, forming the Honolulu Volcanics. Around Maun- Iki (2-WAI), Wailupe (3-WPE), which includes Wiliwilinui (3-WPE- alua Bay, rocks of the Honolulu Volcanics form Koko Head, Koko W), Hawaii Loa Ridge (3.5-HLO), Niu (4-NIU), Kuli'ou'ou (5-KUL), Crater, and Black Point (Clague and Frey, 1982; Stearns and Vaksvik, Ka'alakei (6-KAL), Haha'ione (7-HAH), Kamilo Nui (8-KAN), Kamilo 1935; Winchell, 1947)(Fig. 1). Mantle-incompatible trace element Iki (9-KAI), and Portlock (10-KOK) (Fig. 1, Table 1). Runoff from 6- contents of Makapuu stage lavas can vary widely (Huang and Frey, KAL, 7-HAH, 8-KAN, and 9-KAI enters the Hawaii Kai Marina 2005) but do not vary systematically with age (Jackson et al., 1999). (MAR) and can become trapped there until removal from the sys- Compared to Koolau basalt, rocks of the Honolulu Volcanics are tem by dredging. Because the linkage between runoff into the enriched in alkali and alkali earth elements, light rare earth ele- marina and discharge into Maunalua Bay may be weak, sediment ments, and incompatible trace elements (Clague and Frey, 1982; geochemistry from these basins and the marina will not be dis- Roden et al., 1984). cussed in detail.
3. Approach 3.2. Sediment geochemical analyses
The element aluminum (Al) is a major component of terrestrial Geochemical analyses were performed on the sediment fine sediment (Windom et al., 1989) and a trace component of marine fraction (particle diameter <63 mm) to reduce grain size bias among carbonates (Milliman and Syvitski, 1992). This orders-of-magnitude samples. Bulk sediment was dried at 60 C, disaggregated gently to difference makes Al a sensitive indicator of terrigenous sediment in preserve original grains in an acid-cleaned agate mortar and pestle, environments where carbonate sediment predominates. In other and dry-sieved using stainless steel sieves to obtain the <63 mm coastal environments, Al contents of coastal sediment have been fraction. Fine sediment was decomposed according to EPA Method used to indicate riverine sediment inputs in the Gulf of Papua 3052, a near-total microwave-assisted digestion of siliceous (Brunskill et al., 1995) and the Portuguese shelf (Araújo et al., 2002). matrices using hydrochloric and hydrofluoric acids (USEPA, 1996). Furthermore, Al is not generally enriched by anthropogenic activ- Contents of major, minor, and trace elements including Al, Ba, Ca, ities (Windom et al., 1989) nor is it altered by oxidation-reduction Cd, Co, Cu, Fe, K, Mg, Mn, Na, Pb, Sr, Ti, V, and Zn were determined processes as is the element iron (Fe), another major component on a ThermoFinnigan Element I High Resolution inductively- of basalt. coupled plasma mass spectrometer (ICP-MS) at the Marine To characterize geochemical signatures of drainages contrib- Analytical Laboratories of the University of California at Santa Cruz. uting sediment to Maunalua Bay, immobile and relatively immobile Internal standardization was with germanium (72Ge); external trace element contents were determined in the fine fraction of standardization was with sediment reference materials (SRMs): stream sediment, where possible, to exploit the integrative nature National Institute of Standards and Technology 1646a and 2702 and over space and time of sediment stored in streams (Frissell et al., Canadian Certified Reference Materials Stream Sediment 2 and 3. 1986). Rare earth elements (REE), scandium (Sc), and thorium The reproducibility, expressed by the relative standard deviation (Th) are the most effective sediment provenance indicators (RSD) of a consistency standard analyzed five times, was 6% or (McLennan, 1989; McLennan et al., 1993). Barium (Ba), chromium better for all target elements except strontium (Sr) which had a RSD (Cr), cobalt (Co), nickel (Ni), niobium (Nb), rubidium (Rb), and zir- of 10%. The reproducibility of Fe measurements was better than 3%. conium (Zr) are can also be informative but could be affected by Contents of target elements in 86 soil and sediment samples sorting and weathering (McLennan et al., 1993). The anthropogenic analyzed on the Element I were several orders of magnitude higher trace metals cadmium (Cd), copper (Cu), lead (Pb), and zinc (Zn) than analytical detection limits, which were defined as three times were explored as urban overprints on sediment geochemistry in the standard deviation of blanks. Contents of rare earth elements drainage basins surrounding Maunalua Bay and on the reef flat. (REE), elements in resistant minerals (Cr, Hf, Nb, Y, Zr); elements with inadequate external standards (La, Rb, Th); and nickel (Ni) and 3.1. Sediment collection scandium (Sc), which may have had molecular interferences, were obtained by total digestion of 70 samples that had sufficient fine A low degree of fine-grained sediment storage in the short, material for quantification by SGS Inc., a nationally-recognized hardened stream channels around Maunalua Bay in many cases testing laboratory. SGS used a sodium peroxide sinter for total necessitated sediment collection from retention basins and from digestion and quantified major and minor elements by ICP-AES features such as culverts on the urbanized coastal plain. The upper (atomic emission) and trace elements by ICP-MS. The reproduc- 1e3 cm of sediment, where possible, composed of redeposited soil ibility of three replicate samples analyzed by SGS was better than and organic matter were collected in catchments 7e10 June 2010 10% for all target elements except the rare earth element thulium using acid-cleaned polypropylene sampling tools. Culverts gener- (Tm), which had a reproducibility of 12%. Target element contents ally contained only a veneer of sediment. In situ soil was collected determined by SGS were five or more times higher than analytical 72 R.K. Takesue, C.D. Storlazzi / Estuarine, Coastal and Shelf Science 188 (2017) 69e80
Table 1 Drainage basin codes, names, area estimates, numbers of samples, and numbers of rare earth element (REE) analyses in the Maunalua Bay watershed. ‘TERR’, terrestrial; ‘NS’, nearshore.
Basin code Basin name Drainage area Number of samples Number of REE analyses
(km2)a TERR NSb TERR NS
1-WAN Wai'alae Nui 11.2 5 7 5 3 2-WAI Wai'alae Iki 4.6 2 2 2 1 3-WPE-W Wiliwilinuic 1.0 2 7 2 6 3-WPE Wailupe 9.5 4 6 3 2 3.5-HLO Hawaii Loa Ridge 0.7 4 4 4 3 4-NIU Niu 7.4 5 1 5 1 5-KUL Kuli'ou'ou 4.7 3 7 2 6 MAR Marinad 16.0 2 17 2 17 10-KOK Portlock 1.8 4 4 4 2
a Estimated from ridge to shoreline from Google Earth Pro. b Nearshore includes estuarine sites. c Wiliwilinui area included in Wailupe. d Marina drainages include Ka'alakei, Ha'ha'ione, Kamilo Nui, Kamilo Iki.
detection limits. Major element contents are reported as weight 3.4. Environmental data percent (wt %) and minor and trace element contents as micro- grams per gram (mg/g). Geochemical data are tabulated in the Time series of daily stream discharge and total suspended supplementary material. sediment concentration were obtained from USGS gaging station 16247550 (Wailupe Gulch at E. Hind Drive Bridge) for the two years 3.3. Normalization and criteria for geochemical signatures preceding the study. Annual stream discharge at station 16240500 (Waiakeakua Stream at the head of nearby Manoa Valley) over the Elemental contents of fine sediment were normalized to iron past three decades provided a long-term context for southeast (Fe) to account for basaltic parent rock compositions, a convention Oahu. Wind data were obtained from NOAA National Ocean Service used in Hawaiian soil and sediment (e.g., De Carlo and Spencer, (NOS) Station 1612340 (Honolulu, Hawaii), approximately 15 km 1995) because it can contain primary volcanic minerals (Nelson west of Maunalua Bay. et al., 2013). REE contents were normalized to a North American shale composite (NASC), which represents the composition of average sedimentary rock (McLennan, 1989). Normalized values are 4. Results denoted with subscripts MFe or REEN. Statistics were calculated with StatPlus:mac Pro software. REE ratios were normally distrib- 4.1. Environmental conditions uted in basaltic and in alkalic fine sediment (Shapiro-Wilk). Anal- ysis of variance (ANOVA) was performed on log-transformed data In the 8 months preceding this study, northeast trade winds for non-normally distributed parameters. were prevalent 55% of the time and there were only 16 d when the Geochemical tracers used for sediment-source attribution mean daily wind direction was from the west (Fig. 2). Water year should vary little within each source region so that end members (WY) 2010, which began 1 October 2009, was a dry year. Honolulu are well-constrained, and considerably among source regions to received about half the amount of precipitation as in the preceding allow discrimination (Collins and Walling, 2002). Ratios of 29 major year (Presley and Jamison, 2010). Discharge in Wailupe Stream was and trace elements relative to Al, Fe, Nb, Sc, and Th, and seven REE four times lower than in WY 2009, and discharge in Waiakekua ratios were examined in the initial data review. A high degree of Stream was four times lower than the decadal average (Oki and chemical weathering has been observed in Hawaiian soil and Brasher, 2003). suggests that only immobile element ratios are representative of source compositions (Kurtz et al., 2000; Vitousek et al., 1997). Accordingly, further data exploration focused on immobile ele- 4.2. Terrigenous sediment on the reef flat ments (REE ratios, ScFe,ThFe) and relatively immobile elements (BaFe,CrFe,CoFe,NiFe,NbFe,RbFe,ZrFe) as potential geochemical Aluminum contents of fine-grained sediment on the reef flat signatures in sediment derived from the two types of basalt in the ranged from 0.6 to 6.6 wt % with a mean of 2.7 ± 1.5 wt % (1s) Maunalua Bay watershed, henceforth denoted by THO (sediment compared to a mean of 9.8 ± 1.2 wt % (1s)infine-grained upland derived from tholeiitic basalt, n ¼ 25) and AOB (sediment derived sediment. Assuming fine sediment on the reef flat was a mixture from alkalic olivine basalt, n ¼ 4). One criterion was that the primarily of terrigenous and carbonate material, and that the Al geochemical property should have low coefficients of variation (CV, content of terrigenous sediment on the reef flat was similar to that defined as ratio of the standard deviation to the mean) in AOB and of upland sediment, then 6e66% (mean of 27%) of fine sediment on THO. The other was that the geochemical property should have a the reef flat was land-derived. Fine-sediment Al contents generally large enrichment factor (EF) between AOB and THO. EFs were decreased with distance from shore, except offshore of Wiliwilinui calculated as the ratios of median values (EF ¼ medianAOB/ and Portlock, where Al maxima occurred at 50 m. The three near- medianTHO). A conservative value of 2 was used for the minimum EF shore sites with the largest fractions of terrigenous sediment based criterion, which yielded geochemical signatures that were on on their Al contents were, in decreasing order: the mouth of Wai- average two or more times higher in AOB than THO. EFs of lupe Stream (66%), 50 m offshore of Wiliwilinui Stream (58%), and immobile and relatively immobile elements in AOB and THO were the mouth of Kuli'ou'ou Stream (44%). At two reef sites 400 m from plotted relative to the corresponding CVs in order to identify the the shore, the fine fraction of bed sediment contained only 1.2% and most effective geochemical tracers. 1.4% terrigenous material. R.K. Takesue, C.D. Storlazzi / Estuarine, Coastal and Shelf Science 188 (2017) 69e80 73
Fig. 2. Time series plots of wind direction (top) and wind speed (middle) at NOS Meteorological Station 1612340 in Honolulu. Bottom panel shows stream discharge and total suspended sediment concentration (TSS) in Wailupe Stream at USGS gaging station 16247550 from 1 October 2009 through 30 June 2010.
4.3. Geochemical signatures were evaluated. (La/Yb)N ratios in THO and AOB fell along distinctly different Fine-sediment contents of ScFe,CrFe,CoFe,NiFe, and ZrFe were trends (Fig. 4). Fine sediment at 31 of 36 reef sites had (La/Yb)N relatively similar in AOB and THO (0.8 < EF < 1.3, Fig. 3) and so were not investigated further as geochemical signatures. RbFe had EF ¼ 10.2, however its high CVs precluded its ability to constrain end member compositions (Fig. 3). The median BaFe value was almost five times higher in AOB relative to THO and its CVs were intermediate (Fig. 3). EF values for (La/Yb)N,NbFe, and ThFe were 2.4, 2.5, and 2.4, respectively (Fig. 3). Because (La/Yb)N had the lowest CVs and EF > 2(Fig. 3), it was selected as the geochemical signature by which sources of land-based fine sediment to Maunalua Bay
Fig. 4. Plot of fine-sediment lanthanum (LaN) relative to ytterium (YbN) contents in terrestrial (large symbols) and marine environments (small symbols). The average composition and standard deviation (error bars) of tholeiitic basalt (THO) are shown Fig. 3. Comparisons of enrichment factors (EF) of Fe-normalized immobile and rela- for reference, data from Frey et al. (1994). The dashed line shows the least-squares tively immobile elements in the fine fraction of terrestrial soil and sediment derived regression through HLO (r ¼ 0.9980,n¼ 4) extrapolated from YbN ¼ 0.5 to low YbN from tholeiitic (THO) and alkalic olivine basalt (AOB) relative to their coefficients of values. The compositions of alkalic olivine basalts (AOB) are shown for reference, data variation (CV). Circles show groupings for THO and AOB corresponding to the labeled from Clague and Frey (1982) with a least-squares regression line (r ¼ 0.9949,n¼ 3, element. Dashed line shows EF ¼ 2. Symbols below the line are for ScFe,CrFe,CoFe,NiFe, dotted line). The star shows the composition of Asian dust in North Pacific pelagic and ZrFe. sediment near Hawaii (Nakai et al., 1993). 74 R.K. Takesue, C.D. Storlazzi / Estuarine, Coastal and Shelf Science 188 (2017) 69e80 compositions that more closely followed the trend of tholeiitic orographic rainfall and thus varies spatially and with elevation on basalt than alkalic olivine basalt (Fig. 4). An estuarine sample Oahu (Jackson et al., 1971). In the contributing area of Maunalua collected in Niu Stream (4-NIU) had a (La/Yb)N ratio indicative of an Bay, dust deposition should be greater at high elevations than on AOB source whereas Niu basin consisted of tholeiitic basalt. Stream the coastal plain and in western drainages compared to eastern sediment in Wiliwilinui basin (3-WPE-W) had elevated LaN con- ones (Giambelluca et al., 2013). Following the approach of Kurtz tents relative to YbN values (LaN ¼ 0.5e0.6), representing an et al. (2001) who used soil Th/Nb ratios to estimate the Asian enrichment of 0.1e0.2 units relative to 3-WPE sediment (Fig. 4); dust contribution to soil on the Island of Hawaii, the fine fraction of this difference was significant in a one-way ANOVA (p ¼ 0.02). This soil and sediment collected for this study was estimated to contain enrichment is also apparent at five sites on the reef flat: 0 m and 0e15% Asian dust, assuming Th/Nb ratios of 0.52 for Asian dust 50 m offshore of the mouth of Wiliwilinui Stream (3-WPE-W), at (Kurtz et al., 2001), 0.05 for Koolau basalt (Frey et al., 1994), and the shoreline of 2-WAI, and at the shoreline and 150 m offshore of 0.08 for alkali olivine basalt and basanite (Clague and Frey, 1982). 1-WAN (Fig. 4). One sample from 5-KUL had an elevated Th content and Th/Nb ratio Stream sediment from tholeiitic basins 1e5 and soil from Hawaii that yielded an estimated dust contribution of 36%, but its (La/Yb)N Loa Ridge could not be distinguished individually based on their ratio was not indicative of a large dust fraction so it was considered (La/Yb)N ratios, except for Wiliwilinui. Instead, BaFe,NbFe, and ThFe an outlier and disregarded. The mean ± 1s estimated dust content were examined as potential basin-specific geochemical signatures was lowest at 10-KOK (4 ± 2%), the driest and easternmost basin, because these immobile or relatively immobile elements were and highest on Hawaii Loa Ridge (12 ± 4%), the highest elevation more variable among tholeiitic fine sediment. There were 2e5 site. There were no other longitudinal patterns in the mean esti- samples from individual source areas (Table 1). The fine fraction of mated dust content of fine sediment in drainage basins surround- soil and sediment from basins 1e5 and Hawaii Loa Ridge had ing Maunalua Bay: 1-WAN (7 ± 2%), 2-WAI (7 ± 1%), 3-WPE, generally similar median BaFe and NbFe values with the exception of including 3-WPE-W (7 ± 4%), 4-NIU (8 ± 5%), 5-KUL (7%). NbFe at Wiliwilinui (Fig. 5). Median ThFe contents of fine sediment Wiliwilinui (3-WPE-W) was the only basin whose sediment had appeared to be slightly lower in basins 1 and 2 compared to Hawaii (La/Yb)N ratios that were intermediate between those of AOB and Loa Ridge and basins 4, 5 (Fig. 5). However, log-transformed BaFe, THO (Fig. 5). Asian dust [(La/Yb)N ¼ 1.22 (Nakai et al., 1993),] would NbFe, and ThFe contents of fine sediment were statistically indis- increase the (La/Yb)N ratio of soil developed on THO [0.67 (Frey tinguishable among individual tholeiitic basins except for NbFe at et al., 1994),] and decrease the (La/Yb)N ratio of soil developed on Wiliwilinui (one-way ANOVA with Fisher LSD post-hoc analysis for AOB [1.50 (Clague and Frey, 1982),]. Based on the percentages of p < 0.01), as were other commonly used sediment provenance in- Asian dust estimated from the Th/Nb ratios of two samples from dicators such as log-transformed La/Th, Th/Sc, and Sm/Nd and Eu/ Wiliwilinui Stream, 7% and 14%, Asian dust inputs could explain * Eu . Thus, except in the case of Wiliwilinui, it was not possible to approximately half of the increase of (La/Yb)N ratios relative to THO geochemically discriminate sediment from individual tholeiitic in 3-WPE-W. basins and Hawaii Loa Ridge based on this set of samples and these In the nearshore region, fine sediment with elevated (La/Yb)N geochemical parameters. This result was not wholly unexpected ratios relative to THO were observed at the highway bridge over Niu because immobile elements were found to covary in Makapuu Stream and to the west of 3-WPE-W. No sediment in Niu basin had stage lavas (Huang and Frey, 2005). similarly elevated values, making it likely that the higher ratio re- flected the presence of AOB rather than Asian dust. Estimated dust contents of nearshore fine sediments were insufficient to account 4.4. Geochemical inputs to soil from Asian dust for elevated (La/Yb)N ratios offshore of 3-WPE-W, but could possibly account for those offshore of 1-KAN and 2-KAI. Dust from Asia is transported in the atmosphere across the North Pacific Ocean and deposited on the land surface of Hawaii (Dymond et al., 1974; Rex et al., 1969) where it can be incorporated 4.5. Trace metals in soil (Kurtz et al., 2001; Porder et al., 2007). Because the chemical composition of Asian dust can differ substantially from that of Chromium and Ni contents of Maunalua Bay terrestrial, estua- Hawaiian soil (e.g., Ferrat et al., 2011; Nakai et al., 1993), it can alter rine, and reef fine sediment and soil were correlated with those of geochemical mass balances in the soil column (Kurtz et al., 2001; Fe (Fig. 6). Nickel contents relative to Fe fell along a linear trend for Vitousek et al., 1997). Dust deposition is associated with all basins and depositional environments, whereas Cr contents fell
Fig. 5. Box and whisker plots of Fe-normalized barium (Ba), niobium (Nb), and thorium (Th) contents of fine-sediment from five tholeiitic basins (1e5) and Hawaii Loa Ridge (3.5). The height of the box shows the interquartile range, the whiskers show the maximum and minimum values, and the line crossing the whiskers shows the median value. The composition of sediment derived from alkalic olivine basalt (AOB) in 10-KOK is shown for comparison. R.K. Takesue, C.D. Storlazzi / Estuarine, Coastal and Shelf Science 188 (2017) 69e80 75
Fig. 6. Plots of chromium and nickel content relative to iron content for all samples. The dashed lines show the effects range median values (ERM) for reference. into two populations, one with Cr contents that increased with Wiliwilinui streams. All nearshore fine-sediment contents of Cd, increasing Fe, and one with low Cr and high Fe contents (Fig. 6). All Cu, Pb, and Zn were below marine and estuarine ERM levels. Pb was terrestrial soil and sediment had Cr contents that exceeded the the only anthropogenic trace metal that was consistently elevated probable effects concentration (PEC) for freshwater ecosystems, above background levels in estuarine and reef flat fine sediment. 111 mg/g (MacDonald et al., 2000), as did 24 for Ni [PEC ¼ 48.6 mg/g Total Pb contents ranged from 3 to 114 mg/g (Fig. 7) compared to an (MacDonald et al., 2000),]. At 25 marine and estuarine sites sedi- estuarine background of 14 ± 1 mg/g (Table 2). When normalized to ment had Cr contents that exceeded the level at which adverse Fe to account for Pb in the geologic fraction, patterns of Pb biological effects are probable in marine and estuarine environ- enrichment were less extreme closer to shore (Fig. 7), and the ments [ERM, 370 mg/g (Long et al., 1995)], as did 42 with respect to highest PbFe value was found in carbonate-dominated sediment Ni (51.6 mg/g, Fig. 6). The sites with the highest Cr contents in fine- 400 m offshore of Wiliwilinui (3-WPE-W). The PbFe ratio of that grained sediment were nearshore sites: Wiliwilinui (3-WPE-W) sample, 52, was similar to the ratio in estuarine sediment before the stream mouth (0 m) and 50 m offshore of the stream mouth, and a phase-out of leaded gasoline (De Carlo and Spencer, 1995). In storm drain 125 m to the west of Wiliwilinui stream that carried comparison, the mean Pb content of five other carbonate- runoff from Wailupe basin. dominated sediment samples on the Maunalua Bay reef flat was Among terrestrial soil and sediment samples collected in June 6 ± 1 mg/g, lower than the pre-1927 value from the Ala Wai canal 2010, approximately twice as many had Cu and Pb contents (Table 2). exceeding background levels determined in forested conservation e lands during the 1998 2000 National Water Quality Assessment on 5. Discussion Oahu than for Cd and Zn (Table 2). Trace metal enrichments were up to 6 (Cd), 1.8 (Cu), 65 (Pb), and 6 (Zn) times background levels in High-standing and young islands in the tropics and subtropics fi ne soil and sediment in the contributing area of Maunalua Bay; can undergo high soil erosion rates and contribute large material fi however, median ne-sediment contents of Cd, Cu, Pb, and Zn were fluxes to the coastal ocean (Hilton et al., 2008; Kao and Milliman, at or below background levels (Table 2). Soil and sediment with 2008; Lyons et al., 2002; Milliman and Syvitski, 1992). Such land- elevated Cd, Cu, and Zn generally occurred on land, whereas more derived runoff of sediment, nutrients, carbon, and contaminants marine sites had above-background Pb contents than terrestrial can have large impacts on global biogeochemical cycles (Milliman ones (Table 2). et al., 1999; Nittrouer et al., 1995; Sholkovitz et al., 1999) and ma- Two of 31 terrestrial sites had Pb (1 site) and Zn (2 sites) con- rine and coastal ecosystems (Restrepo et al., 2006; Waycott et al., fi tents of ne sediment that exceeded the PEC (Table 2). These were 2009). Many Hawaiian watersheds are small, steep, urbanized, culverts adjacent to a major highway on the south shore of Hawaii and coupled to downstream ecologic communities, such as coral Kai Marina and Koko District Park. Runoff near the highway had the reefs and seagrasses, and socio-economic activities such as tourism m highest Zn and Pb contents measured in this study (1253 g/g and and recreation, that rely on clean and clear water in coastal areas to m 325 g/g, respectively). The Cu background level in Maunalua soil flourish (GPA, 2006; Nurse et al., 2014). Sediment-geochemical was three times higher than the PEC (Table 2), and all but three signatures that identify areas contributing runoff could aid in the fi terrestrial sites had soil or sediment Cu contents of the ne fraction protection and restoration of coastal waters and ecosystems by m that exceeded the PEC. The highest overall Cu value (391 g/g) was identifying priority areas for management and remediation measured in forested parkland on Hawaii Loa Ridge at 340 m (Bartley et al., 2014). elevation. At no sites did soil or sediment Cd contents of the fine fraction exceed the PEC, which was more than an order of magni- 5.1. Watershed sources and nearshore dispersal of terrigenous tude higher than the Cd background in Maunalua soil (Table 2). sediment In the nearshore region, Cd, Cu, and Zn were elevated above estuarine background levels in fine sediment deposited in estua- Geochemical signatures of tholeiitic and alkalic olivine basalt rine reaches of streams, in marina bottom sediment, at the land-sea were able to distinguish fine sediment from basins in Koolau Range interface, and up to 75 m from shore offshore of Wai'alae Nui and versus a contributing area on the west slope of Koko Head. Based on 76 R.K. Takesue, C.D. Storlazzi / Estuarine, Coastal and Shelf Science 188 (2017) 69e80
Table 2 Levels of anthropogenic trace metals in unimpacted areas, sediment quality criteria, and in Maunalua Bay terrestrial and marine (reef and estuarine) sediment.
Cd (ppm) Cu (ppm) Pb (ppm) Zn (ppm)
Background levels in unimpacted areas Upland forested soila 0.20 ± 0.08 223 ± 17 5 ± 3 210 ± 27 Estuarine sediment, Ala Wai Canalb 0.29 ± 0.03 141 ± 614± 1 122 ± 11
Freshwater sediment quality criteriac
Threshold effects concentration (TEC) 0.99 31.6 35.8 121 Probable effects concentration (PEC) 4.98 149 128 459
Marine and estuarine sediment quality criteriad
Effects range low (ERL) 1.2 34 46.7 150 Effects range median (ERM) 9.6 270 218 410
Terrestrial sediment, June 2010 (n ¼ 31)
Minimum 0.10 126 3 113 Maximum 1.16 391 325 1253 Median 0.18 218 5 167 Mean 0.24 215 23 219 Standard deviation 0.19 50 61 208 # samples exceeding upland background 6 13 14 5
Reef and estuarine sediment, June 2010 (n ¼ 55)
Minimum 0.04 17 3 18 Maximum 2.36 262 114 292 Median 0.09 78 15 86 Mean 0.14 90 21 104 Standard deviation 0.31 55 21 56 # samples exceeding estuarine background 1 11 29 17
a De Carlo et al. (2005), leeward sites (n ¼ 4). b De Carlo and Spencer (1995), core G8B (107e110 cm, n ¼ 3). c Consensus-based values from MacDonald et al. (2000). d Long et al. (1995). its geochemical signature, the majority of fine terrigenous sediment basin Wiliwilinui. Fine sediment from Wiliwilinui had a distinct on the reef flat originated from basins in Koolau Range. Basin- (La/Yb)N and NbFe signature and its dispersal up to 2 km west and specific sources could not be distinguished individually using the 150 m offshore of its source was consistent with the direction of geochemical signatures explored here, except for the small sub- trade-wind-driven sediment transport on south-facing fringing Hawaiian reefs (Ogston et al., 2004; Presto et al., 2006; Storlazzi et al., 2004). The size of Wiliwilinui was small compared to other basins, but the scale of the dispersal of its runoff signature on the reef flat was not. Therefore, runoff mitigation in this small basin could result in a relatively large improvement in land-based runoff impacts on the nearshore ecologic community. The alkalic (La/Yb)N signature in fine sediment near the mouth of Niu Stream was another example of westward terrigenous sediment transport over several km. The predominance of trade- wind-driven westward sediment transport has ecological implica- tions for the reef community in Maunalua Bay in two ways. First, sediment and contaminants entering the bay will be entrained in westward flow and undergo repeated cycles of resuspension, deposition, and interaction with organisms (Ogston et al., 2004) before exiting the bay near Black Point (Presto et al., 2012; Storlazzi et al., 2010). Runoff plumes arising from winter storms, on the other hand, have short residence times in Maunalua Bay, exiting the reef rapidly in the offshore direction (Storlazzi et al., 2010; Wolanski et al., 2009). Second, runoff from the east part of the watershed near Hawaii Kai Marina, which consists of almost 20 km of shore- line with high-density residential and commercial development, had the highest levels of anthropogenic trace metals. Runoff from this region likely contains other compounds associated with ur- banization (Brasher and Wolff, 2004) such as PAHs, pesticides, flame retardants, pharmaceuticals, and personal care products that are growing concerns in urban stormwater (Daughton and Ternes, 1999; Schwarzenbach et al., 2007) and groundwater (Barnes Fig. 7. Plots of total lead (Pb) and iron (Fe)-normalized Pb contents of estuarine and et al., 2008). Urban contaminants that are transported to the reef fine sediment with distance offshore. Estuarine values are shown to the left of 0 m. coastal ocean in surface runoff or groundwater can impact ecologic Dotted lines show marine background levels (Table 2). R.K. Takesue, C.D. Storlazzi / Estuarine, Coastal and Shelf Science 188 (2017) 69e80 77 communities downstream of discharge sites. part of the watershed and a six-fold or higher enrichment over There were few rainfall events during the 2008e2009 winter geologic background levels. This is consistent with previous studies and spring preceding field sampling for this study, and flow in local showing that Cd and Zn in Hawaiian soils are anthropogenically streams was below normal, so runoff to the nearshore was also influenced (e.g., De Carlo and Spencer, 1995; De Carlo et al., 2005; likely below normal. In this context, the relatively large amount of Sutherland, 2000). Elevated fine-sediment contents of Cd and Zn fine-grained terrigenous sediment on the Maunalua Bay reef flat, in the estuarine reaches of streams and near storm drain outfalls which averaged more than 25%, was somewhat surprising, since indicate that these trace metals were entrained in runoff from the the same component averaged less than 20% of the sediment urbanized watershed and deposited at the land-sea interface, as is transported off the reef flat during winter 2008e2009 (Storlazzi typical in estuaries (Li et al., 1984; Turekian, 1977). Although et al., 2010). From the absence of strong Kona storms during nearshore fine sediment with anthropogenic Cd and Zn was winter 2009e2010, it can be inferred that there were fewer large confined to the inner 75 m of the reef flat under the dry conditions wave events that resuspended and transported sediment, resulting preceding this study, sediment and contaminant runoff and its in greater storage of terrigenous material from previous flood dispersal could be higher in wetter years. Furthermore, because events (Draut et al., 2009). Thus calm winter conditions appear to motor vehicle traffic is the primary urban source of Cd and Zn (De have contributed to the storage of land-derived sediment, and by Carlo et al., 2005; Sutherland, 2000), loading of these metals to the association sediment-bound contaminants, on the Maunalua Bay coastal ocean is expected to increase as population increases in the reef flat. state of Hawaii (HOP, 2013). Increasing runoff of sediment and ur- The naturally high and variable contents of Cr and Ni in Ha- ban contaminants is a global concern as coastal regions become waiian basalts are related to the minerals spinel and olivine, more populated and urbanized (Newton et al., 2012). respectively (Frey et al., 1994; Haskins and Garcia, 2004; Jackson Pb contents of urban soils on Oahu have been decreasing since et al., 1999). The strong correlation of Ni and Fe in terrestrial fine the phase-out of leaded gasoline in the 1980s (De Carlo and sediment was indicative of the presence of olivine (De Carlo et al., Anthony, 2002). Fine-grained sediment with high PbFe values 2005; Frey et al., 1994). When olivine is subaerially exposed it is relative to background levels indicate that legacy Pb-contaminated susceptible to alteration; however, Ni was tightly coupled to Fe in soil and sediment are still present in urban watersheds of southeast fine sediment on the reef flat, indicating it was immobile Oahu, though the Pb levels reported here are lower than during the (Marsaglia, 1993; Moberly et al., 1965) over the timescales of USGS National Water Quality Assessment a decade earlier. The Pb erosion and transport in the small drainages surrounding Maunalua content of nearshore carbonate sediment can reflect marine as well Bay and thus not likely to be bioavailable (Sutherland, 2000). The as anthropogenic processes. Carbonate has a strong affinity for Pb high Cr content of upland sediment was also due to a volcanic in seawater (Talbot and Chegwidden, 1983), and Pb enrichment of source (De Carlo et al., 2005; Frey et al., 1994), and its biomodal marine carbonate can occur by passive adsorption (Sturesson, distribution relative to Fe was indicative of two Cr-bearing min- 1976), also called scavenging, particularly when sediment is erals. Cr-spinel is ubiquitous in Hawaiian soil (Frey et al., 1994; Oze resuspended, which likely occurs almost daily on the Maunalua Bay et al., 2004) and its composition would account for sediment with a reef flat for more than 6 months of the year when trade winds high Cr to Fe ratio (Oze et al., 2004), whereas sediment with a lower prevail (Presto et al., 2006). Scavenging can account for nearshore Cr to Fe ratio could have contained Cr-bearing pyroxene (Nelson Pb contents of fine-grained carbonate sediment that were two to et al., 2013; Oze et al., 2004). Spinel and pyroxene are heavy min- three times higher than the marine background, values that were erals with densities approximately two times higher than terrige- comparable to Pb levels in calcareous sediment in other shallow nous alumiosilicates and marine carbonates, so stronger hydraulic water environments in Australia (Esslemont, 2000; Talbot and forcing is required for their suspension and transport than for Chegwidden, 1983) and Central America (Guzman and Jimenez, similarly sized particles of lower density. As a result, heavy mineral 1992). Lead enrichments in nearshore fine sediment in excess of transport can be decoupled from that of other particles in reef Pb-scavenging were attributed to anthropogenic legacy Pb from the environments (Marsaglia, 1993; Moberly et al., 1965). Cr contents of dispersal of land-based runoff across the Maunalua Bay reef flat. about 1000 mg/g found in carbonate-dominated sediment (<3% Legacy Pb was found along three of five transects on the reef flat; terrigenous) offshore of basins Wai'alae Nui and Portlock shows however, it did not occur at levels where biological impacts would that high-Cr mineral grains remained on the Maunalua Bay reef flat be expected. after other terrigenous sediment was winnowed away. Therefore the distribution of Cr in fine-grained sediment on the reef flat is not 5.3. Implications representative of overall fine-sediment transport and deposition patterns. Urbanization and development of coastal zones are increasing worldwide, and concomitant increases in material fluxes from land 5.2. Anthropogenic trace metals in the Maunalua Bay system to sea could irreparably disrupt coastal processes and ecosystems if corrective measures are not undertaken (Hughes et al., 2010; The elevated Cu content of soil and sediment in the Maunalua Jackson, 2008; Steffen et al., 2011). Sediment geochemical signa- Bay watershed, although affected by anthropogenic activities (De tures represent one means to improve the effectiveness of runoff Carlo et al., 2004), is due primarily to a large geologic component reduction and control efforts by identifying runoff-contributing (De Carlo et al., 2005), almost 60% of the total (Sutherland and areas, information that can inform watershed management prac- Tolosa, 2000), which is not biologically available. The occurrence tices (Bartley et al., 2014). In the coastal zone, geochemical signa- of some of the highest soil Cu contents in forested parkland where tures can provide insights about land-based sediment and motor vehicle operation is minimal underscores the geologic con- contaminant transport in relation to ecological communities, in- trol of this trace metal. The three-fold variation of terrestrial soil formation that is important for restoration efforts. and sediment CuFe contents reflects the natural variability of Cu in tholeiitic and alkalic olivine basalt lavas. 6. Conclusions A stronger anthropogenic influence was found for fine sediment Cd and Zn contents in the Maunalua watershed, as evidenced by the The health and resilience of coral reefs has been shown to occurrence of maximum values near the urban center in the east improve when land-based sediment and contaminant runoff is 78 R.K. Takesue, C.D. Storlazzi / Estuarine, Coastal and Shelf Science 188 (2017) 69e80 reduced. Insights about runoff sources gained from sediment 160e162. http://dx.doi.org/10.1007/BF01204458. geochemical signatures would make runoff management more Clague, D.A., Frey, F.A., 1982. Petrology and trace element geochemistry of the Honolulu Volcanics, Oahu: implications for the oceanic mantle below Hawaii. effective by targeting priority contributing areas. A sediment J. Petrol 23, 447e504. http://dx.doi.org/10.1093/petrology/23.3.447. geochemical signature based on (La/Yb)N showed a clear distinction Coles, S.L., DeFelice, R.C., Eldredge, L.G., 2002. Nonindigenous Marine Species at between the two main rock types, tholeiitic and alkalic olivine Waikiki and Hawaii Kai, Oahu, Hawaii. Bishop Musuem Technical Report No. 25, Honolulu, Hawaii. 255 pp.. basalt, in the contributing area of a coral-reef fringed embayment in Collins, A.L., Walling, D.E., 2002. Selecting fingerprint properties for discriminating southeast Oahu, Hawaii. Asian dust was only a minor component of potential suspended sediment sources in river basins. J. Hydrol. 261, 218e244. soil and sediment in southeast Oahu. One of five tholeiitic basins in http://dx.doi.org/10.1016/S0022-1694(02)00011-2. Daughton, C.G., Ternes, T.A., 1999. Pharmaceuticals and personal care products in the overall contributing area could be individually distinguished the environment: agents of subtle change? Environ. Health Perspect. 107 (6), using (La/Yb)N and NbFe. Fine sediment from this small sub-basin 907e938. http://www.justor.org/stable/3434573. was dispersed over a large section of the reef flat, illustrating De Carlo, E.H., Anthony, S.S., 2002. Spatial and temporal variability of trace element concentrations in an urban subtropical watershed, Honolulu. Hawaii. Appl. how hydrodynamic sediment transport can increase the scale of Geochem 17, 475e492. http://dx.doi.org/10.1016/S0883-2927(01)00114-7. runoff impacts in coastal areas, and how geochemical signatures De Carlo, E.H., Beltran, V.L., Tomlinson, M.S., 2004. Composition of water and sus- can identify potential target areas for runoff control. Anthropogenic pended sediment in streams of urbanized subtropical watersheds in Hawaii. e metals Cd, Pb, and Zn were associated with fine sediment from the Appl. Geochem 19, 1011 1037. http://dx.doi.org/10.1016/0025-326X(95)98338- W. most heavily urbanized portion of the watershed, and were De Carlo, E.H., Spencer, K.J., 1995. Sedimentary records of anthropogenic inputs of entrained in runoff to the reef flat. Other trace metals Cr, Cu, and Ni heavy metals to the Ala Wai a small man-made estuary in Honolulu. Hawaii. e were associated with volcanic minerals in Hawaiian soil and sedi- pac. Sci. 49, 471 491. De Carlo, E.H., Tomlinson, M.S., Anthony, S.S., 2005. Trace elements in streambed ment, so were unlikely to pose ecological risks. Sediment-source sediments of small subtropical streams on O'ahu, Hawai'i: results from the attribution with geochemical signatures can be a valuable tool for USGS NAWQA program. Appl. Geochem 20, 2157e2188. http://dx.doi.org/ the implementation of effective runoff control measures, a strategy 10.1016/0025-326X(95)98338-W. Done, T., 1995. Remediation of degraded coral reefs: the need for a broad focus. Mar. that has become critical for the protection and restoration of coastal Pollut. Bull. 30, 686e688. http://dx.doi.org/10.1016/0025-326X(95)98338-W. ecosystems as urbanization and exploitation of coastal areas in- Draut, A.E., Bothner, M.H., Field, M.E., Reynolds, R.L., Cochran, S.A., Logan, J.B., crease worldwide. Storlazzi, C.D., Berg, C.J., 2009. Supply and dispersal of flood sediment from a steep, tropical watershed: Hanalei Bay, Kaua'i, Hawai'i, USA. Geo. Soc. Am. Bull. 121, 574e585. http://dx.doi.org/10.1130/B26367.1. Acknowledgments Dymond, J., Biscaye, P.E., Rex, R.W., 1974. Eolian origin of mica in Hawaiian soils. Geo. Soc. Am. Bull. 85, 37e40. http://dx.doi.org/10.1130/0016-7606(1974) 85<37:EOOMIH>2.0.CO;2. This work was funded by and carried out as part of the USGS Esslemont, G., 2000. Heavy metals in seawater, marine sediments and corals from Pacific Coral Reef Project, an effort in the U.S. and its trust territories the townsville section, great barrier reef marine park. Qld. Mar. Chem. 71, e to better understand the effect of geologic processes on coral reef 215 231. http://dx.doi.org/10.1016/S0304-4203(00)00050-5. Fabricius, K.E., 2005. Effects of terrestrial runoff on the ecology of corals and coral systems. Funding was provided by the USGS Coastal and Marine reefs: review and synthesis. Mar. Pollut. Bull. 50, 125e146. http://dx.doi.org/ Geology Program. The authors thank Cheryl Hapke (USGS) for 10.1016/j.marpolbul.2004.11.028. laying the groundwork for the field study; Alyssa Miller and Jolie Ferrat, M., Weiss, D.J., Strekopytov, S., Dong, S., Chen, H., Najorka, J., Sun, Y., Gupta, S., Tada, R., Sinha, R., 2011. Improved provenance tracing of Asian dust sources Wagner, on staff at the time of this study with Malama Maunalua, using rare earth elements and selected trace elements for paleomonsoon for guidance and logistical support during upland and marina studies on the eastern Tibetian Plateau. Geochim. Cosmochim. Acta 75, sediment collection; Leticia Hallas for assistance processing sam- 6374e6399. http://dx.doi.org/10.1016/j.gca.2011.08.025. Frey, F.A., Garcia, M.O., Roden, M.F., 1994. Geochemical characteristics of Koolau ples; and Rob Franks of the University of California at Santa Cruz for Volcano: implications of intershield geochemical differences among Hawaiian assistance with ICP-MS analyses. John Crusius (USGS) and anony- volcanoes. Geochim. Cosmochim. Acta 58, 1441e1462. http://dx.doi.org/ mous reviewers provided excellent suggestions than improved the 10.1016/0016-7037(94)90548-7. fi Frissell, C.A., Liss, W.J., Warren, C.E., Hurley, M.D., 1986. A heirarchial framework for manuscript. Any use of trade, product, or rm names is for stream habitat classification: viewing streams in a watershed context. Environ. descriptive purposes only and does not imply endorsement by the Manage 10, 199e214. http://dx.doi.org/10.1007/BF01867358. U.S. Government. Giambelluca, T.W., Chen, Q., Frazier, A.G., Price, J.P., Chen, Y.-L., Chu, P.-S., Eischeid, J.K., Delparte, D.M., 2013. Online rainfall atlas of Hawai'i. B. Am. Meteorol. Soc. 94, 313e316. http://dx.doi.org/10.1175/BAMS-D-11-00228.1. Appendix A. Supplementary data GPA, 2006. Protecting Coastal and Marine Environments from Impacts of Land- based Activities: a Guide for National Action. United Nations Environment Programme (UNEP) Global Programme of Action for the Protection of the Ma- Supplementary data related to this article can be found at http:// rine Environment from Land-Based Activities (GPA), The Hague, p. 107. dx.doi.org/10.1016/j.ecss.2017.02.013. Gramlich, J.W., Lewis, V.A., Naughton, J.J., 1971. Potassium-argon dating of Holocene basalts of the Honolulu volcanic series. Geo. Soc. Am. Bull. 82 http://dx.doi.org/ 10.1130/0016-7606(1971)82[1399:PDOHBO]2.0.CO;2, 1,399-391,404. References Guzman, H.M., Jimenez, C.E., 1992. Contamination of coral reefs by heavy metals along the Caribbean coast of Central America. Mar. Pollut. Bull. 24, 554e561. Alloway, B.J., 1995. Heavy Metals in Soils. Blackie Academic & Professional, London. http://dx.doi.org/10.1016/0025-326X(92)90708-E. Araújo, M.F., Jouanneau, J.-M., Valerio, P., Barbosa, T., Gouveia, A., Weber, O., Haskins, E.H., Garcia, M.O., 2004. Scientific drilling reveals geochemical heteroge- Oliveira, A., Rodrigues, A., Dias, J.M.A., 2002. Geochemical tracers of northern neity within the Ko'olau shield, Hawai'i. Contrib. Mineral. Petrol 147, 162e188. Portuguese estuarine sediments on the shelf. Prog. Ocean. 52, 277e297. http:// http://dx.doi.org/10.1007/s00410-003-0546-y. dx.doi.org/10.1016/j.csr.2004.02.010. Hilton, R.G., Galy, A., Hovius, N., Chen, M.-C., Horng, M.-J., Chen, H., 2008. Tropical- Barnes, K.K., Kolpin, D.W., Furlong, E.t., Zaugg, S.D., Meyer, M.T., Barber, L.B., 2008. cyclone-driven erosion of the terrestrial biosphere from mountains. Nat. Geosci. A national reconaissance of pharmaceuticals and other organic wastewater 1, 759e762. http://dx.doi.org/10.1038/ngeo333. compounds in the United States - I). Groundw. Sci. Total Environ. 402, 192e200. HOP, 2013. Hawaii Ocean Resources Management Plan. Hawaii Office of Planning. http://dx.doi.org/10.1016/j.scitotenv.2008.04.028. Coastal Zone Management Program, Honolulu, p. 152. Bartley, R., Bainbridge, Z.T., Lewis, S.E., Kroon, F.J., Wilkinson, S.N., Brodie, J.E., Huang, S., Frey, F.A., 2005. Recycled oceanic crust in the Hawaiian Plume: evidence Silburn, D.M., 2014. Relating sediment impacts on coral reefs to watershed from temporal geochemical variations within the Koolau Shield. Contrib. Min- sources, processes, and management: a review. Sci. Total Environ. 468e469, eral. Petrol 149, 556e575. http://dx.doi.org/10.1007/s00410-005-0664-9. 1138e1153. http://dx.doi.org/10.1016/j.scitotenv.2013.09.030. Hughes, T.P., Graham, N.A.J., Jackson, J.B.C., Mumby, P.J., Steneck, R.S., 2010. Rising to Brasher, A.M.D., Wolff, R.H., 2004. Relations between land use and organochlorine the challenge of sustaining coral reef resilience. Trends Ecol. Evol. 25, 633e642. pesticides, PCBs, and semi-volatile organic compounds in streambed sediment http://dx.doi.org/10.1016/j.tree.2010.07.011. and fish on the Island of Oahu, Hawaii. Archive Environ. Contam. Toxicol. 46, Hunter, C.L., Evans, C.W., 1995. Coral reefs in Kaneohe Bay, Hawaii: two centuries of 385e398. http://dx.doi.org/10.1007/s00244-003-3019-4. western influence and two decades of data. B. Mar. Sci. 57, 501e515. Brunskill, G.J., Woolfe, K.J., Zagorskis, I., 1995. Distribution of riverine sediment Jackson, J.B.C., 2008. Ecological extinction and evolution in the brave new ocean. chemistry on the shelf, slope and rise of the Gulf of Papua. Geo-Mar. Lett. 15, P. Natl. Acad. Sci. U. S. A. 105, 11458e11465. http://dx.doi.org/10.1073/ R.K. Takesue, C.D. Storlazzi / Estuarine, Coastal and Shelf Science 188 (2017) 69e80 79
pnas.0802812105. 559e569. http://dx.doi.org/10.1007/s00338-004-0415-9. Jackson, J.B.C., Kirby, M.X., Berger, W.H., Bjorndal, K.A., Botsford, L.W., Bourque, B.J., Oki, D.S., Brasher, A.M.D., 2003. Environmental Setting and the Effects of Natural Bradbury, R.H., Cooke, R., Erlandson, J., Estes, J.A., Hughes, T.P., Kidwell, S., and Human-related Factors on Water Quality and Aquatic Biota, Oahu, Hawaii. Lange, C.B., Lenihan, H.S., Pandolfi, J.M., Peterson, C.H., Steneck, R.S., Tegner, M.J., U.S. Geological Survey Water-Resources Investigations Report 03-4156, p. 98. Warner, R.R., 2001. Historical overfishing and the recent collapse of coastal Oze, C., Fendorf, S., Bird, D.K., Coleman, R.G., 2004. Chromium geochemistry of ecosystems. Science 293, 629e638. http://dx.doi.org/10.1126/science.1059199. serpentine soils. Int. Geol. Rev. 46, 97e126. http://dx.doi.org/10.2747/0020- Jackson, M.C., Frey, F.A., Garcia, M.O., Wilmoth, R.A., 1999. Geology and geochem- 6814.46.2.97. istry of basaltic lava flows and dikes from the Trans-Koolau tunnel, Oahu. Pandolfi, J.M., Jackson, J.B.C., Baron, N., Bradbury, R.H., Guzman, H.M., Hughes, T.P., Hawaii. b. Volcanol. 60, 381e401. Kappel, C.V., Micheli, F., Ogden, J.C., Possingham, H.P., Sala, E., 2005. Are U.S. Jackson, M.L., Levelt, T.W.M., Syers, J.K., Rex, R.W., Clayton, R.N., Sherman, G.D., coral reefs on the slippery slope to slime? Science 307, 1725e1726. http:// Uehara, G., 1971. Geomorphological relationships of tropospherically derived dx.doi.org/10.1126/science.1104258. quartz in the soils of the Hawaiian Islands. Soil Sci. Soc. Am. J. 35, 515e525. Porder, S., Hilley, G.E., Chadwick, O.A., 2007. Chemical weathering, mass loss, and http://dx.doi.org/10.2136/sssaj1971.03615995003500040015x. dust inputs across a climate by time matrix in the Hawaiian Islands. Earth. Kao, S.J., Milliman, J.D., 2008. Water and sediment discharges from small moun- Planet. Sci. Lett. 258, 414e427. http://dx.doi.org/10.1016/j.epsl.2007.03.047. tainous rivers, Taiwan: the roles of lithology, episodic events, and human ac- Prego, R., Caetano, M., Bernardez, P., Brito, P., Ospina-Alvarez, N., Vale, C., 2012. Rare tivities. J. Geol. 116, 431e448. http://dx.doi.org/10.1086/590921. earth elements in coastal sediment of the northern Galician shelf: influence of Knowlton, N., 2001. The future of coral reefs. P. Natl. Acad. Sci. U.S.A. 98, 5419e5425. geologic features. Cont. Shelf Res. 35, 75e85. http://dx.doi.org/10.1016/ http://dx.doi.org/10.1073/onas.091092998. j.csr.2011.12.010. Kurtz, A.C., Derry, L.A., Chadwick, O.A., 2001. Accretion of Asian dust to Hawaiian Prego, R., Caetano, M., Vale, C., Marmolejo-RodrIguez, J., 2009. Rare earth elements soils: isotopic, elemental, and mineral mass balance. Geochim. Cosmochim. in sediments of the vigo ria, NW iberian peninsula. Cont. Shelf Res. 29, Acta 65, 1971e1983. http://dx.doi.org/10.1016/S0016-7037(01)00575-0. 896e902. http://dx.doi.org/10.1016/j.csr.2009.01.009. Kurtz, A.C., Derry, L.A., Chadwick, O.A., Alfano, M.J., 2000. Refractory element Presley, T.K., Jamison, T.J., 2010. Rainfall, Discharge, and Water-quality Data during mobility in volcanic soils. Geology 28, 683e686. http://dx.doi.org/10.1130/ Stormwater Monitoring, H-1 Storm Drain, Oahu, Hawaii, July 1, 2009 to June 30, 0091-7613(2000)28<683:REMIVS>2.0.CO;2. 2010. U.S. Geological Survey Open-File Report 2010-1161, Reston, Virginia, p. 30. Li, Y.-H., Burkhardt, L., Teraoka, H., 1984. Desorption and coagulation of trace ele- Presto, M.K., Ogston, A.S., Storlazzi, C.D., Field, M.E., 2006. Temporal and spatial ments during estuarine mixing. Geochim. Cosmochim. Acta 48, 1879e1884. variability in the flow and dispersal of suspended-sediment on a fringing reef http://dx.doi.org/10.1016/0016-7037(84)90371-5. flat, Molokai. Hawaii. Estuar. Coast. shelf. Sci. 67, 67e81. http://dx.doi.org/ Long, E.R., MacDonald, D.D., Smith, S.L., Calder, F.D., 1995. Incidence of adverse 10.1016/j.ecss.2005.10.015. biological effects within ranges of chemical concentrations in marine and Presto, M.K., Storlazzi, C.D., Logan, J.B., Reiss, T.E., Rosenberger, K.J., 2012. Coastal estuarine sediments. Environ. Manage 19, 81e97. http://dx.doi.org/10.1007/ Circulation and Potential Coral Larval Dispersal in Maunalua Bay, Oahu, Hawaii. BF02472006. U.S. Geological Survey Open-File report 2012-1040, p. 67. Lotze, H.K., Lenihan, H.S., Bourque, B.J., Bradbury, R.H., Cooke, R.G., Kay, M.C., Restrepo, J.D., Park, E., Aquino, S., Latrubesse, E.M., 2016. Coral reefs chronically Kidwell, S.M., Kirby, M.X., Peterson, C.H., Jackson, J.B.C., 2006. Depletion, exposed to river sediment plumes in the southwest Caribbean: rosario Islands, degradation, and recovery potential of estuaries and coastal seas. Science 312, Colomiba. Sci. Total Environ. 553, 316e329. http://dx.doi.org/10.1016/ 1806e1808. http://dx.doi.org/10.1126/science.1128035. j.scitotenv.2016.02.140. Lyons, W.B., Nezat, C.A., Carey, A.E., Hicks, D.M., 2002. Organic carbon fluxes to the Restrepo, J.D., Zapata, P., Diaz, J.M., Garzon-Ferreira, J., Garcia, C.B., 2006. Fluvial ocean from high-standing islands. Geology 30, 443e446. http://dx.doi.org/ fluxes into the Caribbean Sea and their impact on coastal ecosystems: the 10.1130/0091-7613(2002)030<0443:OCFTTO>2.0.CO;2. Magdalena River, Colombia. Glob. Planet. Change 50, 33e49. http://dx.doi.org/ MacDonald, D.D., Ingersoll, C.G., Berger, T.A., 2000. Development and evaluation of 10.1016/j.gloplacha.2005.09.002. consensus-based sediment quality guidelines for freshwater ecosystems. Arch. Rex, R.W., Syers, J.K., Jackson, M.L., Clayton, R.N., 1969. Eolian origin of quartz in soils Environ. Contam. Toxicol. 39, 20e31. http://dx.doi.org/10.1007/s002440010075. of Hawaiian Islands and in Pacific pelagic sediments. Science 163, 277e279. Marsaglia, K.M., 1993. Basaltic island sand provenance. Geological Society of http://dx.doi.org/10.1126/science.163.3864.277. America Special Papers 284, Boulder, CO. In: Johnsson, M.J., Basu, A. (Eds.), Richmond, R.H., Rongo, T., Golbuu, Y., Victor, S., Idechong, N., Davis, G., Kosta, W., Processes Controlling the Composition of Clastic Sediments, pp. 41e65. Neth, L., Hamnett, M., Wolanski, E., 2007. Watersheds and coral reefs: conser- McLennan, S.M., 1989. Rare earth elements in sedimentary rocks: influence of vation science, policy, and implementation. BioScience 57, 598e607. http:// provenance and sedimentary processes. In: Lipin, B.R., McKay, G.A. (Eds.), Rare dx.doi.org/10.1641/B570710. Earth Elements. Reviews in Mineralogy, Washington D.C, pp. 169e200. Roden, M.F., Frey, F.A., Clague, D.A., 1984. Geochemistry of tholeiitic and alkalic lavas McLennan, S.M., Hemming, S., McDaniel, D.K., Hanson, G.N., 1993. Geochemical from the Koolau Range, Oahu, Hawaii: implications for Hawaiian volcanism. approaches to sedimentation, provenance, and tectonics. Geological Society of Earth. Planet. Sci. Lett. 69, 141e158. http://dx.doi.org/10.1016/0012-821X(84) America Special Paper 284, Boulder. In: Johnsson, M.J., Basu, A. (Eds.), Processes 90079-7. Controlling the Composition of Clastic Sediment, pp. 21e40. Rogers, C.S., 1990. Responses of coral reefs and reef organisms to sedimentation. Mee, L., 2012. Between the devil and the deep blue sea: the coastal zone in an era of Mar. Ecol-Prog. Ser. 62, 185e202. http://dx.doi.org/10.3354/meps062185. globalization. Estuarine. Coast. Shelf Sci. 96, 1e8. http://dx.doi.org/10.1016/ Roussiez, V., Aubert, D., Heussner, S., 2013. Continental sources of particles escaping j.ecss.2010.02.013. the Gulf of Lion evidenced by rare earth elements: Floods vs. normal conditions. Milliman, J.D., Farnsworth, K.L., Albertin, C.S., 1999. Flux and fate of fluvial sedi- Mar. Chem. 153, 31e38. http://dx.doi.org/10.1016/j.marchem.2013.04.010. ments leaving large islands in the East Indies. J. Sea Res. 41, 97e107. http:// Schwarzenbach, R.P., Escher, B.I., Fenner, K., Hofstetter, T.B., Johnson, C.A., von dx.doi.org/10.1016/S1385-1101(98)00040-9. Gunten, U., Wehrli, B., 2007. The challenge of micropollutants in aquatic sys- Milliman, J.D., Syvitski, J.P.M., 1992. Geomorphic/tectonic control of sediment tems. Science 313, 1072e1077. http://dx.doi.org/10.1126/science.1127291. discharge to the ocean: the importance of small mountainous rivers. J. Geol. Sholkovitz, E.R., Elderfield, H., Szymcak, R., Casey, K., 1999. Island weathering: river 100, 525e544. http://www.jstor.org/stable/30068527. sources of rare earth elements tot he Western Pacific Ocean. Mar. Chem. 68, Moberly, R., Baver, L.D., Morrison, A., 1965. Source and variation of Hawaiian littoral 39e57. http://dx.doi.org/10.1016/S0304-4203(99)00064-X. sand. J. Sediment. Petrol 35, 589e598. Short, F.T., Wyllie-Echeverria, S., 1996. Natural and human-induced disturbance of Nakai, S., Halliday, A.N., Rea, D.K., 1993. Provenance of dust in the Pacific Ocean. seagrasses. Environ. Conserv. 23, 17e27. http://dx.doi.org/10.1017/ Earth. Planet. Sci. Lett. 119, 143e157. http://dx.doi.org/10.1016/0012-821X(93) S0376892900038212. 90012-X. Smith, J., Douglas, G.B., Radke, L.C., Palmer, M., Brooke, B.P., 2008. Fitzroy River Nelson, S.T., Tingey, D.G., Selck, B., 2013. The denudation of ocean islands by ground Basin, Queensland, Australia. III. Identification of sediment sources in the and surface waters: the effects of climate, soil thickness, and water contact coastal zone. Environ. Chem. 5, 231e242. http://dx.doi.org/10.1071/EN07094. times on Oahu, Hawaii. Geochim. Cosmochim. Acta 103, 276e294. http:// Stearns, H.T., Vaksvik, K.N., 1935. Geology and Ground-water Resources of the Island dx.doi.org/10.1016/j.gca.2012.09.046. of Oahu, Hawaii. Hawaii Division of Hydrography. Bulletin 1, Honolulu, p. 479. Newton, A., Carruthers, T.J.B., Icely, J., 2012. The coastal syndromes and hotspots on Steffen, W., Persson, A., Deutsch, L., Zalasiewicz, J., Williams, M., Richardson, K., the coast. Estuar. Coast. Shelf. Sci. 96, 39e47. http://dx.doi.org/10.1016/ Crumley, C., Crutzen, P., Folke, C., Gordon, L., Molina, M., Ramanathan, V., j.ecss.2011.07.012. Rockstrom,€ J., Scheffer, M., Schellnhuber, H.J., Svedin, U., 2011. The Anthro- Nittrouer, C.A., Brunskill, G.J., Figueiredo, A.G., 1995. Importance of tropical coastal pocene: from global change to planetary stewardship. Ambio 40, 739e761. environments. Geo-Mar. Lett. 15, 121e126. http://dx.doi.org/10.1007/ http://dx.doi.org/10.1007/s13280-011-0185-x. BF01204452. Storlazzi, C.D., Ogston, A.S., Bothner, M.H., Field, M.E., Presto, M.K., 2004. Wave- and Nurse, L.A., McLean, R.F., Agard, J., Briguglio, L.P., Duvat-Magnan, V., Pelesikoti, N., tidally-driven flow and sediment flux across a fringing coral reef: southern Tompkins, E., Webb, A., 2014. Small islands. In: Barros, V.R., Field, C.B., Molokai. Hawaii. Cont. Shelf Res. 24, 1397e1419. http://dx.doi.org/10.1016/ Dokken, D.J., Mastrandrea, M.D., Mach, K.J., Bilir, T.E., Chatterjee, M., Ebi, K.L., j.csr.2004.02.010. Estrada, Y.O., Genova, R.C., Girma, B., Kissel, E.S., Levy, A.N., MacCracken, S., Storlazzi, C.D., Presto, M.K., Logan, J.B., Field, M.E., 2010. Coastal Circulation and Mastrandrea, P.R., White, L.L. (Eds.), Climate Change 2014: Impacts, Adaptation, Sediment Dynamics in Maunalua Bay, Oahu, Hawaii. U.S. Geological Survey and Vulnerability. Part B: Regional Aspects. Contribution of Working Group II to Open-File Report 2010-1217, p. 59. the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. Sturesson, U., 1976. Lead enrichment in shells of Mytilus edulis. Ambio 5, 253e256. Cambridge University Press, Cambridge, UK, pp. 1613e1654. Sutherland, R.A., 2000. Bed sediment-associated trace metals in an urban stream, Ogston, A.S., Storlazzi, C.D., Field, M.E., Presto, M.K., 2004. Sediment resuspension Oahu. Hawaii. Environ. Geol. 39, 611e627. http://dx.doi.org/10.1007/ and transport patterns on a fringing reef flat, Molokai. Hawaii. Coral Reefs 23, s002540050473. 80 R.K. Takesue, C.D. Storlazzi / Estuarine, Coastal and Shelf Science 188 (2017) 69e80
Sutherland, R.A., Tolosa, C.A., 2000. Multi-element analysis of road-deposited GSA Today 7, 1e8. sediment in an urban drainage basin, Honolulu. Hawaii. Environ. Pollut. 110, Waycott, M., Duarte, C.M., Carruthers, T.J.B., Orth, R.J., Dennison, W.C., Olyarnik, S., 483e495. http://dx.doi.org/10.1016/S0269-7491(99)00311-5. Calladine, A., Fourqurean, J.W., Heck Jr., K.L., Hughes, A.R., Kendrick, G.A., Talbot, V., Chegwidden, A., 1983. Heavy metals in the sediments of Cockburn Sound, Kenworthy, W.J., Short, F.T., Williams, S.L., 2009. Accelerating loss of seagrasses Western Australia, and its surrounding areas. Environ. Pollut. B 5, 187e205. across the globe threatens coastal ecosystems. P. Natl. Acad. Sci. U. S. A. 106, http://dx.doi.org/10.1016/0143-148X(83)90003-4. 12377e12381. http://dx.doi.org/10.1073/pnas.0905620106. Tomlinson, M.S., De Carlo, E.H., 2003. The need for high resolution time series data Wentworth, C.K., Winchell, H., 1947. Koolau basalt series, Oahu, Hawaii. Geo. Soc. to characterize Hawaiian streams. J. Am. Water Resour. As 39, 113e123. http:// Am. Bull. 58, 49e78. http://dx.doi.org/10.1130/0016-7606(1947)58[49:KBSOH] dx.doi.org/10.1111/j.1752-1688.2003.tb01565.x. 2.0.CO;2. Turekian, K.K., 1977. The fate of metals in the coeans. Geochim. Cosmochim. Acta 41, Winchell, H., 1947. Honolulu series, Oahu, Hawaii. Geo. Soc. Am. Bull. 58, 1e48. 1139e1144. http://dx.doi.org/10.1016/0016-7037(77)90109-0. http://dx.doi.org/10.1130/0016-7606(1947)58[1:HSOH]2.0.CO;2. USEPA, 1996. Method 3052: Microwave Assisted Acid Digestion of Siliceous and Windom, H.L., Schropp, S.J., Calder, F.D., Ryan, J.D., Smith, R.G., Burnery, L.C., Organically Based Matrices, Part of the Test Methods for Evaluating Solid Waste, Lewis, F.G., Rawlinson, C.H., 1989. Natural trace metal concentrations in estu- Physical/chemical Methods. U.S. Environmental Protection Agency, p. 20. arine and coastal marine sediments of the southeastern United States. Environ. https://www.epa.gov/hw-sw846/sw-846-test-method-3052-microwave-assis- Sci. Technol. 23, 314e320. http://dx.doi.org/10.1021/es00180a008. ted-acid-digestion-siliceous-and-organically-based. Wolanski, E., Martinez, J.A., Richmond, R.H., 2009. Quantifying the impact of Vitousek, P.M., Chadwick, O.A., Crews, T.E., Fownes, J.H., Hendricks, D.M., watershed urbanization on a coral reef: Maunalua Bay. Hawaii. Estuar. Coast. Herbert, D., 1997. Soil and ecosystem development across the Hawaiian islands. shelf. Sci. 84, 259e268. http://dx.doi.org/10.1021/es00180a008.
View publication stats Marine Pollution Bulletin 71 (2013) 92–100
Contents lists available at SciVerse ScienceDirect
Marine Pollution Bulletin
journal homepage: www.elsevier.com/locate/marpolbul
Can stormwater be detected by algae in an urban reef in Hawai‘i? ⇑ T. Erin Cox a, , Celia M. Smith a, Brian N. Popp b, Michael S. Foster c, Isabella A. Abbott a a University of Hawai‘i, Department of Botany, 3190 Maile Way, Room 101, Honolulu, HI 96822, USA b University of Hawai‘i, Department of Geology and Geophysics, 1680 East–West Road, Honolulu, HI 96822, USA c Moss Landing Marine Laboratories, 8272 Moss Landing Road, Moss Landing, CA 95039, USA article info abstract
Keywords: Nitrogen (N) enrichment of tropical reefs can result in the dominance of invasive algae. The invasive alga ‘Ewa Beach Acanthophora spicifera and the native alga Laurencia nidifica are part of a diverse reef assemblage in ‘Ewa O‘ahu Beach, O‘ahu. Their N contents and d15N values were investigated to determine if N was enriched and to Acanthophora spicifera evaluate potential nitrogenous sources near and removed from storm-drain outlets. d15N values of algae Laurencia nidifica (3.8–17.7‰) were within and above the range for algae around the island (1.9–11.9‰). Elevated algae N Nitrogen isotope isotope values (d15N>+7‰, [N] > 1.6%) and seawater nitrate + nitrite levels (0.59–7.93 M) indicated a Effluent l 15 Nutrient inputs mixed, high nutrient environment. The overlap in d N values with multiple nitrogenous sources pre- cluded identification. However, spatial and temporal patterns did not support stormwater as the domi- nant, nitrogenous source. Patterns were congruent with algal incorporation of terrestrial derived N, subjected to a high degree of biogeochemical cycling. Ó 2013 Elsevier Ltd. All rights reserved.
1. Introduction on algal communities and near shore water quality (Bernardo, 2008). The land in the ‘Ewa Beach area, once used for extensive su- Storm-drains have the potential to act as conduits during rain gar agriculture, is currently being developed into neighborhood events, collecting and focusing nitrogen (N) enriched runoff into subdivisions (Schaefers, 2006). Although climate in ‘Ewa Beach is coastal systems. In contrast, regions with histories of intensive tropical and dry, episodic rainfall events in winter months often re- agricultural production and/or use of septic tanks in high density sults in coastal low-land flooding. The average precipitation be- housing areas could lead to elevated background nutrient levels tween 1949 and 2001 was 508 mm but, 381 mm fell during in groundwater that obscure episodic rain events. Very little is October–March (Otkin and Martin, 2004). To prevent flooding known about the fate of these effluents into oligotrophic waters. and continue neighborhood development, new storm drains have Increased N input into aquatic systems can cause increased macro- been proposed (Schaefers, 2006; Bernardo, 2008) that would dis- algal production, changes in assemblages of organisms, alteration charge into nearby intertidal habitat. The shallow coastal area of food webs, and disruption of nutrient biogeochemical cycling along southwest O‘ahu has historically been an area of high algal (Valiela et al., 1992; Walsh, 2000; Cole, 2003; Choi et al., 2007). abundance and is culturally important to Hawaiians for collection In tropical pristine waters, where nutrient levels are low, addition of edible limu (macroalgae) (Abbott, 1996; Leone, 2004; Ohira, of nutrients can potentially be detrimental. There are several 2005). Longtime residents familiar with ‘Ewa Beach coastline recall examples of blooms of non-indigenous, nuisance seaweeds that the ability to collect burlap bags full of edible macroalgae (Abbott, out-compete native algae and corals in areas adjacent to sewage 1996). The recent macroalgal community appears to be different inputs (Banner, 1974; Russell, 1992; Stimson et al., 2001, 2002, from past descriptions. The local perception is the edible, native 2004, 2005). Also, increase in nutrient supply to oligotrophic marine plants are in decline (Leone, 2004). Terrestrially-derived waters is often cited as a catalyst for community phase shifts from nutrients potentially from past agricultural and waste manage- coral to algal dominated reefs (McCook, 1999). ment may explain the wide availability of limu in the ‘Ewa Beach In ‘Ewa Beach, a growing coastal community on the island of area in the past and the coincident reef algal decline with recent O‘ahu, a primary concern is the impact(s) of storm-drain effluent changes in land use. An investigation (Lapointe and Bedford, 2011) along ‘Ewa Beach coastline used N in algae to identify nitrog- ⇑ Corresponding author. Address: Univ. Pierre et Marie Curie-Paris, Laboratoire enous contributions from a storm-drain. They concluded that d’Oceanographie de Villefranche (UMR7093), 06234 Villefranche-sur-Mer Cedex, storm-drain effluent favors the growth of non-indigenous algal France. Tel.: +33 (0) 4 93 76 38 33; fax: +33 (0) 4 93 76 38 34. species along this shore but did not thoroughly consider all poten- E-mail addresses: [email protected], [email protected] (T. Erin Cox), celia@ tial sources of N or the effects of variation in habitat on algal abun- hawaii.edu (C.M. Smith), [email protected] (B.N. Popp), [email protected] dance and species composition. (M.S. Foster).
0025-326X/$ - see front matter Ó 2013 Elsevier Ltd. All rights reserved. http://dx.doi.org/10.1016/j.marpolbul.2013.03.030 T. Erin Cox et al. / Marine Pollution Bulletin 71 (2013) 92–100 93
Identifying nutrient sources and fates can help predict potential hypothesis of Lapointe and Bedford (2011) that N in storm-drain impacts from storm-drain discharge on water quality and reef algal effluent supports the physiology of introduced and native algae. physiology and abundance. Reef algae are generally N-limited so if We sampled the N content and the d15N values of algae along a gra- a particular N source is dominant, reef algae are well suited for dient leading offshore, in the dry and rainy season, at four sites tracing N in aquatic environments. Many algal species have simple with and six sites without storm-drains that discharge into the morphologies that allow uptake of nutrients from the water col- intertidal zone. Isotopic values in algae are compared to the con- umn and integration into growth or storage (Larned, 1998; centrations and d15N values (nitrate + nitrite) of potential N McCook, 1999; Schaffelke, 1999; Fong et al., 2001, 2003; Lapointe sources. Support for a ‘‘drain effect’’ would be indicated by a spatial et al., 2004; Lin and Fong, 2008; Teichberg et al., 2008; Dailer et al., effect; reef algae with enriched N concentration and elevated d15N 2010, Dailer et al., 2012b). values in near-shore macroalgae at sites with a storm-drain and The N isotopic composition of macroalgae can also be used to after large rain events. Finally, to gain a wider perspective on water identify and track sources of N. The d15N values of macroalgae vary quality and reef physiology, we compare d15N values of algae col- as a function of N sources, the biogeochemical history of the N uti- lected at ‘Ewa Beach to values of algae collected from several other lized, species of N used, and the extent of dissolved N use (Peterson nearshore sites around O‘ahu island. and Fry, 1987; Kendall, 1998; Robinson, 2001). Atmospheric nitro- gen (N ) can be fixed and reduced to ammonium (NHþ, ammonifi- 2 4 2. Methods cation) and converted under low oxygen conditions to nitrite and nitrate (NO and NO , nitrification) by bacteria and archaea. Dur- 2 3 2.1. Description of sites along ‘Ewa Beach ing these transformations the N isotopic composition is altered via kinetic isotope fractionation, or the differential reaction of Ten sites (numbered 1–10 from west to east) were selected along 15N relative to 14N. N fixed in the atmosphere typically has d15N the shore of O‘ahu, from One‘ula Beach Park (21°18036.5900 values of 4‰ to +4.0‰ (Owens, 1987; Macko and Ostrom, N–158°0027.5200W) to ‘Ewa Beach proper (21°18042.1700N–158°00 1994) while dissolved N in groundwater including that derived 15.6600W), to examine the N sources of reef algae and to determine from manure and waste from sewage treatment with increased the extent of N contribution from storm-drain discharge (Fig. 1). bacterial processing has d15N values that can range from about Four sites (sites 4, 6, 8, 9) each contain an existing large storm-drain +6‰ to +22‰ (Macko and Ostrom, 1994; Kendall, 1998). For that directly discharges accumulated urban runoff into the intertidal macroalgae with simple morphologies that completely use the zone. These sites are referred to as ‘‘Drain sites’’. ‘‘Control sites’’ surrounding available pool of nitrate and/or ammonium for growth (sites: 1, 2, 3, 5, 7, 10), were selected for comparison and are inter- the incorporation of 15N occurs with little fractionation (Peterson spersed between sites with large storm-drains’ in areas with similar and Fry, 1987; Gartner et al., 2002; Cohen and Fong, 2005; topography, limestone and sand substrate, and shore direction as Umezawa et al., 2007) and quickly, within 7 days (Gartner et al., the adjacent Drain site. 2002). Hence when a single N source is dominant in the ecosystem, The large existing storm-drains have been in use for over the d15N values of macroalgae can reflect the d15N values of the 30 years (C. Morgan, Planning Solutions, personal communication). source of N. The storm-drains at Drain 4, 6, and 8 drain a watershed of approx- Previous investigations of the d15N values and N content in reef imately 13–31 acres (Hiyakumoto, 2012). Drain 9 is an ocean out- algae have had limited replication of sites of interest, limited sea- let for a larger watershed that extends more than 420 acres (C. sonal data, limited tissue replication, and no investigation of d15N Morgan, Planning Solutions, personal communication). values of nitrogenous sources (Derse et al., 2007; Lapointe and The ‘Ewa Plain consists of a Pleistocene age highly porous lime- Bedford, 2011). As such these earlier studies fail to take into ac- stone caprock, a wedge of sediment overlying volcanic rock (Mink, count the complexity of nutrient dynamics and species specific as- 1989). Highly rugose carbonate platforms occur from the intertidal pects of plant metabolism. Identifying sources of N can be zone into the shallow subtidal zones and are dominated by algae. challenging because the isotopic compositions of nitrogenous The urbanized area surrounding these sites is served by sewer sys- nutrients in land and marine environments is dynamic. N not only tems, but some households use on-site disposal of wastewater de- varies in form but increased N inputs result from upwelling of off- fined as either soil treatment, seepage pits, or cesspools (Whittier shore seawater (nitrates), discharge of wastewater (nitrates and and El-Kadi, 2009). ammonium), fertilizer runoff (nitrates and ammonium) and natu- ral sources of N in runoff. Mixing of these multiple sources occurs under field conditions. Mixing and inputs can also vary temporally 2.2. Species selection for fine scale examination with seasonal changes in the environment and with rainfall. For in- stance, groundwater discharge can vary with rainfall amounts and Two closely related reef algae, Acanthophora spicifera and Lau- deliver N into coastal waters (Lapointe et al., 2004). Temporal as- rencia nidifica (Order Ceramiales, Family Rhodomelaceae) were pects could be particularly important when assessing impact from sampled because they are abundant and have a simple morphology 15 a pulsed N source like a storm-drain. Furthermore, biogeochemical that should allow d N values of these reef algae to quickly reach reactions can occur naturally in groundwater and can be associated steady state with dissolved nitrogenous nutrients. A. spicifera,an with septic systems. Thus, a higher d15N value in algae does not introduced species, in Hawai‘i (Doty, 1961; Russell, 1992), is necessarily mean they have been exposed to elevated N inputs known to grow faster in N enriched waters and has been previ- from drainage effluent. Despite these challenges several studies ously used to identify anthropogenic and natural nitrogenous 2 have successfully used bulk d15N values of macroalgae to identify sources (Lin and Fong, 2008). For all sampling, small (1–2 cm ) dis- land-based N inputs into pristine coastal zones where the differ- crete clumps consisting of several whole plants were collected. ences between end points of nutrient gradients are well resolved Laurencia nidifica was rare at sites 8–10 and few or no individuals (Sammarco et al., 1999; Schaffelke, 1999; Umezawa et al., 2002; were collected. Cohen and Fong, 2005; Garrison et al., 2007; Lin et al., 2007; Teichberg et al., 2008; Dailer et al., 2010, 2012a). 2.3. Sampling methodology We used the N content and the d15N values of common reef al- gae to identify the potential source and fate of N along the coastal Sampling occurred from the upper edge of the intertidal zone or region of ‘Ewa Beach, O‘ahu. Specifically, we rigorously test the storm-drains to a distance of approximately 20 m (1–3 m depth) 94 T. Erin Cox et al. / Marine Pollution Bulletin 71 (2013) 92–100
The HawaiianThe Hawaiian Archipelago Archipelago
N O‘aO‘ahu h u
‘EwaBeach‘Ewa‘E w aBeach B e a c h
Fig. 1. Map of the Hawaiian Islands (top) showing southwest O‘ahu and an aerial image of ‘Ewa Beach (Google EarthÓ, bottom) with selected sites 1–10 (circles are control sites with no drain present, squares are drain sites, i.e. drain present), nearby parks (picnic tables), and golf course (golfer). The and enclosed area marked with are the approximate locations of nearby on-site sewage disposal systems (OSDSs) as identified by Whittier and El-Kadi, 2009. More inland OSDS identified in 2009 within the Hasekeko development (northwest cleared portion of image) have been excluded from this representation as these OSDS are temporary structures and possibly not present when this study was conducted. offshore. One algal sample of each species was collected, when occurred, from 12 shallow water sites around the island of O‘ahu encountered, at set intervals ( 1 m) along a transect line. If both for determination of d15N values in reef algae. Results of analyses species did not occur, the available species was collected. of these specimens provide perspective on whether d15N values Replicate algal samples were collected in months that are typi- from algae within the ‘Ewa Beach area differ from those at other cally dry (July–August) and rainy (November–December) to ac- intertidal shores. Nine sites are located in a different watershed count for seasonality and discharge events. Because of permit from sites 1–10. These watersheds have variable rainfall amounts, requirements, sites 1–6 were sampled in 2007 while sites 7–10 drainage areas, land use, onshore wave action, and nearshore were sampled in 2008. The accumulated rainfall 7 days prior to this depths, and are thus predicted to have different N concentrations collection ranged from 1.09 to 8.66 cm in the rainy season and and varying isotopic compositions (Van Houtan et al., 2010). On trace amounts to 1.04 cm in the dry season (Ewa Kalaeloa Airport the 13th of March 2008 samples (n = 5–7) were collected along Station Id:GHCND:USW00022551). the Wai‘anae (West) coast. On the 9th of October 2011 algae were To determine the nutrient inputs delivered by storm drain flow collected from the six other sites and from Control 2 and Drain 4 and water column nutrients, at time of plant collection in the dry (for a total of 9 sites). Control 2 and Drain 4 were sampled again season, 80 ml water samples (n = 1/site) were collected from on this day in 2011 to account for any temporal variation that near-shore or at the mouth of drains. On a large rain event in may have occurred from 2007. One to two replicate samples of 2007, samples of flowing water were collected from street storm each species were collected when encountered at a site. Samples grates that empty directly into Drain 4, 6, and 8 and were used were from separate, growing plants found within a 100 m dis- as a proxy for storm-drain nutrient source. For comparison, on tance from each other in intertidal habitats. the same rain event, water was collected from the ocean at Control 3 and from a shallow water habitat 100 m to the west of Site 1. In 2.5. Nitrogen content and d15N value determination the dry season in 2008, an 80 ml water sample was collected from standing water in Drain 9. One 80 ml water sample was collected Algal samples were immediately cleaned of epiphytes, rinsed in from a nearby non-potable irrigation well that taps into groundwa- de-ionized water, and placed into an oven at 60 °C until dried. ter as a water source. Dried samples were ground into a fine powder, stored in glass scin- tillation vials until further analyses. Water samples were filtered 2.4. Methodology for island wide d15N value comparison (0.22 lm Millipore filter) and stored frozen ( 20 °C) until analyzed. Reef algae used in similar studies examining N inputs in Hawai‘i Carbon and N isotope compositions of reef algae were deter- (Dailer et al., 2010: A. spicifera, Laurencia mcdermidiae, Astronema mined using an on-line carbon–nitrogen analyzer coupled with breviarticulatum, or Ulva lactuca) were collected, when they an isotope ratio mass spectrometer (Finnigan ConFlo II/Delta-Plus). T. Erin Cox et al. / Marine Pollution Bulletin 71 (2013) 92–100 95
Isotope values are reported in standard d-notation relative to an Lastly, to determine if N inputs in the ‘Ewa Beach area were 15 international standard (V-PDB and atmospheric N2 for carbon similar to other areas in O‘ahu, we spatially mapped the d N and N, respectively). Glycine reference compounds with well-char- values of reef algae collected from around the island, following acterized isotopic compositions were used to ensure accuracy of all methods used by Dailer et al. (2010). isotope measurements. Several samples were measured in dupli- cate or triplicate and the reproducibility associated with these 3. Results measurements was typically 0.2‰ for both carbon and N isotopic measurements. 3.1. Fine scale examination Analyses of water samples for dissolved inorganic nutrients (NHþ,NO ,NO ,PO , SiO ) were performed by the Analytical 4 2 3 4 4 The non-linear relationship (n =7, p = 0.05, r-squared = 0.58) Lab at the Marine Science Institute, University of California, Santa d15 Barbara using a continuous flow technique with a Quick Chem between the ln [nitrate] concentration and N (nitrate + nitrite) values in seawater suggested the d15N values of reef algae in the 800 Flow Injection Analyzer manufactured by Lachat Instruments, Inc. The precision associated with these measurements was typi- ‘Ewa Beach area were the result of a mixing of multiple sources of N with different isotopic compositions (Fig. 2). cally 0.05%. 15 A comparison of nutrient concentrations within the water col- Analyses of water samples for d N (nitrate + nitrite) values were conducted at the University of Washington, Isolab using the umn at sites 1–10 with the concentrations measured in potential sources indicated mixing and dilution of N (Table 1). Silicate con- bacterial denitrifier method (following Sigman et al., 2001) along with an autosampler, PreCon GasBench II assembly coupled to a centrations in groundwater were well above those measured in the seawater at any site. Similarly, the standing water in Drain 9 Finnigan Delta Plus. The isotopic compositions of the international reference materials USGS34 and IAEA-NO-3 were used to ensure was particularly high in silicate. The nutrient concentrations were variable among replicate storm-drain samples and tended towards accuracy. Several samples were measured in duplicate and had a standard deviation from 0.02‰ to 0.18‰. relatively higher phosphate, nitrate, and ammonium concentra- tions than those measured in seawater at nearby sites and from concentrations measured in ocean water collected simultaneously, after the same rain event. 2.6. Analyses The average d15N values of reef algae at sites 1–10 ranged from 4.8‰ to 14.7‰ (Table 2); in comparison the d15N values of nitrate A linear regression was used to examine the relationship be- 15 ‰ ‰ tween the ln [nitrate] and d N (nitrate + nitrite) values of seawa- ranged from 1.5 to 2.4 for open ocean water from 150 m and 7.0–7.1‰ from 500 m (as an estimate of upwelled nitrate) at Sta- ter at sites 1–10. Linearity would indicate a closed system with one source of N while a non-linear relationship indicates an open tion ALOHA located 100 km north of O‘ahu (Casciotti et al., 2008) to 27.9‰ for ‘Ewa Plain groundwater. Reef algae collected during system with mixing of N sources (Robinson, 2001). In addition, the d15N values of reef algae, nearby water sources and the concen- both the dry and rainy season at Drains 4 and 6 had higher mean d15N values (4: 12.0–13.3‰, 6: 10.0–11.9‰) than the d15N values tration of nutrients in the water column and source samples were compared in an attempt to identify N sources for algae. of nitrate+nitrite measured from runoff draining into the associ- 15 ated storm-drains (7.7–7.9‰) but algal isotopic values were much Statistical analyses were performed on total N content and d N 15 values (separately) for each species. Specifically for each species, a lower than the d N values measured for nitrate + nitrite in groundwater (27.9‰). Similarly, the d15N values of reef algae at two-way analysis of variance (ANOVA) was used to compare reef ‰ algal values between drain and control sites in the dry and rainy Drain 9 (5.3–7.9 ) were lower than the isotopic values of nitra- te + nitrite measured for standing water in Drain 9 (19.7‰) and season. Data met the requirements for normality and homogeneity groundwater (27.9‰). However, the reef algae at Drain 8 had sim- of variance. Multiple linear regressions were used to examine the 15 15 ilar mean d N values (8.8–10.4‰) as the storm-water collected relationship between d N values of L. nidifica and A. spicifera and ‰ distance from shore. Each site was examined separately because near Drain 8 (10.0 ). In contrast to reef algae, the two ocean water samples collected after the rain event had d15N values of nitrate + of unique site characteristics. Season and an interaction of Dis- nitrite (22.3‰, 20.2‰) that were similar to the d15N values for tance Season were included in the models as seasonal rainfall groundwater (27.9‰) and the values of water in Drain 9 (19.7‰). could alter delivery of N to the shore. Regressions were not per- formed on %N in algae because of the correlation with d15N values. Because nitrogenous nutrients can be derived from a wide 30 range of sources (e.g., natural runoff, groundwater seeps, submar- y = 13.96 - (3.39 + x) 28 r-squared = 0.58, p > 0.05 ine groundwater discharge, agricultural development, etc.) and 26 have the potential to deliver N to coastal environments on a local scale and large scale, two analyses were conducted. First, to deter- 24 mine if inputs were related to fine-scale location, we statistically 22 examined the geographic affinity of total N content and d15N val- 20
ues for the two reef algae collected at sites 1–10. In these fine-scale N (‰) 18 15 analyses, separate statistical analyses were performed for total N δ 16 15 content and d N values for each species in both seasons. Using 14 Euclidean distance, a similarity matrix was constructed between 12 each value at every site. Second a separate, similarity matrix of geographical distance between sites also was constructed. RELATE 10 in PRIMER-E (Clarke and Warwick, 2001), a Mantel-like test, was 8 -101234 used to examine the statistical relationship between the two types ln NO3- of similarity matrices for both species in both dry and rainy season. RELATE performs Spearman Rank correlation to determine Rho Fig. 2. The non-linear relationship between seawater concentration and d15N with 999 permutations and this is used to determine statistical values (nitrate + nitrite) suggesting an open system with mixing of multiple sources significance. with different d15N values. 96 T. Erin Cox et al. / Marine Pollution Bulletin 71 (2013) 92–100
Table 1 Phosphate, Silicate, Nitrite, Nitrate, and Ammonia as ppm values and d15N values from water samples collected inshore at Sites 1–10, storm-drain grates, and well. n = 1 unless otherwise specified, n > 1 values are expressed as mean ± SE. Values with an are below detectable values, and are samples collected on the same rain event.
Site name Dry season (lM) Rainy season (lM) d15N(‰) dry season 3 SiO NO NO þ NO NH 3 SiO NO NO þ NO NH PO4 2 2 2 3 3 PO4 2 2 2 3 3 Control 1 – – – – – 0.10 7.9 0.20 0.96 1.18 – Control/proposed Drain 2 – – – – – 0.09 12.7 0.38 2.53 2.02 – Control 3 0.15 10.0 0.22 1.03 0.81 0.10 7.0 0.29 1.53 0.94 10.7 Drain 4 0.23 10.9 0.18 2.88 2.71 0.12 6.6 0.29 2.57 1.86 15.9 Control 5 0.15 3.3 0.18 0.72 0.54 0.07 8.2 0.26 1.29 0.34 13.2 Drain 6 0.07 16.5 0.14 7.93 7.79 0.07 13.5 0.37 1.71 0.87 15.2 Control 7 – – – – – – – – – – 9.0 Drain 8 0.08 3.7 0.20 1.02 0.82 0.08 3.7 0.20 0.82 3.10 – Drain 9 0.09 2.6 0.13 0.59 0.46 – – – – - 14.0 Control 10 – – – – – 0.09 12.7 0.38 2.53 2.02 – n 3 SiO NO NO +NO NH n d15N PO4 2 2 2 3 3 Drain water near 4 1 36.2 159.0 6.8 37.4 32.9 1 7.9 Drain water near 6 1 17.5 47.2 5.2 128.0 26.5 1 7.7 Drain water near 8 – – – – – – 1 10.0 Standing water in Drain 9 1 1.4 515.6 3.1 99.2 15.6 1 19.7 Ocean water at 3 1 0.2 4.4 0.4 2.4 3.3 1 22.3 Ocean water at One‘ula 1 – – – – – 1 20.2 Well water, proxy for groundwater 1 2.5 ± 0.1 574.2 ± 20.0 0.8 ± 0.0 23.5 ± 0.4 10.1 ± 0.0 1 27.9 Deep ocean (Casciotti et al., 2008) 7.0–7.1 Open ocean (Casciotti et al., 2008) 2.4–3.0
Table 2 Thalli mean ± SE Total Nitrogen as % dry weight and d15N values for A. spicifera and L. nidifica in dry and rainy seasons at sites 1–10.
Site name Species Dry season Rainy season n Total N (%) d15N n Total N (%) d15N(‰) Control 1 A. spicifera 5 2.6 ± 0.1 14.0 ± 0.5 3 1.6 ± 0.1 8.2 ± 0.7 L. nidifica 2 2.6 ± 0.3 14.7 ± 0.9 3 1.6 ± 0.1 9.5 ± 1.1 Control/proposed drain 2 A. spicifera 6 2.1 ± 0.1 11.6 ± 0.2 1 2.3 10.8 L. nidifica 4 2.1 ± 0.1 12.3 ± 0.1 4 2.1 ± 0.4 11.4 ± 0.3 Control 3 A. spicifera 5 1.8 ± 0.1 6.8 ± 0.1 3 1.8 ± 0.1 7.0 ± 0.0 L. nidifica 4 2.1 ± 0.1 6.9 ± 0.1 4 2.0 ± 0.2 7.3 ± 0.1 Drain 4 A. spicifera 5 2.5 ± 0.2 12.2 ± 0.4 5 2.5 ± 0.1 12.0 ± 0.4 L. nidifica 5 2.6 ± 0.2 13.3 ± 0.7 4 2.1 ± 0.3 12.0 ± 0.7 Control 5 A. spicifera 5 2.1 ± 0.1 9.2 ± 0.2 7 1.9 ± 0.2 8.5 ± 0.1 L. nidifica 7 2.1 ± 0.2 9.5 ± 0.2 6 2.0 ± 0.2 8.7 ± 0.1 Drain 6 A. spicifera 8 2.7 ± 0.1 11.3 ± 0.4 8 2.4 ± 0.1 10.5 ± 0.2 L. nidifica 6 2.5 ± 0.2 11.9 ± 0.2 6 2.5 ± 0.1 11.3 ± 0.3 Control 7 A. spicifera 6 1.8 ± 0.1 11.4 ± 0.4 5 1.7 ± 0.1 12.5 ± 0.7 L. nidifica 4 1.5 ± 0.3 11.0 ± 0.5 – – – Drain 8 A. spicifera 6 1.7 ± 0.2 12.4 ± 1.5 5 2.3 ± 0.1 10.4 ± 1.5 L. nidifica 1 1.3 8.8 – – – Drain 9 A. spicifera 6 1.1 ± 0.1 7.5 ± 0.5 5 1.2 ± 0.0 5.3 ± 0.5 L. nidifica –– – –– – Control 10 A. spicifera 6 1.1 ± 0.0 4.8 ± 0.3 5 1.4 ± 0.1 5.3 ± 0.3 L. nidifica –– – –– –
Although the d15N values and N contents of A. spicifera differed controls, p-value = 0.40, drain x season p-value = 0.53). In addition, among sites, differences were not significant between drain and N contents of L. nidifica did not vary between dry and rainy seasons control sites and did not differ between season (Table 2)(d15N: (Table 2)(p-value = 0.98). Two-Way ANOVA, drain vs control p-value = 0.45, season A. spicifera and L. nidifica at several drain and control sites had p-value = 0.41, drain season p-value = 0.89% N: Two-way ANOVA, higher d15N values nearshore with lower values further from shore. drain vs control p = 0.42, season p = 0.92, drain season p = 0.63). This offshore gradient was seen in both seasons (Fig. 3). For nine of Average values of d15N and% N in L. nidifica from sites 1–8 were 10 sites, distance from shore was a significant factor in a multiple similar to those from A. spicifera at sites 1–10 (Table 2). The d15N regression model used to determine d15N values of A. spicifera values of L. nidifica varied among sites but differences were not sig- (Table 3). Values of d15N tended to be lower in the rainy season nificant (two-way ANOVA) between drain and control sites and for most sites did not alter the slope of the relationship (p-value = 0.21), rainy and dry season (p-value = 0.45), or for the between d15N values and distance from shore. However, an inter- interaction of site type (drain or control) and season action of distance and season was observed in samples of A. spicif- (p-value = 0.33). Although N content of L. nidifica at Drains 4 and era at three sites: Control 1, Control 3, and Control 5 (Table 3). 6 was higher in both seasons, the values at these sites were statis- The d15N values of L. nidifica differed at some sites with distance tically similar to Controls 1, 2, 3, 5, 7 and similar to the sample col- from shore and varied for some sites with season (Fig. 3). Three out lected in the dry season at Drain 8 (Two-way ANOVA, drain vs. of six sites sampled showed significant variation in d15N values of T. Erin Cox et al. / Marine Pollution Bulletin 71 (2013) 92–100 97
20 C1 C/PD2 C3 D4 C5 15
10
5 ) 00 / 0 0
N ( 20
15 D6 C7 D8 D9 C10 δ 15
10
5
0 0 5 10 15 20 25 0 5 10 15 20 25 0 5 10 15 20 25 0 5 10 15 20 25 0 5 10 15 20 25 Distance from shore (meters)
Fig. 3. d15N(‰) values of A. spicifera N and L. nidifica d in the dry (closed) and rainy (open) season as distance from shore (meters) for sites.
Table 3 Multiple linear regressions show the relationship between d15NofA. spicifera and L. nidifica with distance from shore, season, and the interaction of terms at drain (D) and control (C) sites (rel = direction of relationship, NS = not significant at a = 0.05, NA = not applicable).
Site Acanthophora spicifera Laurencia nidifica
Terms rel df Fp-Value R2 ( 100) Terms rel df F p-Value R2 ( 100) C1 Distance, p = 0.010 7 27.8 0.004 95.4 Distance, NS 4 11.3 0.043 79.4 Season, NS – Season, p = 0.043 Distance Season, p = 0.011 + Distance Season, NS C2 Distance, NS 5 0.6 0.488 12.7 Distance, NS 7 1.1 0.433 39.6 Season, NA Season, NS Distance Season, NA Distance Season, NS C3 Distance, p = 0.032 8 20.6 0.085 70.6 Distance, NS 8 9.3 0.021 61.6 Season, NS Season, p = 0.021 + Distance Season, p = 0.041 + Distance Season, NS D4 Distance, p = 0.001 9 22.4 0.001 73.6 Distance, p = 0.002 8 37.1 <0.001 95.7 Season, NS Season, NS Distance Season, NS + Distance Season, p = 0.019 + C5 Distance, p = 0.010 11 6.2 <0.001 85.9 Distance, p = 0.003 12 21.4 <0.001 81.1 Season, p = <0.001 Season, p = 0.045 Distance Season, p = 0.016 + Distance Season, NS D6 Distance, p = 0.017 15 14.3 <0.001 68.8 Distance, p = 0.021 11 9.4 0.021 48.4 Season, p = <0.001 Season, NS Distance Season, NS Distance Season, NS C7 Distance, p = 0.046 10 5.30 0.046 37.4 NA Season, NS Distance Season, NS D8 Distance, p = <0.001 10 18.0 0.001 81.8 NA Season, p = 0.009 Distance Season, NS D9 Distance, p = < 0.040 10 10.8 0.005 73.0 NA Season, p = 0.002 Distance Season, NS C10 Distance, p = 0.044 10 2.7 0.128 53.4 NA Season, NS Distance Season, NS
L. nidifica with distance from the shoreline. Distance was a signifi- related to geographic location for the dry season (RELATE test:, cant predictor of d15N values at Drain 4, Control 5, and Drain 6. For Rho = 18.9%, p-value = 0.02), but not in the rainy season these sites, the d15N values were lower with distance from the (Rho = 3.5%, p-value = 0.20). shoreline. Season was a significant predictor of L. nidifica d15N val- Total N content of A. spicifera was related to fine-scale geo- ues for Control 1, Control 3, and Control 5. At two of these sites, graphic location where collection occurred but for L. nidifica, sam- values were lower during the rainy season while Control 3 values pled only at sites 1–8, the relationship was not evident. There was of d15N were higher in the rainy season. Drain 4 was the only a significant relationship between the similarity matrices of geo- regression model that had a significant interaction between dis- graphic location and total N content in dry (Rho = 36.4%, tance from shore and season (Table 3). p = 0.001) and rainy (Rho = 36.3%, p-value = 0.001) seasons for A. Significant relationships were found between geographic loca- spicifera. However there was no relationship between geographic tion and d15N values of A. spicifera for both seasons (RELATE: dry location and N content for either season for L. nidifica (RELATE: season, Rho = 80.9%, p-value = 0.001; rainy season, Rho = 65.4%, dry season, Rho = 4.0%, p-value = 0.69; wet season, Rho = 3.0%, p-value = .001). The d15N values of L. nidifica were significantly p-value = 0.29). 98 T. Erin Cox et al. / Marine Pollution Bulletin 71 (2013) 92–100
δ15N (‰) Laniakea 1.0 - 4.9 10.1 5.0 - 8.9 6.1
9.0 - 12.9 Kaena Point
>13.0 1.9 Species A. spicifera A. breviarticulatum L. mcdermidiae U. lactuca Kane‘ohe Bay Wai‘anae Coast 9.8 Maili-Kahe
Makapu‘u 11.9
9.6 6.1 8.4 3.8 2.7
N 16.8 9.7 km 7.3 7.7
One‘ula, Control 2, Drain 4 3.1 2.6 Diamond Head
Fig. 4. The island of O‘ahu showing mean d15N values (‰, n = 1–7) of reef algae (A. spicifera, A. breviarticulatum, Laurencia mcdermidae, Laurencia nidifica, and Ulva lactuca) collected along the coast in shallow waters in March 2008 (gray box) and October 2011.
3.2. Island wide d15N value comparison Storm-drain effluent alone does not dominate the input of nutrients supporting nearshore algal growth in the ‘Ewa Plain re- The mean d15N values of reef algae sampled at sites around gion studied. This conclusion is based on (1) the N content and O‘ahu ranged from 1.9‰ to 16.8‰; reef algae collected in 2011 d15N values in algae, which did not vary in a predictable manner from the ‘Ewa Beach area fell within the upper end of this range with drain presence (2) the nutrient concentrations in drain-water (7.3–16.8‰)(Fig. 4). Reef algae with the lowest mean d15N value that differed from concentrations measured in seawater (3) the of 1.9‰ were collected at Ka‘ena Point and the highest mean distinct differences in d15N values between algae collected from d15N value was determined from algal samples collected along drain sites and the d15N values of nitrate + nitrite measured in the ‘Ewa Beach area. The d15N values of A. spicifera collected at Con- storm-drain effluent and (4) the relatively consistent d15N values trol 2 and Drain 4 differed between 2011 and 2007. In 2007, A. spi- regardless of season and rainfall amounts (see online materials). cifera from Control 2 during the dry and rainy season had mean The latter is particularly revealing for a fast-growing species which values of 11.6‰ and 10.8‰ (respectively) but in 2011 had a lower likely uses all available substrate quickly and depletes internal N mean value of 7.7‰. The mean d15N values of A. spicifera at Drain 4 stores. Furthermore, algae located at Drain 9 which services the increased from 2007 to 2011 from 12.0‰ to 12.2‰ to 16.8‰. largest drainage area within the study region had the lowest values of d15N (4.8‰), a value that in other systems is more indicative of 4. Discussion non-anthropogenic sources (Garrison et al., 2007) or anthropo- genic sources that derived N from atmospheric sources (Kendall, The results of this study do not support the hypothesis that ni- 1998). trate in storm-drain effluent dominated the input of nutrients for Our findings are contrary to the conclusions drawn in previous nearshore algal growth over a several year evaluation with two investigations (Derse et al., 2007; Lapointe and Bedford, 2011). The species and end member water sources included in the analyses. discrepancy can be explained by differences in sampling design 15 Instead results suggested that reef algae at ‘Ewa Beach incorporate and interpretation of high d N values in algae. Assuming all avail- 15 N from a mixture of sources that cannot be readily identified, able nitrogenous nutrients are used by algae, the d N values of al- requiring more thorough investigation that goes beyond the typi- gae can provide information about nutrient sources particularly 15 cal approaches used for the measurements of d15N values in algae. when compared to d N values (nitrate + nitrite) of potential N The N content and d15N values in algae were often higher nearest sources (Peterson and Fry, 1987; Kendall, 1998; Robinson, 2001). to shore, tied to geographic location, and were similar or above val- However, elevated N microbial processing in an aquifer can lead 15 ues for algae collected from other nearshore environments. These to high d N values of nitrate (e.g., Lindsey et al., 2003) similar to findings support the conclusion that the nitrogenous nutrients values that often indicate anthropogenic enrichment (Dailer used by these plants are derived predominantly from terrestrial et al., 2012a). For this reason, a suite of additional indicators or sources that may more closely model usual coastal development. incorporation of spatial arrangement and temporal changes of al- 15 In areas with coastal development nutrients are predicted to be gae d N values can aid in interpretations (see Gartner et al., anthropogenic in origin, and undergo N cycling. These results are 2002; Lapointe et al., 2005; Cohen and Fong, 2005; Lin and Fong, in marked contrast to Dailer et al. (2010) but those Maui sites ran- 2008; Dailer et al., 2012a). 15 ged from near pristine to sites with focused delivery of up to Despite the differences between studies, elevated d N values 5 M gal d 1 of wastewater, nutrient loading not observed on O‘ahu. (>7‰) and N contents (>2%) in both species of macroalgae in the T. Erin Cox et al. / Marine Pollution Bulletin 71 (2013) 92–100 99
‘Ewa Beach area fell above the range commonly cited for animal Peterson et al., 2007, 2009; Johnson et al., 2008; Street et al., waste (Kendall, 1998; Macko and Ostrom, 1994) but within the 2008; Peterson et al. 2009). Further investigations are needed to range from studies in other locations that investigated anthropo- identify if d15N values that occur along ‘Ewa Beach are associated genic enrichment in tropical and subtropical waters. Published val- with seep locations and studies including quantification of excess 15 ues of d N of various warm water macroalgae range from +1.0‰ N2 gas in groundwater (Lindsey et al., 2003) are needed to distin- to +5.5‰ in the Florida Keys (interpreted as being indicative of guish whether high d15N values in reef algae result from the incor- agricultural runoff, Lapointe et al., 2004), from +5.7‰ to +12.0‰ poration of a land based anthropogenic or natural nitrogenous in waste water influenced locations of South Florida (Lapointe sources that have been modified significantly by microbial reaction et al., 2005), from +7.7‰ to +11.4‰ in southeastern Gulf of Califor- in aquifers. nia (plausible sources included sewage, agriculture, and shrimp In conclusion, our study revealed a complex pattern of nutrient farms, Piñón-Gimate et al., 2009), a mean of 0.5‰ in Hanalei sources that is most congruent with localized delivery of terrestrial Bay, Kaua‘i island, Hawai‘i (concluded to be indicative of synthetic N source (s) for much of the reef algae along the coast that pro- fertilizers from nearby Resorts, Derse et al., 2007), from +2‰ to duces an onshore–offshore trend. No drain effect was evident +8‰ in Ishigaki, Japan (Umezawa et al., 2002), and +1.0‰ to but, we caution that results are not necessarily applicable to other +4.8‰ in Ofu, Samoa (concluded to be indicative of non-anthropo- areas. For instance, the effects of storm-drain effluent on bloom- genic sources, Garrison et al., 2007). Nonetheless, we cannot distin- forming species could have a strong impact in more closed, oligo- guish unequivocally whether the elevated d15N values (>7‰) trophic waters. The results of this study should be considered in observed for Ewa Beach macroalgae result from microbial cycling the design of future assessment schemes using similar techniques, of natural or anthropogenic N sources in soils or groundwater in decisions about how to manage water flow, and can be used as aquifers. baseline to measure future alterations to nutrient inputs in ‘Ewa The general enrichment in 15N contents of algae collected from Beach. The study also highlights the value of considering physio- ‘Ewa Beach shorelines relative to many other sites around O‘ahu logical differences among the algae sampled, the spatial-temporal suggests that the N source in this area could be anthropogenic in patterns of their d15N values and the biogeochemical environment origin. Recent studies by Dailer et al. (2010) on Maui have shown where they occur. remarkably high macroalgal d15N values (>40‰) for field collected materials, well above those recorded in this study for algal material Acknowledgments collected from O‘ahu. In that study, there was a single clear anthro- pogenic N source that could be linked to a suite of sewage injection We would like to acknowledge K. Boyle, P. Buxton, M. Lurie, E. wells using multiple tracers including elevated N and phosphorous Donham, S. V. Gent, M. nGiriou, and M. Kawachi for their assistance levels and high concentrations of pharmaceuticals (Hunt and Rosa, in collecting or processing samples for this study. Also we thank 1 2009). These wells deliver up to 5 M gal d of partially denitrified the staff of University of Hawai‘is Stable Isotope Biogeochemistry primary treated sewage effluent into groundwater (Dailer et al., Laboratory for guidance in stable isotope analyses. The study was 2010, 2012a). In this study, no single signal was detected. Further, supported by Hasekeko Development, Inc. We also acknowledge there is no obvious relationship between elevated marine macroal- the State of Hawai‘i, Division Aquatic Resources for permitting us 15 gal d N values and proximity to known on-site sewage disposal to collect inside the ‘Ewa Beach Limu Management Area. This is SO- systems (Whittier and El-Kadi, 2009). EST contribution #8904. Our d15N values of reef algae were often related to the fine- scale, geographic location where collection occurred. Average d15N values of reef algae varied among sites from 4.8‰ to 14.7‰, Appendix A. Supplementary material a difference of nearly 10‰ that was greater than the largest varia- tion between seasons (5.8‰ difference), the two species at a site Supplementary data associated with this article can be found, in within a season (1.3‰ difference), and individuals (3.8‰ differ- the online version, at http://dx.doi.org/10.1016/j.marpolbul.2013. ence) at a site within a season. Furthermore, this variation occurred 03.030. over a small geographic area and was also observed along shore- lines in other areas of O‘ahu (e.g., Wai‘anae Coast). The geographic References affinity of d15N values and %N of reef algae reveal that the N sources in these areas act on a local scale. However, obvious point Abbott, I.A., 1996. Limu, An Ethnobotanical Study of Some Hawaiian Seaweeds, sources of N were difficult to identify. fourth ed. National Tropical Botanical Garden, Lawai, Kauai, Hawai‘i. Banner, A.H., 1974. Kane‘ohe Bay, Hawaii: Urban pollution and a coral reef Groundwater seeps are thought to occur along the coast in ‘Ewa ecosystem. In: Proc. 2nd Int. Coral Reef Symp. vol. 2, pp. 685-702. Beach (Spengler et al., 1998; Laws et al., 1999). Groundwater can Bernardo, R., 2008. Limu Delays Project to Ease Ewa Flooding. Star Bulletin, deliver substantial amounts of nutrients into coastal systems Honolulu.
dimensionally model the plume across a coral reef on Maui, Hawai‘i, USA. Mar. Owens, N.P.J., 1987. Natural variations in d15N in the marine environment. Adv. Mar. Pollut. Bull. 64, 207–213. Biol. 24, 389–451. Dailer, M.L., Smith, J.E., Smith, C.M., 2012b. Responses of bloom forming and non- Peterson, B.J., Fry, B., 1987. Stable isotopes in ecosystem studies. Ann. Rev. Ecol. bloom forming macroalgae to nutrient enrichment in Hawai‘i, USA. Harmful Syst. 18, 293–320. Algae 17, 111–125. Peterson, R.N., Burnett, W.C., Glenn, C.R., Johnson, A.J., 2007. A box model to Derse, E., Knee, K.L., Wankel, S.D., Kendall, C., Berg, C., Paytan, A., 2007. Identifying quantify groundwater discharge along the kona coast of hawaii using natural sources of nitrogen to Hanalei Bay, Kauai utilizing the nitrogen isotope tracers. In: Sanford, W., Langevin, C., Polemio, M., Povinec, P. (Eds.), A New Focus signature of macroalgae. Environ. Sci. Technol. 41, 5217–5223. on Groundwater–Seawater Interactions. IAHS Press, pp. 142–149. Doty, M.S., 1961. Acanthophora, a possible invader of the marine flora of Hawaii. Pac. Peterson, R.N., Burnett, W.C., Glenn, C.R., Johnson, A.G., 2009. Quantification of Sci. 15, 547–552. point-source groundwater discharges to the ocean from the shoreline of the Big Dulaiova, H., Gonneea, M.E., Henderson, P.B., Charette, M.A., 2008. Geochemical and Island, Hawaii. Limnol. Oceanogr. 54, 890–904. physical sources of radon variation in a subterranean estuary – implications for Piñón-Gimate, A., Soto-Jiménez, M., Ochoa-Izaguirre, M.J., García-Pagés, E., Páez- radon groundwater end-member activities in submarine groundwater Osuna, F., 2009. Macroalgae blooms and d15N in subtropical coastal lagoons discharge studies. Mar. Chem. 110, 120–127. from the Southeastern Gulf of California: Discrimination among agricultural, Fong, P., Kamer, K., Boyer, K.E., Boyle, K.A., 2001. Nutrient content of macroalgae shrimp farm and sewage effiuents. Mar. Pollut. Bull. 58, 1144–1151. with differing morphologies may indicate sources of nutrients for tropical Robinson, D., 2001. D15N as an integrator of the nitrogen cycle. Trends Ecol. Evol. 16, marine systems. Mar. Ecol. Prog. Ser. 220, 137–152. 153–162. Fong, P., Boyer, K.E., Kamer, K., Boyle, K.A., 2003. Influence of initial tissue nutrient Russell, D.J., 1992. The ecological invasion of Hawaiian reefs by two marine red status of tropical marine algae on response to nitrogen and phosphorus algae, Acanthophora spicifera (Vahl) Boerg. and Hypnea musciformis (Wulfen) J. additions. Mar. Ecol. Prog. Ser. 262, 111–123. Ag., and their association with two native species, Laurencia nidifica J. Ag and Garrison, G.H., Glenn, C.R., McMurtry, G.M., 2003. Measurement of submarine Hypnea cervicornis. J. Agric. ICES Mar. Sci. Symp. 194, 110–125. groundwater discharge in Kahana Bay, O‘ahu, Hawai‘i. Limnol. Oceanogr. 48, Sammarco, P.W., Risk, M.J., Schwarcz, H.P., Heikoop, J.M., 1999. Cross-continental 920–928. shelf trends in coral d15N on the Great Barrier Reef: further consideration of the Garrison, V., Kroeger, K., Fenner, D., Craig, P., 2007. Identifying nutrient sources to reef nutrient paradox. Mar. Ecol. Prog. Ser. 180, 131–138. three lagoons at Ofu and Olosega, American Samoa using d15N of benthic Schaefers, A., 2006. Ocean Pointe Drainage Plan Under Attack. Star Bulletin, macroalgae. Mar. Pollut. Bull. 54, 1813–1838. Honolulu.
Coastal water quality in Hawaii: the importance of buer zones and dilution E.A. Laws a,*, D. Ziemann b, D. Schulman a aUniversity of Hawaii, Department of Oceanography, 1000 Pope Road, Honolulu, HI 96822-2285, USA bOceanic Institute, 41-202 Kalanianaole Highway, Waimanalo, HI 96795-1898, USA
Received 20 June 1998; received in revised form 10 December 1998; accepted 27 December 1998
Abstract
A study of the relationship between point and nonpoint source freshwater discharges and marine water quality were studied during a period of 1 year in Mamala Bay, a coastal inden- tation on the south shore of the island of Oahu, Hawaiian Islands. Despite the fact that 100± 300Â106 m3 year 1 of land runo/groundwater seepage and 150Â106 m3 year 1 of treated sewage euent enter Mamala Bay and its tributaries, coastal water quality as judged by stand- ard chemical and physical parameters is high at virtually all locations in the bay. The expla- nation for the high water quality re¯ects several important factors. First, much of the nonpoint source discharge enters either estuaries or harbors, which function as buer zones by trapping some of the sediment and nutrients that would otherwise enter the coastal ocean. Second, the principal point source discharges are located in water suciently deep that their wastewater plumes are trapped below the surface most of the time. When the plumes surface they are suciently diluted that their impact on parameters, such as nutrient concentrations, is unde- tectable. Third, the coastal current system is greatly diluted by exchange with the oshore ocean. Based on a simple box model, the degree of mixing with the oshore ocean is roughly 40 times the rate of input of fresh water from point and nonpoint sources. The oshore waste- water outfalls have no discernible eect on water quality at any recreational beach along the shoreline. The principal impact on water quality at the recreational beaches comes from non- point source discharges, and with the exception of one beach located directly adjacent to a stream mouth, that impact is on the composition rather than the concentration of the plank- ton. There is a systematic shift from a chlorophyte- to a diatom-dominated phytoplankton community due to the high silicate concentration in groundwater and land runo, and there is a systematic increase in the d15N of suspended particles due to the high d15N of the biologically available nitrogen in groundwater seepage. # 1999 Elsevier Science Ltd. All rights reserved. Keywords: Buer zone; Coastal zone; Dilution; Estuaries; Mixing processes; Ocean disposal; Outfalls; Recreational waters; Sewage; Water quality
* Corresponding author. Tel.: +1-808-956-7633; fax: +1-808-956-9225.
0141-1136/99/$ - see front matter # 1999 Elsevier Science Ltd. All rights reserved. PII: S0141-1136(99)00029-X 2 E.A. Laws et al. / Marine Environmental Research 48 (1999) 1±21
1. Introduction
Mamala Bay is a coastal indentation that extends a distance of about 30 km along the southern shoreline of the island of Oahu in the Hawaiian Islands, from Diamond Head in the east to Barbers Point in the west (Fig. 1). The city of Honolulu borders much of the shoreline of Mamala Bay, and along roughly the eastern third of the bay's shoreline are located the hotels and beaches of Waikiki, one of the major tourist destinations of the world. The attraction of the Hawaiian Islands to tourists depends very much on the quality of the environment, including in particular near- shore water quality. Because of the importance of the tourism industry to the state of Hawaii, the state has been much concerned with maintaining high water quality standards in nearshore recreational waters. Maintenance of high water quality is a particularly sensitive issue in an area such as Mamala Bay, which currently receives wastewater in the form of primary treated sewage from a population of roughly 750,000 persons. During 1993 and 1994 a comprehensive study of water quality in Mamala Bay was carried out under the direction of the Mamala Bay Commission (1996). One of the concerns of that study was the impact on water quality of nutrients introduced into the bay from both point and nonpoint sources. Stream water runo and ground- water seepage to the bay have been estimated to be somewhere between 100Â106 m3 year 1 (Freeman, 1993) and 300Â106 m3 year 1 (Stevenson, O'Connor, & Aldrich,
Fig. 1. Mamala Bay showing locations of principal wastewater treatment plant outfalls. Beaches moni- tored were Diamond Head (DH), Queens Surf (QS), Fort DeRussy (FD), Ala Moana (AM), Sand Island (SI), Keehi Lagoon (KL), Fort Kamehameha (FK), Ewa Beach (EB), Oneula Beach (OB), and Barbers Point (BP). Dashed line indicates approximate path of oshore transects and locations of Stations 1±17, where vertical pro®ling work was done. GPS=Global positioning system. E.A. Laws et al. / Marine Environmental Research 48 (1999) 1±21 3
1996). The associated inputs to Mamala Bay of nutrients such as nitrogen (N) and phosphorus (P) are known with less accuracy, since much of the fresh water enters estuaries and harbors, which can be rather eective nutrient traps (Laws, Hiraoka, Mura, Punu, & Yamamura, 1994). Roughly 65% of the stream runo and ground- water seepage enters Pearl Harbor, about 20% enters Keehi Lagoon/Honolulu Harbor, and 10% enters the Ala Wai Canal (Fig. 1). The nutrient delivery from stream runo and groundwater seepage to these estuaries/harbors has been esti- mated to be 30±70 tonnes year 1 of total phosphorus (TP) and 300±700 tonnes year 1 of total nitrogen (TN) (Freeman, 1993; Stevenson et al., 1996). How much of this N and P reaches Mamala Bay is problematic. Estuaries can be ecient sediment and nutrient traps, but their trapping eciency is a function of many factors, including loading rate, estuarine morphology, and water residence time. They may function as sources of inorganic nutrients due to the remineralization of organic matter, particularly on the bottom (Nixon, 1981; Correll, Jordan, & Weller, 1991), but in general they are sinks rather than sources of TN and TP. The Ala Wai Canal is a type B estuary as de®ned by Pritchard (1967). In such estuaries, ``Very little of the suspended sediment introduced from upland sources escapes through the estuary to the sea'' (Biggs & Cronin, 1981, p. 10). Laws et al. (1994) have estimated that sedimentation removes about 40% of the allochthonous inputs of organic carbon to the Ala Wai Canal. It seems likely that a comparable percentage of the organic N and P inputs are also removed via sedimentation. In the case of N, an additional mechanism of removal is denitri®cation (Nixon, 1981). Pearl Harbor is a type C estuary, with an area of 20.1 km2 and a mean tidal range of 0.37 m (Buske, 1974); its trapping eciency is unknown. Point source inputs are known with more accuracy and are summarized in Table 1. The principal inputs come from sewage treatment plants operated by the City and County of Honolulu, namely the Sand Island wastewater treatment plant (WWTP) and the Honouliuli WWTP. Both plants discharge primary treated euent directly into Mamala Bay (Fig. 1). An additional source of sewage is the Fort Kamehameha WWTP, which discharges secondary treated sewage into the mouth of Pearl Harbor (Fig. 1). The discharge of TN and TP from all point sources is estimated to be 2800
Table 1 Principal point source inputs of total nitrogen (TN) and total phosphorus (TP) to Mamala Bay (Stevenson et al., 1996)
Source Nature of Freshwater discharge TN discharge TP discharge discharge (106 m3 year 1) (tonnes year 1) (tonnes year 1)
Sand Island WWTP Primary treated sewage 100 1890 274 Honouliuli WWTP Primary treated sewage 33 800 123 Fort Kamehameha Secondary treated sewage 10 110 16 WWTP Other 7 Total 150 2800 413
WWTP, wastewater treatment plant. 4 E.A. Laws et al. / Marine Environmental Research 48 (1999) 1±21 and 413 tonnes year 1, respectively (Stevenson et al., 1996). While it seems clear from this analysis that point sources dominate the nutrient loading to Mamala Bay, it is important to bear in mind that the Sand Island and Honouliuli WWTPs dis- charge their euent 2 km from shore at depths of 73 and 61 m, respectively. The discharges from harbors/estuaries occur directly along the shoreline and impact the uppermost portion of the water column. Furthermore, stream runo during storms can easily be 10±100 times dry weather ¯ow (Laws, 1993). Thus, from the standpoint of aecting water quality in recreational areas, the nonpoint source discharges are potentially of comparable or even greater importance than the sewer outfalls. The study reported here was aimed at investigating patterns in nutrient-related water quality parameters on recreational beaches and along the approximate isobath of the two major sewer outfalls in Mamala Bay, with the expectation that some insights could be gained about the mechanisms that control water quality in the bay and about the relative importance of point and nonpoint nutrient sources at locations where impacts from fresh water were discernible.
2. Materials and methods
2.1. Beach sampling
Water samples were collected at the 10 beaches shown in Fig. 1 on a monthly basis from August 1993 through July 1994. The samples were collected from just below the surface in acid-cleaned 20-l plastic carboys at a distance from the shoreline where the water column was 1-m deep. Salinity was measured with an Extech Oyster con- ductivity meter calibrated with Copenhagen seawater. Upon return to the laboratory, aliquots of the water were ®ltered through glass ®ber Whatman GF/F ®lters for fur- ther analysis. The ®ltrates were used to determine inorganic nutrient concentrations. Silicate, molybdate-reactive phosphate (MRP), nitrate+nitrite, and ammonium +ammonia concentrations were determined by colorimetric methods on a Technicon AutoAnalyzer. Concentrations of particulate carbon (PC) and particulate nitrogen (PN) were determined from aliquots ®ltered onto precombusted GF/F ®lters and analyzed using a Perkin-Elmer model 2400 elemental analyzer. In addition, the iso- topic composition of the PN (d15N) was determined using either a Finnigan MAT 252 or Delta-S mass spectrometer. A modi®ed Dumas sealed tube method was used to convert the PN to N2 as described by Minawaga, Winter, and Kaplan (1985). Dried ®lters containing the PN were placed in 9-mm outer diameter Vycor brand quartz tubes along with 3 g of pre-combusted cupric oxide and 1 g of copper wire. The tubes were evacuated, sealed, and combusted at 650C overnight. The tempera- ture was then lowered to 500C for 1 h, after which time they were allowed to cool slowly to room temperature. Slow cooling ensured that the oxides of N created by combustion were all converted to N2.N2 was isolated by cryogenic distillation, using silica gel to adsorb the N2 (Mariotti, 1983). Samples for chlorophyll a (chl a) and diagnostic carotenoid and chlorophyll pigment analyses were collected on GF/F
®lters and placed in liquid N2 prior to extraction. The pigments were extracted by E.A. Laws et al. / Marine Environmental Research 48 (1999) 1±21 5 placing the ®lters in 5 ml of acetone for 12 h. The ®lters were then ground using a Wheaton tissue grinder and the ®lter residue removed by centrifugation. The extract- ed pigment concentrations were measured on a Varian 5000 high-performance liquid chromatograph (HPLC) as described by Laws, DiTullio, Betzer, and Hawes (1990).
2.2. Oshore ®eld work
Six cruises were carried out between December 1993, and December 1994, at approximately 2-month intervals. The cruises consisted of horizontal transects along the 60-m isobath between Diamond Head and Barbers Point and vertical pro®les at selected stations located at depths of 20, 60, and 120 m. Typical horizontal transect tracks and the locations of the vertical pro®le stations are shown in Fig. 1. Hydro- graphic pro®les of temperature, salinity, and beam transmissometry to a maximum depth of 100 m were conducted with an Applied Microsystems CTD12 system ®tted with a 25-cm beam transmissometer (CTTD). The CTTD was ®tted with internal RAM and operated in a stand-alone mode. The CTTD was programmed to record data at 2-s intervals for all vertical pro®les, and data were downloaded to a portable computer at the completion of each cast. Water samples for analysis of chemical and biological parameters were collected by Niskin bottles during the December 1993 cruise and with a submersible pump system on subsequent cruises. The system con- sisted of a 110-volt submersible pool pump, portable generator, 60 m of 2-cm inside- diameter opaque plastic hose, and assorted ®ttings. Samples were transferred to polyethylene bottles and chilled in ice. In the laboratory the samples were split, with one fraction ®ltered onto 0.2-mm Nuclepore membrane ®lters for analysis of extract- ed chl a. The other fraction was analyzed (un®ltered) for dissolved nutrients and turbidity. Turbidity was measured using a Turner Designs nephelometer. Inorganic nutrients were analyzed using colorimetric methods with a Technicon Auto- Analyzer, with the exception of nitrate+nitrite, which was analyzed using the chemiluminescent procedure of Garside (1982) on an Antek model 703C nitrogen oxides analyzer.
2.3. Sewage bioassay experiments
Bioassays were carried out using euent from either the Sand Island or Honouliuli WWTPs. The euent was, in each case, a 24-h composite sample obtained from the City and County of Honolulu's Sand Island wastewater laboratory. The sewage was ®ltered through precombusted GF/F ®lters. The particulate matter on the ®lters was analyzed for PC and PN and the d15N of the PN as previously described. The ®ltrate was analyzed for inorganic nutrient concentrations as previously described using colorimetric methods on a Technicon AutoAnalyzer. Concentrations of total dissolved nitrogen (TDN) and total dissolved phosphorus (TDP) were determined by ultraviolet oxidation of the ®ltrates followed by colorimetric analysis of the inorganic N and MRP. Salinity was measured as previously described with an Extech Oyster conductivity meter. 1-week incubations were conducted to determine the d15Nof the biologically available dissolved N in the sewage euent and to determine the 6 E.A. Laws et al. / Marine Environmental Research 48 (1999) 1±21 impact of sewage enrichment on the composition of the phytoplankton community. For these experiments, 100 ml of sewage ®ltrate were added to each of two clear poly- carbonate bottles containing 900 ml of surface seawater collected 1 km o Diamond Head. The bottles were incubated at a temperature of 25C in front of a bank of day- light ¯uorescent lamps at a continuous irradiance of 150 mmol quanta m 2 s 1 (PAR). One bottle was incubated without further nutrient additions and the cells harvested by ®ltration after 1 week for pigment analysis. The second bottle was enriched with the mixture of trace metals, vitamins, and phosphate recommended for IMR (Institute of Marine Resources) medium (Eppley, Holmes, & Strickland, 1967). The phosphate concentration produced by this enrichment was about 50 mM, and the resultant ratio of dissolved inorganic N (DIN) to phosphate was about 2 on a molar basis. The enrich- ment thus assured that the phytoplankton would strip the water of available N before signi®cantly depleting the phosphate concentration, since the ratio of N to P in phyto- plankton is about 16 and virtually never falls outside the range 3±30 on a molar basis (Ryther & Dunstan, 1971). After 1 week virtually all DIN had been stripped from the water. The phytoplankton cells were harvested by ®ltration onto GF/F ®lters and the particulate material analyzed for d15N as previously described.
2.4. Bioassays with stream water and groundwater
In addition to the 1-week incubations conducted with sewage euent, similar bioassays were conducted with stream runo and groundwater to determine the d15N of the biologically available dissolved N. Seven streams were sampled, ®ve of which ¯ow into Pearl Harbor (Aiea, Kalauao, Waiawa, Waikele, Waimalu streams), one of which ¯ows into Keehi Lagoon/Honolulu Harbor (Manaiki stream), and one of which ¯ows into the Ala Wai Canal (Manoa-Palolo stream). The groundwater samples were taken from two wells in the Ewa Plain. For the stream water assays, several liters of water were collected and transferred to clear polycarbonate bottles, which were then enriched with IMR concentrations of trace metals, vitamins, and phosphate. The natural phytoplankton community in the streams was allowed to grow until virtually all DIN had been stripped from the water. In the case of the groundwater, 200 ml were added to 800 ml of Diamond Head surface seawater, and the mixture enriched with IMR concentrations of trace metals, vitamins, and phos- phate. The Diamond Head phytoplankton community was then allowed to grow until virtually all DIN had been stripped from the water. The phytoplankton cells were then harvested by ®ltration onto GF/F ®lters and the particulate material analyzed for d15N as previously described.
3. Results
3.1. Sewage, stream, and groundwater analyses
A summary of the analyses of the Sand Island and Honouliuli sewage euent is given in Table 2. The silicate concentrations of 700±1100 mM re¯ect interactions E.A. Laws et al. / Marine Environmental Research 48 (1999) 1±21 7
Table 2 Summary of results of analyses of euent from Sand Island and Honouliuli wastewater treatment plant euents
Sand Island Honouliuli
Phosphate (MRP) 548 723 Nitrate+nitrite 5.78.1 7.72.8 Ammonium+ammonia 857115 130258 Silicate 69859 108121 Total dissolved N 109461 1684153 Total dissolved P 589 74.33.4 Dissolved organic N 231179 37695 Dissolved organic P 3.61.2 2.72.6 Particulate C 2670219 3162760 Particulate N 474334 37985 PC/PN 6.32.5 7.10.2 d15N 1.93.1 1.83.5 d15N of biologically available N 6.95.8 10.12.2 Salinity 4.40.2 0.70.1
Nutrient concentrations and concentrations of particulate carbon (PC) and particulate nitrogen (PN) are micromolar. d15N is expressed per mil. Salinity is given in practical salinity units. Error bars are one standard deviation based on three analyses. MRP, molybdate-reactive phosphate; N, nitrogen; P, phos- porous; C, carbon; d15N, isotopic composition of PN. between groundwater and basaltic rock and are typical of groundwater in Hawaii. Virtually all the DIN in the sewage consisted of ammonium+ammonia. The d15Nof the biologically available dissolved N in the sewage averaged 7±10%, about 5±8% higher than the d15N of the PN in the sewage. The d15N of the biologically available dissolved N in stream water and groundwater averaged 2.42.7 and 10.40.7%, respectively.
3.2. Beach analyses
The results of the beach sampling are shown in Figs. 2±4. Salinities were relatively constant and close to the salinity of oshore surface water (34.5%) at beaches from Diamond Head to Sand Island. Salinity was clearly depressed at Keehi Lagoon and recovered only gradually toward more typical seawater values at Oneula Beach and Barbers Point. The much higher variability of the salinity readings at Keehi Lagoon Beach Park re¯ects the highly variable nature of stream runo and the fact that the beach park lies directly adjacent to the mouths of the Moanalua and Kalihi streams. Silicate concentrations were approximately a mirror image of the salinity values due to the high concentration of silicate in Hawaiian streams. The con- centrations from Diamond Head to Sand Island were uniformly low and typical of oshore ocean values (1.5 mM). The values at Keehi Lagoon again showed by far the greatest variability due to the intermittent nature of storm runo events. Mean PN and chl a values were highly correlated (r2=0.97). The highest values and greatest variability occurred at Keehi Lagoon. Chl a values were uniformly low at 8 E.A. Laws et al. / Marine Environmental Research 48 (1999) 1±21
Fig. 2. (A) Mean salinities, (B) concentrations of silicate, (C) chlorophyll a, and (D) particulate nitrogen (PN) at the 10 recreational beaches. Error bars are standard errors. BP, Barbers Point; OB, Onuela Beach; EB, Ewa Beach; FK, Fort Kamehameha; KL, Keehi Lagoon; SI, Sand Island; AM, Ala Moana; FD, Fort DeRussy; QS, Quenns Surf; DH, Diamond Head.
Diamond Head, Queens Surf, Sand Island, Fort Kamehameha, and Barbers Point. Slightly higher and more variable values occurred at Ala Moana, Fort DeRussy, Oneula Beach, and Ewa Beach. PN values followed a similar pattern with the exception of Queens Surf, where the results were slightly higher and more variable compared to the pattern of chl a concentrations. Based on a one-way analysis of variance (ANOVA) there was no signi®cant dierence ( p>0.05) between the MRP concentrations at the 10 beaches. The results were uniformly variable, and all mean values fell within the range 0.150.07 mM. There was a highly signi®cant dierence in nitrate+nitrite concentrations (ANOVA on log-transformed concentrations, p=0.02), but the dierence was due entirely to the high and variable results at Keehi Lagoon Beach Park. Mean concentrations at the remaining beaches were more- or-less uniform, and all were less than 0.4 mM. Concentrations of ammonium +ammonia were uniformly higher than nitrate+nitrite, with all mean values falling in the range 0.80.4 mM. There was again a signi®cant dierence in concentrations between beaches ( p=0.03), the dierence in this case being due to the combina- tion of the high values at Ewa Beach and the low values at Fort Kamehameha. The d15N of the PN diered signi®cantly ( p=10 8) and dramatically between beaches. Values systematically increased from east to west, from a low of about 3.5% at Diamond Head to a high of almost 7% at Barbers Point. Analysis of the diagnostic carotenoid and chlorophyll pigments in the suspended particulate material revealed two general patterns. At Keehi Lagoon and at beaches E.A. Laws et al. / Marine Environmental Research 48 (1999) 1±21 9
Fig. 3. Mean concentrations of (A) phosphate, (B) nitrate+nitrite, (C) ammonium+ammonia, and (D) the d15N of particulate nitrogen (PN) at the 10 recreational beaches. Error bars are standard errors. BP, Barbers Point; OB, Onuela Beach; EB, Ewa Beach; FK, Fort Kamehameha; KL, Keehi Lagoon; SI, Sand Island; AM, Ala Moana; FD, Fort DeRussy; QS, Quenns Surf; DH, Diamond Head. east of Sand Island the dominant diagnostic pigment was chl b, and the ratio of 190- hexanoyloxyfucoxanthin (190-hex or hex) to chl b was consistently about 0.5. At Sand Island and at beaches west of Keehi Lagoon the dominant diagnostic pigment was either fucoxanthin or (at Ewa Beach) 190-hex, and the 190-hex/chl b ratio was consistently greater than 0.5. In Diamond Head water enriched with sewage euent, the dominant diagnostic pigment was fucoxanthin, but the 190-hex/chl b ratio was only about 0.1, much lower than at any of the beaches.
3.3. Oshore ®eld work
The results of the oshore ®eld work are shown in Figs. 5±9. Plots of ammonium +ammonia versus silicate and nitrate+nitrite versus silicate reveal the presence of three water types. In most cases the water contained low concentrations of all three nutrients. However, some water samples contained relatively high concentrations of both ammonium+ammonia and silicate, and other water samples contained rela- tively high concentrations of both nitrate+nitrite and silicate. Examination of the depths from which these high-nutrient samples were collected reveals that the former group came from depths of 10 m or greater, and the latter group came from depths of 5 m or less. The samples with the highest concentrations of ammonium+ammo- nia and silicate came from depths of 60 m, and the samples with the highest con- centrations of nitrate+nitrite and silicate came from depths of either 1 or 2 m. 10 E.A. Laws et al. / Marine Environmental Research 48 (1999) 1±21
Fig. 4. Pro®les of normalized diagnostic pigment concentrations for the ®ve beaches at which (A) chlo- rophyll b was the dominant accessory pigment, (B) either fucoxanthin (fuc) or 190-hexanoyloxyfucox- anthin (hex) was the dominant accessory pigment, and (C) sewage-enriched Diamond Head water. Error bars are standard errors.
Further analysis of the relationship between nutrient concentrations and sample location reveals that the high ammonium+ammonia results came almost exclusively from stations in the vicinity of the Sand Island and Honouliuli WWTP outfalls, whereas the high nitrate concentrations came from stations o the mouths of the Ala Wai Canal and Pearl Harbor. Some additional insight into the nature of oshore water quality can be gained by examining the frequency distribution of nutrient concentrations. For example, when results from samples collected from depths of 1 or 2 m are excluded, the nitrate+ nitrite concentrations give an excellent ®t to a log-normal distribution Fig. 7, with a median of 8 nM. This ®gure is well below the limit of detection by traditional col- orimetric methods (0.03 mM), but is comfortably above the detection limit (0.3 nM) of Garside's (1982) chemiluminescent assay. Determining the background marine distribution of phosphate, silicate, and ammonium+ammonia was not as straightforward since both land runo/groundwater seepage and the euent from the oshore WWTP outfalls contain high concentrations of phosphate and silicate, and the euent from the oshore WWTP outfalls is greatly enriched in ammonium +ammonia. In order to estimate the concentrations associated with strictly marine waters, we excluded samples taken from depths of 1±2 m and from stations in the immediate vicinity of the WWTP outfalls (Stations 1±3 and 13±15). The median phosphate, silicate, and ammonium+ammonia concentrations at the other stations E.A. Laws et al. / Marine Environmental Research 48 (1999) 1±21 11
Fig. 5. Scatter diagrams of (A) silicate versus nitrate+nitrite concentrations, and (B) silicate versus ammonium+ammonia concentrations. Numbers next to symbols indicate depths from which samples with high concentrations were taken. and depths were 0.06, 1.5, and 0.06 mM, respectively. The median phosphate and ammonium+ammonia concentrations approach the limit of detection (0.03 mM) by colorimetric methods. The turbidity of most water samples was less than 0.2 nepholometric turbidity units (NTUs), but exceptions to this pattern were clearly evident o the mouths of 12 E.A. Laws et al. / Marine Environmental Research 48 (1999) 1±21
Fig. 6. Concentrations of (A) nitrate+nitrite, (B) turbidity, and (C) ammonium+ammonia versus o- shore station number. Numbers next to symbols indicate depths from which samples with high con- centrations were taken. NTUs, nepholometric turbidity units. the Ala Wai Canal and Pearl Harbor. Turbidity was closely correlated with the dis- tribution of chl a. By far the highest and most variable chl a and turbidity were found o the mouth of the Ala Wai Canal (Station 4) and Pearl Harbor (Station 9). Comparisons were made between the chl a and turbidity readings at depths of 20 and 60 m in the immediate vicinity of the Sand Island and Honouliuli WWTP out- falls and at other stations in the bay to determine whether the impact of the outfalls E.A. Laws et al. / Marine Environmental Research 48 (1999) 1±21 13
Fig. 7. Frequency distribution of oshore nitrate+nitrite concentrations, excluding results from samples taken at depths of 1±2 m. Asterisks indicate expected frequencies based on a log-normal distribution of values with the same mean and standard deviation. on these parameters was discernible in a statistical sense. Based on a t-test, there was no signi®cant dierence between chl a concentrations at either 20 or 60 m near the Sand Island outfall and at other stations in the bay ( p>0.35). However, chl a at both 20 and 60 m was judged to be higher near the Honouliuli outfall than at other stations ( p<0.05). There was no dierence in turbidity at 20 m ( p>0.55), but at 60 m turbidity was signi®cantly dierent and higher at stations near both the Sand Island and Honouliuli outfalls compared to other stations in the bay ( p<0.04). We estimated the chl a concentration of strictly marine waters in the bay by excluding shallow water (<5 m) samples collected immediately o the mouth of the Ala Wai Canal (Station 4) and Pearl Harbor (Station 9) and samples from the immediate vicinity of the two WWTP outfalls. The median chl a concentration of the remaining samples was 0.12 mg/l. Vertical pro®les of temperature, salinity, and turbidity sometimes revealed very obvious gradients in the upper 60 m of the water column and at other times a rela- tively well-mixed system. The data from August 1994 are illustrative of the former condition. The temperature of the water dropped by about 2C between 40 and 55 m, and this gradient was sucient to trap the sewage plumes. The plumes were clearly evident as increases in turbidity at roughly 50 m depth (Fig. 9B,C). No such behav- ior was apparent in December 1993. 14 E.A. Laws et al. / Marine Environmental Research 48 (1999) 1±21
Fig. 8. Mean turbidity and chlorophyll (chl) a concentrations measured at Stations 4 (Pearl Harbor) and 9 (Ala Wai Canal) and at depths of 20 and 60 m over the Sand Island and Honouliuli outfalls and at depths of 20 and 60 m at other stations in the bay. Error bars are standard errors. NTUs, nepholometric turbidity units.
4. Discussion
4.1. Recreational beaches
Of all the recreational beaches included in this study, Keehi Lagoon Beach Park was clearly the most impacted by fresh water. This conclusion is evidenced by its low mean salinity and high concentrations of silicate, chl a, PN, and nitrate+nitrite. The highly variable nature of all these parameters at Keehi Lagoon re¯ects the fact that the impact comes from stream runo, which is temporally highly variable (Laws, 1993). The temporal variability would be much less if the impact were primarily from groundwater seepage or wastewater discharges, which are relatively constant. Both the magnitude of the impact and the temporal variability would be less extreme if the impact were due to the out¯ow from an estuary as opposed to the direct dis- charge from several streams, since estuaries can be rather ecient traps of sediment and N (Biggs & Cronin, 1981; Nixon, 1981; Laws et al., 1994) and can eectively convolute the time series of fresh water and nutrient inputs on a time scale com- parable to the residence time of water in the estuary. The latter consideration is particularly relevant to Pearl Harbor, which has a narrow mouth and within which the ratio of volume to fresh water discharge is 1.0 year (Freeman, 1993). The impacts of freshwater discharges on water quality at beaches other than Keehi Lagoon appear to be relatively subtle. Of the remaining beaches, the one E.A. Laws et al. / Marine Environmental Research 48 (1999) 1±21 15
Fig. 9. Vertical pro®les of temperature, salinity, and turbidity above the Sand Island outfall measured in August 1994 (A±C) and December 1993 (D±F). experiencing the greatest impact would appear to be Fort Kamehameha Beach, which lies directly adjacent to the mouth of Pearl Harbor and within close proximity to the Fort Kamehameha WWTP. The Fort Kamehameha WWTP discharges 26Â103 m3 of secondary treated sewage directly into the mouth of Pearl Harbor. The treated wastewater contains 360 mM nitrate+nitrite and 55 mM TP. The high phosphate concentration and anomalously low ammonium+ammonia concentra- tion at Fort Kamehameha Beach probably re¯ect microbial uptake of ammonium stimulated by the phosphate in the wastewater and the fact that the N/P ratio in the wastewater is less than half the Red®eld ratio. Hamilton, Singer, and Waddell (1996, p. 1) have noted that, ``Circulation patterns in Mamala Bay are complex, with substantial changes occurring over short distances and times.'' While agreeing with this assessment, Blumberg and Connolly (1996, p. 2) concluded that, ``The mean ¯ow is typically westward at all depths.'' This latter conclusion is not surprising, since the Hawaiian Islands lie within the North Equa- torial Current. Since roughly 65% of the stream ¯ow and groundwater seepage to Mamala Bay enters Pearl Harbor, one would expect Fort Kamehameha Beach and beaches immediately to the west of the mouth of the harbor to experience the greatest impact from that fresh water. These considerations probably explain the low salinities and high silicate concentrations at Fort Kamehameha and Ewa Beach. Following this line of reasoning, one would expect the ammonium+ammonia con- centration at Ewa Beach to be low, since the out¯ow from Pearl Harbor contains only about 0.4 mM ammonium+ammonia. Obviously this expectation is not 16 E.A. Laws et al. / Marine Environmental Research 48 (1999) 1±21 supported by the experimental data. The ammonium+ammonia concentration at Ewa Beach was higher than at any of the other beaches. Relevant to this observation is the fact that during the time of this study a block of approximately 150 homes within about 200 m of the Ewa Beach shoreline were not connected to the Honolulu municipal sewer system but instead were serviced by cesspools. Seepage from those cesspools probably accounts for the high ammonium+ammonia concentrations at Ewa Beach. Of the 10 beaches studied, nutrient concentrations at one (Keehi Lagoon) seem clearly impacted by stream runo, at another (Fort Kamehameha) by discharges of treated wastewater, and at a third (Ewa Beach) perhaps by cesspool seepage. The remaining beaches show very little impact from fresh water on the concentrations of either dissolved or particulate materials. However, there appear to be some very signi®cant changes in the composition of the particulate material. Most striking is the steady increase in the d15N of the PN from Diamond Head to Barbers Point. The approximate twofold increase in the d15N of the PN from Diamond Head to Barbers Point and the absence of any corresponding change in the concentration of PN can be readily understood with the simple model shown in Fig. 10. According to this model the coastal zone receives an input of fresh water from the land and exchanges water with the oshore ocean. The net exchange with the oshore ocean is assumed to exactly balance the input of fresh water from the land and, hence, the westward ¯ux of water in the coastal zone remains constant. The rate of change of the quantity of any substance in the coastal zone is then given by the equation:
d VX F X X F X F X F F X; 1 dt P DH BP L L O O L O where X is the average concentration of the substance in the coastal zone box; V is the volume of the coastal zone box; XDH and XBP are the concentrations of X at Diamond Head and Barbers Point, respectively; XL and XO are the concentrations of X in land runo/seepage and the oshore ocean, respectively; and FP, FL, and FO
Fig. 10. Box model of water movement in the coastal zone used to derive Eq. (1). Arrows indicate direc- tion of water movement. FP, FL, and FO are the ¯uxes of water (volume/time) parallel to the coastline (Diamond Head to Barbers Point), from the land to the coastal zone box, and from the ocean into the coastal zone box, respectively. E.A. Laws et al. / Marine Environmental Research 48 (1999) 1±21 17 are the ¯uxes of water (volume/time) parallel to the coastline (Diamond Head to Barbers Point), from the land to the coastal zone box, and from the ocean into the coastal zone box, respectively. Under steady state conditions Eq. (1) equals zero, and the right-hand side can be rearranged to give:
fP XDH XBP 1 fOXL fOXO X; 2
where fO and fP are dimensionless numbers equal to FO/(FO+FL) and FP/(FO+FL); clearly 0 1 fOXL fOXO X: 3