Towards a Tidal Farm in Banks Strait, Tasmania: Inﬂuence of Tidal Array on Hydrodynamics

: The development of tidal energy in Australia is still a challenge with few studies performed on the characterisation of the resource, due to the difﬁculty to acquire data and uncertainties about the inﬂuence of this anthropogenic activity on the marine environment. Changes in ﬂow could lead to alterations in sediment transport and have further inﬂuence on the marine habitat. A case study in a promising area, Banks Strait (Australia), was created using high resolution 2D and 3D models validated against in situ data to investigate changes to hydrodynamic conditions with two scenarios of tidal farms (100 and 300 turbines). Comparison between 2D and 3D is performed to ﬁnd the best compromise between model accuracy and computational time for preliminary assessment. Changes to current speed and bed shear stress over a 35 day period were found to be localised around the tidal farms and did not extent more than 7 km from the farm (300 turbines) for both 2D and 3D. The results showed that for near ﬁeld and far ﬁeld, 2D models are sufﬁcient to give a ﬁrst approximation of the hydrodynamic inﬂuence of tidal farm deployment on its environment.


Introduction
The development of tidal industry has made significant progress over the last decades: numerous tidal projects are under way, examples include development in United Kingdom (Pentland Firth, Shetland, Orkney), in France (Ushant, Aldernay Race) and in China (Zhoustan Archipelago) [1]. There also have been considerable advances in numerical models in the recent years, to assess the influence of tidal farms on the environment principally in 2D [2][3][4][5][6][7][8][9][10][11][12]. However, the conclusions of all these studies are site-specific and dependent on device design; thus, the effects of array size will need to be investigated in the site of interest. Generally, 2D models are used for regional modelling while 3D models are used to assess the local tidal energy resource. Nowadays researchers tend to use 3D models to assess the resource and the environmental impact, although at higher computational costs [13].
The objectives for 2030 for Australia (reducing CO 2 emissions to 26-28% below 2005 levels) and Tasmania's plan to become 200% in renewables by 2040 are a major challenge. To meet these obligations, the Australian government is considering to speed up the deployment of ocean energy, looking towards blue economy projects to reduce fossil fuels reliance. Fourteen Australian tidal power trials have taken place over the last two decades, they are summarised in Table 1 [1,14]. To assess accurately tidal energy sites, the hydrodynamic conditions that produce effective tidal streams must be determined. Further high resolution projects are required to perform a detailed site assessments in order to evaluate the viable production, environmental influence of Australia's tidal resource and tidal farm development. The Australian Tidal Energy (AUSTEn) project [15] funded by the Australia Renewable Energy Agency (ARENA) aims to map Australian tidal energy resource with three components. One includes the characterization of two promising sites (the Banks Strait in Tasmania and the Clarence Strait in Northern Australia) with a focus on the energy resource potential, the influence on physical parameters and the development of fine-scale models. Banks Strait (Figure 1), a strait located in the north east of Tasmania between Swan Island and Clarke Island, has been identified as suitable for commercial tidal farm developments [16]. This 15 km wide strait connects the Tasman Sea to the south east of Bass strait: which due to its strong currents and proximity to grid-connection is a very attractive site for the development of the tidal arrays in Australia. Preliminary coarse resolution hydrodynamic modelling by the Commonwealth Scientific and Industrial Research Organisation (CSIRO) indicated strong currents of around 2.5 m/s in this area [17] and more recently [18] found an area of 9 km 2 south of Clarke Island with maximum current speed higher than 2 m/s. The surveys performed by the AUSTEn project in Banks Strait in 2018 revealed important coastal features including the presence of large sand dunes near Clarke Island (Figure 1). The deployment of a tidal farm in this area may influence these natural phenomenons.
To assess a tidal energy site, the different steps are the collection of data with surveys, then the numerical modelling of the site to perform Environmental Impact Assessment (EIA). Previously in the Banks Strait region, field measurements such as bathymetry, current velocities, sub-bottom properties and suspended sediment concentration were sparse due to the energetic environment found in this area. Recently the AUSTEn campaigns filled this gap using a large set of instruments and surveys techniques. Most of studies conducted in this region [19][20][21][22][23][24][25][26] have focussed on the current flow through the Bass strait, without focussing specifically on Banks Strait region. This strait was not modelled with a resolution high enough to be reliable for accurate tidal energy site study until recently with new work found in the Banks Strait area: [18] used a 400 m mesh resolution model which was validated against an Acoustic Wave and Current Profiler (AWAC) field data with no implementation of turbines. More fine-scale modelling is therefore needed in Banks Strait in order to reduce the uncertainty of the influence of tidal energy extraction on the physical environment. To reduce uncertainty in sediment transport rates, hydrodynamic models should be calibrated to less than 10% errors in the velocity [27,28]. Recent literature and the comprehensive review of the current understanding of the hydro-environmental impacts [29] revealed the next challenges for tidal energy: more complex physics (3D), models coupling wave and current, very far field effect and one of the most important: field measurements, critical to validate model simulations. The aim of this work was to investigate the influence of two tidal farms on physical parameters in Banks Strait using high-resolution 2D and 3D hydrodynamic ocean models. Highly calibrated unstructured mesh hydrodynamic models were validated against a large set of in situ data with implementation of two tidal farms (100 and 300 turbines). Investigation on whether 2D or 3D model dimensions are required for EIA in the near and far field are also presented. Details of the Banks Strait site and the AUSTEn surveys are presented in Section 2, with model validation results presented in Section 3 using field data from the AUSTEn campaigns. The influence in oceanographic systems of two simulated tidal farms is then described in Section 4. This study is also a solid basis for future work to determine the influence of tidal farms on sediment dynamics.

Modelling Methodology
This section outlines the description of the site, field data, numerical model, forcing inputs and validation process.

The Site and Data Collection
Banks Strait is located in the eastern part of Bass Strait between Clarke Island and the mainland of Tasmania. Principal harmonic constituents in the strait are the semidiurnal M2, S2 and diurnal O1, K1 with the semidiurnal M2 tide being the principal constituent [20,25]. The water mass characteristics of Banks Strait are altered by tidal stirring and local air-sea fluxes: with strong warming in summer and strong cooling in winter [26] and with saltiest water in winter [21]. Currents in the strait are primarily driven by tidal forcing, winds and shelf waves in Bass Strait. In Banks Strait, Sandery [21] described the area as well mixed because of the turbulent system due to strong tidal currents occurring. Sandery [20] also mentioned inferred flow in Banks Strait in winter with the sub-Antarctic Surface Water (SASW) intrusion.
The region of Banks Strait has been surveyed in March, July and December 2018, with the AUSTEn project [15], allowing measurements of seasonal variations. Field measurements at up to five stations using Acoustic Doppler Current Profiler (ADCP) were performed in Banks Strait in March, three in July and one in December ( Figure 2). The details of ADCP deployments are described in Table 2. These instruments were configured to measure currents, to calibrate and validate the numerical model and were processed using the Integrated Marine Observing System (IMOS) Toolbox [30].  Figure 2. Field measurements of conductivity Temperature Depth (CTD) casts found vertical fluctuations of less than 5% [31], confirming the well mixed area in this energetic system and allowing the use of barotropic ocean models without any reduction in simulation accuracy. Data from the AWACS found average wave heights between 1 and 2 m with maximum values reaching over 5 m [32]. Studies of the region found currents higher than 2 m/s during ebb tides in the middle of the channel with a tidal asymmetry of 20 degrees between ebb and flood at the location of the station CW4 [18]. Maximum power density over 3 kW/m 2 was found in the middle of the channel, exceeding the 2.5 KW/m 2 threshold outlined for the economic viability of a site [33].

Model Domain and Forcing
High-resolution numerical models of Banks Strait were developed in 2D and 3D using the unstructured mesh MIKE21 and MIKE3 FM modelling software (DHI A/S, Hørsholm, Denmark). The software solves Reynolds-averaged Navier-Stokes (RANS) equations using the Boussinesq assumption [34]. The numerical domain extends between 38 • S and 44 • S and 141 • E and 150 • E. The mesh resolution varied in the domain from an average of 3 km at the boundary to 100 m at the approach of Banks Strait, in the strait the typical node spacing was 50 m and increased to 20 m in the area suitable for tidal farm installation and key interest features such as the sand dunes near Clarke Island (Figure 3). The model shoreline data were extracted from the Australian Bathymetry and Topography grid [35] and from Australian Hydrography Service AusEnc charts (AHS) [36] for better resolution in the Banks Strait area (north east coast of Tasmania and Furneaux Islands). Bathymetry data were also collected from Geoscience Australia (GA 2009) [35] for the larger domain, from the AHS for more precise features in the Banks Strait area (between the mainland of Tasmania and Flinders Island) and from the AUSTEn Project (multi-beam in March 2018). The collected bathymetry dataset was interpolated onto an unstructured mesh using the natural neighbour method. Three open boundaries were forced with tidal constituents (surface elevation) from the Global tide MIKE model [37]. This tide model represents the major diurnal (S1, K1, O1, P1 and Q1), semidiurnal (M2, S2, N2 and K2) and shallow water/nonlinear (M4) constituents with a spatial resolution of 0.125°× 0.125°based on altimetry data from TOPEX/POSEIDON. The open boundaries were set far beyond Banks Strait to reduce their influence on simulation results as recommended [7,38]. The ERA5 climate reanalysis dataset [39] was used for wind components, with the collated dataset mean sea level and pressure closely matching the values from the Bureau of Meteorology measurements [40] in the Banks Strait at Swan Island (R 2 = 0.9) when compared to the year 2018.

Model Settings
The barotropic models were run for 42 days to follow IEC recommendations [38], with 1 week spin up (to reach a steady state), from 16/03/2018 00:00 to 26/04/2018 00:00, generating a 35-day dataset. The Smagorinsky formulation was used for the horizontal eddy viscosity and the log-law formulation for vertical eddy viscosity. Changes in 3D compared to the 2D setup included the vertical layers and the bed resistance where the roughness height was used and calculated from Manning's number (M). A Manning value of 32 was found to give the best calibration after sensitivity tests in 2D with the results shown in Section 3. This Manning value gave us a bottom drag coefficient of 0.0025 which represents a sand/gravel habitat [41]. Vertical dimension was discretised into 19 equidistant vertical layers.
Tidal farms with 100 and 300 turbines were tested to evaluate their influence on hydrodynamics. The tidal turbines were represented as a sub-grid object via an actuator disc model [34], implemented using Matlab scripts to automate the process [42]. The turbines were assumed to be aligned with the flow direction and were modelled using a fixed drag coefficient C D and thus C D = C T (thrust coefficient). The drag force is defined as: where ρ is the density of water (equal to 1025 kg/m 3 ), α a correction factor (set to 1 in this study), C D is the drag coefficient, A e is the effective area of turbine exposed to current and U 0 is the upstream current speed. A generic turbine design was used with a diameter of 20 m, with a constant C D setting of 0.85 and the hub height was set to 17 m [42] to represent a first generation of tidal stream energy converters. Each turbine was arranged in one grid cell of the numerical model for a better representation of the interactions between the turbines. The array layout was arranged with 2.5D (with D the turbine diameters) lateral spacing, recommended by the European Marine Energy Centre (EMEC) standard [43] and according to [44] inter-device spacing could be less than 3D. For the longitudinal spacing 6D was chosen as a compromise between 10D (recommended by EMEC) who allows better recovery of the flow and the characteristics (bathymetry and area of strong current) of the area.
To enable velocity profile comparisons between the 2D model and the 3D model, the power law defined as: was used to transform the 2D results. This equation relates the velocity profile U z with bed roughness (β = 0.4 to consider tidally energetic environment [45]), a power law coefficient α, h the water depth, U the depth average velocity and z being the height above the seabed. It should be pointed out that this equation does not take into account wind and wave-driven surface effects which could alter the velocity profile [46,47].

Validation
Calibration, validation techniques used for the models and results for the baseline are outlined in this section.

Free Surface and Velocities
The 2D and 3D models were calibrated against five ADCP measurements (C1, CW2, CW3, CW4 and CWTb1 /cf. Table 2) from the first AUSTEn campaign in March 2018. Calibration was performed by adjusting the seabed roughness with the Manning number (M) in 2D and roughness height parameter in 3D, with seven tests to obtain the closest match of surface elevations, depth-averaged velocities (DAV) and velocity profiles between ADCP measurements and numerical model's predictions, as shown in Figure 4. The phase shift between predicted and measured timeseries was removed following the EMEC guidelines [43]. Taylor diagrams were used to rank model skill depending on surface elevation and current speed accuracy against the field measurements at the location of the five ADCPs. These Taylor diagrams allow visual comparison of three statistical metrics (Standard deviation, Correlation coefficient and Centered Root Mean Square) and provide a summary of the performance of each simulation. For both surface elevations and DAV, the runs R4 (M = 32) and R9 (M = 30) in 2D showed good agreement with the five ADCP data (Figure 4). Three dimensional models increase computational cost so only R4 and R9 and R6 (reference for lower quality) were compared, with R4 the closest to the measurements (Figure 5a). Further comparisons between R4 and R9 were performed, examining current speed at the hub height, velocity profiles (Ebb and Flood) and harmonic analysis of the principal constituent. Ebb direction was defined for this work as the tidal current going from Bass Strait to Tasman Sea (eastwards) and flood from Tasman Sea to Bass Strait (westwards). R4 and R9 showed similar results for the velocities at the hub height layer (Figure 5b,c).  For velocity profiles, the data were decomposed in ebb and flood, with results showing good agreement with the five ADCP measurements ( Figure 6). The models correlated well with all stations: both runs are pretty close to the measurements. The highest correlation between in situ data and numerical model was found for simulation run R4 at station CW2 deployed 2.2 km southwest of Clarke Island, CW4 deployed 7.5 km north of Swan Island and CWTb1 in the deep pocket at the eastwards entrance of the channel. Those results confirmed the ability of the model to represent accurately the vertical profile distribution of the velocity in Banks Strait.

Harmonic Analysis
Harmonic constituents were derived from the timeseries of the 35 days simulations, the dominant harmonic constituent M2 was determined using the Matlab program UTide [48] from surface elevation and current velocities for both the MIKE models and the ADCP stations. The M2 constituents found were validated for both amplitude and phase for runs R4 and R9, with similar differences (Table 3). Among all comparisons, R4 had the highest correlation for both 2D and 3D depth average velocities and velocity profiles, therefore it was chosen to represent the conditions in Banks Strait. The run R4 was then run for one year in 2D and validated against measurements for three different periods with different ADCP: Autumn (22/03/2018-09/07/2018) with C1, CW2 and CWTb1; Winter (13/07/2018-06/09/2018) with CW1, CW2 bis and CW4 bis and Summer (06/12/2018-15/02/2019) with C1 bis. These large validation datasets give confidence in the ability of the model to reproduce tidal currents in Banks Strait. All periods showed good correlation for DAV, with slightly lower quality results for Winter period with a correlation coefficient of 0.86 (Figure 7). To evaluate the influence of waves on simulation accuracy in Winter, a small subset of the data was analysed: a period of smaller significant wave height in Winter (31/07/2018 to 06/08/2018) was considered from the measurements of waves for CW4 bis, leading to a better correlation coefficient: R = 0.9. The 2D models were also validated using data available for tidal range from the Australian Baseline Sea Level Monitoring Project [49] for the following stations: Portland

The Hydrodynamics in Banks Strait: Baseline Results
The models were applied to develop a baseline case with no tidal turbine and hence no tidal energy extraction to evaluate the current circulation in the strait and determine suitable locations for the placement of tidal turbine farms. Similar to previous studies where the strongest currents were found in the middle of the channel south of Clarke Island [18], the 3D model results reveal mean current speeds of 1.2 m/s and maximum speeds of 2.5 m/s located at (596878/−4498311 UTM 55) for the period of 22/03/2018 to 26/04/2018 (dataset of 35 days). The average maximum for the mid water column is illustrated in Figure 8 for two spring-neap cycles, with the time series of current speed at the surface for one spring-neap cycle at the location of the maximum of velocity. The modelled currents displayed 2 eddies near the area of strong currents: with one at the vicinity of the sand dunes, the presence of this feature could be strongly linked to this eddy (Figure 9).

Influence of the Tidal Farms on the Hydrodynamics
To evaluate the influence of tidal farms on the hydrodynamics (current speed, bed shear stress), simulations were performed comparing the spatial influence of two tidal farms between 2D and 3D, the influence near the seabed and the hydrodynamic influence on sand dunes near Clarke Island to evaluate potential disturbance in the ecologic system. The tidal farms are located south of Clarke Island where the strongest currents were found and where the bathymetry did not exceed 50 m (operational water depth).

Comparison of Model Spatial Extent 2D/3D
Simulations were performed in 2D and 3D in Banks Strait with and without tidal farms of 100 turbines (TF100) and 300 turbines (TF300). The spatial extent of the difference in mean current speed over 35 days between the baseline (no tidal energy extraction) and the scenario TF100 was limited to less than 5km from the tidal farm location for both 2D and 3D DAV (Figure 10a,c). For the TF300 scenario, changes in mean current speed were found to exceed to more than 7 km (Figure 10b,d).
The change in mean current speed was greatest at the location of the farms due to the removal of energy by tidal turbines. For the TF300 case in 2D and 3D, a slight increase of the current speed (0.04 m/s) is noticed between the farm and Clarke Island where sand dunes are present, and also south of the tidal farm, resulting from the barrier effect of the turbines on the flow. In both models (2D/3D), the scenario TF100 (Figure 10a,c) had less influence on the reduction of current speed than the scenario TF300: a maximum of 0.05 m/s against 0.12 m/s for TF300 (Figure 10b,d). The values are quite similar between 2D and 3D, the major difference is the spatial extent where 2D displays a larger influence which could be considered as a conservative measure. For more insight about the difference between the two arrays, the mean current speed over 35 days was extracted at the mid water column in the 3D scenario ( Figure 11). The influence of tidal energy extraction is higher at the location of the turbines for the TF300 scenario: above 0.2 m/s reduction against around 0.1 m/s for TF100. These results are slightly lower than [5] who found a mean reduction of 0.3 m/s with a farm of 290 turbines in the Alderney Race.  A profile AA' (location on Figure 8) parallel to the direction of the currents was extracted to examine the spatial influence of the tidal farms on depth averaged velocities for the 35 days dataset in 2D and 3D (Figure 12). The drop for the baseline case at 2 km is caused by the presence of a canyon. The 2D model baseline results for current speed were overall 5% higher than the 3D DAV results, corresponding to an increase of approximately of 0.05 m/s along the profile. A significant drop in current speed was noticed for both farms at their location, with a greater reduction for the scenario TF300 as expected due to higher levels of tidal energy extraction, equivalent to a reduction of 7% in the 3D model ( Figure 12). For both dimensions, the differences were lower at approximately 13km from the tidal farm reducing to less than 1% up to 0% at 18 km towards the Tasman Sea for the two tidal farms ( Table 4). The influence of the 100 turbine array was minimal in the near field with less than 1% of difference in mean current speed. For smaller scale tidal farm development, this is a positive outcome.  For 2D models, the difference in the magnitude of velocity between the baseline and scenarios with tidal farms was greater than in 3D DAV at the vicinity of the arrays, however in the far field the average difference was the same. The 2D model results match the 3D DAV results for the TF100 scenario with a maximum difference of 0.4% and are similar in the far field for the TF300 scenario ( Table 4). The 2D models overestimated slightly the current speed compared to the 3D DAV results, with the difference between baseline and scenario with turbines highest at the location of the TF300. This should be considered if a preliminary resource assessment is performed in 2D. At the location of the tidal turbines, the velocities difference between 2D/3D reached an average of 0.1 m/s corresponding to 12%, which is lower than the results of [13] who found a 20% difference in an idealised channel similar to Pentland Firth. However, the difference between 2D and 3D DAV was found to be negligible in the far field of the tidal farms.
Small discrepancies were observed between 2D and 3D models in the near and far field. Given the high computational cost required by 3D models, the authors recommend that preliminary assessments be performed in 2D. These results are in accordance to previous studies in an idealised channel for the difference of dimension with a focus on the bed level change. The 2D models are a good compromise between model accuracy and time management, especially for long temporal scale assessments for a first approximation.

Influence Near the Seabed
To enable comparisons of the 2D model to the 3D model at turbine hub heights, the power law defined in Equation (1) Section 2.3 was used. Previous studies have demonstrated that the 1/7th power law (α = 7) is generally the most suitable to capture vertical processes in coastal oceanographic research; moreover, [45] showed that α = 7 correlated well with tidal energy sites conditions in the Irish Sea. However, guidelines [43,50,51] and field data observation [52] recommend a 1/10th power law in tidal energetic sites. To examine the influence of power law factors on simulation results, two α coefficients (7 and 10) were analysed with the velocity profile for four points (inside the tidal farm then at 0.9, 4.5 and 18 km from the tidal farm towards Bass Strait) for the TF300 model as shown in Figure 13. Results indicated that the 1/10th power law correlated the best with 3D results, as the difference was 14% instead of 25% for the α = 7 model. There was a significant difference of velocity at the location of the tidal farm between 2D and 3D. The 3D model represented well the decrease of the flow in the water column due to the presence of the turbines. This velocity deficit at the hub height cannot be reproduced in the 2D model with the power law and thus lead to a misrepresentation of the current speed in 2D at the location of the tidal farm confirming the findings of [13]. It should also be noted that some tidal energy sites have more complex velocity profiles which do not fit the power law equation [53]. Future research using coupled hydrostatic and non-hydrostatic models will be able to capture flow properties at scale ranging from ocean to tidal turbines blades enabling the accurate modelling of all hydrodynamics flow properties including wake environment and velocity profiles. Friction between the near seabed currents and the seafloor creates a layer of shear and increased turbulence above the seabed called the bottom boundary layer. This is the layer where the flow is the most influenced by the seafloor and alterations in this layer could lead to significant impact in sediment transport patterns. The profile AA' was extracted at the bottom layer for the mean current speed (Figure 14a,b) with a coefficient α equal to 10 for the power law. An increase of current speed was noticed at the location of the tidal farm only for the 3D models: with magnitudes of 1 cm/s for TF100 and 3 cm/s for TF300 (2% and 5%, respectively) found. Away from the tidal farm, the reduction in current speeds are quite similar for 2D and 3D: for TF300 6% and 5.6% (3D, 2D, respectively) and 2% for TF100 at 2.2 km from the tidal farm. The differences fade away in the far field reaching less than 1% at 13.4 km for both dimensions. The 3D models appear to model the acceleration caused by the constriction of the flow beneath the turbines, which the 2D models cannot capture. In the 2D model, the bed shear stress is derived from the DAV velocities, whereas the bed shear stress is determined from the bottom cell velocity in the 3D model; therefore, the profile for the bed shear stress followed the same fluctuation as in the bottom layer in 3D with an increase at the location of tidal farm (Figures 14c,d and 15). Around the tidal farms the bed shear stress was reduced by 7% (0.22 Pa) for TF100 arrays in 3D, 6% for 2D, with the differences not significant in the far field. The reduction of bed shear stress for the TF300 scenario was greater reaching 20% (equivalent to 0.58 Pa) around the tidal farm in 3D and 15% (0.51 Pa) in 2D. The 2D results were similar to the 3D results for the scenario TF100 and TF300 in the far field. These results vary compared to the influence of turbines in Ramsey Sound as shown by [7] who found an average reduction of 2.3 Pa in the Ramsey sound with less devices. The peak reduction found here is 0.9 Pa, 25 times lower than [4] with 25 Pa in the Pentland Firth. These discrepancies can be explained by the faster flows occurring in these two sites and therefore a greater alteration in hydrodynamics. Influence on bed shear stress was observed with the TF100 scenario and was greater with the 300 turbines arrays with a spatial influence of 7 km instead of 4.5 km for TF100 (Figure 15), results close to the spatial extent of 12 km found by [7]. Reduction in the bed shear stress around the tidal farm suggests that accretion of sediments may appear and change the dynamics in the channel as it can be seen on maps on Figure 15. The increase of bed shear stress at the location of the farm and at the top and bottom suggest that additional scour and more erosion could occur at the location of the sand dunes between the tidal farm and Clarke Island. A closer look was given to the region of the sand dunes near Clarke Island given the proximity to a potential tidal farm. Alteration in current speed and bed shear stress could lead to a disturbance of the benthic ecology system. Bed shear stress results in 3D and 2D were compared between the case with no tidal extraction and the two scenarios TF100 and TF300 ( Figure 16). The change in the magnitude for the 100 turbines arrays was less than 0.12 Pa and less than 0.26 Pa for the 300 turbines case for both 2D/3D. The difference between 2D and 3D results were not significant: a maximum of 4.5% difference with an average difference around 0.06 Pa for all the scenarios, but the models did not display the same fluctuations. The slight change in bed shear stress at this location needs to be better understood, with more studies on the dynamics of these sand dunes. To pursue the EIA, the focus needs to be on the dynamics of sediment, especially at the location of the sand dunes which are of economic importance. They may provide sediment for the close beaches of Clarke Island, protect them from erosion during storms and may contribute to the fishing industry being a source of nutrients for many fish species [54,55]. These results reinforces previous findings [7,8] that, for preliminary assessment of oceanographic changes of a potential tidal energy sites, 2D models are suitable for examining far field parameters (current speed and bed shear stress). These 2D studies can give a first insight of the potential influence of arrays on a longer time frame given the reduced computational costs when compared to 3D models. Thus, 2D models are assets for the dynamics of sediment transport and are an effective tool to understand the influence of the removal of energy. However further work is needed to understand the capability of 2D models with more complex and challenging scenarios such as sand dunes migration. 3D models remain essential to have more insight at the location of tidal farm; however, they become less powerful in the near and far field as it gives similar results as the 2D model in terms of hydrodynamic environmental influence.

Conclusions
The hydrodynamics of Banks Strait was examined using two high resolution models in 2D and 3D, with resolution varying from 3 km to 20 m in the area of interest to simulate the presence of tidal farms. The developed models were successfully validated against a large dataset of ADCP field measurements, with results indicating good agreement for surface elevation, depth average velocities and velocity profiles. Validation for a full year was then performed for the 2D model with close agreement to the field data found. Those results gave confidence in the ability of these models to reproduce the hydrodynamics in Banks Strait.
The implementation of tidal turbines in the numerical models showed the footprint effect of tidal turbines on the hydrodynamics (current speed, bed shear stress). Scenario with 100 turbines showed minimal influence on the velocities in the near field, with changes in velocity only found less than 3km from the farm. For the alteration of bed shear stress, this scenario did not display any significant difference compared to the baseline, which is a relatively good outcome for the deployment of a small tidal farm. The spatial influence of the 300 turbines farm was greater reaching 7 km from the farm for the difference in the current speed. A small number of devices will not change the system in a significant way compared to the natural variability; but, as these studies have shown, large scale arrays have the potential to influence the natural process. These alterations can lead to greater changes on the dynamics of sediments.
In the far field, the difference between the two dimensions 2D/3D was negligible, confirming that 2D models are suitable for investigating the influence of tidal arrays for a preliminary environmental assessment. At the location of the tidal farm 2D models were found to have slightly higher values for the influence of tidal energy extraction and thus precaution should be taken for resource assessment. Regarding local features here in Banks Strait, more work is needed to fully understand the effect of the tidal arrays on the sand dunes compared to the natural variability.
The findings of this paper are a first step in the understanding of the influence of tidal energy extraction in Banks Strait and gave a solid basis to reduce the uncertainty in the hydrodynamics for future investigation of sediment transport in the region. Future work will focus on the creation and calibration of a sediment transport model with in situ data (multi-beam, echosounder, bottom grabs, sediment traps and sub-bottom data) collected by the AUSTEn project. The migration of sand dunes will also be of a high interest to better understand the influence of anthropogenic activities on the physical environment.