Browse
Abstract
Pressure anomaly set by the open ocean affects the dynamic topography and associated circulation over the continental shelf, which is explored here on a linearized β-plane arrested topographic wave framework that considers the variation in Coriolis parameter with latitude. It was found that on a meridional shelf, a nondimensional parameter Pe β , termed the β Péclet number, signifies the characteristics of open ocean–shelf interaction. The Pe β ≡ D β /α is determined by the ratio of long-wave-limit planetary to topographic Rossby wave speeds, i.e., the β drift D β , and the linear Ekman number α. On the western boundary shelf, due to the westward planetary Rossby wave, open ocean pressure propagates shoreward as Pe β > 1, and shelf circulation peaks where Pe β drops to 1. At this location, the planetary β effect is balanced by the bottom friction. The Pe β = 1 must occur either on the shelf or on the coastal wall when Pe β > 1 is observed at the shelf edge. On the eastern boundary shelf, however, Pe β < 0, the pressure anomaly is removed from the shelf, and hence the inductive circulation decays rapidly from the shelf edge. This β effect is robust on gently sloping meridional shelves. For zonal shelves, the planetary β increases the effective bottom slope on the northern boundary shelf but decreases it on the southern one, in a sense of potential vorticity conservation. However, this effect could be less significant in reality, given the complex dynamics involved. The above mechanism can explain the dynamics driving the Taiwan Warm Current in the East China Sea and its bifurcation around 28°N.
Abstract
Pressure anomaly set by the open ocean affects the dynamic topography and associated circulation over the continental shelf, which is explored here on a linearized β-plane arrested topographic wave framework that considers the variation in Coriolis parameter with latitude. It was found that on a meridional shelf, a nondimensional parameter Pe β , termed the β Péclet number, signifies the characteristics of open ocean–shelf interaction. The Pe β ≡ D β /α is determined by the ratio of long-wave-limit planetary to topographic Rossby wave speeds, i.e., the β drift D β , and the linear Ekman number α. On the western boundary shelf, due to the westward planetary Rossby wave, open ocean pressure propagates shoreward as Pe β > 1, and shelf circulation peaks where Pe β drops to 1. At this location, the planetary β effect is balanced by the bottom friction. The Pe β = 1 must occur either on the shelf or on the coastal wall when Pe β > 1 is observed at the shelf edge. On the eastern boundary shelf, however, Pe β < 0, the pressure anomaly is removed from the shelf, and hence the inductive circulation decays rapidly from the shelf edge. This β effect is robust on gently sloping meridional shelves. For zonal shelves, the planetary β increases the effective bottom slope on the northern boundary shelf but decreases it on the southern one, in a sense of potential vorticity conservation. However, this effect could be less significant in reality, given the complex dynamics involved. The above mechanism can explain the dynamics driving the Taiwan Warm Current in the East China Sea and its bifurcation around 28°N.
Abstract
The La Jolla Canyon System (LJCS) is a small, steep, shelf-incising canyon offshore of San Diego, California. Observations conducted in the fall of 2016 capture the dynamics of internal tides and turbulence patterns. Semidiurnal (D2) energy flux was oriented up-canyon; 62% ± 20% of the signal was contained in mode 1 at the offshore mooring. The observed mode-1 D2 tide was partly standing based on the ratio of group speed times energy c g E and energy flux F. Enhanced dissipation occurred near the canyon head at middepths associated with elevated strain arising from the standing wave pattern. Modes 2–5 were progressive, and energy fluxes associated with these modes were oriented down-canyon, suggesting that incident mode-1 waves were back-reflected and scattered. Flux integrated over all modes across a given canyon cross section was always onshore and generally decreased moving shoreward (from 240 ± 15 to 5 ± 0.3 kW), with a 50-kW increase in flux occurring on a section inshore of the canyon’s major bend, possibly due to reflection of incident waves from the supercritical sidewalls of the bend. Flux convergence from canyon mouth to head was balanced by the volume-integrated dissipation observed. By comparing energy budgets from a global compendium of canyons with sufficient observations (six in total), a similar balance was found. One exception was Juan de Fuca Canyon, where such a balance was not found, likely due to its nontidal flows. These results suggest that internal tides incident at the mouth of a canyon system are dissipated therein rather than leaking over the sidewalls or siphoning energy to other wave frequencies.
Abstract
The La Jolla Canyon System (LJCS) is a small, steep, shelf-incising canyon offshore of San Diego, California. Observations conducted in the fall of 2016 capture the dynamics of internal tides and turbulence patterns. Semidiurnal (D2) energy flux was oriented up-canyon; 62% ± 20% of the signal was contained in mode 1 at the offshore mooring. The observed mode-1 D2 tide was partly standing based on the ratio of group speed times energy c g E and energy flux F. Enhanced dissipation occurred near the canyon head at middepths associated with elevated strain arising from the standing wave pattern. Modes 2–5 were progressive, and energy fluxes associated with these modes were oriented down-canyon, suggesting that incident mode-1 waves were back-reflected and scattered. Flux integrated over all modes across a given canyon cross section was always onshore and generally decreased moving shoreward (from 240 ± 15 to 5 ± 0.3 kW), with a 50-kW increase in flux occurring on a section inshore of the canyon’s major bend, possibly due to reflection of incident waves from the supercritical sidewalls of the bend. Flux convergence from canyon mouth to head was balanced by the volume-integrated dissipation observed. By comparing energy budgets from a global compendium of canyons with sufficient observations (six in total), a similar balance was found. One exception was Juan de Fuca Canyon, where such a balance was not found, likely due to its nontidal flows. These results suggest that internal tides incident at the mouth of a canyon system are dissipated therein rather than leaking over the sidewalls or siphoning energy to other wave frequencies.
Abstract
As part of the Flow Encountering Abrupt Topography (FLEAT) program, an array of pressure-sensor equipped inverted echo sounders (PIESs) was deployed north of Palau where the westward-flowing North Equatorial Current encounters the southern end of the Kyushu–Palau Ridge in the tropical North Pacific. Capitalizing on concurrent observations from satellite altimetry, FLEAT Spray gliders, and shipboard hydrography, the PIESs’ 10-month duration hourly bottom pressure p and round-trip acoustic travel time τ records are used to examine the magnitude and predictability of sea level and pycnocline depth changes and to track signal propagations through the array. Sea level and pycnocline depth are found to vary in response to a range of ocean processes, with their magnitude and predictability strongly process dependent. Signals characterized here comprise the barotropic tides, semidiurnal and diurnal internal tides, southeastward-propagating superinertial waves, westward-propagating mesoscale eddies, and a strong signature of sea level increase and pycnocline deepening associated with the region’s relaxation from El Niño to La Niña conditions. The presence of a broad band of superinertial waves just above the inertial frequency was unexpected and the FLEAT observations and output from a numerical model suggest that these waves detected near Palau are forced by remote winds east of the Philippines. The PIES-based estimates of pycnocline displacement are found to have large uncertainties relative to overall variability in pycnocline depth, as localized deep current variations arising from interactions of the large-scale currents with the abrupt topography around Palau have significant travel time variability.
Abstract
As part of the Flow Encountering Abrupt Topography (FLEAT) program, an array of pressure-sensor equipped inverted echo sounders (PIESs) was deployed north of Palau where the westward-flowing North Equatorial Current encounters the southern end of the Kyushu–Palau Ridge in the tropical North Pacific. Capitalizing on concurrent observations from satellite altimetry, FLEAT Spray gliders, and shipboard hydrography, the PIESs’ 10-month duration hourly bottom pressure p and round-trip acoustic travel time τ records are used to examine the magnitude and predictability of sea level and pycnocline depth changes and to track signal propagations through the array. Sea level and pycnocline depth are found to vary in response to a range of ocean processes, with their magnitude and predictability strongly process dependent. Signals characterized here comprise the barotropic tides, semidiurnal and diurnal internal tides, southeastward-propagating superinertial waves, westward-propagating mesoscale eddies, and a strong signature of sea level increase and pycnocline deepening associated with the region’s relaxation from El Niño to La Niña conditions. The presence of a broad band of superinertial waves just above the inertial frequency was unexpected and the FLEAT observations and output from a numerical model suggest that these waves detected near Palau are forced by remote winds east of the Philippines. The PIES-based estimates of pycnocline displacement are found to have large uncertainties relative to overall variability in pycnocline depth, as localized deep current variations arising from interactions of the large-scale currents with the abrupt topography around Palau have significant travel time variability.
Abstract
Hydrographic measurements recently acquired along the thalweg of the Lifamatola Passage combined with historical moored velocity measurements immediately downstream of the sill are used to study the hydraulics, transport, mixing, and entrainment in the dense overflow. The observations suggest that the mean overflow is nearly critical at the mooring site, suggesting that a weir formula may be appropriate for estimating the overflow transport. Our assessment suggests that the weir formulas corresponding to a rectangular, triangular, or parabolic cross section all result in transports very close to the observation, suggesting their potential usage in long-term monitoring of the overflow transport or parameterizing the transport in numerical models. Analyses also suggest that deep signals within the overflow layer are blocked by the shear flow from propagating upstream, whereas the shallow wave modes of the full-depth continuously stratified flow are able to propagate upstream from the Banda Sea into the Maluku Sea. Strong mixing is found immediately downstream of the sill crest, with Thorpe-scale-based estimates of the mean dissipation rate within the overflow up to 1.1 × 10−7 W kg−1 and the region-averaged diapycnal diffusivity within the downstream overflow in the range of 2.3 × 10−3 to 10.1 × 10−3 m2 s−1. Mixing in the Lifamatola Passage results in 0.6–1.2-Sv (1 Sv ≡ 106 m3 s−1) entrainment transport added to the overflow, enhancing the deep-water renewal in the Banda Sea. A bulk diffusivity coefficient estimated in the deep Banda Sea yields 1.6 × 10−3 ± 5 × 10−4 m2 s−1, with an associated downward turbulent heat flux of 9 W m−2.
Abstract
Hydrographic measurements recently acquired along the thalweg of the Lifamatola Passage combined with historical moored velocity measurements immediately downstream of the sill are used to study the hydraulics, transport, mixing, and entrainment in the dense overflow. The observations suggest that the mean overflow is nearly critical at the mooring site, suggesting that a weir formula may be appropriate for estimating the overflow transport. Our assessment suggests that the weir formulas corresponding to a rectangular, triangular, or parabolic cross section all result in transports very close to the observation, suggesting their potential usage in long-term monitoring of the overflow transport or parameterizing the transport in numerical models. Analyses also suggest that deep signals within the overflow layer are blocked by the shear flow from propagating upstream, whereas the shallow wave modes of the full-depth continuously stratified flow are able to propagate upstream from the Banda Sea into the Maluku Sea. Strong mixing is found immediately downstream of the sill crest, with Thorpe-scale-based estimates of the mean dissipation rate within the overflow up to 1.1 × 10−7 W kg−1 and the region-averaged diapycnal diffusivity within the downstream overflow in the range of 2.3 × 10−3 to 10.1 × 10−3 m2 s−1. Mixing in the Lifamatola Passage results in 0.6–1.2-Sv (1 Sv ≡ 106 m3 s−1) entrainment transport added to the overflow, enhancing the deep-water renewal in the Banda Sea. A bulk diffusivity coefficient estimated in the deep Banda Sea yields 1.6 × 10−3 ± 5 × 10−4 m2 s−1, with an associated downward turbulent heat flux of 9 W m−2.
Abstract
Ocean lee waves occur on length scales that are smaller than the grid scale of global circulation models (GCMs). Therefore, such models must parameterize the drag associated with launching lee waves. This paper compares the lee wave drag predicted by existing parameterizations with the drag measured in high-resolution nonhydrostatic numerical simulations of a lee wave over periodic sinusoidal bathymetry. The simulations afford a time-varying glimpse at the nonlinear and nonhydrostatic oceanic lee wave spinup process and identify a characteristic time scale to reach steady state. The maximum instantaneous lee wave drag observed during the spinup period is found to be well predicted by linear lee wave theory for all hill heights. In steady state, the simulations demonstrate the applicability of parameterizing the drag based on applying linear theory to the lowest overtopping streamline of the flow (LOTS), as is currently employed in GCMs. However, because existing parameterizations are based only on the height of the LOTS, they implicitly assume hydrostatic flow. For hills tall enough to trap water in their valleys, the simulations identify a set of nonhydrostatic processes that can result in a reduction of the lee wave drag from that given by hydrostatic parameterizations. The simulations suggest implementing a time-dependent nonhydrostatic version of the LOTS-based parameterization of lee wave drag and demonstrate the remarkable applicability of linear lee wave theory to oceanic lee waves.
Abstract
Ocean lee waves occur on length scales that are smaller than the grid scale of global circulation models (GCMs). Therefore, such models must parameterize the drag associated with launching lee waves. This paper compares the lee wave drag predicted by existing parameterizations with the drag measured in high-resolution nonhydrostatic numerical simulations of a lee wave over periodic sinusoidal bathymetry. The simulations afford a time-varying glimpse at the nonlinear and nonhydrostatic oceanic lee wave spinup process and identify a characteristic time scale to reach steady state. The maximum instantaneous lee wave drag observed during the spinup period is found to be well predicted by linear lee wave theory for all hill heights. In steady state, the simulations demonstrate the applicability of parameterizing the drag based on applying linear theory to the lowest overtopping streamline of the flow (LOTS), as is currently employed in GCMs. However, because existing parameterizations are based only on the height of the LOTS, they implicitly assume hydrostatic flow. For hills tall enough to trap water in their valleys, the simulations identify a set of nonhydrostatic processes that can result in a reduction of the lee wave drag from that given by hydrostatic parameterizations. The simulations suggest implementing a time-dependent nonhydrostatic version of the LOTS-based parameterization of lee wave drag and demonstrate the remarkable applicability of linear lee wave theory to oceanic lee waves.
Abstract
Turbulent mixing adjacent to the Velasco Reef and Kyushu–Palau Ridge, off northern Palau in the western equatorial Pacific Ocean, is examined using shipboard and moored observations. The study focuses on a 9-day-long, ship-based microstructure and velocity survey, conducted in November–December 2016. Several sections (9–15 km in length) of microstructure, hydrographic, and velocity fields were acquired over and around the reef, where water depths ranged from 50 to 3000 m. Microstructure profiles were collected while steaming slowly either toward or away from the reef, and underway current surveys were conducted along quasi-rectangular boxes with side lengths of 5–10 km. Near the reef, both tidal and subtidal motions were important, while subtidal motions were stronger away from the reef. Vertical shears of currents and mixing were stronger on the northern and eastern flanks of the reef than on the western flanks. High turbulent kinetic energy dissipation rates, 10−6–10−4 W kg−1, and large values of eddy diffusivities, 10−4–10−2 m2 s−1, with strong turbulent heat fluxes, 100–500 W m−2, were found. Currents flowing along the eastern side separated at the northern tip of the reef and generated submesoscale cyclonic vorticity of about 2–4 times the planetary vorticity. The analysis suggests that a torque, imparted by the turbulent bottom stress, generated the cyclonic vorticity at the northern boundary. The northern reef is associated with high vertical transports resulting from both submesoscale flow convergences and energetic mixing. Even though the area around Palau represents a small footprint of the ocean, vertical velocities and mixing rates are several orders magnitude larger than in the open ocean.
Abstract
Turbulent mixing adjacent to the Velasco Reef and Kyushu–Palau Ridge, off northern Palau in the western equatorial Pacific Ocean, is examined using shipboard and moored observations. The study focuses on a 9-day-long, ship-based microstructure and velocity survey, conducted in November–December 2016. Several sections (9–15 km in length) of microstructure, hydrographic, and velocity fields were acquired over and around the reef, where water depths ranged from 50 to 3000 m. Microstructure profiles were collected while steaming slowly either toward or away from the reef, and underway current surveys were conducted along quasi-rectangular boxes with side lengths of 5–10 km. Near the reef, both tidal and subtidal motions were important, while subtidal motions were stronger away from the reef. Vertical shears of currents and mixing were stronger on the northern and eastern flanks of the reef than on the western flanks. High turbulent kinetic energy dissipation rates, 10−6–10−4 W kg−1, and large values of eddy diffusivities, 10−4–10−2 m2 s−1, with strong turbulent heat fluxes, 100–500 W m−2, were found. Currents flowing along the eastern side separated at the northern tip of the reef and generated submesoscale cyclonic vorticity of about 2–4 times the planetary vorticity. The analysis suggests that a torque, imparted by the turbulent bottom stress, generated the cyclonic vorticity at the northern boundary. The northern reef is associated with high vertical transports resulting from both submesoscale flow convergences and energetic mixing. Even though the area around Palau represents a small footprint of the ocean, vertical velocities and mixing rates are several orders magnitude larger than in the open ocean.
Abstract
Towed shipboard and moored observations show internal gravity waves over a tall, supercritical submarine ridge that reaches to 1000 m below the ocean surface in the tropical western Pacific north of Palau. The lee-wave or topographic Froude number, Nh
0/U
0 (where N is the buoyancy frequency, h
0 the ridge height, and U
0 the farfield velocity), ranged between 25 and 140. The waves were generated by a superposition of tidal and low-frequency flows and thus had two distinct energy sources with combined amplitudes of up to 0.2 m s−1. Local breaking of the waves led to enhanced rates of dissipation of turbulent kinetic energy reaching above 10−6 W kg−1 in the lee of the ridge near topography. Turbulence observations showed a stark contrast between conditions at spring and neap tide. During spring tide, when the tidal flow dominated, turbulence was approximately equally distributed around both sides of the ridge. During neap tide, when the mean flow dominated over tidal oscillations, turbulence was mostly observed on the downstream side of the ridge relative to the mean flow. The drag exerted by the ridge on the flow, estimated to
Abstract
Towed shipboard and moored observations show internal gravity waves over a tall, supercritical submarine ridge that reaches to 1000 m below the ocean surface in the tropical western Pacific north of Palau. The lee-wave or topographic Froude number, Nh
0/U
0 (where N is the buoyancy frequency, h
0 the ridge height, and U
0 the farfield velocity), ranged between 25 and 140. The waves were generated by a superposition of tidal and low-frequency flows and thus had two distinct energy sources with combined amplitudes of up to 0.2 m s−1. Local breaking of the waves led to enhanced rates of dissipation of turbulent kinetic energy reaching above 10−6 W kg−1 in the lee of the ridge near topography. Turbulence observations showed a stark contrast between conditions at spring and neap tide. During spring tide, when the tidal flow dominated, turbulence was approximately equally distributed around both sides of the ridge. During neap tide, when the mean flow dominated over tidal oscillations, turbulence was mostly observed on the downstream side of the ridge relative to the mean flow. The drag exerted by the ridge on the flow, estimated to
Abstract
The effects of topography on the linear stability of both barotropic vortices and two-layer, baroclinic vortices are examined by considering cylindrical topography and vortices with stepwise relative vorticity profiles in the quasigeostrophic approximation. Four vortex configurations are considered, classified by the number of relative vorticity steps in the horizontal and the number of layers in the vertical: barotropic one-step vortex (Rankine vortex), barotropic two-step vortex, and their two-layer, baroclinic counterparts with the vorticity steps in the upper layer. In the barotropic calculation, the vortex is destabilized by topography having an oppositely signed potential vorticity jump while stabilized by topography of same-signed jump, that is, anticyclones are destabilized by seamounts while stabilized by depressions. Further, topography of appropriate sign and magnitude can excite a mode-1 instability for a two-step vortex, especially relevant for topographic encounters of an otherwise stable vortex. The baroclinic calculation is in general consistent with the barotropic calculation except that the growth rate weakens and, for a two-step vortex, becomes less sensitive to topography (sign and magnitude) as baroclinicity increases. The smaller growth rate for a baroclinic vortex is consistent with previous findings that vortices with sufficient baroclinic structure could cross the topography relatively easily. Nonlinear contour dynamics simulations are conducted to confirm the linear stability analysis and to describe the subsequent evolution.
Abstract
The effects of topography on the linear stability of both barotropic vortices and two-layer, baroclinic vortices are examined by considering cylindrical topography and vortices with stepwise relative vorticity profiles in the quasigeostrophic approximation. Four vortex configurations are considered, classified by the number of relative vorticity steps in the horizontal and the number of layers in the vertical: barotropic one-step vortex (Rankine vortex), barotropic two-step vortex, and their two-layer, baroclinic counterparts with the vorticity steps in the upper layer. In the barotropic calculation, the vortex is destabilized by topography having an oppositely signed potential vorticity jump while stabilized by topography of same-signed jump, that is, anticyclones are destabilized by seamounts while stabilized by depressions. Further, topography of appropriate sign and magnitude can excite a mode-1 instability for a two-step vortex, especially relevant for topographic encounters of an otherwise stable vortex. The baroclinic calculation is in general consistent with the barotropic calculation except that the growth rate weakens and, for a two-step vortex, becomes less sensitive to topography (sign and magnitude) as baroclinicity increases. The smaller growth rate for a baroclinic vortex is consistent with previous findings that vortices with sufficient baroclinic structure could cross the topography relatively easily. Nonlinear contour dynamics simulations are conducted to confirm the linear stability analysis and to describe the subsequent evolution.
Abstract
Microstructure measurements in Drake Passage and on the flanks of Kerguelen Plateau find turbulent dissipation rates ε on average factors of 2–3 smaller than linear lee-wave generation predictions, as well as a factor of 3 smaller than the predictions of a well-established parameterization based on finescale shear and strain. Here, the possibility that these discrepancies are a result of conservation of wave action E/ω L = E/|kU| is explored. Conservation of wave action will transfer a fraction of the lee-wave radiation back to the mean flow if the waves encounter weakening currents U, where the intrinsic or Lagrangian frequency ω L = |kU| ↓ |f| and k the along-stream horizontal wavenumber, where kU ≡ k ⋅ V. The dissipative fraction of power that is lost to turbulence depends on the Doppler shift of the intrinsic frequency between generation and breaking, hence on the topographic height spectrum and bandwidth N/f. The partition between dissipation and loss to the mean flow is quantified for typical topographic height spectral shapes and N/f ratios found in the abyssal ocean under the assumption that blocking is local in wavenumber. Although some fraction of lee-wave generation is always dissipated in a rotating fluid, lee waves are not as large a sink for balanced energy or as large a source for turbulence as previously suggested. The dissipative fraction is 0.44–0.56 for topographic spectral slopes and buoyancy frequencies typical of the deep Southern Ocean, insensitive to flow speed U and topographic splitting. Lee waves are also an important mechanism for redistributing balanced energy within their generating bottom current.
Abstract
Microstructure measurements in Drake Passage and on the flanks of Kerguelen Plateau find turbulent dissipation rates ε on average factors of 2–3 smaller than linear lee-wave generation predictions, as well as a factor of 3 smaller than the predictions of a well-established parameterization based on finescale shear and strain. Here, the possibility that these discrepancies are a result of conservation of wave action E/ω L = E/|kU| is explored. Conservation of wave action will transfer a fraction of the lee-wave radiation back to the mean flow if the waves encounter weakening currents U, where the intrinsic or Lagrangian frequency ω L = |kU| ↓ |f| and k the along-stream horizontal wavenumber, where kU ≡ k ⋅ V. The dissipative fraction of power that is lost to turbulence depends on the Doppler shift of the intrinsic frequency between generation and breaking, hence on the topographic height spectrum and bandwidth N/f. The partition between dissipation and loss to the mean flow is quantified for typical topographic height spectral shapes and N/f ratios found in the abyssal ocean under the assumption that blocking is local in wavenumber. Although some fraction of lee-wave generation is always dissipated in a rotating fluid, lee waves are not as large a sink for balanced energy or as large a source for turbulence as previously suggested. The dissipative fraction is 0.44–0.56 for topographic spectral slopes and buoyancy frequencies typical of the deep Southern Ocean, insensitive to flow speed U and topographic splitting. Lee waves are also an important mechanism for redistributing balanced energy within their generating bottom current.
Abstract
In this study, a 2-yr time series of velocity profiles to 1000 m from meridional glider surveys is used to characterize the wake in the lee of a large island in the western tropical North Pacific Ocean, Palau. Surveys were completed along sections to the east and west of the island to capture both upstream and downstream conditions. Objectively mapped in time and space, mean sections of velocity show the incident westward North Equatorial Current accelerating around the island of Palau, increasing from 0.1 to 0.2 m s−1 at the surface. Downstream of the island, elevated velocity variability and return flow in the lee are indicative of boundary layer separation. Isolating for periods of depth-average westward flow reveals a length scale in the wake that reflects local details of the topography. Eastward flow is shown to produce an asymmetric wake. Depth-average velocity time series indicate that energetic events (on time scales from weeks to months) are prevalent. These events are associated with mean vorticity values in the wake up to 0.3f near the surface and with instantaneous values that can exceed f (the local Coriolis frequency) during periods of sustained, anomalously strong westward flow. Thus, ageostrophic effects become important to first order.
Abstract
In this study, a 2-yr time series of velocity profiles to 1000 m from meridional glider surveys is used to characterize the wake in the lee of a large island in the western tropical North Pacific Ocean, Palau. Surveys were completed along sections to the east and west of the island to capture both upstream and downstream conditions. Objectively mapped in time and space, mean sections of velocity show the incident westward North Equatorial Current accelerating around the island of Palau, increasing from 0.1 to 0.2 m s−1 at the surface. Downstream of the island, elevated velocity variability and return flow in the lee are indicative of boundary layer separation. Isolating for periods of depth-average westward flow reveals a length scale in the wake that reflects local details of the topography. Eastward flow is shown to produce an asymmetric wake. Depth-average velocity time series indicate that energetic events (on time scales from weeks to months) are prevalent. These events are associated with mean vorticity values in the wake up to 0.3f near the surface and with instantaneous values that can exceed f (the local Coriolis frequency) during periods of sustained, anomalously strong westward flow. Thus, ageostrophic effects become important to first order.