Hydrodynamic and Sediment Modelling within a Macro Tidal Estuary: Port Curtis Estuary, Australia

An understanding of sediment transport processes and resultant concentration dynamics in estuaries is of great importance to engineering design awareness and the management of these environments. Predictive modelling approaches provide an opportunity to investigate and address potential system responses to nominated events, changes, or conditions of interest, often on high temporal and spatial resolution scales. In this study, a three-dimensional hydrodynamic model and wave model were validated and applied to generate forcing conditions for input into a sediment transport model for the period 7 May 2010–30 October 2010 within a macro tidal estuary, Port Curtis estuary (Australia). The hydrodynamic model was verified against surface and near-bottom current measurements. The model accurately reproduced the variations of surface and near-bottom currents at both a mid-estuary and upper-estuary location. Sediment transport model predictions were performed under varying meteorological conditions and tidal forcing over a 180-day period and were validated against turbidity data collected at six stations within Port Curtis estuary. The sediment transport model was able to predict both the magnitudes of the turbidity levels and the modulation induced by the neap and spring tides and wind-wave variations. The model-predicted (converted) turbidity levels compared favourably with the measured surface water turbidity levels at all six stations. The study results have useful practical application for Port Curtis estuary, including providing predictive capabilities to support the selection of locations for monitoring/compliance sites.


Introduction
Estuarine environments are of great importance providing habitats (and nurseries) for numerous flora and fauna, including commercial and recreational fish species [1][2][3].Estuaries also play a significant role in altering the flow of contaminants and sediments between terrestrial and coastal waters [4], moderating water quality [5], absorbing wave energy [6], providing relief from tidal surges [7,8], and providing cultural and recreational benefits [6].Furthermore, the sheltered nature of estuaries often provides suitable settings for port and harbour infrastructure, which offer significant economic benefits in many parts of the world.However, the associated operations and expansions of large commercial ports and harbours often results in multiple environmental stresses, such as dredging to maintain navigable passage for larger ships, land reclamation, and changes in surrounding land use patterns altering sediment and nutrient runoff regimes [9][10][11].
An understanding of estuarine hydrodynamics and sediment transport dynamics is important since they play a critical role in the functionality and health of these systems.Knowledge of suspended sediment dynamics, influenced by tidal-and wind-driven currents and wave action [12][13][14], is essential for quantifying fluxes of substances and determining the fate of pollutants [15].Re-suspension of bottom sediments have been demonstrated to influence productivity in shallow-water environments [16][17][18] through water column enrichment by nutrients originating from the sediment [19] and impacting light availability [20,21].Additionally, re-suspension events potentially impact water quality through release of trace metals and organic contaminants [18,22], whilst also directly impacting on water column biota [23,24].The re-suspension and transport of sediments can cause increased sedimentation, which is often mitigated by expensive dredging operations.In addition to navigational issues sediment re-susepsnion and sedimentation events may also lead to the mortality of benthic communities through smothering and burial [25].
As such, an understanding of transport processes and concentration dynamics in estuaries is of great importance to engineering design awareness and the management of these environments [22,[26][27][28].Mechanisms controlling the re-suspension, and deposition of fine sediments within tide-dominated estuaries are complex [28,29].To date, considerable efforts have been invested to better understand the re-suspension mechanisms and transportation processes of sediments in estuarine settings (e.g., Ali et al. [29], Brand et al. [30], Giardino et al. [31], Guan et al. [32] and Hayter and Mehta [33]).The use and application of numerical models able to describe these physical processes (e.g., Carniello et al. [34], Cancino and Neves [35], Lumborg [36] and Le Normant [37]) are useful tools in planning and designing effective restoration and management strategies [38].Specifically, predictive models provide an opportunity to investigate and address potential system responses (e.g., suspended sediment concentrations and bed level changes) to nominated events, changes, or conditions of interest, often on high-resolution temporal and spatial scales.
Located between Keppel Bay and Rodds Bay, Port Curtis estuary is an area of significant importance adjacent the World heritage listed Great Barrier Reef Marine Park along the central Queensland coastline, Australia.In addition to being Australia's fifth largest multi-commodity port and the world's fifth largest coal export port in the world [39], the estuary is ecologically important for commercial and recreational fish species [40][41][42].Furthermore, the estuary also supports visiting and local populations of vulnerable dugong (Dugong dugon) and endangered green turtle (Chelonia mydas) [43].Activities surrounding Port Curtis estuary, including extensive mangrove and salt marsh reclamation for port infrastructure, marina, industrial and urban development, have altered the catchment, surrounding land use, and the estuarine conditions.According to the environmental and economic regional importance of the estuary and surrounding zones, the region has been the focus of recent scientific efforts [41][42][43][44][45][46][47].Previously, Herzfeld et al. [47] provided an overview of the circulation characteristics and distributions of dissolved material using a pilot stage three-dimensional hydrodynamic model of the Port Curtis estuary.However, the modelling required further calibration and validation in order to achieve full confidence in solutions.To date a predictive understanding of ambient sediment dynamics within the system is still lacking.
In this paper, we first present a three-dimensional hydrodynamic model, Delft3D-FLOW and a wave model, SWAN, followed by a description of a versatile tailored sediment transport model, DREDGEMAP, used as part of this study.We then present the application of the model for describing the suspended sediment dynamics for the case study region of Port Curtis estuary.Thereafter, we compare the sediment transport model predictions under varying meteorological conditions and tidal forcing over a 180-day period to turbidity data collected at six stations within Port Curtis estuary.This is followed by estuary-wide snapshots of total suspended sediment concentrations during example spring and neap flood tide events.

Study Site
Port Curtis estuary is a macro tidal estuary (Figure 1) with energetic physical forcing, including; spatially-and temporally-variable waves, tidal currents, winds, and freshwater inputs.The estuary covers an area of approximately 200 km 2 and includes a network of rivers, creeks, inlets, shoals, channels, and islands.Depths range from <1 m on the intertidal flats to >20 m within the central channel, and a mean depth of 5.3 m.Tides are semi-diurnal, with a maximum tidal range of 4.69 m at Gladstone and up to 6.00 m within the Narrows [41].The tide propagates into the estuary through the straits separating Facing Island from the mainland (Gatcombe Channel) and Curtis Island (North Channel), and through the Narrows via Keppel Bay in the north (Figure 1).Tides undergo a spring-neap cycle with a period of approximately 14 days, with maximum ranges during spring tides and ~1 m during neap tides.The strongest tidal currents are focused in the main tidal channels.Though a less frequent physical forcing mechanism compared to tidal forcing, which causes most of the estuarine mixing, freshwater inputs are occasionally important during periods of high rainfall (i.e., summer months).Freshwater flows originate from the Narrows as a result of the nearby Fitzroy River flooding, the Calliope River and Auckland Creek which discharges into the estuary through Gladstone (Figure 1).The large tides ensure that the water column is vertically well mixed, and are also responsible for significant re-suspension of fine sediment [47].Combined with very large deposits of silt from the hinterland in times of flood, the estuary maintains a highly turbid character.The estuary is characterised by extensive areas of tidal flats that become exposed at low tide and large areas of mangroves fringing the estuary, which behave as a storage buffer for water at high tide.Port Curtis estuary sediments consists primarily of silts and clays within the shallow intertidal banks, while fine and course sand predominant the more tidally energetic deeper regions of the estuary.Additionally, the estuary is influenced by local wind-driven waves, which contribute to turbulences and sediment suspension events.Wave activity is greatest during periods of prolonged south-easterly trade winds.

Model Description
A hydrodynamic model and wave model in conjunction with a sediment transport model were used during this study.Firstly, a hydrodynamic and wave model was used to simulate complex current flows and wave actions throughout the estuary subject to tidal flows and prevailing metocean conditions within Port Curtis estuary.Thereafter a versatile tailored sediment transport model, DREDGEMAP [48], was coupled with the hydrodynamic and wave models over a 180-day duration to model suspended sediments throughout Port Curtis estuary.

Hydrodynamic Model
The hydrodynamic conditions were simulated using Delft3D-FLOW [49].The model has been successfully applied within oceanic, coastal, estuarine, riverine, and flooding conditions [49][50][51].A curvilinear computational grid with varying cell size was developed for the entire extent of the estuary.The finest grid resolution (25 m) was defined along shorelines and regions of steep depth gradients, to accurately define the current patterns.Whilst the coarsest resolution (~400 m) was applied in the outer region nearby the southern entrances of the estuary and neighbouring coastal waters (Figure 1).Datasets used to represent the bathymetry of the estuary were sourced from commercial nautical charts and bathymetric surveys [52].In order to define the three-dimensional current structure, a σ-layer approach was employed with five equally-proportioned vertical layers throughout the model domain.Tidal variations were driven at the open boundaries of the model domain through the northern entrance using water elevation data from Port Alma tidal station (23.58° S, 150.87°E) and the eastern entrances using data from the Gatcombe Head tidal station (23.88° S, 151.37°E).
The "cyclic" advection schemes for momentum and transport were selected.The κ-ε turbulence closure scheme was selected for the vertical viscosity formulation, with background values set to 0 m 2 /s.The Manning roughness formula was used to account for the bottom friction, with a uniform coefficient of 0.02 selected as representative of the seabed in the region on the basis of previous works by Wolanski et al. [53] and King and Wolanski [54].The horizontal eddy viscosity was set to 1 m 2 /s throughout the model domain.Hourly wind speed and direction data, measured at the Gladstone Airport by the Australian Bureau of Meteorology (station identification: 039326; 23.87° S, 151.22°E), was used to account for the wind forcing across the model domain.The wind drag coefficients, CD, is assumed to vary linearly with wind speed from 0.00063 to 0.00723 over the range of U10 = 0-100 m/s.The model time-step was set to 0.2 min, which produced Courant numbers less than 10.

Wave Model
The SWAN model [55] was used to reproduce the prevailing wave conditions.The model is fully spectral (in all directions and frequencies) and computes the evolution of wind waves in coastal regions with shallow water depths.The physical processes selected for the simulations were: white-capping, depth-induced wave-breaking, bottom friction, and triad wave-wave interactions.Using an unstructured triangular mesh approach, the generated computational grid included the entire Port Curtis estuary and offshore regions (24.00° S, 151.50 °E to 23.65° S, 151.50°E, Figure 1).A fine (<50 m) mesh resolution was defined for tidal creeks and shorelines and a coarser mesh resolution (>500 m) was applied over the estuary region and near the open boundaries of the domain.The bathymetric dataset used for the wave model was sourced, as per the hydrodynamic modelling.Surface waves entering the model boundaries were obtained from the WAVEWATCH III ® (WW3) global model (National Centres for Environmental Prediction National Oceanic and Atmospheric Administration, College Park, MD, USA) [56].Outputs of WW3 (significant wave height (Hs), mean wave period (Tp) and peak wave direction period (Pdir)) were extracted from the computational point closest to the area of interest (24.00° S, 152.50°E) and used for the model simulations.Model surface and open boundary data were updated at hourly intervals.The default SWAN parameterisation of depth-induced wave-breaking was employed (0.73 [57]).Bottom friction was activated using the drag law-based model of Collins [58] and a bottom friction coefficient of 0.01 was applied.In addition, quadruplet wave-wave interaction was also included.Hourly wind speed and direction data, measured at the Gladstone Airport by the Australian Bureau of Meteorology, was used as model input.

Sediment Transport Model
Suspended sediment concentrations (SSC) were modelled using a tailored version of DREDGEMAP, an enhancement of SSFATE [48], a three-dimensional sediment transport and fates model, jointly developed by Applied Science Associates Pty Ltd and the U.S. Army Corps of Engineers (USACE) Environmental Research and Development Center (ERDC).The model has been successfully applied in coastal and estuarine settings [59][60][61].The model represents the total mass of sediments suspended over time by a defined sub-sample of Lagrangian particles, allocating an equal proportion of the mass to each particle, where the transport, dispersion and settling of suspended sediment released to the water column is calculated using a random-walk procedure.Sediment particles are divided into five size classes.For the purpose of this study the following grain size classes were adopted: 5 μm, 75 μm, 200 μm, 450 μm and 1000 μm.
Settling of mixtures of particles is a complex process due to interaction of the different size classes, some of which tend to be cohesive, and thus clump together to form larger particles that have different fall rates than would be expected from their individual sizes.Enhanced settlement rates, due to flocculation and scavenging, are particularly important for clay and fine-silt sized particles and as such these processes are implemented in DREDGEMAP.
Minimum sinking rates are calculated using Stokes equations, based on the size and density of the particle.However, sinking rates of finer classes (representing clay and silt-sized particles) are increased based on the local concentration of the same and larger particles, to account for clumping and entrainment.
The settling velocity of each particle size class ( ), is computed from: where: Culi and Clli are the nominal upper and lower concentration limits, respectively, for enhanced settling of grain size i, ai is a grain-size class average maximum floc settling velocity, C is the total concentration for all grain size classes [61]. ( Selected values for Culi, Clli and ai adopted for this study ranged between 1000-8000 mg/L, 50-400 mg/L and 0.0008-0.1 m/s, respectively.The settling velocity for the largest size class (1000 μm) was assumed constant at 0.1 m/s.
Deposition is calculated as a probability function of the prevailing bottom stress, local sediment concentration and size class.Mixing of re-suspended sediment into the water column is a dynamic balance between estimates of the sinking rate and vertical mixing induced by turbulence (as specified by vertical mixing coefficients).
Sediment deposition flux is computed as: If 0 0.05, then 1 Otherwise, where: Pi is the deposition probability of i-th size class, Ci is the sediment concentration, Wsi is the calculated settling velocity and bi is the empirical parameter that includes all other factors influencing deposition other than shear (0.2-1.0 for clays to fine sand) [61].Deposition probability, Pi, is calculated as follows: P1, for clay sized particles (1st size class) where: τcd is the critical shear stress for deposition for the clay fraction, for this study a value of 0.016 was adopted.The bottom shear stress is based on the combined velocity due to current and wave action using the parametric approximation by Soulsby [62].
Pi, for the other size classes (2nd, 3rd and 4th size classes) where: τuli is the shear stress above which no deposition occurs for grain size class i, and τlli is the shear stress below which deposition probability for grain class i is 1.0.For values between τlli and τuli, linear interpolation was used.Adopted τuli and τlli values for this study ranged between 0.030-0.900and 0.016-0.200,respectively, excluding the 1000 μm size class.
The model can employ two different re-suspension algorithms.The first is based on the work of Sanford and Maa [63] and, subsequently, Lin et al. [64], and applies to material deposited in the last tidal cycle, which accounts for the fact that newly deposited material will not have had time to consolidate and will be re-suspended with less effort (lower shear force) than consolidated bottom material.The second algorithm, which was used in this study, is the established Van Rijn method [65] and applies to all other material that has been deposited prior to the start of the last tidal cycle.This method calculates a constantly varying critical threshold for re-suspension, based on the median local particle-size distribution for settled material.Re-suspension algorithms implemented in DREDGEMAP are presented in detail by Swanson et al. [61].
Knowledge of the initial bed sediment distribution was critical in predicting suspended sediment concentrations.Surface sediment grain size contribution maps were established based on data sourced from the analysis of >650 sediment cores collected throughout the estuary [66].The contribution of five sediment size classes (5 μm, 75 μm, 200 μm, 450 μm and 1000 μm) throughout the estuary was spatially assigned and interpolated using the inverse distance weighting method [67] throughout the model domain.A grid resolution of 40 m was used throughout the model domain.Clays and silts dominate the shallow intertidal regions of the estuary, whilst larger sediments, including fine and coarse sands dominate the more energetic deeper channel regions of the estuary.Figure 2 presents the median grain size (D50) of the surface sediments of Port Curtis estuary used as part of the study.Measured sediment densities were incorporated into the model inputs.The sheltering effects of seagrass on shear stress [68] and biologically-altered erodability of the bed sediments [69,70] was not included as part of the modelling approach.Model predictions output as SSC were converted to turbidity in order to provide comparison with coincident turbidity datasets collected at six stations.Suspended sediment concentrations were converted as follows [71]:

Model Validation
Measured data was used to validate the models and to verify their ability to predict prevailing conditions.

Hydrodynamic Model
Model-predicted currents were verified using measured Acoustic Doppler Current Profiler (ADCP; seafloor mounted Teledyne RD ADCPs) data collected between 7th to 27th October 2010 at a mid-estuary (23.78833° S, 151.20142°E; depth ~7.5 m) and upper-estuary site (23.74817°S, 151.16243°E; ~10.0 m) (Figure 1).The model accurately represented the semi-diurnal mixed spring and neap tides throughout the 21-day deployment period, as well as variations of the surface and near-bottom current speeds and directions during the deployment period.Analysis of the measured data revealed greater current speeds in the upper layers of the water column in comparison to the lower depth layers at both sites with maximum surface current speeds exceeding 1 m/s and near-bottom speeds approaching maximum speeds of 0.8 m/s.Predominant flow directions were northwest and southeast, which corresponds to the orientation of the adjacent coastline, occurring during the flood and ebb tide, respectively.Directional differences were observed between depth layers during the time-series as a result of wind influences and bathymetry-induced steering.To provide a quantitative assessment of the performance of the hydrodynamic model, the correlation coefficient (R), Index of Agreement (IOA; Willmott [72]) and Mean Absolute Error (MAE; Willmott [73]) were calculated for surface and near-seabed currents.
Comparison between the measured and predicted surface elevation at both the mid-estuary and upper-estuary sites demonstrated very good agreements (Figure 3).Statistical analyses between the measured and predicted elevation datasets produced R, IOA, and MAE values of 0.99, 0.99, and 0.17 at the mid-estuary site and 0.99, 0.98, and 0.25 at the upper-estuary site, respectively.Additionally, the current speeds and direction were shown to be satisfactorily reproduced with IOA values for current speed ranging between 0.82-0.90and direction ranging between 0.93-0.97for both sites and depth layers compared.Table 1 provides a summary of the results of the statistical analysis performed.The performance of the model indicates the suitability for input into the sediment transport model.

Wave Model
In the absence of inshore wave data, the model accuracy was verified using measured wave data between December 2009 to March 2010 offshore Port Curtis estuary (23.89562° S, 151.50238°E [74]).Table 2 provides a summary of the results of the statistical analysis performed using the measured and predicted wave parameters offshore Port Curtis estuary.The model was shown capable of replicating the general trend at the offshore wave buoy.Predicted long-period ocean waves ranged from ~1 to 4 m in height, with the predicted peak wave direction being ~90°.Due to the orientation of the estuary entrance, long-period ocean swells are typically blocked by Facing Island and the mainland coastline south of Port Curtis and, as such, ocean wave penetration is negligible.The wave climate within the estuary is, therefore, typically dominated by locally generated wind waves.The maximum fetch within the estuary (~20 km) is achieved by southeast winds, which are one of the predominant wind directions of the region [75].Significant wave heights within the estuary have been documented as being <0.5 m and <0.3, ~96% and ~80% of the time, respectively, with wave periods less than 5 s [76].The modelled estuary wave climate demonstrated Hs typically ranging between <0.01-0.5 m with an average Hs value of ~0.2 m within the mid-estuary region.The mean wave period typically ranged between 1 and 4 s within the model domain.

Sediment Transport Model
The sediment transport model was calibrated and compared against turbidity surface water (~0.5 m depth from surface) measurements collected at six different stations within Port Curtis estuary (see Figure 1 and Table 3).Stations were deployed and maintained as part of a Port Curtis estuary monitoring program using YSI 6920 sondes with 650MDS data loggers (YSI Australia, Brisbane) e.g., [77].The modelling considers an extended period characterised by different tidal and meteorological conditions over a 180-day period.In order to illustrate and quantify the model's performance, the results of the simulations were compared to the measured data, based on both hourly-measured datasets (Figure 4) and based on an exponentially-weighted moving average (EWMA) (6 h) (Figure 5).The model-predicted turbidity levels compare favourably with the measured data.The model predicted both the magnitudes of the turbidity levels and the modulation induced by the neap and spring tides and wind-wave variations.As the six turbidity stations were well distributed within the estuary and positioned in locations with varying initial surface sediment grain size composition, water depth, bottom topography, and wind wave fetches, the comparison between the model-predicted and measured turbidity levels is considered a rigorous assessment.Where HAT represents highest astronomical tide.
Model performance was assessed by determining R, IOA, and MAE values at each station based on the hourly measured and predicted turbidity datasets.Table 4 provides a summary of statistical analysis performed using the measured and predicted turbidity hourly datasets at the six stations within Port Curtis estuary.Statistical analyses between the measured and predicted hourly and EWMA datasets produced R values ranging between 0.07-66 and 0.17-0.75,IOA between 0.22-0.68 and 0.30-0.83,and MAE ranging between 7.07-12.57and 4.96-8.14,respectively.Generally the model performed best in the mid-estuary region.Once suspended from the estuary floor, the fate of the suspended sediments is dependent on factors, including particle size and associated settling velocity, ambient current speed and direction, and water depth.Larger particles typically located within the deeper more energetic regions settle more quickly once re-suspended, while the finer particles that are re-suspended remain suspended and are transported greater distances and settle during decreased current speeds and wave action.
Although the model-predicted turbidity levels compared favourably with the measured surface data, the model does possess limitations in terms of included turbidity sources.Areas identified that require further inclusion into the transport model include; sediment loadings from freshwater inflow, biological production of suspended particulate matter owing to algal growth, bioturbation, and inclusion of seagrass/biological influence on bed shear and re-suspension.Furthermore, dredging and ship traffic-induced sediment re-suspension, which may also affect suspended sediment dynamics, was not considered as part of the model.Despite the exclusion of contributing turbidity sources associated with the model, the modelling illustrates the importance of the tidal-induced currents and wind-wave interactions on sediment re-suspension and the applicability of the model in predicting turbidity levels within Port Curtis estuary.Reduced discrepancies between the measured and predicted turbidity levels may be achieved following the model adoption of established site-specific relationships (conversions) between measured SSC and turbidity levels.Such an understanding would prove very important for the model's performance when reporting turbidity levels and also provides flexibility in the model's application.

Hydrodynamic Conditions
The verified hydrodynamic model (and wave model) was used to generate current data as inputs into the sediment transport model.The modelled period was selected to coincide with data made available from the Western Basin Dredging and Disposal Project water quality monitoring program undertaken.Surface current speeds were greater than near-bottom currents throughout the estuary (Table 5).The maximum and mean current speeds occurred during spring periods, with maximum current speeds occurring in the deeper channels centrally located within the estuary.Reduced current speeds are experienced in the shallower intertidal regions of the estuary.Mean and maximum surface current speeds during spring tides within the mid-estuary region were 0.45 ± 0.27 m/s and 1.25 m/s, respectively, and 0.38 ± 0.21 m/s and 1.02 m/s, respectively, in the upper-estuary.Whilst near-bottom mean and maximum current speeds within the mid-estuary and upper-estuary regions were 0.35 ± 0.19 m/s and 0.77 m/s, and 0.31 ± 0.15 m/s and 0.81 m/s, respectively.In comparison, the mean and maximum surface current speeds occurring during neap tides within the mid-estuary and upper-estuary regions were 0.35 ± 0.23 m/s and 1.13 m/s, and 0.26 ± 0.17 m/s and 0.83 m/s, respectively.Near-bottom mean and maximum current speeds in the mid-estuary and upper-estuary regions were 0.26 ± 0.14 m/s and 0.61 m/s, and 0.21 ± 0.12 m/s and 0.62 m/s, respectively.
Current directions predominantly occurred along the northwest-southeast axis, which corresponds to the orientation of the adjacent coastline (see Figure 3).The northwest and southeast currents occurred during the flood and ebb tides, respectively.In lesser instances, easterly flow components were modelled (and observed) within the estuary, for example within the mid-estuary region when ebb waters drain from South Passage Island (northwest of Station 3).Surface and near-bottom currents also demonstrated divergences in current angle.Vertical differences in current angle and speed are presumably a combined result of the wind influence on surface currents and bathymetry (deep channels amongst shallow and inter-tidal areas).Specifically, the water movement in the shallows (<5 m) would experience less steering than the movement of water in regions >5 m that is steered by the channel walls.Due to this factor there are, seemingly, instances where surface and deeper layers decouple due to the difference in momentum.Furthermore, the momentum influence is enhanced during spring tides characterised by stronger currents and greater tidal prism.The vertical variability support the selection of the three dimensional hydrodynamic model used as part of this study in capturing the depth varying current speeds and directions.
It is acknowledged that temporal extremes in flood events and prolonged periods of increased wind events (e.g., cyclonic weather systems) significantly influence current dynamics which, although irregular, should also be considered to provide a more robust approach for future predictions, especially given the importance of these events in delivering suspended material into the estuary and their influence on the regional hydrodynamics throughout Port Curtis.

Spatial and Temporal Turbidity Variations
Modelling of SSC, converted to turbidity, was undertaken to establish a predictive tool in forecasting suspended sediment concentrations (and/or turbidity levels) throughout Port Curtis estuary.The estuary typically maintains a highly turbid character [47] and during the modelling period (May-November 2010) measured turbidity levels demonstrated temporal and spatial variations, attributable to differences and interactions between factors including strength of tidal currents, water depth, bathymetry variation, proximity to the estuary entrance, grain size, and composition of benthic and suspended sediment, and fetch.Suspended sediment concentrations and turbidity levels within the estuary have been shown to demonstrate variations/trends according to wet and dry seasonal periods with greater values occurring during the wet season (i.e., December to April) [77].
Mean predicted and measured surface water turbidity values within Port Curtis ranged from 4.31 ± 6.74 to 13.60 ± 12.80 NTU and 7.48 ± 15.10 to 12.62 ± 11.80 NTU, respectively.Predicted turbidity 50th and 95th percentiles ranged from 1.92 to 8.68 and 17.99 to 42.50, respectively.Surface turbidity is summarised in Table 6.Where SD represents standard deviation and % RSD represents % relative standard deviation.
A clear spatial and temporal trend was evident as demonstrated by the datasets.Generally, turbidity was greatest in the mid-and upper-estuary stations (i.e., 1-5) in comparison to the lower estuary station (6), which was positioned closest to the estuary entrance.Turbidity demonstrated high variability, particularly in the mid-and upper-estuary stations, with fluctuations in turbidity predominantly influenced by spring-neap tide variations.Turbidity was greatest throughout the estuary during spring tides due to the increased tidal energy input during this period, where peak current speeds are greatest.On a tidal cycle time-scale, greatest turbidity typically occurred during the flood tide period.To a lesser extent periods of elevated wind speeds were also shown to influence increased turbidity.
Predicted and measured turbidity levels at all sites exceeded both the Central Coast Queensland Water Quality Guideline (QWQG) regional turbidity guideline for mid-estuary and tidal canals, constructed estuaries, marinas and boat harbours (8 NTU) and upper estuarine waters (25 NTU) [78], and the Australian Water Quality Guideline (AWQG) turbidity guideline for slightly disturbed ecosystems in tropical Australia estuarine and marine waters (20 NTU) [79].
Given the good performance provided by the model in reproducing the measured turbidity levels at the six different stations and under different current and meteorological conditions within Port Curtis estuary we provide model-predicted "snapshots" of surface water turbidity (0-1 m depth) throughout the estuary.Snapshots are provided for spring and neap tide flood periods, providing example comparisons of spring and neap tide conditions (Figures 6 and 7).The snapshot images illustrate the time-and space-varying dynamics within the estuary, under sample spring and neap tide flood periods.Greatest predicted surface suspended sediment concentrations were associated within regions characterised by increased current speeds exceeding the bed shear stress, inducing suspensions and subsequent transportation of suspended sediments.Shallower intertidal regions characterised by reduced current speeds typically demonstrated decreased suspended sediment concentrations.
The results from the modelling have useful practical applications for the Port Curtis estuary, including providing predictive capabilities used in parallel with modelling of fugitive sediments from operations such as dredging and supporting selecting locations of monitoring/compliance sites.Furthermore, modelling approaches are recognised as a valuable tool for understanding and predicting sediment deposition regimes and morphological developments.Predicted zones of sediment deposition and erosion were well represented within the model during this study (model output not shown).
The model-predicted turbidity levels compared favourably with the measured surface water turbidity levels despite the exclusion of known contributing turbidity factors.The model demonstrates an ability to sufficiently predict surface water conditions (i.e., suspended sediment concentrations and/or turbidity), which may be used as a standalone predictive tool, or used in a complimentary manner with event driven, source-specific modelling (e.g., dredge plume modelling).The combined model approach can be used to provide a greater understanding of future suspended sediment conditions and allow a better understanding and flexibility for compliance based operations.
The model presented clearly represents a simplified approach which does not produce all major terms in the sediment transport budget.The inclusion of such identified contributing factors present an opportunity to further refine the operational capability of the modelling approach within the Port Curtis estuary.Additionally, the inclusion of appropriate modelling routines representing light extinction presents an opportunity to broaden the scope of the model application from an ecological management standpoint by providing photosynthetic active radiation model predictions.

Conclusions
A three-dimensional hydrodynamic model and wave model coupled with a sediment transport model was validated and used to investigate turbidity (suspended sediment) dynamics within a macro tidal estuary, Port Curtis.The hydrodynamic model was verified against surface and near-bottom water column measurements.In general, the simulation results accurately reproduced the variations of surface and near-bottom currents at both a mid-estuary and upper-estuary location.The three-dimensional hydrodynamic model was applied to generate forcing currents for input into the sediment transport model for the period 07 May 2010-30 October 2010.The sediment transport model predicted both the magnitudes of the turbidity levels and the modulation induced by the neap and spring tide tidal phase and wind-wave variations.The model-predicted (converted) turbidity levels compared favourably with the measured surface water turbidity levels.Predicted and measured turbidity levels at all sites exceeded both the QWQG regional turbidity guideline for mid-estuary and tidal canals, constructed estuaries, marinas and boat harbours (8 NTU) and upper estuarine waters (25 NTU), and AWQG turbidity guideline for slightly disturbed ecosystems in tropical Australia estuarine and marine waters (20 NTU).Despite the exclusion of contributing turbidity sources associated with the model, the predicted turbidity (suspended sediment) levels illustrates the importance of the tidal-induced currents and wind-wave interactions on sediment re-suspension and the applicability of the model in predicting turbidity levels within Port Curtis estuary.The inclusion of such identified contributing factors present an opportunity to further refine the operational capability of the modelling approach within the Port Curtis estuary.Following model validation the sediment transport model was used to model the spatial and temporal variability of suspended sediment concentrations within Port Curtis estuary.The results of modelling have useful practical application for the Port Curtis estuary, including providing predictive capabilities which may be used in parallel with modelling of fugitive sediments from operations such as dredging and decision support for selecting the locations for monitoring/compliance sites.

Figure 1 .
Figure 1.Site study map showing location of Port Curtis estuary indicating the model domain extents, ADCP site locations and surface waters turbidity measurement stations.

Figure 3 .
Figure 3.Comparison between measured and predicted near-bottom current speeds (a,c) and direction (b,d) at the mid-estuary and upper-estuary ADCP sites.The data is compiled from the measured and predicted datasets 07 October 2010-27 October 2010.

Figure 4 .
Figure 4. Comparison between the measured and predicted surface turbidity based on turbidity hourly datasets at the six stations (Stations 1-6).The data is compiled from the measured and predicted datasets 07 May 2010-30 October 2010.

Figure 5 .
Figure 5.Comparison between the measured and predicted surface turbidity based on 6 h measured datasets.The data is compiled from exponentially-weighted moving average turbidity hourly datasets at six stations (Stations 1-6) from the measured and predicted datasets 07 May 2010-30 October 2010.

Figure 6 .
Figure 6.Model-predicted "snapshot" of surface water (0-1 m depth) suspended sediment concentrations during 02:00-08:00 from a spring tides sample flood period during 08 August 2010.The cross-section inset depicts the predicted water column suspended sediment concentrations along the transect (-----) between points A and B.

Figure 7 .
Figure 7. Model-predicted "snapshot" of surface water (0-1 m depth) suspended sediment concentrations during 00:00-06:00 from a neap tides sample flood period during 20 August 2010.The cross-section inset depicts the predicted water column suspended sediment concentrations along the transect (-----) between points A and B.

Table 1 .
Hydrodynamic model performance summary.The data is compiled from the measured and predicted datasets 07 October 2010-27 October 2010.

Table 2 .
Wave model performance summary.The data is compiled from the measured and predicted datasets 01 October 2009-30 March 2010.

Table 3 .
Station summary of measured turbidity dataset used as part of sediment transport model validation.

Table 4 .
Sediment transport model performance summary based on the comparison between measured and predicted hourly turbidity datasets and exponentially weighted moving average (6 h).The data is compiled from the measured and predicted datasets 07 May 2010-30 October 2010.Location of measurement stations are shown in Figure1.

Table 5 .
Model-predicted surface and near-bottom current speeds based on hourly datasets.Measured dataset shown in square brackets for the mid-estuary and upper-estuary ADCP.The data is compiled from the measured and predicted datasets 07 May 2010-30 October 2010.

Table 6 .
Model-predicted turbidity summary based on hourly datasets.Measured dataset shown in square brackets.The data is compiled from the measured and predicted datasets 07 May 2010-30 October 2010.