Lateral Circulation in a Partially Stratiﬁed Tidal Inlet

: Using a three-dimensional, hydrostatic, primitive-equation ocean model, this study investigates the dynamics of lateral circulation in a partially stratiﬁed tidal inlet, the Barataria Pass in the Gulf of Mexico, over a 25.6 h diurnal tidal cycle. Model performance is assessed against observational data. During ﬂood tide, the lateral circulation exhibits the characteristics similar to those induced by differential advection, i.e., lateral ﬂow consists of two counter-rotating cells and is convergent at the surface. The analysis of momentum balance indicates that, in addition to the pressure gradient and vertical stress divergence, nonlinear advection and horizontal stress divergence are also important contributors. During ebb phase, the lateral circulation is mostly toward the right shoal (when looking into the estuary) for the whole water column and persisting for almost the whole period. The surface divergence suggested by the differential advection mechanism lasts for a very short period, if it ever exists. The main momentum balance across most of the transect during ebb is between the along-channel advection of cross-channel momentum and pressure gradient. The sectional averaged lateral velocity magnitude during ebb is comparable to that during ﬂood, which is different from the idealized numerical experiment result. show


Introduction
The lateral circulation in tidally dominant estuaries can be driven by various mechanisms. Nunes and Simpson [1] identified the effect of differential advection on the generation of secondary circulation. They pointed out that due to frictional retardation, the along-channel velocity is stronger in the channel than over the shoals. When acting upon an along-channel density gradient, it results in greater (smaller) density at the thalweg than at the shoals during flood (ebb) tides. This produces a cross-channel pressure gradient toward the channel (shoals) on the surface and a pressure gradient toward the shoals (channel) at the bottom during flood (ebb) tides. Thus, the lateral flows are convergent (divergent) at the surface over the deep channel and divergent (convergent) at the bottom during flood (ebb) tides. However, the surface axial convergence was only observed during flood tides in Nunes and Simpson's work.
In an idealized, narrow straight channel with weak stratification, Lerczak and Geyer [2] confirmed that secondary circulation was driven by differential advection. Differential rotation of tidal ellipse was also identified as a mechanism for axial convergence fronts [3]. Interactions between barotropic pressure gradient and bathymetry can generate convergence of lateral flow, producing flows rotating toward the channel from the shoals [4,5]. In curved estuaries, an alternative driving mechanism for lateral circulation is the centrifugal acceleration [6][7][8] and advection [9]. Winds can enhance or degrade the local-curvature-induced, two-layer flow and can drive three-layer flow [10]. Ekman-forced lateral circulation varies with the Ekman number. When the boundary layer is comparable to the channel depth (large Ekman number), lateral flow is a single circulation cell; while for thin tidal boundary layer (small Ekman number), lateral flow is complex and varies over the tidal cycle [2]. Boundary mixing on a no-flux boundary layer was confirmed to be one of the driving mechanisms of lateral circulation [2,11]. Cheng et al. [12] investigated the lateral circulation during stratified ebb tides due to the lateral baroclinic pressure gradient, which is generated by differential diffusion caused by a lateral asymmetry in vertical mixing.
In the same idealized numerical experiment, Lerczak and Geyer [2] found that lateral flow is about four times stronger during flood tides than during ebb tides. They attributed it to the interaction between the along-channel tidal currents and nonlinear advective processes over a tidal cycle. This flood-ebb asymmetry in the lateral circulation strength was also observed by Scully et al. [13] in the Hudson River estuary, where stronger lateral flows were observed during flood tides while lateral flows were suppressed during ebb tides. However, in the numerical modeling of James River Estuary, Li et al. [14] found flood-ebb asymmetry during neap tides with stronger lateral circulation on ebb, while flood-ebb asymmetry was reduced during spring tides.
Lateral circulation plays an important role in estuarine dynamics. Many observational and numerical simulation results [2,7,12,15,16] have demonstrated the existence of secondary currents, and discussed its dependence and feedback on density stratification [7,12,17,18], streamwise momentum budget [2], estuarine circulation [12,13], as well as its impacts on sediment transport and geomorphology [19][20][21]. In Barataria Pass, the area of focus for this study, observations showed that there existed a distinct asymmetry in stratification within the diurnal cycle [22]. However, in that analysis no consideration was given to the contribution of the lateral circulation to density stratification.
In this study, we use a three-dimensional (3-D), high-resolution hydrodynamic model to examine the lateral circulation structure in the partially stratified tidal inlet, Barataria Pass, which connects Barataria Bay with the continental shelf in the southeastern Louisiana. The objectives of this study are to elucidate the tidal evolution of lateral circulation and determine its driving mechanisms. The remainder of this paper is organized as follows: Section 2 describes the study area, and the configuration of the 3-D finite volume numerical ocean model. Section 3 presents the validation of the numerical model and the temporal evolution of lateral circulation over a 25.6 h diurnal tidal cycle. In Section 4, we quantify the 2-D and 3-D momentum balance and examine the driving forcing for lateral circulation. Flood-ebb variations in lateral circulation pattern are also discussed in Section 4. Finally, conclusions are given in Section 5.

Study Area
Barataria Bay ( Figure 1) is located at the southeastern Louisiana, on the western side of the Mississippi Birdfoot Delta. It is connected to the Gulf of Mexico through several tidal inlets including Barataria Pass, which is between two barrier islands, Grand Isle Island and Grand Terre Island (Figure 1b). The most significant freshwater source inside the Barataria Bay is the Davis Pond freshwater diversion. The maximum diversion discharge is~300 m 3 s −1 [23]. Barataria Pass is an 800 m wide narrow channel. It is one of the four main tidal passes of Barataria Bay, accounting for~66% of total water exchange [24]. Tidal currents account for~85% of the total flow variance in the inlet, with equal contributions from the O 1 and K 1 constituents. Tidal amplitudes for both O 1 and K 1 constituents are about 0.5 m s −1 [25]. Maximum tidal currents reach as high as 2 m s −1 during tropic tides.
Barataria Bay is composed of broad shallow waters (average depth of 2 m), islands and a 5 m deep main shipping channel, the Barataria Waterway. The shipping channel ends at Barataria Pass. Around Barataria Pass the main channel has an average depth of~20 m and is being periodically dredged, causing a depression of~50 m deep close to the inlet ( Figure 2). A cross-sectional view of the Barataria Pass, along which all analyses are performed, is shown in Figure 1c. Bathymetrically this transect has one 20 m deep channel in the center bordered by extensive shoal regions. The shoal at the left-hand side (looking upstream) is~2 m deep and 250 m wide, while that at the right-hand side is~3 m deep and 100 m wide.

Study Area
Barataria Bay (Figure 1) is located at the southeastern Louisiana, on the western side of the Mississippi Birdfoot Delta. It is connected to the Gulf of Mexico through several tidal inlets including Barataria Pass, which is between two barrier islands, Grand Isle Island and Grand Terre Island (Figure 1b). The most significant freshwater source inside the Barataria Bay is the Davis Pond freshwater diversion. The maximum diversion discharge is ~300 m 3 s −1 [23]. Barataria Pass is an 800 m wide narrow channel. It is one of the four main tidal passes of Barataria Bay, accounting for ~66% of total water exchange [24]. Tidal currents account for ~85% of the total flow variance in the inlet, with equal contributions from the O1 and K1 constituents. Tidal amplitudes for both O1 and K1 constituents are about 0.5 m s −1 [25]. Maximum tidal currents reach as high as 2 m s −1 during tropic tides.
Barataria Bay is composed of broad shallow waters (average depth of 2 m), islands and a 5 m deep main shipping channel, the Barataria Waterway. The shipping channel ends at Barataria Pass. Around Barataria Pass the main channel has an average depth of ~20 m and is being periodically dredged, causing a depression of ~50 m deep close to the inlet ( Figure 2). A cross-sectional view of the Barataria Pass, along which all analyses are performed, is shown in Figure 1c. Bathymetrically this transect has one 20 m deep channel in the center bordered by extensive shoal regions. The shoal at the left-hand side (looking upstream) is ~2 m deep and 250 m wide, while that at the right-hand side is ~3 m deep and 100 m wide.

Model Description and Configuration
The Finite Volume Coastal Ocean Model (FVCOM) is used in this study to simulate the hydrodynamics of the Barataria Basin and adjacent continental shelf. FVCOM is a 3-D, hydrostatic, free surface, primitive-equation ocean model [26,27]. In the finite volume method, the computational domain is discretized using a mesh of non-overlapping triangles in the horizontal and sigma-coordinate (σ-coordinate) in the vertical. The governing equations are solved in their integral forms in each individual control volume. The triangular grid in the horizontal can resolve complex coastal and bathymetric geometries. It uses a cell-vertex-centered (similar to the finite-difference C-grid) method, which facilitates the enforcement of mass conservation in tracer advection and tracer open boundary conditions. Vertical mixing uses modified Mellor and Yamada level 2.5 turbulence model [28,29] and horizontal diffusion uses Smagorinsky eddy parameterization [30]. The model employs mode split approaches (barotropic 2D (external) mode and baroclinic 3D (internal) mode) to solve the momentum equations with second-order accuracy. The bottom boundary conditions apply an exact form of the no flux boundary conditions.
A flooding-drying scheme is implemented in FVCOM to simulate motions in intertidal zones and wetlands. If vertical water column thickness at the cell center is less than a criterion value (typically 5 cm), then the cell is designated as a dry cell and its velocity is set to zero. Whenever the vertical water column thickness exceeds the criterion value, the cell becomes wet and water level and velocity are computed from control equations. The advantage of triangular mesh to accurately represent complex bathymetry and coastlines makes FVCOM ideally suited for Barataria Pass study.
A high-resolution FVCOM Barataria Pass model was developed by configuring FVCOM version 2.6 to the Northwestern Gulf of Mexico continental shelf region with inclusion of the intertidal zones inside the Barataria Bay. The computational domain extends longitudinally from Mobile Bay (Alabama) to west of Galveston Bay (Texas) and offshore to about 27 • N ( Model bathymetry was obtained from various sources. Using an inverse distance weighted interpolation method, a 5 m by 5 m resolution digital elevation model constructed from Light Detection and Ranging (LIDAR) measurement was interpolated into model wetland region. The water depth in channels, bayous, and lakes was interpolated from Coastal Louisiana Ecosystem Assessment and Restoration Report (CLEAR), US Army Corp Survey, and NOAA nautical charts. Shelf and open ocean water depth was interpolated from a coarse resolution ADCIRC model bathymetry. The water depth values in the inlet were obtained by vessel-based surveys [22,31].
Since temperature difference is minute across much of our modeling region during summer time [32], salinity is considered to be the most important factor that influences water density and vertical stratification in Barataria Bay and adjacent coastal oceans in this study. Thus it is the only prognostic tracer variable in FVCOM simulations. Temperature is kept as a spatial and temporal constant (20 • C). The coefficients for horizontal viscosity and diffusivity are both set to be 0.4 m 2 s −1 . The conventional quadratic bottom friction formulation is applied, with drag coefficient C d determined by matching a logarithmic bottom boundary layer velocity to that of the numerical model at the lowest sigma-layer height. However, bottom drag coefficient over the wetlands is defined as five times greater than that in the estuarine channels, mimicking the vegetation damping effect [33].

Model Forcing, Initial, and Boundary Conditions
FVCOM is driven by winds at the surface, sea level elevation at the open boundary, and freshwater inflows from various Mississippi River and Atchafalaya River passes, and the Davis Pond Diversion.
It is initialized on 1 October 2007 and run until 31 December 2008. We use 3-hourly wind data from NOAA National Centers for Environmental Prediction (NCEP) North American Regional Reanalysis (NARR) products and interpolate it onto the entire computational domain. The initial values of sea level elevation and velocity are specified as zero throughout the computational domain. The initial salinity field over the continental shelf is interpolated from HYCOM Gulf of Mexico 1/25 • reanalysis product, while salinity inside the Barataria Bay is interpolated from observations at the United States Geological Survey (USGS) stations. Experiments show that this technique gives more accurate salinity simulation result than that with linear interpolation from estuarine head to the mouth. At the open boundary, the 6-min interval sea level time series at four stations are downloaded from NOAA tides and currents website. Time series of Dauphin Island, Southwest Pass, Freeport, and Galveston Pier 21 are directly used to prescribe sea surface elevations at the easternmost node, the southeastern node, the southwestern node, and the westernmost node, respectively. Sea level elevation at other open boundary nodes are piecewise linearly interpolated from these four nodes. Observed 15-min freshwater discharge at four locations, Mississippi River at Belle Chasse, David Pond diversion, Atchafalaya River at Morgan City, and Wax Lake Outlet, are injected into the computational domain with flux boundary conditions of zero salinity and specified volume and momentum. All model forcing functions are ramped up from zero over a period of 10 days.

Observations
Simulated water elevation, velocity and salinity are compared with in situ observations. Water elevation data is obtained from National Water Information System of USGS shown in Figure 1, including six stations, Barataria Bay near Grand Terre Island (S1), Barataria Bay North of Grand Isle (S2), Hackberry Bay near NW of Grand Isle (S3), Barataria Waterway (S4), Little Lake near Bay Dosgris (S5), and Little Lake near Bay Cutoff (S6). The unit of water elevation is converted to meter and vertical datum is adjusted to mean sea level in order to be compared with model results. Velocity and salinity data are from the field observations conducted at Barataria Pass from 11:30 31 July to 11:10 1 August 2008 UTC (see Li et al. [31] for details about field observation and data processing).

Analysis Methods
The vertically averaged cross-and along-channel momentum equations are written as: where (u, v) are the vertical integrated cross-and along-axis velocity components. The positive u is pointed to the eastern bank, the positive v to the upstream. Terms from the left to the right are local acceleration (DDT), nonlinear advection (ADV), Coriolis force (COR), barotropic pressure gradient (DPBP), baroclinic pressure gradient (DPBC), wind stress (WIND), bottom friction (FRIC), horizontal diffusion (HDIF), and difference between nonlinear terms of vertically-averaged 2-D variables and vertical integration of 3-D variables (AV2D). The expressions for HDIF and AV2D can be referred to [26,27]. Consistent with 3-D currents converted to cross-and along-channel directions, all terms in Equations (1) and (2) are rotated from the original FVCOM computation x-y coordinates (see Appendix A for details). A 3-D FVCOM momentum equation analysis was used to identify the mechanisms that drive lateral circulation. The equation is written as: where (u, v) are cross-and along-channel velocity components. The advection terms are moved to the same side of the pressure gradient, and each term in Equations (3) and (4) is calculated with its corresponding sign (e.g., −∂(u 2 D)/D∂x).

Two-Month Water Elevation Comparisons
We use the Pearson correlation coefficient [34] as evaluation metrics. Two-month (1 July to 31 August 2008) time series of water elevation for both observation and model simulation are shown in Figure 3. The model reproduces the observed tidal variations, including tropic-equatorial modulation, as well as wind-driven water level set-up and set-down with all correlation coefficients above 0.9.  (1) and (2) are rotated from the original FVCOM computation x-y coordinates (see Appendix A for details). A 3-D FVCOM momentum equation analysis was used to identify the mechanisms that drive lateral circulation. The equation is written as: where ( , ) are cross-and along-channel velocity components. The advection terms are moved to the same side of the pressure gradient, and each term in Equations (3) and (4) is calculated with its corresponding sign (e.g., − ( )⁄ ).

Two-Month Water Elevation Comparisons
We use the Pearson correlation coefficient [34] as evaluation metrics. Two-month (1 July to 31 August 2008) time series of water elevation for both observation and model simulation are shown in Figure 3. The model reproduces the observed tidal variations, including tropic-equatorial modulation, as well as wind-driven water level set-up and set-down with all correlation coefficients above 0.9.

Velocity
Similar to the treatment in Li et al. [22], the xand y-velocity components were counterclockwise rotated by 52.7 • from the east and north direction to obtain the cross-and along-channel velocity components (Figure 1b), which are shown in Figure 4. Positive along-channel velocity is flood current. As shown in Figure 4a, the observed along-channel velocity has a stronger tidal signal, with maximum magnitude of~1.5 m s −1 , than the model simulated velocity, which has a maximum magnitude of 1.0 m s −1 . The tidal phase is in agreement with the observations. Both observed and modeled cross-channel velocities (Figure 4b) are much smaller and noisier, and the tidal signal is not clear compared to the along channel velocity component. The discrepancy between the observed and modeled velocity is, in part, due to the fact that observed velocity data points were chosen along a 530 m long transect within a 90 m band (45 m on each side [22]), while modeled velocity data points are exactly along the chosen transect.

Velocity
Similar to the treatment in Li et al. [22], the x-and y-velocity components were counterclockwise rotated by 52.7° from the east and north direction to obtain the cross-and along-channel velocity components (Figure 1b), which are shown in Figure 4. Positive along-channel velocity is flood current. As shown in Figure 4a, the observed along-channel velocity has a stronger tidal signal, with maximum magnitude of ~1.5 m s −1 , than the model simulated velocity, which has a maximum magnitude of ~1.0 m s −1 . The tidal phase is in agreement with the observations. Both observed and modeled cross-channel velocities (Figure 4b) are much smaller and noisier, and the tidal signal is not clear compared to the along channel velocity component. The discrepancy between the observed and modeled velocity is, in part, due to the fact that observed velocity data points were chosen along a 530 m long transect within a 90 m band (45 m on each side [22]), while modeled velocity data points are exactly along the chosen transect.

Vertical Salinity Profile
A seabird Electronics Conductivity, Temperature, and Depth Sensor (CTD) was used for vertical profiles of salinity. A total of 28 CTD casts were made between 31 July and 1 August 2008 over a 25.6 h period (see Table 1 in Li et al. [31] for details). Vertical salinity profiles from FVCOM simulation at three locations, eastern, middle, and western side of the channel (see Figure 1c for station locations), are compared with the CTD measurements in Figure 5. The magnitude of observed salinity ranges between 19 and 28.5. The maximum vertical salinity difference is about 5.5. The magnitude of the simulated salinity varies between 15 and 27. Generally, the model underestimates salinity, perhaps partially due to the neglect of evaporation. However, the model successfully captures characteristic features in salinity vertical profile. For example, cast 6 ( Figure  5a), which was made 4 h after maximum flood, has a weak stratification at the top of the water column and well-mixed state at the bottom. Our model captures this feature. Other casts, such as casts 19, 26, 18, 8, and 16, have similar vertical profiles with the simulation. The temporal evolution of modeled salinity is also consistent with the observations. For example, the sequence of cast numbers (from low to high) for both the observation and the model show that the salinity tends to decrease in Figure 5f (ebb tide, west station). This gives us confidence to apply the model to do further qualitative dynamic analysis.

Vertical Salinity Profile
A seabird Electronics Conductivity, Temperature, and Depth Sensor (CTD) was used for vertical profiles of salinity. A total of 28 CTD casts were made between 31 July and 1 August 2008 over a 25.6 h period (see Table 1 in Li et al. [31] for details). Vertical salinity profiles from FVCOM simulation at three locations, eastern, middle, and western side of the channel (see Figure 1c for station locations), are compared with the CTD measurements in Figure 5. The magnitude of observed salinity ranges between 19 and 28.5. The maximum vertical salinity difference is about 5.5. The magnitude of the simulated salinity varies between 15 and 27. Generally, the model underestimates salinity, perhaps partially due to the neglect of evaporation. However, the model successfully captures characteristic features in salinity vertical profile. For example, cast 6 (Figure 5a), which was made 4 h after maximum flood, has a weak stratification at the top of the water column and well-mixed state at the bottom. Our model captures this feature. Other casts, such as casts 19, 26, 18, 8, and 16, have similar vertical profiles with the simulation. The temporal evolution of modeled salinity is also consistent with the observations. For example, the sequence of cast numbers (from low to high) for both the observation and the model show that the salinity tends to decrease in Figure 5f (ebb tide, west station). This gives us confidence to apply the model to do further qualitative dynamic analysis.

Temporal Variation of Stratification in the Barataria Pass
In order to study the variation of stratification over a diurnal cycle, a 25.6 h time series from the FVCOM simulation, starting from 06:00 31 July to 07:40 1 August 2008, was chosen to perform further analysis. Time series of water level, depth-averaged along-channel velocity and salinity difference between bottom and surface at three cross-channel locations are shown in Figure 6. The water elevation and depth-averaged velocity are not in-phase. Based on the deep channel station, flood tide lasts for ~12 h with maximum vertically averaged flood velocity ~0.8 m/s, while ebb tide lasts for ~13.6 h and maximum vertically averaged ebb velocity ~1.2 m/s. This tidal asymmetry is mainly caused by upstream Davis Pond diversion discharge, which is ~200 m 3 /s during this period of time. Stratification evolution at these three locations shows distinct cross-channel variation. Within 2 h of early flood tides, stratification at all three locations decreases and reaches a well-mixed condition. Then stratification starts increasing. First the station at the western shoal quickly reaches the maximum (Figure 6a), followed by the deep channel station (Figure 6b). The station at the eastern shoal has a more moderate increase rate (Figure 6c). During the remaining period of flood tide, the western shoal experiences variation between well-mixed and stratified conditions. The deep channel and the eastern shoal are always stratified, and the latter has the largest stratification except near the end of flood. During ebb tide stratification at the western shoal is the weakest and remains almost well-mixed for the whole ebb tide. Stratification at the deep channel and eastern shoal have similar evolution, decreasing in the beginning, reaching well-mixed condition 5-6 h after ebbing, and increasing again 3 h before the slack water.

Temporal Variation of Stratification in the Barataria Pass
In order to study the variation of stratification over a diurnal cycle, a 25.6 h time series from the FVCOM simulation, starting from 06:00 31 July to 07:40 1 August 2008, was chosen to perform further analysis. Time series of water level, depth-averaged along-channel velocity and salinity difference between bottom and surface at three cross-channel locations are shown in Figure 6. The water elevation and depth-averaged velocity are not in-phase. Based on the deep channel station, flood tide lasts for 12 h with maximum vertically averaged flood velocity~0.8 m/s, while ebb tide lasts for~13.6 h and maximum vertically averaged ebb velocity~1.2 m/s. This tidal asymmetry is mainly caused by upstream Davis Pond diversion discharge, which is~200 m 3 /s during this period of time. Stratification evolution at these three locations shows distinct cross-channel variation. Within 2 h of early flood tides, stratification at all three locations decreases and reaches a well-mixed condition. Then stratification starts increasing. First the station at the western shoal quickly reaches the maximum (Figure 6a), followed by the deep channel station (Figure 6b). The station at the eastern shoal has a more moderate increase rate (Figure 6c). During the remaining period of flood tide, the western shoal experiences variation between well-mixed and stratified conditions. The deep channel and the eastern shoal are always stratified, and the latter has the largest stratification except near the end of flood. During ebb tide stratification at the western shoal is the weakest and remains almost well-mixed for the whole ebb tide. Stratification at the deep channel and eastern shoal have similar evolution, decreasing in the beginning, reaching well-mixed condition 5-6 h after ebbing, and increasing again 3 h before the slack water.

Residual Currents in the Barataria Pass over One Tidal Cycle
Tidal, ebb-, and flood-averaged along-channel velocities for the same 25.6 h time period are shown in Figure 7. The location of the transect is shown in Figure 1b and the western shoal is on the left side of Figure 7. The magnitude of cross-sectionally averaged ebb tides is −0.46 m/s, and that of cross-sectionally averaged flood tides is 0.33 m/s. The transverse structure of the along-channel residual current differs significantly between ebb and flood tides. For ebb tide, the maximum outflow is at the surface near the western shoal, inclining against the western bank of the deep channel. For flood tide, the maximum inflow is near the mid-depth in the central deep channel. Vertical velocity shear is much larger along the eastern bank than along the western bank. The magnitude of spatially averaged residual current during the 25.6 h period is −0.09 m/s, which is in the ebb direction. The maximum residual current is near the western shoal with the magnitude close to −0.3 m/s. Using idealized numerical experiments, Cheng and Valle-Levinson [35] studied the sensitivity of estuarine exchange flow pattern on two nondimensional parameters, the Rossby number ( / , where U is the estuarine circulation velocity, is the Coriolis parameter, and B is the width of the channel) and the Ekman number ( / , where is the vertical eddy viscosity, H is water depth, and is the Coriolis parameter). They demonstrated that the exchange flow is vertically sheared at large , and horizontally sheared at large . In our case ( ≈ 1.72, ≈

Residual Currents in the Barataria Pass over One Tidal Cycle
Tidal, ebb-, and flood-averaged along-channel velocities for the same 25.6 h time period are shown in Figure 7. The location of the transect is shown in Figure 1b and the western shoal is on the left side of Figure 7. The magnitude of cross-sectionally averaged ebb tides is −0.46 m/s, and that of cross-sectionally averaged flood tides is 0.33 m/s. The transverse structure of the along-channel residual current differs significantly between ebb and flood tides. For ebb tide, the maximum outflow is at the surface near the western shoal, inclining against the western bank of the deep channel. For flood tide, the maximum inflow is near the mid-depth in the central deep channel. Vertical velocity shear is much larger along the eastern bank than along the western bank. The magnitude of spatially averaged residual current during the 25.6 h period is −0.09 m/s, which is in the ebb direction. The maximum residual current is near the western shoal with the magnitude close to −0.3 m/s.
Using idealized numerical experiments, Cheng and Valle-Levinson [35] studied the sensitivity of estuarine exchange flow pattern on two nondimensional parameters, the Rossby number R 0 (U/ f B, where U is the estuarine circulation velocity, f is the Coriolis parameter, and B is the width of the channel) and the Ekman number E k (A z / f H 2 , where A z is the vertical eddy viscosity, H is water depth, and f is the Coriolis parameter). They demonstrated that the exchange flow is vertically sheared at large R 0 , and horizontally sheared at large E k . In our case (R 0 ≈ 1.72, E k ≈ 2.06), the residual current is both vertically and horizontally sheared, which is similar to their case shown in Figure 5c (R 0 = 2.63, E k = 1) of [35]. For a triangular shaped cross section, Wong [36] showed that the estuarine circulation is outward at the surface and inward at the bottom of the deep channel due to the interaction between baroclinic force and triangular bathymetry. Field observations at Barataria Pass captured this characteristics [31], in which the inflow is very weak. Our model also reproduces a weak inflow near the bottom of the eastern slope (Figure 7c). The asymmetry of estuarine-ocean exchange, i.e., inflow tending to the right side of the channel and outflow to the left side (when looking up-estuary) may be attributed to the Coriolis force [37]. 2.06), the residual current is both vertically and horizontally sheared, which is similar to their case shown in Figure 5c ( = 2.63, = 1) of [35]. For a triangular shaped cross section, Wong [36] showed that the estuarine circulation is outward at the surface and inward at the bottom of the deep channel due to the interaction between baroclinic force and triangular bathymetry. Field observations at Barataria Pass captured this characteristics [31], in which the inflow is very weak. Our model also reproduces a weak inflow near the bottom of the eastern slope ( Figure 7c). The asymmetry of estuarine-ocean exchange, i.e., inflow tending to the right side of the channel and outflow to the left side (when looking up-estuary) may be attributed to the Coriolis force [37].   The transverse structure of cross-channel residual current also differs significantly between ebb and flood tides. The maximum flood-averaged and ebb-averaged velocity can reach 0.2 m/s and 0.16 m/s, respectively. The ebb-averaged velocity exhibits no horizontal divergence or convergence in the deep channel (Figure 8b,e), while the flood-averaged velocity shows a strong convergence close to the surface in the mid-channel (Figure 8a,d). However, no closed cell circulation exists for the flood average (Figure 8d), which is different from that described in Nunes and Simpson [1]. The tidally averaged cross-channel velocity is relatively weaker (Figure 8c), and it displays a clockwise circulation cell in the western channel and an even weaker counterclockwise cell in the eastern channel (Figure 8f). 2.06), the residual current is both vertically and horizontally sheared, which is similar to their case shown in Figure 5c ( = 2.63, = 1) of [35]. For a triangular shaped cross section, Wong [36] showed that the estuarine circulation is outward at the surface and inward at the bottom of the deep channel due to the interaction between baroclinic force and triangular bathymetry. Field observations at Barataria Pass captured this characteristics [31], in which the inflow is very weak. Our model also reproduces a weak inflow near the bottom of the eastern slope (Figure 7c). The asymmetry of estuarine-ocean exchange, i.e., inflow tending to the right side of the channel and outflow to the left side (when looking up-estuary) may be attributed to the Coriolis force [37]. The transverse structure of cross-channel residual current also differs significantly between ebb and flood tides. The maximum flood-averaged and ebb-averaged velocity can reach 0.2 m/s and 0.16 m/s, respectively. The ebb-averaged velocity exhibits no horizontal divergence or convergence in the deep channel (Figure 8b, e), while the flood-averaged velocity shows a strong convergence close to the surface in the mid-channel (Figure 8a, d). However, no closed cell circulation exists for the flood average (Figure 8d), which is different from that described in Nunes and Simpson [1]. The tidally averaged cross-channel velocity is relatively weaker (Figure 8c), and it displays a clockwise circulation cell in the western channel and an even weaker counterclockwise cell in the eastern channel ( Figure 8f).   Figure 9 shows the temporal evolution of along-channel velocity, lateral circulation, salinity, and turbulent vertical eddy viscosity for the Barataria Pass transect during flood tide. Each row represents a time instance indicated by an arrow in Figure 6. One hour after the flood starts (T1 in Figure 6), stratification is weak (bottom-top salinity difference~1) across the deep channel and eastern shoal, and the western shoal is almost well-mixed (Figure 9c). Strong lateral circulation mainly occurs in the mid-depth with a convergence zone below 10 m in the deep channel (Figure 9b). The along-channel velocity is~0.4 m/s, extending almost the whole deep channel. Thus, vertical and lateral velocity shear is weak. The vertical eddy viscosity is mostly smaller than 0.005 m 2 /s, probably due to the small flood current magnitude and its shears.

Time Series of Cross-Sectional Salinity and Currents Structures
Two hours later (T2 in Figure 6), the whole water column becomes more or less well-mixed across the deep channel and eastern shoal, while in the western shoal a sharp salinity stratification develops with bottom-top salinity difference~3 in 2 m water column (Figure 9g). The distribution of salinity, and thus density, in the cross-channel direction is such that lateral baroclinic pressure gradients are directed from the central part of the deep channel towards the shoals. This is similar to the situation pointed out in Nunes and Simpson [1]. Based on their theory, such pressure gradient force will induce a lateral circulation with convergence at the surface and divergence near the bottom. This indeed occurs in our numerical results. Surface lateral flow at the west half of the channel changes from 0.1 m/s westward to~0.2 m/s eastward between T1 and T2. A pair of counter-rotating circulations can be clearly seen at T2 with strong convergence 2 m below the surface. Contrast to the idealized case in Lerczak and Geyer [2], the two circulation cells are not closed at this time. The maximum along-channel tidal velocity reaches~0.8 m/s, confined at the mid-depth of the central channel ( Figure 9e). The vertical shear in along-channel velocity is relatively weak, while the horizontal shear is great, which is, over the western slope,~0.6 m/s within 200 m distance. It follows that Nunes and Simpson's argument is also applicable here in that differential advection of along-channel current is at least one of the mechanisms to generate the lateral salinity gradient. Strong vertical mixing (maximum eddy viscosity~0.05 m 2 /s, Figure 9h) occurs at the mid-depth and bottom boundary layer, where either tidal currents or bottom friction are strong (Figure 9e). Strong turbulence mixing tends to destratify the water column, which explains the relatively uniform salinity distribution in the deep channel and east shoal (Figure 9g). At the maximum flood (T3), the along-channel velocities intensify. Maximum along-channel velocity, located at channel thalweg, reaches~1.2 m/s and extends to the surface (Figure 9i), which is quite different from T2 (Figure 9e). The lateral shear of along-channel velocity, i.e., differential advection, is 1.0 m/s across 300 m distance on both sides. However, salinity distribution changes drastically. The western shoal is completely well-mixed at this time, while the eastern shoal and part of the east channel have surface stratification with a salinity difference of~3 within 5-6 m depth (Figure 9k). The stratification in the deep channel, especially below 10 m, is weak, because the tidal mixing is relatively high (vertical eddy viscosity~0.03 m 2 /s) at the bottom boundary (Figure 9l). The lateral circulation pattern is similar to that at T2, but with intensified strength. The convergent zone rises to the ocean surface and the right circulation cell is now a complete circle (Figure 9j).
Four hours before the end of flood (T4), differential advection still persists, although the maximum along-channel tidal current has reduced to~1 m/s and lowered below the surface (Figure 9m). Turbulence mixing weakens (Figure 9p). Hence, a weak stratification (bottom-surface salinity difference~1.5) develops at the western shoal (Figure 9o). The salinity distribution shows a more symmetric pattern relative to the axis of the channel compared with T1, T2, and T3. As a result, the counter-rotating lateral circulation cells are more symmetric and fully developed (Figure 9n).
At the flood slack (T5), the surface water column becomes stratified in the upper 6 m (salinity difference~1) while the deep channel is almost well-mixed (Figure 9s). The lateral circulation almost completely disappears in the deep channel (Figure 9r). At the beginning of ebb tide (T6), although the along-channel ebb current increases to ~0.4 m/s (Figure 10a), the salinity distribution and vertical stratification (Figure 10c) are almost the same as that of T5. There is a weak (less than 0.1 m/s) eastward lateral flow below 6 m (Figure 10b).
Two hours later (T7), the along-channel ebb current reaches above 1.0 m/s, locating mostly on the western slope and reaching the sea surface. The magnitude of along-channel velocity across the majority of the cross section is ~0.8 m/s (Figure 10e), which induces large tidal mixing (maximum vertical eddy viscosity ~0.16 m 2 /s, Figure 10h). Thus, the whole water column is vertically well-mixed (Figure 10g). However, a horizontal salinity gradient exists, with higher salinity located near the channel axis, fresher water on both shoals. The western shoal is fresher than the eastern shoal (Figure 10g). This is because freshwater is flushed out of the estuary through the western shoal, as shown in Figure 7b,c. The lateral circulation shows mostly eastward currents in the deep channel across the whole water column, while the eastern shoal has a convergence area (Figure 10f).

Figure 9.
Cross-sectional profiles of currents (u, v, w), salinity and vertical viscosity during flood tide. The first column (a,e,i,m,q) is along-channel velocity, the second column (b,f,j,n,r) secondary circulation, the third column (c,g,k,o,s) salinity, and the last column (d,h,l,p,t) vertical viscosity. The velocity contour interval is 0.2 m/s, positive is up-estuary. The salinity contour interval is 0.5. Each row corresponds to a time instance indicated by an arrow in Figure 6.
At the beginning of ebb tide (T6), although the along-channel ebb current increases to~0.4 m/s (Figure 10a), the salinity distribution and vertical stratification (Figure 10c) are almost the same as that of T5. There is a weak (less than 0.1 m/s) eastward lateral flow below 6 m (Figure 10b).
Two hours later (T7), the along-channel ebb current reaches above 1.0 m/s, locating mostly on the western slope and reaching the sea surface. The magnitude of along-channel velocity across the majority of the cross section is~0.8 m/s (Figure 10e), which induces large tidal mixing (maximum vertical eddy viscosity~0.16 m 2 /s, Figure 10h). Thus, the whole water column is vertically well-mixed ( Figure 10g). However, a horizontal salinity gradient exists, with higher salinity located near the channel axis, fresher water on both shoals. The western shoal is fresher than the eastern shoal ( Figure 10g). This is because freshwater is flushed out of the estuary through the western shoal, as shown in Figure 7b,c. The lateral circulation shows mostly eastward currents in the deep channel across the whole water column, while the eastern shoal has a convergence area (Figure 10f).
During the next five hours, this cross section is always vertically well-mixed and salinity decreases constantly due to freshwater outflow. The turbulence mixing remains intense in the deep channel during this period. The maximum vertical eddy viscosity can reach up to 0.2 m 2 /s, which results in the water column in the east half of the channel vertically and horizontally well-mixed (Figure 10k, T8). The water column in the west half of the channel is also vertically well-mixed, but has a weak (~1) horizontal salinity decrease westward. The structure of lateral circulation and along-channel velocity remains the same, although the maximum ebb velocity has decreased from 1.0 m/s at T7 to 0.8 m/s at T8. During the next five hours, this cross section is always vertically well-mixed and salinity decreases constantly due to freshwater outflow. The turbulence mixing remains intense in the deep channel during this period. The maximum vertical eddy viscosity can reach up to 0.2 m 2 /s, which results in the water column in the east half of the channel vertically and horizontally well-mixed (Figure 10k, T8). The water column in the west half of the channel is also vertically well-mixed, but has a weak (~1) horizontal salinity decrease westward. The structure of lateral circulation and along-channel velocity remains the same, although the maximum ebb velocity has decreased from 1.0 m/s at T7 to 0.8 m/s at T8.
Later during the ebb period (T9), the vertical stratification returns (Figure 10o), probably due to greatly reduced turbulence vertical eddy viscosity (Figure 10p). The vertical salinity difference is ~2 in the deep channel and on the eastern shoal. Lateral circulation shows a weak flow divergence close Figure 10. Cross-sectional profiles of currents (u, v, w), salinity and vertical viscosity during ebb tide. The first column (a,e,i,m,q) is along-channel velocity, the second column (b,f,j,n,r) secondary circulation, the third column (c,g,k,o,s) salinity, and the last column (d,h,l,p,t) vertical viscosity. The velocity contour interval is 0.2 m/s, positive is landward. The salinity contour interval is 0.5. Each row corresponds to a time instance indicated by an arrow in Figure 6.
Later during the ebb period (T9), the vertical stratification returns (Figure 10o), probably due to greatly reduced turbulence vertical eddy viscosity (Figure 10p). The vertical salinity difference is~2 in the deep channel and on the eastern shoal. Lateral circulation shows a weak flow divergence close to the surface near the western slope (Figure 10n). The maximum along-channel ebb current moves to the surface central channel and its magnitude decreases to 0.6 m/s (Figure 10m).
At the ebb slack (T10), salinity distribution go back to similar to T1 situation. The western shoal is almost vertically uniform. Strong stratification (salinity difference~2 over 6 m depth) occurs in the upper water column of the deep channel and on the eastern shoal. The lower water column in the deep channel has a weak stratification (Figure 10s). Lateral circulation is greatly reduced compared to other time instances of the ebb tide (Figure 10r).

Depth-Averaged Momentum Balance
The Coriolis force, wind stress, and horizontal diffusion in Equations (1) and (2) are at least one order of magnitude smaller than the other terms during the 25.6 h period. Thus, these three terms are not shown. Figure 11 shows time series of six terms (DDT, ADV, DPBP, DPBC, FRIC, and AV2D) in Equations (1) and (2) at seven locations across the channel (Figure 1c). The momentum balance across this narrow channel is very complex, as various locations have different characteristics.
to the surface near the western slope (Figure 10n). The maximum along-channel ebb current moves to the surface central channel and its magnitude decreases to 0.6 m/s (Figure 10m).
At the ebb slack (T10), salinity distribution go back to similar to T1 situation. The western shoal is almost vertically uniform. Strong stratification (salinity difference ~2 over 6 m depth) occurs in the upper water column of the deep channel and on the eastern shoal. The lower water column in the deep channel has a weak stratification (Figure 10s). Lateral circulation is greatly reduced compared to other time instances of the ebb tide (Figure 10r).

Depth-Averaged Momentum Balance
The Coriolis force, wind stress, and horizontal diffusion in Equations (1) and (2) are at least one order of magnitude smaller than the other terms during the 25.6 h period. Thus, these three terms are not shown. Figure 11 shows time series of six terms (DDT, ADV, DPBP, DPBC, FRIC, and AV2D) in Equations (1) and (2) at seven locations across the channel (Figure 1c). The momentum balance across this narrow channel is very complex, as various locations have different characteristics.  In the along-channel momentum balance, the dominant balance is between the barotropic pressure gradient and the nonlinear advection, especially during ebb tide. The magnitudes of these two terms for stations on the western side (C2 and C3, Figure 11i,j) of the channel are at least twice as greater as that of other stations. This is because maximum ebb currents flush out of the Barataria Bay near these two stations. The sign of barotropic pressure and nonlinear advection at stations on the west channel (C1-C3, Figure 11h-j) is opposite to those on the east channel (C4-C7, Figure 11k-n). There is a spike during ebb tide, similar to that in Huang et al. [33]. The reason for this spike is not clear. During flood tide, the magnitudes of all terms at stations in the western channel (C1-C3, Figure 11h-j) are relatively small. This is because the sign of terms in the upper layer is opposite to that of the lower layer during flood tide. They offset each other after integrating over depth. The balance is among the DDT, nonlinear advection, barotropic pressure gradient, baroclinic pressure gradient, and the AV2D.
In the cross-channel momentum balance, the characteristics are similar to that of along-channel, i.e., the dominant balance is between advection term and the barotropic pressure gradient. Excluding station C1, the signs of advection and barotropic pressure terms at the westmost station (C2, Figure 11b) and the eastmost station (C7, Figure 11g) are the same, but opposite to the stations of the deep channel (C3, C4, Figure 11c,d). During flood tides, the balance is among the DDT, nonlinear advection, barotropic pressure gradient, baroclinic pressure gradient, and the AV2D. Note that, the baroclinic pressure is great at stations in the deep channel, and the magnitude is larger than that of along-channel. This indicates the lateral salinity gradient play an important role in momentum balance.

Driving Mechanism of Lateral Circulation
Lerczak and Geyer [2] pointed out that lateral advection plays an important role in the estuarine dynamics when lateral flows are strong enough to advect water parcels relative to 0.5 times the breadth of the channel (4|v|/σB ≥ 1, where |v| is the absolute value of lateral velocity amplitude, σ is the semidiurnal tidal frequency, and B is the channel width) in a tidal cycle. As shown in previous results, lateral circulation in the Barataria Pass is strong both during maximum flood and maximum ebb. Thus, lateral advection is expected to be an important term in the momentum balance. The 3-D momentum Equations (3) and (4) are used to explore the generation mechanisms of the lateral circulation.
two terms for stations on the western side (C2 and C3, Figure 11i,j) of the channel are at least twice as greater as that of other stations. This is because maximum ebb currents flush out of the Barataria Bay near these two stations. The sign of barotropic pressure and nonlinear advection at stations on the west channel (C1-C3, Figure 11h-j) is opposite to those on the east channel (C4-C7, Figure 11k-n). There is a spike during ebb tide, similar to that in Huang et al. [33]. The reason for this spike is not clear. During flood tide, the magnitudes of all terms at stations in the western channel (C1-C3, Figure 11h-j) are relatively small. This is because the sign of terms in the upper layer is opposite to that of the lower layer during flood tide. They offset each other after integrating over depth. The balance is among the DDT, nonlinear advection, barotropic pressure gradient, baroclinic pressure gradient, and the AV2D.
In the cross-channel momentum balance, the characteristics are similar to that of along-channel, i.e., the dominant balance is between advection term and the barotropic pressure gradient. Excluding station C1, the signs of advection and barotropic pressure terms at the westmost station (C2, Figure 11b) and the eastmost station (C7, Figure 11g) are the same, but opposite to the stations of the deep channel (C3, C4, Figure 11c,d). During flood tides, the balance is among the DDT, nonlinear advection, barotropic pressure gradient, baroclinic pressure gradient, and the AV2D. Note that, the baroclinic pressure is great at stations in the deep channel, and the magnitude is larger than that of along-channel. This indicates the lateral salinity gradient play an important role in momentum balance.

Driving Mechanism of Lateral Circulation
Lerczak and Geyer [2] pointed out that lateral advection plays an important role in the estuarine dynamics when lateral flows are strong enough to advect water parcels relative to 0.5 times the breadth of the channel (4〈| |〉/ ≥ 1, where | | is the absolute value of lateral velocity amplitude, is the semidiurnal tidal frequency, and B is the channel width) in a tidal cycle. As shown in previous results, lateral circulation in the Barataria Pass is strong both during maximum flood and maximum ebb. Thus, lateral advection is expected to be an important term in the momentum balance. The 3-D momentum Equations (3) and (4) are used to explore the generation mechanisms of the lateral circulation.  Transverse distributions of various momentum terms of along-and cross-channel components at T3 (flood tide) are shown in Figures 12 and 13, respectively. In the along-channel momentum equation, baroclinic pressure gradient (Figure 12e), Coriolis force ( Figure 12h) and vertical stress divergence (Figure 12i) are at least two orders of magnitude less than other terms at this time, thus are less important. The total pressure gradient comes from the barotropic pressure gradient (Figure 12d), and the largest terms are the nonlinear advection terms (Figure 12b,c). Lateral advection of along-channel momentum (Figure 12b, −uv x − wv z ) is one of the largest terms in the momentum balance. This term has a similar pattern but opposite sign with the along-channel momentum advection (−vv y ). The main force balance at this time is between non-linear advection and the barotropic pressure gradient, which is consistent with the result in 2-D momentum analysis.
Transverse distributions of various momentum terms of along-and cross-channel components at T3 (flood tide) are shown in Figures 12 and 13, respectively. In the along-channel momentum equation, baroclinic pressure gradient (Figure 12e), Coriolis force ( Figure 12h) and vertical stress divergence (Figure 12i) are at least two orders of magnitude less than other terms at this time, thus are less important. The total pressure gradient comes from the barotropic pressure gradient (Figure  12d), and the largest terms are the nonlinear advection terms (Figure 12b,c). Lateral advection of along-channel momentum (Figure 12b, − − ) is one of the largest terms in the momentum balance. This term has a similar pattern but opposite sign with the along-channel momentum advection (− ). The main force balance at this time is between non-linear advection and the barotropic pressure gradient, which is consistent with the result in 2-D momentum analysis. In the cross-channel momentum balance, Coriolis force (Figure 13h), although small, has exactly the same pattern as the transverse distribution of along-channel velocity (Figure 9i). Differential advection of the flood current at this time is clearly shown in the picture. The baroclinic pressure gradient (Figure 13e) has almost the same magnitude as the barotropic pressure gradient (Figure 13d), especially near the bottom, due to its magnitude increasing with depth. The baroclinic pressure gradient is mostly positive in the whole channel, except near the eastern shoal. Combining with the salinity distribution at the transect (Figure 9k), it indicates that the conceptual model of cross-channel baroclinic pressure gradient induced by differential advection is applicable here. However, the simple diagnostic model proposed by Nunes and Simpson [1] to explain the mechanism of lateral circulation generation, in which the momentum balance is between cross-channel pressure gradient and vertical stress divergence, seems to be oversimplified. Our numerical experiment shows that vertical stress divergence is mostly confined to the near bottom layer of the eastern slope (Figure 13i). Several other terms-nonlinear advection (Figure 13b,c) and horizontal stress divergence (Figure 13g)-are as important as the terms suggested in [1].
During ebb tide (i.e., at T7), the main momentum balance in the along-channel direction is among the nonlinear advection terms and the barotropic pressure gradient (Figure 14), which is the same as during the flood tide ( Figure 12).
In the cross-channel momentum balance, the barotropic pressure gradient (Figure 15d) dominates the baroclinic pressure gradient (Figure 15e). Thus, contour lines of the total cross-channel pressure gradient (Figure 15f) are very similar to that of the former. It seems that at In the cross-channel momentum balance, Coriolis force (Figure 13h), although small, has exactly the same pattern as the transverse distribution of along-channel velocity (Figure 9i). Differential advection of the flood current at this time is clearly shown in the picture. The baroclinic pressure gradient (Figure 13e) has almost the same magnitude as the barotropic pressure gradient (Figure 13d), especially near the bottom, due to its magnitude increasing with depth. The baroclinic pressure gradient is mostly positive in the whole channel, except near the eastern shoal. Combining with the salinity distribution at the transect (Figure 9k), it indicates that the conceptual model of cross-channel baroclinic pressure gradient induced by differential advection is applicable here. However, the simple diagnostic model proposed by Nunes and Simpson [1] to explain the mechanism of lateral circulation generation, in which the momentum balance is between cross-channel pressure gradient and vertical stress divergence, seems to be oversimplified. Our numerical experiment shows that vertical stress divergence is mostly confined to the near bottom layer of the eastern slope (Figure 13i). Several other terms-nonlinear advection (Figure 13b,c) and horizontal stress divergence (Figure 13g)-are as important as the terms suggested in [1].
During ebb tide (i.e., at T7), the main momentum balance in the along-channel direction is among the nonlinear advection terms and the barotropic pressure gradient (Figure 14), which is the same as during the flood tide ( Figure 12).
In the cross-channel momentum balance, the barotropic pressure gradient (Figure 15d) dominates the baroclinic pressure gradient (Figure 15e). Thus, contour lines of the total cross-channel pressure gradient (Figure 15f) are very similar to that of the former. It seems that at this time the main momentum balance is between along-channel advection of cross-channel momentum (Figure 15c) and pressure gradient (Figure 15f) over most of the cross-section area, both of which have alternative positive and negative vertical stripes and, when added together, almost cancel out each other. At the eastern shoal, all terms, excluding Coriolis force (Figure 15h) and baroclinic pressure gradient (Figure 15e), are needed in order to balance the cross-channel momentum equation. Because of the insignificant contribution from the baroclinic pressure gradient, no obvious convergence or divergence occurs in this transect during the ebb period ( Figure 10, second column). This is in clear contrast to the near surface convergence that happens during flood tide (Figure 9, second column).
this time the main momentum balance is between along-channel advection of cross-channel momentum ( Figure 15c) and pressure gradient (Figure 15f) over most of the cross-section area, both of which have alternative positive and negative vertical stripes and, when added together, almost cancel out each other. At the eastern shoal, all terms, excluding Coriolis force (Figure 15h) and baroclinic pressure gradient (Figure 15e), are needed in order to balance the cross-channel momentum equation. Because of the insignificant contribution from the baroclinic pressure gradient, no obvious convergence or divergence occurs in this transect during the ebb period ( Figure  10, second column). This is in clear contrast to the near surface convergence that happens during flood tide (Figure 9, second column).   momentum ( Figure 15c) and pressure gradient (Figure 15f) over most of the cross-section area, both of which have alternative positive and negative vertical stripes and, when added together, almost cancel out each other. At the eastern shoal, all terms, excluding Coriolis force (Figure 15h) and baroclinic pressure gradient (Figure 15e), are needed in order to balance the cross-channel momentum equation. Because of the insignificant contribution from the baroclinic pressure gradient, no obvious convergence or divergence occurs in this transect during the ebb period ( Figure  10, second column). This is in clear contrast to the near surface convergence that happens during flood tide (Figure 9, second column).

Flood-Ebb Asymmetry
Here, we follow Lerczak and Geyer [2] using the cross-channel average of depth-averaged velocity amplitude < |u| > = 1 A |u|dA (|u| is the absolute value of depth-averaged cross-channel velocity, A is the cross-channel area) to represent the strength of the lateral flow. The result is shown in Figure 16. The maximum < |u| > during flood is 0.32 m/s, while the maximum during ebb is 0.34 m/s. Generally, cross-channel current amplitude during ebb is comparable, even slightly greater than that during flood, which is different from the idealized case [2], in which lateral circulation is about four times stronger during flood tides than during ebb tides. However, this inconsistence is also observed in James River estuary [14], where the lateral circulation shows no flood-ebb asymmetry during spring tides, and a reversed asymmetry during neap tides, that is stronger during ebb than during flood. Li et al. [14] attribute negligible flood-ebb variations during spring tides to turbulent mixing, which in their numerical experiment was simulated via a turbulence closure model. It reduces the vertical shear and the flood-ebb asymmetry in the vorticity generation. In comparison, the idealized case of Lerczak and Geyer [2] used a constant eddy viscosity. In our numerical experiment, turbulent mixing is simulated using the Mellor-Yamada 2.5 turbulence closure model and the 25.6 h diurnal tidal cycle is close to tropic tide. Maximum turbulent vertical viscosity is twice as great during ebb ( Figure 10, last column) as during flood (Figure 9, last column). Therefore, our simulation result is more similar to Li et al. [14] than to Lerczak and Geyer [2].

Flood-Ebb Asymmetry
Here, we follow Lerczak and Geyer [2] using the cross-channel average of depth-averaged velocity amplitude | | = ∬ | | (| | is the absolute value of depth-averaged cross-channel velocity, A is the cross-channel area) to represent the strength of the lateral flow. The result is shown in Figure 16. The maximum | | during flood is 0.32 m/s, while the maximum during ebb is 0.34 m/s. Generally, cross-channel current amplitude during ebb is comparable, even slightly greater than that during flood, which is different from the idealized case [2], in which lateral circulation is about four times stronger during flood tides than during ebb tides. However, this inconsistence is also observed in James River estuary [14], where the lateral circulation shows no flood-ebb asymmetry during spring tides, and a reversed asymmetry during neap tides, that is stronger during ebb than during flood. Li et al. [14] attribute negligible flood-ebb variations during spring tides to turbulent mixing, which in their numerical experiment was simulated via a turbulence closure model. It reduces the vertical shear and the flood-ebb asymmetry in the vorticity generation. In comparison, the idealized case of Lerczak and Geyer [2] used a constant eddy viscosity. In our numerical experiment, turbulent mixing is simulated using the Mellor-Yamada 2.5 turbulence closure model and the 25.6 h diurnal tidal cycle is close to tropic tide. Maximum turbulent vertical viscosity is twice as great during ebb ( Figure 10, last column) as during flood (Figure 9, last column). Therefore, our simulation result is more similar to Li et al. [14] than to Lerczak and Geyer [2]. During flood tide, lateral circulation shows asymmetry across the section (Figure 9f,j,n). The counterclockwise circulation on the east side is stronger, which is due to the presence of lower salinity water close to and over the eastern shoal (Figure 9k,o). This fresher water increases the lateral baroclinic pressure gradient on the east side, and thus enhances the strength of the lateral circulation cell at this side. The low salinity water is from the Mississippi River plume based on analyses of in-situ observational data and numerical simulation result [22].
Another asymmetry between flood and ebb tides lies in that maximum flood currents are always located, more or less, at the central part of the deep channel ( Figure 9, first column), while maximum ebb currents, when they are greater than 0.8 m/s, are located closer to the western shoal or over the western slope (Figure 10e,i). Only when maximum ebb currents are below 0.6 m/s, they move back to the middle of the channel (Figure 10a,m,q).

Conclusions
Barataria Pass is a tidal inlet that connects the Barataria Bay to the continental shelf. Previous investigation has introduced the tidal straining effect on density stratification during the same 25.6 h period along the same transect [22]. In this study, we conduct a numerical model simulation and During flood tide, lateral circulation shows asymmetry across the section (Figure 9f,j,n). The counterclockwise circulation on the east side is stronger, which is due to the presence of lower salinity water close to and over the eastern shoal (Figure 9k,o). This fresher water increases the lateral baroclinic pressure gradient on the east side, and thus enhances the strength of the lateral circulation cell at this side. The low salinity water is from the Mississippi River plume based on analyses of in-situ observational data and numerical simulation result [22].
Another asymmetry between flood and ebb tides lies in that maximum flood currents are always located, more or less, at the central part of the deep channel ( Figure 9, first column), while maximum ebb currents, when they are greater than 0.8 m/s, are located closer to the western shoal or over the western slope (Figure 10e,i). Only when maximum ebb currents are below 0.6 m/s, they move back to the middle of the channel (Figure 10a,m,q).

Conclusions
Barataria Pass is a tidal inlet that connects the Barataria Bay to the continental shelf. Previous investigation has introduced the tidal straining effect on density stratification during the same 25.6 h period along the same transect [22]. In this study, we conduct a numerical model simulation and illustrate that the lateral variations in the salinity and velocity fields are comparable to or even larger than the vertical variations within a diurnal tidal cycle.
The density distribution within any estuary is a result of both advective and mixing processes. In Barataria Pass, the turbulent mixing is closely related to the magnitude of ebb/flood current and the strength of the tidal bottom boundary layer. Characteristics of horizontal advection processes in the inlet are that maximum flood currents are located at the central part of the deep channel for a large part of the flood period. This differential advection [1], when acting upon the along-channel density gradient, produces a distinct density difference between the shoal and channel waters. In addition, the advection of Mississippi River water to the eastern channel during part of the flood period further enhances the density difference. On the contrary, maximum ebb currents swing between the western slope and central surface of the channel during the ebb. When maximum ebb flows are at the western slope, the differential advection mechanism does not work. When they go back to the channel center, the salinity contour lines are mostly horizontal due to weak vertical turbulence mixing. Thus, both situations are not favorable to produce an extreme density near the middle of the channel.
During flood period, when density distribution is high near the channel center and low at both shoals, the horizontal pressure gradient drives a lateral circulation with two counter-rotating cells and surface or near surface convergence. This result from the Barataria Pass is similar to that reported by Nunes and Simpson [1]. However, detailed analysis of momentum equations indicates that, in addition to the pressure gradient and vertical stress divergence, nonlinear advection and horizontal stress divergence are also important terms.
During ebb period, the lateral circulation is mostly eastward for the whole water column and persisting for almost the whole period. The surface divergence suggested by the differential advection mechanism is either non-existent or lasting for a very short period. The main momentum balance across most of the transect is between the along-channel advection of cross-channel momentum and pressure gradient. In addition, the sectional averaged lateral velocity magnitude during ebb is comparable to that during flood, which is different from the idealized numerical experiment [2].
Interactions among lateral circulation, along-channel tidal currents, and density stratification are complex processes. Lateral circulation can play a critical role in the dynamics of estuaries. It can act as a driving term for the estuarine exchange flow [2], and thus alter the along-channel momentum budget [38]. During flood tide, lateral circulation transported near-bed sediment out of the channel toward the shoal [37]. During ebb tide, it created a convergence in suspended sediment at the transition between channel and shoal [39]. Trapping by lateral circulation alters the lateral distribution of bed sediment. The idealized numerical experiment with constant eddy coefficients [2] demonstrated that the lateral circulation can be significantly different over a spring-neap cycle and density stratification can inhibit lateral circulation. Our study can be a starting point for further investigations of interactions among lateral circulation, estuarine circulation, and estuarine stratification in the partially stratified tidal inlet, the Barataria Pass. where (u, v) are the velocity components in the (x, y) directions, σ is the vertical coordinate.
In order to quantify along-and cross-channel momentum balance, we choose the along-channel direction (y ) to be aligned with the channel. It is positive when pointing into the estuary. The cross-channel direction (x ) is defined to be positive when pointing to the eastern boundary ( Figure A1). The relationship between (u , v ) and (u, v), (x , y ) and (x, y) are as following: where θ is the angle between the x -direction (cross-channel) and the x-direction. where ( , ) are the velocity components in the ( , ) directions, is the vertical coordinate. In order to quantify along-and cross-channel momentum balance, we choose the along-channel direction ( ′ ) to be aligned with the channel. It is positive when pointing into the estuary. The cross-channel direction ( ′ ) is defined to be positive when pointing to the eastern boundary ( Figure  A1). The relationship between ( ′ , ′ ) and ( , ), ( ′ , ′) and ( , ) are as following: where is the angle between the ′ -direction (cross-channel) and the -direction. To project the momentum equations into the cross-and along-channel directions, we treat each term in the momentum equations as a vector in the (x, y) direction and then apply the same decomposition as Equations (A3) and (A4). Thus terms in x -y coordinates can be calculated by corresponding terms in x-y coordinates as follows: To project the momentum equations into the cross-and along-channel directions, we treat each term in the momentum equations as a vector in the ( , ) direction and then apply the same decomposition as Equations (A3) and (A4). Thus terms in ′-′ coordinates can be calculated by corresponding terms in -coordinates as follows: where , 1' and 2' are shown in Figure A2. = ∂u 2 D ∂y sinθcos 2 θ + ∂uvD ∂y sin 2 θcosθ + ∂vuD ∂y sin 2 θcosθ + ∂v 2 D ∂y sin 3 θ + ∂u 2 D ∂x cos 3 θ + ∂uvD ∂x sinθcos 2 θ + ∂vuD ∂x sinθcos 2 θ + ∂v 2 D ∂x sin 2 θcosθ With the same method, ADVVX , ADVUY , and ADVVY are given as: = − ∂u 2 D ∂y sinθcos 2 θ − ∂uvD ∂y sin 2 θcosθ + ∂vuD ∂y cos 3 θ + ∂v 2 D ∂y sinθcos 2 θ + ∂u 2 D ∂x sin 2 θcosθ + ∂uvD ∂x sin 3 θ − ∂vuD ∂x sinθcos 2 θ − ∂v 2 D ∂x sin 2 θcosθ ADVUY = ∂u v D ∂x = − ∂u 2 D ∂y sin 2 θcosθ + ∂uvD ∂y sinθcos 2 θ − ∂vuD ∂y sin 3 θ + ∂v 2 D ∂y sin 2 θcosθ − ∂u 2 D ∂x sinθcos 2 θ + ∂uvD ∂x cos 3 θ − ∂vuD ∂x sin 2 θcosθ + ∂v 2 D ∂x sinθcos 2 θ