INTRODUCTIONTop
Upwelling of deep water, caused by surface wind stress along the coast, occurs naturally in some oceanic regions and favours high primary production. High primary production regions require both large solar radiation input and a high inorganic nutrient content. The above two requirements seldom occur naturally, and either solar radiation or nutrient concentration becomes a limiting factor inhibiting primary production. High solar radiation and nutrient availability occur simultaneously in surface tropical and subtropical ocean regions where natural upwelling along coastal waters occurs. The scarcity of these natural upwelling regions in the whole ocean and the need to obtain fish resources led to the start of aquaculture activities, first on land and, more recently, in coastal marine environments (CostaPierce 2002CostaPierce B.A. 2002. Ecological Aquaculture: The Evolution of the Blue Revolution. Blackwell Science, 382 pp.). However, these activities do not seem to be sufficient to satisfy the fish requirements of an increasing world population (FAO 2010FAO. 2010. The state of the world fisheries and aquaculture. Rome: Food and Agriculture Organization of the United Nations.). The message that comes out most clearly from FAO (2010)FAO. 2010. The state of the world fisheries and aquaculture. Rome: Food and Agriculture Organization of the United Nations. is the concern that the world fisheries are on a dangerous course (Pauly and Froese 2012Pauly D., Froese R. 2012. Comments on FAO’s State of Fisheries and Aquaculture, or ‘SOFIA 2010’. Mar. Pol. 36: 746752.).
Scientists are now looking for ways to avoid these limiting factors and fertilize the ocean artificially at a large spatial scale (UNESCO 2010UNESCO 2010. Ocean fertilization. A scientific summary for policy makers. 17 pp.). The activity is known as mariculture (e.g. McKinley and Takahashi 1991McKinley K.R., Takahashi P.K. 1991. Deep ocean water, artificial upwelling, and open ocean mariculture: A promise for the future. IEEE 1: 195199.), in close relation to aquaculture, which refers to similar fish production activities carried out at smaller spatial scales in land factories. Several attempts at ocean fertilization have been tested. Experiments based on the addition of iron or phosphorus to ocean sites seem to have been concluded negatively (de Baar et al. 2008de Baar H.J.W., Gerringa L.J.A., Laan P., et al. 2008. Efficiency of carbon removal per added iron in ocean iron fertilization. Mar. Ecol. Prog. Ser. 364: 269282. , UNESCO 2010UNESCO 2010. Ocean fertilization. A scientific summary for policy makers. 17 pp.). Alternative ways to fertilize the ocean are based on the artificial upwelling of nutrientrich deep waters. This paper proposes an energetically efficient mechanism to upwell deep waters using offshore wind energy.
Several artificial upwelling mechanisms have already been proposed. One of them is to cause upwelling using surface wave energy (Kirke 2003Kirke B. 2003. Enhancing fish stocks with wavepowered artificial upwelling. Ocean Coast. Manag. 46: 901–915.). Wavepowered inertia pump theory (Liu 1999Liu C. 1999. Research on artificial upwelling and mixing at the University of Hawaii at Manoa. IOA Newslett. 10: 18.), using a simple pump design originally conceived by Isaacs et al. (1976)Isaacs J.D., Castel D., Wick G.L., 1976. Utilization of the energy in ocean waves. Ocean Eng. 3: 175187., predicts a volume transport of 0.45 m^{3} s^{−}^{1} from a depth of 300 m employing a 1.2m diameter rigid tube. Using a “Bristol Cylinder” (Davis 1990Davis J. 1990. Wave energy absorption by the Bristol Cylinderlinear and nonlinear effects. Proc. Inst. Civil Eng. 89: 317340.) of 21.5 m diameter and 40 m length, Kirke (2003)Kirke B. 2003. Enhancing fish stocks with wavepowered artificial upwelling. Ocean Coast. Manag. 46: 901–915. estimated that a power of 1328 kW could be obtained in 1.9m Hawaiian swells and pointed out that such a power would be able to raise 50 m^{3} s^{−}^{1} of water from a depth of 500 m. A recent attempt (White et al. 2010White A., Björkman K., Grabowski E., et al. 2010. An open ocean trial of controlled upwelling using wave pump technology. J. Atmos. Ocean. Tech. 27: 385396.) used an Isaacs wave pump mechanism consisting of a 300mlong flexible tube with a oneway valve at the deeper end and the upper end attached by a 15m cable to a surface buoy. However, this experiment was interrupted early due to material failure caused by the vertical shear and stress of water currents on the tube’s flexible material (Kithil 2011Kithil P.W. 2011. Comments on “An open ocean trial of controlled upwelling using wave pump technology”. J. Atmos. Ocea. Tech., 28: 847849. , White et al. 2011White A.E., Letelier R.M., Björkman K.M., et al. 2011. Reply. J. Atmos. Ocean. Tech. 28: 850851.). Horizontal pressure differences inside and outside the tube may also have played a role in the material failure. A more complex and robust design was proposed by Salter (2009)Salter S.H. 2009. A 20 GW thermal 300metre^{3}/sec waveenergised, surgemode nutrient pump for removing atmospheric carbon dioxide, increasing fish stocks and suppressing hurricanes. Proceedings of the 8th European Wave and Tidal Energy Conference, Uppsala, Sweden.. This device uses sea wave energy to pump surface warm water downwards through a downtube reaching the thermocline. The warm water will mix with cold, nutrientrich water, rise to the neutral density layers, and spread sideways. For a wave climate of typical tradewind sites, and a design including a 100mdiameter, 200mlong tube could generate a volumetric flow rate of about 300 m^{3} s^{−}^{1}.
An upwelling mechanism based on the particular temperature and salinity vertical gradients of some tropical and subtropical ocean regions was proposed by Stommel et al. (1956)Stommel H., Arons A.B., Blanchard D. 1956. An oceanographical curiosity: the perpetual salt fountain. DeepSea Res. 3: 152153.. In these regions the temperature and salinity of sea water are higher in surface layers than in deep layers. The density profile remains statically stable due to the large temperature differences, thus counteracting the effect of salinity gradients on the density profile. The mechanism, named the perpetual salt fountain, consists in initially raising cold and lowsalinity water along a vertical tube that allows heat exchange. As water raises diabatically, its temperature increases due to the warmer environmental water outside the tube. Since the low salinity remains unchanged during the ascent, the water becomes lighter than the environmental water at every depth z, resulting in a perpetual flux of deep water at the upper end of the tube, close to the sea surface. Due to its low efficiency, the perpetual salt fountain was initially considered only an “oceanographic curiosity” (Stommel et al. 1956Stommel H., Arons A.B., Blanchard D. 1956. An oceanographical curiosity: the perpetual salt fountain. DeepSea Res. 3: 152153., Groves 1958Groves G.W. 1958. Flow estimate for the perpetual salt fountain. Deep Sea Res. 5: 209214.). However, this mechanism was recently tested in the ocean (Tsubaki et al. 2007Tsubaki K., Maruyama S., Komiya A., et al. 2007. Continuous measurement of an artificial upwelling of deep sea water induced by the perpetual salt fountain. DeepSea Res. I 54: 7584., Maruyama et al. 2011Maruyama S., Yabuki T., Sato T., et al. 2011. Evidences of increasing primary production in the ocean by Stommel’s perpetual salt fountain. DeepSea Res. I 58: 567574.) following an earlier experiment by Maruyama et al. (2004)Maruyama S., Tsubaki K., Taira K., et al. 2004. Artificial upwelling of deep seawater using the perpetual salt fountain for cultivation of ocean desert. J. Oceanogr. 60: 563568. using a 300mlong, 0.3mdiameter tube whose upper end was located at about 60 m depth. The estimated volume transport was approximately 5.2×10^{−}^{4} m^{3} s^{−}^{1} ≃ 45 m^{3} day^{−}^{1}.
Artificial upwelling has also been proposed and developed in Norwegian fjords. Upwelling in fjords has been used to prevent the development of toxic dinoflagellates in the upper brackish layer during summer and to favour toxinfree, highquality mussels (Berntsen et al. 2002Berntsen J., Aksnes D.L., Foldvik A. 2002. Production enhancement by artificial upwelling: a simulation study. Hydrobiology 484: 177190., Aure et al. 2007Aure J., Strand Ø., Erga S.R., et al. 2007. Primary production enhancement by artificial upwelling in a western Norwegian fjord. Mar. Ecol. Prog. Ser. 352: 3952., McClimans et al. 2010McClimans T.A., Handå H., Fredheim A., et al., 2010. Controlled artificial upwelling in a fjord to stimulate nontoxic algae. Aquac. Eng. 42: 140147. ). In one case, artificial upwelling was caused by forcing down through a vertical pipe (27 m long and 1.25 m in diameter) surface brackish water using an electrical pump (60 kW), which produced a downward volume transport of about 2 m^{3} s^{−}^{1} (Aure et al. 2007Aure J., Strand Ø., Erga S.R., et al. 2007. Primary production enhancement by artificial upwelling in a western Norwegian fjord. Mar. Ecol. Prog. Ser. 352: 3952.). Artificial upwelling starts outside the pipe due to the turbulent entrainment of the ascending buoyant plume with the environmental nutrientrich deep water. In another case, artificial upwelling was caused by a similar vertical plume but forced this time by a 100mlong air bubble curtain submerged at about 40 m depth (McClimans et al. 2010McClimans T.A., Handå H., Fredheim A., et al., 2010. Controlled artificial upwelling in a fjord to stimulate nontoxic algae. Aquac. Eng. 42: 140147. ). The estimated vertical flux of deep seawater was 65 m^{3} s^{−}^{1}. A third approach to force the upwelling of the unstable buoyant plume (McClimans et al. 2010McClimans T.A., Handå H., Fredheim A., et al., 2010. Controlled artificial upwelling in a fjord to stimulate nontoxic algae. Aquac. Eng. 42: 140147. ) used a 26 m^{3} s^{−}^{1} discharged of freshwater at 40 m depth from a hydropower plant. A diffuser plate at the lower end of the pipe was employed to increase turbulent entrainment. In this case the field experiment showed an entrainment of 117 m^{3} s^{−}^{1} of nutrientrich deep seawater. These experiments show an increase in nutrients in the euphotic zone resulting in an increase in phytoplankton and a decrease in toxic algae.
Related also to the upwelling mechanism proposed here is the Japanese TAKUMI project, which was set up in Sagami Bay, Japan in 2003, a prototype of the Ocean Nutrient Enhancer (Ouchi et al. 1998Ouchi K., Yamatogi T., Kobayashi K., et al. 1998. Density Current Generator—A new concept machine for agitating and upwelling a stratified water. Ocean Community Conference, 1, Baltimore, MD, USA, The Marine Technology Society, 129136.). This device consists of a submersed spartype floating structure and a steel riser pipe that upwells dense water from 200 m depth and mixes it with light surface water so that the mixed water spreads at its neutral buoyancy depth within the euphotic layer. The resulting volume transport is about 10^{5} m^{3} day^{−}^{1}≃1.2 m^{3} s^{−}^{1} (Ouchi 2007Ouchi K. 2007. Three years operation of ocean nutrient enhancer TAKUMI in Sagami Bay. Seventh ISOPE Ocean Mining Symposium, 16 July, Lisbon, Portugal, International Society of Offshore and Polar Engineers.). The source of energy of this Ocean Nutrient Enhancer is a diesel generator of about 100 Kw, which makes it unsustainable, though other design possibilities based on ocean thermal energy could overcome this problem (Otsuka 2011Otsuka K. 2011. Inclusive Impact Index Triple I and its application for Ocean Nutrient Enhancer. Proceedings of the Twentyfirst International Offshore and Polar Engineering Conference, International Society of Offshore and Polar Engineers.).
In this article we propose the use of offshore wind as the power source for artificial upwelling. Wind turbine farms are stably operating on land and coastal areas. However, offshore wind energy (Henderson et al. 2003Henderson A.R., Morgan C., Smith B., et al. 2003. Offshore wind energy in Europe – a review of the stateoftheart. Wind Energ. 6: 3552.) is being seriously considered as an emerging sustainable complementary source of energy. From the first offshore wind turbine deployed in 1990, cumulative installed power reached 1500 MW in 2009 mainly in northwest European waters (Zaaijer 2009Zaaijer M.B. 2009. Review of knowledge development for the design of offshore wind energy technology. Wind Energ. 12: 411430.). In the first six months of 2015, 584 commercial offshore wind turbines with a combined capacity totalling 2343 MW were fully grid connected in Europe. Once the 15 commercial wind farms under construction have been completed, the total capacity will be about 4268 MW (EWEA 2015EWEA. 2015. The European offshore wind industry  key trends and statistics 1st half 2015. Tech. rep. European Wind Energy Association.). Offshore wind farms have a higher cost in construction, deployment and maintenance than those on land or in coastal areas. However, the possibility of installing a large number of turbines with a minimum visual impact and the persistence of strong winds in some ocean regions make offshore wind farms an attractive solution to the electric energy demand. Interest in offshore wind farms may increase considering the possibility of using these turbines to favour ocean primary production.
The design proposals for offshore floating wind turbines depend on the way the turbine is stabilized (Jonkman 2009Jonkman J.M. 2009. Dynamics of offshore floating wind turbines – model development and verification. Wind Energ. 12: 459492., Jonkman and Matha 2011Jonkman J.M., Matha D. 2011. Dynamics of offshore floating wind turbines – analysis of three concepts. Wind Energ. 14: 557569. ). One of them, the sparbuoy floating wind turbine, incorporates a long vertical tube (typically 150 m long) submerged below the turbine’s tower that deepens the centre of buoyancy and stabilizes the turbine structure (Roddier and Weinstein 2010Roddier D., Weinstein J. 2010. Floating wind turbines. Mech. Eng. 132: 2832., Viré 2012Viré A. 2012. How to float a wind turbine. Rev. Environ. Sci. Biotechnol. 11: 223226.). It seems possible to combine a sparbuoy floating wind turbine with a wave energy converter, with the shape of a torus attached at the basis of the turbine’s tower so that wind and wave energy are simultaneously obtained at the same location (Muliawan et al. 2013Muliawan M.J., Karimirad M., Moan T. 2013. Dynamic response and power performance of a combined Spartype floating wind turbine and coaxial floating wave energy converter. Renew. Energy 50: 4757.). The idea that we explore in this work is to use the power of a sparbuoy turbine to raise deep water through an extended spar converted into a long (say 300 m deep) hollow tube (Fig. 1).
Dense deep water (with density ρ_{b}) inside a vertical submerged hollow tube open at the bottom and top ends will sink and, in the state of rest, will cause a negative hydraulic jump (say h) at the free surface, where pressure equals the atmospheric pressure, because the water column inside the tube is heavier than that outside. The amount of energy E required to pump a mass M of deep water to the surface is therefore the potential energy E=Mgh, where g is the acceleration due to gravity. The power P_{w} of the wind turbine and the upward volume transport F_{v} are therefore related through P_{w}=F_{v}ρ_{b}gh, or F_{v}=P_{w}/(ρ_{b}gh), which provides the upwelling transport through the tube as a function of the turbine’s power P_{w} and the surface hydraulic jump h. In the next section we calculate the surface jump h, using several background density vertical profiles ρ(z) and therefore provide an estimate of the upwelling transport. It is found that the order of magnitude of h is only 0.1 m, so, due to the large power of wind turbines (typically P_{w}=2 MW) the volume transport could be very large (of the order of 10^{3} m^{3} s^{−}^{1}). With such a large volume transport, pipe friction can only be neglected if very wide pipes, of the order of 15 m radius, are used. For thinner pipes, the head loss h_{f} becomes as important as the hydraulic head and needs to be taken into account. Still, considering pipe friction, the volume transport is very large, of the order of 500 m^{3} s^{−}^{1}. Using the turbulent entrainment theory of buoyant plumes it is found that once the dense water has upwelled to the sea surface, it will sink to a neutral height which, depending on the relative depths of the photic layer and thermocline layers, will usually be located below the photic layer, thus making the mechanism useless for primary production enhancement. To avoid this problem, some form of water dilution at the surface must be introduced. Water sprinkling, as a way to increase turbulent entrainment and water aeration (as a fountain) is considered, but this process is very energyconsuming. Other mechanical ways requiring no energy consumption, for example introducing a diffuser plate or multiport diffusers, as used for ocean outfall dilution in reverse osmosis desalination plants to reject brine, are also possible. We conclude that, though several important physical, ecological, and engineering questions must be addressed before its effectiveness is proven, artificial upwelling using offshore wind power is a promising approach for enhancing ocean primary production on a large scale.
THE SURFACE JUMP hTop
In this section we calculate the free surface height difference, or hydraulic head, h, between the water column inside and outside the tube for a number of density profiles. To do so we assume hydrostatic balance,
$\frac{\partial p}{\partial z}=g\rho$,

(1) 
where p(z) is the pressure field, g>0 is the acceleration due to gravity, ρ(z) is the density field and z is the vertical coordinate. We further assume that, in the steady state, the horizontal pressure gradient at the lower base of the tube is zero, so there is no horizontal flow acceleration in the exterior fluid as it enters the tube. Therefore, at the tube’s bottom depth, z=H (subscript b), the interior pressure p_{i}(H)=p_{ib} equals the exterior pressure p_{e}(H)=p_{eb} below the pump. At the free surface (top, subscript t) interior and exterior pressures are equal to the atmospheric pressure p_{a}, that is, p_{e}(0)=p_{i}(h)=p_{a}. We assume that in the steady state the water inside the tube is homogeneous and has a density equal to the bottom density ρ(H). We therefore neglect density changes as the fluid ascends along the tube. Under these assumptions integration of the hydrostatic condition (1) leads to
$\rho \left(H\right)\left(hH\right)={\int}_{H}^{0}\rho \left(z\right)dz$
, 
(2) 
which means that water masses per unit of horizontal area of the water column inside and outside the tube are equal. From (2) we obtain the solution for the jump h,
$h=H+\frac{1}{\rho \left(H\right)}{\int}_{H}^{0}\rho \left(z\right)dz$
,

(3) 
which depends on the length of the tube and on the vertical distribution of background density, ρ(z). Notice that H<0, and since 0<ρ(z)≤ρ(H) for z∈[H, 0], we have H≤h<0. We describe next three different solutions for h depending on three different simplified vertical density profiles, ρ(z), and as a fourth case we consider an experimental density profile.
Theoretical models of stratification
Twolayer model
The simplest case is that of a twolayer model (see Fig. 2). In the upper layer, where the tube is located, and in the lower layer, below the tube, densities are homogeneous and equal to ρ_{t} and ρ_{b}, respectively. In this case (see Appendix A) the solution for h is:
$h}_{2}=H\frac{{\rho}_{\text{b}}{\rho}_{\text{t}}}{{\rho}_{\text{b}}$.

(4) 
Thus, in the simplest case h_{2} depends linearly both on H and on the density difference ρ_{b}−ρ_{t}. We note that, as expected, h_{2}=0 when H=0 and/or ρ_{b}=ρ_{t}. Assuming values of H=−300 m, ρ_{b}=1029 kg m^{−}^{3}, and ρ_{t}=1026 kg m^{−}^{3}, typical of the western Mediterranean Sea in the summer season, we obtain h_{2}≃−87 cm. Notice that this profile is an extreme case in the sense that it maximizes the jump h. An upper density layer deeper than the tube’s depth will cause a zero jump, while a shallower layer will cause the weight difference between exterior and interior water column decrease, and therefore will cause a smaller h.
Linear density profile
This is the simplest case of a continuously stratified ocean, in which density ρ_{l} is assumed to be linearly dependent with depth, that is, ρ_{l}(z)≡ρ_{t}+ϑ_{0}z, where ρ_{t}=ρ_{l}(0) is the surface density, and ϑ_{0}<0 is a constant density stratification (see Fig. 2). Following the approach described above (see Appendix B) we obtain the surface jump h_{l},
$h}_{1}=\frac{{\vartheta}_{0}{H}^{2}}{2\left({\rho}_{\text{t}}+{\vartheta}_{0}H\right)}=\frac{H}{2}\frac{{\rho}_{\text{b}}{\rho}_{\text{t}}}{{\rho}_{\text{b}}}=\frac{{h}_{2}}{2$
. 
(5) 
The jump h_{l}=h_{2}/2 because the mass of the background water column (per unit area) in the linear density profile is half that in the twolayer profile. Thus, for the same typical values used previously, we obtain the jump h_{l}=h_{2}/2≃−44 cm.
Arctan density profile
Here we consider a density profile ρ_{a}(z) that is more realistic than the previous ones in the sense that it reproduces the ocean pycnocline. This profile may be modelled using the arctangent function
ρ_{a}(z)= c_{1} + c_{2} arctan(c_{3} + c_{4}z) , 
(6) 
where c_{1}, c_{2}, c_{3}, and c_{4} are given constants used to match the ρ_{a} asymptotes to the boundary values ρ_{t} and ρ_{b}, and to adjust the pycnocline depth (see Fig. 2 and Appendix C). Integrating (3) we obtain
$\begin{array}{c}{\rho}_{\text{b}}\left({h}_{\text{a}}H\right)=\frac{{c}_{2}}{{c}_{4}}\left[{c}_{3}\mathrm{arctan}\left({c}_{3}\right)\frac{1}{2}\mathrm{ln}\left({c}_{3}^{2}+1\right)\right]\\ {c}_{1}H\frac{{c}_{2}}{{c}_{4}}\left({c}_{3}+{c}_{4}H\right)\mathrm{arctan}\left({c}_{3}+{c}_{4}H\right)+\\ +\frac{{c}_{2}}{2{c}_{4}}\mathrm{ln}\left[{\left({c}_{3}+{c}_{4}H\right)}^{2}+1\right]\equiv \xi \left(H\right)\text{\hspace{0.17em}},\end{array}$ 
(7) 
where ρ_{b}(h_{a}−H) and ξ(H) stand for the masses per unit area of the water column inside and outside the tube, respectively. Solving for h_{a} we get
$h}_{\text{a}}=H+\frac{\xi \left(H\right)}{{\rho}_{\text{b}}$.

(8) 
For the same typical values in the western Mediterranean used previously we obtain h_{a}≃−18 cm, which has a magnitude smaller than the previous h_{2}=87 cm and h_{l}=44 cm.
Experimental density data
Finally, we obtain the water jump h_{e} using experimental in situ density profiles ρ_{j} (t, z_{i}) representative of the Catalan Sea (western Mediterranean). No offshore wind farm has yet been erected in the Mediterranean, but some studies have already addressed the suitability of climatological Mediterranean offshore winds for profitable energy applications (Lavagnini et al. 2006Lavagnini A., Sempreviva A.M., Transerici C., et al. 2006. Offshore wind climatology over the Mediterranean Basin. Wind Energ. 9: 251266.). The profiles were spatially averaged to produce a representative profile of the region, and timeaveraged to produce a representative profile of each season,
${\overline{\rho}}_{A}\left({z}_{i}\right)$
from January to March, ${\overline{\rho}}_{B}\left({z}_{i}\right)$
from April to June, ${\overline{\rho}}_{C}\left({z}_{i}\right)$
from July to September, and
${\overline{\rho}}_{D}\left({z}_{i}\right)$
from October to December. The sampling depths are z_{i}={−0, −10, −20, −30, −50, −75, −100, −125, −150, −200, −250, −300} m , i = 1, . . . , N = 12. These mean density profiles are plotted in Figure 3. In order to compute h_{e} for the four mean density profiles (subscript e={A, B, C, D}), we use the trapezoidal rule to numerically integrate (3), that is,
$h}_{\text{e}}=H+\frac{1}{{\overline{\rho}}_{\text{e}}\left(N\right)}{\displaystyle \sum _{i=1}^{N1}\frac{{\overline{\rho}}_{\text{e}}\left({z}_{i}\right)+{\overline{\rho}}_{\text{e}}\left({z}_{i}+1\right)}{2}}\Delta {z}_{i$,

(9) 
where the vertical grid size is nonuniform ∆z_{i}≡z_{i}−z_{i}+1. Using (9) and the data ρ_{e} (Fig. 3), we obtain the jumps h_{A}=−4.9 cm, h_{B}=−8.1 cm, h_{C}=−16 cm, and h_{D}=−12 cm, corresponding approximately to the four seasons.
In winter the reduced solar radiation causes a cooling of the upper ocean and therefore a vertical homogenization of the water (ρ_{A} in Fig. 3). The surface jump h_{A}=4.9 cm therefore reaches a minimum value. In spring the surface layer becomes warmer and therefore lighter, so vertical stratification increases. As a consequence, the surface jump increases to h_{B}=8.1 cm. In summer, with the pycnocline completely developed, vertical stratification reaches maximum values, and ρ_{C} becomes more similar to the arctangent profile ρ_{a} (Fig. 2). The surface jump increases to h_{C}=16 cm and approaches the jump size of the arctangent profile h_{a}≃ 18 cm. In autumn, surface cooling decreases vertical stratification and therefore the surface jump decreases to h_{D} = 12 cm. The cases described above show that the surface water jump h is only of the order of a few centimetres. This means that very large water volumes inside the tube, of the order of thousands of m^{3} s^{−}^{1}, will be transported up if the complete turbine power (also assuming a 100% pump efficiency) is used to raise the water. Such a huge amount of dense water on the surface runs a risk of sinking fast to aphotic layers due to vertical static instability (the so called fountain effect), ruining the enhancement of primary production. This possibility is analysed in a section below. We notice that the power required to raise the deep water is the potential energy associated with the water surface jump h only in the steady state, in which the vertical velocity inside the tube is steady. However, the water volume is initially forced upward from a state of rest, and therefore part of the wind energy must be used to increase the kinetic energy of the water inside the tube. This is a transient effect that, though interesting from the energy balance, is not addressed in this study, which only considers steady state energy balances. We have also assumed that the increase in kinetic energy is negligible and that the water withdrawn within the vertical tube is instantaneously replaced with deep water to keep the hydrostatic balance. This assumes that friction effects in the pipe are negligible, which in turn assumes, given the large volume flux involved, a very large pipe radius, R, so the vertical velocity inside the pipe is small. In the next section we consider the loss of power due to the increase in kinetic energy and the loss of power due to friction, which depend on the size of the pipe radius R.
LOSS OF POWER DUE TO THE INCREASE IN KINETIC ENERGY AND ENERGY DISSIPATION IN THE PIPETop
The loss of power due to the increase in kinetic energy is simply related to the velocity head h_{u}=u^{2}/(2g)= F_{v}^{2}/(2gπ^{2}R^{4}). We obtain the loss of power due to friction in the pipe, related to the head loss, h_{f} in hydraulic engineering, using the DarcyWeisbach equation
$h}_{\text{f}}={f}_{\text{D}}\frac{L}{D}\frac{{u}^{2}}{2g$,

(10) 
where f_{D} is the DarcyWeisbach friction factor, u is the fluid velocity, and L and D are the length and diameter of the pipe, respectively. Velocity head and head loss can be combined
${h}_{\text{u}}+{h}_{\text{f}}=\frac{{u}^{2}}{2g}\left(1+{f}_{\text{D}}\frac{L}{D}\right)$.

(11) 
The DarcyWeisbach friction factor f_{D} is obtained from the ColebrookWhite equation for Reynolds numbers Re>4000,
$\frac{1}{\sqrt{{f}_{\text{D}}}}=2{\mathrm{log}}_{10}\left(\frac{\epsilon}{3.7D}+\frac{2.51}{\mathrm{Re}\sqrt{{f}_{\text{D}}}}\right)$,

(12) 
where ε is the roughness height of the pipe. For very large Reynolds numbers, Re>10^{8}, corresponding to the expected highly turbulent flow inside the pipe, (12) simplifies to
$f}_{\text{D}}={\left[2{\mathrm{log}}_{10}\left(\frac{1}{3.7}\frac{\epsilon}{D}\right)\right]}^{2$.

(13) 
As a conservative value we take f_{D} = 0.03, which approximately corresponds to a very large ratio ε/D= 0.005, or a roughness height ε=5 cm for a diameter D=10 m. Typical roughness height of clean materials are two orders of magnitude smaller than this magnitude, but macro fouling caused by coarse matter, such as algae or mussels, will substantially increase ε.
When changes in kinetic energy and friction are considered, the power provided by the turbine P_{w} is used to save the hydraulic head h (now taken as a positive quantity), the velocity head h_{u}, and the head loss h_{f}, so the volume flux is
$F}_{\text{V}}=\frac{{P}_{\text{W}}}{g{\rho}_{\text{b}}\left(h+{h}_{\text{u}}+{h}_{\text{f}}\right)$

(14) 
Replacing h_{u}+h_{f} given by (11) in (14) and using F_{v} = πR^{2}u, we obtain the depressed cubic equation for the volume flux F_{v},
${F}_{\text{V}}^{3}\frac{{\rho}_{\text{b}}}{2{\pi}^{2}{R}^{4}}\left(1+{f}_{\text{D}}\frac{L}{2R}\right)+{F}_{\text{V}}g{\rho}_{\text{b}}h{P}_{\text{W}}=0$.

(15) 
For large pipe radius R, the first term in (15) may be neglected, and F_{v}(R) tends to the frictionless solution for slow flows F_{v}=P_{w}/(gρ_{b}h_{f}). The solution of F_{v}(R) is plotted in Figure 4.
As expected, the increase in kinetic energy and the loss of energy due to friction can only be neglected for a large radius R. Given a hydraulic jump h=0.15 m, a pipe radius R≃10 m, and a wind turbine providing a power F_{W}≃2 MW, the resulting volume transport is F_{v}≃547 m^{3} s^{−}^{1}, which is still a very large volumetric flow. In this case the velocity head h_{u}≃0.15 cm is similar to the hydraulic head h, and the head loss h_{f}≃0.06 cm is about h/2. Given the large values of R required to reduce the velocity head and pipe resistance head, we may instead think of the vertical tube or pipe as a wide cylindric container, with open bottom and surface ends, filled with dense deep water, the upper end located above the ocean free surface to avoid entrance of light water. The energy delivered by the wind turbine is used to withdraw the heavy water from this container. In the next section we calculate the horizontal pressure jump across the tube in order to have an insight into the resistance requirements of the tube’s materials.
HORIZONTAL PRESSURE JUMPTop
The horizontal pressure jump ∆p(z) across the tube is
$\Delta p\left(z\right)=\{\begin{array}{cc}{p}_{\text{e}}(z),& z\ge {h}_{\text{a}};\\ {p}_{\text{e}}(z)+{p}_{\text{i}}(z),& z<{h}_{\text{a}},\end{array}$

(16) 
where p_{e}(z) and p_{i}(z) are the hydrostatic pressure of the water column at depth z outside and inside the tube, respectively. In the case of the arctangent density profile, which is the most realistic profile, ∆p(z) is computed from the integration of (6) from z to z=0 for the tube’s exterior, and from z to z=h for the tube’s interior, giving
$\begin{array}{c}{\displaystyle {p}_{\text{e}}(z)=g\frac{{c}_{2}}{{c}_{4}}\left[{c}_{3}\mathrm{arctan}({c}_{3})\frac{1}{2}\mathrm{ln}\left({c}_{3}^{2}1\right)\right]g{c}_{1}z}\\ {\displaystyle g\frac{{c}_{2}}{{c}_{4}}({c}_{3}+{c}_{4}z)\mathrm{arctan}({c}_{3}+{c}_{4}z)+}\\ {\displaystyle +g\frac{{c}_{2}}{2{c}_{4}}\mathrm{ln}\left[{({c}_{3}+{c}_{4}z)}^{2}1\right]}\end{array}$

(17) 
and
p_{i}(z) = gρ_{b}(h_{a} – z). 
(18) 
The resulting pressure jump ∆p(z) is shown in Figure 5. At the ocean surface ∆p(0)=0 since both pressures are equal to the atmospheric pressure. At z≥h_{a}, ∆p(z) equals the (minus) weight (per unit area) of the exterior water column and therefore decreases almost linearly (density differences are very small) along the first h_{a}=18 cm. At h_{a}=−18 cm we have ∆p(h_{a})≃18×10^{2}P_{a}=0.18 dbar. One dbar approximately corresponds to the hydrostatic pressure of a 1mdepth water column. At depths z<h_{a} the pressure jump increases due to the larger weight per unit depth of the interior water column, and reaches ∆p=0 at z=H, which was the boundary condition chosen to compute h_{a}. Note that ∆p increases rapidly from z=h_{a} to z≃−50 m since is in this layer where density ρ_{a}(z) increases fast with depth (Fig. 2). Though this pressure seems to be small, it may be enough to deform flexible material tubes made of very thin (0.5 mm thickness) polyvinyl chloride or polyethylene, such as those used in some upwelling experiments. Thus, some periodic reinforcement of a flexible tube like that done in Maruyama et al. (2011)Maruyama S., Yabuki T., Sato T., et al. 2011. Evidences of increasing primary production in the ocean by Stommel’s perpetual salt fountain. DeepSea Res. I 58: 567574. seems to be necessary. Deployment of a flexible tube in the open ocean, assuring that is initially full of water, is therefore essential for the correct operation of the device. Deformation, including torsion, in the tube caused by vertical shear currents seems to be a more serious problem to the operation of flexible tubes. The structure supporting the pipe will presumably also have an important effect on the pipe. Strain measurements taken in the riser pipe of TAKUMI (Maeda et al. 2007Maeda K., Takahashi I., Miyabe H., et al. 2007. Strain measurement on riser pipe of TAKUMI. Proceedings of the Seventh ISOPE Ocean Mining Symposium, Lisbon, Portugal, The International Society of Offshore and Polar Engineers.) indicated that the main factor affecting the riser pipe behaviour was not the direct effect of wave and current fields on the pipe, but their indirect effects through the motion of the floating structure.
For a rigid vessel such as that proposed in this study, the circumferential hoop stress σ_{θ} can be obtained using the Barlow thin wall method,
$\sigma}_{\theta}=\frac{\Delta p*R}{t$,

(19) 
where above ∆p* is the pressure difference due to the hydraulic head and head loss, R is the pipe radius, and t is the wall thickness. Thus, an order of magnitude of the maximum permissible radiustowall thickness (R/t )_{max} required for a stainless steel pipe (maximum allowable stess σ_{θ} max of the order of 10^{2} MPa) subject to a pressure difference ∆p*≃2×∆p≈40×10^{2}Pa, is (R/t)_{max}∼σ_{θmax}/∆p*∼25×10^{3}. For a circular pipe of radius R=10 m, this implies a minimum wall thickness t_{min}∼4 mm. Furthermore, horizontal currents will cause a shear force on the tube, which will in turn tend to bend it. Weighting the bottom may help to reduce this problem. Thus, as mentioned before, deformation, including torsion, in the tube by bending moments caused by vertical shear currents seems to be more serious a problem for the operation of the container than the circumferential hoop stress caused by the hydraulic head and head loss.
Neutral height of the sinking plume
As we have seen in the previous sections, the surface water jump h is only of a few centimetres. The power generated by offshore wind turbines is of the order of several MW. If the complete wind power, say P_{w}, is used to uplift water from z=h to the surface z=0, that is, any other energy consuming processes is neglected, the volume flux F_{v} obtained is
$F}_{\text{v}}=\frac{{P}_{\text{w}}}{g{\rho}_{\text{b}}h$.

(20) 
For typical values P_{w}=2 MW, ρ_{b}=1029 kg m^{−}^{3} and h=0.16 m, we have a water volume transport estimate F_{v}≃10^{3} m^{3} s^{−}^{1}. Obviously this is a huge volume flux. Such a volume of water per unit of time with density ρ_{b} brought to the sea surface, where density ratio ρ_{b}/ρ_{t}≃1.003>1, will have no time for a complete mixing, and will inevitably lead to a fast sinking of the dense deep water developing a downwelling plume, leaving little time for primary production to develop. Due to the horizontal turbulent entrainment, this vertical plume will later reach a neutral height H_{N} beyond which no further sinking occurs. Following the classical results of (Morton et al. 1956Morton B.R., Taylor G.I., Turner J.S. 1956. Turbulent gravitational convection from maintained and instantaneous sources. Proc. R. Soc. Lond. 234: 123.) for point buoyancy sources, an order of magnitude of H_{N} is given by
$H}_{\text{N}}=2.1\sqrt[4]{\frac{\pi {f}_{0}}{{N}^{3}}$
, 
(21) 
where f_{0} is the surface buoyancy flux and N is the environment BruntVäisälä or buoyancy frequency,
$N}^{2}\equiv \frac{g}{{\rho}_{0}}\frac{d{\rho}_{\text{B}}}{dz$
, 
(22) 
where ρ_{0} is a constant reference density and ρ_{B}(z) is the environment (background) density profile. In our case, f_{0} is related to the volume flux F_{v} by
$f}_{\text{0}}={g}^{\prime}\frac{{F}_{\text{v}}}{\pi$
, 
(23) 
where gʹ is the reduced gravity
$g}^{\prime}\equiv g\frac{{\rho}_{\text{b}}{\rho}_{\text{t}}}{{\rho}_{\text{0}}$.

(24) 
Finally, the neutral height H_{N} is related to the volume flux F_{v} through
$H}_{\text{N}}=2.1\sqrt[4]{\frac{{g}^{\prime}{F}_{\text{v}}}{{N}^{3}}$
.

(25) 
For the values gʹ=10(1029−1026)/1028 m s^{−}^{2}, F_{v}=10^{3} m^{3} s^{−}^{1}, and N=1/50 s^{−}^{1}, we obtain a neutral height H_{N}≃90 m, or H_{N}≃77 m for F_{v}=500 m^{3} s^{−}^{1}. Thus, due to the very large volume flux F_{v} the nutrientrich deep water will probably sink to the bottom limit of the photic layer, and primary production will not be considerably enhanced.
To avoid this large vertical sinking, the upwelled water must be dispersed in some way at the surface to increase turbulent entrainment. The point source theory for vertical plumes used above must therefore be modified to include finite area sources. A number of corrections have been developed (Hunt and Kaye 2001Hunt G.R., Kaye N.G. 2001. Virtual origin correction for lazy turbulent plumes. J. Fluid Mech. 435: 377396.) to take into account finite area sources in the theory of turbulent entrainment. A simple estimate may be given by noticing that the angle θ of the plume’s cone (Fig. 6) is related to the diameter D of the source and to the virtual height z_{v} of the point source by
$\mathrm{tan}\frac{\theta}{2}=\frac{D}{2{z}_{\text{v}}}$.

(26) 
Experimental results show that the total angle θ enclosed is approximately 18°20°. If the mixed upwelled water is to remain close to the surface, we must impose z_{v}≃H_{N}. Thus, the upwelled water must be dispersed to the surface over a circle of radius r∼H_{N} tan(θ/2)≃15 m. This water dispersion could be achieved by a number of mechanical, nonpowerconsuming ways. For example, a horizontal perforated disc platform of 15 m radius located at the surface but below the deep water outflow could work well as a diffuser plate, or multiport diffusers such as those used in ocean outfall dilution could be installed. Note also that this radius r≃15 m has an order of magnitude similar to the pipe radius R required to keep power losses due to friction below reasonable low values. However, if oxygen enrichment, or aeration, of the deep water is desirable, a portion of the turbine’s energy could be used for water spreading (a fountain aerator) and only the remaining portion will be consumed for water raising. A gross estimation of the order of magnitude of the power required for water sprinkling is given in the next section.
POWER BUDGET CONSIDERING POWER CONSUMPTION REQUIRED FOR WATER SPRINKLINGTop
In this case the wind power P_{w} must be employed both for water upwelling (P_{u}) and water spreading (P_{s}). The larger the water volume transport F_{v}, the larger the amount of energy required for water spreading P_{s }and the smaller the amount of energy for water upwelling P_{u}, which means the smaller the amount of water transport F_{v} achieved. An estimate of the power P_{s }(R, F_{v}) consumed for homogeneously sprinkling, over a circle of radius R, a volume transport of water F_{v} leaving a point source located at the circle’s centre is given in Appendix D. The basic result is that P_{s}(R, F_{v})=(ρg/3)RF_{v}. We stress that the above relation is only a rough estimate of the minimum bound for the ideal power required to spread a volume flux of water F_{v} over a radius R. The practical power P_{s }(R, F_{v}) might be an order of magnitude larger. In order to take this uncertainty into account, we include a factor κ_{0} in P_{s}, that is,
$P}_{\text{s}}\left(R,{F}_{\text{v}}\right)={\kappa}_{0}\frac{{\rho}_{\text{b}}g}{3}R{F}_{\text{v}$
, 
(27) 
and assume that typically κ_{0}∼10. Thus, the variables P_{u}, P_{s }, F_{v} and R are coupled through the following system of equations
P_{w} = P_{u} + P_{s} , 
(28) 
$F}_{\text{v}}=\frac{{P}_{\text{u}}}{g{\rho}_{\text{b}}h$,

(29) 
$R=2.1\sqrt[4]{\frac{{g}^{\prime}{F}_{\text{v}}}{{N}^{3}}}\mathrm{tan}\theta$
, 
(30) 
$P}_{\text{s}}={\kappa}_{0}\frac{{\rho}_{\text{b}}g}{3}R{F}_{\text{v}$
, 
(31) 
where the unknowns are (P_{u}, P_{s}, F_{v}, R). By simple substitution we arrive at a single algebraic equation for F_{v},
$F}_{\text{v}}\left({\kappa}_{0}\frac{2.1}{3}\sqrt[4]{\frac{{g}^{\prime}}{{N}^{3}}}\mathrm{tan}\theta {F}_{\text{v}}^{1/4}+h\right)=\frac{{P}_{\text{w}}}{{\rho}_{\text{b}}g$
, 
(32) 
which can be solved numerically noticing that the second term on the left hand side of (32), F_{v}h, is small compared to the first one. For typical values P_{w}=2 MW, gʹ=0.028 m s^{−}^{2}, N=1/50 s^{−}^{1}, and θ=20°, we obtain F_{v}=6.3 m^{3} s^{−}^{1}, R≃9 m, P_{s}=1.99 MW, and P_{u}=0.01 MW. This solution implies that 99.4% of the wind power will be employed for water sprinkling while only the remaining 0.6% will be used for raising water at a very reduced volume flux of about 6 m^{3} s^{−}^{1}. These results show that water sprinkling, as a way to favour both the turbulent entrainment and oxygen enrichment of the deep water, is very power consuming and, as stated above, use of nonpowerconsuming diffusion devices such as multiport diffusers or a perforated disc platform is more advantageous.
NUTRIENTS TRANSPORTTop
Once we have an estimation of the order of magnitude of the upwelling volume transport (say F_{v}∼500 m^{3} s^{−}^{1}), we can also provide an estimation for the transport of nutrients to the photic layer. Time and space averages of vertical profiles of nitrates, phosphates and silicates, characteristics of the Catalan Sea, are given in Figure 7. As a first approximation we can assume that at a depth z=−300 m the mass concentration of nitrates, phosphates and silicates is, respectively, C_{N}=7 µ mol l^{–1}, C_{P}=0.4 µ mol l^{–1} and C_{S}=5 µ mol l^{–1}. The transport of nutrients to the photic layer is therefore F_{N}= C_{N}F_{v} ∼3.5 mol s^{−}^{1}, F_{P}=C_{P}F_{v}∼0.2 mol s^{−}^{1} and F_{S}=C_{S}F_{v}∼2.5 mol s^{−}^{1}. This means mass transports T_{N}∼217 g s^{−}^{1}=19 t day^{−}^{1}, T_{P}∼19 g s^{−}^{1}=1.6 t day^{−}^{1} and T_{S}∼150 g s^{−}^{1}=13 t day^{–1}, where we have used the molecular masses M_{N}≃62 g mol^{–1}, M_{P}≃95 g mol^{–1} and M_{S}≃60 g mol^{–1}. Typical values of offshore water transport caused by wind stress in natural upwelling coastal areas are of the order of 10^{3} m^{3} s^{−}^{1} per km of coastline (upwelling index U=10^{3} m^{3} s^{−}^{1} km^{−}^{1}). Thus, an upwelling F_{v}∼500 m^{3} s^{−}^{1} produced by a wind turbine is roughly equivalent to the natural upwelling along a 500 m coastline.
We notice that large artificial upwelling in an ocean region with poor biological productivity can change the trophic cycle positively or negatively. The sudden growth of primary producers can generate the expected positive response, i.e. a higher productivity with the consequent increase of economically important marine resources. It can also generate a negative response, favouring proliferation of damaging or simply nonautochthonous species that will put an end to the flora and fauna of the region. Short, medium, and large time scale responses of the ecosystem to large artificial upwelling are therefore a very important issue that has been addressed little (McAndrew et al. 2007McAndrew P.M., Björkman K.M., Church M.J., et al. 2007. Metabolic response of oligotrophic plankton communities to deep water nutrient enrichment. Mar. Ecol. Prog. Ser. 332: 6375., Williamson et al. 2009Williamson N., Komiya A., Maruyama S., et al. 2009. Nutrient transport from an artificial upwelling of deep sea water. J. Oceanogr. 65: 349359.) and needs to be seriously taken into account. The relation between edible fish production and upward flow rate of deep nutrientrich water also needs to be studied (Groves 1958Groves G.W. 1958. Flow estimate for the perpetual salt fountain. Deep Sea Res. 5: 209214., Thorson 1961Thorson J. 1961. Flow estimate for an artificial upwelling in the sea. Deep Sea Res. (1953) 7(4): 296297., Matsuda et al. 1999Matsuda F., Szyper J., Takahashi P., et al. 1999. The ultimate ocean ranch. Sea Technol. 40(8): 1726.). The ecosystem response is a problem beyond the scope of this work, which only develops, in terms of power requirements, the idea of artificial upwelling using wind power. Artificial upwelling using wind power is, however adjustable, and can be easily reduced, or turned off, if there is any risk of adverse effects to the ecosystem.
DISCUSSION AND CONCLUDING REMARKSTop
In this work we have put forward the idea of taking advantage of in situ offshore wind energy to cause upwelling of nutrientrich deep water to favour primary production in the ocean. With a wind turbine power of several MW, the order of magnitude of the upward water transport is several hundreds of m^{3} s^{−}^{1}, a transport much larger than that achieved using surface wave energy or the perpetual salt fountain mechanism. Here we propose enlarging the spar of a sparbuoy wind turbine to channel the deep water up to the sea surface. This long rigid tube also acts as a ballast, providing additional stability to sparbuoy offshore wind turbines. This work is, however, only a first approximation to the general problem. We basically have found that the upward volume transport could be so large that turbulent water entrainment must be enhanced in some way to avoid sinking of the dense water to a neutral depth below the photic layer. An energetically cheap way to dilute the dense upwelled water is therefore required.
There are still a large number of basic questions that need to be properly addressed. One of them is the ecological response. Though a water rich in nutrients is in theory beneficial to primary production, a large and concentrated source of nutrients may be harmful in an oceanic region where the ecosystem is not yet well adapted to a large availability of nutrients. For example, at the Mississippi mouth excess nutrients carried down the river in floods cause massive algal blooms. A slow startup of the fertilization could therefore be desirable. Since it is planned to obtain offshore wind energy in wind farms comprising a large number of wind turbines, it could be possible to use only a part of the electric power of several wind turbines to raise deep water. This will reduce the upward volume transport assigned to each wind turbine, therefore favouring turbulent entrainment, though it would require an expensive 300m sparbuoy on each wind turbine. Thus, the benefits of artificial upwelling might largely counteract minor potential negative environmental impacts of offshore wind farms (Wilson et al. 2010Wilson J.C., Elliott M., Cutts N.D., et al. 2010. Coastal and offshore wind energy generation: Is it environmentally benign? Energies 3(7): 13831422.).
Another important problem concerns engineering issues. Nutrientrich water below the photic layer is located in many ocean regions at a minimum depth of about 300 m. This implies the design, manufacturing, testing, and deployment of vertical spars almost twice as large as the ones currently designed for sparbuoy offshore wind turbines (EWEA 2013EWEA. 2013. Deep water  The next step for offshore wind energy. Tech. rep. European Wind Energy Association.). Also, since energy losses due to friction in the pipe become large for a reasonable pipe radius, a compromise between wall thickness, pipe radius and upward volume transport must be reached through very careful engineering analysis. Fatigue, stress and load analysis of these new structures need to be carried out.
Some physical oceanographic questions also need to be addressed with great care. Pointlike localized upwelling of large water volume transport cannot be maintained for long periods of time unless ocean currents advect the surface mixed water far from the wind turbines and replace it with new, lighter water. Otherwise the complete water column will homogenize into dense water and deep sinking of the upwelled water will occur despite turbulent entrainment. However, while some amount of ocean currents is beneficial to maintain the nutrients in the photic layer, large amplitude currents, and in particular the large vertical shear associated with baroclinic currents, might be detrimental to the stability of the long tube and the turbine’s tower. It is therefore important for mariculture applications to select the offshore wind farm location taking into account the wind, wave and ocean current climatology.
The influence of surface waves and shear currents favouring turbulent entrainment needs to be properly investigated. Also, in the absence of significant background ocean shear currents, the upwelled water mixed with the surface water could form large bowls of dense water which, on a time scale of a few days, will start rotating cyclonically, due to the centripetal and Coriolis accelerations, and will approach the cyclogeostrophic balance that is characteristic of mesoscale and submesoscale dynamics. An order of magnitude of the time required to generate a cyclogeostrophic cyclone from a source of deep water can be roughly estimated assuming that in a completely mixed and cyclogeostrophic balanced state the cyclone develops a Gaussian radial (negative) free surface height anomaly h(r)=h_{0} exp(−r/r_{0}), where r is the radial distance, h_{0}<0 is the height at the vortex centre, and r_{0}≃20 km is the horizontal length scale of the vortex. The volume integral of h(r) between r=0 and r=∞ is the volume anomaly Vʹ=πh_{0}r_{0}^{2}. The mass anomaly of this vortex is therefore Mʹ=−ρ_{0}Vʹ>0. A constant flux of mass anomaly Fʹ_{m}=ρʹF_{v} will take a time t=Mʹ/Fʹ_{m}= ρ_{0}Vʹ/(ρʹF_{v}) to complete the vortex mass anomaly Mʹ. For ρ_{0}=10^{3} kg m^{−}^{3}, ρʹ=3 kg m^{−}^{3}, h_{0}=1 cm, r_{0}=20 km and F_{v}=500 m^{3} s^{−}^{1}, we obtain an order of magnitude of the time period t∼10^{7} s∼10^{2} days. This is still too long time a period to neglect the effect of background ocean currents, but if exceptionally quiescent conditions are found, or if several wind turbines are employed, the high volume transports achieved open up the possibility of artificial submesoscale eddy formation. Mesoscale and submesoscale eddies are very coherent stable structures in which mixing with surrounding waters is inhibited. When these vortices leave the upwelling site, they retain the nutrients and other sea water properties for a long time period. In the extreme case of very quiescent waters, an anticyclone may be artificially generated by downwelling light water using a nearby wind turbine. If a cycloneanticyclone pair is generated close enough, the two vortices will form a dipole vortex which will start translating coherently along the dipole axis. That scenario seems quite simplistic because the upwelling produced by a turbine only occurs when wind is present, so wind stress on the surface ocean layer will likely provide, via turbulent momentum diffusion, enough horizontal momentum to the upper layer waters to avoid generation of shallow eddies, or at least to strongly modify its surface structure. However, we mention this artificial dipole formation because the large vertical transport of water that wind turbines might cause opens up the opportunity to modify the regional submesoscale ocean circulation. On the other hand, large wind farms may produce a significant disturbance on the surface wind stress, generating upwelling/downwelling velocities in the wind wake that might affect the local ecosystem (Broström 2008Broström G. 2008. On the influence of large wind farms on the upper circulation. J. Mar. Syst. 74: 585591.).
As an overall conclusion we might say that, from the energetic point of view, artificial upwelling using offshore wind energy seems to be a promising way to enhance primary production in the ocean. Mariculture application of this approach implies the fertilization of large regions in the open ocean, and is therefore severely subjected to atmosphere and ocean climatology, as well as to ecological dynamics. The political issues also need to be addressed: if productivity is enhanced in the open ocean, costs and benefits must be properly distributed among the interested countries. The general problem is a multidisciplinary one, and we have noticed that some important physical, engineering, and ecological questions need to be seriously addressed to obtain a more complete confidence in the approach presented here.
ACKNOWLEDGEMENTSTop
We thank two anonymous reviewers for their very useful comments. Partial support for this study was obtained through projects CTM201128867 and CTM201456987P (Spanish Ministry of Science and Innovation).
REFERENCESTop
Aure J., Strand Ø., Erga S.R., et al. 2007. Primary production enhancement by artificial upwelling in a western Norwegian fjord. Mar. Ecol. Prog. Ser. 352: 3952.
http://dx.doi.org/10.3354/meps07139
Berntsen J., Aksnes D.L., Foldvik A. 2002. Production enhancement by artificial upwelling: a simulation study. Hydrobiology 484: 177190.
http://dx.doi.org/10.1023/A:1021317424276
Broström G. 2008. On the influence of large wind farms on the upper circulation. J. Mar. Syst. 74: 585591.
http://dx.doi.org/10.1016/j.jmarsys.2008.05.001
CostaPierce B.A. 2002. Ecological Aquaculture: The Evolution of the Blue Revolution. Blackwell Science, 382 pp.
Davis J. 1990. Wave energy absorption by the Bristol Cylinderlinear and nonlinear effects. Proc. Inst. Civil Eng. 89: 317340.
http://dx.doi.org/10.1680/iicep.1990.9394
de Baar H.J.W., Gerringa L.J.A., Laan P., et al. 2008. Efficiency of carbon removal per added iron in ocean iron fertilization. Mar. Ecol. Prog. Ser. 364: 269282.
http://dx.doi.org/10.3354/meps07548
EWEA. 2013. Deep water  The next step for offshore wind energy. Tech. rep. European Wind Energy Association.
EWEA. 2015. The European offshore wind industry  key trends and statistics 1st half 2015. Tech. rep. European Wind Energy Association.
FAO. 2010. The state of the world fisheries and aquaculture. Rome: Food and Agriculture Organization of the United Nations.
Groves G.W. 1958. Flow estimate for the perpetual salt fountain. Deep Sea Res. 5: 209214.
http://dx.doi.org/10.1016/01466313(58)900133
Henderson A.R., Morgan C., Smith B., et al. 2003. Offshore wind energy in Europe – a review of the stateoftheart. Wind Energ. 6: 3552.
http://dx.doi.org/10.1002/we.82
Hunt G.R., Kaye N.G. 2001. Virtual origin correction for lazy turbulent plumes. J. Fluid Mech. 435: 377396.
http://dx.doi.org/10.1017/S0022112001003871
Isaacs J.D., Castel D., Wick G.L., 1976. Utilization of the energy in ocean waves. Ocean Eng. 3: 175187.
http://dx.doi.org/10.1016/00298018(76)900226
Jonkman J.M. 2009. Dynamics of offshore floating wind turbines – model development and verification. Wind Energ. 12: 459492.
http://dx.doi.org/10.1002/we.347
Jonkman J.M., Matha D. 2011. Dynamics of offshore floating wind turbines – analysis of three concepts. Wind Energ. 14: 557569.
http://dx.doi.org/10.1002/we.442
Kirke B. 2003. Enhancing fish stocks with wavepowered artificial upwelling. Ocean Coast. Manag. 46: 901–915.
http://dx.doi.org/10.1016/S09645691(03)00067X
Kithil P.W. 2011. Comments on “An open ocean trial of controlled upwelling using wave pump technology”. J. Atmos. Ocea. Tech., 28: 847849.
http://dx.doi.org/10.1175/2010JTECHO778.1
Lavagnini A., Sempreviva A.M., Transerici C., et al. 2006. Offshore wind climatology over the Mediterranean Basin. Wind Energ. 9: 251266.
http://dx.doi.org/10.1002/we.169
Liu C. 1999. Research on artificial upwelling and mixing at the University of Hawaii at Manoa. IOA Newslett. 10: 18.
Maeda K., Takahashi I., Miyabe H., et al. 2007. Strain measurement on riser pipe of TAKUMI. Proceedings of the Seventh ISOPE Ocean Mining Symposium, Lisbon, Portugal, The International Society of Offshore and Polar Engineers.
Maruyama S., Tsubaki K., Taira K., et al. 2004. Artificial upwelling of deep seawater using the perpetual salt fountain for cultivation of ocean desert. J. Oceanogr. 60: 563568.
http://dx.doi.org/10.1023/B:JOCE.0000038349.56399.09
Maruyama S., Yabuki T., Sato T., et al. 2011. Evidences of increasing primary production in the ocean by Stommel’s perpetual salt fountain. DeepSea Res. I 58: 567574.
http://dx.doi.org/10.1016/j.dsr.2011.02.012
Matsuda F., Szyper J., Takahashi P., et al. 1999. The ultimate ocean ranch. Sea Technol. 40(8): 1726.
McAndrew P.M., Björkman K.M., Church M.J., et al. 2007. Metabolic response of oligotrophic plankton communities to deep water nutrient enrichment. Mar. Ecol. Prog. Ser. 332: 6375.
http://dx.doi.org/10.3354/meps332063
McClimans T.A., Handå H., Fredheim A., et al., 2010. Controlled artificial upwelling in a fjord to stimulate nontoxic algae. Aquac. Eng. 42: 140147.
http://dx.doi.org/10.1016/j.aquaeng.2010.02.002
McKinley K.R., Takahashi P.K. 1991. Deep ocean water, artificial upwelling, and open ocean mariculture: A promise for the future. IEEE 1: 195199.
http://dx.doi.org/10.1109/OCEANS.1991.613930
Morton B.R., Taylor G.I., Turner J.S. 1956. Turbulent gravitational convection from maintained and instantaneous sources. Proc. R. Soc. Lond. 234: 123.
http://dx.doi.org/10.1098/rspa.1956.0011
Muliawan M.J., Karimirad M., Moan T. 2013. Dynamic response and power performance of a combined Spartype floating wind turbine and coaxial floating wave energy converter. Renew. Energy 50: 4757.
http://dx.doi.org/10.1016/j.renene.2012.05.025
Otsuka K. 2011. Inclusive Impact Index Triple I and its application for Ocean Nutrient Enhancer. Proceedings of the Twentyfirst International Offshore and Polar Engineering Conference, International Society of Offshore and Polar Engineers.
Ouchi K. 2007. Three years operation of ocean nutrient enhancer TAKUMI in Sagami Bay. Seventh ISOPE Ocean Mining Symposium, 16 July, Lisbon, Portugal, International Society of Offshore and Polar Engineers.
Ouchi K., Yamatogi T., Kobayashi K., et al. 1998. Density Current Generator—A new concept machine for agitating and upwelling a stratified water. Ocean Community Conference, 1, Baltimore, MD, USA, The Marine Technology Society, 129136.
Pauly D., Froese R. 2012. Comments on FAO’s State of Fisheries and Aquaculture, or ‘SOFIA 2010’. Mar. Pol. 36: 746752.
http://dx.doi.org/10.1016/j.marpol.2011.10.021
Roddier D., Weinstein J. 2010. Floating wind turbines. Mech. Eng. 132: 2832.
Salter S.H. 2009. A 20 GW thermal 300metre^{3}/sec waveenergised, surgemode nutrient pump for removing atmospheric carbon dioxide, increasing fish stocks and suppressing hurricanes. Proceedings of the 8th European Wave and Tidal Energy Conference, Uppsala, Sweden.
Stommel H., Arons A.B., Blanchard D. 1956. An oceanographical curiosity: the perpetual salt fountain. DeepSea Res. 3: 152153.
http://dx.doi.org/10.1016/01466313(56)900958
Thorson J. 1961. Flow estimate for an artificial upwelling in the sea. Deep Sea Res. (1953) 7(4): 296297.
http://dx.doi.org/10.1016/01466313(61)900508
Tsubaki K., Maruyama S., Komiya A., et al. 2007. Continuous measurement of an artificial upwelling of deep sea water induced by the perpetual salt fountain. DeepSea Res. I 54: 7584.
http://dx.doi.org/10.1016/j.dsr.2006.10.002
UNESCO 2010. Ocean fertilization. A scientific summary for policy makers. 17 pp.
Viré A. 2012. How to float a wind turbine. Rev. Environ. Sci. Biotechnol. 11: 223226.
http://dx.doi.org/10.1007/s1115701292929
White A., Björkman K., Grabowski E., et al. 2010. An open ocean trial of controlled upwelling using wave pump technology. J. Atmos. Ocean. Tech. 27: 385396.
http://dx.doi.org/10.1175/2009JTECHO679.1
White A.E., Letelier R.M., Björkman K.M., et al. 2011. Reply. J. Atmos. Ocean. Tech. 28: 850851.
http://dx.doi.org/10.1175/2010JTECHO796.1
Williamson N., Komiya A., Maruyama S., et al. 2009. Nutrient transport from an artificial upwelling of deep sea water. J. Oceanogr. 65: 349359.
http://dx.doi.org/10.1007/s108720090032x
Wilson J.C., Elliott M., Cutts N.D., et al. 2010. Coastal and offshore wind energy generation: Is it environmentally benign? Energies 3(7): 13831422.
http://dx.doi.org/10.3390/en3071383
Zaaijer M.B. 2009. Review of knowledge development for the design of offshore wind energy technology. Wind Energ. 12: 411430.
http://dx.doi.org/10.1002/we.349