Observation of the Dynamics and Horizontal Dispersion in a Shallow Intermittently Closed and Open Lake and Lagoon (ICOLL)

: Among the 135 or more estuaries on the southeast Australian coastline, 45% are intermittently open. A number of others, which would normally have an intermittently open entrance, are kept permanently open to improve the ﬂushing and water quality that are functions of the horizontal dispersion of the estuarine system. Field studies that investigate the dynamics and horizontal dispersion processes of open and closed states of Intermittenly Closed and Open Lakes and Lagoons (ICOLL) are conducted. Clusters of Lagrangian drifters were used to characterise the surface ﬂow. Horizontal dispersion coefﬁcients ( K ) under study varied signiﬁcantly for different tidal phases; the mean ebb tide K (~4.7 ± 3.8 m 2 s − 1 ) is an order of magnitude higher than that of the ﬂood tide. During the closed state, K values were up to two orders of magnitude lower than those observed in the open state. These results highlight the contribution of tidal dispersion to transport and mixing processes within an ICOLL. The results showed that tidal pumping effects and tidal-induced horizontal velocity gradients are responsible for the horizontal surface transport and dispersion occurring during the open inlet state. The horizontal surface dispersion processes identiﬁed at different inlets and environmental conditions within this site are valuable for managing transport of particles, freshwater mixing, larvae transport, waste, and pest control.


Introduction
Intermittently Closed and Open Lakes and Lagoons (ICOLLs) are a subcategory of coastal lakes and lagoons, which exhibit limited exchange with the ocean and river flows [1]. Of the more than 135 estuaries along the southeast Australian coast, about 45% of these are intermittently open and a number of others, which would normally have intermittently open entrances, are kept permanently open. The permanent open state is facilitated by mechanical means to maximize the flushing and water quality that are functions of the horizontal dispersion of the system [2,3]. Many prior studies on ICOLLs have focused on the change in ecological processes caused by their dynamic opening and closing [4,5]. Hydrodynamic processes including vertical mixing, and horizontal dispersion, flushing, and recirculation patterns control the ecological processes [6,7]. Previous research demonstrates (iii) identifying the dominant dispersion mechanism for the inlet conditions. The purpose of this article is to present estimates of dispersion coefficients and dominant dispersion mechanisms from field studies in an ICOLL that provide an insight into transport of material in the systems during open and closed conditions, for time scale less than a tidal cycle and length scale in the order O [1 m].

Study Area
Currimundi Lake (Longitude 153.13 • E, Latitude 26.763 • S) is a coastal lagoon located on Queensland's Sunshine Coast in eastern Australia (Figure 1). The local climate is subtropical with an annual rainfall greater than 1200 mm with most rainfall (>60%), in the summer and part of autumn, between December and March. The location is characterized with average wind speeds between 5 m s −1 from the NE-E direction in the second half of a year and 6.8 m s −1 from the S-SE directions the first half of a year. The average daily global solar exposure is 5 kWh/m 2 (208.3 W/m 2 ) with the lowest of 2.9 kWh/m 2 (121 W/m 2 ) in June and highest of 7.6 kWh/m 2 (317 W/m 2 ) in January [24]. The catchment area for Currimundi Lake is highly urbanised having experienced significant growth commencing around the 1950s. The lower main channel of Currimudi Lake is connected to the ocean and upstream links to a range of stormwater infrastructure, canal networks, and a constructed brackish water body (Lake Kawana). To ensure that water circulation and turnover rates within Lake Kawana are maintained, an exchange system draws freshwater from the Mooloolah River System through Lake Kawana, discharging into the Currimundi Lake system via a weir (Figure 1a) at an average rate of 80 mL/day when the mouth is open. During the closed state, the discharge rate varies between 0-80 mL/day depending on the need to supplement the lake water level inundation of biting midge breeding sites located within the intertidal areas.

Study Area
Currimundi Lake (Longitude 153.13° E, Latitude 26.763° S) is a coastal lagoon located on Queensland's Sunshine Coast in eastern Australia (Figure 1). The local climate is subtropical with an annual rainfall greater than 1200 mm with most rainfall (>60%), in the summer and part of autumn, between December and March. The location is characterized with average wind speeds between 5 m s −1 from the NE-E direction in the second half of a year and 6.8 m s −1 from the S-SE directions the first half of a year. The average daily global solar exposure is 5 kWh/m 2 (208.3 W/m 2 ) with the lowest of 2.9 kWh/m 2 (121 W/m 2 ) in June and highest of 7.6 kWh/m 2 (317 W/m 2 ) in January [24]. The catchment area for Currimundi Lake is highly urbanised having experienced significant growth commencing around the 1950s. The lower main channel of Currimudi Lake is connected to the ocean and upstream links to a range of stormwater infrastructure, canal networks, and a constructed brackish water body (Lake Kawana). To ensure that water circulation and turnover rates within Lake Kawana are maintained, an exchange system draws freshwater from the Mooloolah River System through Lake Kawana, discharging into the Currimundi Lake system via a weir (Figure 1a) at an average rate of 80 mL/day when the mouth is open. During the closed state, the discharge rate varies between 0-80 mL/day depending on the need to supplement the lake water level inundation of biting midge breeding sites located within the intertidal areas.  The Currimundi Lake main channel is characterised by a mixed tide with a predominantly semi-diurnal tidal pattern and is a micro-tidal estuary with a spring tidal range limited to 0.8 m. The depth of the channel ranges between 3-5 m mid-estuary. The system is considered an ICOLL with the recorded natural opening and closing events prior to the current field experiments, in the early 2000s and was naturally closed recently in July 2017 [25]. However, from the mid-2000s it has been artificially disconnected from the ocean once per year, by blocking the inlet with sand, usually for about six weeks in September/October. Three different drifter experiments were conducted along the main channel during tidal (open state-E1 in 27-28 April 2015) and non-tidal (closed state-E2 in 9 September 2015) conditions, and during the opening of the inlet (E3 in October 2015). Here we focus on the effect of the inlet condition on the dynamics and the horizontal surface dispersion of channel flow by contrasting observations during the open (E1) and closed state (E2) experiments. The forcing conditions for the two experiments are discussed in details in the next Section 2.2.

Environmental and Inlet Conditions during the Field Experiments
The major forcing conditions in Currimundi Lake during an open condition include tide, river discharge, wind and rainfall. In the absence of tidal forcing in a closed state, the solar radiation and heat convection potentially become important. A range of tide, wind, solar exposure and flow conditions was encountered during the field studies (Table 1). During the open state (E1), the average tidal range was 0.6 m. The average fresh water input into the channel from Lake Kawana was 80 mL/day for the open state. During the closed state, the water level indicated a gradual increase of 0.1 m between 31 August 2015 and 2 September 2015, likely from the combination of storm water run-off and freshwater from Lake Kawana, and the water level then remained approximately constant at 0.78 m AHD through the closed state experiment (Figure 2c). The total rainfall within the catchment seven days prior to the measurements was less than 10 mm for both open and closed states. Therefore, runoff from storm water drains was assumed to be insignificant. The wind velocity peaked at 4 and 7 m s −1 from the S-SE directions for the open and closed experiments, respectively. Table 1. Summary of the experimental conditions and instruments deployed at Currimundi Lake; flow velocity and wind speed were measured using ADV, weather station and sonic anemometer deployed at the Pontoon (Figure 1b). Water level, salinity, and temperature of the channel are continuously monitored by the Sunshine Coast Council (SCC) at the bridge marked X in Figure 1b. The water level in the channel was between 0.08-0.78 m Australian Height Datum (m AHD) during the open state (E1), while the water level was constant at 0.78 m AHD during the closed state (E2) experiment ( Figure 2). The water levels are shown relative to the surveyed bed profile near the pontoon in Figure 3. A Conductivity, Temperature and Depth meter, CTD3100, was used to sample water quality at a temporal resolution of 5 min placed at −0.32 m AHD elevation. Figure 2 shows the time series of water level and salinity during the open state (E1) and the closed state (E2) experiments. During the open state (E1), the water level and salinity exhibited tidal trends. The diurnal trend in the salinity observed during the closed state (E2) is likely associated with daily thermal and wind mixing in the water column.

Instrumentation
To examine the horizontal surface dispersion within the channel, Lagrangian drifters were deployed. The drifters used in this study are similar in shape and design, but different in dimension and electronic components to high resolution (HR) drifters described by [13]. These low cost, low resolution (LR) drifters are made of PVC cylindrical pipes with 4 cm diameters and lengths of 25 and 50 cm for the short and long types, respectively ( Figure 1d). The drifters contain off-the-shelf Holux GPS data loggers with absolute position accuracy of 2-3 m, sampled at 1 Hz. The effect of this error on analysed parameters is discussed in Section 2.3. The drifters were positively buoyant for continuous satellite position fixation with <3 cm height unsubmerged to limit wind effects [13,26]. The current configuration results in direct wind slip estimated as less than 1% of the ambient wind for the short drifter version [13,[26][27][28]. The drifters showed good agreement (R 2 ~ 0.9) with surface velocity measurements obtained from an Acoustic Doppler Current Profiler (ADCP) within a 60 m radius [26,29]. However, unfiltered drifter data are not suitable for measuring small-scale processes in low flow applications because of the inherent position error that manifests as a large speed variance; σ 2 ~ 0.0005 m 2 s −2 in stationary tests at frequencies F > 0.01 Hz [23].
Although a Lagrangian drifter has been shown to be applicable to studying shallow water systems, there are clear limitations to its application to ICOLL due to stratification and ensuing baroclinic motions, particularly in closed states. The drifter used in the current study is designed to follow the upper 0.5 m layer of the water column; therefore, the observed dispersion coefficients and dynamics are representative of horizontal surface dispersion. In addition, the drifter only captures motion on a scale greater than its radius; therefore, the size limits the range of eddies it captures. Further limitations of drifter application to examining dispersion in shallow waters are discussed in [13].
During the open state (E1), a fleet containing 10 short and 10 long drifters was deployed in the ebb and flood tides on 27 April 2015. Eighteen drifters were retrieved at the end of the experiment; two drifters were lost. The ebb deployment occurred at Site A (Figure 1b), adjacent to the main channel, while the flood deployment was made at Site C, upstream of the bridge marked X. The drifters were deployed in a line abreast configuration with approximately 1 m separation [23]. Drifters trapped along banks were collected and redeployed at the adjacent centerline. Data quality control is described in Section 2.3. During the 4-h ebb tide experiment, the number of times an individual drifter reached the banks varied from two to nine with an average of six times per drifter. Conversely, only five of the 18 drifters were redeployed during the 3-h flood tide experiment. This suggests that the ebb surface current induced many recirculation zones close to the banks, trapping the drifters. During the closed state (E2), all deployments occurred at Site B. Drifters were deployed in four clusters (C1, C2, C3 & C4), each containing 3-5 drifters at the centre of the channel. Two

Instrumentation
To examine the horizontal surface dispersion within the channel, Lagrangian drifters were deployed. The drifters used in this study are similar in shape and design, but different in dimension and electronic components to high resolution (HR) drifters described by [13]. These low cost, low resolution (LR) drifters are made of PVC cylindrical pipes with 4 cm diameters and lengths of 25 and 50 cm for the short and long types, respectively ( Figure 1d). The drifters contain off-the-shelf Holux GPS data loggers with absolute position accuracy of 2-3 m, sampled at 1 Hz. The effect of this error on analysed parameters is discussed in Section 2.3. The drifters were positively buoyant for continuous satellite position fixation with <3 cm height unsubmerged to limit wind effects [13,26]. The current configuration results in direct wind slip estimated as less than 1% of the ambient wind for the short drifter version [13,[26][27][28]. The drifters showed good agreement (R 2~0 .9) with surface velocity measurements obtained from an Acoustic Doppler Current Profiler (ADCP) within a 60 m radius [26,29]. However, unfiltered drifter data are not suitable for measuring small-scale processes in low flow applications because of the inherent position error that manifests as a large speed variance; σ 2~0 .0005 m 2 s −2 in stationary tests at frequencies F > 0.01 Hz [23].
Although a Lagrangian drifter has been shown to be applicable to studying shallow water systems, there are clear limitations to its application to ICOLL due to stratification and ensuing baroclinic motions, particularly in closed states. The drifter used in the current study is designed to follow the upper 0.5 m layer of the water column; therefore, the observed dispersion coefficients and dynamics are representative of horizontal surface dispersion. In addition, the drifter only captures motion on a scale greater than its radius; therefore, the size limits the range of eddies it captures. Further limitations of drifter application to examining dispersion in shallow waters are discussed in [13].
During the open state (E1), a fleet containing 10 short and 10 long drifters was deployed in the ebb and flood tides on 27 April 2015. Eighteen drifters were retrieved at the end of the experiment; two drifters were lost. The ebb deployment occurred at Site A (Figure 1b), adjacent to the main channel, while the flood deployment was made at Site C, upstream of the bridge marked X. The drifters were deployed in a line abreast configuration with approximately 1 m separation [23]. Drifters trapped along banks were collected and redeployed at the adjacent centerline. Data quality control is described in Section 2.3. During the 4-h ebb tide experiment, the number of times an individual drifter reached the banks varied from two to nine with an average of six times per drifter. Conversely, only five of the 18 drifters were redeployed during the 3-h flood tide experiment. This suggests that the ebb surface current induced many recirculation zones close to the banks, trapping the drifters. During the closed state (E2), all deployments occurred at Site B. Drifters were deployed in four clusters (C1, C2, C3 & C4), each containing 3-5 drifters at the centre of the channel. Two separate deployments were made for each cluster during the 6-h field experiment. The drifters all moved upstream toward the pontoon (Figure 1), due to wind forcing.
A Sontek ® Acoustic Doppler Velocimeter (ADV) sampling at 50 Hz was mounted rigidly to the pontoon 0.5 m below the water surface to capture surface flow velocity. Figure 3 shows the profile of the transect next to the pontoon and the corresponding location of the ADV. Coinciding with the drifter experiments during the closed state (E2), vertical profiles of DO, salinity, and temperature were collected by the Sunshine Coast Council at the pontoon and the bridge (X) (Figure 1) using Hydrolab DS5X. In addition, a two-dimensional wind anemometer was mounted next to the pontoon about 0.5 m above the water surface to capture the wind speeds and directions at the free surface.

Data Analysis
Raw drifter data were converted to local East-North-Up (ENU) coordinates from a World Geodetic System (WGS84). The drifter position data were quality controlled using velocity and acceleration de-spiking [13]. Sections of the trajectories where the drifters did not float in the channel, where they were found next to the bank or affected by occasional proximity to boats were removed using MATLAB scripts developed to work with the experimental event log. A low-pass filter with cut off frequency, F c = 0.01 Hz, was then applied to the position time series. Drifter velocities relative to East and North were obtained by combining the quality controlled position and speed time series, which are directly logged by holux data-loggers: where V E and V N are the Easting and Northing velocities, respectively, θ is the direction based on the position time series (E,N), and Sp is drifter speed at time t.
The objective of the cluster analysis herein was to estimate the time variation mean values of the dispersion coefficient (K c ) and differential kinematic properties (DKP) for individual drifter deployments. The dispersion coefficient was obtained using the method described by [30]. The key aspect of the method is the calculation of the rate of change of a patch size with the assumption that the drifters are 'locked' to the 'patch' of fluid within which they were deployed.
Cluster analysis only commenced once data from all drifters in the cluster are simultaneously logged. However, the initial 300 s of measurement, when drifter motion was affected by boat movement and initial deployment, was excluded. For the open state (E1), where drifters were frequently beached and re-introduced into the channel, times with drifter numbers n < 5 were noisy and were not included in the results; therefore, at least five drifters were used in the calculation. For the closed state (E2), calculations were terminated immediately when one drifter in a single cluster was beached. The dispersion coefficient (K c ) estimates for long and short drifter clusters were similar during most of the experiment. Therefore, both drifter types were combined into a single cluster to improve the stability of the estimate. Using the local ENU coordinates, the centroid (represented with overbar) of a cluster is defined as: where i is the drifter counter and n is the total number of active drifters in a cluster at time, t. The variance of an individual drifter from the centroid of the cluster is then defined as: Water 2018, 10, 776 8 of 21 The cluster relative dispersion coefficient is calculated from the averaged variance such that: DKP values were calculated from velocity gradients obtained via a Taylor's series expansion about the centre of mass of a cluster, following the method described in [31]. This method assumes that the fluid domain is small and finite in size, such that the velocity gradient is uniform across a cluster. Therefore, the calculations were not conducted when a member of a cluster was within close proximity of a bank where strong velocity gradients are expected. The method further assumes that cluster velocity is correctly represented in the linear (first order) term of a Taylor's series, while the second and higher order terms are considered as turbulence. Individual drifter velocity can be described as: where ∂V E ∂E , ∂V E ∂N , ∂V N ∂E and ∂V N ∂N are linear centroid velocity gradient terms, and v N and v E are non-linear turbulence velocity terms (second and higher order terms). The velocity gradient and non-linear terms are estimated using a least squares approach which converts the equation into a matrix format [31]. The non-linear turbulence terms are not presented here because their magnitudes are not significantly higher than the standard deviation of the speed of the drifter placed in a fixed location (i.e., during a stationary test). DKPs were then described in terms of the resulting velocity gradients such that: The horizontal divergence (δ) is a measure of change in area of the cluster without a change in orientation or shape, while relative vorticity (ζ) is a measure of the change in the cluster orientation without a change in area or shape. Shearing (b) is the measure of change in shape produced by boundaries parallel to the cluster boundary without change in the area or orientation, while stretching (a) is a change in shape without a change in area and orientation, but is caused by boundaries perpendicular to the cluster boundary [32]. Truesdell's kinematic vorticity number, T K , is defined to identify the periods (regions) of strain and vorticity field [33]: T K > 1 corresponds to vorticity dominance or the presence of stronger eddy-like structures, whilst T K < 1 corresponds to strain field dominance or periods (regions) of convergence or divergence where dispersion is stronger [17].
The growth of a cloud of particles (such as tracers, larvae and buoyant pollutants) in the surface layer is governed by several mechanisms which include small scale turbulent eddies at length scale smaller than the size of the cloud and large scale dispersion processes at length scale larger the size of the cloud. The dispersion coefficient measures the change in the size of a cloud with time. The dispersion coefficient estimated using this cluster method is an apparent dispersion coefficient because it includes the effect of both small scale eddies diffusion and horizontal shear dispersion [34]. For a shallow tidal channel, within a period less than a tidal cycle, the dispersion coefficient varies as a function of length scale of the patch with a power weaker than the Richardson's 4/3 relation [35]. Therefore, the average dispersion coefficient estimate for a period is biased towards the largest Water 2018, 10, 776 9 of 21 dominant eddies encountered by the patches while the standard deviation reflects the range of eddy length scales encountered during the period.

Measurement Uncertainty from Fixed Stationary Drifter
The estimate of the centroid, velocity, dispersion coefficient and DKP associated with inherent errors due to GPS position fixing and hardware noise is examined. Assuming the GPS position fixing is independent of drifter motion and location, measurement taken at a fixed location is representative of the inherent errors [13,36]. In order to estimate the magnitude of this errors on the estimates of the key analysis parameters, three independent observations were obtained at fixed locations using the GPS data logger in the low resolution drifters. Four dataloggers separated by~1 m apart were sampled continuously for 4-8 h. The dataset was quality controlled to synchronise the drifter position time series, lowpass filtered with cut-off frequency, F c = 0.01 Hz and analysed using relevant equations in Section 2.3. The results are summarised in Appendix A. For a fixed measurement, theoretical K c is zero if the GPS position fixing and hardware noise are zero. These measurements showed that while the inherent position accuracy of the drifters was between 2-3 m, the estimates of the K c based on the measured noise were at least an order of magnitude less than the values presented in Section 3. Correspondingly, magnitudes of DKP based on the noise were at least an order of magnitude less than those reported in Section 3. The few exceptions to this condition of minimum of an order of magnitude difference between inherent error and estimate of DKP, from Currimundi Lake, are flagged in Section 3.3.

Basic Flow and Physiochemical Observations
During the open state experiment (E1) when the channel was connected to the ocean, measured salinity varied between 25 and 33 ppt. Horizontal transport was governed by river and tidal inflow. Comparison of water levels within the channel with tidal data from tidal gauge in the ocean indicated that the high and low waters in the channel lag the ocean by 2-3 h (not shown). The time series of mean flow velocity magnitude and the water level at Currimundi Lake is presented in Figure 4. Flow direction is indicated by the direction of the mean streamwise velocity component, which is positive in the upstream direction. The mean flow velocities were obtained using a moving average based on a 200 s average every 10 s along the entire dataset which is appropriate for estuaries of this type [37]. The channel exhibits a rapid and short flood tide compared to the ebb tide. The flood period occurred within 4-h with a peak velocity of 0.25 m s −1 , while ebb tide extended to a period up to 9-h with peak velocity limited to 0.18 m s −1 . The short-lived flood tide and the tidal asymmetry within the channel were likely influenced by the net outflow from Lake Kawana. The end (−3600-3700 m North; Figure 5) of the main channel exhibited secondary flows during the ebb, and higher flow velocity magnitudes >0.4 m s −1 during the peak flood tide due to constriction of the channel. This accelerated flow in these narrow reaches, combined with the short-lived flood flow (Figure 5a), could accelerate streambank erosion in this portion of the channel, especially in areas not protected by woody riparian vegetation. As such, changes in channel morphology could occur in these upstream reaches.
Except for the locations in the constricted end of the main channel, variations in instantaneous velocity within the channel were limited to 0-0. run-off and from Lake Kawana, via the weir, after the inlet was closed (Figure 2c). The wind strength (speed < 5 m s −1 ) observed near the surface was weak and unable to disrupt the stratified layer resulting in salinity gradient in the channel system. This inhibition of vertical mixing probably explains the depletion of oxygen near the channel bed through the consumption of oxygen by organisms in sediment-water layer ( Figure 6).   run-off and from Lake Kawana, via the weir, after the inlet was closed (Figure 2c). The wind strength (speed < 5 m s −1 ) observed near the surface was weak and unable to disrupt the stratified layer resulting in salinity gradient in the channel system. This inhibition of vertical mixing probably explains the depletion of oxygen near the channel bed through the consumption of oxygen by organisms in sediment-water layer ( Figure 6).  Daytime vertical profiles of the physiochemical properties of the channel during the closed state (E2) indicated that the water column was strongly stratified in density from the salinity gradient ( Figure 6b). With an insignificant vertical temperature gradient as observed during the closed state (E2), the source of stratification was likely connected with the freshwater inflow from stormwater run-off and from Lake Kawana, via the weir, after the inlet was closed (Figure 2c). The wind strength (speed < 5 m s −1 ) observed near the surface was weak and unable to disrupt the stratified layer resulting in salinity gradient in the channel system. This inhibition of vertical mixing probably explains the depletion of oxygen near the channel bed through the consumption of oxygen by organisms in sediment-water layer ( Figure 6). Water 2018, 10, x FOR PEER REVIEW 11 of 21

Horizontal Surface Dispersion Coefficient
In the open state (E1) during the flood time, the time variation of the cluster dispersion coefficient, Kc and centroid horizontal resultant velocity, VH, varied throughout, which is a reflection of the tidal influence ( Figure 7). However, temporal variation of Kc reflected different stages characterised by divergence (Kc > 0; Stages 1 and 2) and convergence (Kc < 0; Stages 3 and 4) of the cluster. Regions in the channel corresponding to these stages are identified in Figure 7b; drifter trajectories are colored by time. The magnitude of mean Kc for the stages was not significantly different from the overall mean. Of interest is the pronounced oscillation apparent in Kc, VH, and cluster area manifested as two divergence/convergence events (Figures 8 and 9a). This modulation of divergence and convergence has a period ~3000 s (e.g., Stage 2). Similar modulation of period ~3000 s can be observed in the time series of velocity of the cluster centroid and ADV (Figure 4). This corresponds to an internal resonance/wave of length 3.3 km in shallow water with an average depth, d ~ 0.5 m. This is probably due to resonance between the mouth and the fixed weir located 3.6 km from the channel inlet. The resonance effect is a reversible process manifested as cluster expansion/contraction [17,38]. However, this effect is expected to have some impact on the dynamics of the surface layer of the channel by destroying stratification and enhancing mixing.

Horizontal Surface Dispersion Coefficient
In the open state (E1) during the flood time, the time variation of the cluster dispersion coefficient, K c and centroid horizontal resultant velocity, V H , varied throughout, which is a reflection of the tidal influence ( Figure 7). However, temporal variation of K c reflected different stages characterised by divergence (K c > 0; Stages 1 and 2) and convergence (K c < 0; Stages 3 and 4) of the cluster. Regions in the channel corresponding to these stages are identified in Figure 7b; drifter trajectories are colored by time. The magnitude of mean K c for the stages was not significantly different from the overall mean. Of interest is the pronounced oscillation apparent in K c , V H , and cluster area manifested as two divergence/convergence events (Figures 8 and 9a). This modulation of divergence and convergence has a period~3000 s (e.g., Stage 2). Similar modulation of period~3000 s can be observed in the time series of velocity of the cluster centroid and ADV (Figure 4). This corresponds to an internal resonance/wave of length 3.3 km in shallow water with an average depth, d~0.5 m. This is probably due to resonance between the mouth and the fixed weir located 3.6 km from the channel inlet. The resonance effect is a reversible process manifested as cluster expansion/contraction [17,38]. However, this effect is expected to have some impact on the dynamics of the surface layer of the channel by destroying stratification and enhancing mixing.

Horizontal Surface Dispersion Coefficient
In the open state (E1) during the flood time, the time variation of the cluster dispersion coefficient, Kc and centroid horizontal resultant velocity, VH, varied throughout, which is a reflection of the tidal influence ( Figure 7). However, temporal variation of Kc reflected different stages characterised by divergence (Kc > 0; Stages 1 and 2) and convergence (Kc < 0; Stages 3 and 4) of the cluster. Regions in the channel corresponding to these stages are identified in Figure 7b; drifter trajectories are colored by time. The magnitude of mean Kc for the stages was not significantly different from the overall mean. Of interest is the pronounced oscillation apparent in Kc, VH, and cluster area manifested as two divergence/convergence events (Figures 8 and 9a). This modulation of divergence and convergence has a period ~3000 s (e.g., Stage 2). Similar modulation of period ~3000 s can be observed in the time series of velocity of the cluster centroid and ADV (Figure 4). This corresponds to an internal resonance/wave of length 3.3 km in shallow water with an average depth, d ~ 0.5 m. This is probably due to resonance between the mouth and the fixed weir located 3.6 km from the channel inlet. The resonance effect is a reversible process manifested as cluster expansion/contraction [17,38]. However, this effect is expected to have some impact on the dynamics of the surface layer of the channel by destroying stratification and enhancing mixing.    (7)). Note: the number of active drifters 5-15; periods with drifter numbers, n < 5 are not included.   (7)). Note: the number of active drifters 5-15; periods with drifter numbers, n < 5 are not included.

Effect of Inlet on the Dispersion Coefficient
A single deployment resulted in a range of cluster dispersion coefficients, K c values and the data varied with the instantaneous effective length scale which is equivalent to the size of the patch. Due to the limited number of drifters in each individual cluster in some deployments, K c values obtained from Equation (4) were noisy. Therefore, only the mean values over cluster deployments are reported in Table 2. The definition of a cluster allows for negative values of dispersion coefficient to be calculated, indicating cluster contraction, i.e., rather than spreading. Dispersion coefficients resulting from spreading and contraction are separately averaged for individual clusters and are determined by taking the mean over a deployment duration. Contraction of clusters of buoyant particles as opposed to spreading has been observed in environmental flows such as lakes, estuarine embayments and nearshores [35,38]. This phenomenon is usually related to a combination of physical processes-e.g., secondary flow, 3D contraction effect on quasi 2D flow, proximity to fronts and banks and internal waves (as discussed in Section 3.2), and dynamics of particles in a turbulent flow influenced by their inertia and drag [39]. The effective dispersion coefficients as a function of the mean cluster velocity are compared for the two inlet conditions (Figure 8). The results indicate that the ebb K c = 4.73 ± 3.8 m 2 s −1 is an order of magnitude greater than that of the flood tide and up to two orders of magnitude greater than values obtained when the inlet was closed. The ebb cluster diffusivities were similar to those (K = 6.5-9.9 m 2 s −1 ) observed from dye tracer studies in natural rivers (depth = 0.58-1.56 m; width = 20-40 m) such as Green-Dumanish River and Powell River in the USA (Table 5.3; Fischer et al. [10]). A mean K c value of 0.42 ± 0.34 m 2 s −1 was obtained during the flood tide. This estimate is similar to K c = 0.5 m 2 s −1 , obtained in a similar independent drifter experiment at a tidal inlet [40]. The estimate is also similar to the dispersion coefficient of 0.57 m 2 s −1 obtained from absolute dispersion of high resolution drifters at Eprapah Creek, Australia, a smaller tidal channel [13].
The 1-2 orders of magnitude difference between the open and closed state dispersion coefficients highlight the role of tidal inflows in the horizontal transport and mixing within the channel. Mechanisms responsible for dispersion are those related to freshwater input and local wind effects. During the closed state (E2), wind intensity next to the channel surface was generally weak, less than 30% of the local wind velocity for the catchment (obtained at 4.5 m AHD). This is consistent with the boundary layer effect of mangroves and buildings [41] and explains the low dispersion coefficient.

Differential Kinematic Properties
Differential kinematic properties (DKP) are particularly useful for describing the structure of relative motion of clustered particles and for relating the velocity field to dynamic processes. This is important for studying passive particle transport such as larvae, litter, and logs in response to the velocity gradient in an underlying flow field. Divergence provides an indication of vertical motion whilst vorticity is a diagnostic of external forcing such as interaction of tide with bottom topography or wind-induced forcing in the absence of tide [32]. Shearing and stretching deformation in the context of semi-enclosed waters are a function of horizontal boundary influences on cluster motion.
We first present the variation of DKPs observed during the flood tide for the open inlet condition. DKP varied temporally and spatially across the channel (Figure 9). Mean DKP values were in the order of 10 −4 s −1 . Divergence and vorticity were predominantly positive, whilst shearing and stretching deformation were negative. Observed values were an order of magnitude larger than those observed in a non-tidal lake [17] and up to ten orders of magnitude lower than observed in the ocean [32]. The variation of T K with time suggests that the channel was systematically dominated by strain fields for most of the flood tide during the baseline tidal condition (Figure 9f). Mixing during these periods was likely linked with interaction of the tide with irregular bathymetry, boundaries, and channel roughness elements such as riparian vegetation. However, some periods of vorticity dominance were observed. These periods corresponded to divergence-convergence interphase linked with internal resonance within the channel as identified from the cluster dispersion.
Mean DKP values were in the order of 10 −4 s −1 ( Table 3). The standard deviation of DKP was, on average, an order of magnitude greater than the mean, indicating large variability and a limited number of particles in the clusters. This result indicates that on average the DKP of the channel was higher during the closed state (E2). Large values of DKP during the closed state (E2) are indicative of small spatial separation of the clusters. The drifter clusters tended not to separate great distances even 2 to 3-h after deployment, consistent with the observed lower dispersion coefficient. Consistent with wind driven transport, the vorticity and stretching effects were typically amplified during the closed state (E2) compared to the open state (E1). For all clusters and deployments, the time variation of T K was predominantly less than 1, suggesting that dispersion of the clusters was dominated by strain fields rather than eddy-like structures. The dominance of the strain field within the channel implies that water parcels are stretched and sheared horizontally. This effect was probably associated with horizontal velocity shearing resulting from the interactions of the tide with topographic structures within the channel (e.g., bars, banks, and vegetation) during the open state and wind-driven flow forced against the banks during the closed state.

Open State (E1)
During the open state (E1), horizontal dispersion in the natural channel is influenced by the combination of the dispersive effect of wind, tide, and river currents. In a long, narrow channel such as Currimundi Lake, dispersion resulting from wind-driven currents at moderate wind speeds (e.g., 0-5 m s −1 ) observed during the field study, is probably insignificant compared with the influences of the tide and river. The physical mechanisms for horizontal dispersion based on the competition between the tide and freshwater inflow are related to the additive effects of turbulence, shear dispersion, tidal trapping, and pumping. Turbulent mixing resulting from the flow boundaries such as beds and banks is parameterised as eddy diffusivity, K turb . Assuming well-mixed conditions, K turb could be expressed as [8]: where C D is the drag coefficient relating the square of depth-averaged velocity to the bottom stress, U t is the amplitude of the tidal velocity, and h is the averaged water depth. Using a typical range of U t = 0.05-0.3 m s −1 and h = 1-2 m observed during the open state (E1) and assuming C D = 0.003 after [42], the vertical mixing within Currimundi Lake can be estimated as  [27], who assumed a logarithmic turbulent velocity. The concept has been extended to oscillatory tidal flows [10]. The resulting horizontal dispersion, K shear , can be approximated by the tidal velocity and width, W as: Using the observed field tidal velocity and an average channel width, W~50 m, K shear is estimated between 0.25-3 m 2 s −1 during the open state (E1).
Tidal trapping occurs as a result of geometrical variation that leads to a difference in the phase of tidal flow between the main channel and branching channels [8]. Okubo [43] used an analytical approach to estimate the dispersion coefficient associated with this phenomenon based on the trapped volume in the branching channels, the tidal excursion length, L, and the tidal velocity amplitude: Assuming a trapped volume r~0.2 based on the ratio of the volume of the branching and main channels in Currimundi Lake, a tidal excursion length L t~2 000 m based on the tidal range, and depth and length of the main channel, K trap is estimated to vary between 5-30 m 2 s −1 . The tidal pumping effect, on the other hand, describes the net transport that is caused by the difference between the ebb and flood tidal currents [8]. This effect may significantly influence the horizontal dispersion in Currimundi Lake considering the difference in the flow behaviour between the flood and ebb tides observed with the ADV and drifters. However, a resulting dispersion coefficient from tidal pumping was not estimated due to difficulty in estimating the residual recirculation owing to the complexity of the channel network and dead ends.
By comparing the observed K c values for different experiments with those predicted using approximate solutions for individual mechanisms, the dominant underlying mechanisms can be explained. The K values observed during the open state (E1) were larger than the inherent eddy diffusivity within the channel, indicating dominance of larger scale processes. The flood tide mean K = 0.42 m 2 s −1 was similar in magnitude with that predicted by the shear dispersion mechanism. The ebb tide mean K = 4.73 m 2 s −1 , on the other hand, was similar in magnitude to that predicted by a trapping mechanism. This suggests that the horizontal dispersion within the channel during an open state is dominated by large-scale mechanisms such as horizontal shear dispersion and tidal trapping. Trapping within the channel could be linked to the influence of the narrow and shallow channel branches compared with the main channel. This geometrical effect reduces the resistance of the channel branches to inertia forces leading to quicker ebbing in the channel branches and potentially increases in erosive forces along banks. In addition, Currimundi Lake is connected to Lake Kawana via a weir head located at 0.65 m AHD. This structure restricts the over-passage of flow at tidal elevations less than 0.65 m AHD while freshwater is supplied into the system creating an induced trapping and tidal asymmetry.

Closed State (E2)
During the closed state (E2), the tidal effect was absent. Therefore, the driving forces for horizontal transport were expected to be from horizontal density gradients and wind forcing. The water level was fairly constant and there was no significant stormwater runoff into the channel during the study period. There was no significant difference between the density between the pontoon and the bridge. Therefore, the influence of horizontal density gradient currents within the channel was negligible. To examine the influence of wind on horizontal transport, drifter motion and wind velocity measured 0.5 m above the water surface during the closed state (E2) are presented using a rose diagram ( Figure 10). To obtain a better visual match with the drifter motion, the wind rose is presented in terms of the direction to where the wind was blowing (Figure 10), i.e., rotated by 180 o compared to conventional wind rose diagrams, and transformed to a local coordinate system. The directions illustrated in the wind rose approximately represent the ensemble average drifter direction and speed in a Eulerian perspective at a particular interval and do not necessarily represent the Lagrangian trajectory [44]. A comparison of the drifter rose with the wind rose should therefore be interpreted cautiously. The near surface wind speed was predominantly between 1-5 m s −1 blowing in a N-NW direction. The drifter speeds were mainly between 0-0.15 m s −1 and drifted towards N-NW. The result showed that the drifter speeds were between 0-5% of the wind speed and within 20 • of the wind direction. This suggests a significant contribution of wind to the dynamics of the surface flow during the closed state (E2) at low to moderate wind speeds observed within the channel. The drifter motion did display wind driven transport, such as net deflection to the left of the wind direction as expected in the southern hemisphere. However, the net angle of deflection over~0.5 m of the channel was limited to 20 • compared to the expected 45 • Ekman deflection at the surface. This was largely due to the small size and elongation of the main channel that substantially modified the Ekman influence. During the closed state (E2), the tidal effect was absent. Therefore, the driving forces for horizontal transport were expected to be from horizontal density gradients and wind forcing. The water level was fairly constant and there was no significant stormwater runoff into the channel during the study period. There was no significant difference between the density between the pontoon and the bridge. Therefore, the influence of horizontal density gradient currents within the channel was negligible. To examine the influence of wind on horizontal transport, drifter motion and wind velocity measured 0.5 m above the water surface during the closed state (E2) are presented using a rose diagram ( Figure 10). To obtain a better visual match with the drifter motion, the wind rose is presented in terms of the direction to where the wind was blowing (Figure 10), i.e., rotated by 180 o compared to conventional wind rose diagrams, and transformed to a local coordinate system. The directions illustrated in the wind rose approximately represent the ensemble average drifter direction and speed in a Eulerian perspective at a particular interval and do not necessarily represent the Lagrangian trajectory [44]. A comparison of the drifter rose with the wind rose should therefore be interpreted cautiously. The near surface wind speed was predominantly between 1-5 m s −1 blowing in a N-NW direction. The drifter speeds were mainly between 0-0.15 m s −1 and drifted towards N-NW. The result showed that the drifter speeds were between 0-5% of the wind speed and within 20° of the wind direction. This suggests a significant contribution of wind to the dynamics of the surface flow during the closed state (E2) at low to moderate wind speeds observed within the channel. The drifter motion did display wind driven transport, such as net deflection to the left of the wind direction as expected in the southern hemisphere. However, the net angle of deflection over ~0.5 m of the channel was limited to 20° compared to the expected 45° Ekman deflection at the surface. This was largely due to the small size and elongation of the main channel that substantially modified the Ekman influence.

Implications for Current Lake Management Strategies
Several lake management issues can be informed by these results: the horizontal dynamics of particles in the channel litter accumulation, and water quality in the upstream reaches, whilst pest control remains of concern throughout the ICOLL.
Biting midge management is a major concern in Australian coastal areas and across many other tropical and sub-tropical areas. Midge larvae are difficult to control by direct application of larvicides in their breeding sites because they are deposited in mud and wet soils high in organic matter [45]. While tidal cycles have been observed to aid breeding of intertidal biting midge species [19], recent research has shown that the larvae of many species are not able to survive inundation for a long period [20], thus creating an opportunity for larvae control at sites next to channel banks. Midge larvae management conducted by the local government in Currimundi Lake entails the annual

Implications for Current Lake Management Strategies
Several lake management issues can be informed by these results: the horizontal dynamics of particles in the channel litter accumulation, and water quality in the upstream reaches, whilst pest control remains of concern throughout the ICOLL.
Biting midge management is a major concern in Australian coastal areas and across many other tropical and sub-tropical areas. Midge larvae are difficult to control by direct application of larvicides in their breeding sites because they are deposited in mud and wet soils high in organic matter [45]. While tidal cycles have been observed to aid breeding of intertidal biting midge species [19], recent research has shown that the larvae of many species are not able to survive inundation for a long period [20], thus creating an opportunity for larvae control at sites next to channel banks. Midge larvae management conducted by the local government in Currimundi Lake entails the annual drowning of the larvae laid in sandy banks by blocking Currimundi Lake inlet for a period of approximately 14 days; this is the typical lifespan of biting midge in mid-summer [19]. The results presented herein support the practice of closing the channel inlet at a period coinciding with the peak midge breeding period as an effective control strategy. This is because dispersion was lower and significant flow fluctuations were absent during the closed inlet experiment; therefore, larvae in the banks of the channel likely remain inundated as the water level rises. In addition, the vertical profile water quality parameters showed a depletion of oxygen closer to the bed, thus reducing larval access to oxygen.
Litter in southeast Queensland is dominated by cigarette butts [28,46]. Currimundi Lake area attracts a large number of tourists because of the combination of a sheltered inlet in close proximity to a surfing beach. To reduce the impact of litter on the water quality, catchment community groups actively pursue litter clean ups. The strong direction coherence of the motion of drifters with wind during the closed state (E2) indicates that wind information can provide some guidance for litter management. During the closed condition, the flow is dominated by the wind and dispersion coefficients are an order of magnitude lower, resulting in less spreading of litter (see Figures 5,8 and 10) during the closed state (E2). The dominant wind direction for Lake Currimundi is SE-S and litter will therefore accumulate along the NW-N banks of the lake.

Summary and Conclusions
The dynamics and horizontal dispersion of Currimundi Lake, a shallow intermittently closed and open lake and lagoon, were studied to examine the dispersion coefficients and dispersion mechanisms for different inlet conditions and their implications on some key management issues. The analysis of drifter deployments, ADV observation, and physiochemical properties highlight a significant contribution of tide in the flow dynamics of the channel. A fixed ADV mounted near the surface and a GPS-tracked drifter during the open and closed states were used to observe the flow dynamics of the channel. The analysis of the water levels within the channel and on the ocean side indicated that the channel exhibited high and low tides that lagged the ocean within 2-3 h. Generally, the inflow and outflow velocities were weak with speeds <0.25 m s −1 . The channel exhibited a rapid flood tide flow compared to the ebb tide. This accelerated flow, combined with short-lived flood flow, may modify the upstream channel morphology. The channel was stratified during the closed system, likely linked with freshwater inflow, low wind forcing, and low solar insolation [6]. The freshwater pumped into the lake likely keeps the channel connected to the ocean [21] and also induced salinity gradient during the closed system. The results here further suggest that the freshwater input induced tidal asymmetry and enhanced horizontal dispersion within the channel.
Importantly, the closed state exhibited a dispersion coefficient up to two orders of magnitude less than those observed during the open state (E1), indicating that the absence of tide due to the blocked inlet reduced transport and horizontal mixing processes within the channel. During the open state (E1), it was found that the dispersion coefficient varied drastically for different tidal phases with the ebb tide average dispersion coefficient K~4.7 m 2 s −1 an order of magnitude greater than flood estimates.
The comparison of the observed dispersion coefficients with theoretical estimates showed that the tidal trapping effect and the tidal-induced horizontal velocity gradient likely explained the horizontal transport and dispersion during the open inlet state. During the closed state (E2), the drifters followed the wind driven surface flow with drifter speeds between 0-5% of the near water surface wind speed and aligned to the wind direction within 20 • . This indicated a modified Ekman surface transport during the closed inlet. This highlights that accurate wind information is required for modelling the transport of particles and the physiochemical properties of the channel during the closed inlet condition. An ongoing analysis examines the response of channel flow dynamics and the processes involved in restoration of water in the channel immediately after a mechanical opening operation.
Author Contributions: All authors (K.S., R.B., R.S., H.F., A.M., N.M and C.A) designed the experiment. K.S., R.B., R.S., H.F., A.M. and N.M. performed the experiments. K.S. analyzed the data and wrote the manuscript. All authors interpreted the data, discussed and edited the manuscript. In order to examine the measurement uncertainty associated with the key parameters used in this study, stationary tests using the GPS data-logger included in the low resolution drifters were carried out. The estimate of the centroid, velocity, dispersion coefficient and DKP associated with inherent errors due to GPS position fixing and hardware noise is examined. In order to estimate the magnitude of this errors on the estimates of the key analysis parameters, three independent observations were obtained at fixed locations using the GPS data logger in the low resolution drifters. Four dataloggers separated by~1 m apart were sampled continuously for 4−8 h. The dataset was quality controlled to synchronise the drifter position time series, lowpass filtered with cut-off frequency, F c = 0.01 Hz and analysed using relevant equations in Section 2.3. Table A1 above presents the 95% confidence interval around the mean estimates of the parameters using the measurements made at fixed points from three independent tests. For a given fixed cluster, theoretical values for centroid (with axes at 0,0), velocities, K c and DKP are zero if the GPS position fixing and hardware noise is zero. These measurements showed that while the inherent position accuracy of the drifters was between 2-3 m, the post-processed estimate of the centroid is in a decimeter accuracy while the velocity the centroid velocity is accurate up to, 10 −4 m s −1 . The estimates of the K c based on the measured noise were at least an order of magnitude lower than the values presented in Section 3. Correspondingly, magnitudes of DKP based on the noise were at least an order of magnitude lower than those reported in Section 3. A few exceptions to this condition of a minimum of order of magnitude different between inherent error and estimate of DKP, from Currimundi Lake, are flagged in results shown in Table 3 in the Section 3.3.