AMERICAN METEOROLOGICAL SOCIETY - AMS Journals

3 downloads 124 Views 6MB Size Report
derived surface currents, that in Bodega Bay (northern California) the ..... observed close to Huntington Beach, California (~33.63° N) by Nam and Send (2013) ...
 ET AL. MIHANOVIC

NOVEMBER 2016

3487

Diurnal Sea Breezes Force Near-Inertial Waves along Rottnest Continental Shelf, Southwestern Australia HRVOJE MIHANOVIC´ School of Civil, Environmental and Mining Engineering, and The University of Western Australia Oceans Institute, The University of Western Australia, Crawley, Western Australia, Australia, and Institute of Oceanography and Fisheries, Split, Croatia

CHARITHA PATTIARATCHI AND FLORENCE VERSPECHT School of Civil, Environmental and Mining Engineering, and The University of Western Australia Oceans Institute, The University of Western Australia, Crawley, Western Australia, Australia (Manuscript received 21 January 2016, in final form 17 August 2016) ABSTRACT Observations of upper-ocean dynamics close to the critical latitude (ratio of the local inertial to diurnal frequency is 0.94) from a range of platforms (surface currents using HF radar, moored instruments, and satellite remote sensing data) off southwest Australia indicated the presence of energetic, near-inertial waves generated through the diurnal–inertial resonance. During the austral summer, when southerly winds and land–sea breeze (LSB) system dominated the wind regime, strong counterclockwise diurnal motions (amplitudes surpassing 0.3 m s21) penetrated to 300-m depth with diurnal vertical isotherm fluctuations up to 60 m. The upward phase propagation speed of ;140 m day21, deep penetration of diurnal currents below the mixed layer, and the ;1808 phase difference between the upper and lower water column suggested that the local LSB system caused the resonant diurnal motions. Relative vorticity fluctuations along two cross-shore transects indicated changes to the local effective Coriolis frequency by more than 50% (60.5f ). In the presence of strong and relatively consistent cross-shore diurnal wind forcing in the study area, the main factors that controlled the observed energetic but sporadic near-inertial oscillations were the Leeuwin Current strength and spatial–temporal variations. These variations controlled the effective Coriolis frequency and enabled the effective pumping of diurnal wind energy into the ocean particularly when the effective Coriolis frequency was ;24 h.

1. Introduction Interaction between the atmosphere and ocean is of fundamental importance for understanding the dynamics and mixing in the upper ocean. In combination with the Coriolis force, surface currents generated by wind stress have a maximum response at the local, inertial frequency (Ekman 1905; Simpson et al. 2002; Kim and Crawford 2014). A resonance condition occurs when the period of wind forcing is close to the local, inertial frequency and occurs in coastal regions close to the critical latitude (308N or 308S 6 108). These resonant wind current

Denotes Open Access content.

Corresponding author address: C. Pattiaratchi, The University of Western Australia, 35 Stirling Highway, Crawley, WA 6009, Australia. E-mail: [email protected] DOI: 10.1175/JPO-D-16-0022.1 Ó 2016 American Meteorological Society

responses close to the critical latitude have been addressed in a variety of field and theoretical studies and have been shown to enhance the upper-ocean velocity field and vertical mixing (e.g., Simpson et al. 2002; Hyder et al. 2002, 2011; Zhang et al. 2009; Kim and Crawford 2014). In this paper, we present data obtained from a range of platforms: surface currents using HF radar (HFR), from moored instruments (bottom-mounted ADCP and thermistor strings), and remotely sensed sea surface temperature data from the Rottnest continental shelf, southwest Australia, to examine the influence of the strong land–sea breeze (LSB) cycle on the upperocean dynamics through the diurnal–inertial resonance in a region close to the critical latitude, where the local inertial period is 22.6 h. In microtidal coastal regions, where tidal forcing is negligible, vertical mixing is strongly influenced by wind events such as extratropical storms (e.g., Ferré et al. 2005; Verspecht and Pattiaratchi 2010; Hopkins et al.

3488

JOURNAL OF PHYSICAL OCEANOGRAPHY

2014), tropical cyclones/hurricanes (e.g., Jacob et al. 2000; Zhang et al. 2014), and strong sea breezes (e.g., Pattiaratchi et al. 1997; Gallop et al. 2012). In deep waters off the continental shelf, vertical mixing contributes to the transport of nutrients across the pycnocline to the primary producers at the surface layer (Sharples and Tett 1994; Sharples et al. 2001; Burchard and Rippeth 2009). In strongly stratified systems, the pycnocline prevents diffusion of nutrients, oxygen, and phytoplankton, and in the absence of vertical mixing, this can cause the detrimental depletion of nutrients in the surface boundary layer (e.g., Sharples et al. 2001). In oligotrophic environments, processes that are able to transport nutrients from below the pycnocline to the upper ocean are important for the local primary production. Wind-induced resonance occurs when either the broadbanded, open-ocean winds excite the near-inertial frequency more effectively than other frequencies or when the main wind forcing is at a single frequency, which causes the ocean to respond at the same frequency. The latter situation is proposed to occur off southwest Australia, where the diurnal LSB forcing is close to the local inertial frequency. Resonance of this kind has also been observed at the critical latitude off Namibia at 28.68S (Simpson et al. 2002); the Thermaikos Gulf, northwest Aegean Sea (Hyder et al. 2002); the Texas–Louisiana shelf (Chen et al. 1996; DiMarco et al. 2000; Zhang et al. 2009, 2010; Gough et al. 2016); and at other critical latitudes around the world (Rosenfeld 1988; Lerczak et al. 2001; Hunter et al. 2007; Edwards 2008). Kim and Crawford (2014) used HFR-derived surface current data obtained along the west coast of the United States to assess the relation between LSB forced currents and latitudes between 328 and 478N and concluded that the broadbanded, wind-driven, near-inertial currents at locations closer to critical latitude were 8 to 12 times stronger than at locations farther away from the critical latitude. Alongshore coastal wind forcing at a frequency slightly lower than the inertial frequency (the ‘‘subinertial’’ case) causes strong pycnocline oscillations dominated by coastal-trapped waves (Clancy et al. 1979; Wolanski 1986; Rosenfeld 1988; Cushman-Roisin 1994). In contrast, for the superinertial case (forcing frequency is higher than the inertial frequency) waves mostly radiate away from the coast during the periodic coastal upwelling, and the cross-shore wind produces more pycnocline oscillations than the alongshore wind (Orlic´ and Pasaric´ 2011). Inertial currents oscillating in an unbounded ocean display no vertical variability, but when a coastal boundary is introduced, the zero-flux condition at the boundary changes the behavior of the inertial currents, resulting in a phase shift between the

VOLUME 46

surface and bottom layers (Chen et al. 1996). The zeroflux condition induces an opposing cross-shelf sea surface slope at the coast, which acts as a pressure gradient through the water column. The wind forcing acts only at the surface, which causes a cross-shelf flow with a 1808 phase shift at the interface between the upper mixed layer and the lower stratified layer to drive a reverse flow in the lower stratified layer (Craig 1989; Simpson et al. 2002). The different phases in the surface and bottom layer motion produce vertical shear, which enhances vertical water column mixing and turbulent dissipation (Simpson et al. 2002). Low-frequency background currents can also change the local relative vorticity and local effective Coriolis frequency (Kunze 1985; Lerczak et al. 2001). Kunze (1985) derived a dispersion relation for near-inertial internal waves propagating in geostrophic current shear and concluded that local relative vorticity z shifted the internal wave band’s lower bound from the planetary Coriolis frequency to the effective Coriolis frequency: feff 5 f 1 z/2. For the Southern Hemisphere, the effective frequency is shifted below (above) the Coriolis frequency in regions where low-frequency motions have positive (negative) vorticity. Counterclockwise (anticyclonic) eddies in the Southern Hemisphere have positive vorticity and therefore decrease the local Coriolis frequency, whereas clockwise (cyclonic) eddies increase the local Coriolis frequency. This paper uses data from different platforms (HF radar, moored instruments, and satellite remote sensing) to examine the generation and influence of nearinertial oscillations along the Rottnest continental shelf, southwest Australia. In particular the following are examined: 1) the forcing, response, and the extent of the diurnal variability of the winds and the surface current field; 2) the vertical structure and horizontal extent of the rotatory motion; 3) the influence of the wind forcing and horizontal current shear on near-inertial motions and isotherm fluctuations in the vertical; and 4) the relationship between the wind forcing and the local effective Coriolis frequency in the presence of low-frequency motion. The paper is arranged as follows: Section 2 describes the study area’s oceanographic and meteorological characteristics. Data collection and analysis techniques are outlined in section 3, with the results and discussion presented in section 4. The conclusions and implications of the results are presented in section 5.

2. Study site The study site is the Rottnest continental shelf, southwest Australia (Fig. 1a). This microtidal region has a mean spring tidal range of ;0.6 m and diurnal tides

NOVEMBER 2016

 ET AL. MIHANOVIC

3489

FIG. 1. (a) Locations of the Rottnest continental shelf, Western Australia measurement sites. The major current systems Leeuwin and Capes Currents are shown schematically as southward and northward flowing, respectively. The black squares denote the mooring sites off Two Rocks and in the Perth Canyon (Fig. 1b). The black triangle denotes the Rottnest Island meteorological station. HF radar shore stations were located at Guilderton and Fremantle. Bathymetry contours are in meters. (b) HF radar (WERA) operational coverage between 10 Mar 2010 and 30 Apr 2012. The color scale indicates the percentage of data recovery. The small black dots denote the HFR grid, and the white squares indicate the mooring sites. The larger black dots in the grid denote the two cross-shore transects where the effective Coriolis frequency (period) was estimated.

(Pattiaratchi and Eliot 2009). Continental shelf processes, which are mainly wind driven, include mixing, circulation, and particulate resuspension dominated by LSB activity, particularly during the summer months (Pattiaratchi et al. 1997; Zaker et al. 2007; Verspecht and Pattiaratchi 2010; Gallop et al. 2012). Local sea breezes, superimposed upon synoptic southerly winds, with speeds often .15 m s21, occur almost daily during the austral spring, summer, and early autumn (October–April). The unusually strong alongshore sea breeze results from the interaction of the local sea breeze generated due to the land–sea temperature gradient and the synoptic pressure system (Pattiaratchi et al. 1997). The Leeuwin Current system, which includes the Leeuwin Current (LC), the Leeuwin Undercurrent, and the Capes Current on the shelf (Pattiaratchi and Woo 2009), controls the circulation along the Australian west coast (Fig. 1a). The LC is a warm, low-salinity eastern boundary current, which flows poleward along the continental shelf edge. This is an anomalous direction for an eastern boundary current; although the equatorward winds provide conditions conducive for coastal upwelling, downwelling conditions are a feature of the region (Smith et al. 1991). The LC transports tropical waters southward from northern Australia and flows strongest (5–6 Sv) during the austral autumn/winter (April– September) when the opposing southerlies and sea breezes are weakest (Smith et al. 1991; Pattiaratchi and Woo 2009). The Leeuwin Current system also responds to the El Niño and La Niña cycles with the current being

stronger (weaker) during La Niña (El Niño) events (Pattiaratchi and Buchan 1991). Of particular interest to this study, the region experienced a marine heat wave in February and March 2011, which was associated with the warming related to the La Niña event defined as Ningaloo Niña (Feng et al. 2013). This event increased the Leeuwin Current’s volume transport in February—an unusual event at this time of the year that resulted in unprecedented warm sea surface temperature anomalies (;58C higher than normal) off Australia’s west coast, with widespread coral bleaching and fish kills (Feng et al. 2013; Doi et al. 2013; Kataoka et al. 2014). The Capes Current (Fig. 1a) flows northward along the continental shelf, inshore of the Leeuwin Current (Pearce and Pattiaratchi 1999). The current originates with winddriven upwelling of cooler, saltier waters around 348S, which flows northward particularly during the summer (Gersbach et al. 1999). This water flows along the inner continental shelf and broadens as it moves north from 348S. The Capes Current is strongest in spring and summer when sea breezes and southerly winds prevail. From November to March, the northward wind stresses push the Leeuwin Current offshore and reinforce the Capes Current (Pearce and Pattiaratchi 1999). The Perth Canyon drops from the 200-m isobath, west of Rottnest Island, to ;4000 m within ;50 km (Fig. 1a). Rennie et al. (2006), using thermistor strings, observed the impact of sea breeze on the thermal structure in the Perth Canyon and observed a periodic variation at the diurnal frequency, which manifested in the upper 200 m

511.5 31.78678S, 114.94008E WATR500

206 31.72418S, 115.03828E WATR200

151.5 31.68628S, 115.12878E WATR150

108

24 Sep 2010–25 Mar 2011; 10 min; 32.7, 37.8, 42.9, 48.0, 54.2, 74.3, 94.0, 99.0 m; SBE 39 and SBE 39plus 24 Sep 2010–25 Mar 2011; 10 min; 26.3, 31.4, 41.4, 51.3, 69.6, 89.6, 99.4, 124.2, 147.7 m; SBE 39 and SBE 39plus 22 Dec 2010–18 May 2011; 10 min; 30.0, 39.9, 55.0, 74.0, 103.9, 128.7, 153.5, 178.3, 200.0 m; SBE 39, SBE 39plus, XR420, and Nortek

ADCP data dates; temporal resolution; bin depths; instrument type 31.63338S, 115.19578E

Meteorological data were collected at the Rottnest Island meteorological station (32.00698S, 115.50228E), located 10–30 km from the Perth Canyon sites and 50–60 km from the mooring locations along Two Rocks (Fig. 1a). Wind speed and direction were recorded 10 m above the ground every 30 min, with short (almost insignificant) gaps. These data were interpolated and resampled into 1-h intervals to correspond with the radar data and other analyses. The coordinate system was rotated 258 counterclockwise to align with the coast to provide cross-shore and alongshore winds. Wind vectors are presented in subsequent figures in the direction the

WATR100

1) METEOROLOGICAL DATA

Station depth (m)

Data presented here were collected as part of Australia’s Integrated Marine Observing System (IMOS) and include surface currents from HFR, vertical current profiles, and thermal structure from moored instruments and sea surface temperature from satellite remote sensing augmented by meteorological data. The HFR surface current data from February 2011 and February 2012 were used to examine the horizontal extent of the surface highfrequency current patterns. ADCP currents measured at 200 and 500 m (WATR200 and WATR500, respectively) during the 2010/11 austral summer were analyzed in detail, with a focus on the high-frequency variability observed in February and March 2011 (Table 1); this period was compared with a period during the 2011/12 austral summer (February–March 2012; Table 2), which also showed strong diurnal oscillations.

Latitude, Longitude

a. Data

Station

3. Data collection and analysis

Temperature data dates; temporal resolution; mean sensor depths; instrument types

of the water column. They also identified a phase relationship between the wind and ocean temperature during sea-breeze conditions. The sea breeze is highlighted in their study as a source of the strongest diurnal periodicities and isotherm displacements of up to 80 m over a 24-h period that could not be accounted for by the tide. The 1 yr of current measurements from an ADCP mooring within our study area was used by Cresswell (2009) to highlight the current reversal and variability due to wind forcing and transport off the shelf. The ADCP data revealed periods of several days when crossshelf flow dominated. Because of the coastal boundary effect, the surface and bottom layers were often out of phase. Cresswell (2009) used progressive vector diagrams to show the generation of current oscillations as loops and cusps and found that the LSB cycle produced baroclinic inertial waves.

VOLUME 46 22 Dec 2010–18 May 2011; 10 min; 6-m bins from 30 to 192 m; Nortek Continental 190 kHz 26 Oct 2010–8 Jun 2011; 15 min; 16-m bins from 48 to 480 m; Teledyne RDI Workhorse Long Ranger 75 kHz

JOURNAL OF PHYSICAL OCEANOGRAPHY

TABLE 1. Mooring deployment details for austral summer 2010/11.

3490

 ET AL. MIHANOVIC

NOVEMBER 2016

21 Dec 2011–4 Jul 2012; 10 min; 6 m bins from 30.7 to 192.7 m; Nortek Continental 190 kHz

2) MOORED TIME SERIES DATA

12 Jan 2012–8 Jun 2012; 5 min; 26.8, 31.8, 36.8, 41.8, 50.7, 71.3, 92.3, 99.4 m; SBE 39 and SBE 39plus 21 Dec 2011–4 Jul 2012; 5 min; 25.7, 30.7, 40.8, 50.7, 70.4, 90.7, 101.1, 126.0, 146.6 m; SBE 39, SBE 39plus, and XR420 21 Dec 2011–4 Jul 2012; 5 min; 30.6, 40.5, 55.5, 74.8, 104.7, 129.5, 154.4, 179.2, 200.8 m; SBE 39, SBE 39plus, XR420, and Nortek 20 Oct 2011–11 Apr 2012; 10 min; 25.9, 35.9, 50.9, 70.4, 100.5, 125.7, 150.8, 175.9 m; SBE 39, SBE 39plus, and XR420

Temperature data dates, temporal resolution, mean sensor depths, instrument types

ADCP data dates, temporal resolution, bin depths, instrument type

wind is blowing to but referred to in the text using meteorological convention.

The moorings were located ;40–80 km offshore, 200 km south of the critical latitude for diurnal–inertial resonance at 308S (Fig. 1). The positions of the four northern moorings on the continental shelf off Two Rocks (Fig. 1; Tables 1 and 2) were in water depths of 100, 150, 200, and 500 m, respectively (WATR100–WATR500) and a southern mooring in the Perth Canyon was in water depth of 200 m (WACA200). Data from nearshore (;50-m water depth) moorings WATR050 and NRSROT were not used here; however, the mooring locations are shown in all relevant figures for reference. Bottom-mounted, upward-looking acoustic Doppler current profilers (ADCP) were moored on the continental shelf off Two Rocks in 200- and 500-m water depths (WATR200 and WATR500, respectively, during the 2010/11 austral summer and at WATR200 in 200-m water depth during the 2011/12 austral summer; Fig. 1a). Mooring configurations, ADCP instrument type, and sampling configurations are presented in Tables 1 and 2. The coordinate system was rotated 258 counterclockwise so that u and y represent the cross-shore and alongshore components, respectively. The thermistors (Seabird Electronics SBE 39, SBE 39plus, and RBR XR420) were fixed onto the mooring lines at the Perth Canyon site and at the four sites on the continental shelf off Two Rocks (Fig. 1; see also Tables 1 and 2). Changes in depths with time due to mooring inclinations were taken into account using the pressure sensors associated with the thermistors.

204.5 31.98338S, 115.22778E WACA200

206.5 31.66678S, 115.00008E WATR200

152 31.68678S, 115.12838E WATR150

104 31.64388S, 115.20318E WATR100

Latitude, Longitude

Station depth (m)

3) HF RADAR DATA

Station

TABLE 2. Mooring deployment details for austral summer 2011/12.

3491

Australian Coastal Ocean Radar Network (ACORN) operates surface current mapping HF radar systems as part of the Australian IMOS. This facility provides hourly, time-resolved coastal ocean surface current maps using a high-frequency, phased array wave radar— the Wellen Radar (WERA)—which is shore based and provides a reliable dataset of surface currents (Gurgel et al. 1999). The Rottnest shelf WERA HF radar system transmits at a frequency of 8.5125 MHz with a bandwidth of 33 kHz. High-frequency radar systems were first used to measure surface currents more than three decades ago (Barrick et al. 1977; Hammond et al. 1987; Shay et al. 2007) and are now widely used around the world with a high level of accuracy (Paduan and Washburn 2013). Two WERA HFR shore stations were located at Guilderton and Leighton Beach in Fremantle (Fig. 1a) and the radar coverage area overlapped the mooring sites (Fig. 1b). The WERA instrument system had a maximum

3492

JOURNAL OF PHYSICAL OCEANOGRAPHY

offshore range of ;180 km. The temporal resolution for the HFR measurements was 1 h, and the surface current vectors were derived over a regular grid with a horizontal resolution of 4 3 4 km. WERA HFR data were reprocessed using ACORN-developed software, which included an improved analysis algorithm to step from raw spectra to surface current components and provide quality control flags on the data points (Heron and Prytz 2011).

b. Analysis techniques The Rottnest Island wind vector W time series for February 2011 and February and March 2012 were used to calculate the wind stress vector t. The wind stress was calculated using the quadratic law t 5 raCDjWjW, with air density ra 5 1.25 kg m23 and the drag coefficient CD calculated according to Anderson (1993). The unfiltered wind vectors W and current vectors u were used to estimate the rate of work done by the local wind at WATR200, t  u 5 t x ux 1 ty uy (surface energy flux, W m22), and KE per unit volume KE 5 (1/2)r0juj2 (J m23), where the density r0 was set to 1023.25 kg m23, which represented the average density at 30-m depth at WATR200 for the study period. Each vector time series (winds and currents) was subjected to a 48-h, Butterworth filter (half-power point 5 48 h) to separate the time series into high- and low-frequency components. This allowed for the examination of the high-frequency signal close to the diurnal band that is addressed in this paper. Rotary energy spectra were calculated for the Rottnest Island wind data (22 December 2010–22 December 2012) following Gonella (1972) and using the fast Fourier transform (FFT) method and a 150-day, half-overlapping Hamming window. Similarly, rotary spectra were calculated for the ADCP currents (22 December 2010–15 April 2011) with a 19-day overlapping Hamming window. Wavelet analysis was used for the wind data from the Rottnest Island meteorological station and all available current data (30-m depth) at the WATR200 ADCP station between 22 December 2010 and 22 December 2012 (four ADCP deployments). Wavelet analysis examined spectra moving in time, rather than one spectrum for the whole dataset, which revealed the seasonal changes in the dominant frequencies. The analyses were based on wavelet methods given in Torrence and Compo (1998) and Grinsted et al. (2004), with wavelet bias corrections from Liu et al. (2007) and Veleda et al. (2012). The WERA HFR-derived surface current vectors with a large geometrical dilution of precision (GDOP) caused by poor intersecting beam geometry were excluded from the analysis (Chapman et al. 1997; Shay et al. 2007). The resulting HFR data were examined for spikes, with anomalous values removed from the dataset. The data

VOLUME 46

were then spatially interpolated with the condition that the missing grid point within the area had to have a small GDOP and at least three contemporaneous hourly measurements within the 7-km radius. Short, nonoperational periods (maximum of 3 h) were also interpolated in time. The coordinates were then rotated 258 counterclockwise so that u and y represent cross-shore and alongshore surface components, respectively. The percentage of HFR data operational coverage over the Perth Canyon and Two Rocks mooring sites was ;99% (Fig. 1b) from 10 March 2010 to 30 April 2012, based on 14 691 operational temporal points: that is, temporal points with at least one gridpoint measurement (out of an overall 18 784 hourly points available during the analyzed period; the radars were not operational through the entire period). HFR surface currents were presented as a temporal evolution of surface currents at positions corresponding to the three mooring sites (WACA200, WATR200, and WATR500; Fig. 1b) in the 2010/11 and 2011/12 austral summers. The HFR surface current data were also used to estimate the surface KE per unit volume and the rate of work done by the local wind at WATR200 plus the local effective Coriolis frequency (period) along two transects corresponding to mooring sites (see Fig. 1b). Complex demodulation (Emery and Thomson 1998) was used to obtain the clockwise and counterclockwise phases and amplitudes of the Rottnest Island wind data, ADCP current measurements, and the HFR currents data. We used 15–35-h, bandpassed filtered winds and currents (Nam and Send 2013) to extract the phases and amplitudes of the 24-h harmonics. Two diurnal periods of 48 h were overlapped and slid in time with 6-h time steps (Emery and Thomson 1998; Jarosz et al. 2007). Inertia–gravity waves can be generated when the forcing frequency is larger than the lowest frequency allowed. If the vertical shear is negligible the effective Coriolis frequency is the lowest frequency allowed (Federiuk and Allen 1996). However, if vertical shear is large the lowest frequency is (Federiuk and Allen 1996) vM 5 feff 2

f (j›u/›zj)2 . 2N 2

(1)

Therefore, vertical shear and stratification can decrease the lowest allowed frequency (increase the longest allowed period). The vertical shear and buoyancy frequency were estimated using HFR and ADCP data and the thermistors.

4. Results and discussion a. Wind forcing The diurnal LSB system off the Western Australian coast is among the globally strongest (frequently . 10 m s21

NOVEMBER 2016

 ET AL. MIHANOVIC

3493

FIG. 2. Normalized wavelet spectra for wind from Rottnest Island meteorological station and currents at WATR200 (;30 m) for 22 Dec 2010 to 22 Dec 2012: (a) u (cross shore) wind component, (b) u (cross shore) current component, (c) y (alongshore) wind component, (d) y (alongshore) current component. The thick black line indicates the 95% confidence level, with red noise as the background spectrum. The blue line denotes the diurnal period. Four individual ADCP deployments are presented in (b) and (d). Lighter shadings on either end denote regions where edge effects become important.

with maxima reaching 20 m s21) and most consistent, resulting from a combination of forcing through the synoptic system and the diurnal heating–cooling cycle (Pattiaratchi et al. 1997; Masselink and Pattiaratchi 2001). Wavelet analysis of 2 yr of measured wind data (22 December 2010–22 December 2012) from the Rottnest

Island meteorological station revealed energy, with a period of 1 day, was present in austral spring, summer, and early autumn but absent in winter (Figs. 2a,c). Energy in the weather band (periods 5–20 days) was present during winter. The energy at the diurnal frequency was higher in the u (cross shore) component than in the y component

3494

JOURNAL OF PHYSICAL OCEANOGRAPHY

VOLUME 46

FIG. 3. Rotary energy spectra calculated for the period 22 Dec 2010 to 22 Dec 2012 (wind) and the period 22 Dec 2010 to 15 Apr 2011 (currents): (a),(b) wind data from Rottnest Island meteorological station, (c),(d) currents at 30-m depth at WATR200, and (e),(f) currents at 132-m depth at WATR200. Clockwise components are denoted in (a),(c), and (e) and counterclockwise components are denoted in (b),(d), and (f).

due to cross-shore nature of the LSB in the region. Similarly the ADCP currents at 30-m water depth at WATR200 (Fig. 1b) indicated higher energy in the diurnal band during the summer months with energy at the diurnal frequency being higher in the u (cross shore) component than in the y component (Figs. 2b,d). The LSB cycle indicated counterclockwise rotation (Figs. 3, 4) with the offshore extent exceeding 150 km as revealed by the surface current response that extended to the limit of the HFR coverage (Figs. 4, 5). This compares with the offshore extent of sea breezes to 250 km reported along the east and west coasts of the

United States by Edwards (2008) and Lerczak et al. (2001), respectively.

b. HF radar surface currents To examine the mean response of the coastal surface ocean to the LSB cycle, the 48-h, high-pass filtered HFR surface currents for the entire month of February 2011 were used to construct mean values over a diurnal cycle: that is, the mean values at a particular time stamp. For example, 0130 local time (LT) is the mean value recorded at that time stamp for the whole month of February 2011 (Fig. 4a). Individual mean hourly

NOVEMBER 2016

 ET AL. MIHANOVIC

3495

FIG. 4. The 3-hourly snapshots of the 48-h, high-pass filtered, high-frequency HFR surface currents representing the mean daily cycle in February 2011. A vector corresponding to wind at Rottnest Island for each snapshot is shown in the bottom-left corner of each panel. Time flags represent local time (LT, 18 h UTC).

snapshots contained points that had at least 10 out of 28 measurements representing that particular hour (Fig. 4). The mean daily wind cycle starts at 0130 LT with a weak land breeze, that is, northeasterly winds (offshore) that reach a maximum speed of 3.4 m s21 at 0730 LT (Fig. 4c). The land breeze decreases in strength by 1030 LT (Fig. 4d) and then switches direction to the sea breeze by 1330 LT (Fig. 4e). The sea breeze (from the southwest) then increases in speed with a maximum speed 5.4 m s21 at 1630 LT (Fig. 4f) and then decreases by the end of the cycle (Fig. 4h). Note that the wind speeds were means of highpassed winds over 1 month with the instantaneous wind speeds being much higher for a particular day. The mean daily cycle of high-frequency surface currents reveals an almost instantaneous response of the surface currents to changes in the wind direction and a strong counterclockwise rotation that corresponds to the diurnal rotation of the wind vector. There is a clear offshore intensification of the surface currents, with speeds . 0.3 m s21 in the western part of the radar domain (Fig. 4). The diurnal ellipse characteristics of the Rottnest Island winds in February 2011 indicated a rectilinear ellipse with counterclockwise rotation with the major axis tilted 358 counterclockwise from the east (Fig. 5), reflecting the direction of the LSB system. The diurnal surface current ellipses for the same month indicated circular motion in the deeper water above the Perth Canyon (Fig. 1a), while on the continental shelf region in water depths , 50 m, the circulation was less circular (Fig. 5). The surface currents rotated counterclockwise

across the whole region. The current amplitudes were up to a factor of 5 larger in the deeper offshore regions compared to those on the continental shelf (Fig. 5). Numerical model simulations by Davies and Xing (2005) indicated that the near-surface response to wind

FIG. 5. Diurnal wind ellipse (black) and surface current ellipses (gray) calculated from the February 2011 data. Wind and current vectors at local time midnight are shown within each ellipse. For clarity, current ellipses denote every third HFR point that had at least two-thirds of temporal coverage in February 2011. All ellipses rotated counterclockwise.

3496

JOURNAL OF PHYSICAL OCEANOGRAPHY

VOLUME 46

FIG. 6. High-frequency currents (m s21) between 1 Feb and 28 Feb 2011: WATR200 (a) crossshore u component and (b) alongshore y component; WATR500 (c) cross-shore u component and (d) alongshore y component. Note that the panels have different depth scales.

forcing is independent of the magnitude of forcing but rather increases with distance offshore. Similarly, Kaplan et al. (2005) reported, through the analysis of HFRderived surface currents, that in Bodega Bay (northern California) the near-inertial variance was at maximum offshore, with significant decrease in near-inertial energy over the inner shelf region. The more polarized, decreased surface current magnitudes in the inner shelf region are a feature of the interaction of inertial motions with the coastal boundary (Kaplan et al. 2005).

c. Two Rocks ADCP currents The rotary energy spectra of the currents at WATR200 indicated a narrow diurnal peak in the kinetic energy (KE) density of the positive, counterclockwise, rotary component at 30 and 132 m below the surface (Figs. 3d,f) and a negligible peak in the negative, clockwise component (Figs. 3c,e). These results indicated that the diurnal currents were enhanced at the surface and could be due to direct wind forcing. Another mechanism that may enhance the surface currents is subharmonic parametric instability (SHPI) in the tides, where the instability grows if

the energy level of the internal waves is sufficiently strong (Simmons 2008). Here, when the internal tide crosses the turning latitude (for frequencies lower than the M2, the turning latitudes are located at 628.88) of its subharmonic, instabilities will develop where energy is transmitted to a subharmonic with high-mode vertical structure and to super harmonics (Simmons 2008). Off the coast of Namibia (;308S) both wind forcing (Simpson et al. 2002) and the SHPI mechanism (Simmons 2008) indicated a strong maximum. The results of Simmons (2008; Fig. 3) indicate an absence of SHPI energy along the east coast of Australia, most likely due to the dominance of diurnal tides in this region (Pattiaratchi and Eliot 2009). Thus, the diurnal currents in the study region results were enhanced at the surface due to direct wind forcing. Vertical profiles of the high-frequency currents from the shelf mooring (WATR200; Figs. 6a,b) and deep mooring (WATR500; Figs. 6c,d) at Two Rocks in February 2011 revealed a diagonal banding in the upper 200–300 m. At each depth level, a daily change in direction (alternating between positive and negative) was observed in both cross-shore u and alongshore y components

NOVEMBER 2016

 ET AL. MIHANOVIC

(Fig. 6). The banding strength in the time series varied; however, the high-frequency current amplitudes were .0.3 m s21 at both sites (at 200- and 500-m water depths). A counterclockwise rotation in the upper 200–300 m, where the u component leads the y component by p/2 (or 908), can be seen in Fig. 6. A counterclockwise rotation in the Southern Hemisphere (with the u component leading y; Simpson et al. 2002) indicates upward phase propagation (Lerczak et al. 2001; Nam and Send 2013). This banding feature, which is caused by obliquely propagating waves in a continuously stratified fluid, defined as inertia–gravity waves, is often observed around critical latitudes (Kundu 1990; Sutherland 2010) and was observed to 300-m water depth at WATR500 (Figs. 6c,d)—a much greater depth than the Ekman depth of ;70 m at this location. The low-pass filtered temperature measurements were used to define the mixed layer depth (MLD) as the depth where the difference between the uppermost available layer and lower layers was .0.58C (Monterey and Levitus 1997). The MLD at WATR200 in February 2011 fluctuated between 30 and 100 m and was deepest around 15 February 2011 due to the marine heat wave associated with Ningaloo Niña (Feng et al. 2013). This warming was strongest between 17 and 25 February 2011, and the water column was less stratified than it was at the beginning of the month. Deep diurnal banding (between 200 and 300 m) was observed at WATR500 during this time (Figs. 6c,d). The rotary energy spectra for the upper (30 m) and lower (132 m) layers at WATR200 both revealed a peak at the diurnal frequency for the counterclockwise component, but not for the clockwise component (Figs. 3c–f), indicating the presence of the diurnal signal through the water column. However, the cross-shore and alongshore components showed a pronounced phase shift of 1808 below the mixed layer when compared to the surface. Downward energy flux and upward phase propagation were visible, and the vertical phase speed around 10 February 2011 was estimated to be 130–140 m day21, which resembled the phase speed observed in De Soto Canyon, Gulf of Mexico (138 m day21), in summer (Jarosz et al. 2007), but was about twice as fast as the speed observed close to Mission Beach, California (68 m day21; Lerczak et al. 2001). The LSB system is modulated by the west to east propagation of the high pressure system (Verspecht and Pattiaratchi 2010), and during the period 1 February to 19 March 2011, the LSB system reflected a periodicity of around 5–6 days (Fig. 7a). During strong diurnal wind forcing periods, corresponding to larger amplitudes in Fig. 7a, the counterclockwise wind phase was relatively constant around 2908 and only fluctuated about 308 from this value when the diurnal wind energy was low

3497

(Fig. 7c). The current phases from the HFR and ADCP measurements varied slightly. They gradually decreased from 1 to 26 February 2011 (Fig. 7c) and then remained around a relatively constant value. The counterclockwise phase at two depths (30 and 132 m) from the ADCP data at WATR200 revealed a distinct phase difference between the water column’s upper and lower layers (Figs. 7c,d), which was close to 1808 for most of the study period. These observations are consistent with the mechanism of resonant response to the LSB with a coastal boundary, proposed for regions close to the critical latitude where it was also found that the directions and phases of surface currents match those of wind stress at the diurnal frequency at the surface with an 1808 phase difference between upper and lower layers (e.g., Simpson et al. 2002; Edwards 2008; Zhang et al. 2009; Nam and Send 2013).

d. Local effective Coriolis frequency The wind-forcing frequency must be near the local effective Coriolis frequency to generate near-inertial oscillations (Kunze 1985; Nam and Send 2013). Data collected from austral summers (February 2011 and February–March 2012) were used to examine the relation between the wind forcing and the local effective Coriolis frequency. In February and March 2011, the region experienced a marine heat wave associated with the Ningaloo Niña event resulting in ocean temperatures ;58C higher than normal (Feng et al. 2013). Time series maps combining sea surface temperatures, altimetry data, and the daily averaged surface currents obtained from HFR (Fig. 8) were used to examine the synoptic-scale dynamics in February 2011. In early February 2011, a system of two counterrotating eddies dominated the synoptic-scale dynamics (Fig. 8a). A clockwise eddy, located to the west of the study region, impacted the deeper Two Rocks moorings (WATR150–WATR500) with southward-flowing currents under the influence of the Leeuwin Current (eastern rim of the clockwise eddy). The continental shelf region zone was under the influence of the northward-flowing Capes Current (WATR050 and WATR100). A counterclockwise eddy was present to the northwest, and the warmer Leeuwin Current water was advected along the western and southern perimeters of this eddy and then entrained into the southern clockwise eddy and advected south (Fig. 8a). The clockwise eddy was dominant in the Perth Canyon until about 12 February 2011. After this date, the strengthening Leeuwin Current (Figs. 8b–d) transported warm water southward across the study region (Figs. 9g,h,i), which reduced the buoyancy frequency through the water column at WATR200 (Fig. 9f) and dissipated the clockwise and counterclockwise eddies.

3498

JOURNAL OF PHYSICAL OCEANOGRAPHY

VOLUME 46

FIG. 7. (a) Counterclockwise (gray) and clockwise (black) amplitudes of diurnal winds from the Rottnest Island meteorological station for 1 Feb to 19 Mar 2011; (b) counterclockwise amplitudes of HFR-measured diurnal surface currents at WATR200 (thick black line) and ADCP-measured diurnal currents at WATR200 (dashed line: 30-m depth; thin black line: 132-m depth) for 1 Feb to 19 Mar 2011; and (c) counterclockwise diurnal phases for the Rottnest wind data (gray line), HFR-measured surface currents at WATR200 (thick black line), and ADCP-measured currents at WATR200 (dashed line: 30-m depth; thin black line: 132-m depth). As the HFR and ADCP currents were mostly counterclockwise, only counterclockwise amplitudes and phases [(b) and (c), respectively] are shown; (d) the difference in phases between counterclockwise diurnal winds at Rottnest and counterclockwise diurnal HFR-measured currents (thick black line), and ADCP-measured currents at WATR200 (dashed line: 30-m depth; thin black line: 132-m depth).

The influence of the subinertial shear on the effective Coriolis period was substantial. At the start of February 2011, the clockwise (cyclonic) eddy increased (in a negative sense) the relative vorticity at WATR200 and decreased the local Coriolis period to less than 17 h on 8 February 2011 (Fig. 9c). As the clockwise eddy moved southward, the vorticity became positive, which increased the local Coriolis period to around 24 h on 10 February 2011. As the diurnal wind was also strong at the time (Fig. 9a), energetic

counterclockwise diurnal oscillations below the MLD at WATR200 were also observed (Fig. 9d). The change in the local vorticity shifted the internal wave band’s lower bound, which allowed the intermittent, resonant generation of near-inertial waves and their deep penetration below the MLD. The effective Coriolis frequency’s temporal variation likely affected the drift of phase variation of counterclockwise current at WATR200 (Fig. 7c); however, the moorings

NOVEMBER 2016

 ET AL. MIHANOVIC

3499

FIG. 8. Sea surface temperature (color), altimetry (m; gray contours), daily averaged HFR surface currents (black arrows), and daily averaged winds from the Rottnest Island meteorological station (red arrows) for 4 days during the February 2011 Leeuwin Current marine heat wave. For clarity, only every second surface current vector is shown. The temperature scale is the same for all panels.

at WATR050 and WATR150 were mostly in superinertial regime during February 2011 (and also April 2011; see Fig. 10) because of the Leeuwin Current’s more onshore position and the weaker and narrower Capes Current (Fig. 8b). The resulting counterclockwise shear decreased the local effective Coriolis frequency in that part of the Two Rocks transect. The land–sea breeze system dominated the wind regime throughout February 2011 (Fig. 9a). The surface currents measured by HFR indicated large, counterclockwise oscillations in the surface currents at the three ADCP mooring sites (Fig. 9b). From 9 February 2011 onward, banding was seen extending through the whole water column in the 200-m ADCP current profiles (Figs. 6a,b) and down to 250–300 m in the 500-m ADCP current profiles (Figs. 6c,d). Counterclockwise motions were dominant at WATR200,

and the upper mixed layer and lower layer were clearly separated (Fig. 9d). The 100- and 150-m moorings (Figs. 9g,h, respectively) measured diurnal isotherm oscillations, with brief vertical isotherm fluctuations reaching almost 50 m at the 100-m mooring (WATR100) on 11 February 2011. As the thermocline deepened (because of the warmer Leeuwin Current waters flowing south), the 200-m Two Rocks mooring measured isotherm oscillations of about 15–20 m around 22–23 February 2011 (Fig. 9i). The above discussion considered the effect of horizontal shear on the effective Coriolis frequency on diurnal–inertial motion. Previous studies (Kunze 1985; Federiuk and Allen 1996; Nam and Send 2013) have demonstrated that vertical shear may also play a significant role in modulating the inertial–gravity wave energy. The low-frequency vertical shear and buoyancy

3500

JOURNAL OF PHYSICAL OCEANOGRAPHY

FIG. 9. February 2011. (a) Counterclockwise (red) and clockwise (blue) amplitudes of the diurnal winds for Rottnest Island; (b) counterclockwise (red) and clockwise (blue) amplitudes of the diurnal HFR surface currents at WACA200, WATR500 and WATR200; and (c) effective Coriolis period at WATR200 (blue line) and the corresponding planetary Coriolis period (red line). The black line denotes the diurnal period representing the wind forcing. (d) Counterclockwise amplitudes (m s21) of the diurnal ADCP currents at WATR200. The white line denotes the mixed layer depth. (e) Vertical shear squared (s22) for low-passed ADCP currents at WATR200. The white line denotes the mixed layer depth. HFR measurements were used to calculate the shear between the surface and 30-m depth when available. (f) Brunt–Väisälä frequency squared (s22) at WATR200. The black line denotes the mixed layer depth. Temperature profile (8C) at (g) WATR100, (h) WATR150, and (i) WATR200. In (g),(h), and (i), the black line denotes the 23.38C isotherm.

VOLUME 46

NOVEMBER 2016

 ET AL. MIHANOVIC

3501

FIG. 10. The effective Coriolis period (h) estimated along the Two Rocks transect (upper panel) and the Perth Canyon transect (lower panel) between January and April 2011. The black contour lines denotes the 24-h period. See Fig. 1b for the transect positions. The mooring locations are also indicated.

frequency estimated from the HFR and ADCP data and thermistor measurements indicated that the vertical shear squared was much smaller than buoyancy frequency squared for majority of the time (Figs. 9e,f). Hence, the lowest frequency allowed was the effective Coriolis frequency [see Eq. (1)]. However, on 9–10 February 2011 and just below relatively shallow mixed layer on 23–24 February 2011, the ratio between vertical shear squared and buoyancy frequency squared was relatively high (Figs. 9e,f): in the first case (9–10 February 2011), it enabled resonance in the middle layer (decreasing the lowest frequency allowed toward diurnal frequency), while on 23–24 February 2011 (when the effective Coriolis frequency was close to the diurnal; Fig. 9c), it also decreased the lowest frequency allowed, thus shifting it away from the diurnal resonance, which therefore was mostly pronounced at the surface and in the mixed layer. The surface energy flux during February 2011 was mostly negative due to opposing directions of southwardflowing Leeuwin Current and southerly winds (Fig. 11b). The influence of the Leeuwin Current may be seen as the negative component in the low-passed energy flux (Fig. 11c). Diurnal oscillations of surface KE per unit volume were highest between 12 and 15 February 2011 (oscillating between 40 and 250 J m 23 ) and around 23–25 February 2011, when they fluctuated between 0 and 60 J m23 (Fig. 11e). The energetic oscillations occurred when the oscillatory component of the energy flux (Fig. 11d) and the effective Coriolis period before/ during those two periods was close to 24 h (Fig. 11a).

The conditions were markedly different in February– March 2012 (Fig. 12). In early February 2012, a single counterclockwise eddy, centered on the Perth Canyon, dominated the synoptic-scale dynamics (Fig. 12a). The continental shelf region zone was under the influence of the northward-flowing Capes Current. With time, the counterclockwise eddy dissipated, and the southward Leeuwin Current was established in the offshore region (Fig. 12d). Several episodes of energetic diurnal motions extending through the whole water column occurred at WATR200 in February and March 2012 (Figs. 13d,e). The vorticity along the Two Rocks transect was positive at deeper stations more often (resulting in a lower effective Coriolis frequency; Fig. 14) than it was in February 2011 because of the Leeuwin Current’s more offshore position (cf. Fig. 8 and Fig. 12). Such conditions were more favorable for effective pumping of the diurnal wind energy into the ocean. The counterclockwise diurnal amplitudes of the HFR surface currents at WATR200 indicated three pronounced peaks: around 8 February, 3–4 March, and 24 March 2012 (Fig. 13b). The effective Coriolis period at station WATR200 was .24 h at the beginning and the end of February 2012 and through most of March 2012 and was close to 24 h during the three pronounced peaks mentioned above (Fig. 13c). For these three episodes the ratio between vertical shear and buoyancy frequency was insignificant (not shown), and the effective Coriolis frequency was the lowest frequency allowed. Therefore, when the counterclockwise diurnal wind amplitude peaks corresponded to the 24-h effective Coriolis period, resonance occurred. Strong diagonal diurnal banding was

3502

JOURNAL OF PHYSICAL OCEANOGRAPHY

VOLUME 46

FIG. 11. February 2011. (a) the effective Coriolis period at WATR200 (blue) and the corresponding planetary Coriolis period (red). The black line indicates the diurnal period representing the wind forcing. (b) The rate of work done by the local wind (surface energy flux, W m–2) at WATR200 calculated using the Rottnest Island wind data and the HFR surface currents data; (c) as in (b) but low pass filtered (‘‘mean’’); (d) as in (b) but high pass filtered. (e) Kinetic energy per unit volume (J m–3) at WATR200 estimated from the HFR surface currents data, and (f) vertical profile of kinetic energy per unit volume (J m–3) at WATR200 calculated from the ADCP data. The black line in (f) denotes the mixed layer depth. The periods when the effective Coriolis period was close to the diurnal period are shaded (see Figure 13c).

observed for the cross-shore and alongshore current components at WATR200 (Figs. 13d,e), with nearbottom, cross-shore currents flowing in the opposite direction to the cross-shore wind stress. LSB circulation was well developed in February and March 2012 but was especially energetic when maximum diurnal currents occurred (Fig. 13a). Between 2 and 6 March 2012, the diurnal wind ellipse was less polarized than usual, with the counterclockwise wind following a more circular path (not shown). The LSB also induced near-inertial oscillations more effectively when the MLD was smaller (e.g., at the start of March 2012 when the MLD was around 50–55 m). Diurnal temperature oscillations were observed in February and March 2012 and were vigorous at WATR150 around 4 March 2012 and from 18 to 20 March 2012 (the isotherm range reached 60 m; Fig. 13g). Sporadic diurnal temperature oscillations

were present at all Two Rocks sites. A prolonged period of diurnal temperature variability was also observed at WACA200 in March 2012, with the vertical isotherm fluctuations reaching 20 m (Fig. 13i). The surface energy flux during this period was mostly positive, with diurnal bursts surpassing 0.04 W m22 during the following periods: 4–8 February 2012, 13–15 February 2012, 29 February–5 March 2012, and 22–24 March 2012 (Fig. 15). Maximum value of diurnal bursts reached 0.1 W m22 (Fig. 15b), more than 5 times larger than the flux observed close to Huntington Beach, California (;33.638N), by Nam and Send (2013), who also suggested that wind stress enhanced the currents. The diurnal oscillations of surface KE per unit volume were also higher on these dates, except between 13 and 15 February 2012 (Fig. 15e). The analysis of the effective Coriolis period (Fig. 15a) indicated large diurnal surface KE peaks occurred when the effective period was close

NOVEMBER 2016

 ET AL. MIHANOVIC

3503

FIG. 12. Sea surface temperature (color), altimetry (m; gray contours), daily averaged HFR surface currents (black arrows), and daily averaged winds from the Rottnest Island meteorological station (red arrows) for 4 days in February and March 2012. For clarity, only every second surface vector is shown. The temperature scale is the same for all panels.

to 24 h. The peaks were not observed between 13 and 15 February 2012 when the effective Coriolis period was much shorter (,20 h; Fig. 15a), representing subinertial conditions. These results are similar to those observed in February 2011 (Fig. 11) where the diurnal oscillations occurred when the oscillatory component of the energy flux (Fig. 15d) and the effective Coriolis period before/ during those two periods was close to 24 h (Fig. 15a). The surface KE per unit volume showed strong diurnal fluctuations between 1 and 4 March 2012. On 4 March 2012, the fluctuations surpassed 400 J m23, which was more than twice as large as the near-surface KE that Nam and Send (2013) observed. Similar surface KE values were observed on 4 February 2012 and somewhat smaller between 18 and 20 March 2012; however, the diurnal variability was less and the rate of work done by the local wind was much larger (in

absolute sense) than it was on 4 March 2012. The surface energy flux was high and positive (up to 0.3 W m22) on 4 February 2012 due to the constructive influence of strong southerly winds and a large, counterclockwise eddy to the west of WATR200 (Fig. 12a); however, the rate of wind work at the surface was high but negative on 19 March 2012 because of the opposing effect of the strong southerly winds and the southward-flowing Leeuwin Current (Fig. 12c). The KE’s temporal evolution at WATR200 showed the KE increased in depth during periods with diurnal banding, especially between 2 and 4 March 2012 (Fig. 15c) when the effective Coriolis period was almost exactly 24 h (Fig. 15a). The low-frequency surface current systems observed during the February 2011 were different from those observed during February–March 2012. In 2011, the Ningaloo Niña event strengthened the Leeuwin Current

3504

JOURNAL OF PHYSICAL OCEANOGRAPHY

FIG. 13. February–March 2012. (a) Counterclockwise (red) and clockwise (blue) amplitudes of the diurnal winds for Rottnest Island; (b) counterclockwise (red) and clockwise (blue) amplitudes of the diurnal HFR surface currents at WACA200, WATR500, and WATR200; and (c) effective Coriolis period at WATR200 (blue) and the corresponding planetary Coriolis period (red). The black line denotes the diurnal period representing the wind forcing. (d) Highfrequency cross-shore u component at WATR200; (e) high-frequency alongshore y component at WATR200. The black line in (d) and (e) denotes the mixed layer depth. Temperature profile (8C) at (f) WATR100, (g) WATR150, (h) WATR200, and (i) WACA200. In (f),(g),(h), and (i), the black line denotes the 228C isotherm.

VOLUME 46

NOVEMBER 2016

 ET AL. MIHANOVIC

3505

FIG. 14. The effective Coriolis period (h) estimated along the Two Rocks transect (upper panel) and the Perth Canyon transect (lower panel) between January and April 2012. The black contour lines denotes the 24-h period. See Fig. 1b for the transect positions. The mooring locations are also indicated.

(Feng et al. 2013). The mean low-frequency southerly wind component was weaker in February 2011 (4.1 m s21) than it was in February 2012 (5.4 m s21), which caused a weaker and narrower Capes Current, fewer counterclockwise mesoscale eddies, and a more onshore position of the Leeuwin Current (Fig. 8). Such surface dynamics resulted in a negative relative vorticity and a higher effective Coriolis frequency at depths greater than 150 m along the Two Rocks transect and at 200 m along the Perth Canyon transect in February 2011 (Fig. 10); however, at depths between 50 and 150 m, the system was mainly in the superinertial regime (especially along the Two Rocks transect), where resonance or freely propagating waves may be expected to occur (Fig. 10). The region with the superinertial regime was shifted farther offshore in February–March 2012 (Fig. 14), and when the effective Coriolis period was close to 24 h, strong, deep diurnal banding was observed at WATR200 (Figs. 13d,e), especially between 2 and 4 March 2012 when the effective period was almost exactly 24 h. The downward energy and upward phase propagation (Figs. 13d,e) confirmed the surface was the energy source for diurnal motions, which was consistent with the resonant response to LSB forcing at the time. The resonant excitation was strongest at the start of March 2012 when the ratio between the surface KE and the work done by the local wind on the surface water was the largest (Fig. 15); this finding supports the notion that low-frequency shears created the conditions that allowed a surface current resonant response to the LSB.

5. Conclusions and implications This paper presented analysis of a set of multiplatform data with wide spatial and temporal resolution to examine the diurnal LSB influence on the surface and subsurface currents and temperature structure along the Rottnest continental shelf in southwestern Australia. The main conclusions were as follows: 1) The study region is subjected to strong LSB forcing along an axis oriented ;358–458 counterclockwise from the east, with counterclockwise rotation extending .150 km offshore. 2) The surface diurnal current response in deeper water was circular and counterclockwise, with current amplitudes . 0.3 m s21. The current amplitudes were stronger in the offshore region when compared to the inner shelf that had weaker currents and more eccentric amplitudes due to the interaction of inertial motions with the coastal boundary. 3) Energy in the diurnal band was observed up to 300-m depth and penetrated deep below the mixed layer. Vertical isotherm fluctuations up to 60 m accompanied these motions with an upward phase propagation speed of ;140 m day21. 4) The upward phase propagation and deep penetration of diurnal currents below the mixed layer and the ;1808 phase difference between the upper and lower water column suggested the local LSB system caused the resonant diurnal motions.

3506

JOURNAL OF PHYSICAL OCEANOGRAPHY

VOLUME 46

FIG. 15. As in Fig. 11 but for February–March 2012.

5) The LSB forcing in the study region area is slightly subinertial, and the rate of wind energy transfer at the surface became more efficient when the effective Coriolis period was ;24 h and the mixed layer was shallower (e.g., at the start of March 2012). The effective Coriolis frequency was also highly dependent on the wider-scale ocean dynamics (such as the marine heat wave in 2011) and the low-frequency southerly winds. In the absence of strong tidal mixing, the intermittent, resonant LSB forcing of near-inertial waves and isotherm oscillations can result in vertical mixing across the pycnocline. The Rottnest continental shelf is oligotrophic due to the dominance of the Leeuwin Current that inhibits large-scale upwelling, with low productivity despite episodic upwelling events during the summer months (Hanson et al. 2005). There is considerable conjecture with regard to the supply of nutrients to the upper ocean in this region (Feng et al. 2009). Vertical mixing has been proposed as an important mechanism for delivering nutrients, from below the Leeuwin Current, to the surface mixed layer (Thompson et al. 2011). Lucas et al. (2014) also studied the diurnal–inertial resonance in the Benguela upwelling system off Namibia and demonstrated that wind forcing enhanced phytoplankton

productivity through a diurnal resonance mechanism embedded in the upwelling’s low-frequency vertical flux forced by Ekman dynamics. Thus, enhanced surface currents’ shear and mixing through the water column through diurnal–inertial resonance could be an important mechanism for nutrient supply to the upper ocean during the summer months in the study region, where low-frequency upwelling episodes driven by strong southerly winds occur often in summer concurrently with land–sea breezes. Acknowledgments. Hrvoje Mihanovic´ acknowledges support received through the Group of Eight (Go8) European Fellowship during his residence at The University of Western Australia (Go8 European Fellowship 2011). We are indebted to Arnstein Prytz for his work on the ACORN analysis algorithm and quality control for the WERA data. The authors thank John Simpson for his early comments on the work and review. The oceanographic data (moorings, HF radar, and satellite SST and altimetry) used in this paper were collected as part of the Integrated Marine Observing System (IMOS), funded by the Australian Government through the National Collaborative Research Infrastructure Strategy and the Super Science Initiative. The meteorological data were obtained from the Bureau of Meteorology.

 ET AL. MIHANOVIC

NOVEMBER 2016 REFERENCES

Anderson, R. J., 1993: A study of wind stress and heat flux over the open ocean by the inertial-dissipation method. J. Phys. Oceanogr., 23, 2153–2161, doi:10.1175/1520-0485(1993)023,2153: ASOWSA.2.0.CO;2. Barrick, D. E., M. W. Evans, and B. L. Weber, 1977: Ocean surface currents mapped by radar. Science, 198, 138–144, doi:10.1126/ science.198.4313.138. Burchard, H., and T. P. Rippeth, 2009: Generation of bulk shear spikes in shallow stratified tidal seas. J. Phys. Oceanogr., 39, 969–985, doi:10.1175/2008JPO4074.1. Chapman, R. D., L. K. Shay, H. C. Graber, J. B. Edson, A. Karachintsev, C. L. Trump, and D. B. Ross, 1997: On the accuracy of HF radar surface current measurements: Intercomparisons with ship-based sensors. J. Geophys. Res., 102, 18 737–18 748, doi:10.1029/97JC00049. Chen, C., R. O. Reid, and W. D. Nowlin Jr., 1996: Near-inertial oscillations over the Texas-Louisiana shelf. J. Geophys. Res., 101, 3509–3524, doi:10.1029/95JC03395. Clancy, R. M., J. D. Thompson, H. E. Hurlburt, and J. D. Lee, 1979: A model of mesoscale air-sea interaction in a sea breeze-coastal upwelling regime. Mon. Wea. Rev., 107, 1476–1505, doi:10.1175/ 1520-0493(1979)107,1476:AMOMAS.2.0.CO;2. Craig, P. D., 1989: A model of diurnally forced vertical current structure near 308 latitude. Cont. Shelf Res., 9, 965–980, doi:10.1016/0278-4343(89)90002-2. Cresswell, G. R., 2009: Current measurements at a year-long continental shelf mooring near Perth, Western Australia. J. Roy. Soc. West. Aust., 92, 165–178. Cushman-Roisin, B., 1994: Introduction to Geophysical Fluid Mechanics. Prentice Hall, 320 pp. Davies, A. M., and J. Xing, 2005: The effect of a bottom shelf front upon the generation and propagation of near-inertial internal waves in the coastal ocean. J. Phys. Oceanogr., 35, 976–990, doi:10.1175/JPO2732.1. DiMarco, S. F., M. K. Howard, and R. O. Reid, 2000: Seasonal variation of wind-driven diurnal current cycling on the TexasLouisiana continental shelf. Geophys. Res. Lett., 27, 1017– 1020, doi:10.1029/1999GL010491. Doi, T., S. K. Behera, and T. Yamagata, 2013: Predictability of the Ningaloo Niño/Niña. Sci. Rep., 3, 2892, doi:10.1038/srep02892. Edwards, C. R., 2008: Coastal ocean response to near-resonant sea breeze/land breeze near the critical latitude in the Georgia Bight. Ph.D. thesis, University of North Carolina-Chapel Hill, 124 pp. Ekman, V. W., 1905: On the influence of the Earth’s rotation on ocean-currents. Arkiv Mat. Astron. Fys., 2, 1–52. Emery, W. J., and R. E. Thomson, 1998: Data Analysis Methods in Physical Oceanography. Pergamon, 634 pp. Federiuk, J., and J. S. Allen, 1996: Model studies of near-inertial waves in flow over the Oregon continental shelf. J. Phys. Oceanogr., 26, 2053–2075, doi:10.1175/1520-0485(1996)026,2053: MSONIW.2.0.CO;2. Feng, M., A. Waite, and P. A. Thompson, 2009: Climate variability and ocean production in the Leeuwin Current System off the west coast of Western Australia. J. Roy. Soc. West. Aust., 92, 67–81. ——, M. J. McPhaden, S.-P. Xie, and J. Hafner, 2013: La Niña forces unprecedented Leeuwin Current warming in 2011. Sci. Rep., 3, 1277, doi:10.1038/srep01277. Ferré, B., K. Guizien, X. Durrieu de Madron, A. Palanques, J. Guillén, and A. Grémare, 2005: Fine-grained sediment dynamics during a strong storm event in the inner-shelf of the

3507

Gulf of Lion (NW Mediterranean). Cont. Shelf Res., 25, 2410– 2427, doi:10.1016/j.csr.2005.08.017. Gallop, S. L., F. Verspecht, and C. B. Pattiaratchi, 2012: Sea breezes drive currents on the inner continental shelf off southwest Western Australia. Ocean Dyn., 62, 569–583, doi:10.1007/s10236-011-0515-3. Gersbach, G. H., C. B. Pattiaratchi, G. N. Ivey, and G. R. Cresswell, 1999: Upwelling on the south-west coast of Australia—Source of the Capes Current? Cont. Shelf Res., 19, 363–400, doi:10.1016/S0278-4343(98)00088-0. Gonella, J., 1972: A rotary-component method for analysing meteorological and oceanographic vector time series. Deep-Sea Res. Oceanogr. Abstr., 19, 833–846, doi:10.1016/0011-7471(72)90002-2. Gough, M. K., A. J. H. M. Reniers, J. H. MacMahan, and S. D. Howden, 2016: Resonant near-surface inertial oscillations in the northeastern Gulf of Mexico. J. Geophys. Res. Oceans, 121, 2163–2182, doi:10.1002/2015JC011372. Grinsted, A., J. C. Moore, and S. Jevrejeva, 2004: Application of the cross wavelet transform and wavelet coherence to geophysical time series. Nonlinear Processes Geophys., 11, 561– 566, doi:10.5194/npg-11-561-2004. Gurgel, K.-W., G. Antonischki, H.-H. Essen, and T. Schlick, 1999: Wellen Radar (WERA): A new ground-wave HF radar for ocean remote sensing. Coastal Eng., 37, 219–234, doi:10.1016/ S0378-3839(99)00027-7. Hammond, T. M., C. B. Pattiaratchi, D. Eccles, M. J. Osborne, L. A. Nash, and M. B. Collins, 1987: Ocean Surface Current Radar (OSCR) vector measurements on the inner continental shelf. Cont. Shelf Res., 7, 411–431, doi:10.1016/0278-4343(87)90108-7. Hanson, C. E., C. B. Pattiaratchi, and A. M. Waite, 2005: Seasonal production regimes off south-western Australia: Influence of the Capes and Leeuwin Currents on phytoplankton dynamics. Mar. Freshwater Res., 56, 1011–1026, doi:10.1071/MF04288. Heron, M. L., and A. Prytz, 2011: The data archive for the phased array HF radars in the Australian Coastal Ocean Radar Network. OCEANS 2011 IEEE - Spain, Santander, Spain, Institute of Electrical and Electronics Engineers, 12207384, doi:10.1109/Oceans-Spain.2011.6003440. Hopkins, J. E., G. R. Stephenson, J. A. M. Green, M. E. Inall, and M. R. Palmer, 2014: Storms modify baroclinic energy fluxes in a seasonally stratified shelf sea: Inertial-tidal interaction. J. Geophys. Res. Oceans, 119, 6863–6883, doi:10.1002/2014JC010011. Hunter, E., R. Chant, L. Bowers, S. Glenn, and J. Kohut, 2007: Spatial and temporal variability of diurnal wind forcing in the coastal ocean. Geophys. Res. Lett., 34, L03607, doi:10.1029/2006GL028945. Hyder, P., J. H. Simpson, and S. Christopoulos, 2002: Sea-breeze forced diurnal surface currents in the Thermaikos Gulf, northwest Aegean. Cont. Shelf Res., 22, 585–601, doi:10.1016/ S0278-4343(01)00080-2. ——, ——, J. Xing, and S. T. Gille, 2011: Observations over an annual cycle and simulations of wind-forced oscillations near the critical latitude for diurnal–inertial resonance. Cont. Shelf Res., 31, 1576–1591, doi:10.1016/j.csr.2011.06.001. Jacob, S. D., L. K. Shay, and A. J. Mariano, 2000: The 3D oceanic mixed layer response to Hurricane Gilbert. J. Phys. Oceanogr., 30, 1407–1429, doi:10.1175/1520-0485(2000)030,1407: TOMLRT.2.0.CO;2. Jarosz, E., Z. R. Hallock, and W. J. Teague, 2007: Near-inertial currents in the DeSoto Canyon region. Cont. Shelf Res., 27, 2407–2426, doi:10.1016/j.csr.2007.06.014. Kaplan, D. M., J. Largier, and L. W. Botsford, 2005: HF radar observations of surface circulation off Bodega Bay (northern

3508

JOURNAL OF PHYSICAL OCEANOGRAPHY

California, USA). J. Geophys. Res., 110, C10020, doi:10.1029/ 2005JC002959. Kataoka, T., T. Tozuka, S. Behera, and T. Yamagata, 2014: On the Ningaloo Niño/Niña. Climate Dyn., 43, 1463–1482, doi:10.1007/ s00382-013-1961-z. Kim, S. Y., and G. Crawford, 2014: Resonant ocean current responses driven by coastal winds near the critical latitude. Geophys. Res. Lett., 41, 5581–5587, doi:10.1002/2014GL060402. Kundu, P. K., 1990: Fluid Mechanics. Academic Press, 638 pp. Kunze, E., 1985: Near-inertial wave propagation in geostrophic shear. J. Phys. Oceanogr., 15, 544–565, doi:10.1175/ 1520-0485(1985)015,0544:NIWPIG.2.0.CO;2. Lerczak, J. A., M. C. Hendershott, and C. D. Winant, 2001: Observations and modeling of coastal internal waves driven by a diurnal sea breeze. J. Geophys. Res., 106, 19 715–19 729, doi:10.1029/2001JC000811. Liu, Y., X. S. Liang, and R. H. Weisberg, 2007: Rectification of the bias in the wavelet power spectrum. J. Atmos. Oceanic Technol., 24, 2093–2102, doi:10.1175/2007JTECHO511.1. Lucas, A. J., G. C. Pitcher, T. A. Probyn, and R. M. Kudela, 2014: The influence of diurnal winds on phytoplankton dynamics in a coastal upwelling system off southwestern Africa. DeepSea Res. II, 101, 50–62, doi:10.1016/j.dsr2.2013.01.016. Masselink, G., and C. B. Pattiaratchi, 2001: Characteristics of the sea breeze system in Perth, Western Australia, and its effect on the nearshore wave climate. J. Coastal Res., 17, 173–187. [Available online at http://www.jstor.org/stable/4300161.] Monterey, G., and S. Levitus, 1997: Seasonal Variability of Mixed Layer Depth for the World Ocean. NOAA Atlas NESDIS 14, 102 pp. Nam, S., and U. Send, 2013: Resonant diurnal oscillations and mean alongshore flows driven by sea/land breeze forcing in the coastal Southern California Bight. J. Phys. Oceanogr., 43, 616– 630, doi:10.1175/JPO-D-11-0148.1. Orlic´, M., and Z. Pasaric´, 2011: A simple analytical model of periodic coastal upwelling. J. Phys. Oceanogr., 41, 1271–1276, doi:10.1175/JPO-D-10-05000.1. Paduan, J. D., and L. Washburn, 2013: High-frequency radar observations of ocean surface currents. Annu. Rev. Mar. Sci., 5, 115–136, doi:10.1146/annurev-marine-121211-172315. Pattiaratchi, C. B., and S. J. Buchan, 1991: Implications of longterm climate change for the Leeuwin Current. J. Roy. Soc. West. Aust., 74, 133–140. ——, and M. Eliot, 2009: Sea level variability in south-west Australia: From hours to decades. Coastal Engineering: Proceedings of the Thirty-First International Coastal Engineering Conference, J. McKee Smith, Ed., World Scientific, 1186–1198. ——, and M. Woo, 2009: The mean state of the Leeuwin Current system between North West Cape and Cape Leeuwin. J. Roy. Soc. West. Aust., 92, 221–241. ——, B. Hegge, J. Gould, and I. Eliot, 1997: Impact of sea-breeze activity on nearshore and foreshore processes in southwestern Australia. Cont. Shelf Res., 17, 1539–1560, doi:10.1016/ S0278-4343(97)00016-2. Pearce, A., and C. B. Pattiaratchi, 1999: The Capes Current: A summer countercurrent flowing past Cape Leeuwin and Cape Naturaliste, Western Australia. Cont. Shelf Res., 19, 401–420, doi:10.1016/S0278-4343(98)00089-2. Rennie, S. J., R. D. McCauley, and C. B. Pattiaratchi, 2006: Thermal structure above the Perth Canyon reveals Leeuwin Current, undercurrent and weather influences and the potential for upwelling. Mar. Freshwater Res., 57, 849–861, doi:10.1071/ MF05247.

VOLUME 46

Rosenfeld, L. K., 1988: Diurnal period wind stress and current fluctuations over the continental shelf off northern California. J. Geophys. Res., 93, 2257–2276, doi:10.1029/JC093iC03p02257. Sharples, J., and P. Tett, 1994: Modelling the effect of physical variability on the midwater chlorophyll maximum. J. Mar. Res., 52, 219–238, doi:10.1357/0022240943077109. ——, C. M. Moore, T. P. Rippeth, P. M. Holligan, D. J. Hydes, N. R. Fisher, and J. H. Simpson, 2001: Phytoplankton distribution and survival in the thermocline. Limnol. Oceanogr., 46, 486–496, doi:10.4319/lo.2001.46.3.0486. Shay, L. K., J. Martinez-Pedraja, T. M. Cook, B. K. Haus, and R. H. Weisberg, 2007: High-frequency radar mapping of surface currents using WERA. J. Atmos. Oceanic Technol., 24, 484– 503, doi:10.1175/JTECH1985.1. Simmons, H. L., 2008: Spectral modification and geographic redistribution of the semi-diurnal internal tide. Ocean Modell., 21, 126–138, doi:10.1016/j.ocemod.2008.01.002. Simpson, J. H., P. Hyder, and T. P. Rippeth, 2002: Forced oscillations near the critical latitude for diurnal-inertial resonance. J. Phys. Oceanogr., 32, 177–187, doi:10.1175/1520-0485(2002)032,0177: FONTCL.2.0.CO;2. Smith, R. L., A. Huyer, J. S. Godfrey, and J. A. Church, 1991: The Leeuwin Current off Western Australia, 1986–1987. J. Phys. Oceanogr., 21, 323–345, doi:10.1175/1520-0485(1991)021,0323: TLCOWA.2.0.CO;2. Sutherland, B. R., 2010: Internal Gravity Waves. Cambridge University Press, 377 pp. Thompson, P. A., K. Wild-Allen, M. Lourey, C. Rousseaux, A. M. Waite, M. Feng, and L. E. Beckley, 2011: Nutrients in an oligotrophic boundary current: Evidence of a new role for the Leeuwin Current. Prog. Oceanogr., 91, 345–359, doi:10.1016/ j.pocean.2011.02.011. Torrence, C., and G. P. Compo, 1998: A practical guide to wavelet analysis. Bull. Amer. Meteor. Soc., 79, 61–78, doi:10.1175/ 1520-0477(1998)079,0061:APGTWA.2.0.CO;2. Veleda, D., R. Montagne, and M. Araujo, 2012: Cross-wavelet bias corrected by normalizing scales. J. Atmos. Oceanic Technol., 29, 1401–1408, doi:10.1175/JTECH-D-11-00140.1. Verspecht, F., and C. B. Pattiaratchi, 2010: On the significance of wind event frequency for particulate resuspension and light attenuation in coastal waters. Cont. Shelf Res., 30, 1971–1982, doi:10.1016/j.csr.2010.09.008. Wolanski, E., 1986: A simple analytical model of low-frequency wind-driven upwelling on a continental slope. J. Phys. Oceanogr., 16, 1694–1702, doi:10.1175/1520-0485(1986)016,1694: ASAMOL.2.0.CO;2. Zaker, N. H., J. Imberger, and C. B. Pattiaratchi, 2007: Dynamics of the coastal boundary layer off Perth, Western Australia. J. Coastal Res., 23, 1112–1130, doi:10.2112/04-0374.1. Zhang, S., L. Xie, Y. Hou, H. Zhao, Y. Qi, and X. Yi, 2014: Tropical storm-induced turbulent mixing and chlorophyll-a enhancement in the continental shelf southeast of Hainan Island. J. Mar. Syst., 129, 405–414, doi:10.1016/j.jmarsys.2013.09.002. Zhang, X., S. F. Di Marco, D. C. Smith IV, M. K. Howard, and R. D. Hetland, 2009: Near-resonant ocean response to sea breeze on a stratified continental shelf. J. Phys. Oceanogr., 39, 2137–2155, doi:10.1175/2009JPO4054.1. ——, D. C. Smith IV, S. F. Di Marco, and R. D. Hetland, 2010: A numerical study of sea-breeze-driven ocean Poincare wave propagation and mixing near the critical latitude. J. Phys. Oceanogr., 40, 48–66, doi:10.1175/2009JPO4216.1.