Wind Speed Measured from Underwater Gliders Using Passive Acoustics

Wind speed measurements are needed to understand ocean–atmosphere coupling processes and their ef- fects on climate. Satellite observations provide sufﬁcient spatial and temporal coverage but are lacking ad-equate calibration, while ship- and mooring-based observations are spatially limited and have technical shortcomings. However, wind-generated underwater noise can be used to measure wind speed, a method known as Weather Observations Through Ambient Noise (WOTAN). Here, we adapt the WOTAN tech- nique for application to ocean gliders, enabling calibrated wind speed measurements to be combined with contemporaneous oceanographic proﬁles over extended spatial and temporal scales. We demonstrate the methodology in three glider surveys in the Mediterranean Sea during winter 2012/13. Wind speeds ranged from 2 to 21.5ms 2 1 , and the relationship to underwater ambient noise measured from the glider was quantiﬁed. A two-regime linear model is proposed, which validates a previous linear model for light winds (below 12ms 2 1 ) and identiﬁes a regime change in the noise generation mechanism at higher wind speeds. This proposed model improves on previous work by extending the validated model range to strong winds of up to 21.5ms 2 1 . The acquisition, data processing, and calibration steps are described. Future applications for glider-based wind speed observations and the development of a global wind speed estimation model are discussed.


Introduction
Quantifying air-sea fluxes is critical in understanding the weather-ocean-climate system. Numerical models need forcing by in situ measurements at an increasingly higher spatial and temporal resolution (Zhang et al. 2006). Sea Denotes content that is immediately available upon publication as open access.
Supplemental information related to this paper is available at the Journals Online website: https://doi.org/10.1175/JTECH-D-17-0209.s1. surface wind speed is a key parameter in forcing numerical models, as well as in quantifying turbulent air-sea fluxes and gas exchanges (Wanninkhof 2014).
Observations made from satellites can provide wind speed data with near-global coverage over the ice-free oceans with a spatial and temporal resolution of about 0.258 and 24 h (e.g., QuikScat, ASCAT), respectively. However, there is a lack of calibration of in situ observations outside the tropics and away from coasts, or in high wind speed (.18 m s 21 ) conditions (Bourassa et al. 2010). In situ monitoring of weather conditions over the ocean is difficult to achieve. Ship-based observations are affected by airflow distortion (Moat et al. 2005) and are sparse (Kent 1998). Making observations from moored meteorological buoys in the long term (years to decades) presents difficulties, such as damage or loss caused by the roughness of the sea surface environment (storms, ice, fishing activity, vandalism, etc.). The spatial coverage offered by moored meteorological buoys and shipbased observations remains limited, as is their ability to provide observations of extreme wind events, because of sensor and platform limitations (Weller et al. 2008).
The Weather Observations Through Ambient Noise (WOTAN) technique (Vagle et al. 1990) enables monitoring of the sea surface weather conditions from underwater, away from the rough sea-air interface, with no difficulties induced by extreme weather events. This approach relies on the analysis of underwater ambient noise, generated by the excitation of the sea surface by the weather conditions. Surface-generated noise can be recorded up to at least 6 km deep (Barclay and Buckingham 2013). The unique characteristics of the main underwater sound sources (e.g., spectrum shape, time variability) allow wind-generated noise to be isolated and quantified (Fig. 1).
Underwater noise generated by surface weather conditions was first studied in the mid-twentieth century, because of its effect on the performance of submarine detection systems (Urick and Kuperman 1989). In the open ocean, and in the absence of sound from marine life or nearby ships, the main source contributing to underwater sound in the frequency band from 500 Hz to 50 kHz is the sound produced by surface weather conditions (Black et al. 1997;Vagle et al. 1990;Wenz 1962). The action of the wind on the surface of the ocean induces air bubbles, spray, splash, and turbulence noise, which contribute to underwater ambient noise (Carey et al. 1993). These complex surface processes are influenced by multiple parameters, such as wind speed, wind duration, fetch, and marine atmospheric boundary layer (MABL) stability. In this study, we will focus on the instantaneous wind speed only, as a first-order approximation. A linear relationship between the logarithm of the surface wind speed and the sound pressure level (decibels) was first reported in shallow water (;40 m) (Piggott 1964), and then extended to deep water (5000 m) in the 1-10-kHz frequency range (Crouch and Burt 1972;Shaw et al. 1978). An empirical linear relationship between the surface wind speed and the sound pressure (mPa) was then proposed in the 4-15 m s 21 wind speed range (Vagle et al. 1990), which is now widely used (Nystuen and Ma 2002;Riser et al. 2008;Vakkayil et al. 1996).
The wind speed measured using this technique is relative to the sea surface, so it is applicable to estimation of heat and moisture fluxes and wind stress (Bourassa et al. 2010). The WOTAN technique spatially averages over an area dependent on the frequency used, the depth of the measurement, and the sound speed profile at the measurement site. At 3 kHz, it varies from 0.1 km 2 for a measurement depth of 100 m to 10 km 2 for a measurement depth of 1000 m (Vagle et al. 1990). For increasing frequencies, the sound absorption coefficient increases, thus the listening area decreases. This spatial scale corresponds to the scales considered in air-sea interaction studies (Bourassa et al. 2010).
The linear relationship between the sound pressure p (mPa) and the surface wind speed U (m s 21 ), where b (mPa) and s (mPa m 21 s) are respectively the offset and slope of the linear regression, is widely used in most recent studies (Nystuen and Ma 2002;Riser et al. 2008;Vakkayil et al. 1996). These studies  Wenz (1962)]. In the absence of heavy rain events or nearby biotic activity, wind-generated noise is predominant in the 500 Hz-20 kHz frequency range. agree on a low wind speed limit of around 2 m s 21 , below which the wind-generated sound is below the background noise level. They also agree on a high wind speed limit around 15 m s 21 , above which measurements are scattered and the correlations are poor, and propose an explanation based on the hypothesis of bubble-layer attenuation of surf noise (Black et al. 1997;. It is worth noting that these previous studies contained few observations of high wind speeds. Ocean gliders are autonomous underwater vehicles, carrying sensors to monitor the ocean. They perform long autonomous missions (several months to a year, and several thousand kilometers) and provide highresolution (;2 h, 2 km) hydrographic profiles (Testor et al. 2010;Rudnick 2016). Glider measurements are not affected by extreme weather events. Their unique way of moving through the water column (buoyancy driven with no propellers) makes them extremely quiet and therefore very suitable for passive acoustic monitoring. Gliders also measure temperature and salinity profiles, from which sound velocity profiles can be derived. This collocated information on the acoustical properties of the water column is of considerable value for soundscape studies. Passive acoustic monitoring (PAM) sensors have been successfully deployed on ocean gliders for cetacean monitoring purposes (Baumgartner and Fratantoni 2008;Baumgartner et al. 2013;Klinck et al. 2012;Moore et al. 2007). The combination of hydrographic profiles with surface weather measurements on submesoscales (i.e., kilometer scale), tracking a weather event or monitoring a selected area for many months, would be highly valuable for air-sea interaction studies.
This paper presents a novel method for measuring surface wind speed using glider-borne underwater ambient noise measurements, from the subsurface to 1 km deep. We deployed gliders equipped with PAM sensors in the Mediterranean Sea during winter 2012/13 in the framework of the Deep Water formation Experiment (DEWEX) experiment (Testor et al. 2017). The gliders recorded 4 months of acoustic data, with recurring opportunities to compare our wind speed estimates with Météo-France meteorological buoys in the area. Focusing on the sound pressure level in the 3-kHz third octave band, which shows the most dynamic response to wind speed, we estimate the surface wind speed around the glider's position (0.1-10 km 2 ) throughout the 2-21.5 m s 21 wind speed range. Section 2 describes the experiment and the associated datasets. Section 3 presents the data acquisition and processing methods, and the results. Section 4 presents the wind speed derivation model and its performance. In section 5, we consider the broader application of the model to different experiments or regions, and we discuss future improvements and the contribution to wider monitoring activities.

Instrumentation and field measurements
Passive acoustic measurements were made using an Acousonde B003A-HF datalogger (Fig. 2), developed by Greeneridge Sciences, Inc. The Acousonde is a self-contained underwater acoustic recorder comprising two hydrophones, sensors for attitude, orientation, depth, and temperature, a digital recorder, and a field-replaceable battery (Burgess 2010). The core of the sensor consists of a high-frequency hydrophone (capable of sampling up to 232 kHz), with a sensitivity of 2204 dB re 1 V mPa 21 . A six-pole linear phase antialiasing filter is used with cutoff at 10 kHz (23 dB) and 42 kHz (222 dB), and a linear high-pass filter with cutoff at 10 kHz (23 dB) (see Fig. S1 in the online supplemental material). Data are stored on a 128-GB flash memory, with a 16-bit sampling resolution. An external three-D-cell tethered battery pack allows up to 200 h of recording. The Acousonde operates autonomously, and has its own battery, memory, and programmed mission. Data processing is undertaken after the sensor is recovered. Initially developed to be attached to marine mammals (Cazau et al. 2017), it has also been used on ocean gliders (Nott 2015). Details of the sampling method used in this study are given in section 3.
The platform we used in this study is the Slocum glider, developed by Teledyne Webb Research. It is driven by buoyancy changes, controlled by 500 cm 3 of oil pumped into or out of a swim bladder, inducing a vertical motion in the water column, from the surface down to 1000-m depth. Fixed wings convert the vertical velocity into forward velocity, internal battery displacements enable pitch and roll management, and a moving rudder enables direction changes. This novel way of propulsion makes it a very quiet platform between the oil pumping phases that occur at the apogee and perigee of each dive (;2 h for 1000-m dives). Water turbulence around the sensor induces flow noise, proportional to the glider's speed in the 5-50-Hz frequency band (Dos Santos et al. 2016;Erbe et al. 2015), with no effects regarding the windgenerated sound levels over 1 kHz. The version of the Slocum gliders we used can complete autonomous missions up to 3 months and 1000 km long. Along with the external PAM sensor, the gliders were equipped with integrated temperature, salinity, and pressure sensors.
The experiment took place within the Mediterranean Ocean Observing System for the Environment (MOOSE) framework, which aims to monitor the oceanographic variability of the northwestern Mediterranean Basin over a continuum of spatial and temporal scales (http:// www.moose-network.fr), and the DEWEX experiment (Testor et al. 2017 (Fig. 3). These lines closely pass the two Météo-France meteorological buoys Lion and Azur, defining two distinct experiment sites in which continuous in situ surface wind measurements are available, thus allowing recurring calibration and validation of the glider measurements.
These two experiment sites present different geographical characteristics, hence different spatial variability scales in their wind fields. The Lion buoy (42.068N, 4.648E) is in the middle of the Gulf of Lion, 100 km from the nearest shore, in waters 2300 m deep, in an open sea area affected by the mistral and tramontane winds, and away from the main shipping lanes. The Azur buoy (43.388N, 7.838E) is 50 km from the shore, in waters 2300 m deep, in an area affected by cyclogenesis in the Gulf of Genoa (Rainaud et al. 2016), and close to the alongshore shipping lanes. Each buoy measurement is considered to represent the conditions for an area around the buoy's position, whose size depends on the spatial variability of the measured parameters and the geographical position of the buoy. To better estimate the area represented by each buoy wind speed measurement, we used output from the atmospheric model Application of Research to Operations at Mesoscale-Western Mediterranean (AROME-WMED; Nuret and Fourrié 2011) to compare the wind speed at the buoy position with the wind speed for the rest of the Gulf of Lion, for the two sites, Lion and Azur (Fig. 4). Around the Lion buoy, the correlation (0.93), root-mean-square (RMS) difference (2.1 m s 21 ), and mean difference (0.2 m s 21 ) with the wind speed at the buoy position are good within 80 km of the buoy, with a rapid deterioration outside. We define this as the Lion confidence area, in which we use the buoy 1-min average measurements of wind speed at 10 m as ground truth. Around the Azur buoy, this radius decreases to 40 km. We defined the Azur confidence area as a 40-km radius around the Azur buoy that provides hourly average measurements of wind speed at 10 m, with mean correlation, RMS difference, and difference of 0.86, 2.5 m s 21 , and 0.2 m s 21 , respectively ( Fig. 4). We undertook four PAM glider deployments around the two sites, Lion (missions Asicsmed and MistralsT02_01) and Azur (missions MooseT00_23 and MooseT00_25). These cover 138 days of data, 37 of which are within the confidence area of one of the buoys and can be used to derive a relationship between sound pressure and surface wind speed (Table 1). Multiple factors can induce differences between our measurements in the different experiments: the two buoys have different measurement sampling frequencies; shipping activity, and therefore its generated noise, is different at the two sites; and our PAM sensors have different sensitivities per experiment, because of the varying mounting position on the glider, and the difficulty of using Acousonde to produce absolute sound level measurements (Wiggins 2013). For these reasons, we decided to process each experiment separately. We merged the two Lion datasets, Asicsmed and MistralsT02_01 (referred to as Lion dataset), as both deployments used the same setup in the same area without any maintenance operations in between.

Acoustic data sampling and preprocessing procedures
We designed the data acquisition protocol to evenly distribute the 200 h of recording time over the 3-month glider mission. The shortest sample allowed by the PAM sensor, 1 min, is sufficient for analysis of wind-generated noise (Nystuen and Ma 2002). The PAM sensor recorded 1 min every 10 min, to allow a mission duration of 80 days. The sampling rate was set at 50 kHz, constrained by the PAM sensor's onboard memory capacity. This sampling method produces 27 GB of data every month.
We adapted the WOTAN technique (Vagle et al. 1990). We processed each 1-min sample individually, extracting the sound level on third octave frequency bands [third octave level (TOL)], according to the standards [American National Standards Institute (ANSI) S1.11-2004] and using the software PAMGuide (Merchant et al. 2015) on 1-s-long non overlapping Hanning windows.
To focus on the wind contribution to the underwater ambient noise, any other sound sources have to be detected and removed from the dataset. The noise generated by the glider's oil pump is a loud broadband signal, lasting for several minutes. Slocum gliders generally activate their pump only during the apogee (surfacing) and perigee (bottom inflexion) phases of a dive. As we programmed only 1000-m-deep dives, we removed data acquired during apogee (depth , 20 m) and perigee (depth . 950 m) to focus on the quiet gliding phase of each dive (Fig. 5). The Lion dataset contains 4021 remaining 1-min samples, the equivalent of 40 days (Table 1).
To eliminate transient sounds (e.g., clicks, whistles, and glider's motors), we assumed that wind speed is constant over a 1-min sample and kept only the minimum TOL values. We therefore condense a 1-min sample spectrogram to a single spectrum in which remains no signal of the transient sounds (Fig. 6).
The pitch of the glider varies significantly between ascending (approximately 268) and descending (approximately 2268) phases. To assess the effects of this variation  on the measurements, we compared the estimated wind speed during climbing and diving profiles. There is an average overestimation from the ascending profiles of 0.3 m s 21 . An ANOVA test on these time series gives a p value of 0.3. Therefore, we neglected the effect of the variation of the pitch in this study. The glider profiling behavior implies significant variations in the depth of the measurements, which affect the received sound level. The effects of refraction, attenuation, spreading loss, and directionality of the sound sources are accounted for, for each depth and frequency, depending on the sound attenuation and velocity profile (Vagle et al. 1990). Therefore, TOL(h, f), measured at depth h and frequency f, can be corrected to an equivalent TOL 0 ( f ) at the surface: where where r is the horizontal distance from the source to the point directly above a hydrophone at h; l is the pathlength from the source to the receiver, including refraction effects; u is the angle between the surface and the path to the receiver; a is the frequency-dependent attenuation coefficient for bubble-free water. For each deployment, we calculated the b profile using the average sound velocity profile, and attenuation, obtained from the glider's temperature and salinity measurements. The average sound velocity and b profiles are shown in Fig. S2 in the online supplement. The applied correction for a depth of 1000 m is 1.3 dB, which means that the received sound pressure (mPa) is 14% lower that the sound pressure at the surface. The associated wind speed estimation correction varies from 0 to 2.5 m s 21 with the sound speed in the 0-25 m s 21 range. The nonlinearity observed at 9-10 m s 21 is due to the regime change. Because of the change of the derivation slope, a similar sound pressure correction will lead to a wind speed correction divided by 4 in the high wind speed regime than in the low wind speed regime.

Wind speed derivation a. Wind speed derivation model
Previous studies Nystuen and Ma 2002;Vagle et al. 1990;Vakkayil et al. 1996) were conducted on wind speeds limited to the 3-15 m s 21 range (Table 2). These studies found a linear relationship [(1)] between sound pressure at 8 kHz (p 8kHz , mPa) and 10-m wind speed (U 10 , m s 21 ). They discussed limitations at high wind speeds (Vagle et al. 1990).
In this study, we use where p fc is the sound pressure for the third octave band of central frequency f c , relative to the minimal pressure p min (f c ) observed on the deployment. The correlation between p fc measured by the glider and U 10 measured by the buoy depends on f c . The best correlations are in the 2-10-kHz frequency band, where the wind-generated sound is predominant (Fig. 7). The poor correlations in frequency bands below 1 kHz is due to the predominance of non-wind-dependent sound sources (e.g., distant shipping). The poor correlations in frequency bands above 10 kHz can be explained by the attenuation of sound in the high-frequency bands observed during high wind speed events (.15 m s 21 ) (Fig. 7). We chose to use the relative pressure in the 3-kHz third octave band, which shows the most dynamic response  to wind speed and shows no attenuation effect during high wind speed events (Fig. 7). Similar results are obtained if other frequency bands (within 2-10-kHz range) are used.
To train our wind speed derivation model, we consider the median values of the acoustic parameter and the measured wind speed over each glider profile. As a result, contaminating sounds of duration of several 1-min samples (e.g., nearby ship) will be filtered out, and the time and length scales are the usual glider profiling scales (;2 h, 2 km for a 1000-m vertical profile). The number of observations is unevenly distributed over the wind speed range (Table 3). Therefore, to give an equal weight to each wind speed during the regressions between p fc and U 10 , we calculated the median sound pressure in each 2 m s 21 wind speed bin.
During our experiment, we observed several high wind speed events of 6-48 h of wind speed above 15 m s 21 , with the glider in the confidence area around the buoy. Our dataset covers a broader wind range (2-21.5 m s 21 ) than previous studies, allowing us to investigate the relationship between p fc and U 10 for wind speed regimes above 15 m s 21 . To assess the validity of the historical single linear model [(1)] for our observations, we fitted a linear model, of parameters b single and s single (Table 4), to the wind speed range below 15 m s 21 , simulating the wind speed range of the previous studies. Its coefficient of determination R 2 5 0:86 indicates a good fit for the 2-15 m s 21 range. We extrapolated this to the high wind speed range and tested it against the data available in the Lion dataset, revealing a tendency of the single linear model to overestimate wind speed in the higher wind speed range (15-21.5 m s 21 ) (Fig. 8).
To improve the agreement at high wind speed we investigated the hypothesis of two distinct sound production regimes, related to the physical processes of wind-generated underwater noise (bubbles, spray, and splash). Whitecaps can be observed on sea surface images at wind speeds as low as 3.7 m s 21 (Callaghan et al. 2008), with this threshold varying with air and sea temperature (Monahan and O'Muircheartaigh 1986). Previous studies using ambient noise show a low wind speed detection threshold at 1.5-4 m s 21 (Table 2). A wind speed threshold in the range of 9-11 m s 21 is commonly used in the description of physical sea surface processes, such as the generation of spume drops, torn off wave crests by high winds (Monahan et al. 1983;Monahan et al. 2017), spray generation function (Large and Pond 1981), and drag coefficient (Foreman and Emeis 2010). As these processes all impact wind noise production, it is likely we will observe regime change in the underwater noise production at wind speeds above this threshold. We therefore defined a low wind speed range (U 10 , 11 m s 21 ), where whitecaps are responsible only for the sound production, and a high wind speed range (U 10 . 9 m s 21 ), where additional physical processes contribute to the sound production. We fitted two linear regressions for low wind speed parameters b low and s low , and for high wind speed parameters b high and s high ( Table 4). The low wind speed model is very close to the single linear model, despite being fitted on a dataset of narrower wind speed range. The two-regime linear model allows a good description of the observations for the complete measured wind speed range (Fig. 8). As an alternative to the two-regime linear model, we investigated the possibility of fitting a power-law (p 5 a 1 U b ) relation. We obtained good results with the quadratic modelp 5 459:6 U 2 -but with a tendency to overestimate low wind speeds. We elected to use a two-regime model instead, as the fit was better and can be explained by a change in physical processes at elevated wind speeds.
The physical processes explaining the regime change appear in a continuous and smooth way, around a wind speed limit that is believed to depend on multiple  parameters (sea state, current, fetch, temperature,. . .).
To better describe this smooth transition between two different linear relationships on both sides of an undetermined joint point, we can use a hyperbola as a transition model (Watts and Bacon 1974). Therefore, we estimate U 3kHz , the wind speed derived from the sound pressure p 3kHz , as U 3kHz 5 U lim 1 s low 3 s high 2(s low 1 s high ) (p 3kHz 2 p lim ) 1 s high 3 s low 2(s high 2 s low ) where › 5 10 000 is the radius of curvature at the joint point (p lim 5 4:5 3 10 4 mPa; U lim 5 10:4 m s 21 ) between the low (s low , b low ) and high (s high , b high ) wind speed linear models.

b. Performance of the wind speed derivation model
Prior to the estimation of the performance of the wind speed derivation model, we detected and removed some data points from the analysis that we suspected to be corrupted by an anemometer failure. They all occur during the 23-28 January 2013 period, during which the buoy measurements are sparse, and disagree with both the glider estimation and the model output (Figs. 9, 11). A comparison between the wind speed derived from the sound pressure using (5), U 3kHz , and surface wind speed measured by the buoy, U 10 , shows a good agreement for the 2-20 m s 21 range. The root-mean-square of the error (U 10 2 U 3kHz ) between the estimation and the FIG. 9. Wind speed U 10 measured at the buoy vs U 3kHz derived from acoustic measurements using the two-regime linear model with hyperbolic transition [(3)]. The RMSE is calculated on the low and high wind speed regimes. Corrupted data points, excluded from the analysis, are shown in red. observation is 2 and 1.4 m s 21 for the low and high wind speed regimes, respectively (Fig. 9). Because the glider traveled repeatedly toward and away from the buoy, the wind speed derivation model, calibrated against U 10 using the data collected within the buoy's confidence area, can be used when the glider is away from the buoy, extending the spatial coverage of the wind speed measurements. The time series of U 10 measured at the buoy position, U 3kHz derived from the glider measurements, and U AROME from the AROME_WMED model output collocated with the glider position show a good fit between U 3kHz and U AROME when the glider is away from the buoy confidence area (light blue color coding) and when the buoy's anemometer failed to provide data (Fig. 10).

Application to other glider campaigns and implications
The same technique was applied to the MooseT00_23 and Moose T00_25 campaigns, around the Azur meteorological buoy. For the MooseT00_23 experiment, the wind speed U 3kHz derived using the associated hyperbolic transition model [(3)] fits U 10 measured by the buoy throughout the 2-20 m s 21 wind speed range, with a root-mean-square error of about 2.2 m s 21 in the low wind speed regime and 1.6 m s 21 in the high wind speed regime (Fig. 11). This long experiment (7 weeks), with eight occurrences of the glider surveying the Azur buoy area, illustrates the opportunity of repeated calibrations using the buoy measurements to estimate the surface wind speed along the glider path (Fig. 11).
For the MooseT00_25 experiment, the wind speed U 3kHz derived using the associated hyperbolic transition model [(3)] fits U 10 measured by the buoy throughout the 2-17 m s 21 wind speed range, with a root-mean-square error of about 1.6 m s 21 (Fig. 12). The poor results in the 16-18 m s 21 wind speed bin can be explained by the low amount of data (six dives), all of which were recorded during one single event on 1 April. The buoy measurements during this event show a 6-h-long plateau at 17.0 m s 21 exactly that could suggest an anemometer failure (Fig. 12).
The noticeable differences in the absolute values of the parameters (Table 5)  between the two Acousonde sensors used, and by the unavoidable variation in sound level with the orientation of the recorder in each deployment. The internal layout of the Acousonde makes it asymmetrical, with the electronic board likely to affect the measurements of the sounds coming from one side. Its cylindrical shape makes it difficult to accurately position the electronic board between the transducer and the glider's hull. This interdeployment sensitivity variability is believed to be responsible for the quantitative differences between our three experiments. Despite this variability in absolute sound level, these results show that once calibrated with buoy data, a relative sound level can be used and produces consistent results. The recent availability of PAM sensors integrated into ocean gliders should reduce this interdeployment sensitivity variability. Being able to calibrate the whole PAM/ glider system in situ once, knowing that the sensitivity will be the same for each following deployment, will allow using the PAM glider for wind speed estimation in areas where in situ calibration data are not available.
The method we propose here is based on relative sound pressure levels. Building a similar model using absolute sound pressure levels would improve the method, since it would negate the need for a specific in situ calibration, although the consistency of these absolute measurements (e.g., for different sensor orientations) would need to be tested.

Discussion
The adaptation of the WOTAN technique (Vagle et al. 1990) to recordings from ocean gliders provides wind speed measurements with long endurance, large spatial coverage, and resilience to severe weather conditions. During this study, we were able to collect data without maintaining ship presence during several high wind speed events. Previous studies were mostly limited to wind speeds below 15 m s 21 (Table 2), so they were unable to successfully characterize higher wind speeds. The significant amount of high wind speed conditions encountered during this experiment allows us to propose a model extending the WOTAN technique measurement range to strong wind conditions, up to 21.5 m s 21 . No evidence of an upper limit of our model, where the sound at 3 kHz would be attenuated, has been found here. Further experiments are necessary to determine whether this critical wind speed exists. The positions of the meteorological buoys relative to the coast and the direction of the main winds limits the fetch to 100 km at the Lion buoy, and 50 km at the Azur buoy. This does not allow a fully developed sea at wind speed higher than 10 and 7 m s 21 at the sites Lion and Azur, respectively. Future experiments in the open ocean (e.g., Southern Ocean) will provide more data at high wind speed and with unlimited fetch, which will allow studying the impact of sea state on the wind noise and assessing the validity of our model at high wind speed. The use of sound pressure level at 3 kHz in this study was driven by the instrument self noise at ;2 kHz and the possible contamination by traffic noise at lower frequencies on one side, and the limitations observed for the high frequencies at high wind speeds on the other side. It is however not necessary to focus on this specific third octave band. Moreover, monitoring higher wind speed conditions than the ones presented here may necessitate the use of multiple frequencies. At high wind speeds, lower frequencies maintain a linear response to wind speed, while higher frequencies show decreasing amplitude/intensity/power with increasing wind speed. Further work is required to constrain the relation between higher frequencies and high wind speeds in relation to sea surface physical parameters (sea state, wave age, energy flux. . .). In the particular case of surface wind speed estimation, the use of p 3kHz for wind speeds up to 25 m s 21 , and p 500Hz for higher wind speeds, when wind-generated noise is loud enough to prevent masking by shipping noise, could be considered, allowing a wide wind speed range to be monitored.
This paper focuses on providing submesoscale measurements of surface wind speed combined with oceanographic profiles. We recorded 1 min every 10 min, and considered the median of the acoustic measurements to provide a robust wind speed estimation for the duration of each glider profile-in this case measured every 2 h and 2 km. This method allows the acquisition of high-resolution in situ forcing data, in weather conditions or locations of interest (e.g., storms, remote regions) for numerical models and satellite calibration (Schmidt et al. 2017). Recent progresses in the integration of PAM sensors into ocean gliders now enables continuous recording during a monthlong glider campaign (limited by the glider's battery). The variation of the listening area with the depth of measurement (Vagle et al. 1990) can be used to choose to measure wind speed over different spatial scales. For example, limiting the glider measurements to the upper 100-200 m would enable monitoring of a sea surface area of 0.1-0.5 km 2 , hence studying of smaller-scale processes than with typical 1000-m glider dives (e.g., land sea breezes, transient events). Continuous sampling throughout the glider dive, combined with sound speed and attenuation profiles measured by the glider, will allow better evaluation of the vertical distribution of ambient noise.
The combination of surface wind speed measurements with glider oceanographic profiles allows better observation and quantification of air-sea interaction processes, and monitoring of the associated processes in the ocean. Surface wind speed measurements allow correction of the bulk SST to obtain skin SST values, which are fundamental for the quantification of air-sea interaction processes (Alappattu et al. 2017). Heat and freshwater transport can be monitored from the air-sea interface to the water column, linking heat content changes in the ocean to meteorological events (Grist et al. 2016).
Wind noise is caused by complex physical processes, such as wave breaking and bubble inclusion (Wenz 1962), driven by weather, sea surface conditions, and history. As a first approximation, we estimate the instantaneous surface wind speed from measurement of wind noise. We believe that further experiments (e.g., glider campaigns together with surface processes monitoring) could improve the wind speed estimation, by
For future deployments, we suggest that in situ calibration of the PAM glider should be performed at the beginning of each deployment. Deploying the PAM glider at the surface from a small boat, next to a calibrated hydrophone, in various noise conditions (e.g., engine on/ off) should be sufficient to allow the acquisition of absolute sound levels. Also, the use of an integrated PAM sensor will reduce the interdeployment variability in the positioning of the sensor, and therefore allow use of the same in situ calibration for multiple deployments.

Conclusions
Surface wind speed can be measured remotely, derived from underwater acoustic measurements. The proposed two-regime linear model yields improved results in high wind speed conditions, extending the wind speed range to 2-21.5 m s 21 .
The PAM sensor can be mounted on an ocean glider, diving from the surface to 1000 m deep, thus allowing the surface wind speed measurements to be combined with the oceanographic profiles (e.g., temperature, salinity). These PAM glider observations provide high-resolution, frequent, and localized in situ data for forcing numerical models and improving the understanding of the air-seaclimate system.