Decay Mechanisms of Near-Inertial Mixed Layer Oscillations in the Bay of Bengal BAY OF BENGAL : FROM MONSOONS TO MIXING

USAGE Permission is granted to copy this article for use in teaching and research. Republication, systematic reproduction, or collective redistribution of any portion of this article by photocopy machine, reposting, or other means is permitted only with the approval of The Oceanography Society. Send all correspondence to: info@tos.org or The Oceanography Society, PO Box 1931, Rockville, MD 20849-1931, USA. Oceanography THE OFFICIAL MAGAZINE OF THE OCEANOGRAPHY SOCIETY

of near-inertial ML currents have been of much interest in the last 30 years.Simmons and Alford (2012) discuss NIW generation and the propagation of the longest wavelengths, while Alford et al. (2016) review a broad range of topics and open questions, including the uncertain fate of short-wavelength NIWs.Our observations and analyses (a) show that groups of short-wavelength NIWs propagate rapidly into the interior and can remove energy from ML inertial oscillations; (b) suggest that mesoscale circulation aids the rapid downward propagation of short-wavelength NIWs; and (c) indicate that alignment of near-inertial current shear and wind leads to spikes in shear generation and presumably mixing and energy extraction from the ML.

Formation and Propagation of Near-Inertial Internal Waves
The classical mechanism for the generation and evolution of near-inertial motions can be explained by considering the response to a rapid increase in wind forcing on a horizontally uniform ocean.After the onset of winds, ML inertial oscillations begin.Due to Earth's rotation, the near-inertial velocity rotates clockwise in time in the Northern Hemisphere at the Coriolis frequency f = 2π/T i = 2Ω sin θ, where the inertial period is T i (= 43-39 hours from 16°N-18°N in our study area), Earth's rotation rate is Ω, and the latitude is θ (Ekman, 1905).Under broad storms in mid-latitudes with length scales reaching 1,000 km, flow converges because the inertial period changes with latitude (D' Asaro et al., 1995), which is referred to as the β effect.As the ML inertial currents converge/diverge, the base of the ML heaves downward/upward, which in turn produces heaving of isopycnals in the stratified interior and generates NIWs there (Gill, 1984).This process explains how long-wavelength NIWs are generated, which propagate far from their source regions, carry much of the energy, and are well reproduced by models (Simmons and Alford, 2012).
The small-scale structure, however, is not well described by this approach based on oceanic vertical normal modes (Gill, 1984;D' Asaro et al., 1995).These vertical modes arise from considering that vertical displacement of fluid particles in the water column is fixed to zero at the surface and bottom.Therefore, only certain wavelengths can be excited depending on the stratification and height of the water column in a finite depth ocean.This construction is analogous to the notes generated on a plucked string on a musical instrument.Alternatively, the propagation of NIWs in a field of mesoscale, geostrophic eddies can be treated as an initial value problem (Young and Ben Jelloul, 1997).The local rotation of the flow in the vicinity of strong fronts and eddies in Earth's rotating frame, also known as the relative vorticity, can be comparable in magnitude to the planetary vorticity (or rotation denoted by f as described in the previous paragraph).Both the relative vorticity associated with the mesoscale eddies and the change in planetary vorticity with latitude (β effect) can enhance the vertical transmission of near-inertial energy out of the ML (D' Asaro et al., 1995;Young and Ben Jelloul, 1997).
In this contribution, we focus on short-wavelength NIWs, which propagate vertically, display greater current shear in the vertical, and can contribute to local mixing.In the interior, typical NIWs have frequencies of 1-1.2 f, vertical wavelengths of 100-400 m, and horizontal wavelengths of order 10-100 km (Alford et al., 2016).In contrast, the lowest modes or longest waves have vertical wavelengths of twice the ocean depth and horizontal wavelengths of hundreds of kilometers.Typical vertical propagation at mid-latitudes is about 10 m day -1 (Alford et al., 2016).The downward group velocity is much faster in some of our observations in the BoB, which may be due to active mesoscale circulation.

INTRODUCTION Importance of Near-Inertial Internal Waves
Steady surface winds produce frictional transport to the right of the wind in the Northern Hemisphere, while an impulsive change in wind speed or direction excites oscillations at the local inertial period in the surface mixed layer (ML; Ekman, 1905).Steady winds produce vertical shear in horizontal currents (shear, hereafter) at the base of the ML and subsequently deepen the ML via mechanical mixing.However, with inertially rotating winds, the similarly rotating currents can propagate into the oceanic interior without much effect on ML depth (MLD; Dohan and Davis, 2011).Also, unidirectional winds with fluctuating magnitude excite such inertial motions (Weller et al., 2014;Majumder et al., 2015), which may have greater effects in the northern Bay of Bengal (BoB) due to its thin ML (Wijesekera et al., in press).The inertial currents heave the base of the ML to produce propagating near-inertial internal waves (NIWs) in the stratified ocean (Simmons and Alford, 2012).Longwavelength NIWs may propagate horizontally far from their source regions, while short-wavelength NIWs propagate vertically and may contribute to local mixing when shear overcomes the buoyancy restoring force in stratified waters or when strain (vertical stretching apart of isopycnals) reduces the stratification (Hebert and Moum, 1994;Alford and Gregg, 2001).Inertial and near-inertial currents often dominate the kinetic energy and contain most of the current shear (Alford et al., 2016).Thus, where and how this episodic near-inertial energy dissipates affects not only the upper ocean but the stratified interior as well.When global coupled atmosphere-ocean models include NIW-driven mixing, substantial changes occur in the MLD, tropical sea surface temperature, and precipitation (Jochum et al., 2013).
The mechanisms leading to the decay

Near-Inertial Internal Waves in the Bay of Bengal
The Ocean Mixing and Monsoon (OMM) program of India's National Monsoon Mission in the Ministry of Earth Sciences and the Air-Sea Interactions Regional Initiative (ASIRI) of the US Office of Naval Research catalyze a large collaboration between Indian and US scientists.
The aims are to understand small-scale processes affecting air-sea interaction in the BoB, to ultimately improve monsoon prediction on hourly to subseasonal time scales (Goswami et al., 2016, in this issue), and to train the next generation of scientists and build capacity (Tandon et al., 2016, in this issue).Climate models of the BoB are biased compared to observations, producing cold sea surface temperatures and deeper MLs (Fousiya et al., 2015;Parekh et al., 2015;Chowdary et al., 2016a), perhaps due to negative biases in surface heat flux (Chowdary et al., 2016a).Also, unrealistically weak biases in stratification allow stronger mixing, which, in turn, allow deep biases in MLD and cold biases in sea surface temperature (Fousiya et al., 2015;Parekh et al., 2015;Chowdary et al., 2016b, in this issue).These biases may be due to a lack of restratification by prevalent but unresolved submesoscale features in the BoB (Ferrari and Boccaletti, 2004;MacKinnon et al., 2016, in this issue).
In another coupled global model, NIW-related mixing improves tropical sea surface temperature and precipitation (Jochum et al., 2013).In the models mentioned above, we speculate that such mixing might exacerbate the biases and produce still deeper MLs and colder sea surface temperatures.On the other hand, we speculate further that if considerable near-inertial energy propagates quickly into the interior in the BoB, the biases could be reduced.The questions of what fraction of the ML near-inertial oscillations contributes to local mixing and what fraction propagates away as NIWs into the ocean interior are, therefore, important for coupled global models (Jochum et al., 2013).How does NIW generation in the BoB compare with that of mid-latitudes?At mid-latitudes, inertially rotating winds force inertial rotation in the ocean.In the BoB, winds can be unidirectional (see later discussion), but NIWs can also be generated by changing wind magnitude rather than inertial rotation (Weller et al., 2014;Majumder et al., 2015).The annual mean wind work and the near-inertial ML energy are almost as large in the BoB as in mid-latitude storm tracks (Chaigneau et al., 2008;Simmons and Alford, 2012).In mid-latitudes, MLs are deeper in winter and thus will project on longer vertical wavelengths, as observed by D' Asaro et al. (1995).Shortwavelength NIWs are expected in the BoB due to shallower MLs.
As noted above, decay of ML nearinertial motions and short-wavelength NIW propagation are not well understood in realistic ocean settings.Further complexity comes from the horizontal and vertical structure of the upper ocean in the BoB (Gordon et al., Hormann et al., MacKinnon et al., Sarkar et al., Shroyer et al., and Thangaprakash et al., 2016, all in this issue).For example, enhanced vertical fluxes are noted at fronts in the BoB and elsewhere (Johnston et al., 2011;Lucas et al., 2016, in this issue).In the BoB, the ML is thin (~10 m) and controlled by freshwater input from major river systems and precipitation (Figure 1; Goswami et al., and Gordon et al., 2016, both in this issue).Further below, a stratification maximum at ~50 m depth lies over the weakly stratified, deep, interior ocean.This transition layer between the ML and stratified interior mediates fluxes from the surface to the interior (Johnston and Rudnick, 2009;Dohan and Davis, 2011;Majumder et al., 2015) and inhibits vertical fluxes in the BoB (Lucas et al., and Sarma et al., 2016, both in this issue).If shear (S) is comparable to stratification (commonly measured with the frequency of vertical oscillations in a stratified fluid or the buoyancy frequency, N ), then mixing is possible when their squared ratio or Richardson number, denoted Ri = N 2 /S 2 , is less than 0.25.Ri is often order one in the transition layer, suggestive of instability.Both pre-existing shear from remotely generated NIWs in the transition layer and downward propagation of locally generated NIWs contribute to the near-inertial ML energy budget (Plueddemann and Farrar, 2006).Wind and near-inertial shear alignment is important for shear generation, subsequent mixing, and, thus, the ML energy budget (Brannigan et al., 2013;Majumder et al., 2015).

Mesoscale Effects on Vertical Energy Propagation
NIW vertical energy flux depends on the product of the energy density of the waves and their vertical group speed, c gz ≈ N 2 α 3 λ x /(2πω), where the aspect ratio of vertical to horizontal wavelengths is α = λ z /λ x and the wave's radial frequency is ω (Gill, 1982).The aspect ratio depends sensitively on the difference between wave frequency and inertial frequency.The horizontal length scales of ML inertial motions are imposed by atmospheric forcing, which has a length scale of order 1,000 km.To reduce their horizontal wavelength or horizontal coherence, oceanic inertial motions in the ML must be disrupted to allow vertical propagation into the interior (D' Asaro et al., 1995).The β effect explains the decay of surface drifters' inertial motions (Park et al., 2009), but neither observations of the predicted rapid decay e-folding scale of less than four days equatorward of ~20°N nor subsurface data are noted in Park et al. (2009).
Eddies can influence vertical propagation in at least five ways.(While we discuss eddies in this paragraph, identical processes occur at fronts, current jets, or other meso-and submesoscale features.)(1) They may bring stratified water closer to the base of the ML.
(2) A key feature is their relative vorticity or, in other words, the local rotation of a fluid parcel (Weller, 1982;Kunze, 1985).The vertical component of eddy vorticity can either be in the same (counterclockwise, positive, cyclonic) or opposite (clockwise, negative, anticyclonic) direction of Earth's rotation.Thus, the effective inertial frequency decreases when the eddies rotate clockwise, which in turn increases the aspect ratio of the waves and increases their vertical group speed.Internal waves may also become trapped in regions of negative vorticity.
(3) As the mesoscale flow converges, the rotational radius of near-inertial motions decreases.Also, the near-inertial current amplitude then increases to conserve angular momentum (Weller, 1982).As the radius/wavelength of near-inertial motions decreases, vertical group speed increases.( 4) As nearinertial motions encounter a field of mesoscale eddies of smaller length scale, the near-inertial motions sample a number of eddies and lose coherence proportional to the kinetic energy in the spatially varying flow field (Young and Ben Jelloul, 1997).This process reduces the length scale of the near-inertial motions.
(5) The shear of the background currents causes internal waves to refract and reduce their vertical wavelengths (Kunze, 1985).When the waves' propagation speed is equal and opposite to the current speed, a critical level is reached, at which the energy cannot propagate vertically any further, and the wave breaks (Munk, 1981).
In summary, the vertical group speed can be increased by: (1) reducing horizontal wavelength from mesoscale currents dephasing inertial motions (Young and Ben Jelloul, 1997), from flow convergence (Weller, 1982), or from the stronger β effect (D' Asaro et al., 1995;Park et al., 2009); and (2) lowering the effective inertial frequency due to the relative vorticity of anticyclones or clockwise eddies and making the NIWs effectively more superinertial (Weller, 1982;Kunze, 1985;Elipot et al., 2010).Also, vertical propagation can be augmented at fronts or eddies, where stratified water is brought closer to the ML.In anticyclones, internal waves may become trapped and amplified (Munk, 1981;Kunze, 1985).

METHODS Mooring
The National Institute of Ocean Technology (NIOT, Chennai, India) runs the Ocean Moored Buoy Network for the Northern Indian Ocean (OMNI), which was set up to investigate the strong coupling between monsoons and cyclones within the northern BoB (Venkatesan et al., 2013).The program has deployed moorings with meteorological and subsurface sensors at six different locations in the BoB.Here, we present an analysis of the data at hourly resolution acquired from August to September 2013 at mooring BD09, located at 17.88°N, 89.67°E (Figure 1a).The meteorological sensors are at 3 m above sea level and measure air temperature, surface pressure, relative humidity, downwelling long-and shortwave radiation, and wind speed and direction.Subsurface conductivity and temperature sensors are at 10 different depths: 5,10,15,20,30,50,75,100,200, and 500 m.The mooring also includes one current meter at 1.2 m depth and one downward-looking acoustic Doppler current profiler (ADCP) at 7.5 m, which measures ocean current from 10-105 m at 5 m resolution.

Float
A SOLO-II profiling float (serial number 8119; Argo platform number 5904302) was deployed at a front near 16.5°N, 85.5°E in the BoB in conjunction with  ship surveys by R/V Roger Revelle in 2013 (Figure 1c).The float profiled from 0-250 m and back every 1.7 hours.After completing 237 profiles, the float entered the typical Argo mode with a profile every 10 days, which does not resolve the inertial frequency.The rapid profiling lasted from November 21 to December 7, 2013.
The float has a pumped Sea-Bird Electronics (SBE) 41CP conductivitytemperature-depth (CTD) instrument from which potential temperature (θ ), salinity (S), in situ density (ρ), and potential density (σ θ ) are obtained on ascents, when the CTD encounters clean flow.Float data are averaged in 10 m bins from 10-250 m depth.The float moves upward by inflating an external bladder to gain buoyancy.Downward motion is achieved by shifting mineral oil back into the pressure case and thereby deflating the bladder and decreasing the vehicle's buoyancy (Davis et al., 2001).A float is assumed to drift with the current, and the depth-mean current over a dive cycle is obtained from GPS fixes at consecutive surface intervals.
Signals with frequencies higher than the profiling frequency are smeared and aliased and cannot be subsequently disentangled.Due to its slow motion of ~20 km day -1 , the drifting float smears this high-frequency temporal variability into spatial variability (Figure 1c; Rudnick and Cole, 2011).Aliasing arises because the profiling frequency is in the internal wave band.By using the analysis method described below, the effects of smearing and aliasing on near-inertial motions are mitigated.
For analysis of NIWs, the data from a slowly moving platform are analyzed as slowly moving time series (Johnston et al., 2015).Since the platform is moving and NIWs are episodic, the expected signal is well suited to a wavelet analysis, which resolves variance as functions of time and wave period.The analysis is done by convolving the time series with a wavelet, a compact function in time and frequency.(Note the contrast with Fourier spectral analysis using sine waves, which may produce precise frequency lines, but has no resolution in time.)To analyze the expected signal of NIW groups, we use the Morlet wavelet, which has a similar character-namely, a sinusoidal oscillation within a Gaussian envelope.Our averaging window for inertial motions is determined from this envelope, which has a full width at the e-folding level of 2 √ -T i or about five days.Near-inertial and semidiurnal motions are band passed within ±25% of the central frequency using wavelets.Mesoscale motions are obtained from the float via a three-day low-pass filter using wavelets.Since NIWs have flat aspect ratios (long horizontal and short vertical wavelengths), they are better detected with velocity rather than in isopycnal displacements.However, both shear (vertical gradient of horizontal currents) and strain (vertical gradient of isopycnal displacements) are implicated in mixing (Alford and Gregg, 2001).The float does not measure depth-dependent currents and so to obtain information on the vertical propagation of NIWs, we use the band-passed potential density fluctuations.
To help understand the in situ float data, surface winds are obtained from a global data assimilation by the National Center for Environmental Prediction (Kalnay et al., 1996).

Model
We use the Regional Ocean Modeling System (ROMS) to understand the submesoscale and mesoscale structure during the float mission (Figure 1c).To allow the model to evolve freely, it was initialized and given lateral boundary conditions from an oceanic state estimate and run with a reanalysis atmospheric state and a prognostic ocean for 20 days at a time.The purpose of this setup is to limit model drift, but at the same time allow the ocean to evolve without constraints.ROMS is a free-surface, hydrostatic primitive equation ocean circulation model.Horizontal resolution is 3.5 × 3.5 km (1/32° Mercator) with 50 sigma layers.We use the K-profile parameterization (KPP) turbulence closure scheme (Large et al., 1994) in the surface and bottom boundary layer and Richardson number-dependent mixing in the ocean interior.The model is initialized from the HYbrid Coordinate Ocean Model (HYCOM) GLBa0.08 1/12° Navy Coupled Ocean Data Assimilation (NCODA) system (Metzger et al., 2014), and lateral boundary conditions are a radiation/nudging scheme with radiation on outflow and nudging on inflow to HYCOM/NCODA, with a three-day relaxation time scale (Marchesiello et al., 2001).The vertical grid is 50 sigma layers, which in the absence of topography appears as depth coordinates in the upper ocean with levels at 0, 5, 10, 15, 25, 30, 40, 50 m, and so on.Flux forcing is computed using the model state and atmospheric state updated every three hours, with turbulent fluxes calculated using bulk formulae (Large and Pond, 1981).Tidal forcing is implemented through the open boundary conditions with semidiurnal (M 2 and S 2 constituents) and diurnal (O 1 and K 1 constituents) tidal velocities and sea level from TPXO7.2 (Egbert and Erofeeva, 2002).River discharge is from the Dai and Trenberth (2002) monthly river climatology.Horizontal tracer mixing is via a harmonic eddy viscosity coefficient of 2.0 m 2 s -1 and momentum is via a harmonic eddy viscosity coefficient of 10.0 m 2 s -1 .A linear bottom drag of 3 × 10 -4 m s -1 and a nonlinear drag coefficient of 3 × 10 -3 are used.

Enhanced Shear Generation at the Mixed Layer Base
Figure 2 presents an analysis of the data collected at the BD09 mooring from August to September 2013.The wind stress magnitude changes with time (Figure 2a,b) and excites inertial currents in the ML (Figure 2c).The MLD (red line; Figure 2b) corresponds closely with the wind stress magnitude (Figures 2a and 3a).MLD is defined as the depth at which density exceeds the surface value (here z = 5 m is the shallowest measurement) by 0.03 kg m −3 .The mean MLD is 15 m and is modulated by the active-break spell of the summer monsoon winds.After August 28, the ML is relatively shallow with a MLD hovering around 5 m due to very fresh water at the ocean surface and the consequent strong salinity stratification (Figure 2d,e).
The MLD and the strength of stratification below the ML are important parameters that determine the extent to which the ML inertial currents can propagate into the deeper ocean.On August 6-7, the MLD approaches the depth of the stratification maximum (Figure 2e), which may be related to the downward NIW group propagation on August 9-11 below 40 m (note the upward phase propagation or wave crest orientation in time; Figure 2c).
In September, even when both ML inertial currents are large and the depth difference between the ML and the stratification maximum is small, it is difficult to relate bursts of near-inertial currents in the interior to surface currents.Much of the shear remains at the base of the ML.
To better understand the near-inertial motions, we examine the vertical structure of the base of the ML and the time evolution of shear there.This interface between the stratified interior and unstratified surface layer of the ocean is known as the transition layer (TL; Ferrari and Boccaletti, 2004;Johnston and Rudnick, 2009).We calculate the transition layer thickness (TLT) as the depth of the shear maximum (D) below the MLD: TLT = D − MLD.Similar calculations can be made with the stratification maximum.Our estimates of TLT are 3-40 m, which are consistent with a previous study in the Arabian Sea (Majumder et al., 2015).Due to the wider vertical spacing of temperature and salinity measurements at the mooring, we use TLT based on shear for the rest of our calculations because the ADCP bin sizes are 5 m.The TL is the region between the seven-day low-passed lines in Figure 2f.On August 6-7, when the TL thins for two days (TLT is difference between MLD in gray, which in not low pass filtered, and the depth of the stratification maximum; Figure 2e), there is some downward near-inertial energy propagation.Otherwise, until August 15, the TL is thick and, despite some inertial motions in the ML, little energy propagates downward (Figure 2c).The TL acts as a barrier.After August 22, the stratification maximum near 40 m depth disappears and is replaced by the shallow stratification maximum due to fresher surface water (Figure 2e).Then, when wind changes excite ML near-inertial motions, there is downward propagation beneath the ML (Figure 2b,c).ML near-inertial motions decay in about three to seven days (Figure 2c).For the idealized scenario of an infinitely deep ocean without background flow and with an abrupt jump in the stratification at the base of the ML from zero in the ML to a finite value in the ocean interior, Moehlis and Llewellyn-Smith (2001) gave an analytical estimate for the time scale over which the near-inertial energy in the ML decays via downward propagation of NIWs.Specifically, they calculate that over a time of T d and 2T d , 50% and 82% of the ML energy is transferred to regions below the ML via downward propagation of NIWs.Here, T d is given by , where H is the MLD, and N 2 the uniform deep ocean stratification.Typical BoB values are: f = 4.47 × 10 -5 s -1 , β = 2.17 × 10 -11 s -1 m -1 , and N = 0.02 s -1 .For MLDs of 10 and 50 m, T d = 31 and 11 days, respectively, thus insufficiently explaining the decay time scales found in our measurements.We note, however, that either a sheared background flow or an equivalent β effect via the presence of mesoscale eddies can affect the decay time scales (see Introduction).
Next, the time evolution of bulk shear of the TL is estimated from the mooring data.We use a simple two-layer model of the upper ocean forced only by wind for the calculation (Brannigan et al., 2013;Majumder et al., 2015).The model approximates (1) the ML as a slab that accelerates due to wind and inter facial stress, and (2) the lower layer as quiescent.Given the observed wind stress, MLD, and the TLT, the model predicts current shear between the upper and lower layers.The shear is the velocity difference between the layers divided by the TLT.Shear generation is maximized (1) for thin ML and TL, and (2) when the wind stress, which varies relatively little in direction, aligns with the rotating inertial currents.Multiple events of mean bulk-shearsquared with amplitudes varying from 0.1-2 × 10 -4 s -2 are present during the summer monsoon in 2013 (Figure 3b).The shear magnitude increases in the presence of the thin freshwater layer during September.The shear generation term, P(S 2 ), roughly follows the time evolution term (∂S 2 /∂t; Figure 3c).When the phases of the bulk shear and the relative stress align, shear generation occurs (0° = direction of relative shear in Figure 3d).Note the main shear generation event is September 16-18 when this alignment occurs and is coincident with a wind event.Despite somewhat larger wind stresses in mid-August, the lack of alignment leads to little shear generation.
Since there are other sources of observed shear besides near-inertial currents, the time series are band passed around the inertial frequency (Figure 4).In the TL, there is little near-inertial shear in early August due to weaker winds and thicker TL.However, from late August onward, the ML shoals (Figure 2).ML inertial currents and near-inertial shear in the TL are elevated in four episodes centered on August 26 and September 3, 11, and 23 (Figures 2c and 4a).During and after the wind event on September 19 (Figure 2a),   x 10 -8 x 10 -9 x 10 -5 separately resolves counterclockwise and clockwise rotation.This particular wind shift projects onto a clockwise-rotating inertial signal from November 25-28 (Figure 5b).The magnitude of the inertial signal in the wind is about 1 m s -1 (blue line; Figure 5f), which is a fraction of the total wind speed of over 10 m s -1 (Figure 5a).Until November 25, 2013, the winds have a diurnal component, which weakens and then grows again on November 29 (Figure 5b).
Initially, the depth-mean current oscillations are semidiurnal (Figure 5c).Also on November 27, a distinct change to near-inertial oscillations in currents is visible without any filtering (Figure 5c).The vector components of the current are denoted u and v, which are positive eastward and northward.The nearinertial currents decay by a factor of two by December 4.
NIW currents rotate in the clockwise direction due to Earth's rotation.Thus, the wind shift's projection onto clockwise-rotating winds favors the shear generation and time evolution agree well (Figure 4b).Agreement is also good under moderate winds from the beginning of the record until August 18.From August 16-23, the ML deepens under strong winds (Figure 2a,b), which indicates an energy sink for the shear that is being produced.From August 25 to September 7, shear generation is larger than shear evolution, which suggests energy is being supplied to NIWs and dissipating in the TL (Figure 4b).This is consistent with a thin ML and a TL thickening (August 31 to September 7; Figure 2f).These results, along with previous work by Brannigan et al. (2013) and Majumder et al. (2015), suggest this twolayer model is a useful tool for evaluating the near-inertial ML energy budget.

Rapid Downward NIW Propagation at a Front
The float was released from R/V Revelle at a freshwater front in the northern BoB and advected westward in a looping trajectory, which is a combination of mean flow, semidiurnal tide (period is 0.5 day), and near-inertial oscillations (period is 1.7 days; Figures 1c and 5c).Here, we examine how these inertial motions decay in the surface ML and propagate downward into the ocean interior.Based on observations from this float and two other glider missions, clear and rapid downward propagation is found near fronts, where the distance between the ML and the stratified interior is reduced.Otherwise, when wind events occur, the inertial motions are concentrated in the near-surface fresh layer and do not cross the barrier layer beneath it.
On November 27, there is a sudden wind shift (Figure 5a).To determine the frequency content of the vector wind time series, we examine the scalogram obtained from the wavelet analysis.The scalogram is analogous to a power spectrum from Fourier analysis, but displays variance as a function of time and the oscillation's period (Figure 5b).Furthermore, a wavelet transform of the vector winds or currents excitation of clockwise currents.The currents and winds in the near-inertial and semidiurnal frequency bands are reconstructed from the wavelet transform by using wavelet coefficients from a pass band within 25% of each frequency.A near-inertial peak (ω = 1.1 f ) is found from November 28 to December 3 (Figure 5e,f).Semidiurnal currents are initially present, abruptly disappear during the wind shift from November 27 to December 2, and then reappear (Figure 5e,f).The reason for the abrupt change in semidiurnal currents is unclear, but could be related to the effect of stratification on semidiurnal internal wave propagation, to the spring neap cycle, or to the position of the float with respect to the front.
The float encounters water on either side of the front multiple times due to surface-intensified inertial motions, which move surface waters more with respect to the float as it profiles to 250 m.With southward near-inertial currents on November 30 (Figure 5f), for example, the freshwater front is advected southward (Figure 6b), the mixed layer deepens (black line; Figure 6a) and becomes less stratified (Figure 6a), and denser water is moved over the float (Figure 6b).This description is qualitatively consistent with the model, where the observed float trajectory is along a modeled freshwater front (Figure 1c).
Since the velocity measurements are a depth mean, we use near-inertial density fluctuations to understand the depth-dependent NIW propagation into the interior (Figure 6b-d).An example scalogram from 150 m depth has similar structure to that of the depth-mean currents-namely, semidiurnal fluctuations are interrupted on November 28 and resume two days later, and nearinertial fluctuations are maximum from November 30 to December 3 (Figure 6c).From similar scalograms, density fluctuations in the near-inertial frequency band are obtained from a wavelet band pass at each depth from the surface to 250 m.To visualize the fluctuations better, they are scaled by a factor mainly related to the stratification, and the plotted quantity is energy, once it is squared (Figure 6d).
The fluctuations are largest in the ML and then propagate downward reaching 250 m on December 2, which is within four to six days (or two to three inertial periods) of the wind shift.This downward energy propagation of about 30 m day -1 is rapid compared to other available measurements.The energy propagates downward in a group of several waves, which is expected for a forcing event of limited duration.A feature of internal waves with downward energy or group propagation is that phase propagation is upward (i.e., wave crests slope upward with time).Also, these observations provide an estimate of the vertical wavelength of the waves in the interior.The wavelength changes with stratification, but on December 4 it appears to be greater than 200 m and probably around 300 m.
The internal wave dispersion relation relates wave frequency to other quantities, such as wavelength.We use the observed vertical wavelength and vertical group speed to establish limits for the horizontal wavelength (Figure 7) as follows.
(1) The observed vertical group speed is 33 m day -1 , but lower and upper limits of 20 m day -1 and 40 m day -1 are not unreasonable (Figure 6d).( 2) The observed vertical wavelength is over 200 m and likely at least 300 m (Figure 6d).(3) We use the observed frequencies from nearinertial velocity and density fluctuations as lower and upper limits on frequencies (ω = 1.1 f and 1.2 f ).Combining these three limits produces a range of possible consistent values (yellow stippled area, Figure 7).However, we emphasize that as the waves approach the inertial frequency and the horizontal wavelength increases, wave properties become increasingly sensitive to wave frequency.Based on a frequency of 1.1 f, a vertical group velocity FIGURE 6.Since the float was released at a front, near-inertial oscillations move (a) lessstratified and (b) denser upperocean water from the north (i.e., southward current in Figure 5f) over the float as it profiles deeper.Panel (a) shows the mixed layer depth (black line) and the depth of maximum stratification (white line).(c) A scalogram shows variance in density fluctuations at semidiurnal and inertial periods.(d) Near-inertial density fluctuations from a band pass (within 25% of the inertial frequency as in Figure 5f) are scaled to equal potential energy when squared.Downward energy propagation to 250 m in six days is found, while mixed layer fluctuations decay.
of 33 m day -1 , and a vertical wavelength of 300 m, the horizontal wavelength is about 300 km.This value is comparable to the spacing of the adjacent cyclonic (clockwise) and anticyclonic (counterclockwise) eddies, which lie south and north of the float's trajectory (Figure 1c).

Mixed Layer Energy Decay Estimate
At the same time as the NIWs propagate downward in the interior ocean, the near-inertial density fluctuations in the ML decay substantially by December 4 (Figure 6d) as do the depth-mean currents (Figure 5c,f).Thus, the observed decay scale of the currents and density fluctuations is about seven days.Here, we estimate the ML energy removal by the downward-propagating NIW group.First, the ML kinetic energy is estimated.The depth-mean currents over the float's dive profile are obtained from adjacent position fixes during surface intervals.However, since these currents are initiated by a sudden wind event, they are concentrated in the 10-50 m deep ML (black line denotes the MLD in Figure 6a and is from a 0.1 kg m -3 density difference from the surface).To estimate the near-inertial currents in the ML, we use the observed near-inertial depthmean velocity (0.08 m s -1 , Figure 5f) and assume near-inertial fluctuations over the profiling range beneath the ML average to zero.For a 50 m deep ML, this yields an ML near-inertial current of 0.4 m s -1 .This corresponds to an ML kinetic energy of about 4 kJ m -2 .This value decreases if near-inertial fluctuations contribute to a non-zero depth-mean current over the profiling range.
Next, we estimate the downward energy flux.For this, we need the frequency of the NIWs seen in density fluctuations (ω = 1.2 f from a wavelet analysis; Figure 6c), the kinetic to potential energy ratio of these waves at this frequency (about 6), and the observed potential energy (about 1 J m -3 ; Figure 6d).The downward energy flux is then estimated as 2.5 mW m -2 from the product of the observed downward group velocity of 33 m day -1 and the total energy density of about 7 J m -3 , which is obtained from the sum of the observed potential energy and the estimated kinetic energy (potential energy times six, as noted above).
The ML energy decay time scale is estimated as the ML kinetic energy divided by the downward energy flux, yielding a value of about 19 days.We emphasize that this calculation is approximate, with large uncertainty, because (1) the depthvarying currents were not actually measured, and (2) the kinetic to potential energy ratio is especially sensitive to the wave frequency near the inertial frequency.If the frequency of the NIWs is instead 1.1 f (as seen for the depthmean currents; Figure 5e), then the kinetic to potential energy ratio almost doubles, the vertical energy flux estimate almost doubles, and the decay scale decreases to 11 days.Given the observed ML near-inertial decay scale of seven days, the downward-propagating NIW group may remove about half of the ML kinetic energy.
The remaining near-inertial energy may either propagate away as longwavelength waves not detected by our methods or dissipate via turbulence in the highly sheared and stratified transition layer (region between black lines, Figure 6d).The change of near-inertial density fluctuations with depth in the transition layer implies current shear of about 0.2 m s -1 over 10 m (or S 2 = 4 × 10 -4 s -2 ), which depends on the kinetic to potential energy ratio.Since the Richardson number is about 1 (N 2 in the transition layer is about 10 -4 s -2 ; Figure 6a), this result suggests shear instability and turbulence are possible.
In summary, our interpretation of the data is that an impulsive wind event generates ML inertial motions, and a considerable fraction of the energy radiates rapidly downward in a group of short-wavelength NIWs near a density front.Some caveats are noted.The short vertical wavelength motions are likely well captured due to their slow horizontal propagation speeds, but waves with long vertical wavelengths have rapid horizontal propagation speeds and move rapidly out of the study area.On the other hand, long waves may propagate into the study area.Such waves could impact the ML near-inertial energy budget.No attempt has been made to identify such waves because of the float's limited profiling depth.However, the float traversed about 3° of longitude, moving with the depth-mean flow in the upper ocean, and does not provide a time series at a single point.This distance is comparable to the FIGURE 7. The internal wave dispersion relation is used to calculate the vertical (a) group velocity and (b) wavelength as functions of frequency and horizontal wavelengths.The horizontal wavelength is estimated by using observed limits on the vertical group velocity, the vertical wavelength, and the frequency-the yellow stippled area with frequencies of 1.1-1.2f , group velocities of 20-40 m day -1 (black lines), and vertical wavelengths of 200-300 m (gray lines).estimated horizontal wavelength of the near-inertial waves but is much smaller than the scale of the atmospheric forcing.Thus, the float placement with respect to atmospheric forcing is not special.However, the float moved along a front sampling a coherent short-wavelength near-inertial signal.

SUMMARY
Shear at the base of the ML contributes to mixing, which affects the heat content and sea surface temperature of the ML and thus precipitation.Such effects may be more pronounced with a thin ML as in the BoB.Thus, how energy is removed from inertial motions, by either downward NIW energy propagation or shear generation at the base of the ML, is consequential.Rapid downward propagation of short vertical wavelength NIW groups occurs near fronts, as observed by a float described here (and two additional glider missions, which were omitted for brevity).The NIW group propagates down to 200 m depth in about six days, yielding a vertical group speed of 30 m day -1 .The NIW vertical energy flux is sufficient to remove energy from the ML inertial motions in 11-19 days and is sensitive to the wave frequency.The observed ML inertial decay is about seven days, which suggests about half of the observed decay is due to the downward-propagating NIW group.No attempt has been made to estimate the contribution of waves with long wavelengths, which may radiate energy in the horizontal.
Elevated shear at the base of the ML likely dissipates some of the remaining energy through shear instability.Mooring data show shear generation is maximized beneath the ML in a two-layer model of the upper ocean, (1) when near-inertial currents align with the wind stress, and (2) for thin ML and transition layers.

FIGURE 1 .
FIGURE 1. Sea surface salinity from the Argo climatology shows (a) large freshwater input to the northern part of the Bay of Bengal in August, and (b) extension of the freshwater down the boundaries of the bay in November.The position of mooring BD09 of the Ocean Moored Buoy Network for the Northern Indian Ocean (OMNI) program is indicated with a plus sign in (a).The broad picture from Argo is complemented by model results on December 1, 2013, averaged over the upper 23 m, which shows (c) freshwater and salty water strained into narrow filaments by cyclonic (counterclockwise) and anticyclonic (clockwise) eddies south and north of the float trajectory (line).The float was released on November 21, 2013 (white circle), at the front and completed rapid profiling on December 7, 2013 (black circle).The position of the float is indicated on December 1, 2013, by a plus sign.Current vectors are lowpassed beyond three days using wavelets.A scale vector of 0.5 m s -1 is plotted in the upper left.

FIGURE 2 .
FIGURE 2. Analysis of the data obtained during August to September 2013 from OMNI mooring BD09.(a) Wind stress as a function of time (plotted as a vector with zonal and meridional components).(b) Temporal rate of change of zonal (blue) and meridional (green) components of wind stress as a function of time; the red curve shows the variation of the mixed layer depth with time.(c) Net magnitude of the near-inertial currents (obtained after filtering the zonal and meridional velocity components around the local inertial frequency) plotted as a function of time and depth.(d) Salinity plotted as a function of time and depth.(e) Square of the buoyancy frequency plotted as a function of time and depth.(f) Seven-day low-passed mixed layer depth (red) and maximum shear-squared depth (green).The gray curves in (c) and (e) and the black curve in (d) show the location of the bottom of the mixed layer.

FIGURE 4 .
FIGURE 4. (a) Near-inertial shear in the transition layer is obtained by band passing the total shear (Figure3b).(b) The near-inertial shear evolution (blue) and near-inertial shear generation (red) are in better agreement than the unfiltered terms in Figure3c.

FIGURE 5 .
FIGURE 5. (a) Wind vectors show a sudden shift on November 27.(b) The shift projects onto a clockwise near-inertial rotation in a scalogram, which displays wind variance (or wavelet energy density, WED) as a function of time and oscillation period.T i and T D2 denote the inertial and semidiurnal periods.(c) Depth-mean currents also show a shift from semidiurnal to near-inertial currents after the wind shift.(d) Current vector variance is shown in scalograms for counterclockwise and clockwise (e) rotation.Clockwise variance is dominant due to internal waves with a sudden onset of nearinertial currents on November 27.(f) Northward currents and winds are band passed within 25% of the semidiurnal and inertial frequencies (compare to Figure 5c).Gray, red, and blue lines indicate semidiurnal current, near-inertial current, and near-inertial wind, respectively.