Seasonal Variability of Forces Controlling Sedimentation in the Sundarbans National Forest, Bangladesh

Southwest Bangladesh, located on the Ganges-Brahmaputra-Meghna delta, is experiencing the impacts of sea level rise (SLR) due to processes at both the local and global scale. In particular, regional alterations of the hydrodynamic network due to embankment construction has drastically altered effective SLR, placing millions of inhabitants at risk of prolonged inundation, and threatening the world’s largest continuous mangrove stand, the Sundarbans National Forest (SNF). In order to effectively employ landscape recovery solutions, an understanding of local sediment transport and deposition is critical. This field-based study investigates the sediment dynamics between the mangrove platform and tidal channels of the SNF using data from a variety of instruments and sediment samples collected within a forested sub-basin (~20 km2) fed by a major tidal channel. We observe profound seasonal variability within the sub-basin, with the wet season exhibiting deeper and longer inundation of the mangrove platform and greater suspended sediment concentrations (SSC). Further, there exists a trend of decreasing SSC and median grain size from the perimeter of the SNF to the interior, and decreasing SSC from the tidal channel to the platform at both locations. We project seasonal platform sedimentation rates ranging from 0.17 ±0.16 cm in the dry season to 1.8 ± 0.35 cm in the wet season. Importantly, the annual deposition rate measured at either location is sufficiently rapid to keep pace with observed rates of effective SLR published in other studies (~1.0-1.7 cm/yr). Based on our results, it appears that many controls on sedimentation are both covariant and of similar importance to land aggradation in the SNF. While inundation depth and frequency will likely increase under future SLR scenarios, sediment supply is threatened by India’s proposed River Linking Project, which could decrease the sediment loads of the Ganges and Brahmaputra Rivers by as much as 75% and 25%, respectively. These rivers provide the sediment for the entire delta, and we predict that with decreasing SSC, some regions—particularly interior sediment depleted regions—may begin to deteriorate and become submerged, including within the SNF.


INTRODUCTION
Mangrove forests provide a variety of societal benefits, including storm-surge buffering (e.g., Badola and Hussain, 2005), economic stimulus from ecotourism (e.g., Satyanarayana et al., 2012), land stability and erosion minimization (e.g., Carlton, 1974;Saenger and Siddiqi, 1993), and carbon sequestration (e.g., Mcleod et al., 2011;Ray et al., 2011;Alongi, 2012;Twilley et al., 2018). As populations around the world move toward coastlines and either alter landscapes or eradicate mangroves altogether, it is imperative that we understand how these systems behave under pseudo-natural conditions to help protect these valuable resources. This study will focus on a portion of the Sundarbans National Forest (SNF), a relatively pristine mangrove forest located in southwestern Bangladesh (Figure 1). The total extent of the Sundarbans reaches across political boundaries into India, expressing a continuum of fluvial and tidal forces with increased distance from the modern Ganges-Brahmaputra-Meghna River mouth (Allison and Kepple, 2001;Rogers et al., 2013;Flood et al., 2015Flood et al., , 2018. Immediately to the north of the SNF is an inhabited region where the construction of earthen embankments ("polders") has resulted in a substantial modification of the tidal prism, leading to tidal amplification and extensive network reorganization (e.g., Pethick and Orford, 2013;Wilson et al., 2017). While attributed to anthropogenic modification, this tidal amplification also impacts the SNF, andcoupled with eustatic sea-level change and compaction-has resulted in an effective sea level rise of 0.7 cm/year at the southern coast, increasing to 1.7 cm/year near the city of Khulna, which enhances flood risk across the entire tidal deltaplain (Figure 1; Steckler et al., 2010;Pethick and Orford, 2013). Furthermore, the National River Linking Project (NRLP) proposed by India is expected to reduce freshwater discharge of the Ganges River by 24%, with a corresponding reduction in sediment load by 39-75% if fully implemented (Higgins et al., 2018). Those authors further predict a reduction in delta-wide aggradation from 0.36 to 0.25 cm/year under this scenario. Given that >60% of the sediment deposited in the SNF is recently derived from the mainstem Ganges-Brahmaputra-Meghna (Allison and Kepple, 2001;Rogers et al., 2013), we expect the impacts of this decrease in sediment supply to be particularly acute in the SNF.
Several recent investigations have investigated the health and stability of the SNF platform elevation in the face of locally accelerated relative sea level rise (e.g., Rogers et al., 2013;Auerbach et al., 2015;Bomer et al., in review). These previous efforts have addressed the timing of seasonal sedimentation, and estimated deposition rates of 1.0 ± 0.9, 1.1 ± 1.0, and 1.32-2.16 ± 0.17-0.26 cm/year, respectively. Collectively, these studies demonstrate that at present, there is sufficient sediment being delivered to the SNF for the vertical accretion to keep pace with rates of local sea level rise, and that the platform is maintaining its elevation relative to mean high water. However, knowledge gaps remain concerning what controls sediment delivery and accumulation within the SNF, and how this varies spatially. To address this need, this study investigates sediment delivery and deposition within a sub-basin in the SNF, and how we might expect similar regions to respond given changes to these functions. Specifically, we combine data from acoustic and optical instrumentation with physical sediment samples to address seasonal changes in: (i) tidal channel suspended sediment concentrations (SSC), water velocity, and grain size; and (ii) platform inundation depth, duration, and grain size. These specific parameters have been identified as especially important in controlling wetland sedimentation in zero-dimensional, massbalance models of marsh growth and mangrove sediment dynamics (e.g., Krone, 1987;Allen, 1990;French, 1993;Furukawa and Wolanski, 1996;Furukawa et al., 1997;Temmerman et al., 2004;Capo et al., 2006;Horstman et al., 2015).

STUDY AREA
The SNF is the largest continuous mangrove stand in the world, with a total area of ∼10,000 km 2 comprised of several major tidal channels separated by vegetated sub-basins (Figure 1). This project focuses on a ∼20-km 2 sub-basin in north-central SNF, immediately adjacent to the human-modified agricultural region, ∼100 km inland from the Bay of Bengal (Figure 1). Water and sediment are primarily delivered to this region by a tidal channel extending from the Sutarkhali River, called the Sorbathkhali, which lies between two primary conduit channels that connect this region to the broader Ganges-Brahmaputra-Meghna river system: the Pussur and Shibsa Rivers (Figure 1). The majority of freshwater and sediment are delivered to the region by inland advection from the Bay of Bengal during the southwest monsoon season (July to September; e.g., Barua, 1990;Barua et al., 1994;Allison and Kepple, 2001;Rogers et al., 2013;Shaha and Cho, 2016;Hale et al., 2019). As a result, the entire SNF is characterized by fresh water (∼0 PSU) from July to late November, and brackish to saline waters (∼20-30 PSU) for the remainder of the year (Shaha and Cho, 2016;Ayers et al., 2017). Accompanying the fresh water is a substantial change in local suspended sediment concentration (SSC), which varies by a factor of five or more from dry-season minima (∼0.2 g/l) in June to Monsoon maxima (>1.5 g/l) in September (Hale et al., 2019). These previous studies provide an excellent backdrop for the research presented here, which focuses on the specific processes responsible for delivering sediment from the major tidal channels onto the mangrove forest platform within the SNF.
The specific focus of this study is on two areas within a representative ∼20-km 2 sub-basin in the SNF, supplied with sediment and water by a single channel, the Sorbothkhali (Figure 1). The Sorbathkhali is ∼90 m wide at its diversion from the Sutarkhali River, with typical thalweg depths ranging from 7 to 10 m. The channel narrows and shoals with distance away from this junction, to a width of ∼15 m and depth of ∼0 m at most interior reaches during spring low tides. While the absolute platform elevation is difficult to constrain owing to seasonal variability and limited benchmark availability, the average inundation depth has previously been reported as ∼0.3 m (Auerbach et al., 2015). Assuming that the Sorbathkhali channel is the conduit for the majority of flooding water (Auerbach et al., 2015), approximately 6 × 10 6 m 3 of water is conveyed during spring tides (e.g., Auerbach et al., 2015;Bomer et al., in review The primary vegetation here is the endangered mangrove Heritiera fomes, locally known as the Sundri, which can grow to a maximum height of 25 m and prefers freshwater. Other species found in this region, presented using their local names, include Gewa, Bain, Pussur, and Kakra (Ghosh et al., 2016).

MATERIALS AND METHODS
This study benefits from the use of a variety of datasets. Longterm seasonal changes were measured using a Schlumberger conductivity, temperature, depth (CTD) sensor continuously deployed in a channel adjacent to the SNF at the Sutarkhali Forest Ranger Station (SFRS) from October 2017 to March 2018 (Figure 1). In addition, suspended sediment concentration was recorded ∼30 cm above bed (cmab) using a Fisher Scientific optical backscatter sensor 3 + (OBS). The OBS was calibrated to SSC using a method described by Hale et al. (2019). In short, water samples were collected at the same time and depth as OBS readings, prior to the long-term instrument deployment. Known sample volumes were then filtered using 0.4 mm glass-fiber filters to determine SSC. A linear relationship between measured OBS response and filtered SSC was then derived. Within the SNF, discrete observations were made during individual spring tide cycles during both the wet and dry seasons at two locations within the sub-basin: an exterior site approximately 1 km from the confluence with the Sutarkhali River (Site 1; Figure 1), and an interior site approximately 5 km further into the forest (Site 2; Figure 1). These locations were selected based on their relative proximity to the sediment source (the Sutarkhali River), and correspond to the locations used in a companion study of platform accumulation (Bomer et al., in review). Wet season measurements were collected in October 2017 and dry season observations took place in March 2018. On each survey day, measurements began as close to low tide as possible, and continued through the entire period of platform inundation. Due to daylight constraints, our results are focused primarily during the rising limb, slack high, and onset of the ebb portions of the tidal cycle, with limited observations at the termination of the falling limb. While the specifics of each survey day varied slightly due to instrumentation and environmental conditions, the general approach included: near-bottom measurements of water depth, velocity profiles (50-250 cmab), temperature, and SSC (5 cmab); water-surface measurements of temperature, conductivity, and SSC, and; measurements of water depth and SSC from the forest platform during periods of inundation (Figure 1). Velocity profiles were measured using an acoustic Doppler current profiler (ADCP), which measured velocity in 10-cm bins extending from 50 to 250 cmab. Measurements were collected at 1 Hz in 512-second bursts every 10 min, with the burst-averaged velocities reported herein. During the dry season, additional water velocity measurements were made further up the bank using a tilt current meter deployed ∼1.25 m below the platform elevation. Platform measurements occurred at distances of ∼50 m and ∼90 m from the creek bank for Sites 1 and 2, respectively, at the locations of the sediment elevation tables used in a concurrent study by Bomer et al. (in review; Figure 1). Efforts were made to ensure that the exact locations of instrument deployment were consistent across seasons. Surface water samples were also collected every 30-60 min to calibrate optical instrumentation measuring SSC, using the method described previously in this section. We used the observed reduction in platform SSC to estimate the deposition (dz) associated with individual tidal cycles, according to the following equation: where SSC1 and SSC2 are the pre-slack-water maximum and post-slack-water minimum SSC (kg/m 3 ), h is the maximum inundation depth (m), and ρ p is the bulk density of the platform (kg/m 3 ). All values used for this calculation are described in Table 1.
In addition to the observations made using acoustic and optical instrumentation, sediment samples were collected from the platform at Sites 1 and 2, and at locations spaced throughout the sub-basin at approximately 400-m intervals in between sites, with repeat sampling in March 2015 and October 2015 (Figure 1). Platform samples used to measure grain size were composed of the seasonally deposited sediment recovered from the sediment tiles deployed as part of a companion study by Bomer et al. (in review). Channel samples (n = 40) were collected at or near low tide, and an effort was made to collect sediment from the mid/upper bank, as well as the channel center (presumptive thalweg) at each station. To the best of our ability, we scraped the upper ∼1 cm of these samples to represent the recently transported material. All samples were measured for their particle size distribution using a Malvern Mastersizer 2000 laser diffraction particle size analyzer, after a ∼0.5 g subsample was immersed in sodium hexametaphosphate and sonified to disaggregate particles. For simplicity, the median grain size (D50) is presented in this study. Finally, sediment dry bulk density was calculated as the quotient of dry sediment weight over volume, using sediment from the upper 10 cm of cores collected at sites 1 and 2.

Multi-Year Observations From Channel Adjacent to SNF
Tidal-channel water depth at Suterkhali Reserve Forestry Station (SRFS) demonstrates a strong semi-diurnal signal, and pronounced spring-neap variability (Figure 2A). Peak water levels were higher in the wet season than the dry by ∼25 cm, while the water level minimum during the dryseason observations was ∼33 cm lower than that of the wet season (Figures 2A,D).
Conductivity during our study period began to increase in early November, reaching a maximum of ∼6 PSU in early January, before apparently dropping to ∼3 PSU until recovery in March ( Figure 2B). An important caveat here is that the CTD was covered in encrusting worms upon recovery. Conductivity during the dry-season field work, as measured by replacement sensors not subjected to long-term deployment, was closer to 15 PSU ( Figure 2B). We suggest that the observed salinity decrease that began in early 2018 is the product of biofouling.
Although the OBS record is incomplete because of deployment challenges in this harsh environment, we can see that SSC at SFRS ranged from 0.05 to ∼1.0 g/l during the wet season field work, and decreased to a minimum of 0.01-0.15 g/l  in early January, at which point the sensor appears to have been completely buried before being re-exposed in March. SSC during the dry season field work ranged from 0.01 to 0.3 g/l ( Figure 2C).

Single-Day Hydrodynamics and Sediment Dynamics Within the SNF
Tidal range remained relatively constant during the days during which discrete sampling took place (Figures 2A,D). As measured at the SRFS, the tidal range during Site 1 sampling days was 4.14 and 4.48 m for wet and dry seasons, respectively (Table 1 and Figure 3). Similarly, Site 2 sampling days exhibited tidal ranges of 4.29 and 4.19 m for the wet and dry seasons, respectively. Despite the similarity of tidal range, the platform inundation depth and duration were drastically different as a result of the overall drop in absolute water levels during the dry season (Figure 2A). Site 1 was inundated for 266 min during the wet season survey day, and to a maximum depth of 0.39 m ( Table 1). The same location was inundated for 122 min during the dry season, to a maximum depth of 0.06 m, which correspond to reductions of 54 and 84%, respectively (Table 1 and Figure 3 columns I, III). At Site 2, the platform was inundated for 321 min, to a maximum depth of 0.43 m during the wet season. During the dry season, on the other hand, the inundation duration reduced 67% to 105 min, and inundation depth reduced 81% to 0.08 m (Table 1 and Figure 3 columns II, IV). On the platform, average SSC was substantially greater during the wet season than the dry season at both locations. At Site 1, we observed an average SSC of 0.27 g/l during the wet season, versus 0.08 g/l during the dry season. Site 2 was similarly dramatic, with average SSC dropping from 0.22 to 0.05 g/l between seasons (Table 1 and Figure 3).
Water velocities were faster during the wet season than the dry season, although not substantially. At Site 1 within the channel, the maximum depth-averaged flooding velocity was 0.64 m/s during the wet season, compared to 0.45 m/s during the dry season, a 30% reduction (Table 1 and Figure 3). Moving further into the sub-basin at Site 2, the depth-averaged velocity reduced ∼47% from 0.44 m/s in the wet season to 0.30 m/s in the dry season (Table 1 and Figure 3). One complication with the burstaveraged velocities presented here is that during ebbing tides at Site 1, the interplay between two channels (Figure 1) sets up ∼30-m-wide eddies, resulting in rapid flow reversals on ∼1minute time scales. Because the burst length is ∼8.5 min, the burst-averaged velocity (Figure 3) is considerably lower than the instantaneous velocity. We do not observe this complication during flooding tides at Site 1, because flow is divergent rather than convergent, preventing eddy formation. At Site 2, the single channel results in uniformly directed currents, with ebb tide burst-averaged velocities that are faster than at Site 1 (Figure 3).
Average near-bed SSC reduced 33% from 0.62 to 0.41 g/l from wet to dry seasons at Site 1 (Table 1 and Figure 3). Unfortunately, near-bed SSC was not measured at Site 2 in the wet season. During conditions of reduced water velocity observed during the dry season, near-bed SSC was ∼0.11 g/l. SSC at the water surface was substantially greater during the wet season than the dry season at both locations. From the OBS, wet season mean SSC was 0.52 and 0.50 g/l at Sites 1 and 2, respectively. During FIGURE 3 | Observations from the single-day experiments, including depth (relative to SRFS), SSC, and water velocity within the tidal channel. In general, we observe a decrease in both flood-oriented current velocity and SSC within the tidal channel from Site 1 to Site 2 (exterior to interior of the sub-basin), and from wet season to dry season. Vertical black line corresponds to the timing of platform inundation, as measured by the platform OBS.
the dry season, the same locations exhibited mean SSC of 0.32 and 0.10 g/l, representing reductions of 38 and 80% (Table 1 and Figure 3). SSC maxima at the channel surface were greatest during the wet season, although the relative decrease during the dry season was not as drastic as for the mean. Max SSC during the wet season were 0.86 g/l (Site 1) and 0.79 g/l (Site 2), compared to dry season SSC of 0.64 and 0.55 g/l at the same locations, for respective reductions of 26 and 30% (Table 1 and Figure 3). In each case, the decrease in SSC within the SNF is smaller than at the longer-term observation station SRFS, where SSC maxima during the wet season field work were ∼0.8 g/l, vs. ∼0.4 g/l dry season (Figure 1).

Sediment Grain Size
Median grain size on the platform decreases from 31 µm at the confluence with the Sutarkhali River to 15 µm near Site 2 (Figure 4) collected in the thalweg vs. the channel flank, nor is there a discernible change between seasons (Figure 4).

DISCUSSION
The field instrument deployment strategy used in this study is admittedly limited in both temporal and spatial coverage due to access issues, safety concerns, and available resources. That said, we are able to provide measurements of sediment in the tidal channels at two locations within this sub-basin, paired with observations on the mangrove platform 50-90 m inland (Figure 1). In doing so, we can address sedimentation conditions in the best-case scenario (wet season at Site 1), and a substantially less-optimal example (dry season at Site 2). Finally, extrapolating from a handful of tidal cycles to a year's worth of deposition will never capture the full range of sedimentary or diagenetic processes, even if our measurements occur during periods corresponding to the annual minima and maxima for conditions related to sedimentation (Figure 2). Instead, as we demonstrate in the following sections, we are encouraged that based on the similarity of our extrapolated deposition rates, and those observed by previous and concurrent studies (Rogers et al., 2013;Bomer et al., in review), we appear to be measuring and comparing the appropriate controlling forces. Future research in this area should build upon these initial observations to develop numerical simulations (sensu Horstman et al., 2015) that will help explain both spatial and temporal variability that this study can not address.

Sedimentation Controls
A persistent challenge in untangling the controls on sedimentation is that they are often covariant. This is clearly the case on the platform in the SNF, as the dry season features lessfrequent and shallower inundation depths, shorter inundation durations, and lower SSC (Table 1 and Figure 3), all of which likely combine to reduce seasonal deposition. One parameter that remains consistent across seasons is the grain size distribution, which controls settling velocity and therefore the efficiency of sediment delivery to and across the platform (Figure 4). Our laboratory analysis of particle sizes, however, ignores the potential for sediment aggregation as a result of either biological processes or the electrochemical properties of fine-grained particles in salt water (e.g., Mietta et al., 2009;Manning et al., 2013). Aggregates (loosely cohesive packets of fine-grained sediment and organic matter) can behave like coarser particles in the water column, requiring additional turbulence to be advected upward, and settling more rapidly during periods of slack water (Manning et al., 2013). We might expect, therefore, that surface SSC in the dry-season channel is reduced more rapidly at slack water than during the wet-season. This effect can be observed most clearly in the time series observations of SSC at Site 1 during the wet and dry seasons (Figure 3). During the monsoon, there is relatively little difference between SSC at the channel surface and bottom, suggesting a well-mixed water column, and relatively little settling by what we assume to be disaggregated sediment particles. While the lowest values of SSC are observed near slack-high tide (Figure 3), there is generally poor correlation between velocity and SSC ( Figure 5).
During the dry season, on the other hand, not only is there a stronger coupling between SSC and velocity ( Figure 5), but while near-bed SSC reduces by ∼60% to ∼0.2 g/l, surface SSC decreases to effectively 0 g/l (Figure 3). Similarly, wet season SSC on the platform at Site 1 is 61% lower than SSC in the channel at the onset of inundation (0.26 g/l vs. 0.67 g/l; Figure 3). During the dry season, the reduction from channel to platform is closer to 90% (0.50 to 0.05 g/l; Figure 3). In each case, the implication is that not only is the overall sediment supply reduced during the dry season, but the sediment in suspension is actually largely unavailable to nourish platform interiors within vegetated sub-basins (here extending ∼6 km from the main Sutarkhali channel). This observation does not discount the presence of biological aggregation, however, we interpret the enhanced dry-season settling as indicative of the relatively greater importance of electrochemical processes here. Based on our observations of environmental conditions, we expect substantially more sediment to be deposited on the platform during the wet season than the dry season. At a fundamental level, we observe that the platform floods more frequently and to greater depth because of the seasonal waterlevel setup (∼25 cm; Figure 2), and the water flooding the platform is carrying dramatically more sediment in suspension (0.22-0.27 vs. 0.05-0.08 g/l; Table 1 and Figures 3, 5). These conditions should promote enhanced deposition, as described by Rogers et al. (2013). In their study of deposition rates at the same locations as this study, Bomer et al. (in review) also observe substantially more frequent and prolonged inundation during the wet season. Those authors measured an average vertical accretion of 1.44 ± 0.35 cm during the wet season (typically measured in October), compared to 0.94 ± 0.27 cm during  the dry season (typically measured in March-May), with no significant difference between the groups.
As described in the section "Materials and Methods, " we estimate the amount of sediment deposited on each tidal cycle, assuming that the difference between the maximum flood-limb SSC and the minimum ebb-limb SSC is the product of deposition. This assumption is supported by a comparison of the settling velocities of the particles at each site, to the depth of the water column at each location. Based on Gibbs et al. (1971), we would expect the median grain size particles at Site 1 (31 µm) and Site 2 (14 µm) to settle at rates of 0.085 and 0.023 cm/s, respectively. These values are likely conservative, as mud particles often aggregate or flocculate as described above, settling at faster rates than the constituent particles (Mietta et al., 2009). Based on these conservative settling rates and the inundation depths described in Table 1, we would expect that during the wet season, the water column at Site 1 could be evacuated of sediment in as few as 7 min during slack high water, while Site 2 might require as long as 30 min of placid conditions. According to our velocity measurements in the tidal channel -a much more dynamic environment -near-slack conditions (U < 0.10 m/s) persist for ∼40 min at each location, offering ample time for sediment to settle from suspension (Figure 3). Previous research has demonstrated the effectiveness of pneumatophores at trapping sediment, with the mangrove platform often serving as a oneway pump for fine-grained material (Wolanski, 1995;Furukawa et al., 1997;Victor et al., 2004). As such, conditions promoting sediment accumulation are unlikely to improve with increased distance from the channel. Our observations should therefore be considered a "best-case" scenario. Locations further from the sediment source are unlikely to demonstrate deposition rates faster than those estimated here, given the well-understood process of reduced SSC with distance from source both from the literature (e.g., Delaune et al., 1978;Reed, 1988;Christiansen et al., 2000), and between Site 1 and Site 2 in this study (Figure 3).
To estimate deposition on the platform during a single tidal cycle, we can begin with an examination of the observations made at Site 1 during the wet season, where we observe a maximum SSC of 0.39 g/l during the flood limb, and a minimum SSC of 0.14 g/l on the ebb. Interestingly, this reduction in SSC 0.4 0.6 0.8 SSC (g/1) (0.25 g/l) is almost identical to the difference back-calculated by Auerbach et al. (2015), when considering sediment deposition in the adjacent poldered region. Assuming that the SSC we measured near the platform bed was consistent throughout the ∼0.39 m water column, 97.5 g/m 2 of sediment would come out of suspension at slack high water. Dry bulk density for surface samples collected at Site 1 is 1100 kg/m 3 in the wet season, and the resulting quotient suggests deposition of 8.7 × 10 −3 cm per inundation cycle ( Table 1). While this value may appear unimpressive, the significance is gained when we extrapolate based on the number of inundation periods that occur under similar conditions of SSC and depth. If we assume, for example, that this location is influenced by similar peak monsoon forces for 90 days (180 inundation periods) per year, we would expect deposition of ∼1.6 cm during that time frame (Table 1); similar to the average deposition of 1.1 cm/wet season observed by Rogers et al. (2013). The same location we would expect an order of magnitude lower deposition (0.3 cm) based on the conditions observed during the dry season ( Table 1). Predicted deposition rates at Site 2 are similar to Site 1, ranging from 2.1 to 0.06 cm/season for the wet and dry seasons, respectively ( Table 1). Of course, these estimates assume that the entire decrease in SSC at slack high tide is a result of deposition, that there is no subsequent erosion, and that this accounts for all of the deposition that occurs during a given inundation cycle. These assumptions are justified based on our calculations of settling times (<10 min), our observations of little to no sediment resuspension on the platform (Figure 3), and the observations of Furukawa et al. (1997) who describe sediment trapping by pneumatophores. We concede here, however, that these backof-the-envelope calculations only apply directly to the interior of the vegetated platform proximal to (50-90 m) the channel margin, as it is readily accepted that sediment accumulation decreases -often exponentially -with distance from the channel edge (Delaune et al., 1978;Reed, 1988;Christiansen et al., 2000). For the wet season, at least, our average predicted deposition (1.8 cm; Table 2) compares favorably with the measurements made by Bomer et al. (in review; 1.4 cm; Table 2). Our dry-season prediction (0.3 cm) on the other hand, is much lower than was observed by Bomer et al. (in review; 0.9 cm; Table 2). Again, our observations are from individual days extrapolated across an entire season, so a one-for-one comparison is likely not reasonable. Furthermore, Bomer et al.'s (in review) measurements of seasonal sediment deposition were typically made in September-October (wet season) and March-May (dry season); periods of time that were sensibly based on the Ganges-Brahmaputra river discharge. As we can see in this study, however, monsoon-like conditions (nearly fresh water, elevated SSC, higher water levels) persist in the tidal channels well into November (Figure 2), allowing for the potential of additional deposition that would be attributed to the dry season in that study. Additionally, in this analysis we are assuming that conditions remain similar throughout either the monsoon or dry season, which we know to be a simplification (Figure 2). We are nevertheless encouraged by the similarities between our calculation and measurements from previous research, and look forward to future colleagues extending these results with a fullscale numerical simulation, the likes of which was outside the scope of the present research.

Comparison Between Locations
In addition to the profound seasonal changes observed at our study locations, there are also important distinctions from within our sub-basin that help us better understand how sediment dynamics evolve within Ganges Brahmaputra tidal sub-basins as a whole. One factor is platform elevation, which is not only different between Sites 1 and 2, but the magnitude of this difference appears to change across seasons. To assess this variability, we adjust the inundation depth by the difference in the absolute magnitude of tide height as measured at SFRS (Figure 5). During the wet season, the platform elevation at Site 2 is approximately 4 cm lower than Site 1. During the dry season, this difference extends to 15 cm ( Figure 5). This calculation relies on the assumptions that the instrument deployment locations are identical across seasons, as is the modification of the tidal waveform between SRFS and the study sites. Given the similarities between the relative inundation depths at each site across seasons, we consider these assumptions to be valid. We suggest that sediment texture plays an important role in controlling platform elevation. In short, the finer grains present at Site 2 (15 vs. 30 µm; Figure 4) are more susceptible to compaction from dewatering, which would be particularly important during the dry season when inundation frequency and duration are reduced (e.g., Knott et al., 1987;Nuttle and Hemand, 1988;Nuttle et al., 1990;Bomer et al., in review). This is further supported based on our observations of sediment dry bulk density at the platform surface, which is ∼35% greater at Site 1 than Site 2 (means of 935 and 685 kg/m 3 , respectively, Table 1), and which also increases by ∼10% from the wet to dry seasons at Site 2 (876 and 968 kg/m 3 , respectively), suggesting that compaction is taking place. Another factor driving the sediment dynamics is proximity to large conduit channels, particularly during the dry season when SSC delta-wide are lower (Barua, 1990;Hale et al., 2019). As described above, we observe relatively little difference between the in-channel SSC at Sites 1 and 2 during the wet season (Figure 3). In contrast, there is a dramatic difference between the two sites during the dry season, with substantially less SSC Frontiers in Earth Science | www.frontiersin.org 9 August 2019 | Volume 7 | Article 211 at the interior location (Site 2). The reason for this may be similar to the difference in the observed reduction of SSC from channel to platform described in section "Sedimentation Controls" above, where an overall decrease in material availability, and increased settling velocities due to salinity-induced flocculation combine to reduce the travel distance for a given particle. While this may seem relatively intuitive, this provides the first empirical evidence for this phenomenon within the SNF. This has very important implications as tidal river management (TRM) projects are being implemented in the poldered region with the goal of combating some of the deleterious side effects of the embankments. With TRM, periodic breaches of the polder wall reconnect sub-basin land with the tidal channel to accelerate land accretion (Khadim et al., 2013;Hossain et al., 2015). Based on our observations, TRM effectiveness would change dramatically based not only on the season, but also the distance from the larger tidal channels. We expect this concept to be explored in greater detail in the near future.
Beyond the obvious elevation difference observed in the sub-basin between Sites 1 and 2, and their relative distances from the Sutarkhali channel, they also exhibit fundamentally different behavior between channel hydrodynamics and sediment transport. In neither location is velocity a strong predictor of SSC ( Figure 5). Indeed, the strongest correlation occurs at Site 2 during the dry season, where velocity explains ∼35% of the in-channel SSC variance. This weak relationship might suggest that at each site, local resuspension plays a less important role in controlling in-channel SSC than advection from up-or down-stream sources. While this may be true, the roles of local resuspension and advection are likely much more nuanced.
A comparison of the Shields parameter (θ t ) and roughness Reynolds number (Re * ) reveals that with the exception of slack high water, the critical shear stress for resuspension is always exceeded ( Table 1; Shields, 1936;Miller et al., 1977). Here, shear velocity was determined calculated as the quotient of von Kármán's constant divided by the slope of the curve of velocity versus the natural log of depth, as described by the Law of the Wall. We performed these calculations using data from the relatively quiescent dry season; the stronger velocities observed during the wet season would only serve to enhance shear in the water column and stress on the bed. Of course, evaluating the Shields parameter in fine-grained systems can be fraught, however, our primary reason for using it here is to highlight an interesting difference in the timing of SSC maxima observed at our two field locations. At Site 1, the maximum flood-limb SSC occurs immediately prior to platform inundation during both the wet and dry season (Figure 5). This region is characterized by broad (∼15 m-wide), muddy banks, which can continuously provide sediment for resuspension as they inundate, assuming sufficiently strong shear velocities. In contrast, peak flood-limb SSC at Site 2 is observed near the lowest water levels, long before any material in suspension could be delivered to the platform. The channel has narrowed substantially at this distance inland, and is completely drained during spring low tides during the wet and dry seasons. The banks here are not as wide as at Site 1 and are heavily vegetated, which can baffle near-bed flow as the water elevation rises, thereby preventing additional erosion (Furukawa and Wolanski, 1996). Figure 3 compares the average current speed from the upward-looking ADCP (located on the channel bottom), with the velocity on the channel flank, ∼1.25 m below the platform elevation, during the dry season. While the flood velocities in the channel were similar at the two locations, the velocities on the bank were dramatically different. At Site 2, the maximum flooding velocity on the bank was 0.03 m/s, a full order of magnitude lower than the fastest channel velocity, and unlikely to resuspend bank sediment. In contrast, Site 1 demonstrated bank velocities as fast as 0.15 m/s, which is 34% as fast as the channel velocity, and certainly capable of resuspension (Figure 3). This additional source of sediment is one of several mechanisms that help to explain the increased deposition seen at Site 1 relative to Site 2, as observed by Bomer et al. (in review; Table 1).

Spatio-Temporal Variability and Future Implications
Interestingly, the deposition rates determined by this and other recent studies (e.g., Rogers et al., 2013;Auerbach et al., 2015;Bomer et al., in review) are similar to the local deposition rates observed across the entire Holocene (0.7-1.0 cm/year; Ayers et al., 2016). Further to the west, where sedimentation has been reduced as a result of the eastward migration of the Ganges River over the past 4000 years, the longterm deposition rates are closer to 0.1-0.2 cm/year (Flood et al., 2015). It would be an oversimplification to imply that this spatial variability stems exclusively from a moving sediment supply. Many boundary conditions, including upstream sediment supply, reservoir construction, sea level, and polder construction have all played a role in effecting delta-wide sedimentation. Nevertheless, the area that has been subject to the most-dramatic decrease in sediment supply (western delta plain) also demonstrates a reduced accumulation rate. This provides important context for the effects of a reduced sediment supply to the SNF-which we can augment with our observed differences between wet season and dry season SSC and sedimentation. For example, consider the possible impacts of a reduced sediment discharge associated with India's planned dams and diversions through the NRLP on sedimentation in the SNF: if fully implemented, one result would be a reduction in the annual suspended sediment load of the Ganges and Brahmaputra rivers by 39-75% and 9-25%, respectively, with the majority of this reduction occurring during the wet season (Higgins et al., 2018). Collectively, this amounts to a reduction in the total suspended sediment load by 14-18%. Rogers et al. (2013) estimate that approximately two thirds of the sediment deposited each monsoon season is of recent (<6 month) Ganges-Brahmaputra origin, with the remainder being sourced from antecedent shelf and floodplain deposits. Using the more optimistic estimate of sediment load reduction (14%), average wet-season deposition is reduced from 1.8 to ∼1.6 cm ( Table 2). Using instead the high-end estimate for total sediment load reduction, wet-season deposition decreases further to ∼1.4 cm ( Table 2). In both scenarios, dry season deposition is <0.2 cm ( Table 2). Of course, by reducing the freshwater discharge of these rivers by 24 and 6% respectively, we can also expect the effective duration of monsoon-like conditions to be reduced by a month or more (Higgins et al., 2018). Indeed, if we rework our initial estimate of deposition to reflect a substantially shorter monsoon season (i.e., from 3 to 2 months, or 180 platform-inundating events to 120), wetseason deposition drops to 1.2 cm. With less sediment being delivered to the region, a smaller freshwater contribution to the wet season water levels, and a longer duration of conditions that do not favor sedimentation, we expect sedimentation rates within the SNF to be reduced as a result of a fully implemented IRLP. Further, we expect the combined effects of these impacts to result in a greater reduction in sedimentation than any factor acting in isolation. Using the combined worstcase scenario of maximum suspended load reduction and a shortened wet season, our estimate of deposition drops to <1 cm ( Table 2). Without adequate sedimentation to offset the regional relative and effective rates of SLR, the SNF may drown in place as predicted by Loucks et al. (2010). Importantly, our analysis neglects the possibility of sediment supplied by channel scouring downstream of the dams associated with the IRLP, similar to the scenario on the Mississippi River described by Nittrouer and Viparelli (2014). That said, as we discussed in the previous sections, the majority of the platform is being supplied with silt and finer-grained sediment, as coarser grains are incapable of being transported to the platform interior. As such, this remains an important area of future research.

CONCLUSION
Sedimentation in the SNF is controlled by a combination of independent factors, including maximum water depth, sediment particle size, inundation duration, and SSC. Understanding the interaction of each of these parameters, as well as their relative importance, is critical to understanding the health and sustainability of the SNF and tidal deltaplain as a whole. Furthermore, the SNF stands out as a natural case study to contrast with the adjacent, heavily modified poldered region, where significant change is necessary to ensure continued habitability. In this study, we observe two overarching trends from a tidal sub-basin located in the SNF. First, we see that environmental conditions during the wet season are substantially more conducive to sediment deposition, with morefrequent and deeper inundation of the platform by fresh water enriched with suspended sediment. The apparent importance of salinity on controlling the formation of aggregates was a particularly interesting find with delta-wide implications in the face of future changes to the fresh water supply. Second, we see that the location proximal to the larger tidal channel (Site 1) exhibits greater SSC values, coarser grain size, and sits at a higher elevation than the location deeper in the forest (Site 2), regardless of season. Despite this, predicted monsoon deposition is actually greater at Site 2 than Site 1, largely because of the lower sediment bulk density (Table 1). This combination is particularly interesting, as it suggests that the coarser grains at Site 1 are building a more-stable platform capable of maintaining a higher elevation along the perimeter of SNF sub-basins, while the finer-grained interior is subject to more-intense compaction and dewatering. We must be cautious, therefore, in extrapolating our measurements from channel margins to span the entire lower delta, as the inaccessible interior regions may be more vulnerable to changes in sediment supply. We should also consider these results in practical applications like TRM, and ensure that sites selected for sluice gates have a sediment supply that is both robust and sufficiently coarse to build lasting platforms. Perhaps more important than the observations of the modern conditions are the implications of these observations in the context of scenarios for reduced sediment supply in response to upstream river damming and diversions. In the worst case, where sediment delivery is reduced by 14-18%, and the duration of fresh water conditions is shortened by one month or more, we predict deposition rates ∼50% lower than are currently observed. Not only will this diminish the effectiveness of land reclamation strategies such as TRM, these rates are below those of local effective sea level rise, resulting in a threat to the continued viability of the SNF. Elevated SSC and fresh water conditions not only appear critical to the delivery of sediment from the tidal channels to the mangrove platform, but are also most likely to be affected by changes to the fresh water or sediment discharges of the larger Ganges and Brahmaputra rivers.

AUTHOR CONTRIBUTIONS
RH was responsible for study design and acquiring funding. RH, CW, and EB conducted the field work. RH drafted the initial report, which was subsequently revised by all authors, each of who contributed to the intellectual content.