Browse
Abstract
Topographic form stress (TFS) plays a central role in constraining the transport of the Antarctic Circumpolar Current (ACC), and thus the rate of exchange between the major ocean basins. Topographic form stress generation in the ACC has been linked to the formation of standing Rossby waves, which occur because the current is retrograde (opposing the direction of Rossby wave propagation). However, it is unclear whether TFS similarly retards current systems that are prograde (in the direction of Rossby wave propagation), which cannot arrest Rossby waves. An isopycnal model is used to investigate the momentum balance of wind-driven prograde and retrograde flows in a zonal channel, with bathymetry consisting of either a single ridge or a continental shelf and slope with a meridional excursion. Consistent with previous studies, retrograde flows are almost entirely impeded by TFS, except in the limit of flat bathymetry, whereas prograde flows are typically impeded by a combination of TFS and bottom friction. A barotropic theory for standing waves shows that bottom friction serves to shift the phase of the standing wave’s pressure field from that of the bathymetry, which is necessary to produce TFS. The mechanism is the same in prograde and retrograde flows, but is most efficient when the mean flow arrests a Rossby wave with a wavelength comparable to that of the bathymetry. The asymmetry between prograde and retrograde momentum balances implies that prograde current systems may be more sensitive to changes in wind forcing, for example associated with climate shifts.
Abstract
Topographic form stress (TFS) plays a central role in constraining the transport of the Antarctic Circumpolar Current (ACC), and thus the rate of exchange between the major ocean basins. Topographic form stress generation in the ACC has been linked to the formation of standing Rossby waves, which occur because the current is retrograde (opposing the direction of Rossby wave propagation). However, it is unclear whether TFS similarly retards current systems that are prograde (in the direction of Rossby wave propagation), which cannot arrest Rossby waves. An isopycnal model is used to investigate the momentum balance of wind-driven prograde and retrograde flows in a zonal channel, with bathymetry consisting of either a single ridge or a continental shelf and slope with a meridional excursion. Consistent with previous studies, retrograde flows are almost entirely impeded by TFS, except in the limit of flat bathymetry, whereas prograde flows are typically impeded by a combination of TFS and bottom friction. A barotropic theory for standing waves shows that bottom friction serves to shift the phase of the standing wave’s pressure field from that of the bathymetry, which is necessary to produce TFS. The mechanism is the same in prograde and retrograde flows, but is most efficient when the mean flow arrests a Rossby wave with a wavelength comparable to that of the bathymetry. The asymmetry between prograde and retrograde momentum balances implies that prograde current systems may be more sensitive to changes in wind forcing, for example associated with climate shifts.
Abstract
Previous studies have concluded that the wind-input vorticity in ocean gyres is balanced by bottom pressure torques (BPT), when integrated over latitude bands. However, the BPT must vanish when integrated over any area enclosed by an isobath. This constraint raises ambiguities regarding the regions over which BPT should close the vorticity budget, and implies that BPT generated to balance a local wind stress curl necessitates the generation of a compensating, nonlocal BPT and thus nonlocal circulation. This study aims to clarify the role of BPT in wind-driven gyres using an idealized isopycnal model. Experiments performed with a single-signed wind stress curl in an enclosed, sloped basin reveal that BPT balances the winds only when integrated over latitude bands. Integrating over other, dynamically motivated definitions of the gyre, such as barotropic streamlines, yields a balance between wind stress curl and bottom frictional torques. This implies that bottom friction plays a nonnegligible role in structuring the gyre circulation. Nonlocal bottom pressure torques manifest in the form of along-slope pressure gradients associated with a weak basin-scale circulation, and are associated with a transition to a balance between wind stress and bottom friction around the coasts. Finally, a suite of perturbation experiments is used to investigate the dynamics of BPT. To predict the BPT, the authors extend a previous theory that describes propagation of surface pressure signals from the gyre interior toward the coast along planetary potential vorticity contours. This theory is shown to agree closely with the diagnosed contributions to the vorticity budget across the suite of model experiments.
Abstract
Previous studies have concluded that the wind-input vorticity in ocean gyres is balanced by bottom pressure torques (BPT), when integrated over latitude bands. However, the BPT must vanish when integrated over any area enclosed by an isobath. This constraint raises ambiguities regarding the regions over which BPT should close the vorticity budget, and implies that BPT generated to balance a local wind stress curl necessitates the generation of a compensating, nonlocal BPT and thus nonlocal circulation. This study aims to clarify the role of BPT in wind-driven gyres using an idealized isopycnal model. Experiments performed with a single-signed wind stress curl in an enclosed, sloped basin reveal that BPT balances the winds only when integrated over latitude bands. Integrating over other, dynamically motivated definitions of the gyre, such as barotropic streamlines, yields a balance between wind stress curl and bottom frictional torques. This implies that bottom friction plays a nonnegligible role in structuring the gyre circulation. Nonlocal bottom pressure torques manifest in the form of along-slope pressure gradients associated with a weak basin-scale circulation, and are associated with a transition to a balance between wind stress and bottom friction around the coasts. Finally, a suite of perturbation experiments is used to investigate the dynamics of BPT. To predict the BPT, the authors extend a previous theory that describes propagation of surface pressure signals from the gyre interior toward the coast along planetary potential vorticity contours. This theory is shown to agree closely with the diagnosed contributions to the vorticity budget across the suite of model experiments.
Abstract
Slowly evolving stratified flow over rough topography is subject to substantial drag due to internal motions, but often numerical simulations are carried out at resolutions where this “wave” drag must be parameterized. Here we highlight the importance of internal drag from topography with scales that cannot radiate internal waves, but may be highly nonlinear, and we propose a simple parameterization of this drag that has a minimum of fit parameters compared to existing schemes. The parameterization smoothly transitions from a quadratic drag law (
Abstract
Slowly evolving stratified flow over rough topography is subject to substantial drag due to internal motions, but often numerical simulations are carried out at resolutions where this “wave” drag must be parameterized. Here we highlight the importance of internal drag from topography with scales that cannot radiate internal waves, but may be highly nonlinear, and we propose a simple parameterization of this drag that has a minimum of fit parameters compared to existing schemes. The parameterization smoothly transitions from a quadratic drag law (
Abstract
An array of moorings deployed off the coast of Palau is used to characterize submesoscale vorticity generated by broadband upper-ocean flows around the island. Palau is a steep-sided archipelago lying in the path of strong zonal geostrophic currents, but tides and inertial oscillations are energetic as well. Vorticity is correspondingly broadband, with both mean and variance O(f) in a surface and subsurface layer (where f is the local Coriolis frequency). However, while subinertial vorticity is linearly related to the incident subinertial current, the relationship between superinertial velocity and superinertial vorticity is weak. Instead, there is a strong nonlinear relationship between subinertial velocity and superinertial vorticity. A key observation of this study is that during periods of strong westward flow, vorticity in the tidal bands increases by an order of magnitude. Empirical orthogonal functions (EOFs) of velocity show this nonstationary, superinertial vorticity variance is due to eddy motion at the scale of the array. Comparison of kinetic energy and vorticity time series suggest that lateral shear against the island varies with the subinertial flow, while tidal currents lead to flow reversals inshore of the recirculating wake and possibly eddy shedding. This is a departure from the idealized analog typically drawn on in island wake studies: a cylinder in a steady flow. In that case, eddy formation occurs at a frequency dependent on the scale of the obstacle and strength of the flow alone. The observed tidal formation frequency likely modulates the strength of submesoscale wake eddies and thus their dynamic relationship to the mesoscale wake downstream of Palau.
Abstract
An array of moorings deployed off the coast of Palau is used to characterize submesoscale vorticity generated by broadband upper-ocean flows around the island. Palau is a steep-sided archipelago lying in the path of strong zonal geostrophic currents, but tides and inertial oscillations are energetic as well. Vorticity is correspondingly broadband, with both mean and variance O(f) in a surface and subsurface layer (where f is the local Coriolis frequency). However, while subinertial vorticity is linearly related to the incident subinertial current, the relationship between superinertial velocity and superinertial vorticity is weak. Instead, there is a strong nonlinear relationship between subinertial velocity and superinertial vorticity. A key observation of this study is that during periods of strong westward flow, vorticity in the tidal bands increases by an order of magnitude. Empirical orthogonal functions (EOFs) of velocity show this nonstationary, superinertial vorticity variance is due to eddy motion at the scale of the array. Comparison of kinetic energy and vorticity time series suggest that lateral shear against the island varies with the subinertial flow, while tidal currents lead to flow reversals inshore of the recirculating wake and possibly eddy shedding. This is a departure from the idealized analog typically drawn on in island wake studies: a cylinder in a steady flow. In that case, eddy formation occurs at a frequency dependent on the scale of the obstacle and strength of the flow alone. The observed tidal formation frequency likely modulates the strength of submesoscale wake eddies and thus their dynamic relationship to the mesoscale wake downstream of Palau.
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.
