Abstract

Quantitative benthic samples were taken during spring and summer at three locations on the Louisiana continental shelf from 2004 to 2012 to assess the influence of hypoxia on the mean sizes (wet weight) of polychaete annelid worms. While the mean body size over the entire study of 64 samples was 3.99 ± 4.66 mg wet weight per individual, the mean ranged from 2.97 ± 2.87 mg during consistently hypoxic conditions (<2 mg/L) to a high of 7.13 ± 7.60 mg () under oxic conditions (>2 mg/L). The variations in size within assemblages were estimated from conventional biomass size spectra (BSS) and normalized biomass size spectra (NBSS) across a broad range of oxygen concentrations. The decline in size was due to the elimination of large species under hypoxic conditions (<2 mg/L), not a reduction in size within species. At “severe” levels of hypoxia (<1 mg/L), the smallest species also declined in abundance, whereas the ubiquitous “medium-sized” Paraprionospio pinnata flourished. These results suggest that there will be enhanced selection for small sizes and species with enlarged branchial palps such as those in P. pinnata if, as predicted, hypoxia becomes more commonplace in time and space worldwide.

1. Introduction

The northern Gulf of Mexico continental shelf just to the west of the Mississippi River Delta experiences widespread hypoxia seasonally [1, 2]. Hypoxia (dissolved oxygen < 2 mg/L) results from two synergistic factors: the summer vertical stratification of the water column that prevents the mixing of oxygen into deep water [35] and the demand for oxygen by the decomposition of organic matter below the pycnocline [1, 6]. Since the benthic infauna is less mobile than demersal crustaceans and fishes, they are “trapped” in deep and near-bottom water and thus susceptible to low oxygen stress. The negative effects of hypoxia on coastal benthos (decline of biodiversity, abundance, biomass, and bioturbation) have been broadly recognized [79]. Many of the few species that can survive intermittent low oxygen stress are small polychaete annelid worms, the most common taxon among all benthic macrofauna [10, 11], but we know little about how these organisms are able to persist under such stress [12].

Body size is a quantifiable trait that reflects a species’ evolution and ecological relationships [1216]. Food availability [17], temperature [18], depth [19, 20], latitude [21], and even oxygen levels [22], among others, are assumed to affect body size. The potential relationship between oxygen and body size has been investigated across a wide range of taxa [2326].

It is reasonable that smaller animals with large surface area to volume ratios would be better able to absorb oxygen in low concentrations [27]. It has also been suggested that small organisms may be less vulnerable to disturbance because of their high reproductive rate and diminished demand for resources [28]. Thus we are suggesting that small size in the polychaetes may be one trait that allows them to persist during hypoxic conditions. To pursue this idea, we have generated a set of data on polychaete sizes encountered across an extreme range of interannual and seasonal oxygen concentrations.

The relationship between size and low oxygen stress is not necessarily simple and direct. The stress could eliminate large species in general, leaving only small species to survive. Alternatively, size could decline within a species, or, given wide and divergent functional morphology within the entire taxa, the survival may be independent of size. The survivors may possess some morphological or physiological talents that make them capable of utilizing oxygen at very low concentrations [29, 30].

Generation of biomass size spectra (BSS) is a simple, timesaving way to compare the structure of assemblages or populations [31], a method first proposed by Sheldon and Parsons [32] for particulate matter. Size spectra of populations [33] and assemblages [34] have historically provided significant insights into growth rates and secondary production, especially within the benthos. Size spectra have provided insights on survival in stressful habitats, such as oxygen minimum zones (OMZ) [34, 35].

Variations in configuration and shape of simple size distribution spectra are sometimes difficult to quantify and compare [36]. To obviate this problem, Platt and Denman [37, 38] suggested constructing “normalized biomass size spectra” (NBSS) that would facilitate comparisons of size distributions between different communities [36]. We have therefore included NBSS plots as an alternative option for characterizing size structure within our assemblages.

For most animal communities, NBSS is a linear plot of the total biomass as a function of narrow size categories of NBSS, from small to large organisms on the -axis. A negative slope represents abundant but small organisms with relatively high total biomass down to large organisms that are fewer in number and lower in total biomass. The greater the dominance of the small size classes, the steeper the negative slope.

The slope and intercept of an NBSS regression line have been utilized previously in comparisons of the size structure of communities from different habitats [14, 34]. For example, a steep slope of the regression illustrates an assemblage dominated by small-sized organisms compared to a modest slope that reflects a higher frequency of larger individuals. Thus, according to the Pearson and Rosenberg [39] model, a community dominated by small opportunists following a disturbance would have a steep slope but the slope would rise during recovery as the small, fast growing opportunists are replaced by larger, more slowly growing species. In general, the intercept and the height of the regression line above the -axis are a function of total biomass. It has been suggested that the departure from a statistically significant negative regression is a departure from steady state [31, 36, 40, 41].

In this paper, we will (1) compare bottom water oxygen concentrations in spring and summer at three standard locations along a 50 km long transect within a narrow depth range (11 to ca. 20 m) in the northern Gulf of Mexico; (2) construct polychaete size spectra from benthic samples subjected to different oxygen concentrations at each of the locations; and (3) compare biomass size spectra (BSS) and normalized biomass size spectra (NBSS) among assemblages at different oxygen concentrations. This comparison is intended to provide insights into the influence of seasonal hypoxia on polychaete body size structure under (a) oxic (>2 mg/L), (b) hypoxic (<2 mg/L), and (c) severely hypoxic (<1 mg/L) conditions.

2. Methods

2.1. Sampling Strategies

The sampling in this study began in 2004 and extended somewhat irregularly until 2012. Stations were placed close to the Mississippi River Birdfoot Delta on the continental shelf (Figure 1, A) and extended to the west between the 10 to 20 m isobaths (Figure 1, C). This pattern followed sampling sites occupied by the Mechanisms Controlling Hypoxia (MCH) project that has attempted to make consistent, repeated sampling of the same locations from early spring, with the onset of stratification, through to late summer, when hypoxia and stratification are often most intense (Figure 1 and Table 1). The locations were designed to compare the intensity and the causes of hypoxia at what has been termed the “brown” water zone close to the river mouth (A), the “green” water zone with high concentrations of plant biomass just to the west of the brown water (B), and the “blue” water zone further to the west of B, where the water column is comparatively more transparent (C); these zones were more or less evenly placed along the extension of the river plume as it flows and disperses from east to west [42], as proposed by Rowe and Chapman [5]. The purpose of this abbreviated pattern in time and space was to observe the setup of the physical stratification in spring that prevents deep water from becoming oxygenated [43, 44] and thus leads to eventual hypoxia by summer. A somewhat similar line of sampling locations has been used to assess the effects of hypoxia on polychaete trophic structure [45], biogenic structures within the sediments [46], and sediment erodibility [47], but in slightly deeper water.

Oxygen concentrations in the water column were profiled using a SeaBird 911 CTD from the surface to 1 to 2 meters above bottom. The CTD was equipped with a SB43 oxygen electrode, as well as conductivity and thermistor probes. A rosette of 12 Niskin bottles caught individual water samples that were analyzed for plant pigments, salinity, dissolved oxygen, and inorganic nutrients [48]. Near-bottom (ca. < 0.5 m) water samples were taken using a four-bottle sampler that tripped when a trigger weight encountered the mud. MCH data are available through NODC under Accession number 0088164.

Three subcores (pseudoreplicates) with a diameter of 13.9 cm each were taken from a single 0.2 m2 GOMEX box corer at each location on each cruise. These subcores were initially incubated for several hours in a water bath aboard ship to determine sediment community oxygen consumption and nutrient regeneration [49, 50]. The entire volume of sediment of each subcore was then sieved through 0.5 mm stainless steel mesh to remove the macrofauna. In 2012 three replicate 0.1 m2 GOMEX box cores were taken at sites B, C, and AB, an intermediate site added in 2009 and 2012 (Figure 1 and Table 1). The surface 15 cm of each of these was also sieved through 0.5 mm stainless steel mesh. Following the sieving at sea, all samples were fixed aboard ship with 10% buffered formalin in filtered sea water.

2.2. Sample and Data Processing

In the lab, the sieved and formalin-fixed biota of the subcore samples were stained with several mL of 5% Rose Bengal for 24 hours and the stained invertebrates were sorted to the lowest possible taxonomic level using a dissecting microscope. Samples were preserved in 70% ethyl alcohol following sorting. The sorted polychaete annelid (segmented worm) fractions were then sorted again to the identified species (species list available in Supplemental Information or from the authors, on request) (see Supplementary Material available online at http://dx.doi.org/10.1155/2015/983521). The 64 samples were subcores (pseudoreplicates) from 25 box cores and thus the three subcores were aggregated into a single sample from a location. Thus, the species-sample matrix contained 25 samples containing a total of 146 separate polychaete species.

Total wet weight (mg) of each polychaete species was obtained by weighing the sorted samples on a Sartorius CP2P and OHAUS AS260D balance after blotting on a paper towel for 30 seconds. Mean body weight of individuals of each “small” species in each subcore was then calculated by dividing the total weight of the species by the number of individuals in the subcore. Individuals of large species were weighed separately.

A species-sample matrix was constructed and the species abundance data were fourth-root transformed prior to multivariate analyses. CLUSTER was carried out using PRIMER 6.0. SIMPROF (Similarity Profile) was used simultaneously with CLUSTER to test the significance () of differences between clusters [51]. SIMPER estimated the contribution of each species to a faunal group and the contribution of each species to the dissimilarity between faunal groups. The species that contributed the most to a group was considered a “characteristic” species [52]. Regressions were drawn by SigmaPlot 10.0 and Excel 2007. Best-fit regression models were chosen based on the adjusted and numbers of variables. Significant differences between the mean abundances were assessed with a one-way analysis of variance (ANOVA).

Biomass size spectra (BSS) and normalized biomass size spectra (NBSS) were constructed in 14 logarithmic (log2) biomass (mg wet weight) size classes. These size classes are aligned along the -axis as the “independent” variable. NBSS is a graph of the biomass of a size category [] divided by the “width” of that size category ( is the difference between the upper and lower boundaries of the size class) in log2 units on the dependent or -axis as a function of the size categories in log2 units on the independent or -axis. In the figures, we represented as ; that is, [37, 38]. Analysis of covariance (ANCOVA) was used to test the significance of differences of the slope and intercept of the NBSS linear regression lines. All statistical tests were performed with SPSS 17.0.

3. Results

3.1. Oxygen Concentrations

Oxygen in the bottom water varied appreciably over the course of the sampling (Figure 2): 12 were oxic and 12 were hypoxic. In April 2004 and March 2007 all three sites were oxic. On the other hand, in June 2004 and August 2004 all three sites were hypoxic. At the other locations and times the sites were mixed, but not in a predictable fashion. In March 2005 the middle site B was hypoxic but A and C were not; by May 2005 the site closest to the river (A) had dropped to the hypoxic boundary but mid-site B had risen into oxic concentrations. The last three samplings were separated by several years each and exhibited a wide range of conditions: totally oxic in March 2007, three of four locations hypoxic in August 2009, and two of three sites hypoxic in August 2012. A nearshore site, AB5, added in late summer of 2009 and 2012, was always hypoxic; the most distal location, C, was always oxic from 2007 to 2012. The highest single concentration overall was 6.76 mg/L at the most distant site in early spring (C, March 2007), whereas the lowest was 0.12 mg/L in late summer (B and AB5, August 2009). In general, as expected, higher values were located at the greatest distance from the river (C) during spring (March and April), whereas the lowest were in mid- to late summer and close to the river (A, AB5, and B). There were however some persistent exceptions. Likewise, there was not a consistent spatial pattern. (It should be noted that variations in oxygen concentration in time and space have been the subject of numerous publications, especially near sites B and AB5, for several decades [1, 2, 46, 47, 53, 54].) All MCH data are available through NODC under Accession number 0088164.

3.2. Distribution of Mean Animal Sizes

The mean abundance of polychaetes in all 64 subcores was 2746 ind. polychaetes/m2 ( = ±2258), and these amounted to 6.9 gm wet weight/m2 ( = ±5.3). There was a significant relationship between mean body size and bottom oxygen concentration (Figure 3): higher oxygen meant bigger worms. Mean body size ranged from 0.6 mg to 30 mg, with an average value of 3.99 mg per worm ( = ±4.7). The highest mean size (30 mg) was sampled in August 2009 and included the large bamboo worms Clymenella torquata and Asychis elongatus.

Throughout most of the early sampling period (2004-05), the average body size of the polychaete populations sampled was 2.97 mg wet weight ( = ±2.87); however in 2007–09 the total mean body size rose to 7.13 mg ( = ±7.60) (). By 2012, this value dropped back to 2.85 mg ( = ±2.97), but this latter mean was not significantly different () from any of the prior values.

The overall mean biomass was 6.9 gm wet weight/m2 ( = ±5.2 gm/m2; ) for the entire study. If the entire set of samples is arbitrarily divided into three time blocks (2004-05, 2007–09, and 2012), the mean of the total wet weight biomass for all the polychaetes was 7.6 ( = ±5.9), 6.3 ( = ±4.4), and 4.7 ( = ±2.7) g/m2, but these were not statistically different (). Thus, while sizes may have decreased in hypoxic conditions, the total biomass remained relatively high at 6.9 gm wet weight/m2 and did not vary consistently.

The overall mean abundance was 2,747 ind./m2 ( = ±2259; ). During the early intense hypoxia this was 3321 ( = ±2447) inds./m2, whereas abundance went down to  inds./m2 in the sampling two years later (2007 to 2009, ). Seven years later (2012) abundance was at an intermediate level (2226 inds./m2, = ±1049 inds./m2).

3.3. Biomass Size Spectra (BSS, Figure 4)

The biomass size spectra have been plotted at the three different categories of oxygen concentration: oxic, hypoxic, and severely hypoxic (Figure 4). When oxygen was higher than 2 mg/L, large worms were present in size classes 7 (64–128 mg) and 8 (128–256 mg). When oxygen concentration was lower than 2 mg/L, species in size classes 7 and 8 disappeared, but when oxygen dropped down below 1 mg/L, another two size classes (5 and 6) also disappeared. That is, ultimately, there were no big worms when oxygen was severely low. On the other hand, the smallest worms in the tiniest sizes were eliminated too (size classes −1 down to −5), leaving only intermediate sized species (size classes 0 to 4).

The three spectra had some similar peaks within the same size classes (0, 3, and 5, Figure 4); the exception was the large numbers in very small size class (−3) at oxygen lower than 2 mg/L (red line, Figure 4). Species in that size class (−3) were mostly Cossura soyeri, Levinsenia gracilis, Prionospio cristata, and Magelona sp. H. This peak indicates that these relatively small species flourished under moderately hypoxic conditions, but then, as the hypoxia worsened, they declined too (green line).

3.4. Normalized Biomass Size Spectra (NBSS, Figure 5)

The slopes of the normalized biomass size spectra (NBSS) of both oxic and hypoxic (oxygen < 2 mg/L) conditions were significant () (Figure 5). The slope of hypoxic (oxygen < 2 mg/L) NBSS appeared to be slightly more negative than the slope for oxic conditions, but the difference was not significant (−0.68 versus −0.51; ). When oxygen was below 1 mg/L, the NBSS regression fit was very low () with no apparent slope (). This reflects the elimination of both the largest and the smallest species under these severe conditions.

P. pinnata (size class 0) is a medium-sized polychaete and the highest point in the “poor” NBSS plot was in size class 0, dominated by P. pinnata (Table 2). When oxygen concentration was high, small-sized species and P. pinnata were almost equal in abundance. When oxygen was <1 mg/L, the percentage of small-sized and large-sized worms both declined, leaving P. pinnata to flourish, with abundances four to five times more than the small-sized species. This peak in the middle attributed to P. pinnata produced the peculiar nonlinear shape of the NBSS.

3.5. Species Groups Associated with Location, Time, and Hypoxia

Similarity analysis of the entire list of 146 species from all 25 locations was conducted to determine the relationships between species groups, hypoxia, and location along the east to west transect. The 25 samples indicate that there were five faunal groups (>40% similarity; ). The largest group of species (beginning on the left of Figure 6) appeared 18 “times” at sites A, B, and C from 2004 through to 2009 in both spring and summer. This demonstrates that the fauna was remarkably uniform in time and space across most of the entire study. The other four smaller groups were encountered at only seven times and locations (all on the right in Figure 6). Within those seven, the three samples (in red) taken in August 2012 stood out by themselves. The other three groups included only one or two sites taken in August, with the exception of a single sampling in April at site A adjacent to the river’s entrance onto the continental shelf (Figure 6).

The large grouping of species assemblages at all locations and seasons into Group 1 demonstrates that the fauna in general was remarkably uniform; however there is some suggestion that the fauna in Group 1 can be partitioned into subgroups that conform to “zones A, B, and C” as suggested in Rowe and Chapman [5]. The A’s were clustered on the left subgroup (green triangles), whereas the C’s were clustered just to their right, but again as green triangles. Finally the B’s tended to be clustered together next as green triangles on the right of subgroup one (Figure 6). On the left, most of the group was close to the river mouth (A), whereas many in the next subgroup were at C, the most distant location, with B in the middle of the transect. These species subgroups thus tended to reflect location rather than season or oxygen concentration. The distinctiveness of August 2012 at AB5, B, and C indicates that this period was also slightly different in species composition, as was inferred from the size data.

There was little evidence that the groups or subgroups were related to oxygen concentrations at the time of the sampling. For example, Group 2 (red dots) from August 2009 all grouped together but spanned a wide range of oxygen concentrations. The same can be stated for all the groups, on close inspection. Thus, species composition varied very little in time and space, but what variation did occur was not related to the oxygen concentration at the time of the sampling.

P. pinnata has long been recognized as the most common polychaete inhabiting the Gulf of Mexico hypoxic zone [53] and likewise its remarkable abundance dominated Groups 1 (18 branches) and 2 (three red branches) (Figure 6). P. pinnata and Mediomastus californiensis appeared in nearly every sample throughout the study and their abundance accounted for more than half of the total abundance everywhere. They contributed the most to the similarities within Group 1 and between Groups 1 and 2 (Table 2). The characteristic species of Group 1 (18 branches on the left of Figure 6) were all deposit feeders (both selective and nonselective), while Group 2 (branches with red dots) included carnivores, omnivores, and deposit feeders [55].

Nereis micromma contributed most to the dissimilarity between Groups 1 and 2 (3.7%). It was not common in the 2004-05 samples; only two individuals were found during those first five cruises. In 2012 (the red dots of Figure 6) however 98 individuals were sampled (Table 3).

Abundance of P. pinnata increased from  inds./m2 to and  inds./m2 along the decreasing oxygen gradient. Abundance of P. pinnata under hypoxic conditions (<2 mg/L and <1 mg/L) was significantly higher than that in oxic conditions (). Abundance however of M. californiensis was not different within the 3 oxygen levels.

4. Discussion

4.1. Effect of Hypoxia on Mean Animal Size

The results in terms of (a) the mean size of individual polychaetes (mg wet weight), (b) the biomass size spectra (BSS), and (c) the normalized biomass size spectra (NBSS) all demonstrate that animals were smaller under hypoxic conditions. Clearly, however, the northern Gulf of Mexico hypoxic area was not a “dead zone” as is often suggested in the popular media. Animals were small but not gone; in fact some were in great abundance when oxygen was low. The total biomass of the worms did not differ statistically throughout the wide range of conditions, even though mean sizes did vary.

A principal reason for the shift in mean size was exclusion of large-sized species that may be more vulnerable to low oxygen physiological stress [28], not a change in size within species [29]. Smaller species with a greater surface to volume ratio may be better able to cope with low oxygen than large species, with the exception of the medium-sized P. pinnata that may be able to survive because of its large branchial palps.

The positive relationship between size and oxygen agrees with the patterns observed by McClain and Rex [26] in turrid gastropods and “polar gigantism” in amphipod crustaceans [22]. Those environmental conditions however were at the high end of oxygen concentrations, so comparisons are not particularly relevant to physiological stress.

The general negative slope of an NBSS reflects biomass dominance by small size classes. The relative steepness (negativity) of the slope may be a reflection of stress if low oxygen stress favors small sizes. In our study, the NBSS slope ranged from −0.56 to −0.48, but Quiroga et al. [34] reported a much broader range (−0.481 to −0.908) in the oxygen minimum zone (OMZ) off Chile. The OMZ however is a more permanent condition that often approaches anoxia [35]. On the other hand, Saiz-Salinas and Romas [31] tabulated an even more negative NBSS slope along a depth gradient off Antarctica (−0.76 to −1.31), but the latter reflects the commonly observed decrease in mean size of the macrobenthos which accompanies the decrease in available food as water depth increases [17, 19, 20, 56, 57] and high mean sizes at shallow depths at high latitudes [18, 22].

4.2. The Ubiquity of Paraprionospio pinnata and the “Altered” NBSS

The lack of a significant negative regression in the NBSS when oxygen was lower than 1 mg/L was related to the impact on very small species but not on the ubiquitous P. pinnata. Generally, small-sized species tend to be relatively more abundant under stressful conditions [13, 36] and thus NBSS regression lines steepen as the larger animals are eliminated due to low oxygen stress; the smallest-sized polychaetes in size class −5 (Mediomastus californiensis, Aricidea (Acmira) sp., and Aricidea taylori) thus theoretically would be expected to be dominant [13]. But this was not the case. Low oxygen did depress the very small species but not P. pinnata. On the contrary, an increase in abundance of P. pinnata was encountered under moderate and severe hypoxia, reinforcing the findings of Baustian and Rabalais [54] at their station C6B.

Lamont and Gage [29] suggested that the enlargement in size and branching of the branchiae of the polychaete Paraprionospio sp. A on the continental slope off Oman were an adaptation to low oxygen in the oxygen minimum zone of the Arabian Sea. Morphological adaption of P. pinnata to low oxygen might be inferred from the inverse relationship with oxygen concentration noted by Rabalais et al. [53] and Baustian and Rabalais [54]. Our data agree that P. pinnata densities actually increased when subjected to hypoxia. It is reasonable to assume that their large branchial palps provide an enhanced ability to better scavenge oxygen at very low concentrations, but whether this is a phenotypic selection or genetic differentiation as implied by Lamont and Gage [29] remains as yet unknown [58].

The results of this study have implications for the oceanic biota as a whole. While the overall polychaete biomass (total grams, wet weight m−2) remained relatively high, mean size of assemblages declined in a continuum as oxygen decreased. Likewise, the overall ecological impact of low oxygen was to diminish both species richness and the diversity of individual size together at the same time at extreme conditions. While size in marine animals has increased over geologic time scales [59] in those habitats with a surfeit of oxygen, it must be acknowledged that size and diversity may also be decreasing as oxygen declines in many environments worldwide [8].

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This study was supported by the Mechanisms Causing Hypoxia (MCH) project at Texas A&M University with funding from NOAA, Steve DiMarco, Principal Investigator (NOAA Center for Sponsored Coastal Research NA03N0S4780039 and NA06N0 S4780198). The authors are grateful to Dr. Matthew Howard for constructing the map. Fangyuan Qu was supported through the joint graduate education program agreement between the Department of Oceanography at Texas A&M and Ocean University of China.

Supplementary Materials

The SM contains a complete list of polychaete species in each of the 64 sample.

  1. Supplementary Material