Quantifying the residence time and flushing characteristics of a shallow, back-barrier estuary : application of hydrodynamic and particle tracking models

This paper is not subject to U.S. copyright. The definitive version was published in Estuaries and Coasts 38 (2015): 1719-1734, doi:10.1007/s12237-014-9885-3.Estuarine residence time is a major driver of eutrophication and water quality. Barnegat Bay-Little Egg Harbor (BB-LEH), New Jersey, is a lagoonal back-barrier estuary that is subject to anthropogenic pressures including nutrient loading, eutrophication, and subsequent declines in water quality. A combination of hydrodynamic and particle tracking modeling was used to identify the mechanisms controlling flushing, residence time, and spatial variability of particle retention. The models demonstrated a pronounced northward subtidal flow from Little Egg Inlet in the south to Pt. Pleasant Canal in the north due to frictional effects in the inlets, leading to better flushing of the southern half of the estuary and particle retention in the northern estuary. Mean residence time for BB-LEH was 13 days but spatial variability was between ∼0 and 30 days depending on the initial particle location. Mean residence time with tidal forcing alone was 24 days (spatial variability between ∼0 and 50 days); the tides were relatively inefficient in flushing the northern end of the Bay. Scenarios with successive exclusion of physical processes from the models revealed that meteorological and remote offshore forcing were stronger drivers of exchange than riverine inflow. Investigations of water quality and eutrophication should take into account spatial variability in hydrodynamics and residence time in order to better quantify the roles of nutrient loading, production, and flushing.Funding was provided by the New Jersey Department of Environmental Protection and the Coastal and Marine Geology Program of the U.S. Geological Survey

Are elevation and open-water conversion of salt marshes connected?

Author Posting. © American Geophysical Union, 2020. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geophysical Research Letters 47(3), (2020): e2019GL086703, doi:10.1029/2019GL086703.Salt marsh assessments focus on vertical metrics such as accretion or lateral metrics such as open‐water conversion, without exploration of how the dimensions are related. We exploited a novel geospatial data set to explore how elevation is related to the unvegetated‐vegetated marsh ratio (UVVR), a lateral metric, across individual marsh “units” within four estuarine‐marsh systems. We find that elevation scales consistently with the UVVR across systems, with lower elevation units demonstrating more open‐water conversion and higher UVVRs. A normalized elevation‐UVVR relationship converges across systems near the system‐mean elevation and a UVVR of 0.1, a critical threshold identified by prior studies. This indicates that open‐water conversion becomes a dominant lateral instability process at a relatively conservative elevation threshold. We then integrate the UVVR and elevation to yield lifespan estimates, which demonstrate that higher elevation marshes are more resilient to internal deterioration, with an order‐of‐magnitude longer lifespan than predicted for lower elevation marshes.This study was supported by the USGS through the Coastal Marine Hazards/Resources Program, the National Park Service through the Natural Resource Preservation Program, and the U.S. Fish and Wildlife Service through the Science Support Partnership. Erika Lentz, Elizabeth Pendleton, Meagan Gonneea, Joel Carr, and two anonymous reviewers provided constructive advice on the study. S.F. was partly supported by US National Science Foundation award 1637630 (PIE LTER), 1832221 (VCR LTER). The geospatial data used in this study are published in the Coastal Wetlands Synthesis Products catalog on ScienceBase (https://www.sciencebase.gov/catalog/item/5b73325ee4b0f5d5787c5ff3)

Identifying salt marsh shorelines from remotely sensed elevation data and imagery

© The Author(s), 2019. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Farris, A. S., Defne, Z., & Ganju, N. K. Identifying salt marsh shorelines from remotely sensed elevation data and imagery. Remote Sensing, 11(15), (2019): 1795, doi: 10.3390/rs11151795.Salt marshes are valuable ecosystems that are vulnerable to lateral erosion, submergence, and internal disintegration due to sea level rise, storms, and sediment deficits. Because many salt marshes are losing area in response to these factors, it is important to monitor their lateral extent at high resolution over multiple timescales. In this study we describe two methods to calculate the location of the salt marsh shoreline. The marsh edge from elevation data (MEED) method uses remotely sensed elevation data to calculate an objective proxy for the shoreline of a salt marsh. This proxy is the abrupt change in elevation that usually characterizes the seaward edge of a salt marsh, designated the “marsh scarp.” It is detected as the maximum slope along a cross-shore transect between mean high water and mean tide level. The method was tested using lidar topobathymetric and photogrammetric elevation data from Massachusetts, USA. The other method to calculate the salt marsh shoreline is the marsh edge by image processing (MEIP) method which finds the unvegetated/vegetated line. This method applies image classification techniques to multispectral imagery and elevation datasets for edge detection. The method was tested using aerial imagery and coastal elevation data from the Plum Island Estuary in Massachusetts, USA. Both methods calculate a line that closely follows the edge of vegetation seen in imagery. The two methods were compared to each other using high resolution unmanned aircraft systems (UAS) data, and to a heads-up digitized shoreline. The root-mean-square deviation was 0.6 meters between the two methods, and less than 0.43 meters from the digitized shoreline. The MEIP method was also applied to a lower resolution dataset to investigate the effect of horizontal resolution on the results. Both methods provide an accurate, efficient, and objective way to track salt marsh shorelines with spatially intensive data over large spatial scales, which is necessary to evaluate geomorphic change and wetland vulnerability.This project was supported by the U.S. Geological Survey (USGS) Coastal/Marine Natural Hazards and Resources Program as well as the Massachusetts O ce of Coastal Zone Management under interagency agreement 16ENMALQ006000

Estimating time-dependent connectivity in marine systems

This paper is not subject to U.S. copyright. The definitive version was published in Geophysical Research Letters 43 (2016): 1193–1201, doi:10.1002/2015GL066888.Hydrodynamic connectivity describes the sources and destinations of water parcels within a domain over a given time. When combined with biological models, it can be a powerful concept to explain the patterns of constituent dispersal within marine ecosystems. However, providing connectivity metrics for a given domain is a three-dimensional problem: two dimensions in space to define the sources and destinations and a time dimension to evaluate connectivity at varying temporal scales. If the time scale of interest is not predefined, then a general approach is required to describe connectivity over different time scales. For this purpose, we have introduced the concept of a “retention clock” that highlights the change in connectivity through time. Using the example of connectivity between protected areas within Barnegat Bay, New Jersey, we show that a retention clock matrix is an informative tool for multitemporal analysis of connectivity.New Jersey Department of Environmental Protectio

Development and application of landsat-based wetland vegetation cover and unvegetated-vegetated marsh ratio (UVVR) for the conterminous United States

© The Author(s), 2022. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Ganju, N. K., Couvillion, B. R., Defne, Z., & Ackerman, K. Development and application of landsat-based wetland vegetation cover and unvegetated-vegetated marsh ratio (UVVR) for the conterminous United States. Estuaries and Coasts, (2022), https://doi.org/10.1007/s12237-022-01081-x.Effective management and restoration of salt marshes and other vegetated intertidal habitats require objective and spatially integrated metrics of geomorphic status and vulnerability. The unvegetated-vegetated marsh ratio (UVVR), a recently developed metric, can be used to establish present-day vegetative cover, identify stability thresholds, and quantify vulnerability to open-water conversion over a range of spatial scales. We developed a Landsat-based approach to quantify the within-pixel vegetated fraction and UVVR for coastal wetlands of the conterminous United States, at 30-m resolution for 2014–2018. Here we present the methodology used to generate the UVVR from spectral indices, along with calibration, validation, and spatial autocorrelation assessments. We then demonstrate multiple applications of the data across varying spatial scales: first, we aggregate the UVVR across individual states and estuaries to quantify total vegetated wetland area for the nation. On the state level, Louisiana and Florida account for over 50% of the nation’s total, while on the estuarine level, the Chesapeake Bay Estuary and selected Louisiana coastal areas each account for over 6% of the nation’s total vegetated wetland area. Second, we present cases where this dataset can be used to track wetland change (e.g., expansion due to restoration and loss due to stressors). Lastly, we propose a classification methodology that delineates areas vulnerable to open-water expansion based on the 5-year mean and standard deviation of the UVVR. Calculating the UVVR for the period-of-record back to 1985, as well as regular updating, will fill a critical gap for tracking national status of salt marshes and other vegetated habitats through time and space.This work was supported by the U.S. Geological Survey’s Coastal and Marine Hazards/Resources Program

Salt marsh erosion rates and boundary features in a shallow Bay

Author Posting. © American Geophysical Union, 2016. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Earth Surface 121 (2016): 1861–1875, doi:10.1002/2016JF003975.Herein, we investigate the relationship between wind waves, salt marsh erosion rates, and the planar shape of marsh boundaries by using aerial images and the numerical model Coupled-Ocean-Atmosphere-Wave-Sediment-Transport Modeling System (COAWST). Using Barnegat Bay, New Jersey, as a test site, we found that salt marsh erosion rates maintain a similar trend in time. We also found a significant relationship between salt marsh erosion rates and the shape of marsh boundaries which could be used as a geomorphic indicator of the degradation level of the marsh. Slowly eroding salt marshes are irregularly shaped with fractal dimension higher than rapidly deteriorating marshes. Moreover, for low-wave energy conditions, there is a high probability of isolated and significantly larger than average failures of marsh portions causing a long-tailed distribution of localized erosion rates. Finally, we confirm the existence of a significant relationship between salt marsh erosion rate and wind waves exposure. Results suggest that variations in time in the morphology of salt marsh boundaries could be used to infer changes in frequency and magnitude of external agents.Department of the Interior Hurricane Sandy Recovery program Grant Number: GS2-2D; NSF DEB Grant Number: 0621014; OCE Grant Number: 12382122017-04-2

Physical response of a back-barrier estuary to a post-tropical cyclone

© The Author(s), 2017. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Journal of Geophysical Research: Oceans 122 (2017): 5888–5904, doi:10.1002/2016JC012344.This paper presents a modeling investigation of the hydrodynamic and sediment transport response of Chincoteague Bay (VA/MD, USA) to Hurricane Sandy using the Coupled Ocean-Atmosphere-Wave-Sediment-Transport (COAWST) modeling system. Several simulation scenarios with different combinations of remote and local forces were conducted to identify the dominant physical processes. While 80% of the water level increase in the bay was due to coastal sea level at the peak of the storm, a rich spatial and temporal variability in water surface slope was induced by local winds and waves. Local wind increased vertical mixing, horizontal exchanges, and flushing through the inlets. Remote waves (swell) enhanced southward flow through wave setup gradients between the inlets, and increased locally generated wave heights. Locally generated waves had a negligible effect on water level but reduced the residual flow up to 70% due to enhanced apparent roughness and breaking-induced forces. Locally generated waves dominated bed shear stress and sediment resuspension in the bay. Sediment transport patterns mirrored the interior coastline shape and generated deposition on inundated areas. The bay served as a source of fine sediment to the inner shelf, and the ocean-facing barrier island accumulated sand from landward-directed overwash. Despite the intensity of the storm forcing, the bathymetric changes in the bay were on the order of centimeters. This work demonstrates the spectrum of responses to storm forcing, and highlights the importance of local and remote processes on back-barrier estuarine function.Department of Interior Hurricane Sandy Recovery progra

Toward a comprehensive water-quality modeling of Barnegat Bay : development of ROMS to WASP coupler

Author Posting. © Coastal Education and Research Foundation, 2017. This article is posted here by permission of Coastal Education and Research Foundation for personal use, not for redistribution. The definitive version was published in Journal of Coastal Research SI78 (2017): 34-45, doi:10.2112/SI78-004.1.The Regional Ocean Modeling System (ROMS) has been coupled with the Water Quality Analysis Simulation Program (WASP) to be used in a comprehensive analysis of water quality in Barnegat Bay, New Jersey. The coupler can spatially aggregate hydrodynamic information in ROMS cells into larger WASP segments. It can also be used to resample ROMS output at a finer temporal scale to meet WASP time-stepping requirements. The coupler aggregates flow components, temperature, and salinity in ROMS output for input to WASP via a hydrodynamic linkage file. The coupler was tested initially with idealized cases designed to verify the water mass balance and conservation of constituent mass using one-to-one and one-to-many connectivity options between segments. A realistic example from the Toms River embayment, a subdomain of Barnegat Bay, was used to demonstrate the functionality of the coupling. A WASP eutrophication model accounting for dissolved oxygen (DO), nitrogen, and constant phytoplankton concentrations was applied to explore the distribution and trends in DO and nitrogen in the embayment for the period of July–August 2012. Results of DO modeling indicate satisfactory agreement with measurements collected at in-bay stations and also indicate that this coupled approach, despite substantial differences in spatiotemporal discretization between the models, provides adequate predictive capabilities.Funding was provided by the NJDEP and the Coastal and Marine Geology Program of the USGS

Spatial distribution of water level impacting back-barrier bays

© The Author(s), 2019. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Aretxabaleta, A. L., Ganju, N. K., Defne, Z., & Signell, R. P. Spatial distribution of water level impacting back-barrier bays. Natural Hazards and Earth System Sciences, 19(8), (2019): 1823-1838, doi: 10.5194/nhess-19-1823-2019.Water level in semi-enclosed bays, landward of barrier islands, is mainly driven by offshore sea level fluctuations that are modulated by bay geometry and bathymetry, causing spatial variability in the ensuing response (transfer). Local wind setup can have a complementary role that depends on wind speed, fetch, and relative orientation of the wind direction and the bay. Bay area and inlet geometry and bathymetry primarily regulate the magnitude of the transfer between open ocean and bay. Tides and short-period offshore oscillations are more damped in the bays than longer-lasting offshore fluctuations, such as a storm surge and sea level rise. We compare observed and modeled water levels at stations in a mid-Atlantic bay (Barnegat Bay) with offshore water level proxies. Observed water levels in Barnegat Bay are compared and combined with model results from the Coupled Ocean–Atmosphere–Wave–Sediment Transport (COAWST) modeling system to evaluate the spatial structure of the water level transfer. Analytical models based on the dimensional characteristics of the bay are used to combine the observed data and the numerical model results in a physically consistent approach. Model water level transfers match observed values at locations inside the bay in the storm frequency band (transfers ranging from 50 %–100 %) and tidal frequencies (10 %–55 %). The contribution of frequency-dependent local setup caused by wind acting along the bay is also considered. The wind setup effect can be comparable in magnitude to the offshore transfer forcing during intense storms. The approach provides transfer estimates for locations inside the bay where observations were not available, resulting in a complete spatial characterization. An extension of the methodology that takes advantage of the ADCIRC tidal database for the east coast of the United States allows for the expansion of the approach to other bay systems. Detailed spatial estimates of water level transfer can inform decisions on inlet management and contribute to the assessment of current and future flooding hazard in back-barrier bays and along mainland shorelines.This work was supported by the US Geological Survey, Coastal and Marine Hazards/Resources Program

Role of tidal wetland stability in lateral fluxes of particulate organic matter and carbon

Author Posting. © American Geophysical Union, 2019. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research-Biogeosciences 124(5), (2019): 1265-1277, doi:10.1029/2018JG004920.Tidal wetland fluxes of particulate organic matter and carbon (POM, POC) are important terms in global budgets but remain poorly constrained. Given the link between sediment fluxes and wetland stability, POM and POC fluxes should also be related to stability. We measured POM and POC fluxes in eight microtidal salt marsh channels, with net POM fluxes ranging between −121 ± 33 (export) and 102 ± 28 (import) g OM·m−2·year−1 and net POC fluxes ranging between −52 ± 14 and 43 ± 12 g C·m−2·year−1. A regression employing two measures of stability, the unvegetated‐vegetated marsh ratio (UVVR) and elevation, explained >95% of the variation in net fluxes. The regression indicates that marshes with lower elevation and UVVR import POM and POC while higher elevation marshes with high UVVR export POM and POC. We applied these relationships to marsh units within Barnegat Bay, New Jersey, USA, finding a net POM import of 2,355 ± 1,570 Mg OM/year (15 ± 10 g OM·m−2·year−1) and a net POC import of 1,263 ± 632 Mg C/year (8 ± 4 g C·m−2·year−1). The magnitude of this import was similar to an estimate of POM and POC export due to edge erosion (−2,535 Mg OM/year and − 1,291 Mg C/year), suggesting that this system may be neutral from a POM and POC perspective. In terms of a net budget, a disintegrating wetland should release organic material, while a stable wetland should trap material. This study quantifies that concept and demonstrates a linkage between POM/POC flux and geomorphic stability.Use of brand names is for identification purposes only and does not constitute endorsement by the U.S. Government. This study was supported by the USGS Coastal and Marine Geology Program, the Department of the Interior Hurricane Sandy Recovery program (GS2‐2D), and the USGS Mendenhall Post‐doctoral Research Program. Viktoria Unger and Paula Zelanko are acknowledged for field and lab assistance. Core collection was funded under NJ SeaGrant/NOAA Grant 6210‐0011. Gil Pontius provided helpful feedback on statistical measures. Kevin Kroeger and two anonymous reviewers provided constructive reviews of the manuscript. All time series and water sample data are available at the U.S. Geological Survey's Oceanographic Time‐Series Data Collection (at https://stellwagen.er.usgs.gov/).2019-10-2