Submesoscale wake formation at Green Island (∼7 km) in the Kuroshio is examined by the three‐dimensional numerical simulations, which are validated by field observations. On the basis of geophysical (rotating and stratified) flow, the wake exhibits sequentially detached recirculation, containing upwelling of cold water, propagates downstream via advection, forming an along‐stream oscillating wake, resembling to the von K árm án vortex streets (VKVS). Evidence includes (1) the shedding frequency as a function of the horizontal eddy viscosity shows a trend analogous with classical wakes; (2) the wake behaviors depend on the Reynolds number (Re), where the turbulent transition regime is determined; and (3) the aspect ratio of the island wakes is similar to the ratio of the VKVS. Unlike classical wakes, the vortex street features are adapted by inertial and barotropic instabilities. The inertial instability has large growth rate and tends to slightly destabilize the anticyclonic recirculation. The barotropic instability could be a secondary process to generate eddy kinetic energy at downstream. Finally, our model suggests the hotspot of the turbulent mixing in the wake is located at the plane free shear layer as a result of the vertical shear instability, which is induced by the island‐shelf effect and the tilting of the vertical vorticity.

Oceanic vortex evolution on the lee side of Taiwan’s Green Island (~7 km in diameter), where the Kuroshio flows at a speed of 1–1.5 m s−1, is observationally examined and compared to theories and the preceding results of laboratory experiments. In the near wake, recirculation occurs with a relative vorticity of ζ ~ 20 f (where f is the planetary vorticity) and subsequently sheds at a combination of periods resulting from the tidal oscillations and the intrinsic time scale of eddy evolution. The tidal oscillations are the predominant processes. Our analysis suggests that an island positioned in the Kuroshio with periodic and cross-stream tidal excursions is analogous to a cross-stream oscillating cylinder. Consequently, the shedding period of the vortex is synchronized to a tidal period occurring close to the intrinsic period. The free shear layer, which is characterized by an ~30 f relative vorticity band (2 km wide) and a wavy thermal front, develops between the Kuroshio and recirculation. The frontal wave occurring over a time period of 0.5–2 h resembles Kelvin–Helmholtz instability corresponding to high Re values. For the far wake, repeated cross-wake surveys suggest that cyclonic and anticyclonic vortices are alternatively present at a period close to the period of M2 tides in agreement with near-wake measurements. Repeated along-wake surveys reveal a cyclonic eddy shedding downstream at a speed of 0.35 m s−1, 1/3 of the upstream current speed, from the near wake. In comparing our observations with the results of previous water tank experiments, an Re value of O(103) for the submesoscale wake regime is expected.

In a stably stratified shear layer, multiple competing instabilities produce sensitivity to small changes in initial conditions, popularly called the butterfly effect (as a flapping wing may alter the weather). Three ensembles of 15 simulated mixing events, identical but for small perturbations to the initial state, are used to explore differences in the route to turbulence, the maximum turbulence level and the total amount and efficiency of mixing accomplished by each event. Comparisons show that a small change in the initial state alters the strength and timing of the primary Kelvin–Helmholtz instability, the subharmonic pairing instability and the various three-dimensional secondary instabilities that lead to turbulence. The effect is greatest in, but not limited to, the parameter regime where pairing and the three-dimensional secondary instabilities are in strong competition. Pairing may be accelerated or prevented; maximum turbulence kinetic energy may vary by up to a factor of 4.6, flux Richardson number by 12 %–15 % and net mixing by a factor of 2.

Studies of Kelvin–Helmholtz (KH) instability have typically modelled the initial flow as an isolated shear layer. In geophysical cases, however, the instability often occurs near boundaries and may therefore be influenced by boundary proximity effects. Ensembles of direct numerical simulations are conducted to understand the effect of boundary proximity on the evolution of the instability and the resulting turbulence. Ensemble averages are used to reduce sensitivity to small variations in initial conditions. Both the transition to turbulence and the resulting turbulent mixing are modified when the shear layer is near a boundary: the time scales for the onset of instability and turbulence are longer, and the height of the KH billow is reduced. Subharmonic instability is suppressed by the boundary because phase lock is prevented due to the diverging phase speeds of the KH and subharmonic modes. In addition, the disruptive influence of three-dimensional secondary instabilities on pairing is more profound as the two events coincide more closely. When the shear layer is far from the boundary, the shear-aligned convective instability is dominant; however, secondary central-core instability takes over when the shear layer is close to the boundary, providing an alternate route for the transition to turbulence. Both the efficiency of the resulting mixing and the turbulent diffusivity are dramatically reduced by boundary proximity effects.

Studies of Kelvin-Helmholtz (KH) instability have typically modelled the initial flow as an isolated shear layer.In geophysical cases, however, the instability often occurs near boundaries and may therefore be influenced by boundary proximity effects. Ensembles of direct numerical simulations are conducted to understand the effect of boundary proximity on the evolution of the instability and the resulting turbulence. Ensemble averages are used to reduce sensitivity to small variations in initial conditions. Both the transition to turbulence and the resulting turbulent mixing are modified when the shear layer is near a boundary: the time scales for the onset of instability and turbulence are longer, and the height of the KH billow is reduced. Subharmonic instability is suppressed by the boundary because phase-lock is prevented due to the diverging phase speeds of the KH and subharmonic modes. In addition, the disruptive influence of three-dimensional secondary instabilities on pairing is more profound as the two events coincide more closely. When the shear layer is far from the boundary, the shear-aligned convective instability is dominant; however, secondary central core instability takes over when the shear layer is close to the boundary, providing an alternate route for the transition to turbulence. Both the efficiency of the resulting mixing and the turbulent diffusivity are dramatically reduced by boundary proximity effects.