Barotropic Instability Research Articles

The Role of the Subtropical Jet in Deficient Winter Precipitation Across the Mid‐Holocene Indus Basin

AbstractThe mid‐Holocene (7–5 ka) was a period with an increased seasonal insolation cycle, resulting from decreased insolation during northern hemisphere winter. Here, a set of six CMIP5 models is used to show that the decreased insolation reduced the upper‐tropospheric meridional temperature gradient, producing a weaker subtropical jet with less horizontal shear. These effects work to reduce the baroclinic and barotropic instability available for perturbations to grow, and in consequence, storm‐tracking results show that there are fewer winter storms over India and Pakistan (known as western disturbances). These western disturbances are weaker, resulting in a reduction in winter precipitation of around 15% in the north Indus Basin. Combined with previous work showing greater northwestward extent of the Indian monsoon during the mid‐Holocene, our Global climate models‐derived results are consistent with the Indus Basin changing from a summer‐growing season in the mid‐Holocene to a winter‐growing season in the present day.

Barotropic Instability across the Moat and Inner Eyewall Dissipation: A Numerical Study of Hurricane Wilma (2005)

Abstract Radar imagery of some double-eyewall tropical cyclones shows that the inner eyewalls become elliptical prior to their dissipation. These elliptical features indicate that the barotropic instability (BI) across the moat (aka, type-2 BI) may play a role in the process. To investigate the mechanism for dissipation, a WRF simulation of Hurricane Wilma (2005) is performed. The results reveal an elliptical elongation of the inner eyewall and a change in the structure of the radial flow from wavenumber (WN) 1 to WN 2 at the lower levels. A linear stability analysis as well as idealized nonlinear experiments using a nondivergent barotropic vorticity model initialized with the vorticity fields before the change in the dominant wavenumber of the radial flow are presented with the results supporting the presence of a type-2 BI at the lower levels. The accompanying WN-2 radial flow is also found to dilute the vorticity within the inner eyewall and the eye. However, this dilution is not seen at higher levels as the type-2 BI becomes weak and short lived at the middle levels and reaches its weakest strength at the upper levels. This phenomenon is traced to the fact that a higher growth rate comes with a narrower moat for type-2 BI. As the outward slope of the outer eyewall is larger than that of the inner eyewall, the moat width increases with height so that the growth rate decreases with height. The results presented here thus highlight the potential role played by the barotropic instability across the moat in inner eyewall dissipation.

Optimal precursors of double-gyre regime transitions with an adjoint-free method

In this paper, we find the optimal precursors which can cause double-gyre regime transitions based on conditional nonlinear optimal perturbation (CNOP) method with Regional Ocean Modeling System (ROMS). Firstly, we simulate the multiple-equilibria regimes of double-gyre circulation under different viscosity coefficient and obtain the bifurcation diagram, then choose two equilibrium states (called jet-up state and jet-down state) as reference states respectively, propose Principal Component Analysisbased Simulated Annealing (PCASA) algorithm to solve CNOP-type initial perturbations which can induce double-gyre regime transitions between jet-up state and jet-down state. PCASA algorithm is an adjoint-free method which searches optimal solution randomly in the whole solution space. In addition, we investigate CNOP-type initial perturbations how to evolve with time. The experiments results show: (1) the CNOP-type perturbations present a two-cell structure, and gradually evolves into a three-cell structure at predictive time; (2) by superimposing CNOP-type perturbations on the jet-up state and integrating ROMS, doublegyre circulation transfers from jet-up state to jet-down state, and vice versa, and random initial perturbations don’t cause the transitions, which means CNOP-type perturbations are the optimal precursors of doublegyre regime transitions; (3) by analyzing the transition process of double-gyre regime transitions, we find that CNOP-type initial perturbations obtain energy from the background state through both barotropic and baroclinic instabilities, and barotropic instability contributes more significantly to the fast-growth of the perturbations. The optimal precursors and the dynamic mechanism of double-gyre regime transitions revealed in this paper have an important significance to enhance the predictability of double-gyre circulation.

Energetics of Eddy–Mean Flow Interactions along the Western Boundary Currents in the North Pacific

AbstractThe energetics of eddy–mean flow interactions along two western boundary currents of the North Pacific, the Kuroshio and Ryukyu Currents, are systematically investigated using 22 years of numerical data from the Ocean General Circulation Model for the Earth Simulator (OFES). For the time-mean and time-varying flow fields, all the energy components and conversions exhibit inhomogeneous spatial distributions. In the two currents, complex cross-stream and along-stream variations are seen in the eddy–mean flow energy conversions. East of Taiwan, the kinetic energy is mainly transferred from the mean flow to the eddy field through barotropic instability, whereas the baroclinic energy conversions form a meridional dipole structure caused by the topographic constraint. In the northern area, particularly, the eddy field drains 2.25 × 108 W of kinetic energy and releases 2.82 × 108 W of available potential energy when interacting with the mean flow, indicating that mesoscale eddies impinging on the Kuroshio decay with baroclinic inverse energy cascades. In the Ryukyu Current, inverse energy conversions from the eddy field to the mean flow also dominate the power transfer in the subsurface layer. The eddy field transfers 0.16 × 108 W of kinetic energy and 1.89 × 108 W of available potential energy to the mean flow, suggesting that meososcale eddies play an important role in maintaining the velocity and hydrographic structure of the current. In other areas, both barotropic and baroclinic instabilities contribute to the generation of eddy kinetic energy with the latter one providing more than 3 times as much power as the former one.

Observed Storm Track Dynamics in Drake Passage

AbstractThe dynamics of an oceanic storm track—where energy and enstrophy transfer between the mean flow and eddies—are investigated using observations from an eddy-rich region of the Antarctic Circumpolar Current downstream of the Shackleton Fracture Zone (SFZ) in Drake Passage. Four years of measurements by an array of current- and pressure-recording inverted echo sounders deployed between November 2007 and November 2011 are used to diagnose eddy–mean flow interactions and provide insight into physical mechanisms for these transfers. Averaged within the upper to mid-water column (400–1000-m depth) and over the 4-yr-record mean field, eddy potential energy is highest in the western part of the storm track and maximum eddy kinetic energy occurs farther away from the SFZ, shifting the proportion of eddy energies from to about 1 along the storm track. There are enhanced mean 3D wave activity fluxes immediately downstream of SFZ with strong horizontal flux vectors emanating northeast from this region. Similar patterns across composites of Polar Front and Subantarctic Front meander intrusions suggest the dynamics are set more so by the presence of the SFZ than by the eddy’s sign. A case study showing the evolution of a single eddy event, from 15 to 23 July 2010, highlights the storm-track dynamics in a series of snapshots. Consistently, explaining the eddy energetics pattern requires both horizontal and vertical components of W, implying the importance of barotropic and baroclinic processes and instabilities in controlling storm-track dynamics in Drake Passage.

Implementation of the Vector Vorticity Dynamical Core on Cubed Sphere for Use in the Quasi‐3‐D Multiscale Modeling Framework

AbstractThe dynamical core that predicts the three‐dimensional vorticity rather than the momentum, which is called Vector‐Vorticity Model (VVM), is implemented on a cubed sphere. Its horizontal coordinate system is not restricted to orthogonal, while the vertical coordinate is orthogonal to the horizontal surface. Accordingly, all the governing equations of the VVM, which are originally developed with Cartesian coordinates, are rewritten in terms of general curvilinear coordinates. The local coordinates on each cube surface are constructed with the gnomonic equiangular projection. Using global channel domains, the VVM on the cubed sphere has been evaluated by (1) advecting a passive tracer with a bell‐shaped initial perturbation along an east‐west latitude circle and along a north‐south meridional circle and (2) simulating the evolution of barotropic and baroclinic instabilities. The simulated results with the cubed‐sphere grids are compared to analytic solutions or those with the regular longitude‐latitude grids. The convergence with increasing spatial resolution is also quantified using standard error norms. The comparison shows that the solutions with the cubed‐sphere grids are quite reasonable for both linear and nonlinear problems when high resolutions are used. With coarse resolution, degeneracy appears in the solutions of the nonlinear problems such as spurious wave growth; however, it is effectively reduced with increased resolution. Based on the encouraging results in this study, we intend to use this model as the cloud‐resolving component in a global Quasi‐Three‐Dimensional Multiscale Modeling Framework.

Jet Instability over Smooth, Corrugated, and Realistic Bathymetry

AbstractThe stability of a horizontally and vertically sheared surface jet is examined, with a focus on the vertical structure of the resultant eddies. Over a flat bottom, the instability is mixed baroclinic/barotropic, producing strong eddies at depth that are characteristically shifted downstream relative to the surface eddies. Baroclinic instability is suppressed over a large slope for retrograde jets (with a flow antiparallel to topographic wave propagation) and to a lesser extent for prograde jets (with flow parallel to topographic wave propagation), as seen previously. In such cases, barotropic (lateral) instability dominates if the jet is sufficiently narrow. This yields surface eddies whose size is independent of the slope but proportional to the jet width. Deep eddies still form, forced by interfacial motion associated with the surface eddies, but they are weaker than under baroclinic instability and are vertically aligned with the surface eddies. A sinusoidal ridge acts similarly, suppressing baroclinic instability and favoring lateral instability in the upper layer. A ridge with a 1-km wavelength and an amplitude of roughly 10 m is sufficient to suppress baroclinic instability. Surveys of bottom roughness from bathymetry acquired with shipboard multibeam echo sounding reveal that such heights are common beneath the Kuroshio, the Antarctic Circumpolar Current, and, to a lesser extent, the Gulf Stream. Consistent with this, vorticity and velocity cross sections from a 1/50° HYCOM simulation suggest that Gulf Stream eddies are vertically aligned, as in the linear stability calculations with strong topography. Thus, lateral instability may be more common than previously thought, owing to topography hindering vertical energy transfer.

Barotropic growth of monsoon depressions

Although monsoon depressions are a principal synoptic‐scale element of the South Asian monsoon, producing extreme rainfall over India and surrounding regions, there exists no widely accepted mechanism explaining their occurrence. This study presents a hierarchy of numerical experiments aimed at finding such an explanation. Using a perturbation‐basic state decomposition, we derive an anelastic system of equations that can represent disturbances growing in the complex, three‐dimensional monsoon basic state. We find that modal solutions to these equations linearized about this basic state can explain many features of observed monsoon depressions, including their warm‐over‐cold core structure, westward propagation, and lower‐tropospheric wind maximum. For the zonally symmetric case, these modes are barotropically unstable, drawing energy from the meridional shear of the monsoon trough. For the zonally varying basic state, modal solutions still derive energy from barotropic conversion, but fail to achieve positive net growth rates when dissipative processes are included. For the nonlinear equation set, these modes can be excited by a heating impulse, and their energy then remains roughly constant over several days as barotropic energy transfers oppose dissipative losses. Our results support the idea that the general concept of barotropic instability can explain the structure, propagation, and geographic distribution of monsoon depressions, but not their rapid growth rates. We speculate that condensational heating coupled to these waves is needed to obtain a positive growth rate.

A Multiscale Energetic Diagnosis of the Response of Mokpo Sea to Typhoon Bolaven

The ocean response to typhoon is usually characterized by a cooling on the sea surface. In August 2012, however, a warming (instead of cooling) event occurs in the Yellow Sea outside Mokpo, South Korea, as the typhoon Bolaven (2012) passes. This study gives a brief introduction to this abnormal sea surface warming. It also provides a multiscale energetic diagnosis of the oceanic response to Typhoon Bolaven. We used a recently developed analysis tool named “multiscale window transform” (MWT). Based on the MWT, we also expanded a localized multiscale energy and vorticity analysis (MS-EVA). The fields are reconstructed onto three scale windows: large-scale, abnormal warming-scale, and high frequency tide-scale windows. The results show that the kinetic energy (KE) in the abnormal warming-scale window of the Mokpo area is obviously enhanced during the passage of Bolaven, which can be attributed to three processes: transfer, transport process of KE and wind stress work. At the same time, the large-scale window in the Mokpo area experiences barotropic instabilities with KE transfers from large-scale window to warming-scale window. Besides, the strong wind stress bought by the passage of Bolaven not only inputs a large amount of KE into warming-scale window, but also causes the increase of KE flux convergence.

Revisiting the Charney Baroclinic Instability Problem and Point-jet Barotropic Instability Problem, Part II: Matched Asymptotic Expansions & Overreflection Without Delta-Functions

Abstract Baroclinic instability generates the cyclones and anticyclones of midlatitude weather. Charney developed the first effective theory for the infancy of this cyclogenesis in 1947. His linear eigenproblem is analytically solvable by confluent hypergeometric functions. It is also, with extension of the domain of the coordinate from [0,∞] to [−∞,∞] by reflection about the origin, the point-jet model of barotropic instability, important for tropical cyclogenesis. (Note that the coordinate is height z in the Charney model, but latitude y for the point-jet bartropic instability. It is a great simplification that the Charney and point-jet instability problems are mathematically identical, but it also is confusing that the mathematical analysis in y also applies to the Charney problem with the substitution of z for y.) Unfortunately, the theory is full of distributions like the Dirac delta-function and the reflected Charney eigenfunction has a discontinuous first derivative at y = 0. Here we regularize the Charney problem by replacing a linear mean current, U = |y|, by either U = є log(cosh(y/є)) or U = є y erf(y/є), followed by matched asymptotic perturbation expansions in powers of the small regularization parameter є. The series is carried to third order because the lowest nonzero correction to the phase speed is O(є2) and this correction is determined simultaneously with the third order approximation to the eigenfunction. The result is both an explicit, analytic regularization of a problem important in atmospheric and ocean dynamics, but also a good school problem because the series is explicit with nothing worse than polylogarithms and confluent hypergeometric functions. The primary meteorological conclusion is that the delta functions in the Charney problem are harmless as demonstrated both by third order perturbation theory and by spectrally-accurate numerical solutions. The physics of the regularized Charney problem is not significantly changed from that of the original Charney problem.

Oceanic Eddy Characteristics and Generation Mechanisms in the Kuroshio Extension Region

AbstractThe Kuroshio Extension region is well known for its strong eddy activity. In this paper, using satellite altimetry‐measured sea surface height anomaly data from 1993 to 2012 in an extended Kuroshio Extension region (140–180°E, 25–45°N), we analyze eddy characteristics: eddy size, polarity, lifetime, intensity, trajectory, and spatial and temporal distributions. Using temperature and salinity vertical profiles measured by Argo floats, we examine the eddy impact on vertical stratification. During the 20‐year period, 7,574 eddies are identified (based on following complete eddy trajectories) with a lifetime equal to or longer than 4 weeks. The numbers of cyclonic and anticyclonic eddies are found to be approximately the same. The distribution of eddy sizes peaks at a radius of about 40 km. The radius at the peak is at the same order as the first baroclinic deformation radius or the horizontal shear scale of the Kuroshio flow. The normalized eddy statistical characteristics show that eddies have different characteristics at different stages of their lifetimes. Among eddies with lifetimes longer than 50 weeks, more anticyclonic (cyclonic) eddies are found north (south) of 35°N. In contrast, among eddies with lifetimes shorter than 20 weeks, more cyclonic (anticyclonic) eddies are found north (south) of 35°N. The asymmetric distribution of eddies suggests two different eddy generation mechanisms: (1) the development of meanders in the Kuroshio path leading to the pinch off of eddies with longer lifetime (larger size) and (2) horizontal shear instability (barotropic instability) leading to eddies of shorter life (smaller size). We further apply an eddy‐resolved numerical product to quantitatively investigate the eddy generation mechanisms.

Mesoscale variability of the Black Sea circulation by the simulation results in 2011 and 2016

The paper analyzed the simulation results of the Black Sea circulation in 2011 and 2016. Numerical experiments were carried out with a horizontal resolution of 1.6 km and took into account the real atmospheric forcing SKIRON. The simulated temperature and salinity were compared with observations. It revealed the correspondence to real data except some deviations caused by the movement of intensive mesoscale eddies. Analysis of velocity fields showed that the reducing of the wind influence in 2016 led to the weakening of the Rim Current in contrast to 2011. The Sevastopol anticyclone in the spring-summer of 2011 was weaker than in 2016, but the Batumi anticyclone was stronger in 2011 due to intensification of coastal circulation near the Anatolian coast. The eddy energy associated with mesoscale variability had been growing in 2011 due to both barotropic and baroclinic instability processes. The eddy energy had been decreasing during 2016 and started increasing only at the end of 2016 due to development of baroclinic instability processes.

The Influence of Topography on the Stability of the Norwegian Atlantic Current off Northern Norway

AbstractThe steep continental slope off the Lofoten–Vesterålen islands of northern Norway appears to be the source of the most intense mesoscale eddy field in all of the Nordic Seas. Here we use linearized two-layer shallow-water equations to study the stability of the Norwegian Atlantic Current in this region. The study extends previous works that used linearized quasigeostrophic vertical mode equations to examine the effects of bottom topography on baroclinic instability here. We find evidence of baroclinic instability in the stacked shallow-water model but also of barotropic instability that is associated with the upper-layer lateral shear. The calculations give an indication that growth rates of barotropic instability may be comparable to or larger than those of baroclinic instability over the steepest parts of the continental slope.

Topographic Influence on Baroclinic Instability and the Mesoscale Eddy Field in the Northern North Atlantic Ocean and the Nordic Seas

AbstractA weak planetary vorticity gradient and weak density stratification in the northern North Atlantic and Nordic seas lead to time-mean currents that are strongly guided by bottom topography. The topographic steering sets up distinct boundary currents with strong property fronts that are prone to both baroclinic and barotropic instability. These instability processes generate a macroturbulent eddy field that spreads buoyancy and other tracers out from the boundary currents and into the deep basins. In this paper we investigate the particular role played by baroclinic instability in generating the observed eddy field, comparing predictions from linear stability calculations with diagnostics from a nonlinear eddy-permitting ocean model hindcast. We also look into how the bottom topography impacts instability itself. The calculations suggest that baroclinic instability is a consistent source of the eddy field but that topographic potential vorticity gradients impact unstable growth significantly. We also observe systematic topographic effects on finite-amplitude eddy characteristics, including a general suppression of length scales over the continental slopes. Investigation of the vertical structure of unstable modes reveal that Eady theory, even when modified to account for a bottom slope, is unfit as a lowest-order model for the dynamics taking place in these ocean regions.

The Momentum Budget in the Stratosphere, Mesosphere, and Lower Thermosphere. Part II: The In Situ Generation of Gravity Waves

The contributions of gravity waves to the momentum budget in the mesosphere and lower thermosphere (MLT) is examined using simulation data from the Ground-to-Topside Model of Atmosphere and Ionosphere for Aeronomy (GAIA) whole-atmosphere model. Regardless of the relatively coarse model resolution, gravity waves appear in the MLT region. The resolved gravity waves largely contribute to the MLT momentum budget. A pair of positive and negative Eliassen–Palm flux divergences of the resolved gravity waves are observed in the summer MLT region, suggesting that the resolved gravity waves are likely in situ generated in the MLT region. In the summer MLT region, the mean zonal winds have a strong vertical shear that is likely formed by parameterized gravity wave forcing. The Richardson number sometimes becomes less than a quarter in the strong-shear region, suggesting that the resolved gravity waves are generated by shear instability. In addition, shear instability occurs in the low (middle) latitudes of the summer (winter) MLT region and is associated with diurnal (semidiurnal) migrating tides. Resolved gravity waves are also radiated from these regions. In Part I of this paper, it was shown that Rossby waves in the MLT region are also radiated by the barotropic and/or baroclinic instability formed by parameterized gravity wave forcing. These results strongly suggest that the forcing by gravity waves originating from the lower atmosphere causes the barotropic/baroclinic and shear instabilities in the mesosphere that, respectively, generate Rossby and gravity waves and suggest that the in situ generation and dissipation of these waves play important roles in the momentum budget of the MLT region.

Instability-Driven Benthic Storms below the Separated Gulf Stream and the North Atlantic Current in a High-Resolution Ocean Model

AbstractBenthic storms are important for both the energy budget of the ocean and for sediment resuspension and transport. Using 30 years of output from a high-resolution model of the North Atlantic, it is found that most of the benthic storms in the model occur near the western boundary in association with the Gulf Stream and the North Atlantic Current, in regions that are generally collocated with the peak near-bottom eddy kinetic energy. A common feature is meander troughs in the near-surface jets that are accompanied by deep low pressure anomalies spinning up deep cyclones with near-bottom velocities of up to more than 0.5 m s−1. A case study of one of these events shows the importance of both baroclinic and barotropic instability of the jet, with energy being extracted from the jet in the upstream part of the meander trough and partly returned to the jet in the downstream part of the meander trough. This motivates examining the 30-yr time mean of the energy transfer from the (annual mean) background flow into the eddy kinetic energy. This quantity is shown to be collocated well with the region in which benthic storms and large increases in deep cyclonic relative vorticity occur most frequently, suggesting an important role for mixed barotropic–baroclinic instability-driven cyclogenesis in generating benthic storms throughout the model simulation. Regions of the largest energy transfer and most frequent benthic storms are found to be the Gulf Stream west of the New England Seamounts and the North Atlantic Current near Flemish Cap.

Optimal Initial Error Growth in the Prediction of the Kuroshio Large Meander Based on a High-resolution Regional Ocean Model

Based on the high-resolution Regional Ocean Modeling System (ROMS) and the conditional nonlinear optimal perturbation (CNOP) method, this study explored the effects of optimal initial errors on the prediction of the Kuroshio large meander (LM) path, and the growth mechanism of optimal initial errors was revealed. For each LM event, two types of initial error (denoted as CNOP1 and CNOP2) were obtained. Their large amplitudes were found located mainly in the upper 2500 m in the upstream region of the LM, i.e., southeast of Kyushu. Furthermore, we analyzed the patterns and nonlinear evolution of the two types of CNOP. We found CNOP1 tends to strengthen the LM path through southwestward extension. Conversely, CNOP2 has almost the opposite pattern to CNOP1, and it tends to weaken the LM path through northeastward contraction. The growth mechanism of optimal initial errors was clarified through eddy-energetics analysis. The results indicated that energy from the background field is transferred to the error field because of barotropic and baroclinic instabilities. Thus, it is inferred that both barotropic and baroclinic processes play important roles in the growth of CNOP-type optimal initial errors.

Decadal Variability of Eddy Characteristics and Energetics in the Kuroshio Extension: Unstable Versus Stable States

AbstractUsing the Estimating Circulation and Climate of the Ocean Phase II product, this study investigates the eddy characteristics and energetics in two dynamical regimes within the Kuroshio Extension (KE) region. Based on the empirical orthogonal function analysis, it is found that the decadal evolution of eddy kinetic energy in the KE region is characterized by a delayed negative correlation between the upstream and downstream. Besides the out‐of‐phase change in eddy activity levels, eddy characteristics and energetics in the upstream KE are also different from those in the downstream. In the upstream region, eddies are stretched in the zonal direction and the unstable state is characterized by strong eddy‐shedding processes. During these processes, cyclonic eddies are found to be formed in the first quasi‐stationary meander and finally pinched off from the KE jet. Energy budget illustrates that the eddy shedding processes are triggered by baroclinic instability, while barotropic instability becomes the dominated energy source after the meander is fully developed. Accompanied by the generation of cyclonic eddies, significant inverse energy cascade is detected. When the upstream KE is stable, the eddy activity is dominated by baroclinic instability and the inverse energy cascade occurs similarly. Distinct from the upstream, eddies in the downstream KE tend to be stretched in the meridional direction and the decadal variability of eddy kinetic energy in this region is mainly regulated by the upstream through lateral advection.

Interannual Eddy Kinetic Energy Modulations in the Agulhas Return Current

AbstractInterannual variability in the mesoscale eddy field in the Agulhas Return Current (ARC) of 32–42°S and 15–35°E is investigated based on satellite altimeter observations and state estimate from the Estimating the Circulation and Climate of the Ocean, Phase II from 1993 to 2016. It is found that the interannual modulation of eddy kinetic energy in the ARC region is externally mediated by the wind stress forcing that generates the westward propagating sea surface height anomalies across the South Indian Ocean subtropical gyre. The wind‐forced sea surface height anomalies influence the upstream Agulhas Current volume transports. By modulating the intensity of barotropic instability of the ARC mean flow centered around the retroflection region, the Agulhas Current inflow variability leads to the downstream interannual eddy kinetic energy fluctuations in the ARC region.

Numerical Studies of Submesoscale Island Wakes in the Kuroshio

AbstractSubmesoscale 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.