Figure 5. 3-hour rainfall (top) and ENT concentrations (bottom) from Taylor’s Creek during perigean tide events from June 12-17, 2022 (a) and July 11-16, 2022 (b). The red line represents the EPA’s single sample maximum threshold concentration for safe public use (104 MPN 100 mL-1), and the grayed area represents the minimum detection limit (below 10 MPN 100 mL-1 for samples diluted 10:1). Point shapes correspond to the tidal stages of high, ebb, and low, while the colors correspond to sampling locations.
Although the highest ENT concentrations in the waterway during the perigean tide events were observed during ebb and low tide, the majority of concentrations were below the EPA’s threshold for safe public use, demonstrating that the observed conditions did not pose a public health hazard. Yet, the observed increase in ENT concentrations during ebb tide in a large and dynamic waterway like Taylor’s Creek demonstrates that tidal floods and stormwater network inundation can serve as a source of fecal contamination to surface waters and are worthy of further investigation as an emerging public health concern, particularly in coastal communities that drain to relatively small tidal creeks with minimal potential to flush or dilute pollutants.
From the results of our robust linear mixed effect models for the perigean tides (Table 1), tidal height at sample collection and time to last high tide were associated with the greatest and second greatest marginal R2 values (0.20 and 0.13, respectively). Though they cannot explain the majority of ENT variance (i.e., they do not have marginal R2 > 0.5), the marginal R2 values of tidal height at sample collection and time to last high tide support prior recommendations by Price et al. (2021) that tidal variables be included and documented as part of coastal ENT monitoring efforts, as tidal effects can play a role in transporting ENT within coastal systems.
For perigean tide conditions, height of the last high tide and rainfall variables had low marginal R2 values and coefficients that were not significantly different from zero (Table 1). Because height of last high tide was not a significant descriptor of ENT concentration variance, we conclude that maximum inundation prior to sampling was not the most important factor explaining the elevated ENT concentrations during ebb tide for our study site. However, height of last high tide should not be discounted as a potential variable to monitor, as other coastal stormwater network configurations may exhibit increase in ENT concentrations with higher tidal heights. Likewise, our results do not necessarily indicate that rainfall is not an important descriptor to ENT variability during the ebb tide and drainage of the stormwater network, as other studies have demonstrated a direct link between rainfall and ENT concentration in coastal waterways (Mallin et al., 2000). The low marginal R2 values of our tested rainfall variables is likely from our sampling design which targeted the tidal dynamics of the 6-7 hour period between the higher high tide and the higher low tide, and from the paucity of rainfall events that occurred during our sampling window.
3.5 Stormwater runoff drives higher ENT concentrations than tidal inundation
During baseline conditions, stormwater runoff drove higher ENT concentrations in the waterway than tidal inundation. Many studies have reported increased FIB concentrations in waterways following rainfall events (Converse et al., 2011; Gonzalez et al., 2012; Parker et al., 2010; Stumpf et al., 2010), and our results support this finding for baseline conditions (i.e., when there is no roadway flooding from complete stormwater inundation by tides). Prior rainfall totals were significantly related to daily ENT variance in the waterway during baseline conditions (Table 1). Specifically, antecedent rainfall over 3, 6, 12, 24, 48, and 72 hours prior to sample collection had larger marginal R2 values (0.10-0.35) in the baseline model than the significant tidal variables (0.04-0.06). While antecedent rainfall predominantly explained the baseline ENT concentrations, rain variables did not significantly explain ENT variance during perigean tides, as discussed above. However, we observed that recent rainfall events coincided with the maximum ENT concentrations observed in Taylor’s Creek during the perigean tide periods (Fig 5). Low tide samples collected on July 15, 2022, which coincided with a compound flood event (rainfall plus reduced capacity in the stormwater network due to tidal inundation), displayed the maximum ENT concentrations observed in the waterway during both perigean tides. These samples were collected 20 hours after a 22-mm, 8-hour rainfall event occurred the previous day, and during the first hour of a 17-mm rainstorm that lasted over 9 hours. Approximately 3.6 mm of rainfall occurred during the first hour of the event.
4 Conclusions
Our findings indicate that chronic stormwater network inundation and tidal flooding caused by SLR present pathways for fecal contamination of coastal surface waters and roadway floodwaters, but research is needed to fully understand whether network inundation and floods create public health hazards. In this study, we consistently observed high ENT concentrations above the US EPA’s single sample maximum threshold for safe public use of recreational waters in tidal and compound floodwaters on a roadway in Beaufort, NC. These measurements were taken from minor floods with minimal overland footprints, leading us to believe the sources of the elevated ENT in the floodwaters were within the stormwater network and could include undrained stormwater runoff, pipe biofilms and sediment, and exfiltrated sewage. However, we did not observe similarly concerning ENT concentrations in Taylor’s Creek, the adjacent coastal waterway, as a function of tide-driven stormwater network inundation. Rainfall-driven runoff produced higher ENT concentrations in Taylor’s Creek than stormwater network inundation and tidal flooding, reinforcing that stormwater runoff is the dominant driver of increased ENT concentrations in the coastal waterway. However, we observed short-lived increases in ENT concentrations in Taylor’s Creek as perigean tides receded from the stormwater network, and we hypothesize the temporary increase in ENT concentrations is due to dilution or flushing given that Taylor’s Creek is a large coastal waterway. In relatively small tidal creeks draining urban centers, the effects of dilution and flushing may not be as apparent, and the role of stormwater network inundation and tidal flooding on ENT concentrations may be more pronounced.
This study provided needed longitudinal water quality observations across a two-month period that included nine tidal floods from two perigean tide events, which allowed us to assess potential public health hazards from tidal flooding and make inferences on the relative roles of stormwater network inundation, tidal flooding, and rain on fecal contamination both in a coastal waterway and roadway floodwaters. However, the small size of the tidal floods we observed prevents us from concluding whether more extensive floods are associated with higher or lower ENT concentrations. Additionally, while our results demonstrate that tidal floods are associated with problematic water quality in both roadway floodwaters and coastal surface waters as was shown in previous studies (Macías-Tapia et al., 2021, 2023; Hart et al., 2020; Price et al., 2021; McKenzie et al., 2021), similar studies as the one presented here should be conducted across more locations to better understand site-specific effects related to stormwater infrastructure and land use. We suggest that future work incorporate sampling across multiple seasons, communities, and flood extents to better understand the prevalence of poor water quality from chronic stormwater network inundation and tidal flooding. Seasonal differences may contribute to the persistence of ENT within the stormwater network, while community differences in infrastructure, topography, land use, and waterway size could contribute to the type and number of ENT sources as well as the effect of tidal flushing on ENT concentrations in the waterway.
Ultimately, the water quality effects of chronic stormwater network inundation and tidal flooding are evolving given the non-stationarity of global climate change and SLR. SLR brings uncertainty to whether tidally-driven inundation and associated pollutant loading will remain at their current levels, as the magnitude of the loading will likely change with increased flooding frequency and damage to infrastructure. For Beaufort and other cities on the Atlantic coast of the US, rates of SLR are 30% higher than the global average of 3 mm/year (Ezer & Atkinson, 2014), meaning that future repetition of this study could produce differing results as tidally-driven inundation becomes more intense and widespread.
Acknowledgments
The authors report no financial conflicts of interest nor affiliations for any author that may be perceived as having a conflict of interest with respect to the results of this paper. We thank David Bennett and Grant Caraway of the North Carolina Maritime Museum for access to their dock, and Dave Eggleston, Melissa LaCroce, and the NC State University Center for Marine Sciences and Technology for use of lab facilities. This work was supported by the U.S. National Science Foundation award number 2047609.
Open Research
 All data and R code are provided as supporting materials with the manuscript for peer review purposes. Following review, the data and code will be archived with Dryad and a permanent DOI will be included in this section.
 References
Allen, T. R., Crawford, T., Montz, B., Whitehead, J., Lovelace, S., Hanks, A. D., Christensen, A. R., & Kearney, G. D. (2019). Linking Water Infrastructure, Public Health, and Sea Level Rise: Integrated Assessment of Flood Resilience in Coastal Cities. Public Works Management & Policy, 24(1), 110–139. https://doi.org/10.1177/1087724X18798380
Befus, K. M., Kurnizki, A. P. D., Kroeger, K. D., Eagle, M. J., & Smith, T. P. (2023). Forecasting Sea Level Rise-driven Inundation in Diked and Tidally Restricted Coastal Lowlands. Estuaries and Coasts, 46(5), 1157–1169. https://doi.org/10.1007/s12237-023-01174-1
Boehm, A. B. (2007). Enterococci Concentrations in Diverse Coastal Environments Exhibit Extreme Variability. Environmental Science & Technology, 41(24), 8227–8232. https://doi.org/10.1021/es071807v
Burkhart, T. H. (Tsung H. S. (2013). Biofilms as sources of fecal bacteria contamination in the stormwater drainage system in Singapore [Thesis, Massachusetts Institute of Technology]. https://dspace.mit.edu/handle/1721.1/82805
Center for Operational Oceanographic Products and Services. (2022). Tide & Currents for Beaufort, Duke Marine Lab, NC - Station ID: 8656483. National Oceanic and Atmospheric Administration. https://tidesandcurrents.noaa.gov/stationhome.html?id=8656483
Converse, R. R., Piehler, M. F., & Noble, R. T. (2011). Contrasts in concentrations and loads of conventional and alternative indicators of fecal contamination in coastal stormwater. Water Research, 45(16), 5229–5240. https://doi.org/10.1016/j.watres.2011.07.029
Division of Water Resources. (2023, March). 2022 North Carolina 303(d) List. North Carolina Department of Environmental Quality. https://edocs.deq.nc.gov/WaterResources/DocView.aspx?dbid=0&id=2738821
Ezer, T., & Atkinson, L. P. (2014). Accelerated flooding along the U.S. East Coast: On the impact of sea-level rise, tides, storms, the Gulf Stream, and the North Atlantic Oscillations. Earth’s Future, 2(8), 362–382. https://doi.org/10.1002/2014EF000252
Gold, A., Anarde, K., Grimley, L., Neve, R., Srebnik, E. R., Thelen, T., Whipple, A., & Hino, M. (2023). Data From the Drain: A Sensor Framework That Captures Multiple Drivers of Chronic Coastal Floods. Water Resources Research, 59(4), e2022WR032392. https://doi.org/10.1029/2022WR032392
Gold, A. C., Brown, C. M., Thompson, S. P., & Piehler, M. F. (2022). Inundation of Stormwater Infrastructure Is Common and Increases Risk of Flooding in Coastal Urban Areas Along the US Atlantic Coast. Earth’s Future, 10(3), e2021EF002139. https://doi.org/10.1029/2021EF002139
Gonzalez, R. A., Conn, K. E., Crosswell, J. R., & Noble, R. T. (2012). Application of empirical predictive modeling using conventional and alternative fecal indicator bacteria in eastern North Carolina waters. Water Research, 46(18), 5871–5882. https://doi.org/10.1016/j.watres.2012.07.050
Habel, S., Fletcher, C. H., Anderson, T. R., & Thompson, P. R. (2020). Sea-Level Rise Induced Multi-Mechanism Flooding and Contribution to Urban Infrastructure Failure. Scientific Reports, 10, 3796. https://doi.org/10.1038/s41598-020-60762-4
Hart, J. D., Blackwood, A. D., & Noble, R. T. (2020). Examining coastal dynamics and recreational water quality by quantifying multiple sewage specific markers in a North Carolina estuary. Science of The Total Environment, 747, 141124. https://doi.org/10.1016/j.scitotenv.2020.141124
Health and Ecological Criteria Division, Office of Science and Technology. (2012, December). Recreational Water Quality Criteria. (EPA-820-F-12-061). Environmental Protection Agency. Office of Water. https://www.epa.gov/sites/default/files/2015-10/documents/rwqc2012.pdf
IDEXX. (2019, May). Standard Test Method for Enterococci in Water Using Enterolert. (ASTM D6503-19). IDEXX. https://www.astm.org/d6503-19.html
Koller, M. (2016). robustlmm: An R Package for Robust Estimation of Linear Mixed-Effects Models. Journal of Statistical Software, 75, 1–24. https://doi.org/10.18637/jss.v075.i06
Lewis, D. J., Atwill, E. R., Pereira, M. das G. C., & Bond, R. (2013). Spatial and Temporal Dynamics of Fecal Coliform and Escherichia coli Associated with Suspended Solids and Water within Five Northern California Estuaries. Journal of Environmental Quality, 42(1), 229–238. https://doi.org/10.2134/jeq2011.0479
Lüdecke, D., Ben-Shachar, M., Patil, I., & Makowski, D. (2020). Extracting, Computing and Exploring the Parameters of Statistical Models using R. Journal of Open Source Software, 5(53), 2445. https://doi.org/10.21105/joss.02445
Lüdecke, D., Ben-Shachar, M. S., Patil, I., Waggoner, P., & Makowski, D. (2021). performance: An R Package for Assessment, Comparison and Testing of Statistical Models. Journal of Open Source Software, 6(60), 3139. https://doi.org/10.21105/joss.03139
Macías-Tapia, A., Mulholland, M. R., Selden, C. R., Loftis, J. D., & Bernhardt, P. W. (2021). Effects of tidal flooding on estuarine biogeochemistry: Quantifying flood-driven nitrogen inputs in an urban, lower Chesapeake Bay sub-tributary. Water Research, 201, 117329. https://doi.org/10.1016/j.watres.2021.117329
Macías-Tapia, A., Mulholland, M. R., Selden, C. R., Loftis, J. D., & Bernhardt, P. W. (2023). Five Years Measuring the Muck: Evaluating Interannual Variability of Nutrient Loads From Tidal Flooding. Estuaries and Coasts, 46(7), 1756–1776. https://doi.org/10.1007/s12237-023-01245-3
Mallin, M. A., Williams, K. E., Esham, E. C., & Lowe, R. P. (2000). Effect of Human Development on Bacteriological Water Quality in Coastal Watersheds. Ecological Applications, 10(4), 1047–1056. https://doi.org/10.1890/1051-0761(2000)010[1047:EOHDOB]2.0.CO;2
McKenzie, T., Habel, S., & Dulai, H. (2021). Sea-level rise drives wastewater leakage to coastal waters and storm drains. Limnology and Oceanography Letters, 6(3), 154–163. https://doi.org/10.1002/lol2.10186
Olds, H. T., Corsi, S. R., Dila, D. K., Halmo, K. M., Bootsma, M. J., & McLellan, S. L. (2018). High levels of sewage contamination released from urban areas after storm events: A quantitative survey with sewage specific bacterial indicators. PLoS Medicine, 15(7), e1002614. https://doi.org/10.1371/journal.pmed.1002614
Parker, J. K., McIntyre, D., & Noble, R. T. (2010). Characterizing fecal contamination in stormwater runoff in coastal North Carolina, USA. Water Research, 44(14), 4186–4194. https://doi.org/10.1016/j.watres.2010.05.018
Price, M. T., Blackwood, A. D., & Noble, R. T. (2021). Integrating culture and molecular quantification of microbial contaminants into a predictive modeling framework in a low-lying, tidally-influenced coastal watershed. Science of The Total Environment, 792, 148232. https://doi.org/10.1016/j.scitotenv.2021.148232
R Core Team (2013). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. ISBN 3-900051-07-0, URL http://www.R-project.org/
Sercu, B., Van De Werfhorst, L. C., Murray, J. L., & Holden, P. A. (2011). Sewage exfiltration as a source of storm drain contamination during dry weather in urban watersheds. Environmental science & technology, 45(17), 7151-7157.
State Climate Office of North Carolina. (2022) Michael J Smith Field Airport (KMRH) hourly and daily rainfall, June-August 2022 [Data set]. NOAA National Weather Service. https://products.climate.ncsu.edu/data/
Strauss, B. H., Ziemlinski, R., Weiss, J. L., & Overpeck, J. T. (2012). Tidally adjusted estimates of topographic vulnerability to sea level rise and flooding for the contiguous United States. Environmental Research Letters, 7(1), 014033. https://doi.org/10.1088/1748-9326/7/1/014033
Stumpf, C. H., Piehler, M. F., Thompson, S., & Noble, R. T. (2010). Loading of fecal indicator bacteria in North Carolina tidal creek headwaters: Hydrographic patterns and terrestrial runoff relationships. Water Research, 44(16), 4704–4715. https://doi.org/10.1016/j.watres.2010.07.004
Sweet, W., Hamlington, B.D., Kopp, R. E., & Zuzak, C. (2022). Global and Regional Sea Level Rise Scenarios for the United States: Updated Mean Projections and Extreme Water Level Probabilities Along U.S. Coastlines (p. 111) [Technical Report]. National Oceanic and Atmospheric Administration. https://aambpublicoceanservice.blob.core.windows.net/oceanserviceprod/hazards/sealevelrise/noaa-nos-techrpt01-global-regional-SLR-scenarios-US.pdf
University of North Carolina - Chapel Hill. (2022). 2022 King Tides Calendar. North Carolina King Tides Project. https://nckingtides.web.unc.edu/