Anthropogenic Mixing of Seasonally Stratified Shelf Seas by Offshore Wind Farm Infrastructure

The offshore wind energy sector has rapidly expanded over the past two decades, providing a renewable energy solution for coastal nations. Sector development has been led in Europe, but is growing globally. Most developments to date have been in well-mixed, i.e. unstratified, shallow-waters near to shore. Sector growth is, for the first time, pushing developments to deep water, into a brand new environment: seasonally stratified shelf seas. Seasonally stratified shelf seas, where water density varies with depth, have a disproportionately key role in primary production, marine ecosystem and biochemically cycles. Infrastructure will directly mix stratified shelf seas. The magnitude of this mixing, additional to natural background processes, has yet to be fully quantified. If large enough it may erode shelf sea stratification. Therefore, offshore wind growth may destabilize and fundamentally change shelf sea systems. However, enhanced mixing may also positively impact some marine ecosystems. This paper sets the scene for sector development into this new environment, reviews the potential physical and environmental benefits and impacts of large scale industrialization of seasonally stratified shelf seas and identifies areas where research is required to best utilise, manage and mitigate environmental change.


INTRODUCTION
Renewable energy solutions, including offshore wind, are prerequisite for clean growth and thus the reduction of greenhouse gas emissions needed to mitigate against climate change. Offshore wind energy in shelf seas has seen a rapid increase over the past decade (Xu et al., 2020;Díaz and Soares, 2020), motivated by: high-quality and reliable energy (wind) resources (Esteban et al., 2011); space availability and site accessibility for installation of large, efficient, turbine systems (Sun et al., 2012); rapidly maturing, reliable and energy-efficient technologies (Jansen et al., 2020); and reduced visual impact on populated areas (Wen et al., 2018). Government programmes have helped drive development of renewable offshore wind energy from offshore wind farm arrays, of tens increasing to hundreds, of offshore wind turbines (OWT) supported by various fixed foundation designs with new floating foundations being designed to access deeper water  . Offshore wind and seasonally stratified shelf seas. (a) Shelf sea development; the limited extent of well-mixed waters, defined using the "h/u 3 M 2 " criterion, is highlighted against projected global offshore wind energy growth (Global Wind Energy Council, 2020). (b) NEMO model of the Northwest Europe summer potential energy anomaly, φ, a measure of the amount of stratification (Guihou et al., 2018). (c) Fishing hotspots in the Celtic sea caused by topographically-enhanced mixing of stratified waters . In (a) and (b) seas are partitioned into regions prone to seasonal stratification and those remaining well mixed, based on: (a) h/u 3 M 2 = 220 s 3 m −2 (Simpson and Sharples, 2012); and (b) φ = 20 J m −3 (Gowen et al., 1995). In (b) and (c) offshore wind farms are separated into active sites, sites still in construction or development and identified zones for future development (4C Offshore, 2021).  Figure 2. Existing and emerging offshore wind engineering solutions (Díaz and Soares, 2020;Jiang, 2021b), including fixed monopile and jacket foundations and floating semi-submersible and spar-buoy foundations, in comparison to shelf sea regimes from coastline to open ocean (Simpson and Sharples, 2012).
sites. Northwest Europe has led sector development, with the UK leading in GW operational capacity to date (Global Wind Energy Council, 2020). The sector has grown rapidly, with technological advances reducing the Levelised Cost of Electricity to a point where price is competitive with alternative energy solutions (Shen et al., 2020). Thus, to meet demand, global development of offshore wind energy in shelf seas is predicted to grow from 35 GW operational in 2020 to 243 GW operational by 2030 ( Figure 1a).
With over 80% of the global population living within 100 km of the ocean, shelf seas have significant economic and social value, including fishing, shipping, carbon storage ("blue" carbon) and recreation. Despite comprising 8% of the total area of the global ocean (Figure 1a), shelf seas support 15 -30% of global ocean biological production (Wollast, 1998). This biological production ultimately supports > 90% of the world's fish landings (Pauly et al., 2002) and plays a disproportionately important role in the absorption of CO 2 from the atmosphere (Roobaert et al., 2019). Thus, high biological productivity means shelf seas are key components of global biogeochemical cycles, supporting societally important bioresources and also the biological uptake and storage of carbon in the marine environment. However, interplay of social and economic drivers already places significant stress on shelf seas (Kröger et al., 2018). Further industrialization of shelf seas will enhance these stresses, with the potential for significant long term environmental impact. Shelf sea dynamics directly control primary production: the growth of microscopic marine plankton. However, from OWT scale to coastal scale, the impact of offshore wind development on shelf seas has yet to be fully considered. Therefore, future offshore wind energy development must be grounded in advanced understanding of impact on shelf sea dynamics. This is critical to enable balance of key global Societal Goals, i.e. to "ensure access to affordable, reliable, sustainable and modern energy" and to "conserve and sustainably use the oceans, seas and marine resources" (United Nations, 2015).
DRAFT R. M. Dorrell et al.

Infrastructure Mixing of Shelf Seas
To date, most offshore wind farms have been installed in the near-shore shallow water regions, up to 50 m depth, of shelf-seas ( Figure 1b). Near-shore shallow-water installations have been preferred due to the cost reduction from ease of access for installation, grid connection and operation and maintenance (Jacobsen et al., 2019). With sector plans for an additional 208 GW of operational capacity in the next decade, and targets of 1.4 TW total by 2050 (Offshore Renewable Energy Action Coalition, 2020), near-shore and shallow-water sites are rapidly becoming limited. The scale of expansion of offshore wind energy means the sector is now expanding into deeper water sites further from shore (Soares-Ramos et al., 2020). The transition from near-shore and shallow-water environments to deeper water further from shore marks a fundamental change in the marine environment. Shallow waters are typically well-mixed; however deeper waters may be subject to seasonal stratification, where density varies vertically with depth (Figure 1a and b). Stratified waters are a vital part of shelf seas, controlling primary production and biogeochemical cycling (Simpson and Sharples, 2012). Expansion into this new environment means that offshore wind farms will increasingly come into conflict with its environmental functioning, controlled by natural mixing of water column stratification (Figure 1c).
Addressing engineering challenges, both fixed and floating foundations are being developed to enable expansion into deeper waters. Fixed foundations, which span the entire water depth, include monopiles, gravity bases and jacket constructions (see, e.g, Figure 2 and Esteban et al., 2019;Díaz and Soares, 2020;Jiang, 2021b). However, floating foundations are crucial to deep water, > 50 metres, deployment. Learning from the petroleum industry (Schneider and Senders, 2010), designs include tension-leg platforms (Uzunoglu and Soares, 2020), spar-pendulum (Cottura et al., 2021) and spar-buoy platforms (Jacobsen and Godvik, 2021), and semi-submersible platforms (Castro-Santos et al., 2020). Using the submerged structural buoyancy and mooring forces to balance atmospheric thrust and wave loads, floating foundations typically have large draft, e.g. spar platforms, or large cross sectional area, e.g. semi-submersible platforms (Butterfield et al., 2007). Thus, with sector development requiring larger of turbines that need bigger rotors, which are subject to greater atmospheric loads, the draft and diameter of fixed and floating foundations will need to increase. However, the > 20 m draft of current floating foundations is already large enough to penetrate the thermocline and directly mix seasonally stratified shelf seas ( Figure 2). For the first time, large scale industrialisation of seasonally stratified marine environments is planned. Over two decades of research has already focused on the direct impacts of offshore wind farm development on well-mixed shallow water marine ecosystems, from: benthic habitats (Dannheim et al., 2020), fisheries (Gray et al., 2005) to seabirds (Exo et al., 2003). Whilst this research is translatable with sector growth, the seasonally stratified regime offers a fundamentally new challenge: the introduction of infrastructure will lead to enhanced 'anthropogenic' mixing of stratified waters. Enhanced mixing may lead to profound impacts on shelf sea dynamics and thus marine ecosystem functioning. The aim of this paper is to investigate the scope of these potential impacts. Section 2 reviews the ecosystem and physical functioning of stratified shelf seas. Section 3 then describes our current understanding of the impact of offshore infrastructure on unstratified and stratified waters. Section 4 discusses current research challenges, the potential impact of offshore wind on stratified shelf seas and the sector requirements needed to ensure acceleration of renewable energy and its sustainable development. It is concluded that offshore wind farm infrastructure may have significant, and long lasting, effects on fragile shelf sea ecosystems. Criteria for Environmental Impact Assessments, must therefore be revised and updated to enable the sustainable growth, and acceleration, of renewable offshore wind energy development.

OCEANOGRAPHY OF STRATIFICATION IN SHELF SEAS
Shelf seas lie on the continental shelf between the coast and the continental slope, where at 200 m water depth the sea floor slopes down to the deep ocean. Despite only accounting for 0.5% of ocean volume shelf seas play a key role in the Earth system, dissipating > 70% of the tidal energy (Egbert and Ray, 2000) and are disproportionately important in supporting ocean biological production (Wollast, 1998), fish landings (Pauly et al., 2002) and the absorption of CO 2 from the atmosphere (Roobaert et al., 2019). Biological production is underpinned by the growth of microscopic marine phytoplankton, which is tightly controlled by the timing and strength of seasonal stratification. The seasonally stratified zones in shelf seas act as an important net sink of carbon (Thomas et al., 2004). This makes the physical and biogeochemical processes described here a key dynamic component of the global carbon cycle (Bauer et al., 2013) linking the atmospheric, terrestrial, and oceanic carbon pools. Figure 3. Seasonal time-series of mooring observations from 2014 in the Celtic Sea (49 • 24 N 8 • 36 W), collected as part of the UK Natural Environment Research Council (NERC) funded CaNDyFloSS project (NE/K00168X/1) (Scannell, 2020). a) Measured thermocline current speed over 20-50 m mean depth; average speed (blue line), and range (shaded region) over 2 tidal periods. b) Surface heat exchange. c) Vertical temperature structure. d) Surface and bottom layer temperatures.

Distribution and Seasonal Cycle of Stratification
Away from coastal regions of freshwater influence, the water column structure in temperate shelf seas undergoes a seasonal cycle in response to changes in heat exchange at the surface. In spring and summer, some areas of the temperate shelf seas become thermally stratified whilst neighbouring areas remain well mixed. The first order control on the water column structure has been established as a balance between the stratifying influence of surface heating and turbulent mixing due to the tides (Simpson and Hunter, 1974). In regions of shallow water and strong tidal currents (of order metres per second), the rate of buoyancy input due to surface heating is insufficient for the establishment of persistent stratification, and in consequence the water column remains homogeneous. However, in regions of deeper water and/or weaker tidal currents (and associated lower levels of turbulence), surface heating dominates and seasonal stratification develops. Away from coastal regions, shelf seas generally exhibit tidal currents of order of tens of cm per second, meaning that stratification are typically found when water depths are greater than 80 m. In stratified regions, a warmer surface water layer 5-40 m thick overlies a deeper cooler water layer. The two layers are separated by a region of strong vertical temperature gradient, the thermocline, which forms a barrier to vertical exchange of heat, salt, nutrients and momentum.
Tidal mixing fronts separate regions of seasonal stratification from well mixed regions. Simpson and Hunter (1974) use an energetics argument to derive a single parameter to predict the positions of these fronts. By considering only vertical exchange processes and assuming the surface input of heat was the only stratifying influence, and that tidal currents are the only source of energy driving mixing, they showed that the first order determinate for the position of shelf sea fronts is given by the ratio h/u 3 M 2 , where u 3 M 2 is the principle lunar M2 tidal current amplitude, and h is water depth (Figure 1a). In Figure (1a) h/u 3 M 2 is calculated from bathymetry and M 2 tide data taken from TPX09 global tidal atlas (Egbert and Erofeeva, 2002), which applies a generalised inverse method, assimilating satellite altimeter data, into a global barotropic tidal model; here model resolution limits precise location of stratified fronts from DRAFT global data. In terms of area, regions of seasonal stratification dominate the continental shelf seas ( Figure  1a). Whilst the critical value for the ratio characterising the position of tidal mixing fronts was initially estimated for the Irish Sea (Simpson and Hunter, 1974;Simpson and Bowers, 1981), consistent values have subsequently been estimated for a range of shelf seas globally [e.g. the Gulf of Maine and Bay of Fundy, Garrett et al. (1978); Loder and Greenberg (1986); the Yellow Sea, Lie (1989); the Patagonian Shelf, Glorioso and Flather (1995); the northwest European Shelf Seas, Pingree and Griffiths (1978) and the Bering Sea, Schumacher et al. (1979)]. The robustness of the critical value highlights the key role of the tides in determining the position of shelf sea fronts and provided the first quantitative link between the dissipation of tidal energy and ocean mixing.
The strength of stratification may be quantified in terms of the potential energy anomaly, φ (J m −3 ), which describes the energy required to fully mix a stratified water column (Simpson and Bowers, 1981) where h is the water depth and ρ the water density, ρ denotes the density calculated using the mean water temperature and salinity (Holt and Proctor, 2008). Geographical variation of φ across the NW European Shelf Seas is plotted on the map in Figure 1b.
A typical seasonal stratification cycle shows a time-series of warming and cooling, varying with water depth (Figure 3). Surface mixed layer temperatures warm from April into the summer in response to a net positive buoyancy input due to surface heating. The observed surface temperatures of 15 − 20°C are typical for temperate shelf seas in midsummer, while the deep water remains close to its winter temperature of 10°C. Over this period the strength of stratification grows with a surface to bed temperature difference exceeding 10°C by August, which gives rise to a density difference of 0.3 kg m −3 . The stratification weakens into autumn as surface cooling leads to convection and storms drive turbulent mixing, such that the water column becomes well mixed during the winter.
During the stratified period the deeper water is isolated from the surface layer by a thermocline. A slow warming of the deep-water results from mixing down of heat from the sea surface. The rate of warming, set by thermocline mixing, varies geographically and is important as it determines the timing of the autumnal breakdown of stratification Rippeth (2005), and transport of nutrient rich deep water up to the surface layer.

Ecosystem Response to Stratification
Primary production of organic matter by phytoplankton forms, directly or indirectly, the primary food source for almost all marine organisms. Phytoplankton growth requires CO 2 , sunlight and nutrients, the availability of which are determined by water column structure, with profound implications for the biological functioning of the shelf seas.
In turbulent unstratified regions, primary production occurs mainly during summer months when sunlight is strong. However, plankton are continuously mixed from sea surface to bed by turbulence, and spend much of their time below a depth where light intensity is sufficient for growth. In contrast, stratified waters provide ideal conditions for phytoplankton growth in spring. As stratification forms phytoplankton become trapped in the well lit surface layer, with the thermocline acting as a barrier to mixing. Phytoplankton retained in the surface layer enjoy the abundance of light, and exhibit rapid growth forming the annual "spring bloom", a biological abundance visible from space that forms the year's first supply of significant new organic fuel. As the phytoplankton grow, they fix inorganic carbon in the surface water into organic carbon, which causes the sea surface to replenish its dissolved carbon concentration by absorbing CO 2 from the atmosphere. The timing of the spring bloom is so significant that zooplankton and fish larvae have evolved to use it as a food source (Platt et al., 2003), with further implications higher up the food-web, e.g. for shrimp survival (Ouellet et al., 2011) and seabird breeding success (Frederiksen et al., 2006).
During the spring bloom, the availability of nutrients in the surface layer becomes exhausted, and further production is limited by nutrient supply. Despite this limitation on plankton growth, a persistent and significant level of primary production is sustained at depth, throughout the period of seasonal stratification. This sub-surface phytoplankton layer located in the stratified thermocline water is a ubiquitous feature and is known as the "subsurface chlorophyll maximum" (SCM) (Pingree et al., 1982). In shelf seas, the SCM occupies a 10-30 m thick layer, the depth of which varies from 10-40 m across the hundreds of kms DRAFT Figure 4. Observational data from the Celtic Sea collected in summer 2008, supported by NERC's Oceans2025 Programme . Section of temperature (line contours) and chlorophyll concentration (colours) measured using a Scanfish CTD towed along the path show on the map. High chlorophyll concentration in the SCM (subsurface chlorophyll maximum) indicate phytoplankton production extending hundreds of kms across the shelf. Enhancements in concentration over rough topography, such as Jones Bank, are a result of elevated turbulent mixing, driving nutrient fluxes which correspond directly to hotspots of marine biodiversity and thus fisheries ( Figure 1c). that it extends over the shelf (Figure 4). The SCM plays a vital role in supporting the pelagic food web during summer. Estimates based on observations of primary production rates within the SCM suggest that subsurface carbon fixation accounts for up to 50% of annual primary production in the seasonally stratified North Sea (Richardson et al., 2000;Weston et al., 2005). An extrapolation using microstructure-based nitrate flux estimates also gives the same approximate figure . The persistence of production in the SCM is dependent on a vertical flux of nutrients from the deep nutrient-rich water below (Sharples and Tett, 1994). In consequence, the processes responsible for mixing across the thermocline, discussed in section 2.5, are key to delivering the limiting nutrients to the euphotic zone and sustaining the SCM (Sharples et al., 2001b(Sharples et al., ,a, 2007Williams et al., 2013a). Episodic mixing associated with storm initiated inertial oscillations (Burchard and Rippeth, 2009;Lincoln et al., 2016) have been shown to drive significantly enhanced nutrient fluxes (Williams et al., 2013b). Local enhancement of primary productivity is evident in regions of steep topography, where tidally induced internal waves elevate mixing and nutrient fluxes. Chlorophyll concentration in the SCM is greatly elevated at the shelf break (Sharples et al., 2007), and mid-shelf sand banks . The dependence of the marine food web and fisheries on upward flux of nutrients is evident in the distribution of seasonal fishing hot-spots ( Figure 1c).
The organic products of the spring and summer primary production sink into the deeper waters, where bacteria remineralise the organic material (nutrients and carbon) back to the inorganic components. Remineralisation removes oxygen from the deep water. In addition, the barrier role of the thermocline limits the replenishment of that oxygen from the atmosphere (Mahaffey et al., 2020). Together both processes determine the dissolved oxygen concentrations available to benthic and pelagic organisms. Thus, high shelf sea biological productivity means shelf seas are key components of global biogeochemical cycles, Marginal stability and energy dissipation in seasonally stratified shelf seas. a) The equivalence between buoyancy frequency N 2 (blue) and vertical current shear S 2 (red), temporally averaged (over two tidal cycles). b) The coincident temporally averaged profile of TKE dissipation rate ε. Measurements from direct observations in the Western Irish Sea in June 2002 (Rippeth, 2005).
supporting societally important bioresources, and also the biological uptake and storage of carbon in the marine environment.

Shelf Sea Mixing Processes
Currents in shelf seas provide energy for stirring the water column and are generally dominated by tidal motions, with episodic contributions by the wind in the upper water column. An example of the current variability from the Celtic Sea, a typical shelf sea location, is presented in Figure 3a, where velocities vary from 0.1-0.7 m s −1 . For this location, the semi-diurnal tide lunar M 2 produces two high and low tides a day with 4 peaks in current speed. The interaction with the principle solar tidal component, S 2 , produces the 14 day spring-neap cycle. Wind driven currents are also observed in the top 50 m, and take the form of inertial oscillations, which have a latitude dependent period, which is 14.9 hours at the mooring location.
Friction at the seabed and the sea surface generate vertical current shear in the flow, and turbulent eddies which "cascade" to ever smaller scales until their energy is dissipated either to heat, or to potential energy via mixing.
In the absence of convection, a three-way local balance is assumed (assuming both at least quasistationarity and that transport terms can be ignored): where P is the (total) production of turbulent kinetic energy (TKE), B is the buoyancy production (mixing), and ε is the TKE dissipation rate (heat).
The efficiency of mixing by turbulence can be quantified by the flux Richardson number Rf = B/P and is widely assumed to have a value Rf 1 in a stratified fluid. Since ε is a commonly measured turbulence metric it is often used to infer the rate of mixing using the closely related flux dissipation coefficient, defined in terms of the buoyancy production as Γ = B/ε. A value of Γ ≈ 0.2 (i.e. Rf ≈ 1/6) is DRAFT routinely applied, and has been verified for the shelf sea thermocline by a number of different observational approaches (Oakey and Greenan, 2004;Inall et al., 2000;Palmer et al., 2008;Bluteau et al., 2013), though it has been found to vary in other regimes (Monismith et al., 2018).
A consequence of shear production at the seabed by barotropic tidal currents is that measured rates of turbulence are extremely low in the seasonal thermocline, orders of magnitude lower than at the boundaries ( Figure 5b). Mean dissipation rates range from ε = 10 −7 − 10 −5 W m −3 (Rippeth, 2005), 2-3 orders of magnitude smaller than rates commonly found in shallow well mixed waters (Simpson et al., 1996). Empirical estimates of the bulk mixing efficiency of the barotropic tide in stratified waters, are very low, Rf BT ≈ 0.0037 (Simpson et al., 1978), as most turbulence is produced in the well mixed bottom layer, so that no mixing is possible. In addition, strong density gradients in the thermocline inhibit vertical mixing, and as such rates of vertical mixing observed in shelf seas are comparable to using a hand mixer in a swimming pool.
The production of turbulence in stratified waters is inhibited by buoyancy forces arising from vertical density gradients, which are quantified by the Brunt-Väisälä, or buoyancy, frequency, N , where describes the frequency at which a displaced parcel of fluid will oscillate in a stratified system and is thus a measure of the stability of stratified waters. Conversely, the vertical current shear, S = ∂u/∂z, is a measure of the extraction of energy from the mean flow, and therefore power available to overcome buoyancy forces and generate turbulence. The generation of instabilities in stratified water is quantified using measurements of the buoyancy frequency, N , and vertical current shear, S, to calculate the gradient Richardson Number Both theory (Miles, 1961;Howard, 1961) and observations (Silvester et al., 2014;Lincoln et al., 2016) show that shear instability and internal wave breaking occurs when Ri g 0.25. Average of measurements, over two tidal cycles, show N 2 and S 2 are approximately equal implying a (close to) marginally stable water column where Ri g ≈ 1 ( Figure 5a). Therefore, any additional sources of shear can generate instabilities and mixing. This marginal stability has been widely observed across shelf seas (van Haren et al., 1999;Rippeth, 2005;Palmer et al., 2008) but it does not universally explain the observations (Palmer et al., 2013) suggesting additional sources of shear may, at times, control system dynamics.
One dimensional vertical exchange models using turbulence closure schemes (also commonly referred to as RANS models) fail to represent the measured rates of mixing through tidal and wind driven boundary processes (Simpson et al., 1996). The deficit in the predicted mid-water ε points to either an incorrect parameterisation of the small scale physics away from the boundaries or to the absence of key physical processes in the model. A recent study (Luneva et al., 2019) evaluated a range of alternative one dimensional turbulence closure schemes for the Northwest European shelf seas, packaged within the Generic Length Scale two-equation formulation (Umlauf and Burchard, 2005). Evaluating the schemes against profile data (28000 profile in total) confirmed that there was no outright winner, with all schemes under representing the thermocline properties and suggested that physical processes are still missing. Candidate mechanisms to account for the deficit in mid water mixing include internal waves generated by stratified flow over topography and wind generated inertial currents.
Internal tides propagate at the thermocline in response to tidal currents flowing over steep topography (Rippeth, 2005;Inall et al., 2021). These waves generate strongly sheared currents about the thermocline which can lead to the development of shear instability, draining energy to turbulence which then supports mixing in the quiescent mid-water region. Although the energy in the internal tide is much less than the barotropic tide, an empirical estimate of bulk mixing Rf IT ≈ 0.056 (Stigebrandt and Aure, 1989) shows it to be highly efficient compared with the barotropic tide Rf BT ≈ 0.0037 (Simpson et al., 1978), on account of the turbulence coinciding with the stratification. A consequence of the differing bulk mixing efficiencies is that whilst the barotropic tide dominates the budget for the dissipation of barotropic tidal energy in the DRAFT 9 R. M. Dorrell et al. seasonally stratified shelf seas to the west of the UK, the shelf break generated internal tide is estimated to dominate diapcynal mixing . With advances in computational performance and associated resolution increase, simulation of internal tides are now routine for European Shelf simulations (Guihou et al., 2018) implemented operationally (Graham et al., 2018;Tonani et al., 2019). Though, as hydrostatic simulations, the non-hydrostatic mixing processes must still be parameterised.

Submesoscale Mixing
The size of many existing and proposed offshore wind farms coincides with the ocean submesoscale (∼ 1 − 20 km). The imprint of submesoscale eddies can be seen in satellite images of chlorophyll concentration in the North Sea as seen in Figure 6 (note that the future development zones shown for reference will likely consist of multiple smaller windfarms). The submesoscale range is characterised by Rossby numbers, Ro ≡ u/(f L) ∼ 1, where u and L are characteristic horizontal velocity and length scales associated with submesoscales and f is the Coriolis parameter. This makes submesocales dynamically distinct from mesoscales (∼ 100 km) where Ro 1. While the Earth's rotation is important on submesoscales, it does not constrain the flow as strongly as it does for larger scales, leading to a unique set of physical processes. Submesoscale currents are generally more energetic in regions with large horizontal density contrasts which includes coastal waters.
It is not clear if, or how, offshore wind farms might interact with submesoscale currents. However, since the spatial extent of wind farms typically lie within the submesoscale range, array scale flow patterns generated by enhanced drag or mixing within offshore wind farms could be subject to the physical processes that are active on the submesoscale. Below we very briefly consider some of these processes. See Thomas et al. (2008); McWilliams (2016); Gula et al. (2021) for in-depth reviews of submesoscale dynamics.
Submesoscales are typically generated by dynamical instabilities or flow/topography interactions. The instability mechanisms include ageostrophic versions of instabilities that also exist at larger scales (e.g. baroclinic instability) and instabilities that are unique to the submesoscale (e.g. inertial/centrifugal and symmetric instability). Although the details and energy pathways differ, the net effect of submesoscale DRAFT instabilities is to convert potential energy associated with horizontal density gradients into kinetic energy. In the process, submesoscales increase the water column stratification and tend to reduce the depth of the surface mixed layer. Submesoscales can also be generated as currents move around topographic features. For example, submesoscale eddies can be generated by flow past islands (Marmorino et al., 2018), along continental slopes Gula et al., 2015), and over seamounts (Srinivasan et al., 2019). It is possible that the enhanced drag experienced by the flow through a submesoscale offshore wind farm could similarly generate submesoscale eddies. The enhanced mixing and drag within wind farms could also influence submesoscale instabilities and eddies, although these hypotheses remain untested.
Submesoscales are important partly because they interact with vertical mixing processes. Submesoscales are characterised by timescales that range from hours to ∼ 1 day. As a result, the tendency for submesoscales to increase the vertical stratification of the water column is fast enough to compete with vertical mixing driven by winds and tides. At the same time, submesoscales can induce strong vertical circulations and locally enhance the exchange between the surface mixed layer and thermocline (Mahadevan and Tandon, 2006). Gula et al. (2015) provide an in-depth review of the connections between submesoscales and ocean mixing.
Submesoscales can have a strong impact on biogeochemistry (Lévy et al., 2012;Mahadevan, 2016). For example, shoaling of the surface mixed layer depth and suppression of turbulent mixing induced by submesoscales can trigger phytoplankton blooms in otherwise light-limited conditions (Taylor and Ferrari, 2011;Mahadevan et al., 2012;Taylor, 2016). In nutrient-limited conditions, submesoscales can upwell nutrient-rich waters to the euphotic zone, enhancing primary production (Lévy et al., 2001;Mahadevan and Archer, 2000). Finally, downwelling circulation and the suppression of turbulent mixing can enhance the export of particulate organic matter from the surface mixed layer (Omand et al., 2015;Taylor et al., 2020).

Mixing Across Density Interfaces
As noted, stratified shelf seas sometimes exhibit a two-layer density structure with relatively homogeneous mixed layers at the top and bottom separated by a stratified thermocline (Figure 3). There are two paradigms for mixing, at the small scales of density interfaces; 'scouring' and 'overturning' (see Caulfield (2021) and references therein). The mixing across the density interfaces and the ultimate fate of the density interface depends strongly on which regime the turbulence is in. The mixing regimes are controlled by the relative strength of stratification and turbulence. This can be quantified by the kinetic energy associated with three-dimensional turbulent eddies and the potential energy associated with the density interface. When the stratification at the density interface is sufficiently strong, turbulent eddies do not have kinetic energy to overturn the density interface. Instead, if there is a source of turbulence in the surrounding mixed layers, turbulent eddies will 'scour' the density interface, pulling characteristic wisps of fluid into the mixed layers ( Figure 7b). On the other hand, when the kinetic energy associated with the turbulent eddies is large enough, turbulence is able to 'overturn' the density interface ( Figure 7a).
The fate of the density interface is intimately tied to the mixing regime; scouring tends to sharpen density interfaces, while overturning tends to mix the interface into a more diffuse state. The mixing regime can be quantified by considering the budget for the buoyancy frequency in sorted density coordinates Taylor and Zhou, 2017): Here z * is the height of a vertically sorted isopycnal, adiabatically sorted to monotonically decrease with depth; N 2 * is the buoyancy frequency, Eq.
(3), in sorted height coordinates; and ρ is the potential density, where ρ 0 is a reference density. The 'effective' diffusivity is (Winters and D'Asaro, 1996;Nakamura, 1996) where κ is the molecular diffusivity and · z * denotes a spatial average for a fixed z * . Note that this formulation assumes that density is controlled by one variable (temperature or salinity). The key advantage . Cartoon illustrating overturning (a) and scouring (b) regimes at a density interface (Caulfield, 2021).
of using sorted height coordinates is that κ e is strictly positive. A purely laminar flow has κ e = κ, and turbulence will lead to κ e > κ by distorting the isopycnals, increasing the density surface and thereby increasing the density gradients.
The relative sizes of the terms on the right hand side of Eq. (5) can be used to diagnose the mixing regimes. Terms B and C in Eq. (5) represent translation and diffusion of the sorted density profile, respectively, and since κ e > 0 this diffusion acts to spread out density interfaces. Term A in Eq. (5) can be positive or negative, and this term dictates whether turbulence at a density interface is in the scouring or overturning regimes. For example, consider a flow in the overturning regime where shear-driven turbulence is generated at a density interface (Figure 7a). If the flow above and below the interface is relatively quiescent, κ e could exhibit a local maximum at the density interface with ∂ 2 κ e /∂z 2 * < 0. In this case term A will act to reduce the stratification at the interface. On the other hand, in the scouring regime ( Figure 7b) strong stratification will suppress mixing at the interface where κ e will be relatively small. In this case strong mixing on either side of the density interface can result in a flow with ∂ 2 κ e /∂z 2 * > 0, in which case term A will act to increase N 2 * and sharpen the interface.
The scouring regime requires a source of turbulence external to the density interface. In natural stratified shelf seas, turbulence generated by bottom friction and surface forcing (e.g. wind, waves, and/or convection) generate turbulence below and above mid-water density interfaces. Turbulence associated with flow past OWTs could also play this role. Indeed, as discussed in Section 3.2, horizontal shear can be very effective at generating and maintaining layers in stratified flows. On the other hand, energetic turbulence driven by the flow past OWTs could be strong enough to overturn and mix natural density interfaces.
Thus, it is not clear whether OWTs would generate turbulence in the overturning or scouring regime. The distinction will likely depend on the strength of stratification, the speed of flow through offshore wind farms, and the geometry of wind turbines. As will be expanded on in the following section more research is needed to understand mixing generated by offshore wind infrastructure and their impact from the scale of OWTs to the shelf sea scale.

FLOW STRUCTURE INTERACTIONS
The dynamics of flows past a range of structures has seen significant research over the last 50 years, both due to immediate real-world applications and the increase in computational resources and experimental measurement fidelity (Williamson, 1996). Here work studying the dynamics of unstratified and stratified flow past infrastructure, relevant to the offshore wind sector, are integrated and reviewed.

Unstratified Flow
Flow past cylindrical structures are well studied, owing to their geometrical simplicity and vast engineering importance; such studies are directly analogous to flow past OWT foundations, such as monopiles. Unsteady vortices shed by cylinders can lead to vibration, acoustic noise, resonance, and ultimately structural failure. Shed vortices form coherent wakes that are spatially vast, and can be detected several hundred diameters downstream of their source, depending on the background flow conditions (e.g. turbulence properties). Cylinder wakes fundamentally alter flow conditions as a source of anthropogenic turbulence (mixing), particularly coherent along the cylinder axis. Dynamics are complex due to the interaction of at least three shear layers; the boundary layer, shear layer, and wake ( Figure 8). In unstratified waters, flow structure is dependent on the cylinder roughness, end conditions, freestream conditions, and the ratio of inertial to viscous forces, the cylinder Reynolds number Re d , where u ∞ denotes the freestream velocity, d the cylinder diameter and ν the fluid viscosity. The dependence of the flow structure on Re d arises due to transitions in the different shear layers. The most well known instability arising from the flow past a cylinder is the Karman vortex (KV) (Williamson, 1996), which develops for Re d 49 (below which the flow is steady and laminar). The Karman instability is associated with alternating 2D vortices shed from either side of the cylinder, aligned with the cylinder axis, and is a consistent feature of even high Reynolds number flows. The unsteady KV is characterised by the dimensionless frequency, the Strouhal number St = f KV d/u ∞ , where f KV is the frequency of vortex shedding. St varies with Re d depending on shear layer transitions, as shown in Figure 8c. The KV is a dominant feature of flow past a cylinder (Figure 9a).
The typical Re d of shelf sea currents past offshore wind foundations (Figure 2) is at least Re d 10 5 , estimated assuming a small 5 m diameter monopile with a minimum tidal velocities of 0.02 m s −1 (Vindenes et al., 2018). This minimum Reynolds number is sub-critical (see Figure 8), indicating that shear layer instabilities may be present, manifesting as Kelvin-Helmholtz (KH) type instabilities. Three dimensionality is a dominant feature of cylinder wakes at these high Reynolds numbers. Critical transition occurs between 2 × 10 5 < Re d < 5 × 10 5 , and is associated with boundary layer transition to turbulence which causes the separation point to occur further downstream on the cylinder surface. For a smooth cylinder, asymmetric separation-reattachment of the boundary layer either side of the cylinder causes a sudden increase in vortex shedding frequency (Figure 8). Super-critical flow is associated with Re d > 5 × 10 5 , where symmetric separation bubbles and turbulent boundary layers are present on both sides of the cylinder. Roughness effects are felt primarily in the critical transition regime; surface roughness causes earlier transition to turbulence and bypasses the asymmetric regime ( Figure 8b).
With distance x downstream, the far-wake x/d 50 is particularly sensitive to freestream conditions. The spectral energy of the KV (energy associated with f KV ) decays downstream of the cylinder and the wake width grows like √ x (under low levels of freestream turbulence) and is approximately self-similar (Ghosal and Rogers, 1997). A secondary vortex street emerges at a lower frequency than the KV, arising due to the merging of vortex pairs, or via hydrodynamic instability of the mean flow, depending on the flow Reynolds number (Jiang, 2021a). Coherent structures in the far-wake can also arise due to non-linear interactions between freestream structures and the KV. Cimbala and Krein (1990) and Williamson and Prasad (1993) found that the interaction between the KV and freestream waves could lead to resonant peaks in spectral energy associated to their non-linear interaction, indicative of hydrodynamic instability. These peaks can be detected far downstream of the cylinder, x/d > 300, indicating that non-linear dynamics have a large effect on the wake even far downstream. It is therefore unsurprising to observe non-linear interaction between wake effects from multiple monopiles in offshore wind farms (Figure 9). Under high freestream turbulence the wake spreads more rapidly, depending on the background turbulent intensity and the turbulent integral length scale. Eames et al. (2011) demonstrated that the growth rate of the wake increases from ∼ √ x to ∼ x when the wake deficit velocity is approximately equal to the background turbulent intensity, and the integral length scale of the turbulence is comparable to the cylinder diameter. When subject to high turbulence a cylinder wake will dissipate downstream more rapidly and diffuse into the background turbulence. In contrast, when background turbulence levels are low the turbulence generated in the wake of a cylinder can persist hundreds of diameters downstream. In current shallow-water offshore wind farms, where levels of turbulence are high, wakes have been observed at least 1 km in length (e.g. Figure 9). The formation of turbid plumes, where suspended sediment is trapped and transported in KVs and secondary vortex streets, is correlated with high levels of TKE (Grashorn and Stanev, 2016). As the wakes pass through the wind farm plumes are observed to spread and interact (Figure 9b). Wake spreading indicates that the effects of monopile wakes are not limited to a short narrow region downstream. Although KV are expected to decay for x/d ≈ 60 (Jiang and Cheng, 2019), observed plume interaction at large lengths, x/d > 150, suggests wake-wake interactions are important. However, conditions vary considerably between the current shallow coastal offshore wind farms and the deep water future development sites. It is to be expected that in deep water sites, with lower background turbulence, that wakes may be even larger.

Stratified Flow
Relatively few studies have investigated stratified flow interaction with vertically oriented cylinders, analogous to proposed offshore wind foundation deployment in seasonally-stratified shelf seas. This may be in part due to similarity between unstratified and stratified flow past vertical cylinders at very low Reynolds, Re d 45 (Meunier, 2012), where wakes are inherently two-dimensional. However, differences arise when three-dimensionality is present in the wake, which is certainly characteristic of the high Reynolds number flows associated with offshore wind farm infrastructure. Three-dimensionality in the wake of vertical cylinders is important as it can lead to fundamental reorganisation of stratified fluids. Layering can emerge where such a flow can lead to multiple intermittent regions of fairly constant density neighbouring thin rbulence sitting in an otherwise unperturbed and stably stratified environment. We ll such a localized turbulent region a turbulent cloud; we consider stratified turbulent ouds, where the cloud interior and the surroundings are initially linearly stratified. ratified turbulent clouds are governed by inhomogeneous dynamics, in contrast to the namics of homogeneous stratified turbulence. The latter class of problems has been vered extensively, with theoretical (see Billant   interfaces with steep density gradients (Bosco and Meunier, 2014). This layering behaviour may arise from any process that produces spatially periodic mixing in the vertical direction (Thorpe, 2016).
Cylinder wakes are prone to layering due to their increase in both horizontal shear, normal to the direction of stratification, and organised vertical vorticity, parallel to the direction of stratification, through the DRAFT formation of KV in particular. The susceptibility of a stratified flow, perturbed by moving cylinders, to layering is dependent on the mechanisms that generate layers, which vary considerably depending on the type of forcing (e.g. oscillatory versus continuous stirring, Thorpe, 2016). Experiments have demonstrated that layers can develop from the 'cloud' of stratified turbulence that results from dragging vertical cylinders through a stratified flow, equivalent to flow through a wind farm (Figure 10). In the experiments of Maffioli et al. (2014) intrusions developed from local regions of mixed fluid, which grew until their length scales were approximately balanced by l ∼ũ/N , whereũ represents a characteristic velocity scale of the turbulence. The buoyancy frequency may also be used to define a Froude number, the ratio of inertial to gravitational forcing, where L is some length scale over which flow velocity and N are averaged, Thus, in the experiments of Maffioli et al. (2014), packets of turbulent fluid grew until they collapsed under gravity (F r ∼ 1) and spread outwards as pancakes, triggering horizontally propagating internal waves. The outward spread of intrusions generates a layered density profile, where well-mixed intrusions are neighboured by thin regions of strong density gradient. It is important to note here that the cylinders are not critical to the reorganization of stratification; they only act as the source of turbulence, via horizontal shear, in the stratified flow. Addition of purely horizontal shear to simulations of stratified shear flow can lead to coherent vertical vortices (Basak and Sarkar, 2006), analogous to the KV and KH instabilities of a cylinder wake. The coherent vertical vortices exhibit pairing, tearing, and amalgamation, resulting in a complex braided vorticity structure that ultimately leads to vertical variability. Layers can subsequently develop as intrusions, coupled with internal waves (Basak and Sarkar, 2006), once again with a characteristic vertical scale of the order ofũ/N , when these two quantities are estimated appropriately.
Indeed, provided the flow Reynolds number is sufficiently large, there is accumulating evidence, that horizontal shear (and hence vertical vorticity) in vertically stratified fluids inevitably forms layers on this scale. Dating back to the first theoretical analysis of Billant and Chomaz (2000b), there is clear evidence that such flows are prone to a class of 'zig-zag' instabilities. These instabilities imprint theũ/N vertical scale on the flow, and can be connected directly to the inherently nonlinear layered structures that develop at finite amplitude (Lucas et al., 2017). Furthermore, several studies (Deloncle et al., 2008;Waite and Smolarkiewicz, 2008;Augier et al., 2015) have demonstrated that the breakdown of these vertical vortices in a stratified fluid introduce a new, inherently stratified route to turbulence in a stratified fluid, and hence substantially enhanced mixing.
Whilst work on layering has been restricted to small scale experiments and comparatively low Reynolds numbers (at least by oceanographic standards), it has provided oceanographers with vital information on the fine-scale density structure of the ocean. However, it is unclear how turbulence generated at high Reynolds number, by flows past offshore wind foundations, will interact with the essentially two-layer density profile of seasonally stratified shelf seas. In addition to the high Reynolds number, offshore wind foundations are a similar width to the typical thermocline thickness. This contrasts with previous work where turbulence was either generated by structures with length scales two orders of magnitude smaller than the density gradient length scale (Maffioli et al., 2014), or (as in the 'zig-zag' instability studies mentioned above) the unstable horizontal shear/vertical vorticity is embedded initially in a linearly stratified fluid with close to constant buoyancy frequency. Where KV shed from offshore wind foundations are of the same order diameter as the density gradient length scale, flow mixing and density profile reorganisation may be fundamentally different, and this generic stratified flow geometry with a range of key characteristic length scales is very poorly understood, even in highly idealised circumstances.
Expansion of the offshore wind sector to deeper waters is predicated on development of floating foundations, which are finite in depth and anchored to the seabed. Finite depth structures induce complex dynamics within stratified fluids, especially if structures intersect sharp density gradients, i.e the thermocline (Figure 2). A good example case study of the impact of finite depth obstacles on stratified flow, is flow past a horizontal cylinder. At high Froude numbers dynamics are similar to unstratified flow. But as the Froude number decreases several new regimes occur (Boyer et al., 1989). Under stratification internal waves can be generated not only by the structure but also the KV. This in turn can increase drag, and thus the amount of mixing, by up to 100% (Arntsen, 1996). The drag coefficient of spheres in stratified flows also varies by up to 100%, when the vortex shedding frequency tends to the buoyancy frequency (Cocetta et al., 2021). the initial 350 to 750 s, depending on the selected U ∞ , to avoid contaminating the results. Since this period is between 3τ and 8τ , the mean velocity at the end of the observational period is always close to the asymptotic value U ∞ , and we define our topographic Froude number using this velocity, i.e.
We explore the impact of the initial acceleration of the flow further in § 4.1.

Basic observations
We consider first an experiment with N = 1.75 rad s −1 and U ∞ = 14.8 mm s −1 giving F R = 0.13 and consider the internal wave field at t = 388 s, shortly before the fluid initially located by the discs reaches our observation window. By this time the uniform flow velocity is within a few per cent of its asymptotic value. The b x and b z fields (figure 3a,b, determined by the Synthetic Schlieren) give clear evidence of an internal wave field propagating upwards, away from the obstacle, as the waves are swept down stream. These waves are reflected from the free surface and then Figure 11. Stratified flow over a hill, taken from the experiments of Dalziel et al. (2011). Internal waves can be observed in the lee of the structure propagating vertically and horizontally, clear in both the horizontal buoyancy gradient b x = −∂ρ/∂x and vertical buoyancy gradient b z = −∂ρ/∂z. The wake is also particularly apparent in the b z field, demonstrating the existence of significant mixing). N is the background buoyancy frequency and R denotes the height of the hill. not necessarily aligned with the direction of stratification (e.g. jackets and semi-submersible foundations, Figure 2).
Through considering reflectional symmetry, flow past floating foundations may be analogous to oceanic and atmospheric flow past sea-mounts and hills. Stratified flow past such obstacles displaces fluid and leads to both vertically and horizontally propagating internal waves, as well as a significant downstream wake ( Figure 11). If an obstacle is wide compared to its height, most fluid impinging on the obstacle passes over it, otherwise flow can pass around it. For low Froude numbers, where L is scaled with obstacle depth, flow blockage by internal waves creates additional drag, which is a very effective momentum sink on the impinging flow. Flow blockage can increase the drag force by 1-2 orders of magnitude compared to unstratified flow over the same obstacle (Smith, 1978;Castro et al., 1990;Cummins et al., 1994). In all cases understanding drag increase is critical to evaluate both the hydrodynamic loads placed on foundations and anchors, and the mixing in stratified shelf seas, associated with the flow past offshore wind foundations. Such wave generation and propagation also suggests the possibility that the effect of flow structures can be both local and non-local, as the emitted waves may transport momentum flux significant distances until they 'break'.

Mixing of Stratified Shelf Seas by Offshore Wind Foundations
To date there has been only two limited studies observing offshore wind foundation induced mixing of stratified waters (Floeter et al., 2017;Schultze et al., 2020a). These studies have been restricted to developments in aperiodically stratified regions of freshwater influence, in shallow water (depths of approximately 40 m and 24 m, respectively). Neither studies have investigated the potential impacts on seasonally, or permanently, stratified shelf seas. Floeter et al. (2017) performed surveys on two wind farms in the German Bight, North Sea. Water property transects through the wind farm revealed a consistent weakening of stratification near the centre. Effects extended into the surrounding area by approximately half the diameter of an ambient tidal excursion. However, it was unclear how much of this was due to 'infrastructure' turbulence from turbine foundations rather than natural topological effects. In addition to reduced stratification, Floeter et al. (2017) also measured local upwelling at the edges of the wind farms, similar to those observed near islands in stratified waters (e.g. Simpson et al., 1982). Shallow water model studies of the effects of offshore wind farms on local oceanic circulation patterns have shown arrays of infrastructure induce strong horizontal shear in wind stress leading to local regions of upwelling and downwelling (Broström, 2008;Paskyabi and Fer, 2012), consistent with the observations of Floeter et al. (2017).
Field measurements were also taken by Schultze et al. (2020a) who measured the stratified wake from an offshore monopile at the leading edge, with respect to local flow, of the DanTysk wind farm (Figure,  12). The monopile wake spread to a width 10 times the (6 m) diameter of the monopole and had reduced the potential energy anomaly by up to 65 % at a downstream distance of x/d ≈ 20 (Figure 12). The full distance required to return to pre-monopile conditions was not captured by the survey, even after approximately 300 m, a distance over 50 times the 6 m monopile diameter. The survey, clearly demonstrates that turbulence generated by monopiles reduces stratification. Complementing the survey, Schultze et al. (2020a) used Large Eddy Simulations (LES) to model flow past a single monopile which was simulated under different levels of background stratification. TKE dissipation rate was found to be up to two orders of magnitude larger in the thermocline than when monopiles were not present. High TKE dissipation rate persisted far downstream of the cylinder and, at x/d > 40, was still an order of magnitude larger than without monopiles at the thermocline depth. Further, the TKE dissipation rate was greater than that generated by the bottom boundary layer.
During a period of stronger stratification, Schultze et al. (2020a) conducted a second survey at the opposite end of the DanTysk wind farm. The second survey was less conclusive than their previous measurements, finding that no clear signal from the wake could be separated from background variability. This may be because stratification was strong enough to suppress the growth and interactions of the wakes, although Schultze et al. (2020a) only sampled the wake at distances greater than 200 m (or 33.3 diameters) from the monopile. It is therefore unclear if this is a result consistent with the earlier survey in Figure 12.
Leading order arguments for the mixing induced by offshore wind foundations have been explored by Rennau et al. (2012) and Carpenter et al. (2016), here they are reviewed. Models start by assuming that the turbulence produced by foundations is equal to the power lost to drag. Over an arbitrary vertical layer of fluid, L, the power lost to drag is given by where c D denotes a drag coefficient and |u| 3 L the velocity magnitude cubed, averaged over the layer L. Whilst Rennau et al. (2012) and Carpenter et al. (2016) assumed that c D is a constant in reality it is variable and dependent on the Reynolds number, surface roughness and structure geometry. The drag coefficient is also likely a function of depth and time and will vary with background shear. Crucially c D is also dependent on the Froude number, especially for highly stratified flows (see Section 3.2). At DRAFT present there is no parametrisation of the Froude number dependence of the drag coefficient for vertical infrastructure in stratified flows. We therefore assume c D = 1 as per the "high drag" case of Carpenter et al. (2016), although this estimate does not account for the potential effects of the Froude number.
A wind farm comprises many turbines, potentially supported by a range of foundations. Turbines are separated by a distance of approximately 10D (Howland et al., 2019), where D is the rotor diameter. The average production of TKE per unit volume P D over a layer of thickness L is therefore: Here it is assumed that the power lost to drag at a turbine foundation acts to increase TKE production over the full area occupied by a foundation. In reality this production will be localised to the comparatively narrow wake, and will vary by several orders of magnitude in space.
Here the potential production of TKE P D per unit volume, generated by the addition of a 10 m diameter monopile, is estimated using the natural conditions, averaged over the thermocline thickness, at the CaNDYFLOSS experiment site (Figure 3). Power generated by flow past such a potential foundation is strongly dependent on the flow velocity, with low frequency oscillations associated with neap-spring tide cycles and higher frequencies associated to daily tidal cycles ( Figure 13). The potential power added to the thermocline, P D , varies between approximately 2 × 10 −6 W m −3 to 4 × 10 −4 W m −3 with an average value of 5.9 × 10 −5 W m −3 . Over the same period, at the same site, the natural background dissipation rate (ε B ) at the thermocline is 3.5 × 10 −5 W m −3 (Scannell et al., 2021). In this approximate calculation, the additional (average) turbulence production is therefore 69% greater than background dissipation without monopiles present: Under several very strong assumptions, the relation between TKE produced by offshore wind foundations and how it is balanced by viscous dissipation and work performed on the buoyancy field may be expressed, from Eq. (2), by where subscript D represents contributions from structures (drag) and subscript B represents all other (background) contributions to TKE production (P) and dissipation rate (ε).
For simplicity, the turbulent flux coefficient can be assumed equal for each component of B and ε such that Γ = B/ε = B B /ε B = B D /ε D . In deriving (11) it is also assumed that OSW infrastructure contributes linearly to P, and is simply related to the structure-induced TKE dissipation rate by P D = (1 + Γ)ε D . In the present example we have shown that P D ≈ 1.69ε B , such that Under the simplest conventional assumption that Γ = 0.2, as typical at the pycnocline in seasonally stratified waters (see Section 2), we obtain ε D ≈ 1.4ε B , such that the total TKE dissipation rate, ε = ε B + ε D , is at least 140 % higher at the thermocline when OSW infrastructure is present. As noted, several sweeping assumptions have been made to arrive at this estimate, which inevitably has a large amount of implicit uncertainty. It is reasonable to suppose that this estimate is likely to be a conservative lower bound on the effect of OSW infrastructure on turbulent dissipation and mixing at the thermocline for at least two reasons. The first is that there is clear potential for a strongly nonlinear effect of such infrastructure on the dissipation rate. Secondly, setting Γ 0.2 may well be an under-estimate of the vigorous overturning mixing likely to be triggered in the wake of such infrastructure (Caulfield, 2021). Nevertheless, these simple conservative estimates show that foundations produce turbulence at levels that will clearly affect the leading order balance of TKE transport, even when normalised by the total area between offshore wind foundations. Thus, offshore wind has the potential to impact directly the stability of seasonally-stratified shelf seas ( Figure 5). The dissipation rate estimate indicates that OSW infrastructure is likely to affect the leading order balance of TKE transport in vicinity of the windfarm, but the implications of this on regional-scale fluid dynamics is more challenging to evaluate. To address this, Carpenter et al. (2016) derived similar models for TKE production as (10), and constructed arguments based on estimates of the different timescales in the flow. In particular, timescales of OSW-induced mixing, and timescales of flow parcels convecting through DRAFT the windfarm. By comparing different estimates of these timescales Carpenter et al. (2016) theorised the extent that the water column would mix as a parcel passed through the array. It was concluded that small developments in shallow waters (Bard 1 and Global Tech 1, German Bight) were unlikely to affect stratification, where TKE production from the bed and free surface dominate. However, despite the simplification of the models and assumptions made, large scale developments in deep water were recognised to have potential for significant impact (Carpenter et al., 2016). Rennau et al. (2012) adopted regional-scale numerical models to assess the effects of offshore wind development in the Baltic Sea, where tidal currents are minimal and stratification are driven by dense saline currents beneath fresh water. They modified a two-equation (RANS) turbulence closure scheme to capture the enhanced mixing arising from OSW installations and found that current installations did not have a significant effect on regional stratification (though cautioned that future development, covering more of the Western Baltic Sea, could lead to significant impacts, such as reduced bottom salinity). However these closure schemes are known to be overly dissipative in stratified flows (Hewitt et al., 2005), and demonstrated to generate overly diffuse pycnocline structures in realistic simulations (Luneva et al., 2019). These schemes are therefore improperly conditioned, at present, to address the impact of OSW infrastructure in stratified seas.

DISCUSSION
The scale of planned offshore wind energy industry is much greater than past, and existing, sea use. Installed offshore wind capacity will increase by 600% in the next decade (Figure 1a), requiring an extra ∼ 20, 000 10 MW+ turbines. Fixed, and floating, offshore wind infrastructure will penetrate the thermocline, adding 'anthropogenic' mixing on top of natural mixing. The impact of such infrastructure will be fundamentally different from existing sea use. For example, even large surface vessels have comparatively small drafts, ∼ 10 m (Golbraikh and Beegle-Krause, 2020;Nylund et al., 2020). Further, whilst offshore oil and gas platforms are similar to offshore wind infrastructure, spanning the thermocline, only ∼ 6, 500 platforms have been installed, globally, over the last 75 years (Schneider and Senders, 2010). Therefore, planned offshore wind developments will add new and large scale infrastructure sources of turbulent mixing in seasonally stratified seas. This discussion reviews our current knowledge gaps, frames potential impacts on shelf sea dynamics and thus marine ecosystem functioning and highlights routes for sustainable growth of the offshore wind sector.

Infrastructure Mixing
The mixing of stratified waters by offshore wind infrastructure is poorly understood. As evidenced in Section 3.3 there is a dearth of research on high Reynolds number stratified flow past vertical structures, which is vital for understanding and parameterising 'infrastructure' mixing processes in natural environments. The problem is particularly complex due to its scale; laboratory experiments and fully resolved numerical simulations are limited to relatively low Reynolds number flows, certainly by comparison to the real oceanographic flows. In all scenarios of fixed or floating OWTs, fine scale vertical density structure and enhanced mixing is anticipated due to horizontal shear generated by the flow past the infrastructure. The O 10 −2 − 10 0 m length scales associated with the resulting turbulence will be unresolved by low resolution numerical simulations. Wakes from individual structures may persist for O 10 2 − 10 3 m downstream. From individual wind turbines to a single offshore wind farm, infrastructure adds a wide range of length scales of O 10 1 − 10 4 m. Further, multiple multiple wind farms are distributed at shelf-wide scales. The vast range of length scales present in flow past offshore wind farms necessitates a variety of modelling techniques.
Although coarse numerics, i.e. Large Eddy Simulation (LES) or other more sweeping turbulence closure schemes, may be desirable modelling strategies for high Reynolds number flow (Rennau et al., 2012;Carpenter et al., 2016), without a firm understanding of the underlying physics, and robust datasets for validation, low resolution numerical models should be adopted with care, particularly for stratified (and therefore highly anisotropic) flows (Hewitt et al., 2005;Khani and Waite, 2015). Integration of modelling disciplines ranging from Direct Numerical Simulation (DNS) of idealised flows to regional scale ocean circulation models is therefore required to determine the impact of offshore wind farms on seasonally stratified shelf seas. Field surveys are essential to support and validate physical and numerical studies. However, the limited work to date (Floeter et al., 2017;Schultze et al., 2020a) has suffered from DRAFT uncertainties, where either effects of infrastructure were difficult to discern from topographical effects or where wakes were difficult to separate from background variability. While both studies concluded effects of infrastructure on mixing may be large, further surveys are required to support this. In addition there is a clear need for repeat and 'before and after' surveys, as noted by van Berkel et al. (2020).
There are also key research questions regarding the geometry of offshore wind infrastructure, and how this impacts mixing. OWTs shed vortices with horizontal scales that are comparable to the thickness of the thermocline. Such mixing contrasts with previous work studying mixing by thin vertical cylinders, several orders of magnitude smaller in diameter compared to density length scales. It is vital to understand how large scale vortex structures interact with a relatively thin thermocline; mixing processes may well be fundamentally different to those studied before, with complex spatio-temporally variable vortex, turbulence and mixing dynamics. In addition, future floating technology raises further questions regarding the impact of geometry on mixing. It is expected that spar-buoy designs ( Figure 2) will act in a similar way to monopiles, given they penetrate through the thermocline and into the well-mixed deep water. However, semi-submersible (or any other small-draft designs) will introduce non-trivial effects by intersecting the thermocline. Research of stratified flow over finite topography, similar to small draft floating OWT, has demonstrated that baroclinic effects can enhance drag by up to two orders of magnitude (Section 3.2). Here, floating structures will introduce infrastructure mixing via shed lee waves, internal waves, blockage effects, and wake-wake interactions in the case of semi-submersible designs.
Crucially, a lack of insight into key multi-scale mixing processes adds uncertainty to current attempts to quantify the impact of infrastructure mixing on shelf sea dynamics and ecosystem functioning. For example, the turbulent flux coefficient is often assumed constant, Γ = 0.2 (Section 2.3 and 3.3), but it is unknown if this holds in the wake of OWTs where fine scale density structures and strong spatial variability are present. Further, production of TKE due to infrastructure has been assumed constant in time and evenly distributed over the area 'occupied' by the monopile. In reality production is focused in the narrow wake of individual monopiles, and could vary by two orders of magnitude during the spring-neap tidal cycle (Section 3.3). TKE production, and thus mixing, arising from stratified flow past infrastructure is also dependent on the drag coefficient c D , yet little is known about the dependence of c D , particularly at high Reynolds numbers or where mixing length scales are large in comparison to density length scales. Meanwhile, research on drag past other obstacle forms suggest estimates for drag may be incorrect by orders of magnitude. Advancing our understanding of each of these processes is vital for assessing their impact on shelf sea oceanography and ecosystem functioning.

Shelf Sea Dynamics
Offshore wind farms are anticipated in having large local impacts on shelf sea dynamics, in a similar fashion to natural topographically controlled mixing, e.g. driven by internal waves and flow over seafloor sand banks. It is anticipated that broader and more diffuse thermocline would develop as a result of enhanced mixing, weakening it as a barrier to vertical mixing and transport ( Figure 14). Subsequent change to surface water characteristics would likely alter exchange across the ocean-atmosphere interface, with impacts on heat storage, atmospheric CO 2 uptake and benthic resupply of O 2 . The scale of this response will be site and infrastructure specific. At regional scales the water column should re-stratify subject to natural buoyancy forcing. In a more extreme scenario, strongly enhanced mixing could prevent stratification from forming around wind farms. This effect is observed around islands, where reduced surface temperatures and well-mixed waters result from enhanced mixing due to flow acceleration and seabed shoaling (Simpson et al., 1982). Small, but significant residual currents sweep this well-mixed water into an observable downstream wake.
Enhanced mixing from infrastructure may also impact seasonal, and shorter timescale, cycles. The first order response of the vertical density structure to enhanced mixing in a wind farm region would likely be delayed onset and early breakdown of seasonal stratification with weaker stratification throughout the summer season. Development of near surface stratification during periods of low wind stress would no longer be expected to occur as the enhanced mixing would act to persistently stir the normally episodically mixed surface layer.
With infrastructure development from local scale of a single turbine to regional scales of multiple developments offshore wind farms in shelf seas have the potential to have regional impact. Horizontal variation in density may arise from wind farm scale mixing, and these variations may enhance submesoscale DRAFT processes driving additional vertical transport and mixing. The density of offshore wind farms, the regional distribution of mixing and wake-wake interactions between wind farms will be of critical importance in determining shelf sea response to offshore wind development.

Shelf-Sea Ecosystem Functioning
Enhanced mixing rates due to infrastructure, would not only lead to temporal and spatial variations in vertical density structure, but also impact biogeochemical function at a fundamental level that would cascade up the ecosystem.
Stronger mixing in the thermocline will drive more nutrients from the bottom water up into the subsurface chlorophyll maximum and, if the mixing is strong enough, up into the surface layer where it could support additional phytoplankton growth. The rate of turbulent mixing strongly affects the simulation of ecosystem behaviour, as demonstrated by (Luneva et al., 2019), who found that different mixing schemes caused a shift in the spring bloom by 1 month, and change regional chlorophyll differences by order 100% .
Mixing also alters the light experienced by the phytoplankton, with stronger mixing potentially disrupting light sufficiently to hinder photosynthesis. The net effect, i.e. whether or not the extra mixing aids net phytoplankton production, will depend on some balance between the nutrient and light effects. The summer reduction in bottom water oxygen concentrations will also respond to the increased mixing. Oxygen will be supplied from the surface water downward, potentially offsetting some of the normal bacterial demand for oxygen as they recycle the sinking organic detritus from phytoplankton growth. At the same time, however, there may be an increased supply of organic detrital material if the net effect of mixing on primary production was positive, in which case the bacterial demand for oxygen in the bottom water will increase. Understanding the net effect of mixing on bottom water is important, as there are large areas of stratified shelf seas that are currently viewed as being close to experiencing oxygen depletion in late summer (Ciavatta et al., 2016).
Changes to biogeochemical functioning would need to be assessed over several years. For instance, an immediate positive impact on net phytoplankton production because of the extra nutrient supply will mean that the total amount of the water column inventory of nutrients used for that year will have increased. The shelf is not restocked with fresh nutrients from the open ocean every year (Ruiz-Castillo et al., 2019), so the ability of the shelf system to maintain the increased production will depend on how efficient the system is at recycling organic nutrients back to inorganic nutrients particularly over the winter.
A natural analogue to the effects of renewable energy infrastructure in a stratified environment may be found in the central Celtic Sea, where a number of seabed banks interact with tidal flows to inject significant internal wave-driven mixing into the thermocline (Palmer et al., 2013). While the exact biogeochemical and biological reasons are not yet clear, it does appear that the bank-driven mixing increases the overall biological activity of the region, ultimately resulting in commercial fishing (see, e.g., Figures 1c and4 Sharples et al., 2013). The island mixing effect too (Simpson et al., 1982), produces significantly enhanced nutrient fluxes, with a corresponding increase in plankton production and fishing activity. Thus, combined with reduced fishing pressure, appropriate infrastructure in suitable locations could prove a positive impact on regional ecology, with benefits to wildlife and fisheries.
The changes in the timing of stratification expected from increases in mixing also need to be considered. Late development of stratification and the spring bloom has the potential to impact wildlife, which has evolved to take advantage of this abundance, with seabird breeding populations and fish stocks a notable risk. At the other end of the stratified period, additional mixing will lead to earlier re-mixing of the whole water column, shortening the total productive time of the area and also altering the timing of the autumn bloom. Similarly, some seabird colonies are located to take advantage of the enhanced productivity at tidal mixing fronts (Trevail et al., 2019), the location of which may be affected by additional mixing from offshore wind farm infrastructure.

Management and Mitigation Requirements
Environmental Impact Assessments (EIA) are required for all new offshore wind farms to mitigate negative impacts resulting from construction, operation and decommissioning across the lifecycle of use. Surveying also covers geophysical site assessment, focus on seabed composition, current assessments and the potential for seabed scour around anchors and cables. Current EIAs have been developed for mixed  Figure 14. Offshore wind infrastructure adds wake turbulence throughout the upper water column, including directly at the thermocline. Here swirl size indicates turbulence intensity arising from near bed, near surface and flow-structure induced shear. Wake turbulence mixes cold nutrient rich bottom water with warm nutrient poor surface water, reducing the strength of stratification and potentially enhancing plankton growth in the subsurface chlorophyll layer. Changes in the subsurface chlorophyll layer would have further impacts on nutrient pathways, ecosystem functioning and oceanic carbon sequestration.
coastal environments to ensure sustainable growth of the offshore wind sector. Beyond the construction phase, impact has been focused on the benthic habitat and individual species interactions with infrastructure. Impact on marine mammals, seabird colonies, and fisheries in the region are individually assessed, but without consideration of potential alteration of the primary production through enhanced mixing.
To ensure the continued growth of the sector (The Crown Estate, 2020), impacts of the new generation of developments in deeper seasonally stratified regimes will likely require a more fundamental assessment. Baseline surveys must include the natural cycle of water column stratification, biogeochemical fluxes, and primary production. Accurately forecasting the interaction between the flow, infrastructure and stratification will require site, array and design specific observations and model scenarios. Only with a comprehensive understanding of this physical modification to the biogeochemical functioning of shelf seas, can impact throughout the marine web be adequately assessed.

CONCLUSIONS
Previous work has considered the environmental impact that offshore wind energy has on well-mixed shallow-water marine ecosystems, including from benthic habits, fisheries to seabirds. Much of this work remains relevant to enable sector growth. However, sector expansion from well-mixed shallow water to seasonally stratified deeper water represents a fundamental change, where physical and environmental impacts are not understood.
For the first time planned developments of both fixed and floating offshore wind infrastructure will add large scale anthropogenic mixing to seasonally stratified shelf seas. Large scale mixing may force shelf sea physics, establishing a 'new normal' for biogeochemical cycles and shelf sea ecosystem functioning. The potential benefits and risks posed by infrastructure mixing of stratified shelf seas, on top of climate change, represents a combined hazard that has not been considered.

Infrastructure Mixing of Shelf Seas
Locally, flow past offshore wind energy infrastructure results in (barotropic) drag and turbulence that, by itself, dissipates at least 140% more energy than exists naturally at the thermocline (section 3.3). However, in stratified waters, the additional baroclinic (wave) drag can exceed the barotropic drag. For example, the baroclinic drag for a stratified flow past topography can be 1-2 magnitudes larger than the barotropic drag (Smith, 1978). Baroclinic drag from vertical infrastructure in stratified flow, such as OWT, is as yet unquantified. Moreover, the role of horizontal shear on vertical mixing, here produced by obstacles that scale with density gradient length scales, is poorly understood, particularly in complex stratifications like the thermocline. It is of great importance to extend our understanding of the fundamental fluid dynamics of flows past vertical structures in ocean-realistic stratifications, in particular, the onset of turbulence and ensuing mixing associated with the breakdown of induced vortices.
Regionally, the first order paradigm for seasonally stratified shelf seas is the balance between the stratifying influence of surface heating, and the input of mechanical energy to mix the water column at the upper and lower boundaries (due to wind stress and the tidal shear respectively). High resolution shelf sea models have some success in reproducing this (Luneva et al., 2019). However, the addition of mixing from large scale offshore wind farm development, limits our ability to understand the trajectory of shelf sea ecosystems. To address this, research is urgently needed that scales processes from: a single turbine; an array of turbines composing a wind farm; to an entire shelf sea region with multiple farms. Advances in regional ecosystem modelling must then be validated against direct before-and-after observations to skillfully assess the direct and indirect impacts of anthropogenic mixing, and so guide sustainable development.
Growth of the offshore wind energy industry must be accelerated to meet global 2050 Net Zero commitments. Risks posed by offshore wind development in stratified shelf seas should be mitigated against, but potential benefits should be identified and maximised. Research should consider how nature based solutions, foundation design and array layout can enable sector growth, mitigation of risks and maximisation of benefits.