Freshwater Composition and Connectivity of the Connecticut River Plume During Ambient Flood Tides

The Connecticut River plume interacts with the strong tidal currents of the ambient receiving waters in eastern Long Island Sound. The plume formed during ambient flood tides is studied as an example of tidal river plumes entering into energetic ambient tidal environments in estuaries or continental shelves. Conservative passive freshwater tracers within a high-resolution nested hydrodynamic model are applied to determine how source waters from different parts of the tidal cycle contribute to plume composition and interact with bounding plume fronts. The connection to source waters can be cut off only under low-discharge conditions, when tides reverse surface flow through the mouth after max ambient flood. Upstream plume extent is limited because ambient tidal currents arrest the opposing plume propagation, as the tidal internal Froude number exceeds one. The downstream extent of the tidal plume always is within 20 km from the mouth, which is less than twice the ambient tidal excursion. Freshwaters in the river during the preceding ambient ebb are the oldest found in the new flood plume. Connectivity with source waters and plume fronts exhibits a strong upstream-to-downstream asymmetry. The arrested upstream front has high connectivity, as all freshwaters exiting the mouth immediately interact with this boundary. The downstream plume front has the lowest overall connectivity, as interaction is limited to the oldest waters since younger interior waters do not overtake this front. The offshore front and inshore boundary exhibit a downstream progression from younger to older waters and decreasing overall connectivity with source waters. Plume-averaged freshwater tracer concentrations and variances both exhibit an initial growth period followed by a longer decay period for the remainder of the tidal period. The plume-averaged tracer variance is increased by mouth inputs, decreased by entrainment, and destroyed by internal mixing. Peak entrainment velocities for younger waters are higher than values for older waters, indicating stronger entrainment closer to the mouth. Entrainment and mixing time scales (1–4 h at max ambient flood) are both shorter than half a tidal period, indicating entrainment and mixing are vigorous enough to rapidly diminish tracer variance within the plume.

The Connecticut River plume interacts with the strong tidal currents of the ambient receiving waters in eastern Long Island Sound. The plume formed during ambient flood tides is studied as an example of tidal river plumes entering into energetic ambient tidal environments in estuaries or continental shelves. Conservative passive freshwater tracers within a high-resolution nested hydrodynamic model are applied to determine how source waters from different parts of the tidal cycle contribute to plume composition and interact with bounding plume fronts. The connection to source waters can be cut off only under low-discharge conditions, when tides reverse surface flow through the mouth after max ambient flood. Upstream plume extent is limited because ambient tidal currents arrest the opposing plume propagation, as the tidal internal Froude number exceeds one. The downstream extent of the tidal plume always is within 20 km from the mouth, which is less than twice the ambient tidal excursion. Freshwaters in the river during the preceding ambient ebb are the oldest found in the new flood plume. Connectivity with source waters and plume fronts exhibits a strong upstreamto-downstream asymmetry. The arrested upstream front has high connectivity, as all freshwaters exiting the mouth immediately interact with this boundary. The downstream plume front has the lowest overall connectivity, as interaction is limited to the oldest waters since younger interior waters do not overtake this front. The offshore front and inshore boundary exhibit a downstream progression from younger to older waters and decreasing overall connectivity with source waters. Plume-averaged freshwater tracer concentrations and variances both exhibit an initial growth period followed by a longer decay period for the remainder of the tidal period. The plume-averaged tracer variance is increased by mouth inputs, decreased by entrainment, and destroyed by internal mixing. Peak entrainment velocities for younger waters are higher than values for older waters, indicating stronger entrainment closer to the mouth. Entrainment and mixing time scales (1-4 h at max ambient flood) are both shorter than half a tidal period, indicating entrainment and mixing are vigorous enough to rapidly diminish tracer variance within the plume.

INTRODUCTION
Rivers deliver terrestrial freshwater, nutrients, sediments, and contaminants to the marine environment (Meybeck, 2003). River plumes influence the physical dynamics, biogeochemistry, ecosystems and fisheries in coastal and open-ocean waters (e.g., Dittmar and Kattner, 2003;Hickey et al., 2010;Grimes, 2001). Due to these manifold effects on the marine environment, it is critically important to continue studying and improving the understanding of river plumes. The structure, mixing, and transport of river plumes has been studied extensively, as described in several research reviews (O'Donnell, 2010;Chant, 2011;Hetland and Hsu, 2013;Horner-Devine et al., 2015). Particular attention has been devoted to the main plume front, a region characterized by strong surface convergence, downwelling, and mixing of the plume with ambient waters (e.g., Garvine, 1974a;Garvine and Munk, 1974;O'Donnell et al., 1998;Marmorino and Trump, 2000;Orton and Jay, 2005). Flow within the plume can overtake the propagating plume front, thereby creating a large intake region of plume waters that interacts with the front (e.g., Garvine, 1974a;Mazzini and Chant, 2016). This connectivity with the front can extend all the way back to the river mouth, can be short-circuited by developing interior fronts, or can be limited by many other processes including tides, winds, ambient currents, and changing source conditions (e.g., McClimans, 1978;Garvine, 1984;O'Donnell, 2010;Chant, 2011;Huguenard et al., 2016;Cole et al., 2020). The degree of connectivity between river source waters and bounding fronts merits further study. The degree of source-front connectivity likely varies considerably for different plumes and forcing situations. The approach of Cole et al. (2020) provides guidance for plume connectivity studies. The transport, mixing, and trapping of pollutants, sediments, nutrients, and plankton within the plume is influenced by when and where source waters reach the plume front. Thus, identifying times and locations of particularly strong or weak sourcefront connectivity within plumes is important for understanding coastal biogeochemical distributions.
Tides modulate many river plumes. Strong tidal currents within rivers can reverse flow, preventing freshwater from exiting the mouth during river flood tides and intensifying freshwater outflow during ebbs. Two well-studied examples of ebb-pulsed plumes entering ambient environments with weaker tides are the Columbia River (e.g., Cudaback and Jay, 1996;Hickey et al., 1998;Horner-Devine et al., 2009;Nash et al., 2009;Kilcher and Nash, 2010;Akan et al., 2018) and the Merrimack River (e.g., MacDonald et al., 2007;Hetland and MacDonald, 2008;Chen et al., 2009;Kakoulaki et al., 2014;Cole et al., 2020). Ebb-pulse outflows also have been studied for smaller plumes (e.g., Luketina and Imberger, 1987;Pritchard and Huntley, 2006). Other tidal plumes are swept back and forth by strong ambient tidal currents in the receiving waters of the adjacent estuary, sea, or continental shelf (e.g., Spicer et al., 2021). The combination of tidally pulsed outflows through the mouth and swift ambient tidal currents along the coast can cut off plume waters from the river source and create a series of discrete tidal plume pulses (de Ruijter et al., 1997). Two formative examples are the Connecticut River plume (e.g., Garvine, 1974b;Garvine and Munk, 1974;Garvine, 1977;O'Donnell et al., 1998;Ackelson and O'Donnell, 2011;Jia and Whitney, 2019) and the Rhine region of freshwater influence (e.g., van Alphen et al., 1988;Simpson and Souza, 1995;de Ruijter et al., 1997;Hessner et al., 2001). Flow reversals and mixing by strong ambient tides also have been studied for smaller plumes (e.g., Bricker et al., 2006;Basdurak et al., 2020). Plume connectivity and overall behavior in these situations depends on the relative magnitudes of the subtidal mouth outflow velocity, plume propagation speed (scaled by the internal wave speed), and river and ambient tidal current amplitudes as well as the relative phasing of river and ambient tides. These factors influence which portions of the tidal cycle have a plume fed by source waters, a plume starved or cut off from the source, or no plume at all.
The present study investigates the composition of and connectivity within the Connecticut River plume, an important example of plumes influences by strong ambient tides. The Connecticut, named after a Native American word meaning 'on the long tidal river' (Trumbull, 1881), is the largest freshwater source flowing into Long Island Sound (LIS, Figure 1A). Average annual discharge (from 2009 to 2019) is 600 m 3 /s and average monthly discharge ranges from 1294 m 3 /s in April down to 304 m 3 /s in September (USGS 01193050, Middle Haddam, CT). The Connecticut River has a tidal salt-wedge estuary with a salinity intrusion that typically extends 5-15 km from the mouth (Meade, 1966;Garvine, 1975;Ralston et al., 2017). The navigational channel is bounded by two jetties (collectively referred to as the Saybrook Jetty) extending beyond the mouth ( Figure 1B). Planetary vorticity (f = 9.6 × 10 −5 s −1 ) does not exert a strong influence on the mouth outflow, indicated by a large Rossby number associated with the relatively fast outflow velocities O(0.1-1 m/s) through the 1 km wide mouth. Plume waters enter into the ambient macrotidal environment of eastern LIS, where tidal currents are O(1 m/s) with substantial spring-neap variability (O'Donnell et al., 2014). Ambient tides are predominantly rectilinear and shore-parallel near the coast, but have more ellipticity in the deeper waters farther offshore within the Sound (Bennett et al., 2010). Tides advect plume waters eastward during ambient ebbs and westward during flood tides (Garvine, 1974b;Garvine and Munk, 1974). The present study focuses on the plume that forms during ambient flood tides. During this tidal stage, the plume can pass over Long Sand Shoal, which is within 5 km from shore and stretches 10 km westward from the river mouth ( Figure 1B; Ellis and Gen, 1874;Poppe et al., 2000). Higher river discharge decreases salinities and increases the extent of the plume (Figures 2A,B; Garvine, 1974b). After high discharge events (e.g., tropical storm Irene flooding), coherent plumes formed during prior tidal cycles are evident farther offshore ( Figure 2C; Whitney et al., 2014). Connecticut River water persists in LIS after it is no longer part of a dynamically distinct tidal river plume. Some of its freshwater joins the LIS estuarine outflow and progresses through Block Island Sound onto the continental shelf within several days (Jia and Whitney, 2019). The rest of the Connecticut River water takes longer months-long routes through LIS, where it comprises most of the freshwater residing in the large estuary (Deignan-Schmidt and Jia and Whitney, 2019). Despite the many studies on the Connecticut River plume, the degree of connectivity between the front and the rest of the plume is still unknown. Early models imply connectivity extending from the front all the way back to the river source (Garvine, 1974a;1984). This high connectivity may be representative of frontal areas near the mouth, such as the upstream front. As the front travels farther from the mouth, however, increasingly large velocities within the plume would be necessary for water to travel from the mouth to the leading edge front before ambient tidal conditions switch from flood to ebb (or vice versa). Thus it is probable that less plume water interacts with the downstream front as the plume grows during a tidal cycle, resulting in relatively weak plume-front connectivity as seen for the Merrimack plume (Cole et al., 2020).
The main objectives of this study on the Connecticut River plume formed during ambient flood tidal conditions are: 1) determining the contributions of river source waters from different parts of the tidal cycle and 2) quantifying the degree and spatial distribution of connectivity of these source waters with the bounding plume fronts. A high-resolution numerical modeling approach is taken. The character of the tidal plume will be described to provide context for the subsequent analysis. Conservative passive freshwater tracers are employed as a key tool for studying plume composition and connectivity. A tracer variance budget also is applied to describe how the plume entrains waters and is mixed during the tidal cycle. The main case with low discharge during neap tides is compared to other cases with higher discharge and/or spring tides. Results are discussed relative to other plumes in weaker ambient tidal environments to emphasize how strong ambient tides can reshape the plume and patterns of source-front connectivity.

Model Setup
The Regional Ocean Modeling System (ROMS, Haidvogel et al., 2000Haidvogel et al., , 2008 is applied to model the Connecticut River plume and surrounding waters. ROMS is a free-surface, hydrostatic, primitive-equation model that evolves the governing momentum, mass, and tracer conservation equations in finitedifferenced form. The applied advection schemes are 3rd-order upstream for 3D momentum, 4th-order centered for 2D (depthaveraged) and vertical momentum, and the Wu and Zho (2010) 3rd-order scheme for salinity, temperature, and passive tracers. Vertical turbulent viscosity and diffusivity are parameterized with the generic length scale method k-epsilon closure scheme. Note that salinity is calculated with the Practical Salinity Scale, thus reported salinities are dimensionless.
The outer (mother) grid domain includes the tidal extent of the Connecticut River (and other rivers), LIS, Block Island Sound, and a large portion of the adjacent continental shelf ( Figure 1A). Full details of the outer model configuration are included in Jia and Whitney (2019). The horizontal resolution is 500 m within LIS and vertical resolution is supplied by 30 sigma levels evenly distributed through the water column. Outer model bathymetry is derived from the NOAA 3 arc-second United States Coastal Relief Model. The outer model is forced along the open ocean boundaries with eight semi-diurnal and diurnal tidal coefficients from the OSU TOPEX/Poseidon Global Inverse Solution (Egbert and Erofeeva, 2002) and subtidal velocities, temperatures, and salinities from the Hybrid Coordinate Ocean Model (Chassignet et al., 2008). Chapman (1985) and Flather (1976) boundary conditions are applied for surface elevation and 2D velocities along the open boundaries. Radiation with nudging open boundary conditions are applied for 3D velocities, temperature, and salinity. Surface winds, humidity, air temperature, surface pressure, net shortwave radiation, and downward longwave radiation are supplied by the North American Regional Reanalysis (NARR, Mesinger et al., 2006). Bulk flux formulae (Fairall et al., 2003) are applied within ROMS to compute surface fluxes of momentum and sensible, latent, and longwave heat. Jia and Whitney (2019) describe how the NARR temperatures were modified to better align with NOAA buoy observations. Daily river discharge for the Connecticut and 78 other rivers are derived from USGS streamgage records. Jia and Whitney (2019) validated tidal performance with tidal harmonic analysis at NOAA coastal stations along LIS and with depth-averaged tidal current observations throughout LIS (Bennett et al., 2010).
This study adds a nested intermediate grid with 100 m resolution and an inner grid with 20 m horizontal resolution FIGURE 2 | Connecticut River plume observations: (A) surface salinities at high slack (near spring tide) for moderate discharge (720 m 3 /s) conditions (from Figure 15 in Garvine, 1974b), (B) surface salinities at high slack (near neap tide) for high discharge (1580 m 3 /s) conditions (from Figure 13b in Garvine, 1974b), and (C) a true-color Landsat 5 satellite image coincident with max ambient flood tide (near spring tide) during high discharge (1633 m 3 /s) and sediment load from tropical storm Irene (from Figure 2 in Whitney et al., 2014). Discharges in parentheses are daily data from USGS station 01184000 (Thompsonville, CT) multiplied by 1.086 to represent total Connecticut River discharge. The salinity contour interval is 5, as in Garvine (1974b).
covering the lower Connecticut River estuary and the tidal plume region ( Figure 1A). The nested runs are initialized with a mother grid simulation that spans at least 2.5 years prior to nesting. The nested runs each are evolved for only 1.5 days, these short runs are sufficient to capture the intratidal evolution of the plume. Light wind periods are selected and then wind stress is set to zero to remove all immediate wind influence during the short nested runs. Four cases are analyzed that are combinations of low and high discharge and neap and spring tides. Connecticut River discharge (Q R ) is set at 250 and 1500 m 3 /s for low and high discharge cases, respectively. The dates and details of these cases are summarized in Table 1. The low discharge and neap tide case is analyzed in greatest detail as it produces the most compact plume and pairs with new observational efforts. Instantaneous fields (in ROMS 'history' files) output at 20-min intervals are analyzed. Results are mapped on an x-y grid with the origin at center channel of the river mouth, the positive x axis is oriented along the coast in the direction of LIS flood tides (approximately west-southwestward), and the positive y axis is pointed down-river and offshore (approximately southsoutheastward) (Figure 1).

Freshwater Tracers and Plume Characteristics
Freshwaters from the Connecticut River are tracked with conservative passive tracers using the tracer routines included within ROMS. The passive tracers are evolved in the Eulerian framework and are akin to continuous dye fields. Tracer concentrations are normalized so that they represent the fraction of the tracked freshwater. Individual freshwater tracers are released at an hourly interval spanning a tidal cycle. Each tracer is imposed as an hour-long release at the up-river boundary of the inner grid (4 km up-river from the mouth) with the same volume flux and freshwater fraction of the Connecticut River waters flowing through the boundary over the release time interval. There are no other interior or boundary tracer sources. There are no tracer sinks; therefore, concentrations of these conservative passive tracers decrease only through dilution. The hourly tracers are summed to create a combined tracer that tracks all the recently discharged Connecticut River freshwater within the plume. Plume waters are identified as including all locations (both horizontal and vertical) outside of the Connecticut River mouth where the combined freshwater tracer concentration exceeds a low threshold value (C 0 ) set to 10 −6 for standard analysis. In practice, the plume boundaries change little with an order of magnitude increase or decrease of C 0 because of typically sharp concentration gradients at plume fronts. Plume volume (V p ) and area at the surface (A p ) are calculated within the plume boundaries. The average plume depth (H p ) can be calculated as V p /A p . Plume length (L p ) is calculated as the farthest downstream extent of the tracer plume. The plume width (W p ) is reported as the farthest offshore extent of the main plume marked by a sharp surface front and does not include more diffuse tracer waters farther offshore that are included within V p . Plume perimeter (P p ) is calculated as the total length of the plume boundary at the surface. The composition of plume boundary points in terms of individual freshwater tracers is used to quantify connectivity between the plume fronts and source waters released at different times of the tidal cycle.
Plume volume changes due to the volume flux through the mouth (Q m , with positive into the plume) and entrainment volume flux into the plume (Q e , positive into the plume): Q m is calculated as the section integral of the v velocity component through the mouth where there is tracer (with higher concentration than C 0 ), operationally v is multiplied by a binary coefficient that is set to one if the tracer is present above the concentration threshold and zero otherwise. The instantaneous Q m can either be into or out of the plume, but represents a net source to the plume overall. The entrainment flux adds water with zero tracer concentration to the plume. V p , its rate of change dV p /dt, and Q m are calculated directly from model output fields and Q e is solved for as the unknown in (1). The plume-averaged entrainment velocity (w e ) is calculated by dividing Q e by the plume surface area: Note that, despite the suggestive notation, entrainment does not have to occur exclusively through vertical velocities; entrainment velocities are locally normal to plume boundaries. The analysis is completed with the combined tracer to represent conditions for the entire plume. Individual hourly tracer fields identify the plume regions associated with source waters distinguished by river release time. Thus, corresponding analysis of individual hourly tracers indicates the entrainment within different parts of the plume.

Freshwater Tracer Variance Budgets
Analysis of salinity variance has been applied to studying mixing in estuaries and coastal waters (e.g., Burchard and Rennau, 2008;Li et al., 2018;MacCready et al., 2018;Wang and Geyer, 2018;Burchard, 2020;Warner et al., 2020). Analogous tracer variance analysis is applied to the freshwater tracers described above. The first step is decomposing tracer concentration (C) into a spatial-mean value (C) and deviations from the mean (C ) such that C = C + C . The tracer variance is defined as C − C 2 which equals C 2 . The spatial mean may be taken over a fixed volume (e.g., an estuary or part of a continental shelf) or a time-varying volume tracking a water mass (e.g., a river plume). Most prior published examples analyze fixed volumes extending throughout the water column, whereas the present analysis considers expanding control volumes encompassing river plumes. In this application, the plume control volume (V p ) is the region with freshwater tracer concentrations above the threshold value (C 0 = 10 −6 , as described above). The control volume only extends throughout the water column where the plume interacts with the bottom, otherwise there is a subsurface interface between plume and ambient waters. The plumeaveraged concentration (C) is calculated as: Following MacCready et al. (2018), the volume-integrated tracer variance budget is: where, u n is the outward normal velocity through the bounding area (A b ) of the control volume, K is the eddy diffusivity, and ∇C is the concentration gradient. K is written as a scalar in (4) for simplicity (as in MacCready et al., 2018), but a more detailed representation involves a diagonal diffusivity tensor with different horizontal and vertical eddy diffusivities (e.g., Li et al., 2018). Alternately, horizontal diffusion is sometimes considered negligible and only vertical diffusion is included (e.g., Warner et al., 2020). The final integrand in (4) represents the tracer mixing per unit volume; it is symbolized with χ in some prior studies (e.g., Wang and Geyer, 2018;Burchard, 2020). The term on the left-hand side of (4) can change both because of tracer variance changes and plume volume changes. The first term on the right-hand side of (4) is the advective transport of variance through the control volume boundaries. Note that (4), does not include a term representing diffusive transport of tracer variance across the boundary between the plume and ambient waters. This diffusive transport is implicitly considered negligible relative to advective transport in (4), but the diffusive transport of tracer variance is considered in other studies (e.g., Burchard, 2020;Wang et al., 2021). The final term in (4) is the dissipation of variance by mixing inside the control volume; it can be symbolized as -M for 'mixing'; M is positive definite. The M notation follows MacCready et al. (2018) except in that study M is a tidal-averaged quantity, whereas here it is tidally varying. In the present application, M is calculated (as described below) such that it includes vertical and horizontal physical mixing as well as numerical diffusion arising from discretization errors associated with the advection scheme (Burchard and Rennau, 2008). For the plume control volume, transport through the plume boundaries can be partitioned into contributions through the river mouth and by plume entrainment. The mouth will tend to have higher concentrations than the plume and therefore have a strong tracer variance signature. With an outflow through the mouth cross-sectional area (A m ) and into the tracer plume, the velocity through the mouth (v m ) is positive and adds both volume and tracer variance to the plume. Entrained water comes from outside the tracer plume and its tracer concentration equals the low threshold value (C 0 ). The corresponding variance of entrained water is C 2 e = C 0 − C 2 and the entrainment tracer variance flux is C 2 e Q e , a positive flux. Substituting this information into (4), the volume-integrated tracer budget for the tracer plume is: Note that, C 0 C since an exceedingly small C 0 has been selected (i.e., the entrained water has nearly zero tracer concentration) and C 2 e ≈ C 2 . Therefore, the entrainment tracer variance flux is approximately C 2 Q e . This approximation is used when describing results and discussing generalization of tracer variance budgets for other plumes. In summary, (5) indicates that river inputs and entrainment both increase volume-integrated variance, while mixing within the plume reduces it. Analyzing the volume-averaged (rather than volumeintegrated) budget provides a different perspective that isolates changes in tracer variance and facilitates intercomparison of results for different plumes. A first step is defining the volume-averaged variance: With this definition, the term on the left-hand side of (5) can be expressed as the time derivative of the product of the volumeaveraged variance and plume volume (7, the intermediate step). Then the product rule for derivatives can be applied and the plume volume budget (1) can be substituted to separately show the contributions associated with changing volume-averaged tracer variance, mouth volume inputs, and entrained volume flux (7).
It also is useful to define the flux-weighted average tracer variance at the mouth: With this definition, the first term on the right-hand side of (5) can be expressed as: Dividing the volume-integrated tracer variance budget (5) by V p and substituting in expressions from (6), (7), (8), and (9) yields the volume-averaged tracer variance budget: It can be seen from this balance that a volume source from the mouth tends to increase plume-averaged tracer variance as long as there is higher variance at the mouth (C 2 m /C 2 >1), as would be the case for high-concentration source waters. Since C 2 e is approximately C 2 for small C 0 , entrainment will tend to increase plume-averaged tracer variance when the variance is smaller than plume-averaged concentration squared (C 2 /C 2 <1); otherwise C 2 /C 2 >1 and the entrainment term is negative and tends to decrease C 2 . The mixing term is simply divided by V p and remains negative as it represents the internal destruction of tracer variance. The relationship between (5) and (10) is similar to the relationship between heat and temperature budgets in the sense that any net advective inflow will increase heat and volumeintegrated tracer variance yet the same inflow can either increase or decrease temperature and volume-averaged tracer variance depending on whether the inflow has higher or lower scalar values than the control-volume average. Both perspectives are useful, but the volume-averaged budget (10) is analyzed in this study to focus on tracer variance changes.
The model provides the information to directly calculate all terms in tracer variance balances (5) and (10). The computational approach in this study applies model results to directly calculate all but the final term in (5) and (10). Then M is solved as the unknown in (5) and divided by V p for the final term in (10). This calculation method for M ensures that the tracer variance budget closes exactly and that M represents the total mixing in the model, which includes physical mixing and numerical mixing from advection scheme discretization errors. Analyzing the combined tracer fields provides the tracer variance budget associated with the entire active plume. Analyzing individual hourly tracers indicates the tracer variance budgets of plume regions associated with different source waters, distinguished by river release time.

Salinity and Velocity
For the standard case with low discharge and neap tides, the along-estuary currents in LIS (2 km east of the river mouth within 4 km from shore) indicate the ambient LIS tidal amplitude (U Ta ) is 0.75 m/s ( Figure 3A). With the time origin at the onset of ambient flood tide the time-varying ambient velocity (u a ) can be approximated as: where, t is time, T is the tidal period (12.42 h for the dominant semidiurnal lunar tides), and u a is positive during flood tide. The corresponding ambient tidal excursion (L Ta = U Ta T/π) is 10.7 km. The surface flow is outward through the mouth for much of the tidal cycle ( Figure 3A). The tidal-averaged surface velocity at the mouth (V D , following the notation in de Ruijter et al., 1997) is 0.25 m/s and the mouth tidal amplitude (V Tm ) is 0.70 m/s. The V D /V Tm ratio is 0.36 and meets the V D /V Tm < 1 'pinching off criterion' for flow reversal at the mouth and shutoff of plume source waters (de Ruijter et al., 1997). The surface velocity at the mouth (v m ) can be approximated as the combination of the steady flow and sinusoidal tidal variations: where, t lag it the lag time of the mouth tides relative to ambient tides, such that a positive lag indicates ambient flood leads the mouth flood tide. Note that v m is negative while the flood tide reverses the flow to inward through the mouth. The ambient flood tides lead mouth flood tides by approximately 1 4 of the tidal cycle for this model run ( Figure 3A, t lag = 3.1 h). As found in de Ruijter et al. (1997), trigonometric investigation of (12) indicates the mouth flow reversal time span (t rev ) is: For this low-discharge neap-tide case, t rev is 4.8 h; indicating the plume is cut off from source waters for at least 38% of the tidal cycle. Since this study focuses on plumes during ambient flood tides, it is useful to express the time interval when there is a mouth reversal during ambient flood tides. With equations (11), (12), and (13), this mouth-reversal during ambient-flood time interval can be expressed as: ) For this model run, flow at the mouth is reversed throughout late ambient flood beginning at t = (T-t rev )/2 (t = 3.8 h). This indicates the plume is fed by source waters during early and max ambient flood tides, but is cut off from the source during late flood.
Mouth surface-averaged salinity falls from 14.0 to 8.8 throughout mouth ebb tides, then increases to 29.6 near the end of mouth flood (Figure 3B). Bottom-averaged salinities at the mouth remain salty throughout the tidal cycle ( Near the start of ambient flood tide, the plume has expanded 6 km outward from the mouth with a western (downstream, relative to ambient flood tidal flow) bias in lower salinities, higher velocities, and sharper boundaries ( Figure 4A). The flow of lowsalinity source waters is bifurcated by Saybrook Jetty outside the mouth and bifurcation effects are evident out to the plume boundary. The lowest salinities on the eastern (upstream) side are within 2 km from the mouth, but more diffuse plume waters associated with the prior ambient ebb extend to 5 km eastward along the coast. The sharp boundary on the downstream side bends around farther offshore and the highest velocities (0.9 m/s) and lowest salinities (20) are near the middle of this front. The Froude number in the downstream frontal zone (Fr d ) exceeds two, for S d = 8.5 and h d = 2 m, indicating supercritical flow. Two hours later at max ambient flood (Figure 4B), the plume stretches 8 km downstream (L p > L Ta /2), extends 4.5 km offshore, has higher salinities than before, and remains connected to source waters. Velocities within the plume increase downstream from 0.7 to 1.1 m/s and are predominantly westward except the offshore flow near the mouth. The plume is supercritical, the plume internal wave speed (c p ) is 0.27 m/s and the offshore plume flow. The propagation of the upstream front has been arrested by the opposing tidal flow because the Froude number associated with ambient tides (Fr Ta = U Ta /c p ) exceeds one. Plume salinities increase downstream along the offshore boundary, but remain lower than inshore. The downstream front is clearly defined and the plume overtakes the slower ambient waters. Two hours later during late ambient flood (Figure 4C), the plume has extended to 14 km from the mouth (farther than L Ta ) and is bending toward the coast. Flow is reversed at the mouth and velocities increase along the plume to 0.9 m/s near the downstream front, approximately double the weakening ambient tidal flow. The plume, now cut off from source waters, has higher salinities than before and the lowest salinities are in a partially disconnected area (7-14 km downstream), where the minimum salinity is 24.2. The plume remains supercritical, with Fr d exceeding three for S d = 4.5 and h d = 2 m. Two hours later as ambient tides switch to ebb (Figure 4D), the freshwater supply still is cut off and the remnants of the plume are evident. In the plume area, the salinity difference has diminished ( S d < 2.0) and velocities have reversed with the tides.
Velocities and salinities at max ambient flood tide are examined in further detail in Figures 5, 6, respectively. With a closer view of surface conditions (Figures 5A, 6A), the main plume and ambient flow features can still be seen and a secondary more diffuse plume farther offshore (from y = 4.5 to 8 km) with a convergent downstream front is evident. The mouth velocity field (Figure 6B) indicates the outflow is within 1-2 m from the surface where salinities are lowest ( Figure 5B). There is stratification throughout the channel water column. The nearfield section shows the sharp offshore salinity front ( Figure 5C) which has strong convergence of the cross-shore flow ( Figure 6C). The surface-intensified freshwater layer increases thickness from the coast out to the offshore front, contrary to large-scale buoyancy-driven plumes (e.g., Yankovsky and Chapman, 1997). The salinity front is farther offshore and not as sharp in the farfield section ( Figure 5D), but there is still strong flow convergence ( Figure 6D). The freshwater layer plunges deeper just inshore of Long Sand Shoal, likely due to plume interaction with bathymetry along the shoal.

Freshwater Tracers
Surface concentrations of freshwater tracers at the mouth ( Figure 3C) show the pulses associated with hourly tracer releases. Peak concentrations range from 0.20 to 0.65; the hourly tracers with the highest peak concentrations (>0.60) are mixed less within the river estuary before reaching the mouth. The tracers that first appear at the mouth during early mouth ebb (late ambient ebb) were released during late river flood. Concentrations peak at the mouth about an hour after peaking at the release point 4 km upriver. The final tracer appearing at the mouth was released during late river ebb (early ambient flood). Despite all being imposed as hour-long pulses, the time interval of elevated concentrations at the mouth increases substantially for each subsequent tracer. The combined tracer increases during late ambient ebb, remains high through max ambient flood, and drops through late ebb.
By max ambient flood the plume is composed of a variety of source waters, mostly originally released in the river during ambient ebb tide ( Figure 7A). Average freshwater ages are calculated as the time interval between the average tracer release time at each location (determined by the concentration-weighted average of all tracers) and the mapped time (3 h into ambient flood). Freshwater age increases along the plume, with the oldest water (with earliest release time) at the downstream tip. The offshore secondary plume has the oldest waters, indicating it is a relic of the prior ebb plume. Overall, the upstream plume boundary has the youngest freshwater, the downstream boundary has the oldest, and the offshore and inshore boundaries have a downstream progression from younger to older waters. The mouth section (Figure 7B) indicates the outflowing surface layer has the youngest waters. The nearfield section ( Figure 7C) shows age decreases toward the offshore front with some vertical variability, whereas the farfield section ( Figure 7D) shows less spatial variability in water age.
Higher combined tracer concentrations correspond to lower salinity, so surface concentrations highlight the same plume patterns (Figure 8A) as the surface salinity field (Figure 5A). The texture of plume composition is better seen with maps of the hourly tracers. By max ambient flood, the two tracers introduced during late river flood and exiting the mouth in late ambient ebb ( Figure 8B) either are in the secondary offshore plume (a relic of the prior ebb plume) or have been incorporated in the new flood plume. Concentrations increase from <0.001 at the mouth to 0.05 toward the downstream front. These source waters are mostly disconnected from the mouth at this time, but some tracer remains trapped in coves along the river. The two early river ebb tracers, exiting the mouth from late ambient ebb into early ambient flood, (Figure 8C) are distributed throughout the plume and tend to increase in concentration from 0.02 along the coast to 0.12 toward the offshore front. The tracers introduced later in river ebb ( Figure 8D) enter the plume later in ambient flood. These tracers are present in highest concentrations (>0.30) near the mouth, particularly along the arrested upstream front, and extend only halfway downstream into the plume. Overall, the downstream half of the plume only includes source waters introduced from late river flood to early river ebb, while the upstream half is composed of a greater variety of source waters with the late river ebb waters at highest concentrations. It is evident that all tracers interact with some portions of the plume boundaries, but the downstream front only sees older waters.
The plume perimeter grows linearly over the ambient flood tidal cycle (Figure 9A). The perimeter of each individual tracer field also approximately grows linearly with time, but the growth rate decreases for later tracers. How much each tracer interacts with plume fronts is quantified by the percent of the plume boundary with each hourly tracer present ( Figure 9B). Each tracer goes through an initial period (lasting 1.5 to 3 h) of increasing interaction with plume boundaries. Frontal connectivity is highest for the oldest waters (present when the early ambient flood plume forms), with prevalence at the front exceeding 80%. Frontal connectivity is lowest for the youngest waters (entering as the ambient flood progresses), which interact with only 20% of the plume front. Overall, frontal connectivity is proportional to water age.
Tracking the perimeter allows for analysis of downstream front advancement relative to ambient and interior plume waters. At max ambient flood, the downstream front advance velocity (u f ) is 0.96 m/s, which is 0.21 m/s faster than U T a and is 94% of U Ta +c p . The interior water (the leading edge of the −0.3 h tracer release, Figure 8D) advance velocity (u i ) is 1.05 m/s, which is 0.09 m/s faster than the downstream plume front is advancing. Thus, the plume front is overtaking ambient water while interior waters are more gradually closing the gap at relative velocity u iu f . The interior waters, however, are not traveling quickly enough to reach the downstream plume front by the end of ambient flood.
Plume volume (Figure 10A) grows at an increasing rate over time throughout ambient flood tide, reflecting a combination of plume lengthening, widening, and deepening. As for perimeters, the volumes and increase rates of individual tracers are higher for earlier tracer releases than later ones. Following (1), plume volumes increase only due to entering source waters and entraining ambient waters. After source inputs are shutoff (when Q m = 0), plume volume exponentially increases due only to entrainment with an e-folding time scale (T e ) equal to the average plume depth divided by the entrainment velocity: Since H p = V p /A p and w e = Q e /A p , T e also can be expressed as V p /Q e . Note that T e can be time variable, as both the entrainment velocity and plume depth can vary. Plume-averaged entrainment velocities, calculated with (2), are 0.7-2.0 mm/s with variability during the tidal cycle ( Figure 10B). Most of the tracers have entrainment velocities within about the same range (0.7-2.3 mm/s), indicating a similar entrainment efficiency to each other and the plume as a whole. The younger waters included in the latest two tracer releases have notably higher peak entrainment velocities (3.3-3.8 mm/s). Since these younger waters are concentrated in the upstream half of the plume (Figure 8D), the result indicates more vigorous entrainment in the nearfield plume. The linear increase in average plume depth over time (not shown) results in a linear increase in plume entrainment time scale (T e ) from 1 to 6 h from the onset to the end of ambient flood tides, so plume volume exponentially increases at a slower rate over time.  Figure 3). Average freshwater age can be calculated by subtracting the average release times from 3 h (the mapped time). Shaded contouring only includes locations where the combined tracer concentration is above the C 0 = 10 -6 threshold.
Volume-averaged tracer concentrations within the plume (C, Figure 11A), calculated with (3), initially increase (for 0.5-1 h) due to source water supply and then decrease due to dilution by continued entrainment of ambient waters. Peak plume-averaged concentrations range from 0.035 to 0.085, with a 0.060 peak for the combined tracer. Concentrations exponentially decay at the same rate (set by T e ) at which plume volume increases. Note that the decay in concentrations is only due to dilution, as the conservative nature of these tracers precludes imposed exponential tracer loss or other interior tracer sinks. By the end of ambient flood, average tracer concentrations decrease to <0.004, which are <5% of their peak average concentrations and <0.4% of the undiluted unit concentration. Thus, entrainment leads to a voluminous but highly diluted plume by the end of ambient flood tide.

Freshwater Tracer Variance Budgets
The freshwater tracer variance budgets assess how source water inputs, entrainment, and mixing affect plume evolution over time. The analysis is completed over the expanding plume control volume, defined as the region where the combined freshwater tracer concentrations exceed the threshold value (C 0 = 10 −6 ). The plume control volume boundaries are shown in Figures 7, 8A. At the mouth ( Figure 7B) and in the nearfield (Figure 7C), the control volume extends throughout the water column. Farther downstream (Figure 7D), the plume does not reach the bottom in all locations. Volume-averaged tracer variances (C 2 , Figure 11B) follow patterns similar to volume-averaged tracer concentrations ( Figure 11A), with an initial increase pattern and a long decay period through the rest of the tidal cycle. Peak tracer variances range from 0.002 to 0.021 and are two to three times larger than the corresponding squared peak concentration values (i.e., C 2 /C 2 >1). By the end of ambient flood tide, tracer variances decrease to <10 −4 , which are <1% their peak concentrations. The volume-averaged tracer variance budgets (10) track how C 2 changes. The rate of change for the combined tracer variance ( Figure 12A) is positive with a peak during the initial period of concentration and variance increase (Figure 11) and then is negative with smaller magnitude during the longer decay period. Mouth input (Figure 12B) is the single largest term and the only one that acts to increase C 2 , as it adds high concentration water with higher tracer variance than the plume average. Entrainment ( Figure 12C) brings in low concentration water (at the C 0 threshold value) and tends to decrease C 2 . The entrainment term in (10) is negative because C 2 /C 2 e >1. As described earlier, an approximately equivalent condition is C 2 /C 2 >1 because C 2 e is approximately C 2 for small C 0 . The C 2 /C 2 >1 condition is met for this plume, but this does not have to be the case in general. The internal mixing term ( Figure 12D) always acts to destroy tracer variance; the term grows smaller while C 2 decreases and V p increases. This mixing within the plume control volume reflects vertical and/or horizontal mixing that tends to homogenize tracer concentrations within the plume. Dividing the tracer variance by the absolute value of the mixing term in (10) yields a mixing time scale: T mix is 1.3 h and T e is 2.2 h at max ambient flood, indicating mixing ( Figure 12D) is about twice as powerful as entrainment ( Figure 12C) in reducing tracer variance at this stage of plume evolution. Sensitivity testing indicates tracer variance budgets at max ambient flood are similar when higher concentration threshold values (up to C 0 = 10 −3 ) are applied, but the term magnitudes approximately double (increasing by a factor of 2.1-2.2) from the original C 0 = 10 −6 to the C 0 = 10 −3 results. For all tested C 0 , mixing is approximately twice (1.8 to 1.9 times) the entrainment at max ambient flood. During early ambient flood (1 h after the onset of flood tide), however, sensitivity testing indicates the mixing to entrainment ratio increases from 0.6 to 1.5 from the C 0 = 10 −6 to the C 0 = 10 −3 results. The volume-averaged variance budgets for the hourly tracer are broadly similar (Figure 12). The entrainment and mixing terms contribute about equally to reducing variance for the combined tracer and most hourly tracers. Entrainment is more than double the mixing term for the two latest tracer releases; this is consistent with the higher associated entrainment velocities ( Figure 10B, described above) for younger waters in the nearfield plume. Overall, tracer variance inputs from the mouth exceed reductions via entrainment and mixing only during an initial 0.5-to 1-h period; tracer variance is eroded throughout the rest of the ambient flood tide.

Other Discharge and Tide Conditions
The standard case is representative of summer low discharge conditions during neap tides. By comparison, the low-discharge spring-tide case increases ambient tidal currents by 33% ( Table 1). Tides through the river mouth increase by a smaller amount (7%) and t lag is shorter ( Table 1). Surface flow reverses an hour earlier during spring tides; therefore, the plume is cut off from source waters prior to max ambient flood. Consistent with stronger river tides, the bulk vertical salinity difference at the mouth is smaller ( S m = 16.0) and the Froude number at the mouth is larger (Fr m = 2.0) ( Table 1). The plume is somewhat longer (10 km length), narrower (3.5 km width), and has a lower downstream salinity anomaly ( S d = 4.0) during spring tides (Figure 13 and Table 1), but the dimensions and salinity anomaly are within 30% of the neap plume characteristics. The springtide plume is supercritical with Fr d exceeding five. Comparing the high-discharge neap-tide case to the standard case indicates how river mouth and plume conditions change with six times stronger discharge ( Table 1). Even though ambient tides are the same, river tides are weaker with a shorter t lag ( Table 1). The tidalaveraged flow exceeds the river tidal amplitude (V D /V Tm < 1), so currents do not reverse and the 'pinching off criterion' (de Ruijter et al., 1997) is not met. Consequently, the plume is fed by source waters throughout the tidal cycle. Completely freshwater reaches the mouth ( Table 1) and low-salinity water occupies the entire water column in the channel (not shown). The fresher and thicker mouth outflow ( S m = 23.0, h m = 10 m) has a faster c m (1.32 m/s), but source conditions are still at least critical at max ambient flood ( Table 1). The main plume is somewhat longer (9.5 km length) and wider (5.5 km width), but the much higher salinity anomaly ( S d = 14.0) and the more apparent offshore relic plume are the most conspicuous changes relative to the low-discharge plume (Figure 13 and Table 1). The stronger salinity anomaly increases c p to 0.46 m/s and decreases Fr d to 2.4, but the plume still is supercritical at max ambient flood. Even with a stronger c p , ambient tides still arrest the upstream front (Fr Ta > 1). Comparing the high-discharge runs with neap and spring tides indicates a qualitatively similar spring-neap tidal dependence as shown for the two low-discharge runs.
Intercomparison of the four model runs at max ambient flood (Figure 13 and Table 1) indicates the mouth outflow and plume velocities always become supercritical during part of the tidal cycle. Under low-discharge conditions, the tidal plume is fed by source waters until max ambient flood, when flow reversal at the mouth cuts off the source water supply. During high-discharge, the plume maintains a connection with source waters because the mouth surface flow always is outward. The plume at max ambient flood is always faster than ambient tidal currents (i.e., u d > U Ta ), and longer than the time-integral of the tidal flow (i.e., L p > L Ta /2). Plume salinities decrease with weaker tides and particularly with higher discharge. Plume dimensions, however, are less sensitive to changing tidal and discharge forcing. The upstream extent is limited because ambient tidal currents arrest the opposing plume propagation, as Fr Ta > 1 in all runs. The downstream extent of the tidal plume always is within 20 km from the mouth, which is less than twice the ambient tidal excursion.
The hourly freshwater tracers released within each run indicate the composition and connectivity within the plume under various forcing conditions. In all cases, water age increases downstream and offshore to a lesser degree (Figure 14). The downstream front has older waters, the upstream front has younger waters, and the main offshore front and more gradual inshore boundary span young to old waters. The oldest waters are found in the relic offshore plume (Figure 14). At max ambient flood, the plumes are chiefly composed of source waters released in the river during the preceding ambient ebb. The oldest water ages are detected during neap tides, particularly within the relic plume that persists under high-discharge conditions. The low-discharge spring-tide plume has the narrowest age range because the oldest waters have been mixed away and the youngest waters are cut off from the plume by the mouth flow reversal (Figure 14). The plume front advances downstream faster than ambient tidal currents in all runs, but the u f advancement velocity is somewhat less than (64-94%) U Ta +c p . Interior waters in all runs, assessed with the hourly tracer closest to but not at the downstream front, approach the downstream front with a relative velocity (0.09 < u i < 0.19 m/s) at max ambient flood ( Table 1). The interior waters, however, are not traveling quickly enough to reach the downstream plume front by the end of ambient flood. Thus, the downstream front always has the lowest overall connectivity with source waters, as interaction is limited to waters that exited the mouth by early ambient ebb. In stark contrast, the arrested upstream front has high connectivity, as all freshwaters exiting the mouth immediately interact with this boundary. Thus, energetic ambient tides lead to a strong upstream-to-downstream asymmetry in source-front connectivity.
Analysis of the plume volume budgets, identified with the combined tracer fields, indicates how the entrainment velocity and the associated e-folding time scale vary among runs. The w e values for both spring-tide cases (1.0-1.2 mm/s) are approximately the same and are about twice the w e values for neap-tide cases (2.0 mm/s) ( Table 1). The w e tidal dependence, and lack of sensitivity to discharge, point to entrainment driven by tidal shear. The entrainment time scale (T e ) at max ambient flood is about the same (2 h) for all but the high-discharge neaptide plume, which has double the T e . For the combined tracer at max ambient flood, C ≥ 0.026 in both high discharge runs; these values are at least double the concentration in the standard low-discharge neap-tide run (C = 0.013) and several times greater than in the low-discharge spring-tide run (C = 0.004) ( Table 1). Note that the higher concentration values in the highdischarge cases correspond to lower plume salinities. Both high discharge runs also have much higher variance at max ambient flood (C 2 ≥ 0.0069) than the low-discharge neap-tide run (C 2 = 0.0019) and the low-discharge spring-tide run (C 2 = 0.0002) ( Table 1). The ratio C 2 /C 2 at max ambient flood ranges from 9.5 to 12.5 ( Table 1) and C 2 /C 2 >1 during the entire ambient flood tide in all cases. Tracer variance budgets are qualitatively similar, but term magnitudes vary by orders of magnitude among runs. In all runs, C 2 is increased only by mouth inputs and is decreased both by entrainment and mixing ( Table 1). Entrainment decreases tracer variance because C 2 e is approximately C 2 and C 2 /C 2 >1 for all tidal and discharge conditions explored; indicating that the plume remains in the same regime in this regard. By max ambient flood, mixing is about twice as strong as entrainment and T mix (1-2 h) is shorter than T e . Both time scales are less than half the tidal cycle, indicating mixing and entrainment are vigorous enough to rapidly diminish tracer variance within the plume.

DISCUSSION
This study investigates the Connecticut River plume formed during ambient flood tides. Observations (e.g., Garvine, 1974b; Figure 2) and model results indicate strong plume salinity variability with forcing conditions, but the downstream extent of the tidal plume always is within 20 km from the mouth. The plume length is strikingly short in comparison to plumes exiting other estuaries with similar river discharge ranges. For comparison, the buoyant outflow from the Delaware Bay can stretch over 100 km along the coast during buoyancy driven conditions (e.g., Whitney and Garvine, 2006). The key difference is the Connecticut River plume forms in a much stronger ambient tidal environment. Tidal advection during flood increases downstream transport and extent, but plume length ultimately is limited by ebb currents strong enough to arrest and then reverse the plume. Simple estimates are useful for characterizing plume length generally for other tidally reversing plumes. The frontal advancement velocity is approximated as the internal wave speed superimposed on the ambient tidal velocity: u f ≈c p +u a , where u a can be expressed in terms of a sinusoidal tide as in (11). Integrating over half a tidal cycle (either flood or ebb) yields L p ≈(T/2)(c p +2U Ta /π), where 2U Ta /π is the mean tidal current during one phase (flood or ebb) of the tide. Writing c p in terms of the salinity anomaly as before, the ratio of the plume length to the tidal excursion (L Ta , as defined earlier) is approximately: Plumes experiencing substantial tidal reversals are expected to be approximately within two tidal excursions from the mouth (L p < 2L Ta ) because c p should be less than 2U Ta /π in these strong tidal regimes. For the Connecticut River, 2L Ta ranges from 22 km during neap tides up to 36 km during spring tides. The actual plume extent is smaller mostly because c p typically is much less than 2U Ta /π, particularly during spring tides. Results indicate the downstream front has very limited connectivity with interior waters. These waters are approaching the front with an overtake velocity equaling the difference between the interior water velocity and the frontal advancement velocity: u = u i -u f . This overtake velocity does not depend on the ambient tidal velocity. The time available for interior waters FIGURE 13 | Surface velocity (arrows) and salinity (shaded) maps at max ambient flood for (A) low discharge during neap tides (as in Figure 4B), (B) low discharge during spring tides, (C) high discharge during neap tides and (D) high discharge during spring tides. Velocity arrows are shown for every 40th point (800 m spacing) and every 20th point (400 m spacing) in x and y, respectively.
to travel from the mouth to the advancing front position (x f ) is at most the time between when the interior water first exits the mouth (t i ) and when the ambient tidal phase (in this case flood) ends at T/2. The downstream extent of the interior water (x i ) is less than x f while approaching the front and equals x f after reaching the front, if that occurs. Assuming constant u i and u f for simplicity, the condition for reaching the front by the end of the tidal phase is u i (T/2-t i ) ≥ u f T/2. Expressed as a velocity ratio, the condition for interior water interacting with the downstream front is: Note that u f and u i , respectively can be scaled as c p +u a and u+u f , so there is an implicit dependence on ambient tidal currents that tends to reduce u i /u f with stronger tides. Lateflood interior water exits the mouth with t i ≥ T/4 and (18) indicates that these waters can only reach the downstream front if u i /u f ≥ 2. This is a difficult threshold to meet, so it is unlikely that late-stage waters interact with the downstream front of tidally reversing plumes. For the Connecticut River plume, three of the runs have u i /u f ≤ 1.2 (Table 1), so t i ≤ 1 h to reach the downstream front. For the low-discharge spring-tide case, u i /u f = 1.5 and t i ≤ 2 h to overtake the front. Thus, interaction with the downstream front is limited to source waters exiting the mouth by early ambient flood. The tracers are identified by the first time they are released within the river, which can be several hours earlier than when they exit the mouth. For the standard low-discharge neap-tide case, the tracer released at −1.3 h is the youngest tracer with t i ≤ 1 h ( Figure 3C) and, consistent with the u i /u f value and (18), it is the youngest water to interact with the downstream front ( Figure 8C).
The present results are most directly comparable to the Cole et al. (2020) plume connectivity study. Cole et al. (2020) studies the Merrimack River plume, which is fed by a river-ebb pulse of fresh source waters and evolves in a much lower-energy ambient tidal environment than the Connecticut River plume experiences. Simulated tracers injected at the Merrimack River mouth indicate that only source waters released within the first 2 h (t i ≤ 2 h) overtake the front within a half tidal cycle after plume formation. This limited frontal connectivity is similar to results for the downstream front of the Connecticut River plume.
A key difference between the plumes is weak ambient tides allow the Merrimack plume to spread almost radially outward from the mouth, while strong ambient tides arrest the upstream front (Fr Ta > 1) in the Connecticut plume and elongate it in the tidally FIGURE 14 | Surface maps of average freshwater composition (indicated by average release times of tracers present in each location at max ambient flood) for (A) low discharge during neap tides (as in Figure 7A), (B) low discharge during spring tides, (C) high discharge during neap tides and (D) high discharge during spring tides. Tracer release times are expressed relative to the onset of ambient flood tides (as in Figure 7). Average freshwater age can be calculated by subtracting the average release times from 3 h (the mapped time). downstream direction. The Connecticut's arrested upstream front has strong and immediate connectivity with all source waters exiting the mouth. There is no part of the Merrimack plume front that has such strong connectivity. The strong connectivity is not guaranteed for all plumes with an arrested upstream front, as it is possible for later source waters to slip past in the plume interior without interacting with the bounding front. For the Connecticut plume, energetic ambient tides lead to large differences in connectivity and source-water composition along the bounding plume front. It is likely other plumes with a tidally arrested upstream front (where Fr Ta > 1) exhibit a similar transition from youngest to oldest waters progressing from the upstream front, along the offshore and onshore fronts, and to the downstream front. Consequently, such plumes will have more source-front connectivity than plumes in weaker ambient tidal environments, but the downstream front still will have very limited connectivity. Relative to plumes in weak ambient tidal regimes, the Connecticut River plume and other tidally reversing plumes should have key differences with respect to the fate and transport of pollutants, sediments, and biogeochemical substances. The relatively short length of tidally reversing plumes should tend to keep the highest concentrations close to the mouth. Spatial gradients of substances along the river can result in large differences in transport pathways. Material passing through the mouth in late ambient ebb or early flood will tend to be transported the farthest in plumes formed during ambient flood tides; this material can be processed at the downstream front. Substances in late-flood source waters would tend to remain closer to the mouth and only interact with the arrested upstream front. The strong connectivity of the arrested upstream front and its constrained location make it a good location for studying the frontal processing of pollutants, sediments, and biogeochemical substances. Future research incorporating plume biogeochemistry and source-front connectivity characteristics would help advance the understanding of plume effects on coastal ecosystem dynamics. It is also important to connect processes acting within the tidally reversing active plume to tracking substances over longer time scales (e.g., days to months), after plume waters have been extensively mixed with ambient waters.
Entrainment and mixing occur throughout the Connecticut River plume; as indicated by the results that concentrations and variances for all hourly tracers diminish quickly (Figure 11). The plume-averaged tracer variance budgets (10) for the combined tracer are qualitatively similar for all runs: tracer variance is increased by mouth inputs, decreased by entrainment, and destroyed by internal mixing. Mouth inputs increase tracer variance because source waters have higher concentrations than the plume (leading to C 2 m /C 2 >1), a common situation for plumes. Entrainment decreases tracer variance because C 2 /C 2 e >1; this condition is approximately C 2 /C 2 >1 because C 2 e is approximately C 2 for small C 0 . The C 2 /C 2 >1 condition should be met by other strongly mixed low-concentration plumes. The tracer variance budget terms vary by orders of magnitude among runs (Table 1). It is useful to simplify and scale the budget in a way that collapses results for the investigated range of tidal and discharge conditions. For C 2 m /C 2 1 and C 2 /C 2 1, the mouth input term and entrainment term in (10) are approximately C 2 m Q m /V p and −C 2 Q e /V p , respectively. These simplifications are reasonable approximations for the Connecticut River plume results. Analysis of tracer variance budgets of other plumes should include evaluating the C 2 m /C 2 and C 2 /C 2 ratios. The latter ratio is larger for strongly mixed low-concentration (low freshwater fraction) plumes, while the former ratio will be greater than one for any case with much higher concentration (fresher) source waters than the plume. The estimated mouth input term (C 2 m Q m /V p ) at t = 0 h (onset of ambient floods) can be used to scale the timevarying terms in (10). This mouth-based scaling reduces term magnitude differences among runs from orders of magnitude to approximately a factor of two. The successful collapsing of budget results indicates a dose-response relationship, whereby stronger tracer variance inputs by source waters have correspondingly larger variance decreases by entrainment and destruction via internal mixing.
Plume-averaged tracer budgets are valuable for analyzing plume evolution as a whole, but combine together horizontal and vertical processes as well as processes occurring in different regions of the plume. Future analysis of the Connecticut River plume should explore the decomposition of tracer variance into its horizontal and vertical parts, as introduced in Li et al. (2018) and applied in other studies (e.g., Wang and Geyer, 2018;Warner et al., 2020). The vertical tracer variance is particularly useful since it is directly related to stratification and it may be more readily applied to understanding biogeochemical processes in plumes. Analyzing spatial variations of tracer variance and budget terms as in Li et al. (2018) and other works (e.g., Warner et al., 2020) would help highlight important areas for variance advection and dissipation. The multiple tracer approach applied in the present study already allows for distinguishing between near-field and far-field plume regions. The mixing term in the tracer budgets for the present study includes the total mixing acting in the model, which includes physical mixing and numerical mixing associated with the advection scheme. Several studies have shown that numerical mixing can be comparable to physical mixing (e.g., Burchard and Rennau, 2008;Li et al., 2018;Burchard et al., 2021;Wang et al., 2021) and the ratio of numerical to physical diffusivities can have pronounced spatial variability within estuaries and plumes (Burchard et al., 2021).
Mixing in idealized tidal plumes representative of the Connecticut River plume and other situations has been studied in more detail in Spicer et al. (2021), following the turbulent buoyancy flux approach of Pritchard and Huntley (2006). Spicer et al. (2021) finds that bottom-generated tidal mixing dominates over mixing associated with buoyant plume shear if the estuarine Richardson number (Ri E = c m 2 V D /U Ta 3 ; modified from Fischer (1972) to use ambient tides) divided by the mouth Rossby number (Ro m = V D /(fW m ), where W m is mouth width) is less than one (Ri E /Ro m < 1). Provided that c p ≤ c m , this condition is satisfied when the product of the ambient tidal Froude number squared and the corresponding Rossby number exceeds one (Fr Ta 2 Ro Ta > 1, where Ro Ta = U Ta /(fW m )), which is true for all model runs in this study. Thus, the Connecticut River plume belongs in the class of plumes strongly influenced by tidal mixing (Spicer et al., 2021).
Other plumes in energetic ambient tidal environments, whether generated during flood or ebb, are expected to share similarities with the Connecticut River plume. These tidally reversing plumes should be within two tidal excursions from the mouth, be cut off from source waters while flow reverses at the mouth, have limited connectivity at the downstream front, strong connectivity along the arrested upstream front, and rapid entrainment and mixing that reduces the plume buoyancy signature over intratidal time scales. Follow-up work should compare the present results to connectivity within the Connecticut's ambient ebb plume. Future research should extend the connectivity and tracer analysis approach to the transport and processing of pollutants, sediments, and biogeochemical material by the Connecticut River plume and other plumes in energetic ambient tidal environments.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are publicly available. This data can be found here: https://opencommons. uconn.edu/marine_sci/12.

AUTHOR CONTRIBUTIONS
MW led this research, completed analysis, constructed figures, and wrote this manuscript with input from co-authors. YJ completed the hydrodynamic model runs, contributed to analysis, and edited and suggested improvements to the manuscript. KC, DM, and KH contributed to the ideas underpinning this research, helped set the scope of this manuscript, and improved the manuscript through their suggestions. All authors contributed to the article and approved the submitted version.