Feedback between Basin Morphology and Sediment Transport at Tidal Inlets: Implications for Channel Shoaling

<jats:p>Increasing societal pressures (e.g., population growth and urbanization) are driving land use change practices in coastal areas that could potentially alter the hydrodynamics and sediment transport patterns near coastal inlets in ways that might exacerbate existing shoaling conditions. To investigate the potential impact of coastal development, a numerical model is used to predict the long-term evolution of an idealized lagoonal-type barrier island inlet under five different morphological conditions that transitioned from net sediment import to net sediment export. The simulations were designed to address the potential effect of inter-tidal placement and land reclamation on sediment transport and the resulting deposition/erosion patterns. Estuaries that were deeper and devoid of extensive tidal flats tended to promote sediment import and had a greater propensity to exacerbate channel shoaling. Simulations that were characteristic of inter-tidal placement showed net export, yet the likelihood of channel shoaling was increased because some of the material eroded from the tidal flats was deposited in the deeper channels as opposed to being carried out the inlet throat. Alternatively, it was found that regions in which the intertidal area was restricted to elevations higher in the tidal frame, which also showed a net export, produced greater sediment loss in the inter-tidal zone that tended to bypass the deeper sections, reducing the likelihood of channel shoaling.</jats:p>

section_content: Increasing societal pressures (e. g., population growth and urbanization) are driving land use change practices in coastal areas. In response to the greater societal presence, engineering activities, such as armoring, sediment diversions, damming, and reclamation, are altering sediment supply, thus affecting the total volume and distribution of sediments in coastal bays and estuaries [1, 2]. These modifications can alter flow pathways, increase or decrease fetch and average water depth, and change tidal propagation characteristics. The long-term effect of anthropogenic modifications to coastal inlets in a regime of global sea level change, which itself can affect sediment distribution patterns, is not well known. Urbanization and changes in land use practices could have unforeseen consequences to sediment transport and channel shoaling patterns, thereby exacerbating existing conditions through increased dredging requirements and associated costs. 

Numerical modeling can be used to explore the dynamic feedback mechanisms between the hydrodynamics and the morphology. Long-term simulations can identify fundamental shifts in the large-scale sediment transport patterns that arise from changes in the sediment supply. The initial sediment distribution (bathymetry) can be manipulated to explore the role of sediment placement to control the dominant sediment transport flux (import vs. export). Configurations that favor sediment import and increase the likelihood of channel infilling can be categorized in terms of morphology (and to some extent hypsometry), sediment supply, or new sediment placement. Developing a generalized framework applicable to a broad class of inlets could help to assess the potential impact to navigation from large-scale sediment redistributions (e. g., dredge disposal, beneficial placement, land reclamation, and armoring). 

Previous work aimed at understanding tidal dynamics at inlets has focused on tidal wave propagation, which is controlled by friction, ocean tidal characteristics, bay configuration, and the vertical distribution of land cover [3, 4]. Non-linear tidal distortion has been used to explain the behavior of many natural inlets and is a primary metric as to whether a system is an importer or exporter of sediment [3, 5]. Flood dominated systems have stronger currents during the flood portion of the tidal cycle and tend to import sediment, while ebb-dominated systems have stronger currents during ebbs and tend to export sediment. 

Flood or ebb dominance is dictated by the tidal wave propagation characteristics, which are partially determined by the vertical distribution of land mass in the bay, or hypsometry (e. g., [6] [7] [8] [9] ). Open bays with fewer inter-tidal areas tend to be flood dominated, whereas bays with extensive inter-tidal areas tend to be ebb dominated (e. g., [5] ). A system that either imports or exports sediment will also experience a net gain or loss of sediment. This can alter the vertical distribution of land mass in the bay, producing feedback such that an open bay may develop more inter-tidal flats and begin to shift towards ebb dominance. Alternatively, a bay with extensive inter-tidal flats may deepen and experience a net increase in areas with open water and shift towards flood dominance. Bays without substantial inter-tidal and upper marsh habitat may be particularly vulnerable to conversion to open water as relative sea level rises [10]. Inlets that are only weakly flood or ebb dominated may be more susceptible to the effects of sea level change and coastal engineering practices that alter sediment availability and distribution, as small changes in hypsometry could produce a net shift towards import or export. In such cases, less impetus is required to alter the dominant sediment transport direction and could alter the rates of channel infilling. 

The hydrodynamic processes that are responsible for flood or ebb dominance have been studied extensively in the absence of sediment transport and morphological change (e. g., [5, 11, 12] ). By neglecting the feedback between morphology and hydrodynamics, there is no way to determine if a flood (ebb) dominated system is likely to maintain its present configuration or develop an alternative tendency towards ebb (flood) dominance in the future. Long-term morphodynamic modeling can be used to explore the feedback between tidal hydrodynamics and morphology to identify the key processes that trigger a shift in the net sediment transport direction (i. e., import versus export). 

Two-dimensional numerical process-based models have been used successfully to investigate long-term hydrodynamic and sediment transport trends in a variety of examples. While such studies have included realistic basin shapes and nearshore bathymetry [13], many examples have used idealized bathymetry and schematized rectangular basins, including tides only [14 -16], tides plus waves [17], and varying timeframes ranging from 50 to 8000 years [14, 15, 18]. Styles et al. [19] investigated the long-term evolution of tidal inlets using idealized models of nine inlets, representing conditions for the U. S. Atlantic, Pacific, and Gulf coasts. That study demonstrated that the Coastal Modeling System (CMS) was able to correctly predict the large-scale morphological features associated with tidal inlets, i. e., the development and growth of the ebb and flood shoals, the emergence of intertidal flats, the formation of a narrow and deep entrance channel, and the formation and progression of an extensive creek network. This report extends the work of Styles et al. [19] through a series of long-term model simulations to investigate the role of hypsometry in regulating sediment import and export. 

The objective of this study is to investigate the morphodynamic characteristics that determine whether an inlet imports or exports sediment. Hypsometry is used as a singular metric that reduces the complex bathymetric variability associated with real inlets into a simple and efficient modeling framework that can be consistently applied to a broad class of inlets with similar characteristics. In particular, the effect of inter-tidal placement on channel shoaling is investigated through the judicious arrangement of sediment to match specified hypsometry. 

The approach is to utilize the CMS modeling system to conduct 200-year simulations of an idealized barrier island tidal inlet using 5 different initial bathymetric configurations that span a range of hypsometries from an open basin to a system with extensive intertidal flats. It is reasonably assumed that all inlets exist between these two extremes, so that the simulations cover the full range of possible lagoonal type bar-built barrier island systems. This also ensures that the suite of model simulations cover both sediment import and export. The methods section describes the CMS modeling system. The results focus on two major aspects of the model simulations: (1) the hydrodynamics and (2) the morphology. The discussion ties sediment transport patterns to the different hypsometries and how these relationships may impact the carrying capacity of the system, the role of hypsometry in regulating sediment yields and associated changes in channel depth that could affect navigation.
section_name: Introduction

section_content: The methodology used in this study is intended to demonstrate the impact of hypsometry on the century scale morphology response of a lagoonal-type barrier island system. Twenty years of tide and wave forcing were applied to idealized bathymetry in CMS with a morphologic acceleration factor of 10 to simulate 200 years of morphologic evolution for each case. Test cases varied bay morphology between a system which was 79% intertidal flats and a system with featureless relatively deep bathymetry. The morphology factor of 10 was confirmed to be acceptable for this study, as differences introduced were small compared to the other model simplifications. Morphologic and hydrodynamic model output were analyzed to investigate the inlet and bay characteristics related to the transition from net sediment import to export. This approach is adapted from the century scale simulations of Styles et al. [19], but expanded to include an idealized approach to varying bay hypsometry between simulations. While the model simulations were idealized, the boundary conditions, forcing, basin configuration, and hypsometry are representative of lagoonal type barrier island tidal inlets.
section_name: Methods

section_content: Simulated domains and starting bathymetry were idealized based on the simplifications of the characteristics of the Humboldt Inlet following Styles et al. [19]. While the domain characteristics were defined to be representative of Humboldt Inlet, these features were idealized to permit comparisons across a broader range of bay systems, particularly lagoonal type barrier island systems. The domain was simplified to a rectangular inlet, bay, and nearshore oriented with the alongshore and along-bay axes on the Y-axis and the ocean on the left side of the domain. The along-inlet axis was oriented on the X-axis. Inlet and bay geometry was selected by matching the approximate horizontal aspect ratio and area of the real Humboldt Inlet system. The inlet dimensions were selected as approximately 2000 m by 550 m (Figure 1 ). The initial depth in the inlet was set as a constant value of 11. 9 m, representing the average of considered bathymetry within the real inlet. While the real-world Humboldt Inlet is jettied, this characteristic was only partially represented in the numerical model domain. Cells along the sides of the inlet are defined as non-erodible, capturing the shore protection and inlet stabilization function of the jetties. The idealized nature of the approach and relatively coarse resolution was determined to make the jetty impact on ebb flow dynamics and littoral drift inadequately resolved. The bay dimensions were selected as approximately 22 km by 3 km (Figure 1 ). A reference depth of 3. 4 m in the bay was determined from the average of considered bathymetry. Depths in the bay were varied to represent sediment distribution and availability and are described in Section 2. 2. The offshore extent of the domain was established following Styles et al. [19] with the intention of minimizing the interaction between the ebb shoal and the offshore boundary. The dimensions of the ocean section of the domain were selected as 6 km by 22 km (Figure 1 ). Bathymetry outside of the bay and inlet were specified by an equilibrium beach profile [20], using the same profile as Styles et al. [19] (Figure 1 ). A single grain size of 0. 20 mm was selected to match the sediment in the real inlet, following Styles et al. [19]. Sediment transport was calculated using the van Rijn [21] transport algorithm for noncohesive sediments. A common model domain was applied in both CMS Flow and CMS Wave. The model domain was composed of 18,396 active cells, ranging from 50 m by 50 m in the inlet to 200 m by 200 m at the offshore edge. The X-coordinate grid refinement continued throughout the bay at the same cell spacing. et al. [19] with the intention of minimizing the interaction between the ebb shoal and the offshore boundary. The dimensions of the ocean section of the domain were selected as 6 km by 22 km (Figure 1 ). Bathymetry outside of the bay and inlet were specified by an equilibrium beach profile [20], using the same profile as Styles et al. [19] (Figure 1 ). A single grain size of 0. 20 mm was selected to match the sediment in the real inlet, following Styles et al. [19]. Sediment transport was calculated using the van Rijn [21]
section_name: CMS Domain

section_content: Bay bathymetry was varied between simulations to represent the morphologic response to different hypsometric curves. Initial bathymetry within the bay was specified to be along-bay (Y, Figure 1 ) uniform and symmetric about the bay centerline (Xbc, Figure 1 ). A variety of starting bathymetries was analyzed to investigate the impact of sediment distribution within the bay on morphodynamic patterns.
section_name: Initial Bay Bathymetry

section_content: In each case, bay depths were specified to fit a curve defined by the hypsometric curve equation presented in Strahler [6] :

where: y = non-dimensional bed elevation that can range from 0 to 1;
section_name: Systematically Varying Hypsometric Curves between Cases

section_content: Bay bathymetry was varied between simulations to represent the morphologic response to different hypsometric curves. Initial bathymetry within the bay was specified to be along-bay (Y, Figure 1 ) uniform and symmetric about the bay centerline (X bc, Figure 1 ). A variety of starting bathymetries was analyzed to investigate the impact of sediment distribution within the bay on morphodynamic patterns.
section_name: Initial Bay Bathymetry

section_content: In each case, bay depths were specified to fit a curve defined by the hypsometric curve equation presented in Strahler [6] :

where: y = non-dimensional bed elevation that can range from 0 to 1; R = non-dimensional bay area that can range from 0 to 1; r = non-dimensional constant that can range from 0. 01 to 0. 5; z = non-dimensional constant exponent that can range from 0. 03 to 2. 

Bathymetry variations among cases altered bay storage volume and sediment availability. The selected hypsometric curves were determined as variations between a case with large inter-tidal areas and a case with no intertidal area. The extensively inter-tidal case was established with approximately half of the bay area submerged at mean tide. A hypsometric curve was selected by sampling viable parameter space and manually choosing between parameters which satisfy the half-submerged at mean tide condition. The bathymetry formed by the hypsometric curve defined as r = 0. 02 and z = 0. 05 was selected to represent a bay with a large inter-tidal area and is referred to as Hypsometry Case 1 (HC1, Figure 2 ). J. Mar. Sci. Eng. 2022, 10, 442 5 of 32 R = non-dimensional bay area that can range from 0 to 1; r = non-dimensional constant that can range from 0. 01 to 0. 5; z = non-dimensional constant exponent that can range from 0. 03 to 2. 

Bathymetry variations among cases altered bay storage volume and sediment availability. The selected hypsometric curves were determined as variations between a case with large inter-tidal areas and a case with no intertidal area. The extensively inter-tidal case was established with approximately half of the bay area submerged at mean tide. A hypsometric curve was selected by sampling viable parameter space and manually choosing between parameters which satisfy the half-submerged at mean tide condition. The bathymetry formed by the hypsometric curve defined as r = 0. 02 and z = 0. 05 was selected to represent a bay with a large inter-tidal area and is referred to as Hypsometry Case 1 (HC1, Figure 2 ). Additional starting bathymetries were selected as variations between HC1 and a uniform bathymetry (HC4; r = 10 -5, z = 0. 99999; Figure 2 ). Two intermediate cases were established as Hypsometry Case 2 (HC2; r = 0. 0167, z = 0. 2083; Figure 2 ) and Hypsometry Case 3 (HC3; r = 0. 0133, z = 0. 3667; Figure 2 ). Areas above mean low water (MLW) and water storage vary over a wide range between HC1 and HC4 (Figure 2, Table 1 ). An additional case was evaluated in which the inverse of the non-dimensional hypsometric curve defining HC1 was shifted to match the same condition of equal submerged and emergent areas at mean tide (IHC1, Figure 2 ). Additional starting bathymetries were selected as variations between HC1 and a uniform bathymetry (HC4; r = 10 -5, z = 0. 99999; Figure 2 ). Two intermediate cases were established as Hypsometry Case 2 (HC2; r = 0. 0167, z = 0. 2083; Figure 2 ) and Hypsometry Case 3 (HC3; r = 0. 0133, z = 0. 3667; Figure 2 ). Areas above mean low water (MLW) and water storage vary over a wide range between HC1 and HC4 (Figure 2, Table 1 ). An additional case was evaluated in which the inverse of the non-dimensional hypsometric curve defining HC1 was shifted to match the same condition of equal submerged and emergent areas at mean tide (IHC1, Figure 2 ). The selected hypsometric curves were made dimensional with the bay half-width (W bh, Figure 1 ) and reflected across the bay centerline (X bc, Figure 1 ). An along-bay channel was added to the domain to maintain simulation stability (Figure 3 ). Maximum bay depths were constant in all starting bathymetries. Several changes were determined to be necessary to maintain simulation stability. The width of the channel along the center of the bay was iteratively increased for simulation stability. The width of the along-bay channel was set to 400 m because simulations with along-bay channel widths of 200 m became unstable. A limiting of the slope of 1:20 was also imposed on the across-bay elevation profile to minimize avalanching and resolve issues with the release of water trapped in inter-tidal areas. A wetting/drying depth of 0. 2 m was also established to resolve further issues with the release of water trapped in inter-tidal areas. Modified initial bathymetries (Figures 3 and 4 ) are expected to address the questions that guided initial hypsometric curve selection, but are not described by the same curves. The selected hypsometric curves were made dimensional with the bay half-width (Wbh, Figure 1 ) and reflected across the bay centerline (Xbc, Figure 1 ). An along-bay channel was added to the domain to maintain simulation stability (Figure 3 ). Maximum bay depths were constant in all starting bathymetries. Several changes were determined to be necessary to maintain simulation stability. The width of the channel along the center of the bay was iteratively increased for simulation stability. The width of the along-bay channel was set to 400 m because simulations with along-bay channel widths of 200 m became unstable. A limiting of the slope of 1:20 was also imposed on the across-bay elevation profile to minimize avalanching and resolve issues with the release of water trapped in inter-tidal areas. A wetting/drying depth of 0. 2 m was also established to resolve further issues with the release of water trapped in inter-tidal areas. Modified initial bathymetries (Figures 3 and 4 ) are expected to address the questions that guided initial hypsometric curve selection, but are not described by the same curves. Each case represents a different basin configuration between extensive intertidal flats and a relatively deep uniform, featureless bathymetry (Figure 4 ). HC1 has extensive intertidal flats around a central channel, HC2 has many shallow subtidal areas, HC3 is relatively deep, HC4 is relatively deep and featureless, and IHC1 has many intertidal areas that are relatively high in the tidal frame. HC1 and IHC1 begin with similar intertidal storage volumes, but an additional 33% of the bay is intertidal in HC1 and an additional 48 M m 3 of sediment begins the simulation in the bay (Table 1 ). HC2, HC3, and HC4 all begin with little to no intertidal area, resulting in similar intertidal storage volumes, but HC2 begins with 47 M m 3 more sediment than HC3 and 98 M m 3 more sediment than HC4 (Table 1 ). The differences between HC1, HC2, HC3, and HC4 represent a transition from large features, many of which are intertidal, to no features. Sediment availability and features on which sediment can deposit on or erode from were captured in the range of simulations. IHC1 has several similarities to HC1, but while HC1 represents large-scale wetland reclamation with areas at a large range of intertidal elevations, IHC1 instead represents land reclamation in which all intertidal area is high in the tidal frame. 

storage volumes, but an additional 33% of the bay is intertidal in HC1 and an additional 48 M m 3 of sediment begins the simulation in the bay (Table 1 ). HC2, HC3, and HC4 all begin with little to no intertidal area, resulting in similar intertidal storage volumes, but HC2 begins with 47 M m 3 more sediment than HC3 and 98 M m 3 more sediment than HC4 (Table 1 ). The differences between HC1, HC2, HC3, and HC4 represent a transition from large features, many of which are intertidal, to no features. Sediment availability and features on which sediment can deposit on or erode from were captured in the range of simulations. IHC1 has several similarities to HC1, but while HC1 represents large-scale wetland reclamation with areas at a large range of intertidal elevations, IHC1 instead represents land reclamation in which all intertidal area is high in the tidal frame.
section_name: Systematically Varying Hypsometric Curves between Cases

section_content: The ocean boundaries of the domain were forced with 13 tidal constituents from Humboldt (Table 2 ) to match Styles et al. [19]. Tidal constituent amplitude and phase were extracted from the National Oceanic and Atmospheric Administration (NOAA) tide database [22]. Mean tidal range was 1. 49 m, with predominantly semi-diurnal components (primarily M2, S2 and N2), but large diurnal (primarily K1, O1, and P1) components as well. The thirteen selected constituents did not capture all tidal variation over the 20 years of tidal forcing at the physical site, but differences were considered to have a smaller impact than most of the simplifying assumptions. Selected tidal constituents capture multiyear tidal range modulations, with maximum boundary condition water levels of approximately 1. 5 m (Figure 5 ).
section_name: Tidal Forcing

section_content: The ocean boundaries of the domain were forced with 13 tidal constituents from Humboldt (Table 2 ) to match Styles et al. [19]. Tidal constituent amplitude and phase were extracted from the National Oceanic and Atmospheric Administration (NOAA) tide database [22]. Mean tidal range was 1. 49 m, with predominantly semi-diurnal components (primarily M2, S2 and N2), but large diurnal (primarily K1, O1, and P1) components as well. The thirteen selected constituents did not capture all tidal variation over the 20 years of tidal forcing at the physical site, but differences were considered to have a smaller impact than most of the simplifying assumptions. Selected tidal constituents capture multi-year tidal range modulations, with maximum boundary condition water levels of approximately 1. 5 m (Figure 5 ).
section_name: Tidal Forcing

section_content: A representative span of 20 years was selected from WIS station 83,047 to force the long-term inlet model simulation of a simplified representation of Humboldt Bay. Wave spectra characteristics (i. e., significant wave height (H s ), peak period (T p ), and mean wave angle (θ) were determined from the WIS station at a depth of 260 m. The offshore vector was estimated as 298 • clockwise from north. The 20-year span with mean alongshore sediment transport rate closest to the mean over the entire duration was selected to represent a characteristic 20-year period. Alongshore sediment transport rate was estimated with the Coastal Engineering Research Center (CERC) formula as:

where Q is the alongshore sediment transport rate in cubic meters per second, k is a constant (0. 39), γ is the breaker index (0. 78), H br is breaking wave height, and θ br is the breaking wave angle following Rosati et al. [23]. Wave height and angle at breaking were determined by iteratively shoaling the wave with Snell's law and conservation of energy (in a method similar to the description provided by Dean and Dalrymple [20] ) until the depth limited breaking criterion is met. The selected depth limited breaking criterion is described by:

where h br is the depth at breaking. Wave spectra were transformed from the WIS station to depth limited breaking, following a similar procedure to Styles et al. [19]. Wave estimates at breaking were used to calculate CERC alongshore sediment transport rates over the entire record. The CERC alongshore sediment transport rate (ASTR) estimates were sampled with overlapping 20-year averaging windows and one year step size, and distributions of mean alongshore sediment transport estimates and mean alongshore sediment transport magni-tude estimates were determined (Figure 6 ). The 20-year span with the most representative mean alongshore sediment transport rate estimate was determined to begin on 1 January 1990. Selected significant wave heights varied from 0. 205 to 9. 08 m with a mean of 2. 15 m, selected peak periods varied from 2. 64 to 26. 0 s with a mean of 11. 2 s, and selected incident wave angles varied from -32. 7 • to 61. 6 • with an energy weighted mean of 1. 14 • (Figure 7 ). 

where hbr is the depth at breaking. Wave spectra were transformed from the WIS station to depth limited breaking, following a similar procedure to Styles et al. [19]. Wave estimates at breaking were used to calculate CERC alongshore sediment transport rates over the entire record. The CERC alongshore sediment transport rate (ASTR) estimates were sampled with overlapping 20-year averaging windows and one year step size, and distributions of mean alongshore sediment transport estimates and mean alongshore sediment transport magnitude estimates were determined (Figure 6 ). The 20-year span with the most representative mean alongshore sediment transport rate estimate was determined to begin on 1 January 1990. Selected significant wave heights varied from 0. 205 to 9. 08 m with a mean of 2. 15 m, selected peak periods varied from 2. 64 to 26. 0 s with a mean of 11. 2 s, and selected incident wave angles varied from -32. 7° to 61. 6° with an energy weighted mean of 1. 14° (Figure 7 ).
section_name: Wave Record Selection

section_content: Model output from simulations was analyzed with a variety of techniques to compare morphodynamic and hydrodynamic differences associated with the range of the considered starting bathymetries and the factors related to net sediment transport to or from the bay.
section_name: Analysis Techniques

section_content: Morphology was output from the model at intervals of 20 effective days of forcing. This higher frequency data facilitates temporal comparisons, but spatial patterns are only presented as starting and ending conditions to limit extraneous content. Morphology and intertidal area patterns are described in Section 3. 1 to present the full model results that are used to quantify the various differences among the cases in the subsequent sections. Elevation and bathymetry data in the bay are binned into hypsographs in Section 3. 4 to characterize the basic information about changes in elevation distributions among the different cases. Morphology change was calculated by subtracting the morphology at a specific time-step from the initial morphology. This is interpreted as net sediment transport and is integrated over specific regions of the domain in Section 3. 2 to describe sediment and large-scale import and export. The combined morphology information indicates the most important differences among the cases that gain or lose sediment from the bay.
section_name: Morphology and Elevation Data

section_content: Tidal asymmetry was addressed through analysis of the peak velocities and the duration of the flood and ebb cycles in the inlet. The impact of bathymetry on tidal asymmetry was investigated through one-year simulations without wave forcing or sediment transport. Tidal asymmetry associated with bathymetry from 0, 25, 50, 100, and 200 years of simulated morphologic change was investigated for each case, for a total of 25 additional simulations. Water elevation and velocity data were output at 300 s at each point along defined transects in the inlet to improve temporal resolution of tidal signals (Figure 8 ).
section_name: Tidal Asymmetry

section_content: Model output from simulations was analyzed with a variety of techniques to compare morphodynamic and hydrodynamic differences associated with the range of the considered starting bathymetries and the factors related to net sediment transport to or from the bay. ration of the flood and ebb cycles in the inlet. The impact of bathymetry on tidal asymmetry was investigated through one-year simulations without wave forcing or sediment transport. Tidal asymmetry associated with bathymetry from 0, 25, 50, 100, and 200 years of simulated morphologic change was investigated for each case, for a total of 25 additional simulations. Water elevation and velocity data were output at 300 s at each point along defined transects in the inlet to improve temporal resolution of tidal signals (Figure 8 ). A weighted average was applied to along-channel velocity data from the three central transects across the inlet to normalize by water volume. Peak velocities were determined for each ebb and flood cycle from the spatially averaged inlet velocities and compared to evaluate peak velocity asymmetry. Peak velocity asymmetry is presented as the average of the difference between these peak flood and ebb velocities over the year of simulated tidal forcing divided by the average of the non-dominant peak velocities. The durations of ebb and flood cycles were also determined from spatially averaged inlet velocities, with positive velocities corresponding to flood tide and negative velocities corresponding to ebb tide. Ebb and flood durations were compared to evaluate tidal duration asymmetry and are presented as relative differences. Tidal duration asymmetry was calculated as the difference between flood and ebb tide durations divided by the duration of the non-dominant tide (i. e., ebb tide for flood dominant asymmetries and vice versa). Tidal asymmetry informs the interpretation of the morphologic results, as results can generally be grouped between flood-and ebb-dominant asymmetry in the inlet.
section_name: Analysis Techniques

section_name: Results

section_content: Morphology results indicate that the impact of hypsometry on bay morphology change varies between regions of the domain. Differences are particularly evident in the ebb shoal and channel patterns in the bay. These impacts were assessed by comparing bathymetry to morphologic changes after 200 years. 

Morphology changes in several similar ways in all five simulations, indicating the tested range of initial bathymetries only had minor impacts in several regions of the domain. In all cases, the inlet deepens, and relatively deep channels develop through much of the bay. Inlet morphology does vary between simulations, but all cases developed a deep channel along the northern side of the inlet (Y > 0) and a shoal along the southern side (Y < 0). A large, more than 100 M m 3 ebb shoal also forms a small distance offshore of the inlet in all simulations. Ebb shoal volume is compared to volumes of sediment transported to or from other regions of the domain (see Section 3. 2), but specific ebb shoal morphology patterns are not presented in detail. The model also predicts the extensive flattening of the shallower nearshore area within approximately 500 m of the ocean-side shoreline for all tested simulations. Nearshore erosion in tested simulations is not expected to match natural systems, but likely results from limitations in the model setup and the long duration. 

Simulations do predict important differences between initial bathymetry configurations in other regions of the domain. Flood dominant cases (HC2, HC3, and HC4) build extensive flood shoals and deposit sediment in intertidal shoals along the sides of the meandering channels that form in the bay (Table 3 ). Ebb dominant cases (HC1 and IHC1) deposit some sediment across from the inlet (Y ≈ 0), but erode larger volumes of sediment from the intertidal areas (Table 3 ). The bay loses more sediment in simulations where more sediment begins in the bay (Table 3 ), but differences on the order of tens of millions of cubic meters remain in the final volumes of sediment in the bay.
section_name: Morpho-Sedimentary Evolution

section_content: In HC1, morphology change in the bay is dominated by the initial channel deepening (Figure 9 ). The channel connecting the inlet to the along-bay channel also carves through the inter-tidal area near the inlet to shorten the flow path. Accretion is predominantly in an ebb shoal, but some accretion also occurs in flood shoal landward of the inlet. The intertidal area decreases from 81% of the bay to 62% throughout the simulation, with the greatest intertidal area loss near the inlet and adjacent to the along-bay channel (Table 3, Figure 10 ). The morphologic evolution in HC1 simulates the continued prevalence of the singular deep channel, as the channel erodes and intertidal area reduces. The realism of this singular streamlined channel is difficult to assess. Other idealized numerical modeling has suggested that several of the common techniques for incorporating of bed slope effects can adversely impact long time-scale channel behavior predictions specifically leading to over-incision and a lack of branching (e. g., [24] ).
section_name: HC1 Morphology

section_content: In HC2, the shallow subtidal areas adjacent to the main channel along the bay and across from the inlet accumulate sediment (Figure 11 ). The central channel deepens and migrates towards the inlet. Several secondary channels also form off the primary along-bay channel. The intertidal area more than doubles through the course of the simulation with the formation of emergent shoals across from the inlet and on the edges of the channel along the bay (Figure 12 ). The 85 M m 3 more sediment and 71. 8% less of the bay beginning as intertidal area in HC2 are predicted to correspond to shallower, less streamlined channels and more shoal formation throughout the bay than in HC1 (Table 1, Figure 11 ). Many of the depositional features in the bay also become intertidal (Figures 11 and 12 ). The net sediment flux between the bay and the other regions of the domain show that more sediment is eroded from the channels than is deposited in the shallow subtidal areas. 

intertidal area decreases from 81% of the bay to 62% throughout the simulation, with the greatest intertidal area loss near the inlet and adjacent to the along-bay channel (Table 3, Figure 10 ). The morphologic evolution in HC1 simulates the continued prevalence of the singular deep channel, as the channel erodes and intertidal area reduces. The realism of this singular streamlined channel is difficult to assess. Other idealized numerical modeling has suggested that several of the common techniques for incorporating of bed slope effects can adversely impact long time-scale channel behavior predictions specifically leading to over-incision and a lack of branching (e. g., [24] ). intertidal area decreases from 81% of the bay to 62% throughout the simulation, with the greatest intertidal area loss near the inlet and adjacent to the along-bay channel (Table 3, Figure 10 ). The morphologic evolution in HC1 simulates the continued prevalence of the singular deep channel, as the channel erodes and intertidal area reduces. The realism of this singular streamlined channel is difficult to assess. Other idealized numerical modeling has suggested that several of the common techniques for incorporating of bed slope effects can adversely impact long time-scale channel behavior predictions specifically leading to over-incision and a lack of branching (e. g., [24] ). In HC2, the shallow subtidal areas adjacent to the main channel along the bay and across from the inlet accumulate sediment (Figure 11 ). The central channel deepens and migrates towards the inlet. Several secondary channels also form off the primary alongbay channel. The intertidal area more than doubles through the course of the simulation with the formation of emergent shoals across from the inlet and on the edges of the channel along the bay (Figure 12 ). The 85 M m 3 more sediment and 71. 8% less of the bay beginning as intertidal area in HC2 are predicted to correspond to shallower, less streamlined channels and more shoal formation throughout the bay than in HC1 (Table 1, Figure 11 ). Many of the depositional features in the bay also become intertidal (Figures 11 and 12 ). The net sediment flux between the bay and the other regions of the domain show that more sediment is eroded from the channels than is deposited in the shallow subtidal areas.
section_name: HC2 Morphology

section_content: In HC3, channels form and deepen through the relatively deep sub-tidal areas (Figure 13 ). These channels are shallower and more numerous than the channels of HC1 or HC2. Shoals form across from the inlet and along the channels. The pattern of channels and shoals in the bay is also less symmetric about the line X = 0. The intertidal area increases from 2% to 16% of the bay, as many of the shoals that form become intertidal (Figure 14 ). The deeper subtidal areas of HC3, corresponding to 47 M m 3 less sediment beginning in the bay, contribute to 16 M m 3 more sediment being imported into the bay than in HC2 (Table 3 ). Net sediment flux is transported into the bay, so more material is deposited in the shoals than is eroded from the channels.
section_name: HC3 Morphology

section_content: In HC4, channels form and deepen through the initially featureless bay bathymetry (Figure 15 ). Channels in HC4 are wider and less numerous than in HC3 and shallower than the channels of HC1 or HC2. Shoals form across from the inlet and along the channels. Accretion occurs throughout the central two-thirds of the bay (Figure 15 ). The balance of shoal and channel formation result in the net import of 22 M m 3 of sediment (Table 3 ). HC4 imports 6 M m 3 more sediment than HC3, but the initially deeper bay results in 7% less of the bay being intertidal. The intertidal areas that do form are located between the wider channels (Figures 15 and 16 ). These differences from HC3 are related to the additional 51 M m 3 of sediment with which HC3 began in the relatively deep sub-tidal areas (Table 3 ). The morphologic evolution in HC4 simulates the formation of several channels through an initially featureless bay with relatively large intertidal shoals forming between channels (Figure 15 ). In HC4, channels form and deepen through the initially featureless bay bathymetry (Figure 15 ). Channels in HC4 are wider and less numerous than in HC3 and shallower than the channels of HC1 or HC2. Shoals form across from the inlet and along the channels. Accretion occurs throughout the central two-thirds of the bay (Figure 15 ). The balance of shoal and channel formation result in the net import of 22 M m 3 of sediment (Table 3 ). HC4 imports 6 M m 3 more sediment than HC3, but the initially deeper bay results in 7% less of the bay being intertidal. The intertidal areas that do form are located between the wider channels (Figures 15 and 16 ). These differences from HC3 are related to the additional 51 M m 3 of sediment with which HC3 began in the relatively deep sub-tidal areas (Table 3 ). The morphologic evolution in HC4 simulates the formation of several channels through an initially featureless bay with relatively large intertidal shoals forming between channels (Figure 15 ).
section_name: HC4 Morphology

section_content: In IHC1, much of the morphology change in the bay is intertidal erosion, particularly within 1 km of the inlet (Figure 17 ). The intertidal area near the inlet decreases slightly, but the changes only represent 1% of the bay area (Figure 18 ). Material is removed from the intertidal areas in such a way that it is transported out of the bay and only has a minor impact on the fraction of the bay that is intertidal. While all other cases lose material from relatively extensive channels and channel networks, channel features in IHC1 do not extend more than 2 km from the inlet. Approximately 18 M m 3 of sediment is exported from the bay, but the volume of sediment in the bay remains 36 M m 3, 50 M m 3, and 95 M m 3 higher in IHC1 than in HC2, HC3, and HC4, respectively (Table 3 ). The morphologic evolution in IHC1 simulates the erosion of the edges of the elevated intertidal area, particularly near the inlet. In IHC1, much of the morphology change in the bay is intertidal erosion, particularly within 1 km of the inlet (Figure 17 ). The intertidal area near the inlet decreases slightly, but the changes only represent 1% of the bay area (Figure 18 ). Material is removed from the intertidal areas in such a way that it is transported out of the bay and only has a minor impact on the fraction of the bay that is intertidal. While all other cases lose material from relatively extensive channels and channel networks, channel features in IHC1 do not extend more than 2 km from the inlet. Approximately 18 M m 3 of sediment is exported from the bay, but the volume of sediment in the bay remains 36 M m 3, 50 M m 3, and 95 M m 3 higher in IHC1 than in HC2, HC3, and HC4, respectively (Table 3 ). The morphologic evolution in IHC1 simulates the erosion of the edges of the elevated intertidal area, particularly near the inlet.
section_name: IHC1 Morphology

section_content: The domain was divided into four regions to investigate sediment import and export trends in the ebb shoal, the inlet, the bay, and the nearshore (Figure 19 ). The inlet and bay have simple designations as the area between the two sides of the barrier and the area landward of the barrier, respectively. The area indicated as the ebb shoal was selected by hand as the extent of substantial accretion outside of the inlet with a manually selected buffer. The nearshore region was configured to exclude the area designated as the ebb shoal.
section_name: Morphodynamic Change by Region

section_content: Each simulation predicts the formation of a 150+ M m ebb shoal, the removal of 5-10 M m 3 of sediment the inlet, and the loss of 60-100+ M m 3 of sediment from the nearshore (Figure 20 ). Specific volumes vary among cases, but magnitudes and signs are similar. The inlet loses sediment in all cases. The net volume of sediment predicted to leave the inlet appears to be related to tidal prism in these simulations. This aligns with a variety of studies that have found inlet cross-section to be related to tidal prism (i. e., [25] [26] [27] [28] ). After 200 years of morphology change, the tidal prisms of HC2, HC3, and HC4 are similar, HC1 has approximately three quarters of that magnitude, and IHC1 has approximately two thirds of the tidal prism in HC4 (Table 3 ). The volume of sediment that leaves the inlet follows a loosely similar pattern, with HC4 inlet losing the most sediment and the IHC1 inlet the least. All simulations also lose sediment from the nearshore. The transfer of sediment away from the nearshore is likely impacted by some of the limitations of the modeling approach but appears to have some connection to the volume of sediment being imported into or exported from the bay and the volume of ebb shoal. 

Ebb shoals grow in all cases, but the volumes vary. Researchers have related ebb shoal volume to a variety of influences, but the volumes of several ebb shoals have been related to tidal prism with strong influences from the tide and wave energy regime, energy flux through the inlet [29], and the angle between the outflow jet and the shoreline [30]. Similarly, the closely related features that often form at the mouths of rivers (river mouth bars) have been found to be strongly influenced by outflow inertia (e. g., [31, 32] ). Differences between ebb shoal volume predictions do appear to be related to these characteristics, but they are not entirely described by these differences between cases. Ebb shoal volume in HC1, HC2, HC3, and HC4 also appears to be strongly related to volume imported into or exported from the bay, with the sediment supply from HC1 corresponding to a large ebb shoal and the sediment trapping of HC3 and HC4 corresponding to smaller ebb shoals. Large volumes of available sediment have been found to have a strong influence on ebb shoal volume and formation (e. g., [33] [34] [35] ), which these model results also indicate. Differences in sediment availability do not appear to describe the different ebb shoal volume of IHC1. Although the bay loses a similar volume of sediment in HC2 and IHC1, the reduced tidal prism of IHC1 may contribute to the smaller ebb shoal volume, which is consistent with previous findings on the impact of the tidal prism on the ebb shoal volume (i. e., [29, 30] ). 

Sediment import or export from the bay does vary among cases. The net volume of sediment entering and leaving the bay relates to the balance of channel formation and accretion or erosion elsewhere. The ebb dominant cases (HC1 and IHC1) both lose sediment from the bay, but HC1 exports more than twice as much sediment as IHC1 (Figure 20 ). This additional erosion corresponds to the longer channel through the bay and more widespread erosion. The bay does not gain sediment in all flood-dominant cases. The balance between channel formation and sediment deposition results in net export in HC2, but net import in HC3 and HC4. 

In all cases, most regions gain or lose sediment at much faster rates in the first 10 to 25 years of the simulation. Sediment export from the bay in HC2 does not follow the same trend as the other cases. The volume of sediment in the bay remains relatively constant for the first 60 to 70 years and decreases through the rest of the simulation. This transition corresponds to the gradual shift from adding intertidal area to maintaining the same coverage of the bay (Figure 21 ). The rate of intertidal area addition also reduces to near zero over the course of HC3 and HC4, but these changes correspond to less of a change in the volume of sediment in the bay. Sediment export from the bay loosely corresponds to more sediment in the ebb shoal, inlet, and nearshore, but the distribution varies between simulations. The boundary conditions allow sediment to enter and exit the domain, so sediment is not conserved. Ebb grow in all cases, but the volumes vary. Researchers have related shoal volume to a variety of influences, but the volumes of several ebb shoals have related to tidal prism with strong influences from the tide and wave energy regime ergy flux through the inlet [29], and the angle between the outflow jet and the shore corresponds to the gradual shift from adding intertidal area to maintaining the same coverage of the bay (Figure 21 ). The rate of intertidal area addition also reduces to near zero over the course of HC3 and HC4, but these changes correspond to less of a change in the volume of sediment in the bay. Sediment export from the bay loosely corresponds to more sediment in the ebb shoal, inlet, and nearshore, but the distribution varies between simulations. The boundary conditions allow sediment to enter and exit the domain, so sediment is not conserved.
section_name: Net Sediment Transport between Regions of the Domain

section_content: Changes to relative intertidal area differ from changes to sediment volume in the bay, as both shoal formation or erosion and channel formation or deepening impact the sediment budget (Figures 20 and 21 ). The continuous export of sediment from the bay in HC1 only impacts intertidal area abruptly in the first 5 years of simulation and after 75 years. Although IHC1 also loses sediment throughout the simulation, much of this erosion results in lower intertidal area over a relatively constant fraction of the bay. HC2 maintains a similar volume of sediment in the bay as intertidal area increases and loses sediment when the intertidal area stops being added. Net sediment flux into the bay follows a similar trend to intertidal area change in HC3 and HC4, but sediment flux rates into the bay decrease before the rate of intertidal area growth approaches zero. Simulations do not appear to approach a single intertidal area, so initial bathymetry conditions likely impact the achievable relative intertidal area.
section_name: Relative Intertidal Area

section_content: A relationship between intertidal area and ebb-dominant sediment transport has been observed in numerous studies (e. g., [3] [4] [5] [36] [37] [38] ). In these examples, the presence and extent of intertidal flats were related to ebb dominance. The intertidal area is compared to tidal asymmetry in more detail in Section 3. 3, but net sediment transport to or from the bay is also found to have some dependence on intertidal area in the tested simulations. The rate of morphodynamic change was considered between 0, 25, 50, 100, and 200 years and compared with the average intertidal area. Morphodynamic change rates are normalized by the mean tidal flux through the inlet. Intertidal area is normalized by the area of the bay and averaged between the two time-steps represented in the volume change. The volume of sediment imported into the bay generally decreases with increasing intertidal area, with predominantly sediment import where approximately ≤10% of the bay was intertidal (Figure 22 ). The rate of sediment import decreases as cases that import sediment gain intertidal area, with HC3 exporting sediment by the end of the simulation. Cases that export sediment do export less sediment as intertidal decreases. The rates at which sediment is exported from the bay are generally much greater when more than ~60% of the bay is intertidal. The rates at which sediment is imported into the bay are greatest when the bay has no intertidal area and decrease as the bay gains intertidal area until about 10% of the bay is intertidal. 

ized by the mean tidal flux through the inlet. Intertidal area is normalized by the area of the bay and averaged between the two time-steps represented in the volume change. The volume of sediment imported into the bay generally decreases with increasing intertidal area, with predominantly sediment import where approximately ≤10% of the bay was intertidal (Figure 22 ). The rate of sediment import decreases as cases that import sediment gain intertidal area, with HC3 exporting sediment by the end of the simulation. Cases that export sediment do export less sediment as intertidal decreases. The rates at which sediment is exported from the bay are generally much greater when more than ~60% of the bay is intertidal. The rates at which sediment is imported into the bay are greatest when the bay has no intertidal area and decrease as the bay gains intertidal area until about 10% of the bay is intertidal.
section_name: Intertidal Area vs. the Rate of Bay Volume Change

section_content: Tidal asymmetry in the inlet was calculated to compare to net sediment transport direction and erosion/accretion patterns in the bay. Tidal asymmetry was quantified with two parameters: peak velocity asymmetry and tidal duration asymmetry. These quantities were calculated from spatial averages across the inlet in shorter 1-year simulations with 5 min save-point time-steps, but without waves or sediment transport. Simulations were run using the 0-, 25-, 50-, 100-, and 200-year bathymetry for each case. Peak velocity asymmetry indicates how much faster the fastest velocities were in the ebb or flood tide. As an example, a value of 10% represents a scenario in which peak ebb velocities are 10% faster than peak flood velocities. Tidal duration asymmetry indicates how much longer the ebb or flood tide lasted. As an example, a value of 15% corresponds to a scenario in which the flood tide is 15% longer than the ebb tide.
section_name: Tidal Asymmetry

section_content: Peak velocity asymmetries are less than 25% in all cases and, with the exception of the initial conditions for IHC1, do not change direction (Figure 23A ). Both HC1 and IHC1 have ebb dominant peak velocities, but HC1 is marginally more ebb asymmetric. The initial bathymetry in IHC1 makes peak velocities in the inlet flood dominant, but after 25 years, this has reversed and IHC1 remains ebb dominant for the remainder of the simulation. The initially sharp connections between the relatively tall intertidal areas and the channel connecting to the inlet may have some similarities to other findings linking friction between the inlet and bay to flood dominance (e. g., [39, 40] ). After hydrodynamic forcing reorganized a sufficient volume of sediment, IHC1 becomes ebb dominant. Peak ebb velocity asymmetry values range from 7 to 12% during the majority of each simulation. Peak velocities are flood dominant in HC2, HC3, and HC4, but are not clearly ranked until more than 100 years of morphology change have occurred. By the end of the simulation, peak velocities are more flood dominant in HC4 than in HC3 or HC2. Peak velocity asymmetry throughout the simulation matches net sediment transport direction to or from the bay in all cases but HC2.
section_name: Peak Velocity and Tidal Duration Asymmetry

section_content: Tidal asymmetry in the inlet was found to correlate to some degree with the intertidal area after an initial equilibration (Figure 24 ). The relatively angular starting bathymetry of the cases with much more sediment in the bay (i. e., HC1 and IHC1) appears to have had a large impact on tidal asymmetry, which was not reflected in other results. Although tidal asymmetry associated with the initial bathymetry in HC1 and IHC1 was weakly ebb dominant and weakly flood dominant, respectively, the intertidal area decreased and the bay exported sediment quickly. The large change in tidal asymmetry by the 25th year of Tidal duration asymmetries are similar to peak velocity asymmetries but have notable differences (Figure 23B ). Tidal duration asymmetries in HC1 and IHC1 are similar to peak velocity asymmetry, but with larger magnitudes. The faster ebb velocities are related to a larger magnitude increase in the duration of flood tide. Ebb tides begin 10 to 20% longer than flood tides in HC2, HC3, and HC4, but this asymmetry reduces over the course of the simulation. The final morphologies cause 1 to 2% longer flood tides in HC2 and HC3, but 7% longer flood tides in HC4. This results in longer flood tides with faster peak velocities, but HC2 exports sediment from the bay and HC3 and HC4 have little net sediment transport to or from the bay. Net sediment transport must also depend on factors other than tidal asymmetry in the inlet.
section_name: Tidal Asymmetry vs. Intertidal Area

section_content: Tidal asymmetry in the inlet was found to correlate to some degree with the intertidal area after an initial equilibration (Figure 24 ). The relatively angular starting bathymetry of the cases with much more sediment in the bay (i. e., HC1 and IHC1) appears to have had a large impact on tidal asymmetry, which was not reflected in other results. Although tidal asymmetry associated with the initial bathymetry in HC1 and IHC1 was weakly ebb dominant and weakly flood dominant, respectively, the intertidal area decreased and the bay exported sediment quickly. The large change in tidal asymmetry by the 25th year of the simulation indicates that the asymmetry changed with the rapidly changing morphology. After the first 25 years, bathymetry with much intertidal area corresponded to faster and shorter ebb flows through the inlet and cases with less intertidal area had slower and longer ebb flows through the inlet (Figure 24 ). Linear regression skills were evaluated to compare correlations between the intertidal area and the rate of volume change (Figure 22, R 2 = 0. 53), peak velocity asymmetry (Figure 24, R 2 = 0. 85), and tidal duration asymmetry (Figure 24, R 2 = 0. 83). This relationship was not observed in the tidal asymmetry related to the initial bathymetry for HC1 and IHC1. 

the simulation indicates that the asymmetry changed with the rapidly changing morphology. After the first 25 years, bathymetry with much intertidal area corresponded to faster and shorter ebb flows through the inlet and cases with less intertidal area had slower and longer ebb flows through the inlet (Figure 24 ). Linear regression skills were evaluated to compare correlations between the intertidal area and the rate of volume change (Figure 22, R 2 = 0. 53), peak velocity asymmetry (Figure 24, R 2 = 0. 85), and tidal duration asymmetry (Figure 24, R 2 = 0. 83). This relationship was not observed in the tidal asymmetry related to the initial bathymetry for HC1 and IHC1. 

Tidal asymmetry depends on a variety of factors beyond intertidal area, but the results generally match previous findings that shoal formation and tidal flats shorten ebb tides (e. g., [5, 7, [36] [37] [38] ). Peak velocity asymmetry after the first 25 years was found to have the highest correlation with the intertidal area of the considered traits. Peak velocity asymmetry and duration asymmetry are not expected to be entirely independent but were only correlated with an R 2 of 0. 74 across the same data, indicating some independent correlation of each with intertidal area. Tidal asymmetry depends on a variety of factors beyond intertidal area, but the results generally match previous findings that shoal formation and tidal flats shorten ebb tides (e. g., [5, 7, [36] [37] [38] ). Peak velocity asymmetry after the first 25 years was found to have the highest correlation with the intertidal area of the considered traits. Peak velocity asymmetry and duration asymmetry are not expected to be entirely independent but were only correlated with an R 2 of 0. 74 across the same data, indicating some independent correlation of each with intertidal area.
section_name: Tidal Asymmetry vs. Intertidal Area

section_content: Hypsographic curves were developed for each simulation by determining the area above each depth within the bay. Differences between hypsometry at the beginning and end of the simulation describe the changes in bay elevation distribution after 200 years of morphodynamic forcing. Hypsographic curves indicate the distribution of depths throughout the bay, but do not describe horizontal positions. The comparison of initial and final hypsographs demonstrates net gain or loss of area above each elevation. Bay hypsometry changes resulting from the differences in the starting bathymetry follow one of two general patterns. Channels connecting the inlet to the rest of the bay deepen to similar depths in all cases. Flood-dominant cases (HC2, HC3, and HC4) generally deposit sediment in areas outside of channels. Ebb-dominant cases (HC1 and IHC1) generally erode sediment from the intertidal areas. Tidal asymmetry is discussed in greater detail in Section 3. 3. The hypsometry from the beginning and end of each simulation are grouped by tidal asymmetry in the inlet.
section_name: Hypsometry

section_content: Ebb-dominant cases export sediment. Some of the additional sediment in the bay at the start of the simulation is eroded. In HC1, elevations reduced over much of the bay. The intertidal areas lowered, and the channel narrowed and deepened. These morphological changes are demonstrated in Figure 25. In IHC1, elevations within the bay also reduced, but approximately two-thirds of the bay remained at similar elevations (Figure 26 ). Erosion predominantly occurred at intertidal depths and in the area immediately adjacent to the inlet. are compared to (A) relative differences between peak ebb and peak flood velocities, and (B) duration of the tidal cycle that flow enters or leaves the inlet.
section_name: Hypsometry in Ebb-Dominant Simulations

section_content: Ebb-dominant cases export sediment. Some of the additional sediment in the bay at the start of the simulation is eroded. In HC1, elevations reduced over much of the bay. The intertidal areas lowered, and the channel narrowed and deepened. These morphological changes are demonstrated in Figure 25. In IHC1, elevations within the bay also reduced, but approximately two-thirds of the bay remained at similar elevations (Figure 26 ). Erosion predominantly occurred at intertidal depths and in the area immediately adjacent to the inlet.
section_name: Hypsometry in Ebb-Dominant Simulations

section_content: Flood-dominant cases contain large areas of both erosion and accretion within the bay. Erosion in flood-dominant cases is almost exclusively confined to the networks of channels that form. Sediment accretes between these channels. The balance between the channel and shoal formation determines the net sediment transport direction through the inlet. Although elevations in HC2 increase or remain the same across >75% of the bay, channel formation in approximately 12% of the bay drives net sediment export from the bay (Figure 27 ). In HC3, elevations increased or remained the same over approximately 80% of the bay (Figure 28 ), and the intertidal fraction of the bay increased throughout the simulation. The channels that form in HC3 have similar depths to the other cases, but are distributed between multiple somewhat narrower channels, occupying a slightly smaller fraction of the bay. In the Hypsometry Case 4 (HC4) simulation, channel formation and deepening occurred across approximately 12% of the bay, while accretion occurred across ~48% of the bay. The remaining portions showed little change in elevations (Figure 29 ). The volume of sediment in the bay increased more in HC4 than in HC3 over a smaller area, but final elevations were lower and less area was intertidal.
section_name: Hypsometry in Flood-Dominant Cases

section_content: Hypsographs become more similar among cases over the course of each simulation, but important differences remain (Figure 30 ). Similar fractions of the bay are above the 3. 7 m depth contour, but deeper depths vary and shallower depths vary greatly. The channelized flow of HC1 results in deep channels that occupy a similar area to the channels that form in other cases. The lack of a long channel network in IHC1 (Figure 17 ) is also evident in the final hypsography comparisons. Channel dimensions and locations in the bay vary greatly between HC2, HC3, and HC4 (Figures 11, 13 and 15), but the associated
section_name: Hypsometry Convergence and Dependence on Initial Bathymetry

section_content: Hypsographs become more similar among cases over the course of each simulation, but important differences remain (Figure 30 ). Similar fractions of the bay are above the 3. 7 m depth contour, but deeper depths vary and shallower depths vary greatly. The channelized flow of HC1 results in deep channels that occupy a similar area to the channels that form in other cases. The lack of a long channel network in IHC1 (Figure 17 ) is also evident in the final hypsography comparisons. Channel dimensions and locations in the bay vary greatly between HC2, HC3, and HC4 (Figures 11, 13 and 15), but the associated final hypsographs are particularly similar for depths below 3. 7 m, indicating similar elevation distributions. The net import or export of sediment from the bay reduces and stabilizes over time in HC3 and HC4 (Figure 20C ), so hypsometries are not expected to converge noticeably further over a longer time span and variations in the starting bathymetry are indicated to have permanent impacts on the hypsometry. 

final hypsographs are particularly similar for depths below 3. 7 m, indicating similar elevation distributions. The net import or export of sediment from the bay reduces and stabilizes over time in HC3 and HC4 (Figure 20C ), so hypsometries are not expected to converge noticeably further over a longer time span and variations in the starting bathymetry are indicated to have permanent impacts on the hypsometry.
section_name: Hypsometry Convergence and Dependence on Initial Bathymetry

section_name: Discussion

section_content: Initial hypsographic curves labeled HC1-HC4 were constructed to produce morphologies that transition from an estuary with extensive inter-tidal flats (HC1) to an estuary devoid of inter-tidal flats (HC4). These cases were chosen under the assumption that most lagoonal-type bar-built estuaries exist between these two extremes, and that modeling a family of morphological types within this range could offer an insight into a wide class of possible inlet system behavior. The theory implies that tidal asymmetry measured as the difference between the maximum flood/ebb velocities in the inlet throat should shift from ebb dominated (HC1) to flood dominated (HC4) [5]. The model results agree with the theory, as HC1 is ebb dominated, and HC2, HC3, and HC4 progressively become more flood dominated. Correspondingly, the idealized model setup, which is initially unbalanced in terms of total sediment volume within the estuary versus tidal prism, allows for basin-scale adjustments through the redistribution of sediment and corresponding changes in morphology and hypsometry consistent with theoretical predictions. While the final hypsometry, and associated morphology differ, all four cases tend to an equilibrium, in which sediment import or export relaxes asymptotically to a near-zero net gain or loss. 

The majority of basin sediment loss for systems that export sediment is in the subtidal areas and regions that are lower in the tidal frame (HC1). This tends to deepen the basin by the preferential erosion of the sub-tidal areas. This could suggest little change in creek network structure and the size and distribution of tidal flats, while deepening subtidal channels. As the initial bay configuration begins to switch towards sediment import
section_name: Implications of Basin Hypsometry to Navigation

section_content: Initial hypsographic curves labeled HC1-HC4 were constructed to produce morphologies that transition from an estuary with extensive inter-tidal flats (HC1) to an estuary devoid of inter-tidal flats (HC4). These cases were chosen under the assumption that most lagoonal-type bar-built estuaries exist between these two extremes, and that modeling a family of morphological types within this range could offer an insight into a wide class of possible inlet system behavior. The theory implies that tidal asymmetry measured as the difference between the maximum flood/ebb velocities in the inlet throat should shift from ebb dominated (HC1) to flood dominated (HC4) [5]. The model results agree with the theory, as HC1 is ebb dominated, and HC2, HC3, and HC4 progressively become more flood dominated. Correspondingly, the idealized model setup, which is initially unbalanced in terms of total sediment volume within the estuary versus tidal prism, allows for basin-scale adjustments through the redistribution of sediment and corresponding changes in morphology and hypsometry consistent with theoretical predictions. While the final hypsometry, and associated morphology differ, all four cases tend to an equilibrium, in which sediment import or export relaxes asymptotically to a near-zero net gain or loss. 

The majority of basin sediment loss for systems that export sediment is in the sub-tidal areas and regions that are lower in the tidal frame (HC1). This tends to deepen the basin by the preferential erosion of the sub-tidal areas. This could suggest little change in creek network structure and the size and distribution of tidal flats, while deepening sub-tidal channels. As the initial bay configuration begins to switch towards sediment import (HC2), subtidal erosion continues, but sediment also deposits on shallower subtidal areas, building an extensive network of intertidal shoals. The balance between channel erosion and shoal formation tips as subtidal areas are initially deeper (HC3). This becomes more pronounced with greater area increases for bays initially devoid of inter-tidal flats (HC4), but there are also greater gains in the sub-tidal areas decreasing the average bay depth. The increase in inter-tidal areas reduces the tendency to import sediment, producing feedback that supports an equilibrium configuration and a resulting basin with more tidal flats dissected by a dendritic creek network. 

Lagoonal-type estuaries that are deeper and devoid of extensive tidal flats tend to promote sediment import and have a greater propensity to provide material that could exacerbate channel shoaling. Activities, such as wetland mining, channel deepening and other engineering projects that remove sediment and reduce the inter-tidal footprint, could increase the rate of sediment import. By the same token, lagoonal-type estuaries that deepen due to sea level rise would similarly have the potential for greater sediment import and increased shoaling rates. 

On the other hand, lagoonal-type estuaries with extensive inter-tidal flats create ebbdominated currents that promote the export of sediment. Material removed from channels during dredging and placed strategically in other energetic regions within the estuary would have a greater chance of leaving the system as opposed to returning to the channel. The results (HC1) show a greater sediment loss in the deeper section of the system, so any tendency for channel shoaling is at least partially offset by the net erosion and export of sediment.
section_name: Implications of Basin Hypsometry to Navigation

section_content: HC1 may represent several important features of the back-barrier response to systems in which the majority of new material is placed on the existing flats in the inter-tidal zone. Morphology change in HC1 is predominantly channel deepening, with a smaller volume of sediment eroded from the intertidal flats. Additionally, channels that deepen along the extensive intertidal areas of HC1 and the shallow sub-tidal areas of HC2 are deeper and extend farther back into the bay than the channels that form in the relatively deep subtidal areas of HC3 or the flat featureless bathymetry of HC4. Erosion is predominantly from the channel in HC1 and HC2, but the extensive intertidal areas of HC1 continue to erode later in the simulation after the channel has stabilized. 

In real-world systems, this pattern of channel erosion adjacent to large intertidal areas is sometimes explained by a shift in the timing of maximum ebb and flood currents, so that peak vertical shear reduces during flood and increases during ebb (e. g., [41] ). The greater ebb turbulence and shear stress can help to contribute to the increased bottom friction and sediment transport during ebb tide [41]. If peak velocities are shifted late enough in the ebb tide, many intertidal areas may be exposed and erosion may occur in the channels as opposed to occurring on the flats.
section_name: Implications of Sediment Placement in Inter-Tidal Areas to Navigation

section_content: IHC1 can be viewed as a rudimentary model of how a lagoonal-type inlet may respond to a land reclamation project, in which material is used to create new land by raising the bed elevation above tidal and storm influence. Many long-term land reclamation projects reduce overall bay volume [42], primarily in the shallow and inter-tidal areas near adjacent shorelines, while leaving the deeper areas relatively undisturbed. 

Most of the sediment loss occurs in the inter-tidal zone, but in such a way that preserves the overall inter-tidal area, i. e., the inter-tidal zone is nearly constant over the 200-year simulation (30. 5 to 29. 5 km 2 ). Sediment eroded from the inter-tidal zone is exported and bypasses the deeper sections of the bay as indicated by the lack of significant change in the hypsometry in the sub-tidal areas (except in the small region adjacent to the inlet throat). Land reclamation may increase the rate of edge erosion and trend toward steeper side slopes in channels, which could promote unstable conditions. The fact that the erosion is primarily confined to the inter-tidal zone will reduce the area of tidal flats, of which an extreme example would be the removal of all wetlands, leading to open water throughout the basin.
section_name: Implications of Land Reclamation to Navigation

section_content: The net rate of sediment flux through the inlet is a function of not only the temporal evolution of the basin morphology, in which the rate slows as the system asymptotically approaches a neutral configuration, but also the size of the inter-tidal and sub-aqueous areas (Figure 22 ). The highest sediment flux occurs when the inter-tidal area is either less than ~10% (import) or greater than ~60% (export) of the total basin area. If the inter-tidal area lies between 10% and 60%, the import and export rates are reduced and the inlet approaches a stable morphological configuration with a vanishingly small net transport. While the results are limited by the small number of conditions that were modeled, this concept can easily be developed further by modeling many systems and researching inter-tidal coverage in real inlets and correlating that to sediment import and export rates. Further investigation could reinforce this concept, leading to a new metric to determine the relationship between large-scale wetland or land reclamation and changes in channel shoaling patterns that would be very beneficial to future navigation project planning studies.
section_name: Sediment Flux versus Intertidal Area

section_content: Increasing societal pressures (e. g., population growth and urbanization) are driving land use change practices in coastal areas and potentially will affect the hydrodynamics and sediment transport in ways that might exacerbate the existing shoaling issues. To understand the potential impact of coastal development better, a numerical model was used to predict the long-term evolution of a lagoonal-type barrier island inlet under a set range of morphological conditions that transitioned from net sediment import to net sediment export. Two of the model simulations, HC1 and IHC1, were further analyzed to examine the effect of inter-tidal placement, one in which the majority of material was restricted to elevations high in the tidal frame (IHC1). 

The results indicate that the likelihood of channel shoaling may be greater for largescale inter-tidal placement projects. Inter-tidal placement tends to enhance erosion in the sub-aqueous regions and the deeper portion of the inter-tidal zone. This may reduce the rate of new wetland erosion but lead to greater sediment deposits in the deeper channels, where it could increase channel shoaling. Placing material high in the tidal frame leads to greater erosion of the inter-tidal zone, with sediment bypassing the deeper sections. This could reduce the size of tidal flats but produce less sediment deposition in navigation channels. Furthermore, the numerical model results indicated that basins devoid of intertidal areas tended to import sediment, decreasing the overall channel depth and increasing inter-tidal areas. This suggests increased channel shoaling and larger tidal flats. 

A noteworthy result is the relationship between inter-tidal area and sediment transport. If the inter-tidal area comprises a very small portion of the total area (<~10%), then the system tends to import sediment. This effectively decreases basin depth and increases the inter-tidal footprint, which shifts the system towards neutral conditions (zero net sediment import/export). If the inter-tidal area comprises a large portion of the total area (>~60%), the system tends to export sediment. This conclusion is based on a small number of numerical modeling cases, and should be viewed as preliminary. However, increasing the number of inlet types is relatively straightforward, so further examination on this concept is very feasible and could offer an efficient way to examine the impacts on navigation for a wide variety of systems using very simple metrics.
section_name: Conclusions

section_content: Funding: This research was funded by the U. S. Army Corps of Engineers Coastal Inlets Research Program (CIRP) through the Inlet Geomorphic Evolution work unit. Institutional Review Board Statement : Not applicable. Informed Consent Statement: Not applicable.

section_name: Conflicts of Interest:

section_content: Author Contributions: D. R. K. conducted the majority of the analysis and draft preparation. R. S. provided the conceptualization and the decision to use theoretical hypsometries. M. E. B. conducted the setup and the model simulations. All authors have read and agreed to the published version of the manuscript.

section_content: The authors declare no conflict of interest.
section_name: Conflicts of Interest: