Waves and Swells in High Wind and Extreme Fetches, Measurements in the Southern Ocean

The generation and evolution of ocean waves by wind is one of the most complex phenomena in geophysics, and is of great practical significance. Predictive capabilities of respective wave models, however, are impaired by lack of field in situ observations, particularly in extreme Metocean conditions. The paper outlines and highlights important gaps in understanding the Metocean processes and suggests a major observational program in the Southern Ocean. This large, but poorly investigated part of the World Ocean is home to extreme weather around the year. The observational network would include distributed system of buoys (drifting and stationary) and autonomous surface vehicles (ASV), intended for measurements of waves and air-sea fluxes in the Southern Ocean. It would help to resolve the issues of limiting fetches, extreme Extra-Tropical cyclones, swell propagation and attenuation, wave-current interactions, and address the topics of wave-induced dispersal of floating objects, wave-ice interactions in the Marginal Ice Zone, Metocean climatology and its connection with the global climate.


INTRODUCTION
The generation and evolution of ocean waves by wind is one of the most complex phenomena in geophysics. Forecasting skill and understanding of these dynamics is critical across a wide range of oceanic applications, including maritime and coastal engineering, air-sea interactions, ocean dynamics, climate, remote sensing. However, the generation and evolution of waves in high-wind conditions and extreme fetches remains poorly understood.
Models are widely inconsistent for large fetch/duration conditions, to a significant extent due to lack of observational guidance. Extreme storms in the North Atlantic and Pacific are seasonal and diverse in their propagation which hinders systematic in situ observations in these regions (e.g., Meucci et al., 2018;Takbash et al., 2018). Also, they usually do not provide extreme fetches. Observations of high-wind conditions with extreme fetch, however, are possible in the Southern Ocean, where these conditions occur regularly and storms move in the same direction (West to East) throughout the year. Swell, which results from such storms in the Southern Ocean radiates into all major ocean basins, but remains poorly predicted by forecast models, both in magnitude and arrival time.
This paper proposes systematic in situ field observations in the Southern Ocean. Understanding the extreme fetch and forcing conditions, and respective advance of the wave models is possible through such observations conducted by means of deployment of a distributed buoy network (drifting and stationary) and autonomous surface vehicles (ASV) in the Southern Ocean. The buoys and ASVs could be deployed from ships of opportunity and through designated efforts of interested countries (e.g., Schulz et al., 2012).
Apart from the extreme fetch and swell problems, such an observational network in the Southern Ocean would address a number of Metocean topics that remain unresolved for decades. The Sections below cover the following unattended problems of Metocean conditions in the World Ocean: • Wave evolution at extreme fetches; • Severe extra-tropical cyclones at extreme fetches; • Swell dynamics and forecasting, with attention to arrival time; • Wave-driven dispersal of floating objects (search and rescue, transport of microplastic and other pollutants); • Non-linear wave-current interactions; • Wave-ice interactions in Marginal Ice Zone subject to extreme storms and waves; • Metocean climatology in the Southern Ocean.

WAVE EVOLUTION AT EXTREME FETCHES
Following the classical paper by Pierson and Moskowitz (1964), it is commonly accepted that there is a limiting condition for wave development such that, for a given wind speed, the significant wave height H s and peak wavelengths (periods) stop growing. In non-dimensional terms of mean wind speed at standard 10 m height U 10 and phase speed of peak waves c p , the limiting stage of wave development is described by ratio U 10 /c p = 0. 82 (1) which, basically, means that once the dominant waves in a windgenerated field (spectrum) are faster than the wind, the wave development ceases. While intuitively attractive, this concept does not necessarily agree with observations -for example, Young (2006) for Tropical Cyclones and Thomson and Rogers (2014) for lighter winds, demonstrated measurements of windgenerated waves well beyond the PM limit. Thus, 50+ years after Pierson and Moskowitz, such limit is still in need of validation, clarification, understanding and explaining. While it was purely empirical concept originally, now we can speculate on such a limit from a more advanced physical perception of wind-wave evolution. Such perception includes non-linear interactions which have no regard for the wind and maintain the energy flux to low frequencies (i.e., wave periods larger than the spectral peak), provided the energy flux to the high frequencies continues (e.g., Zakharov and Zaslavskii, 1982). In principle, such behavior would signify no full development, but at some stage the very long waves would be so fast that the friction against the air (no matter how small it is), would balance the weakening non-linear energy influx [e.g., the mechanism for wind-wave interactions when waves overtake the wind in Donelan et al. (2012)].
Not surprisingly, in absence of clear physical guidance, performance of wave-forecast models in the context of full development is contradictory and far from being consistent. While typically tuned to the PM saturation in academic tests, the models hardly ever meet the limiting criteria in realistic simulations. In Figure 1 such comparisons are reproduced from Rogers (2002), for three different physics packages: ST1, ST2, and WAM4, with the first two being from the WAVEWATCH-III model (WW3, Tolman, 2002). In all the cases the mean wind speed is U 10 = 15 m/s which, if allowed to persist over unlimited fetch/duration, should lead to ultimate PM wave height of H s = 5.5 m. None of the models do -they do not even come close -and none reaches another asymptote. This is also true of newer physics packages available in recent versions of WW3 (Wavewatch III R Development Group, 2016) -ST3, ST4, and ST6 (unpublished). Models in Figure 1 are dissimilar, so at least two are wrong, and probably all three are wrong, but we note that the behavior simulated here is unvalidated, due to scarce observations. Such extreme fetch/duration is a primary "frontier" area for observations, associated with uncertainty in the models.
In principle, if quasi-full-development exists, it should be easier to reach for lower winds [albeit not in Thomson and Rogers (2014)] than for higher winds. It may never happen for high winds due to very long fetches required, but if it happens anywhere, it would be in the Southern Ocean where, depending on the speed of propagation of extra-tropical cyclones, the fetches can be virtually unlimited. The question of the full development is not hypothetical and/or of pure academic value and interest: lack of understanding of wave evolution at the extreme end of Metocean conditions and wave fetches imposes real limitation on performance of models in circumstances which are most critical for maritime engineering and operations. Therefore, a network or array of wave buoys (or a set of drifting buoys) in east-west direction in the Southern Ocean would be able to prove or disprove the concept of fully developed limiting stage. Most importantly, such concept, intuitively FIGURE 1 | Simulating the PM limit in WW3 (WAM3/ST1, WAM4, ST2 packages from left to right, respectively). U 10 = 15 m/s corresponds to PM H s = 5.5 m. Horizontal scale is duration of wave evolution in days, vertical scale is fetch in kilometers, and color scale is wave height in meters [Figure is reproduced from Rogers (2002)].
attractive and most like correct, needs quantification which is only possible on the basis of solid experimental evidence.

SEVERE EXTRA-TROPICAL CYCLONES AT EXTREME FETCHES
As a reference point for the extreme Metocean conditions, the hurricane-scale classification is often used: that is a tropical storm becomes a hurricane if the wind speed reaches U 10 ∼ 33 m/s. Babanin (2018) argued that such classification is not arbitrary, and indeed signifies change of physical regimes in all the three environments near the air-sea interface: in the atmospheric boundary layer, at the surface, and through the upper ocean. This threshold is approximately the wind speed at which the drag coefficient was found to saturate in the field observations [U 10 ≈ 32-33 m/s, e.g., in Powell et al. (2003)]. This saturation has received a lot of attention lately. Less known are the in situ measurements below the surface, change of the upper-ocean mixing mechanism and of bubble dynamics occur at U 10 > 35 m/s (McNeil and D'Asaro, 2007). Directly at the surface, wave dynamics also undergoes essential transformations, from wave breaking (dissipation) being driven by evolution of nonlinear waves, to the breaking being forced directly by the winds, at U 10 ≈ 34 m/s [Babanin (2011) based on laboratory measurements of Leikin et al. (1995)]. Perhaps related to the wave-breaking change of mechanism is the most striking and abrupt alteration of the gas (CO 2 ) transfer at U 10 = 33.6 m/s in laboratory experiments of Iwano et al. (2013). It is therefore argued that the simultaneous change of physical regime in all the three airsea environments cannot be coincidental. Such change of the regime means that if we extrapolate our parameterisations from regular conditions into the extreme Metocean environments (which is what we usually do), we will obtain biased or even incorrect results.
It is easy to appreciate the significance of understanding and adequate modeling of waves in such conditions, both for practical and academic purposes, and the associated difficulties which to a large extent are due to lack of respective measurements. What is not appreciated, perhaps due to the lack of observations, is how different are the evolution of such waves in extreme Tropical (TC) and Extra-Tropical (ETC) cyclones. While the waves with H s in excess of 15 m are not uncommon in both cases [e.g., Young (2006) for TC and Rapizo et al. (2015) for ETC] their directional spectra are very different. Young (2006) based on a large collection of directional spectra in tropical cyclones demonstrated that direction of peak waves does not follow the local wind and, at some quadrants of TCs can be at 90 and even 180 degrees to the wind, whereas in ETCs Rapizo et al. (2015) did not observe unexpected major deviations between wind and wave propagation angles. This means that, while wave evolution in ETC probably follows the direct wind-forcing pattern, in TCs this evolution is different. Young (2006) argued that the presence of large waves propagating perpendicular or even against the wind can only be explained due to the fact that their growth is controlled by non-linear interactions. If so, this is not just an academic curiosity: both the wave-growth dependences and asymptotic behaviors of respective waves will be different (Badulin et al., 2007).
Needless to say, that such differences are not validated and not even accounted for in the current wave forecast. While measurements in Tropical Cyclones are rare, the detailed and consistent measurements in Extreme-Tropical Cyclones are nearly absent. Dedicated effort in the Southern Ocean, where ETCs are continuously present around the year, would be the best observational ground for such extreme Metocean circumstances. Presently, there is only one flux station available at 47 degrees south, 142 degrees East [south of Tasmania Schulz et al. (2012)], and it is proposed that such stations, or air-sea interaction buoys, are deployed south of New Zealand, South America and South Africa. Autonomous surface vehicles (ASV), deployed from ships of opportunity or as a dedicated effort have also proved an efficient way of investigating the air-sea interactions in extreme weather (Schmidt et al., 2017;Thomson and Girton, 2017).

SWELL DYNAMICS AND FORECASTING, WITH ATTENTION TO ARRIVAL TIME
Swell waves are present in most of ocean wave spectra (e.g., Semedo et al., 2011), and provide significant adverse impact on maritime operations and coastal inundation. Their prediction by wave-forecast models, however, is poor, both in terms of wave amplitude and, particularly, arrival time.
The third-generation models, until recently, have entirely based their physics on dynamics and interactions of windgenerated seas. In phenomenological terms, such models simulate the Radiative Transfer Equation [see, e.g., the state-ofthe-art review by Cavaleri et al. (2007)]: where E is wave spectrum, which changes in space and time and whose integral is the total wave energy, and the right-hand side are sources S in (from the wind), sinks S ds (usually due to wave breaking) and redistribution terms S nl for this energy (more terms are available in specific circumstances, particularly in finite depths). While the forecast based on (2) is applied globally, none of the terms on the right, strictly speaking, applies to the swell: swell is not wind-forced (by definition), swell does not break in deep water because of its low steepness, and the Hasselmann resonant interactions usually employed as S nl are not applicable to swells because they are unidirectional and therefore cannot satisfy the resonance conditions. The very definition of ocean swell is ambiguous: while it is usually perceived as former wind-generated waves, in fact it may reconnect with the local wind through non-linear interactions. The visible swell attenuation is driven by a number of dissipative and non-dissipative processes. The dissipative phenomena include interaction with turbulence on the water and air sides (e.g., Babanin, 2006;Ardhuin et al., 2010), with adverse winds or currents (e.g., Donelan, 1999;Babanin et al., 2017, respectively). Non-dissipative contributions to the gradual decline of wave amplitude come from frequency dispersion and directional spreading, refraction by currents, and lateral diffraction of wave energy (e.g., respectively, Ardhuin et al., 2009;Babanin and Waseda, 2015;Rapizo et al., 2018). The interactions with local winds/waves can, on the contrary, cause swell growth (perhaps some observations by Ardhuin et al. (2009)

fall into this category).
Swell arrival time is the least understood and the most uncertain problem. Joint analysis of buoy observations and model reanalysis shows that swell can be tens of hours early or late by comparison with model predictions (Jiang et al., 2016), see Figure 2. This is where the lack of model performance incurs the worst consequences: many practical applications related to swell depend not so much on swell height and steepness (which is usually low), but on its presence or absence (operating the tankers, dredging, ports).
Obviously, since the arrival-time error can be both negative and positive, no single physical mechanism can be held responsible for such failure to perform. Rather, this is a combination of various mechanisms, particularly as swells propagate very large distances over vast ocean surfaces and hence even a single swell event can be subject to multiple influences (Babanin and Jiang, 2017). Finite frequency resolution of the initial wave spectrum in a model can be a reason, randomly responsible for early/late arrival, albeit small. For the early arrival, swell has to be accelerating as it moves away from the distant source, such acceleration can be perceived (i.e., short wavelength decaying faster than the longer waves) or real. A real acceleration of waves can be caused, for example, by so-called Raman effect -downshift of wave energy due to modulational instability of non-linear waves in dispersive environments. This effect is well known in non-linear optics (e.g., Gordon, 1986), and has been perceived for the surface waves too (Segur et al., 2005). Interactions of swell with local winds, waves, currents, or a combination of those, can bring about a plethora of accelerating/decelerating effects. For example, adverse currents with horizontal velocity gradient instigate modulational instability and may lead to sudden frequency downshift and propagation acceleration . Gradual decrease of wave steepness in the course of wave attenuation should cause slight deceleration of moving swells. Refraction of waves by currents and large-scale eddies, permanent and abundant in the Southern Ocean and at the periphery of other oceans, can bring refraction and sequence of divergence/convergence of swell rays, to result in larger propagation distances and later arrival (e.g., Rapizo et al., 2018). Waves can be trapped by the currents (e.g., Shrira and Slunyaev, 2014) which fact can cause either acceleration or deceleration. Influence of the vertical gradient of surface currents on kinematics of wave orbital motion is likely, but unknown. Shallows and islands, if encountered by swell on its path, should slow it down. Relative to the deep-water value, group velocity is increased in intermediate depth and reduced in shallow depths, and diffraction of waves into the island shade causes reduction of wave energy and hence the velocity.
Field observations of the swell dynamics, however, are even less frequent than those for waves in the tropical cyclones: the three papers by Snodgrass et al. (1966), Ardhuin et al. (2009), andYoung et al. (2013) are perhaps close to the exhaustive list. Only the first paper is based on in situ measurements, and the modern studies are remote sensing. The satellites do provide global coverage in nearly real time and are an effective way of estimating swell decay, but they cross the great circles rather than follow swells and, as far as swell arrival time is concerned, in their measurements have to rely on assumptions on the swell propagation speeds, which fact is not helpful since these speeds obviously change as the swell propagates.
Thus, field in situ observations are critical for unveiling the very complex nature of swell problem. A majority of world's swells are produced by the Southern Ocean storms with its severe weather around the year which radiates swell waves across the Pacific, Indian and South Atlantic Oceans (Aguirre et al., 2017;Portilla-Yandún, 2018, among others). In this regard, it should also be pointed out that, even in its simplest scenario of swell propagation, the main uncertainties in description of swell propagation are within the proximity of ∼4000 km to its source storm (Ardhuin et al., 2009) which fact makes measurements of the Southern Ocean swells close to their origin critical for understanding their nature. Therefore, a network of wave buoys or systematic deployment of drifting wave buoys in the southern parts of the Pacific, Indian and Atlantic Oceans is proposed to address the problem.

WAVE-DRIVEN DISPERSAL OF FLOATING OBJECTS
The spreading of floating objects on the ocean surface is a fundamental problem of fluid mechanics which has a significant practical value, for marine search and rescue, dispersion of pollution. Note that, while mean drift by surface currents and large-scale eddies are well determined nowadays, particularly with implementation of satellite altimetry, random dispersion of surface drifters remains the open problem (e.g., Soomere et al., 2011). In this regard, the impact of ocean waves with random phases and directional spectrum, remains not accounted for or even well-perceived.
The aviation disaster of Malaysian Flight MH370 drew the public attention to the necessity and complexity of oceanic modeling. In particular, it highlighted the fact that, while modeling of the ocean currents is conducted at the top level, there is no coherent and coupled wave-current modeling. Wave orbital velocities can exceed the geostrophic or wind-driven surface currents, and furthermore wave-induced drift and currents can be comparable to the ocean currents, but are unrelated to them both in speed and direction, and therefore search of debris or missing-at-sea people days and even weeks after the incident are essentially impaired without the coupled wave-wind-current approach. Debris (and hence surface pollution and other floating objects) are carried by geostrophic currents, and by wave-induced currents (Stokes drift and momentum passed by wave breaking). The latter cannot be included on average because it is absent if there is no storm in the area and has to be a subject of new modeling development.
Additionally, random waves with directional spectrum would scatter the floating objects. Formally, turbulent dispersion of a passive tracer caused by a random wavefield is similar to the conventional mechanism of the Taylor dispersion (Batchelor and Townsend, 1956), i.e., particle dispersion by a "conventional" turbulent flow, but with the random velocity field is induced by ensemble of random waves (wave turbulence) and not by conventional turbulent flow. This imposes additional analytical and experimental challenges for investigation of this phenomena (e.g., Herterich and Hasselmann, 1982;Balk, 2001;Falkovich, 2009). In some way, this problem is similar to two-dimensional turbulence applications, and hence can borrow from turbulence research, but will also feed back to the fundamental science because the 2D turbulence has received far less attention than its 3D counterpart. And applications of random 2D vorticity, i.e., when vertical scales are much smaller than horizontal scales, range from boundary layers very near the surface to TCs and upper-ocean circulation.
Southern Ocean, if wave buoys with satellite tracking are released as drifters, is the natural body for introducing, developing, investigating, testing and validating the wavedispersal theories, and implementing them in practice. Innovation due to introduction of the wave scattering can be as big as difference between finding and not finding MH370 and other subjects of search and rescue. Since debris of MH370 have been found along the African coast, location of their origin in the Southern Ocean will require solving the problem of inverse scattering, where new methodology can prove an innovation in its own right. Logistically, this important observational issue can be addresses through deployments of drifting wave buoys in the Southern Ocean as suggested in see Sections "Wave Evolution at Extreme Fetches" and "Severe Extra-Tropical Cyclones at Extreme Fetches."

NON-LINEAR WAVE-CURRENT INTERACTIONS
As far as ocean currents are concerned, these conditions are not common, but are not rare either. Major currents such as Gulf Stream, Kuroshio or Agulhas are well known for harsh seas and high likelihood of abnormal (rogue) waves. Tidal inlets with waves on strong and variable currents are a typical feature of shipping routes in coastal areas. Linear effects of currents on waves, such as refraction, Doppler shift or relative speed with respect to the wind are assumed to be implicitly or explicitly included in wave-forecast models. Our review indicated that in the framework of JCOMM/WMO, since 2001, there is a monthly intercomparison of operational wave model with buoys. Furthermore, operational wave model account for wave/currents interactions by using a surface currents forcing, like the global wave system implemented in Copernicus Marine Service.
Still, absolute majority of the buoys are not located in the Southern Ocean and even the central part of the Atlantic and the Pacific Ocean. So, the validation in the large parts of the ocean, especially at high sea states is mostly missing. Moreover, nonlinear effects are usually left out or even unknown. These include changes to non-linear interactions in presence of currents with horizontal or vertical velocity gradients, wave/current energy and momentum exchanges, non-linear modifications of the wave spectrum .
Thus the wave-current interactions, along with the topics discussed above and the wave-ice interactions in section below in this article, join the list of the least well performing physics in wave-forecast models. This is largely due to lack of understanding based on observations rather than because of the lack of will to improve the wave forecast in Metocean community. In the meantime, bias in model predictions due to such deficiencies is not negligible, and perhaps somewhat surprisingly is not limited to the specific circumstances of major currents or tidal inlets.
Even such a simple linear effect as refraction-induced convergence and divergence of wave energy have been shown to be important factors in modulating the spatial distribution of wave height on the mesoscale (e.g., Ardhuin et al., 2017). One of the most evident examples of wave refraction is wave trains propagating over mesoscale ocean eddies (Figure 3, left panels). Due to the inverted horizontal current shear, one side of the eddy diverges the incoming wave rays, whereas the other side converges the rays (e.g., Mathiesen, 1987). Rapizo et al. (2018) demonstrated that eddy scales as observed from global current reanalysis can potentially create this effect on Southern Ocean swells, but the main impact of these current on the wave-height bias globally is due to another linear effect -change of relative wind speed for the waves on currents (Figure 3, right panel).
Therefore, even linear effects due to currents, which are abundant in the Southern Ocean, have global impact on wave climate and bias of wave modeling. Needless to say that the non-linear wave-current exchanges, which for now not accounted for and not even well understood, can potentially have an enormous influence on the waves due to the very large differences between wave and current kinetic energy. This topic requires a dedicated attention of the community through ongoing satellite observations and through the wave buoy network and drifters proposed in this paper (see sections "Wave Evolution at Extreme Fetches, " "Severe Extra-Tropical Cyclones at Extreme Fetches, " and "Swell Dynamics and Forecasting, With Attention to Arrival Time").

WAVE-ICE INTERACTIONS IN MARGINAL ICE ZONE
Ice edge and the Marginal Ice Zone in the Southern Ocean, unlike in the Arctic, is subject to continuous wave forcing and extreme storms all round the year and hence is an ideal environment for studying wave-ice interactions. Metocean dynamics of the Antarctic Marginal Ice Zone (MIZ) is a topic of great scientific challenge and practical significance. Until recently, the wave forecast models did not predict waves in MIZ (due to lack of knowledge and capability), and in the large-scale models the waves are mostly not taken into account until now.
In terms of knowledge, the wave forecast models have to describe physical mechanisms of wave energy growth, decay, spectrum transformation and wave propagation -in presence of ice. Even if the wind-wave and non-linear energy exchanges are neglected, as the first step, by comparison with dominant energy process of wave decay in ice, such decay by itself accommodates multiple physical mechanisms of wave attenuation, both conservative (wave scattering, reflection and refraction) and dissipative (viscoelastic, turbulent, among others) -see, e.g., Thomson et al. (2018a). Speed of wave-energy propagation also changes in the ice, and provides a family of new dispersion relationships depending on the ice thickness and other mechanical properties (e.g., Collins et al., 2018). The sea ice is a porous material which consists of solid and liquid (brine) phases, with complex elastic, viscous and flexural behaviors as a function of temperature and water salinity -these behaviors define the wave dissipation and propagation and hence need to be known (Wang and Shen, 2010;Mosig et al., 2015, among others). Ultimately, ice is brittle and subject to fatigue under circulating wave forcing, and waves can break the ice (von Bock und Polac, 2016;Williams et al., 2017;Dolatshah et al., 2018). Once this happens, the waves enter a very different dynamic regime: (a) dissipation is driven by floe collisions, rafting, overwash, depends on distribution of floe sizes, and overall appears orders of magnitude weaker (e.g., Bennetts and Williams, 2015;Squire and Montiel, 2016); (b) wave dispersion (shoaling) adjustment can change sign (e.g., Peters, 1950), and (c) wind-forced growth becomes (possibly) not negligible (e.g., Rogers et al., 2018). Furthermore, the broken ice can melt (and the wave fetches will increase, promoting the wave growth) or can re-freeze (and ice cover will increase, arresting the wave growth) see, e.g., Liu et al. (2016). The former depends, among other processes, on waveocean mixing, and the latter on air-sea heat exchange, hence wave forecast in MIZ becomes essentially an air-sea-wave coupled problem (which is less pronounced across the rest of the world's oceans), see Khon et al. (2014).
Analytical theories for some of the processes outlined above, albeit not all, are available, but quantitative (experimental) guidance is fairly limited: typically, these are case studies rather than non-dimensional parameterisations suitable for global wave modeling in a general case, which conclusion highlights the fact that this kind of measurements are extremely difficult and rare (e.g., Meylan et al., 2014). Thus, advance and even progress in wave forecast is restricted and requires urgent attention, observational in the first place. It should be stressed that the practical problem of wave forecast cannot be approached incrementally -the global wave models are run automatically and require quantitative knowledge of all the above processes, not just some of them or a selection of them, as well as accurate determination of the regime change (ice breakup) -in order to predict waves in MIZ or in the solid ice. For example, if the viscoelastic behavior of ice is known (which it is to some extent), but turbulent dissipation in the water boundary layer below ice is not (which it is not), the wave attenuation cannot be estimated with a reasonable degree of confidence. And if the ice breakup is missed or misplaced by the model, the wave decay, both in time and in space, will be completely off the scale.
In terms of the capability of wave forecast in MIZ, this is not just computing facilities, but mostly operational knowledge on the ice fields which poses predictive limitations. The highresolution real-time ice information is as essential for modeling wave-ice decay, as good wind fields are for forecasting the windwave growth. In this regard, sophisticated analytical theories or precise experimental parameterisations are not helpful if the operational ocean-ice models or satellite observations are not able to provide the relevant properties of ice. Thus, there will always be a gap between research and operational wave modeling, and the practical applications need to balance between exact knowledge and its realistic implementations.
The coupled nature of wave forecast in Antarctica, furthermore, highlights the fact of reciprocal importance of waves for the oceanographic forecast (and, more generally, for air-sea interaction modeling). If waves break the ice and, as a result, it melts faster in spring/summer, this can have significant impact on air-sea fluxes (and not only heat fluxes), even if the ice refreezes in autumn/winter. Note that the first-year ice will be easier to break next summer, and thus the positive feedback loop may accelerate.
Because the presence of waves is more the rule than the exception at the margins of Antarctic ice, this changes the type of ice that forms during the colder months, which can have a profound impact on the heat fluxes and thus the rate of ice growth. Specifically, new ice in the presence of waves will tend to be frazil and pancake ice, and will tend to be sheet ice (starting as nilas) in the case without waves. With pancake and frazil ice, liquid water is directly exposed to the cold air, allowing faster freezing (Doble, 2009). With sheet ice, heat must pass through the insulating ice (thus, slower freezing). The ice type also affects the albedo (so heat flux, again) and the surface roughness (and thus drag on the atmosphere).
Thus, wave-ice interactions, along with the Metocean topics above in this article, is a poorly understood type of ocean-wave dynamics, which, correspondingly, leads to poor performance of wave-forecast models in respective conditions. Like the other topics, the main problem in advancing the fundamental understanding and practical modeling of such conditions is lack of observations, and the most suitable environment for such observations is the Southern Ocean. Necessary observations, in addition to wave and flux observations proposed in see Sections "Wave Evolution at Extreme Fetches, " "Severe Extra-Tropical Cyclones at Extreme Fetches, " and "Swell Dynamics and Forecasting, With Attention to Arrival Time, " will require measurements of wave and ice properties within the Marginal Ice Zone. It is suggested to use Antarctic-going ships of opportunity for this purpose.

METOCEAN CLIMATOLOGY IN THE SOUTHERN OCEAN
The Southern Ocean is the least studied ocean area in terms of in situ oceanographic and Metocean observations. In the meantime, it demonstrates the fastest growth of winds and waves, both in the mean and in extreme percentiles, by comparison with the other Oceans, at least over the era of satellite remote sensing observations, i.e., since mid 1980s . Seasonally, this most dynamic Metocean region is the only one which demonstrates positive trends for the ocean winds well above mean global values over 2/3rds of the year [except Southern spring, see Table 1 and Zieger et al. (2014)].
Metocean climate, apart from the winds and currents also includes ice-covered area and ice thickness whose trends in the Southern Ocean are different to the Arctic and in need of dedicated investigations. Overall, Metocean characteristics, particularly their consistent trends at large scales in space in time, indicate regional climate changes, which may be also connected to the global climate behaviors. If subject to vigorous measurements, such characteristics and their trends can serve as climate proxies, potentially more robust by comparison with point characteristics (such as temperature) because the nature of Metocean properties is necessarily an integral over large areas.
Thus, major gaps in the global Metocean climatology come from the Southern Ocean where in situ observations, particularly at long-term and systematic basis, are virtually absent. Synergy of the proposed wave and air-sea buoy networks, ice and other Metocean measurements of opportunity (see sections "Wave Evolution at Extreme Fetches, " "Severe Extra-Tropical Cyclones at Extreme Fetches, " "Swell Dynamics and Forecasting, With Attention to Arrival Time, " "Wave-Driven Dispersal of Floating Objects, " "Non-linear Wave-Current Interactions, " and "Wave-Ice Interactions in Marginal Ice Zone") will help to address this topic of practical and research significance.

MEASUREMENTS IN THE SOUTHERN OCEAN
The paper is proposing a dedicated measurement program for the Southern Ocean, and therefore in this Section we will briefly review available in situ observations. Over the last several years, there have been a number of efforts to start Metocean measurements in the Southern Ocean. Most encouraging are attempts of permanent buoy deployments by the Australian Integrated Marine Observing System (Schulz et al., 2012) and by the Ocean Observatories Initiative of the United States National Science Foundation: https://oceanobservatories.org/ array/global-southern-ocean/ (see locations in Figure 4, top panels). A number of moored and drifting buoys were deployed south of New Zealand by the Metocean Solutions and the Royal New Zealand Navy: http://www.metocean.co.nz/southernocean/ (Figure 4, bottom left); Metocean observations are conducted by the University of Cape Town in oceanographic voyages of SA Agulhas II in the Southern African sector of the South Ocean all the way to Marginal Ice Zone (Figure 4, bottom right).
Recently, some investigators have begun using autonomous platforms for Metocean measurements in the Southern Ocean. Thomson and Girton (2017) used a wave glider in the Drake Passage during the austral summer of 2017, with a particular focus on measuring directional wave spectra and wind stress (Thomson et al., 2018b). Schmidt et al. (2017) also use a wave glider to evaluate model winds in the Southern Ocean. These and other mobile platforms continues to operate in 2018-2019. Thus, an expansion of in situ wave observations from autonomous surface platforms in addition to traditional moorings, is also likely.
These few deployments, however, while very promising, is literally a drop in the ocean of the most powerful winds and waves. Hundreds of the wind-wave buoys in the Northern Hemisphere, and very few South of equator such as those off the coast of Brazil in the path of Southern swells (Pereira et al., 2017).   of the Metocean observations, and the emptiness of the Southern Ocean where most of the actual problems of the modern Metocean science and applications can and need to be solved (Young et al., 2017). While permanent buoy deployments can be a substantial challenge, drifting buoys, wave gliders and other moving platforms can prove feasible and valuable solution of this challenge for a dedicated international effort.
For an area as geographically remote as the Southern Ocean, remote sensing offers obvious benefits in providing wind and wave data. In terms of wind measurements, there are three potential platforms (radiometers, scatterometers and altimeters). For wave data there are also three options (altimeters, synthetic aperture radar and CFOSat). A number of studies have already looked at global climatology of wind speed and wave height, including the Southern Ocean (Zieger et al., 2009;Vinoth and Young, 2011;Young et al., 2011Young et al., , 2017Takbash et al., 2018;Young and Donelan, 2018;Ribal and Young, 2019;Young and Ribal, 2019). These studies, however, are limited to wind speed and significant wave height and suffer from very limited possibilities for Southern Ocean Calibrations.
Our capability to measure directional waves, up to the early 1990s with the launch of ERS-1, was restricted to few areas in the world where buoy data, mainly, was available. Synthetic Aperture Radar (SAR) is the only satellite instrument so far capable to measure the directional spectrum, despite some limitations in its high frequency part. SAR data have been available since then, with a myriad of satellites yielding over 25 years of directional spectra with global coverage. Sentinel-1A and its twin 1B are currently operational, sharing the same orbit plane and therefore with a greater revisit rate. Sentinel-1C is scheduled to be launched in the next 3-4 years, which will increase the temporal sample of the constellation. The recent launch of CFOSat, which carries a unique scanning wave scatterometer (SWIM) provides the potential to measure the full directional spectrum for components longer than 80 m. This instrument has great potential to open up a new era of wave measurements in environments such as the Southern Ocean. In the context of the proposed network, directional buoy measurements in the scarcely sampled Southern Hemisphere will contribute to the effort to validate such satellite wave observations.

CONCLUSION
Metocean measurements in the Southern Ocean -marine winds and currents, surface waves and swells, ice cover and thickness, among others -are either critically important or, at the very least, can contribute to solving and addressing problems of major significance. Without the Southern Ocean in situ observations, it is not possible to resolve the issues of limiting fetches, extreme Extra-Tropical cyclones, swell propagation and attenuation, wavecurrent interactions. The topics of wave-induced dispersal of floating objects, wave-ice interactions in the Marginal Ice Zone, Metocean climatology and its connection with the global climate cannot be complete in general case without benchmarking the behaviors of these phenomena against observations in this most dynamic area of the global Ocean.
The paper outlines and highlights important gaps in understanding the Metocean processes and suggests a major observational program for this large, but poorly investigated part of the World Ocean. This would include distributed buoy network (drifting and stationary) and autonomous surface vehicles (ASV), intended for measurements of waves and air-sea fluxes in the Southern Ocean.