A numerical study of island wakes in the Southern California Bight

With the existence of eight substantial islands in the Southern California Bight, the oceanic circulation is significantly affected by island wakes. In this paper a high-resolution numerical model (on a 1 km grid), forced by a high-resolution wind (2 km), is used to study the wakes. Island wakes arise due both to currents moving past islands and to wind wakes that force lee currents in response. A comparison between simulations with and without islands shows the surface enstrophy (i.e., area-integrated square of the vertical component of vorticity at the surface) decreases substantially when the islands in the oceanic model are removed, and the enstrophy decrease mainly takes place in the areas around the islands. Three cases of wake formation and evolution are analyzed for the Channel Islands, San Nicolas Island, and Santa Catalina Island. When flows squeeze through gaps between the Channel Islands, current shears arise, and the bottom drag makes a significant contribution to the vorticity generation. Downstream the vorticity rolls up into submesoscale eddies. When the California Current passes San Nicolas Island from the northwest, a relatively strong flow forms over the shelf break on the northeastern coast and gives rise to a locally large bottom stress that generates anticyclonic vorticity, while on the southwestern side, with an adverse flow pushing the main wake current away from the island, positive vorticity has been generated and a cyclonic eddy detaches into the wake. When the northward Southern California Countercurrent passes the irregular shape of Santa Catalina Island, cyclonic eddies form on the southeastern coast of the island, due primarily to lateral stress rather than bottom stress; they remain coherent as they detach and propagate downstream, and thus they are plausible candidates for the submesoscale “spirals on the sea” seen in many satellite images. Finally, the oceanic response to wind wakes is analyzed in a spin-up experiment with a time-invariant wind that exhibits strips of both positive and negative curl in the island lee. Corresponding vorticity strips in the ocean develop through the mechanism of Ekman pumping.     

Tidally generated island wakes and surface water cooling over Izu Ridge

The Kuroshio current is well known for generating cold wakes behind islands over Izu Ridge in Northwestern Pacific. Observational data from the geostationary Himawari-8 satellite for 2015–2017 revealed the occurrence of cold waters during the period when the Kuroshio current flows away from the islands. With a focus on tidal currents, this study presents an investigation of dynamical processes responsible for the formation of areas with low sea surface temperature (SST) through the adoption of a high-resolution numerical ocean model for an event that happened in July 2017. Areas with cold water emerged only when tidal currents are included in the numerical model. The model results indicate the cold surface waters are formed in the vicinity of the islands because of upwelling and vertical mixing. Qualitative features of the cold water formation for each island are found to depend on its size, topography, and ambient currents. Near Kozu Island, the tidal excursion is large enough to cause eddy shedding. These shed eddies are stirred by tidal currents to extend the surface cooling effect to wider areas. Near Hachijo Island, a persistent wake is formed by the ambient northward current. Inclusion of tidal currents destabilizes the wake, and consequently leads to the formation of a low SST area, although no clear eddy shedding is detected. The flow patterns around the islands are classified using an additional non-dimensional parameter, defined as the ratio between tidal excursion and island diameter.

     

Regional turbulence patterns driven by meso- and submesoscale processes in the Caribbean Sea

The surface ocean circulation in the Caribbean Sea is characterized by the interaction between anticyclonic eddies and the Caribbean Upwelling System (CUS). These interactions lead to instabilities that modulate the transfer of kinetic energy up- or down-cascade. The interaction of North Brazil Current rings with the islands leads to the formation of submesoscale vorticity filaments leeward of the Lesser Antilles, thus transferring kinetic energy from large to small scales. Within the Caribbean, the upper ocean dynamic ranges from large-scale currents to coastal upwelling filaments and allow the vertical exchange of physical properties and supply KE to larger scales. In this study, we use a regional model with different spatial resolutions (6, 3, and 1 km), focusing on the Guajira Peninsula and the Lesser Antilles in the Caribbean Sea, in order to evaluate the impact of submesoscale processes on the regional KE energy cascade. Ageostrophic velocities emerge as the Rossby number becomes O(1). As model resolution is increased submesoscale motions are more energetic, as seen by the flatter KE spectra when compared to the lower resolution run. KE injection at the large scales is greater in the Guajira region than in the others regions, being more effectively transferred to smaller scales, thus showing that submesoscale dynamics is key in modulating eddy kinetic energy and the energy cascade within the Caribbean Sea.     

Storm surge and tidal interaction in the Tjeldsund channel, northern Norway

The aim of this study is to investigate tide–surge interaction in narrow channels with complex and relatively shallow topography. A high-resolution depth-integrated tidal and storm surge model has been implemented for the Tjeldsund channel which is an important sailing lane in northern Norway. A horizontal grid resolution down to 50 m is applied in order to represent the complex bottom topography and the formation of jets and small-scale eddies. Two typically storm surge events in December 2004 have been examined in detail. The tide–surge interaction is found to influence the generation of higher harmonics and the formation of eddies in the current field. In some cases, the magnitude of storm surge currents may reach the same magnitude as the tidal currents enhancing the formation of jets and eddies.     

Residual vorticity induced by the action of tidal currents in combination with bottom topography in the Southern Bight of the North Sea

Abstract

Results are presented of calculations on the generation of residual vorticity by tidal currents over the bottom topography of the Southern Bight of the North Sea. A typical order of magnitude is 10^?6 to 10^?7 s ^?1. This is compared with current measurements on calm days, when similar magnitudes are found. At windspeeds less than about 5 m/s tidal generation of residual vorticity is important; at higher windspeeds wind effects begin to dominate. Our results are relevant in understanding the spatial variability of residual currents, because a non-zero vorticity implies the existence of horizontal gradients in the residual current field.     

Nonlinear interaction of the Tsugaru Warm Current and tide in the Tsugaru Strait

The Tsugaru Strait, which connects the Sea of Japan with the Pacific Ocean, is characterized by the eastward Tsugaru Warm Current (TWC) and oscillating tidal currents of similar magnitude. A 15-day current observation was conducted in one of the two narrow channels in the strait, at the northwest tip of the Shimokita Peninsula. The observation revealed that the spectral energy of the semidiurnal current exceeds that of the diurnal current, contrary to the conventional view. The Tsugaru Strait regional model was developed to study the mechanism of this spectral energy reversal (140–141.5° E, 40.4–42.6° N, 500?m grid resolution). At the eastern and western open boundaries, the model was driven by the constant Tsugaru warm current and tidal elevation, which was adjusted by comparing the model with tidal gauge observations within the channel. The relative magnitude of the spectral energies differed from that of the observation when the model was driven by tide only. However, the spectral energy levels were reversed when the model was driven by both tide and current. The nonlinear interaction of periodic tidal currents and the steady TWC was explained by the vorticity equation, which describes the production and advection of residual currents from tidal currents. According to the model results, flow separation and advection of vorticity by the TWC was the most prominent factor in this phenomenon. Because of the strong nonlinearities, flow separation around the headland occurred during the tidal period with dominant current magnitude and furnished the main difference between the diurnal and semidiurnal interactions. These phenomena were enhanced by the complex topography, and demonstrate the importance of scale interaction, especially when developing high-resolution regional models.     

The influence of tidal currents on the asymmetry of tide-dominated ebb–tidal deltas

An idealized numerical model is developed to study the spatial asymmetry of ebb–tidal deltas under influence of large-scale alongshore tidal currents. It is shown that the asymmetry of the delta depends on the magnitude of the cross-shore and large scale alongshore tidal currents, their phase difference, and on the width of the inlet. Model results are compared with observations of ebb–tidal deltas of the tidal inlet systems of the Dutch Wadden Sea and with the ebb–tidal delta of the Eastern Scheldt, located in the southwestern part of the Netherlands. The modeled current and residual sediment transport patterns agree well with observed ones. The modeled asymmetry of the ebb–tidal delta also agree with observed ones. Furthermore, bottom patterns are consistent with those found with a previous version of the idealized model which focused on the modeling of symmetric ebb–tidal deltas. However, the model is not able to reproduce the observed ebb-dominated channel. The underlying physical processes are explained in terms of vorticity dynamics. The convergence of the mean vorticity flux generates mean vorticity and thereby residual circulation. An analysis shows there is competition between two contributions to the convergence of the mean vorticity flux. This competition explains the sensitivity of the results to the model parameters.     

Observations of submesoscale eddies using high-frequency radar-derived kinematic and dynamic quantities

The spatio-temporal variability of submesoscale eddies off southern San Diego is investigated with two-year observations of subinertial surface currents [O(1) m depth] derived from shore-based high-frequency radars. The kinematic and dynamic quantities — velocity potential, stream function, divergence, vorticity, and deformation rates — are directly estimated from radial velocity maps using optimal interpolation. For eddy detection, the winding-angle approach based on flow geometry is applied to the calculated stream function. A cluster of nearly enclosed streamlines with persistent vorticity in time is identified as an eddy. About 700 eddies were detected for each rotation (clockwise and counter-clockwise). The two rotations show similar statistics with diameters in the range of 5–25 km and Rossby number of 0.2–2. They persist for 1–7 days with weak seasonality and migrate with a translation speed of 4–15 cm s⁻¹ advected by background currents. The horizontal structure of eddies exhibits nearly symmetric tangential velocity with a maximum at the defined radius of the eddy, non-zero radial velocity due to background flows, and Gaussian vorticity with the highest value at the center. In contrast divergence has no consistent spatial shape. Two episodic events are presented with other in situ data (subsurface current and temperature profiles, and local winds) as an example of frontal-scale secondary circulation associated with drifting submesoscale eddies.     

Propagation of barotropic modons over topography

Abstract

This is a broad survey of the interaction of modons with topography in a one-layer, quasigeostrophic model. Numerical simulations of modons interacting with ridges, hills, random topography and other obstacles are presented. The behavior of the modon is compared to numerical simulations of a two-point-vortex model, which proves a useful guide to the basic trajectory deflection mechanism. Under sufficiently strong but quasigeostrophically valid topographic perturbations, the modon is shown to fission into two essentially independent, oppositely-signed vortices. In the breakup of a modon near a hill it is found that the positive vortex migrates to the top of the hill. The resulting correlation between the positive vorticity trapped over the hill and the topography is in sharp contrast with the theories of turbulent flow over topography and generation of flow over topography by large scale forcing, both of which describe the development of vorticity anticorrelated with topography. A heuristic explanation of this new behavior is provided in terms of the dynamics of β bT-plane vortices. Further, it is found that a modon travelling over rough topography homogenizes the field of potential vorticity in its vicinity. This behavior is explained in terms of the induced eddy activity near the modon.     

The impacts of the large-scale hydraulic structures on tidal dynamics in open-type bay: numerical study in Jakarta Bay

The construction of a Giant Sea Wall (GSW) complex in Jakarta Bay has been proposed to protect Jakarta against flood in the Master Plan for National Capital Integrated Coastal Development (NCICD). However, these large-scale hydraulic structures could significantly change the tidal dynamics in Jakarta Bay. This research investigates the potential impacts of a GSW on the tidal dynamics, including tides, currents, and residual currents in Jakarta Bay using a validated numerical model (Finite Volume Coastal Ocean Model (FVCOM)). Results show that the bay is diurnal with a maximum tidal range of ~0.9 m. The flow is mainly in an east-west direction with a maximum depth-mean current speed of up to 0.3 ms^?1. The construction of a GSW would modulate the tidal dynamics by changing the bathymetry, tidal prism, wind effect, and tidal choking effect in the bay. The maximum tidal range would be slightly increased due to the reduced tidal prism of the bay and the increased tidal choking effect. The current would penetrate into the west reservoir through the gates and channels between the artificial islands, with peak speed jets appearing at the gates (~0.3 ms^?1), due to tidal choking. A similar peak current speed appears near the right wing of the GSW due to the pressure gradient would be created by the wing of the GSW. Closing the gates would mainly affect the currents inside the west reservoir. The residual current would be slightly increased after the construction of the GSW. An eddy would be formed at the bottom level near the right wing of the GSW. The direction of the residual current is landward instead of seaward at the surface level outside the GSW. The impact of wind on surface currents would be much reduced due to the decreased water surface area. Although this study is site specific, the findings may have a wider applicability to the impacts of large-scale hydraulic structures on tidal dynamics in open-type bays.     

Submesoscale features and their interaction with fronts and internal tides in a high-resolution coupled atmosphere-ocean-wave model of the Bay of Bengal

Large freshwater fluxes into the Bay of Bengal by rainfall and river discharges result in strong salinity fronts in the bay. In this study, a high-resolution coupled atmosphere-ocean-wave model with comprehensive physics is used to model the weather, ocean circulation, and wave field in the Bay of Bengal. Our objective is to explore the submesoscale activity that occurs in a realistic coupled model that resolves mesoscales and allows part of the submesoscale field. Horizontal resolution in the atmosphere varies from 2 to 6 km and is 13 km for surface waves, while the ocean model is submesoscale permitting with resolutions as high as 1.5 km and a vertical resolution of 0.5 m in the upper 10 m. In this paper, three different cases of oceanic submesoscale features are discussed. In the first case, heavy rainfall and intense downdrafts produced by atmospheric convection are found to force submesoscale currents, temperature, and salinity anomalies in the oceanic mixed layer and impact the mesoscale flow. In a second case, strong solitary-like waves are generated by semidiurnal tides in the Andaman Sea and interact with mesoscale flows and fronts and affect submesoscale features generated along fronts. A third source of submesoscale variability is found further north in the Bay of Bengal where river outflows help maintain strong salinity gradients throughout the year. For that case, a comparison with satellite observations of sea surface height anomalies, sea surface temperature, and chlorophyll shows that the model captures the observed mesoscale eddy features of the flow field, but in addition, submesoscale upwelling and downwelling patterns associated with ageostrophic secondary circulations along density fronts are also captured by the model.     

Tidal and residual currents over abrupt deep-sea topography based on shipboard ADCP data and tidal model solutions for three popular bathymetry grids

The response of tidal and residual currents to small-scale morphological differences over abrupt deep-sea topography (Seine Seamount) was estimated for bathymetry grids of different spatial resolution. Local barotropic tidal model solutions were obtained for three popular and publicly available bathymetry grids (Smith and Sandwell TOPO8.2, ETOPO1, and GEBCO08) to calculate residual currents from vessel-mounted acoustic Doppler current profiler (VM-ADCP) measurements. Currents from each tidal solution were interpolated to match the VM-ADCP ensemble times and locations. Root mean square (RMS) differences of tidal and residual current speeds largely follow topographic deviations and were largest for TOPO8.2-based solutions (up to 2.8 cm?s^?1) in seamount areas shallower than 1,000 m. Maximum RMS differences of currents obtained from higher resolution bathymetry did not exceed 1.7 cm?s^?1. Single depth-dependent maximum residual flow speed differences were up to 8 cm?s^?1 in all cases. Seine Seamount is located within a strong mean flow environment, and RMS residual current speed differences varied between 5 % and 20 % of observed peak velocities of the ambient flow. Residual flow estimates from shipboard ADCP data might be even more sensitive to the choice of bathymetry grids if barotropic tidal models are used to remove tides over deep oceanic topographic features where the mean flow is weak compared to the magnitude of barotropic tidal, or baroclinic currents. Realistic topography and associated flow complexity are also important factors for understanding sedimentary and ecological processes driven and maintained by flow–topography interaction.     

Modeling the formation of undulations of the coastline: The role of tides

An idealized model is developed and analyzed to investigate the relevance of tidal motion for the emergence of undulations of a sandy coastline. The model describes feedbacks between tidal and steady flow on the inner shelf, sand transport in the nearshore zone and an irregular coastline. It is demonstrated that an initially straight coastline can become unstable with respect to perturbations with a rhythmic structure in the alongshore direction. The mechanism causing the growth of these perturbations is explained in terms of vorticity concepts. The relative importance of tide-related and wave-driven sediment fluxes in generating undulations of the coastline is investigated for the Dutch coast. Using parameter values that are appropriate for the Dutch coast it is found that tides can render a straight coastline unstable. The model predicts a fastest growing mode (FGM) with a wavelength that is in the order of the observed length of barrier islands. The mode grows on a time scale of 50 yr and it migrates 200 m per year. The wavelength of the FGM decreases with increasing amplitude of the tidal currents. This result is consistent with data of tides, waves and the lengths of barrier islands that are located along the Dutch and German Wadden coast.     

On the exact shape of the horizontal profile of a topographically rectified tidal flow

Abstract

Starting from the nonlinear shallow water equations of a homogeneous rotating fluid we derive the equation describing the evolution of vorticity by a fluctuating bottom topography of small amplitude, using a multiple scale expansion in a small parameter, which is the topographic length scale relative to the tidal wave length. The exact response functions of residual vorticity for a sinusoidal bottom topography are compared with those obtained by a primitive perturbation series and by harmonic truncation, showing the former to be invalid for small topographic length scales and the latter to be only a fair approximation for vorticity produced by planetary vortex stretching. In deriving the exact shape of the horizontal residual velocity profile at a step-like break in the bottom topography, it is shown that the Lagrangian profile only exists in a strip having the width of the amplitude of the tidal excursion at both sides of the break, and that it vanishes outside that interval. Moreover, in the limit of small amplitude topography at least, it vanishes altogether for the generation mechanism by means of planetary vortex stretching. The Eulerian profile is shown to extend over twice the interval of the Lagrangian profile both for production by vortex stretching and by differential bottom friction. These finite intervals over which the residual velocity profiles exist for a step-like topography are not reproduced by harmonic truncation of the basic equation. This method gives exponentially decaying profiles, indicating spurious horizontal diffusion of vorticity. In terms of orders of magnitude, the method of harmonic truncation is reliable for residual velocity produced by vortex stretching but it overestimates the residual velocity produced by differential bottom friction by a factor 2.     

A study of suspended sediment concentration in Yangshan deep-water port in Shanghai,China

Suspended sediment concentration (SSC) plays an important role in the estuarine environment.Its spatial or temporal variations in coastal zones and estuaries indicate that sediments are suspended,trans...     

Tidal and wind-driven surface currents in the German Bight: HFR observations versus model simulations

Tidal and wind-driven surface currents in the German Bight between shallow mudflats of the North Frisian islands and the island of Helgoland are studied using coastal high-frequency radar (HFR) observations and hindcasts from a primitive equation numerical model. The setup of the observational system is described, and estimates of expected measurement errors are given. A quantitative comparison of numerical model results and observations is performed. The dominant tidal components are extracted from the two data sources using tidal harmonic analysis and the corresponding tidal ellipses are defined. Results show that the spatial patterns of different tidal ellipse parameters are consistent in the two data sets. Model sensitivity studies with constant and variable salinity and temperature distributions are used to study density-related mechanisms of circulation. Furthermore, the role of the surface wind field in driving the German Bight circulation is investigated using the complex correlation between wind and surface current vectors. The observed change of the respective correlation patterns from the coastal to open ocean is shown to be due to a combination of density effects, the coastline and topography. The overall conclusion is that HFR observations resolve the small-scale and rapidly evolving characteristics of coastal currents well in the studied area and could present an important component for regional operational oceanography when combined with numerical modelling. Some unresolved issues associated with the complex circulation and large instability of circulation in front of the Elbe River Estuary justify further considerations of this area using dedicated surveys and modelling efforts.     

Contemporary hydrodynamics and morphological change of a microtidal estuary: a numerical modelling study

Contemporary hydrodynamics and morphological change are examined in a shallow microtidal estuary, located on a wave-dominated coast (Port Stephens, NSW, Australia). Process-based numerical modelling is undertaken by combining modules for hydrodynamics, waves, sediment transport and bathymetry updates. Model results suggest that the complex estuarine bathymetry and geometry give rise to spatial variations in the tidal currents and a marked asymmetry between ebb and flood flows. Sediment transport paths correspond with tidal asymmetry patterns. The SE storms significantly enhance the quantities of sediment transport, while locally generated waves by the westerly strong winds also are capable of causing sediment entrainment and contribute to the delta morphological change. The wave/wind-induced currents are not uniform with flow over shoals driven in the same direction as waves/winds while a reverse flow occurring in the adjacent channel. The conceptual sediment transport model developed in this study shows flood-directed transport occurs on the flood ramp while ebb-directed net transport occurs in the tidal channels and at the estuary entrance. Accretion of the intertidal sand shoals and deepening of tidal channels, as revealed by the model, suggest that sediment-infilling becomes advanced, which may lead to an ebb-dominated estuary. It is likely that a switch from flood- to ebb-dominance occurs during the estuary evolution, and the present-day estuary acts as a sediment source rather than sediment sink to the coastal system. This is conflictive to the expectation drawn from the estuarine morphology; however, it is consistent with previous research suggesting that, in an infilling estuary, an increase in build-up of intertidal flats/shoals can eventually shift an estuary towards ebb dominance. Thus, field data are needed to validate the result presented here, and further study is required to investigate a variety of estuaries in the Australian area.     

A numerical study of island wake generated by an elliptical tidal flow

An idealized numerical study of the influence of a tidal flow around an island has been undertaken with ROMS. The study focusses on coastal island wakes which are mainly controlled by elliptical tidal current flows on shallow shelves. This model is typical of some isolated continental shelf islands. The model is forced by a semi-diurnal barotropic inertia gravity wave imposed on the four open boundaries of a rectangular domain and its propagation results in an elliptical tidal flow within the domain in which the circular island lies. The influence of the surrounding island bathymetry and of the ellipse shape has been studied both in two and three dimensions. In the island vicinity, the residual circulation patterns over a tidal period show alongshore flow divergence along the major axis and convergence along the minor axis. A thin tidal ellipse (i.e. with a large ratio between major and minor axes) leads to strong eddy activity periods in the lee of the island during the flood and ebb phases, with eddy dissipation phases in between. By contrast, an almost round ellipse (axis ratio nearly 1) leads to vorticity filaments which continuously progress around the island without eddy shedding. The presence of a topographic slope in the vicinity of the island strengthens the eddy activity. This study suggests that the tidal current rotation favors the development of the eddy rotating in the same direction and weakens the development of the second eddy. In three dimensions with a surrounding bathymetry, an intense upwelling occurs in a large area in the lee of the island and the vertical velocities are stronger with thinner ellipses. With a flat bottom the vertical motions are almost fully generated by convergence and divergence of the secondary flow. With a varying bottom topography, the vertical motions come from a combination of this mechanism with convergence and divergence of the depth averaged flow.     

Baroclinic and topographic influences on the transport in western boundary currents

Abstract

One of the central unsolved theoretical problems of the large scale ocean circulation is concerned with explaining the very large transports measured in western boundary currents such as the Gulf Stream and the Kuroshio. The only theory up to now that can explain the size of these transports is that of non-linear recirculation in which the advective terms in the momentum equations became important near the western boundary. In this paper an alternative explanation is suggested. When bottom topography and baroclinic effects are included in a wind-driven ocean model it is shown that the western boundary current can have a transport larger than that predicted from the wind stress distribution even when the nonlinear advective terms are ignored. The explanation lies in the presence of pressure torques associated with bottom topography which can contribute to the vorticity balance in the same sense as the wind stress curl.

Three numerical experiments have been carried out to explore the nature of this process using a three dimensional numerical model. The first calculation is done for a baroclinic ocean of constant depth, the second for a homogeneous ocean with an idealized continental slope topography, and the third for a baroclinic ocean with the same continental slope topography. The nature of the vorticity balance and of the circulation around closed paths is examined in each case, and it is shown that bottom pressure torques lead to enhanced transport in the western boundary current only for the baroclinic case with variable depth.     

Wind-induced subduction at the South Atlantic subtropical front

The South Atlantic Subtropical Front, associated with the eastward-flowing South Atlantic Current, separates the colder, nutrient-rich waters of the subpolar gyre from the warmer, nutrient-poor waters of the subtropical gyre. Perturbations to the quasi-geostrophic, eastward flow generate meanders and filaments which induce cross-frontal exchange of water properties. Down-front winds transport denser waters from the South over warm waters from the North, inducing convective instability and subduction. Such processes occur over spatial scales of the order of 1 km and thus require high horizontal spatial resolution. In this modeling study, a high-resolution (4 km) regional grid is embedded in a basin-wide configuration (12 km) of the South Atlantic Ocean in order to test the importance of submesoscale processes in water mass subduction along the subtropical front. Stronger and more numerous eddies obtained in the high-resolution run yield more intense zonal jets along the frontal zone. Such stronger jets are more susceptible to instabilities, frontogenesis, and the generation of submesoscale meanders and filaments with \(\mathcal {O}(1)\) Rossby number. As a consequence, vertical velocities larger than 100 md ¹ are obtained in the high-resolution run, one order of magnitude larger than in the low-resolution run. Wind-driven subduction occurs along the frontal region, associated with negative Ertel potential vorticity in the surface layer. Such processes are not observed in the low-resolution run. A passive tracer experiment shows that waters with density characteristics similar to subtropical mode waters are preferentially subducted along the frontal region. The wind-driven buoyancy flux is shown to be much larger than thermal or haline fluxes during the wintertime, which highlights the importance of the frictional component in extracting PV from the surface ocean and inducing subduction, a process that has been overlooked in subtropical mode water formation in the region.