Open Access
How to translate text using browser tools
1 October 2017 Impacts of Endangered Seabirds on Nutrient Cycling in Montane Forest Ecosystems of Hawai‘i
Julia A. Rowe, Creighton M. Litton, Christopher A. Lepczyk, Brian N. Popp
Author Affiliations +

Allochthonous nutrient flow from marine sources via seabirds to the terrestrial habitats where they nest can impact resident organisms and neighboring ecosystems. Seabird populations are decreasing both in Hawai‘i and globally, yet little is known about what is being lost from the ecosystems where they traditionally nested in large numbers. Given the marked decline in seabirds, we hypothesized that current sparsely populated seabird colonies in wet montane ecosystems of Hawai‘i contribute minimally to nutrient availability, but that this small contribution should still be reflected in vegetative uptake of soil N and in plant community composition. Soil nutrient availability on Kaua‘i was assessed using ion-exchange resin probes. Plant and soil uptake of marine-derived nitrogen was determined using δ15N values in soil and foliage of the two dominant species using a two-end member N isotope mass balance mixing model. To determine if the added nutrients impacted the plant community, we also compared canopy cover (total and by dominant species) and species richness between treatments. Soil in seabird areas had more available ammonium, but nitrate and total inorganic N did not differ between sites. The dominant canopy tree, Metrosideros polymorpha, derived 28% of foliar N from marine sources; this value was 15% for the dominant understory plant, Dicranopteris linearis. Plant species composition was not influenced by presence or absence of seabirds. Because N plays a large role in net primary productivity, use of marine-derived N by native plants under even limited seabird populations is likely important to functioning of these ecosystems.

Globally, seabirds can be important drivers of nutrient cycling via allochthonous deposition of nutrients such as nitrogen (N ) and phosphorus (P) that influence all trophic levels in the ecosystems where they breed and nest (Anderson and Polis 1999, Mulder and Keall 2001, Whelan et al. 2008, Towns et al. 2009, Grant-Hoffman et al. 2010). Guano naturally contains N and P, nutrients that commonly limit ecosystem processes in forest ecosystems globally (Hutchinson 1950, Elser et al. 2007). Guano composition for Whitecapped Noddies (Anous minutus) and Wedgetailed Shearwaters (Puffinus pacificus) is approximately 7.3% N and 1.5% P (Smith and Johnson 1995), but N content of guano increases to 14.9%–28.6% for pelagic seabirds such as Northern Fulmar (Fulmarus glacialis) and Manx Shearwater (Puffinus puffinus) (Burger et al. 1978, Bird et al. 2008). As such, seabird defecation can increase available pools of N and P in terrestrial ecosystems by 100 and 400 times, respectively (Mulder et al. 2011). In turn, marine-derived nutrient subsidies affect a suite of processes including primary productivity (Mulder and Keall 2001), plant community composition (Anderson et al. 2008, Mizota 2009), and the population size of top predators (Rose and Polis 1998).

The role of seabird guano in the transport of nutrients from marine sources to land has been studied in many coastal ecosystems, but tropical montane ecosystems where many burrowing seabirds nest have been poorly studied (Mulder et al. 2011). Nutrient levels in highly weathered tropical ecosystems can be low, particularly for P, and with high rainfall and steep slopes N and P can be readily lost through leaching and erosion (Ehleringer et al. 1986, Posada and Schuur 2011). In these nutrient-poor ecosystems, the addition of guano rich in N and P may be vital to maintaining or enhancing plant communities and ecosystem process rates (Mizota 2009).

Seabirds face many challenges including habitat loss, introduction of predators, collision with man-made structures, light pollution, toxins, change in prey availability, and poisoning (Millenium Ecosystem Assessment 2005, Hebshi et al. 2008, Duffy 2010, Griesemer and Holmes 2011, Loss et al. 2012, 2015, Wiley et al. 2013). These challenges have typically led to severe population declines and in some cases extinctions, resulting in reduced nutrient inputs to the terrestrial habitats that seabirds traditionally occupied. In Hawai‘i, seabirds that nest in montane forest ecosystems have experienced severe population declines, with some seabird species extirpated from the Hawaiian Islands and others remaining in greatly reduced ranges and numbers (Olson and James 1994, Burney et al. 2001, Hearty et al. 2005). The decline of nutrient flux from marine to terrestrial ecosystems due to the reduction of seabirds on the main Hawaiian Islands may have serious effects on plant community dynamics and ecosystem processes. Thus understanding how native plants utilize nutrient subsidies is important to inform the conservation and restoration of native habitats.

Two seabird species that have experienced drastic population declines in the montane regions of Hawai‘i historically, as well as in the last 25 yr, are the Newell's Shearwater (Puffinus newelli ) and Hawaiian Petrel (Pterodroma sandwichensis), which are federally listed as threatened and endangered, respectively (James 1991, Burney et al. 2001, Griesemer and Holmes 2011). Formerly numerous and widespread, these two species are currently limited to remote colonies in hard-to-access locations due to loss of habitat and an increase in introduced predators. We sought to determine whether Newell's Shearwaters and Hawaiian Petrels in low numbers still influence soil nutrient availability and plant nutrient uptake. Specifically, we sought to answer three primary research questions. First, do low numbers of seabirds increase the availability of macronutrients and micronutrients in wet montane forest soils? Second, do the dominant plants in wet montane ecosystems utilize marine-derived N, and if so to what extent? Third, does avian nutrient subsidy influence plant species composition? We hypothesized that (1) soil micronutrient and macronutrient availability would be higher around the seabird colonies than in areas without seabirds but only minimally given the greatly reduced seabird populations (Wainright et al. 1998, Liu et al. 2006); (2) δ15N values would be higher in soil and foliage at seabird nesting sites, reflecting a marinederivednutrient subsidy (Caut et al. 2012, Kazama et al. 2013); and (3) the plant community composition in seabird plots would be biased toward nitrophilic plants adapted to high nutrient levels (Vitousek and Farrington 1997, Martinelli et al. 1999, Bond et al. 2010).

materials and methods

Study Site

Though greatly reduced from historic levels, the island of Kaua‘i is home to the densest populations of montane nesting seabirds in the Hawaiian archipelago. Study sites were located in the montane forests of Upper Limahuli Preserve and Hono O Nā Pali, Kaua‘i. We considered two treatment types in each of these areas: active seabird colonies and nonseabird areas (areas without current seabird colonies and with no evidence of recent nesting) (Figure 1). Notably, it is likely that seabirds historically nested in most, if not all, montane areas in the past, but the control sites contained no burrows (new or old), bird sign (e.g., feathers or guano), or records of bird activity since 2006, when the Kaua‘i Endangered Seabird Restoration Project began working in the area. Furthermore, the density of seabirds in the most heavily used areas is only 0.04 burrows m-2; this is low for colonial nesting seabirds such as shearwaters and petrels, which are known to nest up to 0.76 burrows m-2 for Gray-faced Petrel (Pterodroma gouldi ) (Whitehead et al. 2014) and 0.08 burrows m-2 for Cook's Petrel (Pterodroma cookii ) (Rayner et al. 2007).

We established 24 plots on ridgetops: nine seabird and four nonseabird plots in Upper Limahuli Preserve and eight seabird and three nonseabird plots in Hono O Nā Pali (Figure 1). On the geologically older Hawaiian island of Kaua‘i, ecosystem processes are typically limited by P availability (Crews et al. 1995), which is expected to apply to our study sites as well. Seabird and nonseabird plots were selected opportunistically in areas with and without seabirds, respectively. Steep slopes, lack of helicopter landing locations, and low seabird numbers made random or uniform plot selection unrealistic. Sample size for isotopic comparisons was based on an a priori power analysis (Faul 2009). Effect size was set at 1.27 based on published soil and Metrosideros polymorpha δ15N values (Vitousek et al. 1989); error probability (α) was set at .05, power (1 - β error probability) was .95, and total sample size was calculated as n = 24. Each plot was 5 m in diameter and established in different seabird subcolony clumps. The Upper Limahuli Preserve ranges in elevation from 750 to 980 m, with plots located above 800 m. Hono O Nā Pali is located above 1,200 m, and plots were situated between 1,210 and 1,287 m.


Location of plots in Upper Limahuli Preserve (ULP) and Hono O Nā Pali (Pihea), Kaua‘i, Hawai‘i. Seabird plots are represented with a triangle, and nonseabird plots are represented with a circle. Image courtesy of National Geographic, Esri, DeLorme, HERE, UNEP-WCMC, USGS, NASA, ESA, METI, NRCAN, GEBCO, NOAA, iPC.


All soils in the two study areas were surveyed by reconnaissance survey. The soils in Upper Limahuli are classified as Alakai mucky peat and Waialeale mucky silty clay loam, and the soils in Hono O Nā Pali are classified as rough mountainous land (Soil Report for Island of Kaua‘i, Hawai‘i 2014). Alakai mucky peat taxonomic classification is clayey, ferrihumic, dysic, isomesic Terric Haplosaprists with pH values typically less than 4.0. Waialeale mucky silty clay loam is classified as veryfine, isotic, isothermic Typic Epiaquods with pH values commonly less than 4.4. Based on similar topography and rainfall, it is likely that the rough mountainous land in Hono O Nā Pali has the same classification as the soils in Upper Limahuli.

All plots were located in wet montane forest with the majority of vegetation being native and dominated by Metrosideros polymorpha (‘öhi‘a) in the canopy and the false staghorn fern Dicranopteris linearis (uluhe) in the understory. No nitrogen-fixing plants were present, although there is evidence of N-fixation in the litter of D. linearis (Russell and Vitousek 1997). Mean annual precipitation at both sites ranges between 2,500 and 3,000 mm with rain occurring throughout the year (Giambelluca et al. 2013). Mean annual temperature is 13°C, with warmest temperatures occurring in August and September (Juvik and Juvik 1998). Both study sites are remote with relatively intact forests exposed to introduced pigs (Sus scrofa), cats (Felis catus), mice (Mus musculus), and two rat species (Rattus exulans and Rattus rattus). In recent years, control and restoration measures have been implemented, including ungulate-proof fencing and pig removal in Upper Limahuli, and trapping for invasive animals and nonnative plant removal in both sites (Jon-Carl Watson, Limahuli Preserve Operations Manager at National Tropical Botanical Garden, pers. comm.).

Available Soil Nutrients

To assess the availability of inorganic soil macronutrients and micronutrients, plant root simulator (PRS) probes (Western Ag Innovations Inc., Saskatchewan, Canada) were deployed in the top 5–10 cm of mineral soil. The PRS probes consist of separate anion and cation exchange membranes that assess nutrient supply rates by continuously absorbing charged ions over the period that they are in the soil. Nutrients indexed represent the bioavailable, labile, inorganic pools in the forms of NO3-, NH4+, H2PO43-, SO42-, K+, Mg2+, and Ca2+ (Johnson et al. 2005, Meason and Idol 2008, Beyene and Katzensteiner 2011). To account for soil heterogeneity, eight pairs of PRS probes were deployed per plot (192 total) during peak to late seabird breeding season (8 September–9 October 2013). The probes were retrieved after 4 weeks, when they were presumed to have reached a dynamic equilibrium (Meason and Idol 2008, Beyene and Katzensteiner 2011). Probes were rinsed with deionized water to remove roots and soil and shipped to Western Ag Innovations Inc. for extraction and analysis. Nutrients were extracted by shaking the probe in 35 ml of 0.5 mol liter-1 HCl for 1 hr to remove > 95% of sorbed ions from the membrane. Concentrations of NH4+, NO3-, and PO4 were then analyzed using colormetic analysis with a Technicon AutoAnalyzer, and K+, Ca2+, and Mg2+ were determined using an inductively coupled plasma spectrometer (PerkinElmer Optima 3000-DV ICP, Norwalk, Connecticut) (Johnson et al. 2005, Meason and Idol 2008).

Isotopic Analysis

Nitrogen isotopic ratios were determined for the top 10 cm of mineral soil, sunlit foliage of both M. polymorpha and D. linearis in seabird and nonseabird plots at both sites as well as seabird guano. Samples were composited (five soil samples and five samples each of M. polymorpha and D. linearis foliage) in each plot. Soil samples were collected using a 1.27 cm diameter soil corer, and sunlit, live leaves were collected from the newest fully mature cohort. Hawaiian Petrel and Newell's Shearwater guano samples were opportunistically collected from field sites during the 2014 breeding season. Twelve relatively fresh guano samples were placed in plastic bags (Ziploc) and kept in a cooler on ice until they could be frozen. Organic materials were handpicked out of guano samples before analysis. Freshness could not be determined, but the high rainfall at these sites should preclude guano from remaining on the ground for more than 2 days of rain. The mean carbon:nitrogen ratio for the sampled seabird guano was 1.09 ± 0.21 (n = 8), and δ15N = 8.23‰ ± 1.68‰ (n = 8).

Soil and foliar samples were dried at 70°C, sieved through a 2 mm mesh, homogenized, and powdered in a ball mill (Carter and Gregorich 2006). Guano samples were freezedried and ground using a mortar and pestle. The isotopic composition of all samples was analyzed at the University of Hawai‘i at Mānoa Biogeochemical Stable Isotope Facility using a continuous flow mass spectrometer (ThermoFinnigan Deltaplus XP) coupled with an elemental analyzer (Costech ECS 4010) via a Conflo IV interface. Nitrogen isotopic compositions are expressed as δ15N values in ‰ relative to air:

where Rsample is the isotopic ratio (15N/14N ) of the sample and Rstandard is the 15N/14N of air. Accuracy and precision of δ15N values were <0.2‰, as determined from reference materials analyzed every 10 samples (glycine and a tuna muscle homogenate with δ15N values of 11.25 ± 0.04‰ and 12.97 ± 0.06‰, respectively, which were determined by extensive characterization using National Institute of Standards and Technology certified reference materials, and their δ15N values were verified independently in other laboratories). Percentage N (%N ) was determined for soil and foliage samples from the results of isotopic analysis.

A two-component N isotope mass balance mixing model was used to determine the proportion of marine-derived nutrients in the top 10 cm of soil and the foliage of the two dominant plants in the plots (Phillips et al. 2005, West et al. 2006). As such, we used the mixing model equation from Dawson et al. (2002):

where δ T is the total sampled isotopic value, δA and δ B are the two source values, and ∫A is the portion of the total value that is provided by source A (Dawson et al. 2002). The nonseabird plot data were used to establish terrestrial δ 15N end member, and the δ 15N value of seabird guano was used as the marine end member. Solutions from the mixing model provide the percentage of total N in soil and foliage of M. polymorpha and D. linearis in the seabird plots that was marinederived(Fry 2006, Phillips 2012, Dawson et al. 2002, Boecklen et al. 2011). Uncertainty in the marine-derived fraction of total N was determined by propagation of error using the analytical solution of the partial differentiation of general Taylor series approximation determined using the two-component stable isotope mixing model (see Gelwicks and Hayes 1990, Phillips and Gregg 2001).

Vegetation Community Assessment

Data were collected in the same 5 m plots discussed earlier in Upper Limahuli and Hono O Nā Pali in collaboration with the Kaua‘i Endangered Seabird Restoration Project. The following variables were quantified: percentage cover for each plant species taller than 2 m, percentage cover for each species shorter than 2 m, average vegetation height, and canopy cover. Species richness was measured by number of species per plot. To measure proportional diversity we used the Shannon index (H′): H′ = Σ (pl)\ln pl I, where (pl) is the proportion of the total number of individuals in the population that are in species “l” (Stirling et al. 2001). Percentage cover of each species was used in lieu of number of individuals.

Statistical Analyses

Statistical analyses were performed in SPSS 22 (SPSS Inc. 2007). Levine's test was used to assess homogeneity of variance in soil nutrient concentration as well as the plant community composition. Of all nutrients analyzed (%N, total inorganic N, NO3-, NH4+, SO42-, PO43-, and Ca+2), only NH4+ did not pass Levine's test and these data were log10 transformed for analysis. All plant community composition data passed Levine's test. One-way analysis of variance (ANOVA), with significance set at α = .05, was used to test for differences in available soil nutrients as well as relative differences in δ15N values in soil and M. polymorpha and D. linearis foliage between treatments. All results are presented as means ± SE, unless otherwise noted. We used t tests to determine differences between treatments for %N in M. polymorpha and D. linearis. For plant community composition, we analyzed canopy cover, average vegetation height, total Metrosideros cover, and total number of species present to look for differences between treatments, also using t tests.


Mean and Standard Error (SE) Inorganic Soil Nutrient Availability for Seabird and Nonseabird Plots (µg/10 cm-2/4 weeks)



Available Soil Nutrients

Across all measured inorganic soil nutrients, only the concentration of NH4+ showed higher values in the seabird plots compared with nonseabird plots (F = 4.74; df = 1, 21; P = .04) (Table 1). Total inorganic N availability did not differ between treatments, largely because the availability of NO3- was nearly identical between seabird and nonseabird plots. In addition, PO43-, Ca+2, and SO42- were slightly, but not significantly, lower in seabird plots (Table 1).

Isotopic Analysis

Foliage of M. polymorpha had significantly higher δ15N values in the seabird than in the nonseabird plots (F = 5.07; df = 1, 21; P = .036) (Table 2). Although δ15N values in soil and D. linearis leaves between the two treatments were not statistically different (F = 2.78; df = 1, 21; P = .11; F = 2.82; df = 1, 21; P = .11, respectively), there was a positive trend toward increasing δ15N values in the seabird plots compared with the nonseabird plots (Table 2). Results from the mixing model indicated that 32% of the total soil N was derived from seabirds in the seabird plots. Foliar N of M. polymorpha in seabird plots was 27.9% (8% SE) from a marine source, and D. linearis foliage contained 16.9% (0.08 SE) N from a marine source (Table 3). However, %N did not differ between seabird and nonseabird plots (Table 4) for soil (t = 0.81, df = 32, P = .43), M. polymorpha foliage (t = 0.17, df = 20, P = .26), or D. linearis foliage (t = 0.92, df = 21, P = .37).


Mean and Standard Error (SE) of δ15N Values of Soil and Two Dominant Plant Species in Seabird and Nonseabird Plots on Kaua‘i



Mean Proportions and Propagated Error for a Two-Source, One-Isotope (δ15N ) Model for Soil and Two Dominant Plant Species in Seabird and Nonseabird Plots P = .26), or D. linearis foliage (t = 0.92, df = 21, P = .37).



Mean %N in Soil, M. polymorpha, and D. linearis


Plant Community Composition

Seabird and nonseabird plots had similar species composition and vegetation structure. Specifically, canopy cover (t = -0.13, df = 22, P = .21), average vegetation height (t = -0.20, df = 22, P = .84), total M. polymorpha cover (t = -1.35, df = 22, P = .19), total species recorded (t = 0.48, df = 22, P = .96), and H′ (t = 0.038, df = 22, P = 0.97) were all similar between treatments.


We found support for the first hypothesis because more inorganic N in the form of NH4+ was found where seabirds were present. In support of the second hypothesis, 32% of soil nitrogen was from a marine source, and that marine-derived N accounted for 17%–28% of foliar N in the two dominant native plants in the study system. However, we did not find evidence to support the third hypothesis: that plant community composition would differ with and without allochthonous input of nutrients by seabirds.

Microbial processes can affect the nitrogen isotopic composition of plant-soil systems (Amundson 2003, Szpak 2014). Nitrogen in guano is deposited primarily as uric acid (C5H4O3N4) (Bird et al. 2008), and microbes mineralize this organic N into NH4+ and NO3 -. Inorganic nitrogen as NH4+ and NO3- is available to most plants in soil solution, making this a potentially important addition to the ecosystem. Ammonium can be lost from the system through oxidation (the first step in soil nitrification), which can lead to N loss as nitrates are leached from the soil or through ammonia volatilization. Nitrogen transformations such as ammonia volatilization and denitrification can affect the δ15N values of plants and soil because there are large nitrogen isotope fractionations associated with these processes. We do not, however, believe that nitrogen loss in the systems studied was the main cause of 15N enrichment in seabird sites. The soil types in these areas are acidic (pH less than 4.5), and under these acidic conditions ammonia would be protonated so that the dominant form of reduced N would be ammonium. In Hawaiian rain forests in regions with mean annual rainfall exceeding ∼2,500 mm, soil microbial denitrification completely consumes nitrate in local soil environments, preventing expression of the isotope effect associated with denitrification (Houlton et al. 2006). Under these conditions δ15N values of soils converge on the δ15N values of the N input. We anticipate that these potential losses of N may make seabird subsidies more critical to the ecosystem than currently understood.

The impacts of seabirds on soil and plant characteristics vary across systems and depend at least partially on the life history of the seabirds in question. Durrett et al. (2014) found that trees and shrubs differed in their response to the addition of marine nutrients, with trees increasing slowly in foliar %N and δ15N values with increasing population density while shrubs showed a strong positive response at low densities and negative responses at higher densities. In ecosystems that are N and /or P limited, the addition of guano may increase primary productivity and select for fastgrowingplants that can take advantage of periodic resource subsidies, or seabirds that nest in high density may make the soils toxic to plants with the excessive addition of nutrients (Wainright et al. 1998, Anderson and Polis 1999, Kolb et al. 2010). In similar studies based on coastal seabird colonies, nutrient pulses caused changes that were observed through the food web and even back into the marine environment. These changes included increased plant productivity, increased plant predation, increased arthropod and lizard density, and increased coastal nutrient influx into the nearshore environment leading to increased plankton growth (Barrett et al. 2005, Spiller et al. 2010, McCauley et al. 2012). On Kaua‘i, Newell's Shearwater and Hawaiian Petrel nest colonially but in low densities in burrows. Because burrow-nesting birds are not observed to reach densities that cause nutrient toxicity to plants, seabirds likely played a larger supporting role historically in nutrient cycling in this ecosystem when population numbers were higher.

We estimated the quantity of nutrients potentially added to our study sites by guano input. The two sites total about 200 ha, 160 ha in Upper Limahuli and 40 ha in Hono O Nā Pali, including nonseabird areas as well as areas where the seabirds are nesting. We have no official estimates of seabird population numbers in these areas, but using Griesemer's island-wide estimates for Newell's Shearwater as well as estimates from the field, we estimated 500 pairs (± 250) of Newell's Shear water and Hawaiian Petrel combined between the two sites, with proportionally more Newell's Shearwater in Upper Limahuli and proportionally more Hawaiian Petrels in Hono O Nā Pali (Griesemer and Holmes 2011). No studies have been conducted to indicate how much Newell's Shearwater or Hawaiian Petrel eat or excrete per day. However, the wandering albatross weighs approximately 10 kg and was found to consume 2 kg of food per day (Salamolard and Weimerskirch 1993). The average weight for Newell's Shearwater and Hawaiian Petrel is 0.4 kg (Ainley et al. 1997, Judge et al. 2014). Therefore, because small organisms require more food per unit body mass than larger organisms, we used the allometric relationship:

to estimate food intake, assuming birds are in homeostasis. Doing this we estimated that seabirds in our study sites consume 0.2 kg of marine-based food per day. Based on the seabird intake and calculations of guano production in dovekies (Gabrielsen et al. 1991), we estimated that 500 seabirds could produce 98.5 kg of guano day-1, or 73.9 kg ha-1 year-1. Not all of the guano would end up in the montane ecosystem because one bird of a pair would likely be out to sea, so half of this estimate is 37 kg guano ha-1 year-1. Estimating the nitrogen content of the guano at 22% yields 16 kg N ha-1 year-1 (Bird et al. 2008). For comparison, total N deposition from precipitation was measured as 1 kg N ha-1 year-1 at a site on Hawai‘i Island (P. Vitousek 2004). Acacia koa, a dominant native symbiotic N fixing tree, was not present in our sites, but estimates of N-fixation in dense regenerating A. koa stands range from 23 kg N ha-1 year-1 in 5-yr-old stands to 1.5 kg N ha-1 year-1 in 20-yr-old stands (Pearson and Vitousek 2002).

Isotopic Analysis

Seabirds have been shown to increase N levels in soil and surrounding organisms via marine-derived N (Wainright et al. 1998, Wait et al. 2005, Mizota 2009). However, it was previously unknown if this is also the case in wet tropical montane regions characterized by high rainfall, warm temperatures, and low current population densities of seabirds (Martinelli et al. 1999, Garcia et al. 2002). Although marine N was clearly higher in M. polymorpha, the effect of seabird-added nitrogen may be masked in D. linearis by nitrogen fixation that may occur in the litter (Russell and Vitousek 1997). It should also be kept in mind that although the nonseabird sites had no current evidence for nesting, they were likely colonized in the past and may have a legacy of high δ15N values in soil. Thus, using these nonseabird sites likely resulted in an overestimation of terrestrial δ15N end member values in the isotope mass balance mixing model and thus systematically underestimated the proportion of marine-derived N available in the soil and incorporated into the foliage of both studied species.

Vitousek et al. (1989) and Martinelli et al. (1999) measured δ15N of M. polymorpha foliage in nonseabird areas in Hawai‘i across multiple islands and found mean δ15N values of -3.3 (± 2.3‰ SD) and -2.8 (± 2.6‰ SD), respectively. A comparison of their values with this study is complicated by difference in substrate age between sampling sites in these prior studies and ours on Kaua‘i. Martinelli et al. (1999) also measured a δ15N value of -0.5‰ from a single M. polymorpha on Kaua‘i from a nonseabird area. It is unknown if the δ15N value of -0.5‰ is an outlier or representative of M. polymorpha on Kaua‘i. However, the archipelago averages of -3.3‰ and -2.8‰ in nonseabird areas are slightly lower than the nonseabird values that we measured for M. polymorpha (-2.3‰), and this also indicates a potential underestimate of the importance of current marine-derived N as presented here.

Results of the isotopic mixing models indicate that dominant plants in this ecosystem utilize at least some N derived from a marine source, and soil and foliage of both plant species showed marine influence. Although standard error was high, the amount of marinesourced N was higher across all sampled substrates in seabird plots, and δ15N values were significantly higher in M. polymorpha foliage in seabird plots compared with control plots.

Plant Community Composition

None of the plant species composition measures indicated differences between seabird and nonseabird plots. This is likely due not only to the limited amount of nitrogen being added to the seabird sites, but also to the species-depauperate nature of the Islands and isolation of the study sites. There are a limited number of native species on the Islands to populate these areas, and the isolated nature of the study sites means less influence of the nonnative and invasive plant species that occur in high densities in more disturbed areas in Hawai‘i. Invasive plant species were actively managed in Upper Limahuli and pulled opportunistically in Hono O Nā Pali, though the density of invasive plants is low in both sites due to the remote locations. Another reason for the lack of differences in some of the measured variables is that other drivers besides the addition of nutrients by seabirds may be more influential, particularly in lowdensity seabird sites such as ours. Other researchers have found that on high- and medium-density seabird islands, seabirds drove ecosystem properties such as δ15N values (soil and leaf ), soil and leaf N, NH4+, and NO3-. However, in low-density colonies other ecosystem processes drove these ecosystem properties more than the seabirds (Durrett et al. 2014).


Despite being at historically low population densities, seabirds contribute to the ecosystems where they still nest in montane Hawai‘i via marine-derived nutrient deposition. These study sites in Kaua‘i contain some of the last relatively intact tropical montane ecosystems with native seabirds, yet very little research has occurred there. Studies in the arctic and in coastal systems are abundant (Polis and Hurd 1996, Mulder et al. 2011, Gagnon et al. 2013) and generally show seabirds to increase nutrient availability and biodiversity in the arctic (Keatley et al. 2009, Zmudczyńska et al. 2012) and to fertilize or even create toxic conditions in coastal ecosystems depending on seabird density (Kolb et al. 2010, VanderWerf et al. 2014). However, comprehensive studies about how seabirds and their nutrient subsidies impact tropical montane ecosystems are notably lacking (Hawke and Holdaway 2009).

Our control plots may have a historical legacy of seabirds and thus may still contain nutrients from seabirds. Therefore, our estimate of allochthonous nutrient input by seabirds is likely conservative. The influx of N and P may be more important to the resiliency of these ecosystems, especially in the face of climate change and other stressors, than is currently understood (Perry et al. 2010, Doughty et al. 2015). Historically seabirds in the Hawaiian Islands may have played a major role in controlling soil fertility in areas where they nested. The montane forest seabirds in the Hawaiian Islands have been greatly reduced both in number and in range, yet even in vastly reduced numbers these birds impact soil and vegetation nutrient content. The decrease in nutrient flow and its impacts on the ecosystem are vital for restoration project managers to understand as they attempt to rebuild the ecosystem and restore endemic plants and wildlife. We know that seabird declines have sweeping effects throughout the food web in warm temperate New Zealand (Fukami et al. 2006), but the consequences of seabird declines on tropical vegetation dynamics and ecosystem processes are largely unknown.


Special thanks go to Andre Raine and the Kaua‘i Endangered Seabird Recovery Project (KESRP) for logistics in Kaua‘i and a huge amount of support in the field. We thank National Tropical Botanical Gardens for access to the Upper Limahuli Preserve as well as State of Hawai‘i Division of Forestry and Wildlife's Natural Area Reserves System folks for access to Pihea in Hono O Nā Pali. Thanks also to Trevor Joyce for assistance and advice in the field. We also thank the two anonymous reviewers for pushing us to substantially improve our paper. This is School of Ocean and Earth Science and Technology (SOEST) contribution no. 10059.

Literature Cited


Ainley, D. G., T. C. Telfer, and M. H. Reynolds. 1997. Townsend's and Newell's Shearwater (Puffinus auricularis). In A. Poole, ed. The Birds of North America Online. Google Scholar


Amundson, R. 2003. Global patterns of the isotopic composition of soil and plant nitrogen. Glob. Biogeochem. Cycles 17:1031. Google Scholar


Anderson, W. B., and G. A. Polis. 1999. Nutrient fluxes from water to land: Seabirds affect plant nutrient status on Gulf of California islands. Oecologia (Berl.) 118:324–332. Google Scholar


Anderson, W. B., D. A. Wait, and P. Stapp. 2008. Resources from another place and time: Responses to pulses in a spatially subsidized system. Ecology 89:660–670. Google Scholar


Barrett, K., W. B. Anderson, D. A. Wait, L. L. Grismer, G. A. Polis, and M. D. Rose. 2005. Marine subsidies alter the diet and abundance of insular and coastal lizard populations. Oikos 109:145–153. Google Scholar


Beyene, T. F., and K. Katzensteiner. 2011. Assessment of supply of soil nutrients in different land use types using plant root stimulator probes. Page Conference on International Research on Food Security, Natural Resource Management, and Rural Development. University of Bonn, Bonn, Germany. Google Scholar


Bird, M. I., E. Tait, C. M. Wurster, and R. W. Furness. 2008. Stable carbon and nitrogen isotope analysis of avian uric acid. Rapid Commun. Mass Spectrom. 22:3393–3400. Google Scholar


Boecklen, W. J., C. T. Yarnes, B. A. Cook, and A. C. James. 2011. On the use of stable isotopes in trophic ecology. Annu. Rev. Ecol. Evol. Syst. 42:411–440. Google Scholar


Bond, A. L., G. T. W. McClelland, I. L. Jones, J. L. Lavers, and T. K. Kyser. 2010. Stable isotopes confirm community patterns in foraging among Hawaiian Procellariiformes. Waterbirds 33:50–58. Google Scholar


Burger, A. E., H. J. Lindeboom, and A. J. Williams. 1978. The mineral and energy contributions of guano of selected species of birds to the Marion Island terrestrial ecosystem. S. Afr. J. Antarct. Res. 8:59–70. Google Scholar


Burney, D. A., H. F. James, L. P. Burney, S. L. Olson, W. Kikuchi, W. L. Wagner, M. Burney, D. McCloskey, D. Kikuchi, F. V. Grady, R. Gage, and R. Nishek. 2001. Fossil evidence for a diverse biota from Kaua‘i and its transformation since human arrival. Ecol. Monogr. 71:615–641. Google Scholar


Carter, M. R., and E. G. Gregorich, eds. 2006. Soil sampling and methods of analysis. 2nd ed. CRC Press, Boca Raton, Florida. Google Scholar


Caut, S., E. Angulo, B. Pisanu, L. Ruffino, L. Faulquier, O. Lorvelec, J.-L. Chapuis, M. Pascal, E. Vidal, and F. Courchamp. 2012. Seabird modulations of isotopic nitrogen on islands. PLoS ONE 7:e39125. Google Scholar


Crews, T. E., K. Kitayama, J. H. Fownes, R. H. Riley, D. A. Herbert, D. Mueller-Dombois, and P. M. Vitousek. 1995. Changes in soil phosphorus fractions and ecosystem dynamics across a long chronosequence in Hawaii. Ecology 76:1407–1424. Google Scholar


Dawson, T. E., S. Mambelli, A. H. Plamboeck, P. H. Templer, and K. P. Tu. 2002. Stable isotopes in plant ecology. Annu. Rev. Ecol. Syst. 33:507–559. Google Scholar


Doughty, C. E., J. Roman, S. Faurby, A. Wolf, A. Haque, E. S. Bakker, Y. Malhi, J. B. Dunning, and J.-C. Svenning. 2015. Global nutrient transport in a world of giants. Proc. Natl. Acad. Sci. U.S.A. 113:868–873. Google Scholar


Duffy, D. C. 2010. Changing seabird management in Hawai‘i: From exploitation through management to restoration. Waterbirds 33:193–207. Google Scholar


Durrett, M. S., D. A. Wardle, C. P. H. Mulder, and R. P. Barry. 2014. Seabirds as agents of spatial heterogeneity on New Zealand's offshore islands. Plant Soil 383:139–153. Google Scholar


Ehleringer, J. R., C. B. Field, Z. Lin, and C. Kuo. 1986. Leaf carbon isotope and mineral composition in subtropical plants along an irradiance cline. Oecologia (Berl.) 70:520–526. Google Scholar


Elser, J. J., M. E. S. Bracken, E. E. Cleland, D. S. Gruner, W. S. Harpole, H. Hillebrand, J. T. Ngai, E. W. Seabloom, J. B. Shurin, and J. E. Smith. 2007. Global analysis of nitrogen and phosphorus limitation of primary producers in freshwater, marine and terrestrial ecosystems. Ecol. Lett. 10:1135–1142. Google Scholar


Faul, F., E. Erdfelder, A. Buchner, and A.-G. Lang. 2009. Statistical power analyses using G*Power 3.1: Tests for correlation and regression analyses. Behav. Res. Methods 41:1149–1160. Google Scholar


Fukami, T., D. A. Wardle, P. J. Bellingham, C. P. H. Mulder, D. R. Towns, G. R. W. Yeates, K. I. Bonner, M. S. Durrett, M. N. Grant-Hoffman, and W. M. Williamson. 2006. Above- and below-ground impacts of introduced predators in seabirddominatedisland ecosystems. Ecol. Lett. 9:1299–1307. Google Scholar


Gabrielsen, G. W., J. R. E. Taylor, M. Konarzewski, and F. Mehlum. 1991. Field and laboratory metabolism and thermoregulation in dovekies (Alle alle). Auk 108:71–78. Google Scholar


Gagnon, K., E. Rotha¨usler, A. Syrja¨nen, M. Yli-Renko, and V. Jormalainen. 2013. Seabird guano fertilizes Baltic Sea littoral food webs. PLoS ONE 8:e61281. Google Scholar


Garcia, L. V., T. Maran, F. Ojeda, L. Clemente, and R. Redondo. 2002. Seagull influence on soil properties, chenopod shrub distribution, and leaf nutrient status in semi-arid Mediterranean islands. Oikos 98:75–86. Google Scholar


Gelwicks, J. T., and J. M. Hayes. 1990. Carbon-isotopic analysis of dissolved acetate. Anal. Chem. 62:535–539. Google Scholar


Giambelluca, T. W., Q. Chen, A. G. Frazier, J. P. Price, Y. L. Chen, P. S. Chu, J. K. Eischeid, and D. M. Delparte. 2013. Online rainfall atlas of Hawaii. Bull. Am. Meteorol. Soc. 94:313–316. Google Scholar


Grant-Hoffman, M. N., C. P. H. Mulder, and P. J. Bellingham. 2010. Effects of invasive rats and burrowing seabirds on seeds and seedlings on New Zealand islands. Oecologia (Berl.) 162:1005–1016. Google Scholar


Griesemer, A. M., and N. D. Holmes. 2011. Newell's shearwater population modeling for habitat conservation plan and recovery planning. Tech. Rep. 176. The Hawai‘i — Pacific Islands Cooperative Ecosystem Studies Unit and Pacific Cooperative Studies Unit, University of Hawai‘i at Mānoa, Honolulu. Google Scholar


Hawke, D. J., and R. N. Holdaway. 2009. Nutrient sources for forest birds captured within an undisturbed petrel colony, and management implications. Emu 109:163–169. Google Scholar


Hearty, P. J., H. F. James, and S. L. Olson. 2005. The geological context of Middle Pleistocene crater lake deposits and fossil birds at Ulupau Head, Oahu, Hawaiian Islands. Monogr. Soc. Hist. Nat. Balears 12:113–128. Google Scholar


Hebshi, A. J., D. C. Duffy, and K. D. Hyrenbach. 2008. Associations between seabirds and subsurface predators around Oahu, Hawaii. Aquat. Biol. 4:89–98. Google Scholar


Houlton, B. Z., D. M. Sigman, and L. O. Hedin. 2006. Isotopic evidence for large gaseous nitrogen losses from tropical rainforests. Proc. Natl. Acad. Sci. U.S.A. 103:8745–8750. Google Scholar


Hutchinson, G. E. 1950. The biogeochemistry of vertebrate excretion. Soil Sci. 96:596. Google Scholar


James, H. 1991. The contribution of fossils to knowledge of Hawaiian birds. Pages 420–424 in Acta XX Congr. Int. Ornithol. Allen & Unwin, London. Google Scholar


Johnson, D. W., P. S. J. Verburg, and J. A. Arnone. 2005. Soil extraction, ion exchange resin, and ion exchange membrane measures of soil mineral nitrogen during incubation of a tallgrass prairie soil. Soil Sci. Soc. Am. J. 69:260–265. Google Scholar


Judge, S. W., D. Hu, and C. N. Bailey. 2014. Comparative analyses of Hawaiian petrel Pterodroma sandwichensis morphometrics. Mar. Ornithol. 42:81–84. Google Scholar


Juvik, S. P., and J. O. Juvik, eds. 1998. Atlas of Hawai‘i. 3rd ed. University of Hawai‘i Press, Honolulu. Google Scholar


Kazama, K., H. Murano, K. Tsuzuki, H. Fujii, Y. Niizuma, and C. Mizota. 2013. Input of seabird-derived nitrogen into rice-paddy fields near a breeding/roosting colony of the Great Cormorant (Phalacrocorax carbo), and its effects on wild grass. Appl. Geochem. 28:128–134. Google Scholar


Keatley, B. E., M. S. V Douglas, J. M. Blais, M. L. Mallory, and J. P. Smol. 2009. Impacts of seabird-derived nutrients on water quality and diatom assemblages from Cape Vera, Devon Island, Canadian High Arctic. Hydrobiologia 621:191–205. Google Scholar


Kolb, G. S., L. Jerling, and P. A. Hamback. 2010. The impact of cormorants on plantarthropod food webs on their nesting islands. Ecosystems 13:353–366. Google Scholar


Liu, X. D., S. P. Zhao, L. G. Sun, X. B. Yin, Z. Q. Xie, L. Honghao, and Y. H. Wang. 2006. P and trace metal contents in biomaterials, soils, sediments and plants in colony of red-footed booby (Sula sula) in the Dongdao Island of South China Sea. Chemosphere 65:707–715. Google Scholar


Loss, S. R., T. Will, and P. P. Marra. 2012. Direct human-caused mortality of birds: Improving quantification of magnitude and assessment of population impact. Front. Ecol. Environ. 10:357–364. Google Scholar


Loss, S. R., T. Will, and P. P. Marra. 2015. Direct mortality of birds from anthropogenic causes. Annu. Rev. Ecol. Evol. Syst. 46:99–120. Google Scholar


Martinelli, L. A., M. C. Piccolo, A. R. Townsend, P. M. Vitousek, E. Cuevas, W. McDowell, G. P. Robertson, O. C. Santos, and K. Treseder. 1999. Nitrogen stable isotopic composition of leaves and soil: Tropical versus temperate forests. Biogeochemistry (Dordr.) 46:45–65. Google Scholar


McCauley, D. J., P. A Desalles, H. S. Young, R. B. Dunbar, R. Dirzo, M. M. Mills, and F. Micheli. 2012. From wing to wing: The persistence of long ecological interaction chains in less-disturbed ecosystems. Sci. Rep. 2:1–5. Google Scholar


Meason, D. F., and T. W. Idol. 2008. Nutrient sorption dynamics of resin membranes and resin bags in a tropical forest. Soil Sci. Soc. Am. J. 72:1806. Google Scholar


Millenium Ecosystem Assessment. 2005. Ecosystems and human well-being: Biodiversity synthesis. Page Ecosystems, Washington, D.C. Google Scholar


Mizota, C. 2009. Nitrogen isotopic patterns of vegetation as affected by breeding activity of Black-tailed Gull (Larus crassiostris): A coupled analysis of feces, inorganic soil nitrogen and flora. Appl. Geochem. 24:2027–2033. Google Scholar


Mulder, C. P. H., W. B. Anderson, D. R. Towns, and P. J. Bellingham, eds. 2011. Seabird islands: Ecology, invasion, and restoration. Oxford University Press, New York. Google Scholar


Mulder, C. P. H., and S. N. Keall. 2001. Burrowing seabirds and reptiles: Impacts on seeds, seedlings and soils in an island forest in New Zealand. Oecologia (Berl.) 127:350–360. Google Scholar


Olson, S. L., and H. F. James. 1994. A chronology of ornithological exploration in the Hawaiian Islands, from Cook to Perkins. Stud. Avian Biol. 15:91–102. Google Scholar


Pearson, H. L., and P. M. Vitousek. 2002. Soil phosphorus fractions and symbiotic nitrogen fixation across a substrate-age gradient in Hawaii. Ecosystems 5:587–596. Google Scholar


Perry, G. L. W., J. Ogden, N. J. Enright, and L. V. Davy. 2010. Vegetation patterns and trajectories in disturbed landscapes, Great Barrier Island, northern New Zealand. N. Z. J. Ecol. 34:311–323. Google Scholar


Peterson, B., and B. Fry. 1987. Stable isotopes in ecosystem studies. Annu. Rev. Ecol. Syst. 18:293–320. Google Scholar


Phillips, D. L. 2012. Converting isotope values to diet composition: The use of mixing models. J. Mammal. 93:342–352. Google Scholar


Phillips, D. L., and J. W. Gregg. 2001. Uncertainty in source partitioning using stable isotopes. Oecologia (Berl.) 127:171–179. Google Scholar


Phillips, D. L., S. D. Newsome, and J. W. Gregg. 2005. Combining sources in stable isotope mixing models: Alternative methods. Oecologia (Berlin) 144:520–527. Google Scholar


Polis, G. A., and S. D. Hurd. 1996. Linking marine and terrestrial food webs: Alloch thonous input from the ocean supports high secondary productivity on small islands and coastal land communities. Am. Nat. 147:396–423. Google Scholar


Posada, J. M., and E. A. G. Schuur. 2011. Relationships among precipitation regime, nutrient availability, and carbon turnover in tropical rain forests. Oecologia (Berl.) 165:783–795. Google Scholar


Rayner, M. J., M. N. Clout, R. K. Stamp, M. J. Imber, D. H. Brunton, and M. E. Hauber. 2007. Predictive habitat modelling for the population census of a burrowing seabird: A study of the endangered Cook's petrel. Biol. Conserv. 138:235–247. Google Scholar


Rose, M. D., and G. A. Polis. 1998. The distribution and abundance of coyotes: The effects of allochthonous food subsidies from the sea. Ecology 79:998–1007. Google Scholar


Russell, A. E., and P. M. Vitousek. 1997. Decomposition and potential nitrogen fixation in Dicranopteris linearis litter on Mauna Loa, Hawai‘i. J. Trop. Ecol. 13:579–594. Google Scholar


Salamolard, M., and H. Weimerskirch. 1993. Relationship between foraging effort and energy requirement throughout the breeding season in the wandering albatross. Funct. Ecol. 7:643–652. Google Scholar


Schneider, D. C. 2002. Scaling theory: Application to marine ornithology. Ecosystems 5:736–748. Google Scholar


Smith, J. S., and C. R. Johnson. 1995. Nutrient inputs from seabirds and humans on a populated coral cay. Mar. Ecol. Prog. Ser. 124:189–200. Google Scholar


Soil Report for Island of Kaua‘i, Hawai‘i. 2014. Scholar


Spiller, D. A., J. Piovia-Scott, A. N. Wright, L. H. Yang, G. Takimoto, T. W. Schoener, and T. Iwata. 2010. Marine subsidies have multiple effects on coastal food webs. Ecology 91:1424–1434. Google Scholar


SPSS Inc. 2007. SPSS for Windows, version 16.0. SPSS, Chicago. Google Scholar


Stirling, G., B. J. Wilsey, and B. Wilsey. 2001. Empirical relationships between species richness, evenness, and proportional diversity. Am. Nat. 158:286–299. Google Scholar


Szpak, P. 2014. Complexities of nitrogen isotope biogeochemistry in plant-soil systems: Implications for the study of ancient agricultural and animal management practices. Front. Plant Sci. 5:288. Google Scholar


Towns, D. R., D. A. Wardle, C. P. H. Mulder, G. W. Yeates, B. M. Fitzgerald, G. R. Parrish, P. J. Bellingham, and K. I. Bonner. 2009. Predation of seabirds by invasive rats: Multiple indirect consequences for invertebrate communities. Oikos 118:420–430. Google Scholar


VanderWerf, E. A., L. C. Young, S. E. Crow, E. Opie, H. Yamazaki, C. J. Miller, D. G. Anderson, L. S. Brown, D. G. Smith, and J. Eijzenga. 2014. Increase in wedge-tailed shearwaters and changes in soil nutrients following removal of alien mammalian predators and nitrogen-fixing plants at Kaena Point, Hawaii. Restor. Ecol. 22:676–684. Google Scholar


Vitousek, P. M., and H. Farrington. 1997. Nutrient limitation and soil development: Experimental test of a biogeochemical theory. Biogeochemistry (Dordr.) 37:63–75. Google Scholar


Vitousek, P. M., G. Shearer, and D. H. Kohl. 1989. Foliar 15N natural abundance in Hawaiian rainforest: Patterns and possible mechanisms. Oecologia (Berl.) 78:383–388. Google Scholar


Wainright, S. C., J. C. Haney, C. Kerr, A. N. Golovkin, and M. V. Flint. 1998. Utilization of nitrogen derived from seabird guano by terrestrial and marine plants at St. Paul, Pribilof Islands, Bering Sea, Alaska. Mar. Biol. (Berl.) 131:63–71. Google Scholar


Wait, D. A., D. P. Aubrey, and W. B. Anderson. 2005. Seabird guano influences on desert islands: Soil chemistry and herbaceous species richness and productivity. J. Arid Environ. 60:681–695. Google Scholar


West, J. B., G. J. Bowen, T. E. Cerling, and J. R. Ehleringer. 2006. Stable isotopes as one of nature's ecological recorders. Trends Ecol. Evol. 21:408–414. Google Scholar


Whelan, C. J., D. G. Wenny, and R. J. Marquis. 2008. Ecosystem services provided by birds. Ann. N. Y. Acad. Sci. 1134:25–60. Google Scholar


Whitehead, A. L., P. O. B. Lyver, C. J. Jones, P. J. Bellingham, C. J. MacLeod, M. Coleman, B. J. Karl, K. Drew, D. Pairman, A. M. Gormley, and R. P. Duncan. 2014. Establishing accurate baseline estimates of breeding populations of a burrowing seabird, the grey-faced petrel (Pterodroma macroptera gouldi ) in New Zealand. Biol. Conserv. 169:109–116. Google Scholar


Wiley, A. E., P. H. Ostrom, A. J. Welch, R. C. Fleischer, H. Gandhi, and J. R. Southon. 2013. Millennial-scale isotope records from a wide-ranging predator show evidence of recent human impact to oceanic food webs. Proc. Natl. Acad. Sci. U.S.A. 110:8972–8977. Google Scholar


Zmudczyn´ ska, K., I. Olejniczak, A. Zwolicki, L. Iliszko, P. Convey, and L. Stempniewicz. 2012. Influence of allochtonous nutrients delivered by colonial seabirds on soil collembolan communities on Spitsbergen. Polar Biol. 35:1233–1245. Google Scholar


[1] This research was supported by funding from the National Institute of Food and Agriculture, U.S. Department of Agriculture, under Agreement no. 2010-38420-20381; and the College of Tropical Agriculture and Human Resources, University of Hawai‘i at Mānoa via the Hatch (HAW00132H and HAW01127H to C.M.L.) and McIntyre Stennis (HAW00188M and HAW01123M to C.M.L.) Programs.

Julia A. Rowe, Creighton M. Litton, Christopher A. Lepczyk, and Brian N. Popp "Impacts of Endangered Seabirds on Nutrient Cycling in Montane Forest Ecosystems of Hawai‘i," Pacific Science 71(4), 495-509, (1 October 2017).
Accepted: 1 April 2017; Published: 1 October 2017
Back to Top