Next Article in Journal
An Enhanced Feed-Forward Back Propagation Levenberg–Marquardt Algorithm for Suspended Sediment Yield Modeling
Next Article in Special Issue
Using Single-Species and Whole Community Stream Mesocosm Exposures for Identifying Major Ion Effects in Doses Mimicking Resource Extraction Wastewaters
Previous Article in Journal
Recent Advances in the Reactor Design for Industrial Wastewater Treatment by Electro-Oxidation Process
Previous Article in Special Issue
A Vegetation Assessment of the Kearl Treatment Wetland following Exposure to Oil Sands Process-Affected Water
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evaluation of Nitrogen and Carbon Stable Isotopes in Filter Feeding Bivalves and Surficial Sediment for Assessing Aquatic Condition in Lakes and Estuaries

1
US Environmental Protection Agency, Office of Research and Development, 27 Tarzwell Drive, Narragansett, RI 02882, USA
2
General Dynamics Information Technology, 27 Tarzwell Drive, Narragansett, RI 02882, USA
3
US Environmental Protection Agency, Office of Research and Development, 26 Martin Luther King Drive W., Cincinnati, OH 45268, USA
*
Author to whom correspondence should be addressed.
Water 2022, 14(22), 3712; https://doi.org/10.3390/w14223712
Submission received: 6 October 2022 / Revised: 9 November 2022 / Accepted: 11 November 2022 / Published: 16 November 2022
(This article belongs to the Special Issue Applied Ecology Research for Water Quality Management)

Abstract

:
Excessive inputs of nitrogen from anthropogenic activities in watersheds can cause detrimental effects to aquatic ecosystems, but these effects can be difficult to determine based solely on nitrogen concentrations because of their temporal variability and the need to link human activities to ecological responses. Here, we (1) tested the use of stable isotopes of nitrogen (δ15N) and carbon (δ13C) in benthic organic matter (BOM) as proxies for isotope ratios of filter feeding bivalves in lakes and estuaries, which can be used as indicators but are harder to sample and often spatially sparse, and (2) evaluated if stable isotope ratios in benthic organic matter could be used to assess impacts from anthropogenic land development of watersheds. The δ15N in BOM isolated from surficial sediment (δ15NBOM) was significantly correlated with δ15N in filter feeding unionid mussels (Elliptio complanata, δ15NUN) from lakes and with hard-shell clams (Mercenaria mercenaria, δ15NMM) from estuaries. In lakes, δ13CBOM was significantly correlated with δ13CUN, but δ13CBOM was not significantly correlated with δ13CMM in estuaries. Values of δ15NBOM and δ15NUN were significantly and positively correlated with increasing amounts of impervious surface, urban land cover, and human populations in watersheds surrounding lakes. In estuaries, δ15NBOM was only significantly and positively correlated with greater percent impervious surface in the watersheds. Correlations of δ13CBOM in lakes and estuaries, δ13CUN, and δ13CMM with land use and human population were mostly non-significant or weak. Overall, these results show that δ15NBOM can serve as a proxy for δ15N of filter feeding bivalves in lakes and estuaries and is useful for assessing anthropogenic impacts on aquatic systems and resources. Our study area was limited in size, but our results support further studies to test the application of this sediment stable isotope-based technique for assessing and ranking aquatic resources across broad geographical areas.

1. Introduction

Large monitoring programs have found that nutrient over-enrichment, often resulting from human wastewater and agricultural and urban runoff, is a leading cause of degradation of environmental conditions in aquatic systems across the United States [1,2,3]. The extent and magnitude of nitrogen loading to aquatic ecosystems are of particular concern because they are linked to species loss, significant degradation of water quality, eutrophication, and the development of harmful algal blooms [4,5,6]. These negative effects of nutrient pollution extend from watersheds to coastal and estuarine ecosystems that are commonly nitrogen limited [7,8,9]. The condition of lake and estuarine ecosystems can be difficult to determine based on nitrogen concentrations alone because of their temporal variability and the need to link concentrations to ecological effects [10,11,12,13]. Identifying how organisms are associated with watershed conditions and how they can be used as indicators provides context for understanding and managing nutrient effects, especially because they integrate nutrient conditions over time [14,15,16].
To provide a time integrated metric of nitrogen pollution in aquatic sites, studies have identified relationships between stable isotopes of nitrogen in aquatic organisms and variables indicative of anthropogenic land development in the surrounding watersheds [17,18,19]. Relationships like these provide information on ecological effects of nitrogen, can be used to develop indicators for use in monitoring programs, and can help inform nutrient management efforts. Filter feeding bivalves, such as mussels, are widely used for inferring trophic baselines in streams, rivers, and lakes because of their long lifetimes and because their δ15N ratios change depending on the mix of nitrogen inputs from human wastewater, atmospheric deposition, and agricultural fertilizers entering aquatic sites [20,21,22,23,24]. Recently, the δ15N ratios of chironomids were found to be useful for understanding nitrogen inputs, processing, and transformations in streams and lakes as well, and their δ15N ratios were associated with landscape, water body, and biological factors [25]. In saltwater, δ15N ratios in both soft tissues and shells of hard-shell clams (Mercenaria mercenaria, MM) have been found to directly relate to the percentage of wastewater input to small estuaries [26,27]. Overall, these studies support the use of δ15N in these organisms as a time integrative measure of human impact; however, the use of chironomids, which are almost exclusively freshwater taxa, and the use of filter feeding bivalves as indicators can be limited by the availability of organisms and constrained by the need to compare among species with potentially differing feeding strategies and physiologies [25,28]. These limitations highlight the need for a readily available ecological indicator for assessing aquatic condition in lakes and estuaries.
The δ15N in benthic organic material (BOM) (δ15NBOM) isolated from surficial sediments is one sample type potentially useful as a baseline; especially because of its nearly universal availability in lakes and estuaries. δ15NBOM can be significantly correlated with δ15N of organisms at different trophic positions despite BOM being composed of a mixture of detritus, benthic microbes, and deposited living and dead plankton from the water column. For example, significant relationships were found between δ15NBOM and δ15N of fish and mussels collected from lakes representing a wide range anthropogenic impact [22,28]. Additionally, significant relationships have been reported between δ15N in sedimented organic material from surface sediments of tidal flats and δ15N in manila clams, (Ruditapes philippinarum) among sites spanning a wide range of eutrophication levels [29]. In another study, δ15N ratios in two species of bivalves from estuaries on Cape Cod, USA had significant relationships with δ15N in organic material from sediments [30]. The use of δ15N in sediment as an indicator is further supported by studies that found a correspondence of δ15N ratios between sediment and seston [31] and findings that surface sediment organic matter recorded the weighted mean isotopic ratio produced in surface waters [32]. Carbon stable isotopes (δ13C) provide additional information that complements δ15N results by establishing a baseline value for food chains and for calculation of organismal trophic positions in aquatic systems [21,28,33,34].
With multiple studies showing the potential usefulness of relationships between δ15NBOM and δ15N of other organisms, we surveyed 51 lakes and 28 sites in nine estuaries to document and compare these relationships among sites in these two distinct ecosystems and to examine their possible use as indicators of the effects of anthropogenic disturbance in their watersheds. Our goals were: (1) to determine if δ15NBOM and δ13CBOM sampled from estuaries and lakes represented corresponding isotopic ratios in filter feeding bivalves (primary consumers), and (2) to evaluate the use of δ15N and δ13C in BOM and filter feeding bivalves as indicators of human impact on aquatic systems resulting from anthropogenic land development in watersheds. Documenting these relationships is important to developing the use of δ15NBOM as an indicator of nitrogen effects, which could be particularly useful because of the widespread availability of BOM in lakes and estuaries. Results from this study can be useful to monitoring programs interested in identifying effects of nutrient pollution and for informing decisions and management strategies seeking to protect aquatic resources and to improve those experiencing negative effects of eutrophication from human sources. Results from the present study also may provide support for the use of stable isotopes in BOM for national assessments and rankings of environmental conditions in estuaries [3] and, if appropriate, as an addition to other approaches (e.g., use of chironomids) in future assessments in lakes, streams, and rivers [1,2].

2. Methods

Lakes and estuaries were selected to represent a wide range of anthropogenic land development of surrounding watersheds. Single samples of BOM from surficial sediments were collected from littoral zones of 51 lakes in Rhode Island, USA. These lakes had small surface areas, ranging from 3.4 to 430 ha, but were within watersheds that ranged from 23 ha to 2250 ha for lakes that were part of larger river systems (Supplementary Table S1). Unionid (UN) mussels (Elliptio complanata) were collected from the littoral zones by hand or with a clam rake within ~5 m of where the sediment was sampled for BOM, but they were observed and sampled in only 26 of the 51 lakes. Triplicate samples of BOM from surficial sediments were collected from a total of 28 sampling stations along the shorelines of four Rhode Island and five Massachusetts estuaries, which had watersheds ranging from 250 ha in the smallest salt pond to 275,000 ha for Narragansett Bay. Hard shell clams (Mercenaria mercenaria; MM) were collected with a clam rake within a ~10 m2 area around the location of the BOM sample.
Bivalve tissue samples were from a section of the foot of UN and MM. These individual samples were oven dried at 60° C for >2 days, ground and weighed for isotope ratio mass spectrometry (IRMS) analysis. Each BOM sample and bivalve sample was analyzed individually, and means, or individual δ15NUN and δ13C values where only one bivalve was found, were used in comparisons. We report the number and sizes of bivalve samples found at sample sites in Supplementary Tables S2 and S3.
Sediment samples in both lakes and estuaries were collected with a hand-held piston coring sampler, which obtained a 6.5 cm diameter core. In both lakes and estuaries, samples were taken by wading where water depths were <1 m. The corer was driven into the sediment, which varied from sandy to highly silty. The corer was withdrawn and, while maintained vertically, a core cap was pushed onto the bottom of the core tube, which was then detached from the sampler. Samples of benthic organic material (BOM) were obtained from the surface layer (~the top 1.5 to 2.0 cm) of the intact sediment core. The surface layer was re-suspended by stirring, and the suspended material was removed with a 30 mL pipette. Samples were placed in clean plastic bottles, capped, placed on ice, and transported to the laboratory. Samples were refrigerated (4 °C) until they were processed within 96 h of sampling. Samples were stirred, poured through a 0.5 mm Nytex® screen (to remove large plant material, leaves, twigs, gravel, etc.) and collected in 55 mL centrifuge tubes. Samples were shaken and centrifuged at 1000× g for five minutes. The supernatant was poured off and a spatula was used to take samples from the top ~1.5 cm of the material in the bottom of the tube. Samples were dried in an oven at 40 °C for at least three days, ground with mortar and pestle to a fine powder, weighed to 0.01 mg on a microbalance, and analyzed using an Isotope Ratio Mass Spectrometer (IRMS).
The nitrogen and carbon isotopic composition of the BOM samples and tissue samples were determined by IRMS using a Vario Micro Elemental Analyzer interfaced to an Elementar Isoprime 100 Mass Spectrometer (Elementar Americas, Mt. Laurel, NJ, USA). The nitrogen isotopic composition of the samples is expressed as a part per thousand deviation (δ15N ‰) from that of the reference material, which was N2 in air. Carbon isotopic composition is expressed as a part per thousand deviation (δ13C ‰) from the reference standard Pee Dee Belemnite. All samples were analyzed in duplicate with a typical difference of about 0.1‰ and reported as mean values. An examination of results from wheat flour taken through all steps of the sample preparation procedure with original unprocessed flour showed no appreciable change in δ15N or δ13C ratios demonstrating that no contamination of samples resulted from sample preparation.
We were concerned that collection of BOM samples might also include shell fragments, and if so, that the elevated δ13C of carbonate (CaCO3) in mollusk shells [35] and in sediment might be problematic for using δ13CBOM for comparisons with δ13C in bivalves, and for relating to land cover variables. In lakes, we conducted a pilot study to determine whether removal of inorganic carbon in CaCO3 was necessary prior to analysis of BOM by IRMS. We selected a total of 67 BOM samples from 42 lakes and removed inorganic carbon by acid fumigation [36]. Results showed only minimal changes in δ13CBOM after acid fumigation. The mean difference of δ13CBOM untreated − δ13CBOM fumigated was = −0.05‰, range (−0.23 to + 0.07‰), n = 67); due to this small difference in lakes, we used untreated δ13CBOM samples. In estuaries, δ13CBOM from three stations showed large decreases in δ13CBOM after acid fumigation. Therefore, in estuaries we used results from acid fumigated δ13CBOM samples for comparisons.
In lakes, the means of δ15NUN and δ13CUN or individual values (where only one organism was found) were used for the comparisons with δ15NBOM and δ13CBOM and with land cover variables and population. In estuaries the mean values of the triplicate δ15NBOM and δ13CBOM ratios were compared with the corresponding mean δ15NMM and δ13CMM ratios, or with individual values where only one bivalve was found, at a total of 28 stations from nine estuaries. Only one station was sampled in Narragansett Bay, but the nitrogen and carbon isotope ratios in this station, δ15NBOM = 7.8‰ ± 0.5‰ SD; δ13CBOM = -18.3‰ ± 0.3‰ SD and δ15NMM = 13.0‰ ± 0.2‰ SD; δ13CMM = −17.4‰ ± 0.2‰ SD, agreed with summary values for Narragansett Bay calculated from other sources. These sources measured nitrogen and carbon isotopes in the surface sediment (top 2 cm) of core samples from the axis of the Bay and found δ15N = 7.9‰ ± 0.35‰ SD, and δ13CBOM = −19.9‰ ± 0.84‰ SD [37]. The M. mercenaria samples were from stations taken throughout the Bay and had means of δ15NMM =13.2‰ ± 0.54‰ SD, and δ13CMM −16.76 ± 0.61‰ SD [38]. The grand means of δ15NBOM and δ15NMM from all stations within each of the nine estuaries were used for comparison with land cover and population variables (n = 9).
For comparisons of δ15N ratios in filter feeding bivalves with those in BOM or surface sediment, we developed two relationships: (1) δ15N in mussels (δ15NUN) as a function of δ15NBOM in the 26 lakes, and (2) δ15N in hard-shell clams (δ15NMM) as a function of δ15NBOM from a total of 28 stations from nine estuaries. These regressions were compared using analysis of covariance (ANCOVA). Similar relationships with δ13CUN as a function of δ13CBOM, and δ13CMM as a function of δ13CBOM were developed and analyzed.
Land cover data within the watersheds of lakes and estuaries were obtained from RIGIS (2003–2004) [39] and Mass GIS (2005) [40] imagery with 0.6 m and 0.5 m resolution, respectively. These data were categorized using the land use designations of urban/residential, agricultural, forest, wetland, barren, and water, and percentages of each land group were calculated. Percent impervious surface in watersheds was calculated separately within the watersheds of lakes and estuaries obtained from RIGIS (2007) [39] and Mass GIS (2005) [40]. Population estimates for 2010 were made using dasymetric rasters developed using a USGS tool [41] (Sleeter and Gould 2007, updated 2016) with data from United States Census Bureau [42] (2010). Land use data were from RIGIS (2011) [39] for Rhode Island and from Mass GIS (2005) [40] for Massachusetts.
We compared Spearman and Pearson correlations of land cover variables and population in watersheds with δ15NBOM, δ13CBOM in 51 lakes and in the nine estuaries. Similar comparisons were also made using only data for δ15NBOM, δ13CBOM and with δ15NUN, δ13CUN in the 26 lakes where mussels were found. Comparisons were also made between δ15NBOM and δ15NMM and between δ13CBOM and δ13CMM in the nine estuaries. Correlations used individual δ15NBOM and δ13CBOM values and site mean values of δ15NUN and δ13CUN for lakes and grand means of δ15NBOM, δ13CBOM, δ15NMM, and δ13CMM from all stations within each estuary. We compared the results from both correlation methods to determine the extent to which the presence of non-normality of some variables in the Pearson correlations impacted results.
We compared relationships of land cover variables and population with each other and selected log % impervious surface as the most appropriate variable for representing anthropogenic land development. We developed regression relationships with log % impervious surface as a function of δ15NBOM from the lakes and estuaries.
Correlations, least squares regressions, ANCOVA, and regression model selection were done using SAS software version 9.4 (SAS Institute Cary, NC, USA). Student’s t-tests, F-tests for equality of variance, and regression analyses were done using Microsoft Excel. We reported results of Pearson correlation coefficients with significance indicated for pairwise comparisons to show relationships within the data. The level of significance for statistical tests was p ≤ 0.05.

3. Results

3.1. Comparisons of Bivalves with BOM

Increased δ15N in mussels (δ15NUN) from lakes and increased δ15N in hard shell clams (δ15NMM) from estuaries were associated with greater δ15NBOM (R2 = 0.82 and p < 0.001 for both; Figure 1). The slopes of these relationships were similar in lakes and estuaries (ANCOVA, p = 0.12), but their intercepts were significantly different (p < 0.001). Increased δ13CUN were associated with greater δ13CBOM in lakes (R2 = 0.70, p < 0.001), but δ13CMM was not significantly correlated with δ13CBOM in estuaries (Figure 2).

3.2. Comparisons of BOM with Land Cover and Population

We found a general agreement of the Pearson and Spearman correlations of δ15NBOM with land cover variables and population in watersheds surrounding lakes and estuaries which indicated that the lack of normality in some variables in the Pearson correlations did not adversely impact results (Table 1). In the watersheds of 51 lakes, we found δ15NBOM was significantly positively correlated with % urban residential land, log % impervious surface and log population and significantly negatively correlated with % forest. Lower correlation coefficients of δ15NBOM were found with % wetland, % water, and % agriculture and these correlations were not significant (Table 1). The correlation of δ15NBOM with % barren land, which represented both % barren land and % beaches, showed a significant, but low correlation. A review of the % barren land variable showed 28 of the 51 lakes had % barren land = 0; therefore, the significant correlation was likely an artifact.
In the watersheds of nine estuaries, correlations were not as strong as in the lakes, and significant Pearson correlations were found only for δ15NBOM with log % impervious surface. The variable % forest in watersheds of estuaries was not significant as it was in the lakes, but this may have resulted because of the large extent of salt marsh grasses surrounding several of the estuarine sites and the consequent decreases in % forest. The variables % wetland and % water showed significant correlations with δ15NBOM in estuaries, but this resulted because of the large extent of wetland and ponds surrounding one undeveloped site, Sage Lot Pond, which had the lowest δ15NBOM value and the highest percentages of wetland and water. When the percentages of wetland and water for Sage Lot Pond were removed from the correlations, as in Table 1, their correlations with δ15NBOM became non-significant.
In the 26 lakes with mussels, significant positive Pearson correlations were found for both δ15NBOM and δ15NUN with log % impervious surface, % urban residential and log population, and significant negative correlations were found with % forest (Table 2). Similar correlations of δ15NBOM and δ15NMM with land cover and population in watersheds of the estuaries showed a significant correlation only for δ15NBOM with log % impervious surface (Table 2).
Our comparisons of regressions of land cover and population variables with each other and with δ15NBOM showed the highest coefficients of determination with log % impervious surface, which we used for representing anthropogenic land development and compared it with δ15NBOM in the lakes and estuaries. In lakes we found a significant regression of log % impervious surface as a function of δ15NBOM in the watersheds (Figure 2). An empirical equation was derived from this regression (Equation (1)).
log % impervious surface = 0.16 (δ15N BOM) + 0.16, R2 = 0.62, p < 0.001, n = 51
Within the estuaries we found a significant regression of log % impervious surface as a function of δ15NBOM in the watersheds (Figure 3 and Equation (2)). To determine how the removal of the much larger Narragansett Bay from the regression would impact the relationship we omitted it from a second regression analysis. The regression without the Narragansett Bay watershed was only slightly different (Equation (3)).
log % impervious surface = 0.05 (δ15N BOM) + 0.82, R2 =0.57, p = 0.019, n = 9
log % impervious surface = 0.056 (δ15N BOM) + 0.79, R2 =0.60, p = 0.02, n = 8
Comparisons of δ13CBOM with land cover and population variables showed no significant correlations for watersheds of lakes, or estuaries.

4. Discussion

For environmental assessment of how anthropogenic activities in watersheds affect lakes and estuaries, our results suggest that a sediment-based method presents a useful alternative to mussels and other organisms, which commonly can be absent from sites. In the present study, we found a statistical similarity of slopes and R2 values (0.82) of regression relationships between δ15NUN as a function of δ15NBOM in lakes and δ15NMM as a function of δ15NBOM in estuaries. This level of agreement indicates that δ15NBOM from sediment closely duplicates δ15N found in bivalves, which are often considered primary consumers and have been used in environmental assessments [22,23,24,43]. We found different intercepts of regressions between the lakes and the estuaries. The different intercepts may indicate variation in accumulation through feeding or metabolism by the different species of bivalves present in lakes versus those of estuaries. Studies indicate the dietary components of mussels are poorly understood [44,45], and species differences in δ15N of mussels taken from the same sites have been reported [24,46]. Additionally, species-specific differences of δ15N in hard shell clams and in soft shell clams (Mya Arenaria) collected from the same small estuaries on Cape Cod, USA have been found [26].
In lakes of the present study, the relationship of δ13CUN as a function of δ13CBOM was significant; however, the relationship of δ13CMM from estuaries as a function of δ13CBOM was not significant. The diet of M. Mercenaria is planktonic microalgae [47], whereas the diets of mussels are more varied as cited above [45,46]. These dietary differences may be the reason for the relatively small range of δ13CMM relative to δ13CUN. Our findings indicate that although δ15NBOM could be used as a surrogate for bivalves in assessments of nitrogen concentrations and impacts in lakes and estuaries, δ13CBOM appears to have utility as a surrogate for δ13C only for bivalves in lakes.
Our comparisons of correlations of δ15NBOM with variables indicative of anthropogenic land development of watersheds surrounding lakes showed significant positive correlations with log population, log % impervious surface, % urban residential, and a negative correlation with % forest. These results support the use of δ15NBOM for assessments in watersheds surrounding lakes. Due to the dynamic nature of estuaries, it is often difficult to assess their environmental condition.
Environmental assessments in estuaries have used a variety of methods to develop and refine benthic organism-based indices to evaluate ecological status of sites [48,49,50]. However, collection and identification of species present in estuarine sediment samples are labor intensive and subject to variability depending upon seasonality and site-specific variables (e.g., sediment type). Results of analyses of nutrients and other aquatic variables also have been used to assess and rank aquatic site condition, as in National coastal condition assessments [3]. These assessments have been highly useful for summarizing data, but they do not link nutrients directly to ecological responses and nutrient concentrations are inherently variable over time. Thus, the addition of a time integrative measure such as δ15NBOM may be useful for these assessments.
The correlations of δ15NBOM with land use and population variables in estuarine watersheds were not as strong as in lakes, and δ15NBOM was only significantly correlated with log % impervious surface. These results likely reflect the highly dynamic nature of the estuarine environment and the influence of freshwater and seawater inputs, tides, complex flows, and variable marine algal (food) sources on the composition of the organic material in sediments and therefore on δ15NBOM in estuarine sediments, but still show that δ15NBOM in estuaries is associated with watershed development. One of our estuaries, Narragansett Bay, has a watershed which is ~38 times larger than the next largest estuarine watershed, and has been found to be highly influenced by nutrients from wastewater [38]. To determine whether the inclusion of Narragansett Bay watershed had an over whelming influence on the regression, we removed it from the data set and recalculated the regression. We found almost no change in the regression equation which indicates that the regression is relatively robust and provides further support for using δ15NBOM to indicate effects of anthropogenic land development in watersheds of estuaries.
In other studies, the sources of nitrogen inputs to aquatic sites have been found to consist of human or livestock wastes, inputs of human wastewater, or from application of nitrogen in fertilizers in surrounding watersheds. In the present study, the percentage of agricultural land use in watersheds surrounding lakes and estuaries was low and did not represent a significant variable for predicting δ15N in bivalves or δ15NBOM. Our results showed that variables indicating anthropogenic land development (i.e., increased human population, urban/residential development, percentage of impervious surface and decreased forest) were principal explanatory variables and likely indicated the consequent input of wastewater to sites. One study suggested that drivers of δ15N through nitrogen inputs in lakes may show regional differences [51]. Our study sites were all in New England, and our results coincided with anthropogenic development in watersheds. We could not test regional differences but acknowledge that regional differences in nitrogen sources and inputs to aquatic sites are likely [25,52].
The lakes and estuaries of the present study, except for Narragansett Bay, are relatively small; therefore, land use in the surrounding watersheds was more likely to have been reflected in the δ15NBOM. However, in larger waterbodies comparisons of δ15NBOM with anthropogenic activities in watersheds may be influenced by a large range of activities and inputs which are difficult to quantify. Therefore, in large lakes and estuaries, studies to determine the appropriate sampling strategies for using stable isotopes in sediment will be needed for development of adequate assessments.
Overall, the present study showed highly significant positive relationships between δ15N in bivalves and BOM, which indicated that δ15NBOM may be used as a surrogate for δ15N in bivalve mollusks. Our findings suggest that δ15NBOM provides an effective summary metric for assessing and ranking the environmental conditions of lakes and estuaries resulting from anthropogenic land development of watersheds.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/w14223712/s1, Table S1: Location data for study lakes and watershed size (ha) for lakes and estuaries; Table S2: Number, size data, δ15N and δC13 for Mussels (UN) from 26 Lakes; Table S3: Number, size data, δ15N and δC13 and sample site locations for Mercenaria mercenaria from a total of Estuaries 28 stations from nine Estuaries.

Author Contributions

Conceptualization, J.L.L., J.R.S., A.K. and N.J.S.; Data curation, J.R.S.; Formal analysis, J.L.L.; Investigation, J.R.S.; Methodology, J.L.L. and J.R.S.; Project administration, J.L.L.; Supervision, J.L.L. and A.K.; Visualization, J.R.S., A.K. and M.C.; Writing—original draft, J.L.L.; Writing—review and editing, J.R.S., A.K., M.C. and N.J.S. All authors have read and agreed to the published version of the manuscript.

Funding

Michael Charpentier was funded through a contract with General Dynamics Information Technology. The other authors were funded by the U.S. Environmental Protection Agency.

Institutional Review Board Statement

This manuscript has been through ACESD review procedures and has been cleared for publication.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is located in: https://sciencehub.epa.gov/sciencehub/datasets/3803 (accessed 10 November 2022).

Acknowledgments

The authors thank Cathleen Wigand, Michaela Cashman, Laura Coiro, Betty Kreakie and Marty Chintala for their expertise and helpful technical reviews of this manuscript.

Conflicts of Interest

The authors declare no conflict of interest. This manuscript has been subjected to Agency review and has been approved for publication. The views expressed in this article are those of the authors and do not necessarily represent the views or policies of the U.S. Environmental Protection Agency. Any mention of trade names, products, or services does not imply an endorsement by the U.S. Government or the U.S. Environmental Protection Agency. The EPA does not endorse any commercial products, services, or enterprises.

References

  1. United States Environmental Protection Agency. National Lakes Assessment 2012: A Collaborative Survey of Lakes in the United States; Technical Report EPA 841-R-16-113; Office of Water, Office of Research and Development: Washington, DC, USA, 2016.
  2. United States Environmental Protection Agency. National Rivers and Streams Assessment 2008–2009: A Collaborative Survey; Technical Report EPA/841/R-16/007; Office of Water, Office of Research and Development: Washington, DC, USA, 2016.
  3. United States Environmental Protection Agency. National Coastal Condition Assessment 2015: A Collaborative Survey of the Nation’s Estuaries and Great Lakes Nearshore Waters; Technical Report EPA 841-R-21-001; Office of Water, Office of Research and Development: Washington, DC, USA, 2021.
  4. Sondergaard, M.; Jeppesen, E. Anthropogenic impacts on lake and stream ecosystems, and approaches to restoration. J. Appl. Ecol. 2007, 44, 1089–1094. [Google Scholar] [CrossRef]
  5. Wurtsbaugh, W.A.; Paerl, H.W.; Dodds, W.K. Nutrients, eutrophication and harmful algal blooms along the freshwater to marine continuum. WIREs Water 2019, 6, e1373. [Google Scholar] [CrossRef]
  6. Smucker, N.J.; Beaulieu, J.J.; Nietch, C.T.; Young, J.L. Increasingly severe cyanobacterial blooms and deep-water hypoxia coincide with warming water temperatures in reservoirs. Glob. Chang. Biol. 2021, 27, 2507–2519. [Google Scholar] [CrossRef] [PubMed]
  7. Howarth, R.W.; Sharpley, A.; Walker, D. Sources of nutrient pollution to coastal waters in the United States: Implications for achieving coastal water quality goals. Estuaries 2002, 25, 656–676. [Google Scholar] [CrossRef]
  8. Seitzinger, S.P.; Harrison, J.A.; Dumont, E.; Beusen, A.H.W.; Bouwman, A.F. Sources and delivery of carbon, nitrogen, and phosphorus to the coastal zone: An overview of Global Nutrient Export from Watersheds (NEWS) models and their application. Glob. Biogeochem. Cycles 2005, 19, GB4S01. [Google Scholar] [CrossRef] [Green Version]
  9. Boyer, E.W.; Howarth, R.W.; Galloway, J.N.; Dentener, F.J.; Green, P.A.; Vörösmarty, C.J. Riverine nitrogen export from the continents to the coasts. Glob. Biogeochem. Cycles 2006, 20, GB1S91. [Google Scholar] [CrossRef] [Green Version]
  10. Groffman, P.M.; Law, N.L.; Belt, K.T.; Band, L.E.; Fisher, G.T. Nitrogen fluxes and retention in urban watershed ecosystems. Ecosystems 2020, 7, 393–403. [Google Scholar] [CrossRef]
  11. Alexander, R.B.; Smith, R.A.; Schwarz, G.E.; Boyer, E.W.; Nolan, J.V.; Brakebill, J.W. Differences in phosphorus and nitrogen delivery to the Gulf of Mexico from the Mississippi River Basin. Environ. Sci. Technol. 2008, 42, 822–830. [Google Scholar] [CrossRef] [Green Version]
  12. Kaushal, S.S.; Mayer, P.M.; Vidon, P.G.; Smith, R.M.; Pennino, M.J.; Newcomer, T.A.; Duan, S.; Welty, C.; Belt, K.T. Land use and climate variability amplify carbon, nutrient, and contaminant pulses: A review with management implications. J. Am. Water Resour. Assoc. 2014, 50, 585–614. [Google Scholar] [CrossRef] [Green Version]
  13. Loecke, T.D.; Burgin, A.J.; Riveros-Iregui, D.A.; Ward, A.S.; Thomas, S.A.; Davis, C.A.; St Clair, M.A. Weather whiplash in agricultural regions drives deterioration of water quality. Biogeochemistry 2017, 133, 7–15. [Google Scholar] [CrossRef]
  14. Balint, S.J.; Joyce, E.; Pennino, S.; Oczkowski, A.; McKinney, R.; Hastings, M.G. Identifying sources and impacts of precipitation-derived nitrogen in Narragansett Bay, RI. Estuaries Coasts 2021, 45, 1287–1304. [Google Scholar] [CrossRef]
  15. Smucker, N.J.; Pilgrim, E.M.; Wu, H.; Nietch, C.T.; Darling, J.A.; Molina, M.; Johnson, B.R.; Yuan, L.L. Characterizing temporal variability in streams supports nutrient indicator development using diatom and bacterial DNA metabarcoding. Sci. Total Environ. 2022, 831, 154960. [Google Scholar] [CrossRef] [PubMed]
  16. Yuan, L.L.; Smucker, N.J.; Nietch, C.T.; Pilgrim, E.M. Quantifying spatial and temporal relationships between diatoms and nutrients in streams strengthens evidence of nutrient effects from monitoring data. Freshw. Sci. 2022, 41, 100–112. [Google Scholar] [CrossRef] [PubMed]
  17. Karube, Z.; Sakai, Y.; Takeyama, T.; Okuda, N.; Kohzu, A.; Yoshimizu, C.; Nagata, T.; Tayasu, I. Carbon and nitrogen stable isotope ratios of macroinvertebrates in the littoral zone of Lake Biwa as indicators of anthropogenic activities in the watershed. Ecol. Res. 2010, 25, 847–855. [Google Scholar]
  18. Smucker, N.J.; Kuhn, A.; Cruz-Quinones, C.J.; Serbst, J.R.; Lake, J.L. Stable isotopes of algae and macroinvertebrates in streams respond to watershed urbanization, inform management goals, and indicate food web relationships. Ecol. Indic. 2018, 90, 295–304. [Google Scholar] [CrossRef]
  19. Dalu, T.; Cuthbert, R.N.; Taylor, J.C.; Magoro, M.L.; Weyl, O.L.F.; Froneman, P.W.; Wasserman, R.J. Benthic diatom-based indices and isotopic biomonitoring of nitrogen pollution in a warm temperate Australian river system. Sci. Tot. Environ. 2020, 748, 142452. [Google Scholar] [CrossRef]
  20. Cabana, G.; Rasmussen, J.B. Comparison of aquatic food chains using nitrogen isotopes. Proc. Natl. Acad. Sci. USA 1996, 93, 10844–10847. [Google Scholar] [CrossRef] [Green Version]
  21. Vander Zanden, M.J.; Rasmussen, J.B. Primary consumer δ13C and δ15N and the trophic position of aquatic consumers. Ecology 1999, 80, 1395–1404. [Google Scholar] [CrossRef]
  22. Lake, J.L.; McKinney, R.A.; Osterman, F.A.; Pruell, R.J.; Kiddon, J.; Ryba, S.A.; Libby, A.D. Stable nitrogen isotopes as indicators of anthropogenic activities in small freshwater systems. Can. J. Fish Aquat. Sci. 2001, 58, 870–878. [Google Scholar] [CrossRef]
  23. McKinney, R.A.; Lake, J.L.; Charpentier, M.A.; Ryba, S. Using mussel isotope ratios to assess anthropogenic nitrogen inputs to freshwater ecosystems. Environ. Monit. Assess. 2002, 74, 167–192. [Google Scholar] [CrossRef]
  24. Atkinson, C.L.; Christian, A.D.; Spooner, D.E.; Vaughn, C.C. Long-lived organisms provide an integrative footprint of agricultural land use. Ecol. Appl. 2014, 24, 375–384. [Google Scholar] [CrossRef] [PubMed]
  25. Brooks, J.R.; Compton, J.E.; Lin, J.; Herlihy, A.; Nahlik, A.M.; Rugh, W.; Weber, M. d15N of Chironomidae: An index of nitrogen sources and processing within watersheds for national aquatic monitoring programs. Sci. Total Environ. 2022, 20, 151867. [Google Scholar] [CrossRef] [PubMed]
  26. Carmichael, R.H.; Annett, B.; Valiela, I. Nitrogen loading to pleasant bay, cape cod:application of models and stable isotopes to detect incipient nutrient enrichment of estuaries. Mar. Pollut. Bull. 2004, 48, 137–143. [Google Scholar] [CrossRef]
  27. Carmichael, R.H.; Hattenrath, T.; Ivan Valiela, I.; Michener, R.H. Nitrogen stable isotopes in the shell of Mercenaria mercenaria trace wastewater inputs from watersheds to estuarine ecosystems. Aquat. Biol. 2008, 4, 99–111. [Google Scholar] [CrossRef]
  28. Lake, J.L.; Serbst, J.R.; Kuhn, A.; Smucker, N.J.; Edwards, P.; Libby, A.; Charpentier, M.; Miller, K. Use of stable isotopes in benthic organic material as a baseline for estimating fish trophic positions in lakes. Can. J. Fish Aquat. Sci. 2019, 76, 1227–1237. [Google Scholar] [CrossRef] [Green Version]
  29. Watanabe, S.; Kodama, M.; Fukuda, M. Nitrogen stable isotope ratio in the manila clam, Ruditapes philippinarum, reflects eutrophication levels in tidal flats. Mar. Pollut. Bull. 2009, 58, 1447–1453. [Google Scholar] [CrossRef]
  30. Carmichael, R.H.; Rutecki, D.; Annett, B.; Gaines, E.; Valiela, I. Position of horseshoe crabs in estuarine food webs: N and C stable isotopic study of foraging ranges and diet composition. J. Exp. Mar. Biol. Ecol. 2004, 299, 231–253. [Google Scholar] [CrossRef]
  31. Mitchell, J.S.; Bailey, R.C.; Knapton, R.W. Abundance of Dreissena polymorpha and Dreissena bugensis in a warmwater plume: Effects of depth and temperature. Can. J. Fish. Aquat. Sci. 1996, 53, 1705–1712. [Google Scholar] [CrossRef]
  32. Hodell, D.A.; Schelske, C.L. Production, sedimentation, and isotopic composition of organic matter in Lake Ontario. Limnol. Oceanogr 1998, 43, 200–214. [Google Scholar] [CrossRef]
  33. Post, D.M. Using stable isotopes to estimate trophic position: Models, methods, and assumptions. Ecology 2002, 83, 703–718. [Google Scholar] [CrossRef]
  34. Vander Zanden, M.J.; Shuter, B.J.; Lester, N.; Rasmussen, J.B. Patterns of food chain length in lakes: A stable isotope study. Am. Nat. 1999, 154, 406–416. [Google Scholar] [CrossRef] [PubMed]
  35. McConnaughey, T.A.; Gillikin, D.P. Carbon isotopes in mollusk shell carbonates. Geo. Mar. Lett. 2008, 28, 287–299. [Google Scholar] [CrossRef]
  36. Komada, T.; Anderson, M.R.; Dorfmeier, C.L. Carbonate removal from coastal sediments for the determination of organic carbon and its isotopic signatures, δ13C and δ14C: Comparison of fumigation and direct acidification by hydrochloric acid. Limnol. Oceanogr. Methods. 2008, 6, 254–262. [Google Scholar] [CrossRef] [Green Version]
  37. Boothman, W.; (U.S. Environmental Protection Agency Office of Research and Development, 27 Tarzwell Drive, Narragansett RI, USA. Personal communication, 2021. Manuscript in preparation.
  38. Oczkowski, A.; Nixon, S.; Henry, K.; DiMilla, P.; Pilson, M.; Granger, S.; Buckley, B.; Thornber, C.; McKinney, R.; Chaves, J. Distribution and Trophic Importance of Anthropogenic Nitrogen in Narragansett Bay: An Assessment Using Stable Isotopes, Estuaries and Coasts. J. CERF 2008, 31, 53–69. [Google Scholar] [CrossRef]
  39. RIGIS (2011) Rhode Island Geographic Information System. Land Cover/Land Use for Rhode Island. 1988, 2003–2004, 2011. Accessed on 15 July 2016 for Watershed Population and Lake Data; Accessed on 3 December 2019 for Estuaries. 2011. Available online: www.rigis.org (accessed on 6 October 2022).
  40. Mass GIS (2005) Massachusetts Geographic Information System (Mass-GIS). Massachusetts Land Use 1985, 1999, and 2005. (2005) Massachusetts Geographic Information System (Mass-GIS). Massachusetts Land Use Support/Application-Serv/Office-of-Geographic-Information-Massgis. Accessed on 15 July 2016 for Watershed Population and Lake Data; Accessed on 3 December 2019 for Estuaries. Available online: www.mass.gov/anf/research-and-tech/it-serv-and-MassGIS (accessed on 6 October 2022).
  41. Sleeter, R.; Gould, M. Geographic Information System Software to Remodel Population Data Using Dasymetric Mapping Methods: U.S. Geological Survey Techniques and Methods 11-C2; US Department of the Interior, US Geological Survey: Denver, CO, USA, 2007; 15p.
  42. United States Census Bureau. Historical Data: 2010. 2010. Available online: www.census.gov/programs-surveys/popest.html (accessed on 3 December 2019).
  43. McClelland, J.W.; Valiela, I.; Michener, R.H. Nitrogen-stable isotope signatures in estuarine food webs: A record of increasing urbanization in coastal watersheds. Limnol. Oceanogr. 1997, 42, 930–937. [Google Scholar] [CrossRef]
  44. Vuorio, K.; Tarvainen, M.; Sarvala, J. Unionid mussels as stable isotope baseline indicators for long-lived secondary consumers in pelagic food web comparisons. Fundam. Appl. Limnol. 2007, 169, 237–245. [Google Scholar] [CrossRef]
  45. Yasuno, N.; Shindo, K.; Takagi, Y.; Kanaya, G.; Shikano, S.; Fujimoto, Y.; Shimada, T.; Kikuchi, E. Ontogenetic changes in the trophic position of a freshwater Unionidae mussel. Fundam. Appl. Limnol. 2014, 184, 341–349. [Google Scholar] [CrossRef]
  46. Novais, A.; Dias, E.; Sousa, R. Inter- and intraspecific variation of carbon and nitrogen stable isotope ratios in freshwater bivalves. Hydrobiologia 2016, 765, 149–158. [Google Scholar] [CrossRef]
  47. Roegner, G.C.; Roger, L.M. Hard Clam Mercenaria mercenaria. In VIMS Books and Book Chapters 20; Virginia Institute of Marine Science: Gloucester Point, VA, USA, 1991; Available online: https://scholarworks.wm.edu/vimsbooks/20 (accessed on 8 November 2022).
  48. Pelletier, M.C.; Gillett, D.J.; Hamilton, A.; Grayson, T.; Hansen, V.; Leppo, E.W.Q.; Wiesberg, S.B.; Borja, A. Adaptation and application of multivariate (M-AMBI) in US coastal waters. Ecol. Indic. 2018, 89, 818–827. [Google Scholar] [CrossRef]
  49. Pinto, R.; Patricio, J.; Baeta, A.; Neto, B.D.F.J.M.; Marques, J.C. Review and evaluation of estuarine biotic indices to assess benthic condition. Ecol. Indic. 2009, 9, 1–25. [Google Scholar] [CrossRef] [Green Version]
  50. Pruden, M.J.; Dietl, G.P.; Handley, J.C.; Smith, J.A. Using molluscs to assess ecological quality status of soft-bottom habitats along the Atlantic coastline of the United States. Ecol. Indic. 2021, 129, 107910. [Google Scholar] [CrossRef]
  51. Botrel, M.; Gregory-Eaves, I.; Maranger, R. Defining drivers of nitrogen stable isotopes (δ15N) of surface sediments in temperate lakes. J. Paleolimnol. 2014, 52, 419–433. [Google Scholar] [CrossRef]
  52. Sabo, R.D.; Clark, C.M.; Bash, J.; Sobota, D.; Cooter, E.; Dobrowolski, J.P.; Compton, J.E. Decadal shift in nitrogen inputs and fluxes across the contiguous United States: 2002–2012. J. Geophys. Res. Biogeosciences 2019, 124, 3104–3124. [Google Scholar] [CrossRef]
Figure 1. Regression comparisons of δ15N in unionid (UN) mussels (δ15NUN) from 26 lakes (▲), and in hard-shell clams, M. mercenaria (MM) (δ15NMM) from a total of 28 stations from nine estuaries (●), as functions of δ15N in sediment benthic organic material (BOM), (δ15NBOM). δ15NUN = 1.15 δ15NBOM + 2.24, R2 = 0.82, n = 26, p < 0.001. (2) δ15NMM = 0.92 δ15N + 4.86, R2 = 0.82, n = 28, p < 0.001.
Figure 1. Regression comparisons of δ15N in unionid (UN) mussels (δ15NUN) from 26 lakes (▲), and in hard-shell clams, M. mercenaria (MM) (δ15NMM) from a total of 28 stations from nine estuaries (●), as functions of δ15N in sediment benthic organic material (BOM), (δ15NBOM). δ15NUN = 1.15 δ15NBOM + 2.24, R2 = 0.82, n = 26, p < 0.001. (2) δ15NMM = 0.92 δ15N + 4.86, R2 = 0.82, n = 28, p < 0.001.
Water 14 03712 g001
Figure 2. Regression comparisons of δ13C in unionid (UN) mussels (δ13CUN) from 26 lakes (▲), and in hard-shell clams, M. mercenaria (MM) (δ13CMM) from a total of 28 stations from nine estuaries (●), as functions of δ13C in sediment benthic organic material (BOM), (δ13CBOM). δ13CUN = 1.10 δ13CBOM + 0.74, R2 = 0.70, n = 26, p < 0.001. The δ13CMM regression with δ13CBOM in estuaries was not significant.
Figure 2. Regression comparisons of δ13C in unionid (UN) mussels (δ13CUN) from 26 lakes (▲), and in hard-shell clams, M. mercenaria (MM) (δ13CMM) from a total of 28 stations from nine estuaries (●), as functions of δ13C in sediment benthic organic material (BOM), (δ13CBOM). δ13CUN = 1.10 δ13CBOM + 0.74, R2 = 0.70, n = 26, p < 0.001. The δ13CMM regression with δ13CBOM in estuaries was not significant.
Water 14 03712 g002
Figure 3. Regression comparisons of log % impervious surface in watersheds of 51 lakes (dotted line) and nine estuaries (dashed line) as functions of δ15N in sediment benthic organic material (BOM), δ15NBOM. Equations for lines are log % impervious surface = 0.16 (δ15NBOM) + 0.04, R2 = 0.62, p < 0.001, in 51 lakes (●), and log % impervious surface = 0.05 (δ15NBOM) = 0.82, R2 = 0.57, p = 0.02, in nine estuaries (▲).
Figure 3. Regression comparisons of log % impervious surface in watersheds of 51 lakes (dotted line) and nine estuaries (dashed line) as functions of δ15N in sediment benthic organic material (BOM), δ15NBOM. Equations for lines are log % impervious surface = 0.16 (δ15NBOM) + 0.04, R2 = 0.62, p < 0.001, in 51 lakes (●), and log % impervious surface = 0.05 (δ15NBOM) = 0.82, R2 = 0.57, p = 0.02, in nine estuaries (▲).
Water 14 03712 g003
Table 1. Comparisons of Spearman and Pearson correlations of δ15N in sediment benthic organic material (BOM), δ15NBOM with land cover and population variables in watersheds of 51 lakes and nine estuaries. Land cover variables included log % impervious surface (Log%Impsur) which was calculated separately, % urban residential (% UrRes), log population (Logpop), % Forest, % Agriculture (% Agr), % Wetland (% WetL), % Barren, and % Water.
Table 1. Comparisons of Spearman and Pearson correlations of δ15N in sediment benthic organic material (BOM), δ15NBOM with land cover and population variables in watersheds of 51 lakes and nine estuaries. Land cover variables included log % impervious surface (Log%Impsur) which was calculated separately, % urban residential (% UrRes), log population (Logpop), % Forest, % Agriculture (% Agr), % Wetland (% WetL), % Barren, and % Water.
Lakes
Spearman Correlations
Log%Impsur% UrResLogpop% Forest% Agr% WetL% Barren% Water
δ15NBOM0.810.780.54−0.710.15−0.260.31−0.18
p<0.001<0.001<0.001<0.001NSNS0.03NS
Pearson Correlations
Log%Impsur% UrResLogpop% Forest% Agr% WetL% Barren% Water
δ15NBOM0.790.690.59−0.710.14−0.230.35−0.27
p<0.001<0.001<0.001<0.001NSNS0.01NS
Estuaries
Spearman Correlations
Log%Impsur%UrResLogpop% Forest% Agr% WetL% Barren% Water
δ15NBOM0.730.50.330.020.25−0.17−0.570.05
p0.02NSNSNSNSNSNSNS
Pearson Correlations
Log%Impsur% UrResLogpop% Forest% Agr% WetL% Barren% Water
δ15NBOM0.750.480.670.040.31−0.03−0.13−0.07
p0.02NS0.05NSNSNSNSNS
NS indicates pair- wise correlation was not significant.
Table 2. Pairwise Pearson correlation coefficients for δ15N in sediment benthic organic material (δ15N BOM) and δ15N in the bivalves unionid mussels (δ15NUN) and M. mercenaria15NMM) with log percent impervious surface (Log%Impsur), which was calculated separately, % urban residential land (% UrRes), log population (Logpop) and % Forest in watersheds of 26 lakes and nine estuaries.
Table 2. Pairwise Pearson correlation coefficients for δ15N in sediment benthic organic material (δ15N BOM) and δ15N in the bivalves unionid mussels (δ15NUN) and M. mercenaria15NMM) with log percent impervious surface (Log%Impsur), which was calculated separately, % urban residential land (% UrRes), log population (Logpop) and % Forest in watersheds of 26 lakes and nine estuaries.
Lakes
Log%Impsur% UrResLogpop% Forest
δ15NBOM0.790.740.49−0.73
p<0.001<0.0010.011<0.001
δ15NUN0.750.660.63−0.67
p<0.001<0.001<0.001<0.001
Estuaries
Log%Impsur% UrResLogpop% Forest
δ15NBOM0.750.480.670.04
p0.019NSNSNS
δ15NMM0.610.360.550.19
pNSNSNSNS
NS indicates pair- wise correlation was not significant.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lake, J.L.; Serbst, J.R.; Kuhn, A.; Charpentier, M.; Smucker, N.J. Evaluation of Nitrogen and Carbon Stable Isotopes in Filter Feeding Bivalves and Surficial Sediment for Assessing Aquatic Condition in Lakes and Estuaries. Water 2022, 14, 3712. https://doi.org/10.3390/w14223712

AMA Style

Lake JL, Serbst JR, Kuhn A, Charpentier M, Smucker NJ. Evaluation of Nitrogen and Carbon Stable Isotopes in Filter Feeding Bivalves and Surficial Sediment for Assessing Aquatic Condition in Lakes and Estuaries. Water. 2022; 14(22):3712. https://doi.org/10.3390/w14223712

Chicago/Turabian Style

Lake, James L., Jonathan R. Serbst, Anne Kuhn, Michael Charpentier, and Nathan J. Smucker. 2022. "Evaluation of Nitrogen and Carbon Stable Isotopes in Filter Feeding Bivalves and Surficial Sediment for Assessing Aquatic Condition in Lakes and Estuaries" Water 14, no. 22: 3712. https://doi.org/10.3390/w14223712

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop