Correction: Wetland Elevation Change Following Beneficial Use of Dredged Material Nourishment

Brian Harris^1*Aleksandra Ostojic¹Lenore Tedesco²Keith VanDerSys³Susan Bailey¹Justin Shawler¹Navid Jafari⁴Monica Chasten⁵

• ¹Engineer Research and Development Center (ERDC), Vicksburg, United States
• ²The Wetland Institute, Stone Harbor, United States
• ³University of Pennsylvania, Philadelphia, United States
• ⁴Texas A&M University, College Station, United States
• ⁵US Army Corps of Engineers, Washington, United States

1. Introduction Coastal wetlands are highly productive, dynamic systems that provide diverse ecosystem services like wildlife habitat for migratory bird populations, carbon and nutrient sequestration, and storm protection (Mitsch and Gosselink 1993; Costanza et al. 1998; Day et al. 2008; Reddy et al. 2010; Barbier et al. 2011). They exist at the boundary between land and sea with their distribution limited by inundation and erosion at lower elevations and vegetative competition at the landward boundary (Mendelssohn et al., 1981; Davis et al. 2022). Many of these essential ecosystems are being impacted by the combined effects of hydrologic restrictions, sea level rise, subsidence, and fragmentation that result in reductions in productivity and widespread losses of vast acres as vegetated surfaces transition to open water (Day et al. 2000; Fattet et al. 2011; Harris et al. 2020a; Davis et al. 2022; Cadigan et al. 2023; Harris et al. 2023). For example, wetland loss in the Mississippi River delta continues at an approximate rate of 28 km2 per year due to natural processes such as delta abandonment, subsidence, sea level rise, wave erosion, herbivory, hurricanes, and geologic faults, in addition to anthropogenic activities including levee building, hydrologic isolation, land reclamation for agriculture and industry, canal dredging and associated salt water intrusion, boat-induced shoreline erosion, herbivory by invasive species, and pollution (Boesch et al., 1994; Turner 1997; Day et al., 2000; Day et al. 2020). To rehabilitate these deteriorating environments, hydraulically dredged material (i.e., sediment) can be strategically deposited on degraded wetlands to increase surface elevations and reconnect discontinuous wetlands. This application of uncontaminated dredged material is a form of direct nourishment that could potentially increase the rate of SLR that a marsh can withstand (Davis et al. 2022; Harris et al. 2024). This gain in elevation capital is crucial to long-term wetland stability as it controls the hydroperiod (i.e., flooding frequency and duration) and thereby salinity levels, nutrient availability, plant productivity, and redox potential (Croft et al 2006; Day et al. 2011; Morris et al. 2013; Jafari et al. 2019a; Harris et al. 2020b). Coastal wetlands occupy elevations between mean sea level (MSL) and mean high water (MHW) (McKee and Patrick 1988), with optimum growth occurring near the middle of the range (Morris et al. 2002; Morris et al. 2013) while generally wetlands at lower elevations relative to MHW are in more imminent danger of drowning and converting to open water (Davis et al. 2022). Common sources of dredged materials are waterway and port maintenance events. Historically, this material was taken to upland confined disposal facilities or offshore disposal sites, which permanently removes the sediment from the system. Over the past decades, there has been a shift away from this disposal mindset towards the beneficial use of dredged material (BUDM) to bolster coastal environments. The United States Army Corps of Engineers (USACE) is currently seeking to beneficially use 70% of dredged material by 2030, equating to 100 to 160 million m3within the continental United States annually. Despite the benefits that BUDM offers for both navigation and flood-protection missions, significant knowledge gaps exist on the constructability of wetland nourishments and long-term resiliency of the environment. Coastal managers are faced with a number of construction challenges for BUDM such as: (1) significant uncertainties in dredged sediment characteristics (e.g., slurry density, fines and organic matter concentrations) (Harris et al. 2024; Harris and Shawler 2024); (2) the placement areas commonly being comprised of easily disturbed, highly heterogeneous deltaic soils (Jafari et al. 2018, 2019ab); (3) dredged material delivery methods and positions (Harris and Shawler 2024); and (4) containment efficacy (Singh and Joffrion 2017). Following dredged material placement, wetland surface elevations and biogeochemistry can be drastically altered due to the compression of highly organic sediment that can lead to a temporal loss in elevation (VanZomeren and Piercy 2020; Harris et al. 2024), while silt and clay-dominated dredged material, particularly acid sulfate soils, may hinder revegetation due to oxidation of iron sulfides (Berkowitz and VanZomeren 2020; Puchkoff and Lawrence 2022; Raposa et al. 2022), all resulting in uncertainty in vegetation response relative to natural reference sites (Raposa et al. 2023) or target restoration goals. These uncertainties can be problematic to the long-term stability of the environment and thus the ability to provide desired ecosystem services. To address these knowledge gaps, a multi-phase wetland BUDM nourishment project was investigated through a combination of in-situ, laboratory, and geo-spatial methods pre-and post-placement. Sediment testing and subsurface investigations were performed to develop consolidation models and predict future wetland surface elevations, which were compared to field data, while orthomosaics were used to monitor vegetation cover change across the site. Lastly, these validated consolidation models were combined with the Coastal Wetland Evolution Model (CWEM) to assess future site resiliency to sea level rise. This study investigates the following four questions: (1) what is the surface elevation gained from dredged material applied on a well-drained marsh island, (2) what is the magnitude of dredged sediment and wetland foundation consolidation, (3) how rapidly did vegetation re-establish, and (4) what are the potential long-term responses of the nourished system to rising sea levels? 2. Methodology 2.1 Study Site Description The study site is a 5.4-hectare island located on the bay side of Avalon, New Jersey, USA within a 17-km2 tidal marsh complex and is adjacent to the New Jersey Intracoastal Waterway (NJIWW), Fig. 1 (Harris et al. 2020). Tides in this back-bay environment are semi-diurnal with a normal tidal range of approximately 1.3-m (NOAA 2024). Sturgeon Island is located nearly 7 kilometers from Townsend Inlet to the north and nearly 8 kilometers from Hereford Inlet to the south, likely indicating minimal inlet-derived sediment reaches the site. The backbarrier marsh and channels adjacent to Sturgeon Island experience small waves typically less than <0.25 m, slow currents (~0.1 m/s), and low suspended sediment concentrations (~10-23 mg/L), indicating minimal natural mineral sediment contributions to the island (Fall et al. 2021; Coleman et al. 2022). The island, prior to the BUDM nourishments, ranged in elevation from 0.5-m to 1.3-m NAVD88, with the highest elevations found at the center of the island. The northern half of the island drains to a tidal creek on the eastern side of the island. The island is primarily vegetated by the cordgrass Spartina alterniflora and the invasive common reed Phragmites australis, with tall-form S. alterniflora found between 0.2-m and 0.6-m elevations, short-form S. alterniflora between 0.5 m and 0.7 m, and Phragmites found at elevations >0.9-m NAVD88. The woody shrub Iva frutescens can also be found interspersed from 0.5-m to 0.9-m NAVD88. Sturgeon Island has historically hosted high nesting concentrations of wading birds and is an important nesting site in the state of New Jersey. At least four species of wading birds currently nest there, and the majority of these species are identified as Species of Greatest Conservation Need by New Jersey's State Wildlife Action Plan (NJDEP, 2018). Fig. 1. (a) Location of Sturgeon Island, New Jersey, USA, (b) in relation to the New Jersey Intracoastal Waterway (NJIWW), and (c) dredge pipeline and containment layout. The goal of this BUDM nourishment was to enhance the surface elevation on Sturgeon Island to support nesting bird habitat and stabilize the low-lying northern portion of the island. Dredged material was placed over three phases. Phase 1 utilized the 1,300 horsepower (hp) cutterhead dredge, Fullerton (Barnegat Bay Dredging Co.) with a 0.36-m. (14in.) discharge pipeline (Welp et al. 2024; Perkey et al. 2024a; Harris et al. 2024). Material was dredged from the NJIWW from March 16 to 19, 2020 with the discharge pipe positioned on the northern end of the island and equipped with a wye-valve that directed material across two outfall areas using 0.31-m. (12-in.) pipes. During Phase 2, dredged material was placed directly on the island using the 400-hp cutterhead dredge, Montgomery (Barnegat Bay Dredging Co.), with a 0.31-m. (12in.) discharge pipeline from September 9 to 26, 2020. The discharge pipeline was positioned directly on the wetland surface along the northern end of Sturgeon Island. Lastly, Phase 3 utilized the Fullerton (Barnegat Bay Dredging Co.) to place dredged material on the northwestern portion of the island from September 27 to October 31, 2022, that was contained using hydro-dams to construct a berm feature. However, for this study, this latest nourishment was not considered since it was two years after the initial placements, contained, and impacted only a small portion of the larger nourishment area. 2.2. Sediment Testing To determine dredged material characteristics, 1-m long soil cores were collected within the NJIWW in March 2020, prior to the start of dredging, using a 0.15-m inner diameter corer. Sediment was extruded from the core, placed in buckets, and shipped back to a laboratory to be processed for gravimetric moisture content [w(%)], organic content via loss on ignition (LOI) (ASTM D2974-20e1), grain size distribution (GSD), and Atterberg Limits. A column settling test was conducted with NJIWW material to predict settling behavior during initial placement and dewatering (Thackston et al. 1988) and a combination of self-weight consolidation testing and standard oedometer consolidation testing was used to determine consolidation magnitudes under its own weight (Cargill 1986). Results from the self-weight consolidation test provided the void ratio (e), effective vertical stress (σ'v), and hydraulic conductivity (kv) relationships for use in the long-term consolidation modeling of the dredged material. To determine wetland foundation characteristics (i.e., placement area), soil samples were collected in March 2020, prior to the start of dredging, with 7.62-cm diameter PVC cores and processed for w(%), LOI, and GSD. Sixty (60) cone penetrometer tests (CPTs) were performed across the northern portion of Sturgeon Island to determine site stratigraphy, and the data was processed using a Matlab script to correct the soundings and smooth the data without a distortion of the signal tendency (Savitzky and Golay 1964; Cadigan et al. 2020; Harris et al. 2021; Cadigan et al. 2022). The CPT equipment was specifically developed for manual use in ultra-soft deltaic and estuarine soils and is capable of providing a consistent shear strength profile to determine subsurface micro-stratigraphy, root-depths, and the presence of shells (Jafari et al. 2019a; Harris et al. 2021). An in-depth description of the CPT methodology and data interpolation is provided in Jafari et al. (2019a). Soil cores were collected to perform standard oedometer consolidation tests (ASTM D2435) to determine the void ratio (e), effective stress (σ'v), and hydraulic conductivity (kv) relationships. However, due to laboratory access restrictions during Coronavirus-19, no consolidation tests were performed on cores collected on Sturgeon Island. To supplement this, two wetland foundation consolidation tests from a previous wetland BUDM site ~400 meters northeast of the site were utilized. This was not an ideal solution, but both sites were within the same tidal marsh network and within the same vegetation type and range of elevation. In July 2021, following Phases 1 and 2, a Russian peat corer was used to collect disturbed soil samples and record the thickness of dredged material based on visually distinct horizons between the deposited material and in-situ wetland foundation. A total of 33 samples were collected and analyzed for w(%), LOI, and GSD. Surface elevation at each sample location was determined via RTK-GPS and a current elevation of underlying wetland foundation could be inferred by subtracting dredged material thicknesses. This inferred surface is relatively coarse when compared to surface elevations, but it resolved approximate magnitudes of foundation settlement across the nourishment area. 2.3. Elevation Monitoring Airborne LiDAR (Laser imaging, Detection, And Ranging) for Sturgeon Island was collected by the Joint Airborne LiDAR Bathymetry Technical Center of Expertise (JALBTCX) in July 2020 as part of the 2020 USACE USGS Topobathy LiDAR: Cape May & Atlantic City, NJ collection. Data was collected using the Coastal Zone Mapping and Imaging LiDAR (CZMIL) system in support of the United States Geologic Survey (USGS) to assess elevation changes along the NJ coast. The final product used for analysis was a 1-m resolution digital elevation model (DEM) based on classified point data compiled to meet a 1-m horizontal accuracy at the 95% confidence level and a 0.20-m vertical accuracy at the 95% confidence level for topographic data. Uncrewed Aircraft System (UAS) LiDAR for Sturgeon Island was collected by the University of Washington NSF Natural Hazards Engineering Research Infrastructure (NHERI) Reconnaissance Experimental Facility (RAPID) in July 2021. Data was collected using the Phoenix LiDAR Systems miniRanger with 15 mm vertical accuracy. Flight processing was performed using Inertial Explorer and horizonal accuracy tested and corrected against USACE GNSS ground data. Point cloud processing and classification was performed using Phoenix Lab's Spatial Explorer software. The final product was a 0.3 m resolution bare-earth digital elevation model (DEM) for Sturgeon Island. Table 1 in the appendix provides a summary of survey types, dates, and specific equipment used within this study. 2.4. Land Cover Classification Imagery was acquired annually in late summer over four years using a DJI Mavic 3 RTK-enabled drone equipped with a 5-band multispectral sensor. Exact dates are shown in Table 1. All flights were flown at 61-m elevation with an 80% sidelap and 80% frontlap, resulting in a ground sample distance of approximately 1.4-cm. Automated flight pathing was used for all image collections. Flight times were coincident with solar noon, low cloud cover, and low tides. Prior to the flights, differentially corrected ground control points (GCP) were measured with an RTK GNSS receiver (Emlid RS2+) to increase the spatial accuracy of the resultant multispectral imagery. Photo stitching of the imagery was processed using Pix4Dmapper. Post-calibration of multispectral data was performed using a calibrated reference panel before and after the flight. Additionally, downwelling light sensor (DLS) data was used to improve reflectance calibrations during changing light conditions mid-flight. Subsequent multispectral orthomosaics were classified into 5-class land cover sets in ArcGIS Pro through a pixel-based supervised classification process using Support Vector Machine (SVM) algorithms. Accuracy for all classification types exceeded 85%. SVM-based pixel classifiers have been shown to better discriminate salt marsh plant species than other machine learning algorithms such as Random Forest or Maximum Likelihood (Norris et al. 2024). Classification accuracy was improved by incorporating elevation and Normalized Difference Red Edge (NDRE) bands into 5-band multispectral composites. The additional bands greatly improved the differentiation of species with similar spectral characteristics, like the two different forms of salt marsh cord grass. Comparative changes across the multiyear land cover data were calculated using a customized Zonal Statistics methodology developed in Modelbuilder. 2.5. Consolidation Modeling To predict consolidation within the dredged sediment and underlying foundation material (i.e., native marsh), two coupled models were performed to determine (1) self-weight consolidation of the dredged sediment and (2) consolidation of the in-situ wetland due to the applied load of the dredged material. The dredged material was modeled using USACE's Primary Consolidation Secondary Compression and Desiccation of Dredged Fill (PSDDF) model that is based on finite strain theory because large strains are expected within the dredged fill (Terzaghi et al. 1996; Stark et al. 2005; Jafari et al. 2019). The dredged material was bounded to the top by atmospheric conditions and an assumed impermeable lower boundary to simulate the underlying wetland material which was predominately clay. The wetland foundation layer consisted of a vegetated root mat and underlying clay layer and was modeled using Terzaghi's 1-D consolidation theory, which uses small strain theory because the magnitude of settlement is expected to be negligible compared to the original foundation thickness (Terzaghi et al. 1996). The wetland foundation layer was assumed to be single-drained (i.e., drainage in one direction) and was bounded by an underlying sand layer that was assumed to be incompressible and freely draining. Four specific sites were chosen for consolidation modeling. These sites were selected to cover a range of initial elevations (0.21-m to 0.67-m, NAVD88); depth to sand (1.19-m to 1.65-m), and approximate dredged material thicknesses for both nourishments (low: 0.15-m at Site A; high:0.52-m at Site D). Dredged material thicknesses are approximate values as site access was restricted immediately after dredging concluded. Table 2 in the appendix provides a summary of values for each site. A range of dredged material and foundation surface elevations was modeled to hindcast a highly compressible, quickly draining scenario and less compressible, slowly draining scenario to encompass the uncertainty in laboratory results. Due to limited field data for the foundation layer, only the compression index (Cc) was varied and not the permeability. Permeability values for the dredged material layer ranged by an order or magnitude in both ways to account for uncertainties in laboratory derived measurements (Tavenas et al. 1983; Jafari et al. 2019). The Cc ranged from a less compressible Cc = 0.25, determined from standard oedometer results to a more compressible Cc = 2.10, determined from moisture content relationships (Cc = w(%)/100; Mesri and Ajlouni 2007) collected on Sturgeon Island prior to dredged material placement. The consolidation models were calibrated through an iterative process of matching overall dredged material elevations, which includes consolidation of the dredged material and underlying foundation, to the field survey points and the July 2021 soil core data which provided the elevation of the foundation, allowing for the calibration of the foundation layer. 2.6. Wetland Projections The Coastal Wetland Equilibrium Model (CWEM), formerly known as the Marsh Equilibrium Model (MEM), was used to predict future marsh function (Morris et al. 2002; Morris et al. 2021). Marsh equilibrium theory states that a dynamic equilibrium exists between marsh platform elevation and sea level, and that marshes are continuously moving in the direction of that equilibrium (Alizad et al. 2016). MEM uses a polynomial formulation for salt marsh productivity and accounts for inputs of suspended sediments and implicitly for the in-situ input of organic matter to the accreting salt marsh platform (Alizad et al. 2016). CWEM combines organic (i.e., marsh elevation vegetation thresholds, peak aboveground biomass, belowground turnover rate, root depth) and inorganic (i.e., water levels, sea level rise rates, suspended sediment concentrations, marsh elevation) inputs to model marsh accretion at given elevations throughout time (Morris et al. 2012; Schile et al. 2014 Alizad et al. 2016). Based on long-term measurements of sediment accretion and marsh productivity, Morris et al. (2002) developed CWEM (then MEM) and calibrated the program at a Spartina alterniflora dominated marsh at North Inlet, SC (Schile et al. 2014). The model assumes that plant productivity is constrained by upper and lower elevation limits and there is an optimum elevation for growth within the tidal frame (Morris et al. 2002; Schile et al. 2014). Ultimately, CWEM was used to determine future standing biomass, inundation intervals, and elevations with time. Site C was selected to be used within CWEM as dredged material and compressible foundation thicknesses were close to an average among the four sites modeled for consolidation . An initial SLR rate of 0.4 cm/yr was selected (NOAA 2024) and to account for variations in SLR projections, low, moderate, and high values were used, corresponding to 0.8 cm/yr, 1.0 cm/yr, and 1.2 cm/yr, respectively (Kopp et al. 2019). Table 3 in the appendix provides a summary of physical and biological inputs used and their respective references for the pre-BUDM placement at Site C and moderate SLR scenario. 3. Results 3.1. Sediment Analysis The general subsurface stratigraphy across the site was a coarse sandy foundation (CPT refusal), overlain by a compressible fine-grained layer, and a vegetated root layer of variable thickness (Fig. 2) (Harris 2020; Harris et al. 2024). Figure 2 shows the depth to sand layer, interpolated from CPT soundings, alongside three example CPT shear strength profiles moving from the center of the island northward. The uppermost vegetated mat layer was signified by a bulb-shaped shear strength that peaked at 174-kPa at CPT-A and extended to a depth of 32-cm. This layer was underlain by a weaker homogenous compressible clay layer that extended to a depth of 43-cm (i.e., 11-cm thick). Beneath this clay layer was a sandy foundation, which significantly increased in shear strength (reaching >100 kpa) causing refusal of the CPT. Approximately 35-m north at CPT B, the vegetated mat extended to 31-cm, similar to the previous site, but the clay layer was 79-cm thick (down to 110-cm below ground surface). Further north, (100-m from CPT A), there was a 35-cm thick vegetated layer, underlain by a 102-cm thick layer (down to 137-cm below ground surface), below which was the sand foundation. Fig. 2. Site Geotechnics. (a) Location of CPTs (black dots) and interpolated depth to underlying sand foundation across the placement site. (b) – (d) Example shear strength profiles and interpreted stratigraphy moving northward from the center of the island. Profile letters correspond to stars on the map. (e) Void ratio-vertical effective stress and (f) void ratio-hydraulic conductivity relationships for dredged fill and foundation materials. Point data was derived from laboratory testing and dashed lines represent interpolated laboratory data. The e – σ'v relationship, shown in Figure 2e, controls the magnitude of consolidation. The slope of the linear portion of the relationship is the compression index (Cc). However, a unique Cc cannot be computed for the dredged material sample because the e – σ'v relationship is highly nonlinear over several magnitudes of σ'v, where higher values of Cc lie in the lower σ'v range that decrease continuously with increasing σ'v (Terzaghi et al. 1996). The Cc is influenced by soil composition and structure, and a higher Cc indicates a more compressible soil. For example, a peat deposit can produce Cc between 2 to 12 due to large volumes of water being held within the soil void space (Mesri and Ajlouni 2007). The foundation samples (Foundation 1 and 2) from the previous BUDM wetland site ~400-m northeast of the site produced Cc values of 0.25 and 0.55, which is uncharacteristic for typical fine-grained wetland soils. Moisture contents w(%) determined from soil cores and grab samples across Sturgeon Island prior to the nourishment ranged from 173% and 210% which equates to Cc values of 1.7 and 2.1 using Mesri and Ajlouni (2007) correlation [Cc = w(%)/100]. Ultimately, the dredged material consolidation data produced higher Cc's than the foundation soils which was expected due to large volumes of water entrained within the pore space of dredged material during the dredging process which increases void space (i.e., porosity). The e – kv relationship, shown in Figure 2f, controls the time-rate at which soil consolidates due to the expulsion of excess pore pressure from within the void space of the soil. Permeability is a function of void space, the shape of the voids related to the shape, specific surface, arrangement of the soil particles, and the tortuosity of the flow path which directly represents the distribution of the voids (Mesri and Rokhsar 1974; Mesri and Tavenas 1983). However, this series of variables that control permeability are not easily quantified. The permeability values presented in Figure 2f were not directly measured but inversely calculated from the coefficient of consolidation (cv) and volume compressibility (mv). Previous studies have shown that indirectly evaluated permeabilities can be off by as much as an order of magnitude (Jafari et al. 2018), while others state they should be disqualified due to large errors produced by abusive assumptions in the methods of interpreting the test results on the basis of Terzaghi's consolidation theory (Tavenas et al. 1983a). As a result, the installation of field instrumentation, such as piezometers and cone penetrometer dissipation tests, to verify the in-situ permeability values is recommended. A linear e – kv relationship is referred to as the permeability change index (Ck) and is an acceptable representation of the behavior of most natural soft clay and organic soils for volumetric strains under 20% (Tavenas et al. 1983b). However, the data produced by the dredged material self-weight consolidation test resulted in strains in excess of 50%, thus a linear relationship is only observed in the standard oedometer data from the foundation samples and a portion of the dredged fill data (Fig. 2f). 3.2. Surface Mapping Pre-nourishment surface elevations, vegetation type, root-mat thickness, and root-mat peak shear strength are shown in Figure 3. Within the area of interest, the elevations are highest along the western side of the island at 0.7 m NAVD88 and gently slope downward to 0.3 m on the eastern side (see Fig. 3a). As the site drains from west to east, a tidal creek is located on the eastern edge of the island. Within the nourishment area, there were two dominant types of vegetation pre-nourishment, short-and tall-form S. alterniflora. The short-form is found at the higher elevations along the western part of the island and transitions to tall-form in the lower elevations along the eastern half (Fig. 3b). The root-mat was on average ~32-cm thick but was thicker at lower surface elevations adjacent to the tidal creek. This lower elevation, and thus increased inundation frequency, also resulted in lower shear strengths of ~80 kPa as compared to the 180 to 280-kPa of the short-form S. alterniflora. Surface elevations along four transects are shown in Figure 4 throughout the study period. March 2020 data (blue in Figures 4b-e) captures the surface elevations prior to BUDM. The March 2020 transects indicate a gradual slope from west to east, with the west end of the island at approximately 0.63-m NAVD88 while the east end of the island was approximately 0.27-m NAVD88. July 2020 data (purple in Figures 4b-e) captures conditions following the Spring 2020 BUDM, which also serve as pre-nourishment elevation conditions for the Fall 2020 BUDM. Average wetland surface elevations increased from 0.41 ± 0.15-m NAVD88 pre-placement to 0.60 ± 0.13-m NAVD88 following the Spring 2020 placement, and a gradual slope from 0.74-m to 0.37-m NAVD88 was still noted on the west and east ends of the transect, respectively. July 2021 data (red in Figure 4b-e) captures surface elevations following the Spring and Fall 2020 BUDM placements and the region between this surface and the March 2020 data is shaded gray to represent the placed dredged material. Average wetland surface elevation for the July 2021 transects was 0.67 ± 0.09-m NAVD88, which highlighted a lower average elevation gain following this placement. The standard deviation is also lower, indicating this second placement filled in lower spots and reduced the average slope (i.e., the wetland surface became flatter). An interpolated foundation surface (i.e., previous wetland surface) derived from sediment cores (orange points Figure 4a; vertical black lines in Figures 4b-e) was compared to the March 2020 initial elevation data to calculate the amount of foundation consolidation due to the induced dredged material loading. Transects 1, 2, and 3 saw an average of 0.09-m of consolidation (0.07-m if Transect 4 included) with values as high as 0.17-m on Transect 2. Fig. 4. (a) Map of soil core transects, and (b-e) Elevation transects across the placement area from March 2020, July 2020, and July 2021. Solid black lines and dashed black lines indicate the depth of dredged fill at coring locations and the interpolated wetland foundation from cores, respectively. Brown and gray shaded regions represent native wetland foundation and dredged material, respectively. A summary of orthomosaic surveys and land cover classifications post BUDM are shown in Figure 5. The 2020 data occurred in-between the Spring and Fall nourishments and showed a majority of the nourishment area to be unvegetated. At this point, the majority of the BUDM nourishment area (denoted by the red boundary) was classified as bare earth (65%) while the low marsh portion comprised 34% and occurred at the distal portions of the nourishment area. In 2021, approximately 1 year after the Fall 2020 placement, the bare earth area increased to 77% of the nourishment area while the low marsh portion decreased to under 22% as additional dredged material was deposited the previous September. By 2022, two years post nourishment, the vegetation began to re-establish. The percentage of low marsh area increased to 43% while the bare earth decreased to 56%. In 2023, three years post nourishment, a further increase in low marsh was noticed, which now occupied over half the nourishment area (52%) while the bare earth decreased to 33%. There was a slight increase in open water on the site from 2020 to 2022 (0.1% to ~2%), and a significant increase three years post nourishment to nearly 10% of the nourishment area. Wrack was identified within the nourishment area but never occupied a consistent or significant area (<1%). Fig. 5. Orthomosaic surveys and land cover classifications from September 14, 2020; August 23, 2021; August 26, 2022; and September 9, 2023. Land cover classifications consisted of water, bare earth, wrack, low marsh, and high marsh with a red boundary line denoting the approximate nourishment area. For 2023 UAS Land Cover* the percentages exclude the Phase 3 placement area since additional dredged material was placed here, altering the vegetation cover. 3.3. Consolidation Modeling Survey data of wetland surface elevations are overlayed on estimated ranges of dredged material surface elevations and underlying foundation layers and their respective calibrated surfaces at four (4) sites across the nourishment area in Figure 6. Survey data of the wetland surface show a sharp increase in elevation immediately following the dredged material placements followed by a gradual decrease in elevation. This decrease in elevation is primarily driven by consolidation, where a majority of elevation loss occurs within the first 3 to 6 months. Consolidation occurs because the weight of the dredged material generates excess pore water pressure within the dredged material and underlying foundation. This excess pore pressure dissipates from the soil pore space, reducing soil porosity (i.e., void ratio) until primary consolidation is completed. The shaded ranges of modeled surface and foundation elevation represent the variations of highly compressible dredged material and underlying foundation as the lower band and a less compressible media for the upper band. In general, the field data fell within this potential range. For Sites A-C, the calibrated surface elevations fell within the more compressible side of the band while at D they were on the less compressible end, which could be attributed to a shallower sand confining layer. A summary of base and calibrated conditions for Sites A-D are shown in Table 4 in the appendix. Fig. 6. (Left) Map of sites used for consolidation model. (Right) At each site, survey elevation (average values within 2.5-m radius of each point ± 1 standard deviation) overlayed on modeled wetland surface and foundation consolidation ranges and calibrated surface elevations. The thickness of dredged material was determined during the July 2021 survey via Russian Peat Corer borings. A comparison of approximate dredged material thickness and surface elevation gain derived from field data is shown in Figure 7. The data indicate a positive correlation between approximate dredged material thickness and elevation gain. A linear trend shows that nearly 66% of gained elevation post placement was lost due to self-weight consolidation of the dredged material and primary consolidation of the foundation soils two years post BUDM nourishment. This relationship is dependent on the environment, soil properties of the dredged material and foundation, site geometries (size, slopes, containment), and tidal conditions, among others. In addition, the results displayed are from the unique dredged material placements (e.g., Phase 1 and 2) and not a combination, thus results from the first nourishment may not represent true values of total primary consolidation. Fig. 7. Comparison of approximate dredged fill thickness to elevation gain from field survey data at two years post-nourishment. 4. Discussion 4.1. Dredged Material Geomorphology The application of BUDM is a viable solution for wetlands that are unable to keep pace with sea level rise (Wigand et al. 2017; Raposa et al. 2020; Davis et al. 2022), and these marsh restorations require a three-dimensional understanding of geomorphic processes (Ganju 2019). A Spring and Fall 2020 BUDM nourishment on the northern half of Sturgeon Island placed 15,291 m3 of dredged material and resulted in an elevation increase of 0.19 ± 0.11-m two years after placement. A common construction method for BUDM is via hydraulic placement of dredged material (VanZomeren et al. 2019; Davis et al. 2022, Raposa et al. 2023). The hydraulic placement of dredged material will typically result in mounding of coarser-grained sediment in the immediate discharge area while finer-grained sediment remains in suspension longer and travels farther (Whitbeck et al. 2019, Piercy et al. 2023). In addition, scour paths can form as the flow of dredged material finds a path of least resistance away from the discharge point, resulting in a non-uniform surface elevation (Piercy et al. 2023). These topographic features were identified following the initial placement of dredged material on Sturgeon Island (Figure 4). However, with the additional material deposited during Phase 2, the surface appeared flatter as dredged material filled in the surface depressions, and tides and currents reworked material to more stable configurations, though this may be a by-product of the surveying and interpolation methods. The placement of dredged material can impact wetlands by filling in tidal creeks, which are essential to the function of the ecosystem as they facilitate tidal exchange. These can be temporary impacts though as greater thicknesses of dredged material deposits within these depressions (Tyler and Bailey 2019, Harris et al 2021), resulting in higher magnitudes of consolidation and reformation of the depressional features, so the tidal creek network will reestablish itself within the new tidal prism. These observed geomorphic changes at Sturgeon Island emphasize that the collection of high-resolution topographic data prior to BUDM placement is immensely beneficial. In particular, the initial topography can dictate how dredged material will disperse and how to best contain the placed material, whether through natural means (berms, vegetation) or artificial (coir logs, haybales, earthen containment dikes) options. For BUDM projects, the dynamic nature of geomorphic processes during and immediately (<30 days) post placement should be expected and viewed positively. For example, prior to material placement, the northern portion of Sturgeon Island sloped downward from the west to the east and was drained by a tidal creek network off the eastern edge. This creek made it especially difficult to retain dredged material in the subaerial zones. Despite containment attempts, the dredged slurry consistently scoured a "path of least resistance" underneath or around the containment, which resulted in the nourishment of a mudflat off the eastern edge of the island. In most settings, it is unrealistic to expect complete containment of material. Even with an expansive network of earthen containment dikes and drainage weirs, the complete containment of dredged material is monetarily impossible, and containment methods and their construction may negatively impact the marsh surface and vegetation recovery (Berkowitz and VanZomeren 2020; Harris et al. 2021). The "loss" (i.e., dispersal) of dredged material from the direct placement area should be anticipated and built into construction monitoring and adaptive management plans when possible to provide maximum benefits as this can develop a more natural and stable transition (i.e., slope) between marsh types. In addition, sediment leaving the immediate system is likely to stay within the region. For example, Perkey et al. (2024a) investigated the fate of dredged material leaving a BUDM nourishment via a network of tidal creeks and found that the majority of material remained within adjacent mudflats. Ultimately, these enhanced mudflats can provide a multitude of ecosystem services from expanded feeding grounds to wave reduction (Perkey et al. 2024a). Furthermore, recent flume erosion tests conducted on fine-grained sediment, similar to what was placed on Sturgeon Island, revealed that shear stresses in excess of 0.4 Pa were required for erosion 30-days after deposition (Perkey et al. 2024b). Under typical tidal conditions, an area adjacent to Sturgeon Island experienced approximate current induced shear stresses of only 0.05 Pa (Perkey et al. 2024b) meaning the sediment was relatively immobile one month after placement. However, it is possible to exceed these shear stresses during storm induced currents (Perkey et al. 2024a). 4.2. Wetland Surface Elevation Wetland elevation is a primary driver of vegetation species distribution and health as it impacts hydroperiod (i.e., inundation) (Kirwan et al. 2012; Morris et al. 2013; VanZomeren and Piercy 2020; Jafari et al. 2024). The final elevation of a BUDM nourished wetland depends not only on surface flow and topographic controls but also on consolidation. Specifically, self-weight consolidation of the dredged material and underlying wetland foundation consolidation will produce a reduction of elevation with time (Graham and Mendelssohn 2013; VanZomeren and Piercy 2020; Harris et al. 2020; Harris et al. 2024). For more traditional civil infrastructure settings, the determination of consolidation is often a routine exercise that yields highly accurate results. However, wetland environments are highly dynamic systems driven by biological, geological, geotechnical, and meteorological forcings (French 2006; Fagherazzi et al. 2012; Harris et al. 2024), and the presence of these uncertainties produces a large range of potential outcomes during hydraulic placement of dredged material in these environments (Harris and Shawler 2024). Thus, the collection of actual project outcomes (i.e., elevation with time, soil conditions, etc.) is essential to improving predictability of future BUDM marsh enhancement projects. The repeat elevation data collected following the two BUDM placements on Sturgeon Island indicated that most of the consolidation occurred within the ~6 months following placement. This is consistent with consolidation theory. This follows a logarithmic trend where a majority of elevation loss occurs soon after loading is completed and will decrease as time progresses (Terzaghi et al. 1996). Thus, if planning for elevation surveys post BUDM placement, it is recommended to allocate more resources immediately following completion and lengthening the time in-between surveys as time progresses. A comparison of dredged fill thicknesses to overall elevation gain revealed that approximately two -thirds of gained elevation was lost due to consolidation two years post placement (Figure 7). This consolidation is a combination of self-weight consolidation of dredged material and primary consolidation of the wetland foundation. By contrast, Graham and Mendelssohn (2013) found that 23-cm of increased elevation following a BUDM nourishment in a brackish marsh in Louisiana was largely erased 2.5-years later. The difference between our observations at Sturgeon Island and that of Graham and Mendelssohn (2013) may be attributable to differing subsurface conditions, as Sturgeon Island has a shallow sand foundation (<2-m below ground surface) that is significantly less compressible than typical deltaic soils found within Louisiana. 4.3. Consolidation Modeling Consolidation models can be calibrated to field data using two main adjustments: (1) the time-rate or (2) the magnitude of consolidation (Jafari et al. 2019). These parameters are determined by soil conditions within the field, but it is notoriously difficult to determine via laboratory experiments because these environments commonly comprise easily disturbed, highly heterogeneous deltaic soils (Jafari et al. 2018, 2019ab). Previous studies have shown that laboratory-determined soil permeability can be an order of magnitude off from in-situ conditions (Tavenas et al. 1983; Jafari et al. 2019), which is largely attributed to sample disturbance within the collection, transportation, and extrusion processes that alter the soil matrix. However, laboratory compressibility data can depict in-situ conditions (Jafari et al. 2018). To account for this uncertainty, a range of potential surface elevations were modeled within the investigation. In general, the field data fell within the faster and more compressible ends of this range (Figure 6). The calibrated foundation Cc values were significantly higher than the base values of 0.25 and 0.55, determined via standard oedometer testing (Table 3 in the appendix). This was not unexpected as these Cc values were uncharacteristically low when compared to Cc values determined via w(%) correlations (1.7 to 2.1). An interpolated foundation surface (original wetland surface) was determined during the July 2021 field survey (Figure 2a). Although the core spacings were coarse (between 5-to 20-meters), it was estimated that the foundation had consolidated an average of 0.09-m across the transects with values at some locations experiencing as much as 0.17-m of consolidation. This indicates the need to account for foundation consolidation when considering a BUDM nourishment, especially in wetlands with high fines contents and low bulk densities (0.09 g/cm3) (Mudd et al. 2009; Davis et al. 2022). In addition, the in-situ measurement of soil conditions using cone penetrometer tests is strongly recommended over laboratory methods to limit the effects of sample disturbance. 4.4. Trade-Offs of Thick vs Thin BUDM At Sturgeon Island, dredged material thicknesses averaged 39-cm and 37-cm immediately following Phase 1 and Phase 2, respectively. This resulted in an average gain of 19-cm of elevation capital within the immediate placement area following consolidation (Figure 6). At other recent wetland BUDM nourishments within southern New Jersey, dredged material thicknesses averaged 15-cm at Ring Island in 2014, 17-cm at Fortescue in 2016, and 31-cm at on the marsh platform and >62-cm within the pools at Avalon in 2014 and 2015 (two phases) (NJDEP 2023). Outside of the region, 15-cm to 23-cm of sediment was deposited at a marsh within the Brunswick estuary in Georgia (Mohan et al. 2021), while a marsh at Fowl River in Alabama received 15-to 61-cm within the direct placement area and <31-cm in the surrounding fringe areas (Nester and Rees 1988). It was unclear how long after placement these thicknesses were determined, but most likely the wetland foundation and dredged material would have consolidated to some degree. Regardless, the thicknesses of dredged material deposited at Sturgeon Island, and most full-scale BUDM sites, far exceed natural rates of accretion in the New Jersey region (2–6 mm/yr; Haaf et al., 2021) and storm sedimentation (1–12 cm) observed on many US Atlantic and Gulf Coast marshes (Nyman et al., 1995; Goodbred and Hine, 1995l Cahoon et al., 1995; Baustian and Mendelssohn, 2015; FitzGerald et al., 2020; Moore et al., 2021). This implies that, excepting large-scale overwash events on barrier coasts (e.g. Williams, 2012), burial by BUDM may cause greater disturbance to the marsh systems than natural processes, but the benefits from the gain in elevation capital likely warrant this disturbance. To understand the long-term effects these BUDM nourishments could have on Sturgeon Island, low, moderate, and high SLR scenarios for pre-and post-BUDM at one location (Site C, Figure 6) in CWEM were evaluated (Figure 8). At the time of this work, CWEM did not have the capability to simulate a disturbance, so these model results represent a fully established wetland at two distinct elevations. The model results indicate S. alterniflora aboveground biomass quantities are initially reduced for the post placement conditions because wetland vegetation, specifically S. alterniflora, exhibits a parabolic distribution with elevation (Morris et al. 2013; Davis et al. 2017) and the pre-nourishment elevation at Site C was at a more optimal inundation frequency. After 30 years, the model predicts that the post-nourishment site will have the same aboveground biomass. Furthermore, the model predicts that the post-nourishment scenario will maintain elevation capital relative to sea level for an additional 18-28 years, reducing the percent inundation. Past a certain threshold, the addition of dredged sediment will smother existing vegetation and convert the surface to bare earth (Davis et al. 2012; Harris et al. 2024), as seen at Sturgeon (Figure 5). This concept is consistent with the optimal range for storm-driven overwash deposition on fringing marshes of 5-10 cm, beyond which plant mortality is likely (Walters and Kirwan, 2016). BUDM applications provide immediate elevation gains to at-risk wetlands, but vegetation in nourished wetlands will not resemble control sites for many years post placement and may never approximate reference conditions (Raposa et al. 2023). As for design thicknesses, Raposa et al. (2023) found that vegetation can recolonize faster following a thinner lift of dredged fill when compared to a thicker one, but differences had disappeared at the three-year mark. Similarly, NJDEP (2023) found that vegetation within thinner BUDM colonized faster during the first few years. Douglas et al. (2021) showed vegetation recovery within 3 years at a BUDM with an average thickness of 17-cm. This is consistent with the findings from Sturgeon Island as significant revegetation did not occur until 3 years had passed following an average post-consolidation elevation increase of 19-cm (Figure 5). However, at Ring Island, NJ, vegetation was slow to recover for six years post-placement despite an average placement thickness of only 15-cm (NJDEP 2023). Taken together, these studies indicate that dredged material thickness is not the only factor determining plant recovery. Instead, long-term plant recovery is dictated primarily by the final marsh surface elevation (NJDEP 2023). Thus, because significant site disturbance is likely the outcome regardless of material thickness, it is more advantageous to place a thicker dredged material layer so long as sufficient material is available and the target marsh platform elevation will not be exceeded following consolidation. 5. CONCLUSIONS This study focused on monitoring two BUDM-placements on Sturgeon Island, NJ. Specifically, this investigation determined how much elevation gain occurred, how much loss was attributable to consolidation of the dredged material and foundation soils, how quickly the site revegetated, and the long-term impacts of wetland nourishment when considering projected sea level rise rates. The study concludes that: • With each BUDM placement, the surface elevation became flatter as dredged material filled in the depressions and was reworked by tides and currents. Two-thirds of the immediate elevation gain was eventually lost two-years post nourishment due to a combination of the dredged fill and wetland foundation consolidation and sediment reworking. • The surface was elevated on average by 0.19 ± 0.11-m two years after construction. This value accounts for dredged material placement with subsequent consolidation of soil and sediment reworking. • Consolidation models generally agreed with field conditions, but a reliance upon laboratory derived inputs can contribute large degrees of uncertainty. Instead, modelers should opt for in-situ based methods such as cone penetrometer tests to assess soil conditions. • Post BUDM, the site was converted to bare earth as the vegetation was smothered by the dredged material. More than half of the nourished area revegetated after three growing cycles post-BUDM, with continued revegetation of the low marsh occurring annually. • Sufficient material was placed from the two BUDM placements to offset SLR by 18-28 years based on CWEM results. For any wetland BUDM project, the determination of target elevations is the most critical aspect of the project. This is contingent on the type of soil being dredged, volume of dredged material available, foundation characteristics, topography, tidal ranges, target vegetation species and densities, and the desired project lifespan. It is essential for coastal managers to have realistic expectations for BUDM nourishments and anticipate the site to be bare earth for multiple growing cycles. If sufficient elevations are achieved relative to tidal inundation, more rapid revegetation will occur at lower elevations for S. alterniflora colonization and more slowly at higher elevations with lower tidal inundation frequencies. However, these higher elevations confer increased long-term resiliency of the site that will balance out the short-term disturbance to vegetation. Managers must be prepared for the loss of dredged material from the direct placement area so they can maximize benefits to adjacent areas that may receive material. Lastly, they must account for the consolidation of not only the dredged material, but also the foundation to determine dredged material placement heights and resultant surface elevations. 6. Data availability Data and results will be made available upon request to corresponding author. 7. Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper. Financial support was provided by the US Army Corps of Engineers Dredging Operations and Environmental Research (DOER) Program and the Department of Defense (DoD) through its Science, Mathematics, and Research for Transformation (SMART) program. 8. Acknowledgements The authors would like to acknowledge the U.S. Army Corps of Engineers through the Dredging Operations and Environmental Research (DOER) Program and the Department of Defense (DoD) through its Science, Mathematics, and Research for Transformation (SMART) program. In addition, the USACE Philadelphia District, Joint Airborne Lidar Bathymetry Technical Center of Expertise (JALBTCX), The Wetland Institute, and the University of Washington NSF Natural Hazards Engineering Research Infrastructure (NHERI) Reconnaissance Experimental Facility (RAPID) for support of field data collection. Lastly, the authors would like to thank Mary Bryan Barksdale for noting an error within Figure 7 of the original paper that led to clarifying correction to this article. 9. Appendix Table 1. Timeline of survey types and equipment on Sturgeon Island, NJ utilized in this study and dredged material placement windows. Date Survey Type Equipment March 5, 2020 RTK-Survey Trimble RTX March 16-19, 2020 Phase 1: Spring Dredged Material Nourishment May 3, 2020 Terrestrial LiDAR Leica ScanStation P40 July 3 2020 Airborne LiDAR CZMIIL September 2, 2020 RTK-Profiles Trimble RTK September 9-26, 2020 Phase 2: Fall Dredged Material Nourishment UAS-Orthomosaic DJI Mavic 3 November 18, 2020 RTK-Profiles Trimble receivers/antennas/data collectors April 1, 2021 RTK-Profiles Trimble SX10 scanning total station and Trimble RTK July 14, 2021 UAS-LiDAR miniRANGER-UAV April 20, 2022 RTK-Profiles Trimble RTX August 23, 2021 UAS-Orthomosaic DJI Mavic 3 August 23, 2022 UAS-Orthomosaic DJI Mavic 3 September 27 – October 31, 2022 Phase 3: Fall Dredged Material Nourishment September 9, 2023 UAS-Orthomosaic DJI Mavic 3 *September 14, 2020 occurred during the second dredged material placement. Table 2. Summary of initial elevation, depth to sand layer, and thickness of dredged material additions for Phases 1 and 2 for the modeled Sturgeon Island sites. Site Initial Elevation NAVD88 (m) Depth to Sand Layer (m) Approximate Dredged Material Addition (m) Phase 1 Spring 2020 Phase 2 Fall 2020 A 0.67 1.49 0.15 0.20 B 0.26 1.62 0.51 0.46 C 0.42 1.65 0.37 0.44 D 0.21 1.19 0.52 0.38 Table 3. Coastal Wetland Equilibrium Model inputs for Site C prior to dredged material placement for a moderate sea level rise scenario. Physical Inputs Biological Inputs Name Input Reference Name Input Reference Sea Level Forecast 100 cm Coleman et al. 2022 Kopp et al., Lower Growth Limits 9.2 cm, MSL Davis et al. 2022 Run Time 100 yrs Upper Growth Limits 90.0 cm, MSL TWI data Sea Level at Start -12.2 cm, NAVD88 NOAA Station Optimum Elevation 50 cm, MSL Average, Morris et al. 2002, 2013 Starting SLR 0.42 cm/yr NOAA Station Maximum Root Depth 32 cm Harris et al. 2024 Mean Tidal Amplitude 60.1 cm NOAA Station Peak Aboveground Biomass 1450 g/m2 Zhang et al. 2016; Davis et al. 2022 Marsh Elevation @ t0 (MSL) Site Specific Below Ground Bio to Shoot Ratio 2.0 g/g Morris and Sunberg (2024) Suspended Min. Sed. Conc. 23 mg/l Coleman et al., Below Ground Turnover Rate 0.5/yr Morris and Sunberg (2024) Suspended Org. Conc. 0 mg/l Time to Maturity 1 yrs Morris and Sunberg (2024) Initial Accretion Rate 4.2 mm/yr McGauley (2024) OM decay rate 0.5 g/g/yr Morris and Sunberg (2024) Flood Frequency 704/yr Morris and Sunberg (2024) Sediment LOI above-marsh 20 % of dry weight Morris and Sunberg (2024) Minimum Mineral Input 2.0 mg/cm2/yr Morris and Sunberg (2024) Table 4. Summary of base and calibrated values for relevant permeability (kv) and compression index (Cc) values. Site Dredged Material Foundation Base kv Calibrated kv Base Cc Calibrated Cc A - kv/1 0.55 1.65 B - kv/5 0.55 1.55 C - kv/10 0.55 2.00 D - kv/8 0.55 0.85