Browse
Abstract
The Beaufort Gyre (BG) is hypothesized to be partially equilibrated by those mesoscale eddies that form via baroclinic instabilities of its currents. However, our understanding of the eddy field’s dependence on the mean BG currents and the role of sea ice remains incomplete. This theoretical study explores the scales and vertical structures of eddies forming specifically due to baroclinic instabilities of interior BG flows. An idealized quasigeostrophic model is used to show that flows driven only by the Ekman pumping contain no interior potential vorticity (PV) gradients and generate weak and large eddies,
Abstract
The Beaufort Gyre (BG) is hypothesized to be partially equilibrated by those mesoscale eddies that form via baroclinic instabilities of its currents. However, our understanding of the eddy field’s dependence on the mean BG currents and the role of sea ice remains incomplete. This theoretical study explores the scales and vertical structures of eddies forming specifically due to baroclinic instabilities of interior BG flows. An idealized quasigeostrophic model is used to show that flows driven only by the Ekman pumping contain no interior potential vorticity (PV) gradients and generate weak and large eddies,
Abstract
Previous in situ observations have suggested that bottom water temperature variations in shelf seas can drive significant ocean bottom heat flux (BHF) by heat conduction. The BHF-driven bottom water temperature variations, however, have been overlooked in ocean general circulation models. In this study, we established a sea-sediment fully coupled model through incorporating the BHF. The coupled model included a sediment temperature module/model, and the BHF was calculated based on the sediment heat content variations. Meanwhile, we applied temporally varying BHF in the calculation of the bottom water temperature, which further determined the sediment temperature. The two-way coupled BHF process presents a more complete and physically reasonable heat budget in the ocean model and a synchronously varying sediment temperature profile. The coupled model was validated using a one-dimensional test case, and then it was applied in a domain covering the Bohai and Yellow Seas. The results suggest that when a strong thermocline exists, the BHF can change the bottom water temperature by more than 1°C because the effects of the BHF are limited to within a shallow bottom layer. However, when the water column is well mixed, the BHF changes the temperature of the entire water column, and the heat transported across the bottom boundary is ventilated to the atmosphere. Thus, the BHF has less effect on water temperature and may directly affect air–sea heat flux. The sea-sediment interactions dampen the amplitude of the bottom water temperature variations, which we propose calling the seabed dampening ocean heat content variation mechanism (SDH).
Abstract
Previous in situ observations have suggested that bottom water temperature variations in shelf seas can drive significant ocean bottom heat flux (BHF) by heat conduction. The BHF-driven bottom water temperature variations, however, have been overlooked in ocean general circulation models. In this study, we established a sea-sediment fully coupled model through incorporating the BHF. The coupled model included a sediment temperature module/model, and the BHF was calculated based on the sediment heat content variations. Meanwhile, we applied temporally varying BHF in the calculation of the bottom water temperature, which further determined the sediment temperature. The two-way coupled BHF process presents a more complete and physically reasonable heat budget in the ocean model and a synchronously varying sediment temperature profile. The coupled model was validated using a one-dimensional test case, and then it was applied in a domain covering the Bohai and Yellow Seas. The results suggest that when a strong thermocline exists, the BHF can change the bottom water temperature by more than 1°C because the effects of the BHF are limited to within a shallow bottom layer. However, when the water column is well mixed, the BHF changes the temperature of the entire water column, and the heat transported across the bottom boundary is ventilated to the atmosphere. Thus, the BHF has less effect on water temperature and may directly affect air–sea heat flux. The sea-sediment interactions dampen the amplitude of the bottom water temperature variations, which we propose calling the seabed dampening ocean heat content variation mechanism (SDH).
Abstract
Leading modes of interannual variability in upper-ocean salinity in the tropical Indian Ocean (TIO) and their connections were studied based on 17 years (2002–18) of oceanic historical and reanalysis data. Empirical orthogonal function (EOF) analysis depicted the dominant roles of the first two leading modes in salinity variability in the TIO over a wide range of interannual time scales. Among the rich oscillations of the leading EOF modes, a coherent near-biennial band was identified with basinwide loading of sea surface salinity anomalies (SSSa) (EOF1) leading/lagging the northeast–southwest dipolar mode of SSSa (EOF2) by around 4 months across the TIO, with southwestward migration of SSSa center. The spatial loadings of the SSSa leading modes in the TIO were strongly shaped by sea surface temperature–related freshwater fluxes and wind-driven regional ocean circulation on a near-biennial time scale. Composite analysis of the mixed layer salinity budget reflected characteristic features of basin-scale ocean–atmosphere coupling, both temporally and regionally during the life cycle of the near-biennial fluctuation in anomalous salinity in the TIO. Consistent with the intrinsic oscillation paradigm in the observed Indian Ocean dipole (IOD) variation, the dynamic and thermodynamic feedbacks associated with switches from the positive to negative IOD modes provided the phase-connection mechanisms for the SSSa leading-mode displacement over the TIO.
Significance Statement
This study investigates the leading modes of interannual variability in upper-ocean salinity in the tropical Indian Ocean (TIO). The intrinsic oscillation and associated dynamic and thermodynamic feedbacks over the TIO drive the basinwide connections of upper-ocean salinity variability. Our results show that a coherent near-biennial band is identifiable within the leading modes of sea surface salinity anomalies (SSSa), in which the wind-induced horizontal advections and evaporation-minus-precipitation anomalies associated with the switches from positive to negative Indian Ocean dipole modes mainly provide the phase-transition mechanism of SSSa. This research illustrates substantial evidence for the displacement of basin-scale sea surface temperature anomalies modulating the structures of SSSa and inducing the dynamical connections of leading modes of SSSa on the near-biennial time scale.
Abstract
Leading modes of interannual variability in upper-ocean salinity in the tropical Indian Ocean (TIO) and their connections were studied based on 17 years (2002–18) of oceanic historical and reanalysis data. Empirical orthogonal function (EOF) analysis depicted the dominant roles of the first two leading modes in salinity variability in the TIO over a wide range of interannual time scales. Among the rich oscillations of the leading EOF modes, a coherent near-biennial band was identified with basinwide loading of sea surface salinity anomalies (SSSa) (EOF1) leading/lagging the northeast–southwest dipolar mode of SSSa (EOF2) by around 4 months across the TIO, with southwestward migration of SSSa center. The spatial loadings of the SSSa leading modes in the TIO were strongly shaped by sea surface temperature–related freshwater fluxes and wind-driven regional ocean circulation on a near-biennial time scale. Composite analysis of the mixed layer salinity budget reflected characteristic features of basin-scale ocean–atmosphere coupling, both temporally and regionally during the life cycle of the near-biennial fluctuation in anomalous salinity in the TIO. Consistent with the intrinsic oscillation paradigm in the observed Indian Ocean dipole (IOD) variation, the dynamic and thermodynamic feedbacks associated with switches from the positive to negative IOD modes provided the phase-connection mechanisms for the SSSa leading-mode displacement over the TIO.
Significance Statement
This study investigates the leading modes of interannual variability in upper-ocean salinity in the tropical Indian Ocean (TIO). The intrinsic oscillation and associated dynamic and thermodynamic feedbacks over the TIO drive the basinwide connections of upper-ocean salinity variability. Our results show that a coherent near-biennial band is identifiable within the leading modes of sea surface salinity anomalies (SSSa), in which the wind-induced horizontal advections and evaporation-minus-precipitation anomalies associated with the switches from positive to negative Indian Ocean dipole modes mainly provide the phase-transition mechanism of SSSa. This research illustrates substantial evidence for the displacement of basin-scale sea surface temperature anomalies modulating the structures of SSSa and inducing the dynamical connections of leading modes of SSSa on the near-biennial time scale.
Abstract
Lateral mesoscale eddy-induced tracer transport is traditionally represented in coarse-resolution models by the flux–gradient relation. In its most complete form, the relation assumes the eddy tracer flux as a product of the large-scale tracer concentration gradient and an eddy transport coefficient tensor. However, several recent studies reported that the tensor has significant spatiotemporal complexity and is not uniquely defined, that is, it is sensitive to the tracer distributions and to the presence of nondivergent (“rotational”) components of the eddy flux. These issues could lead to significant biases in the representation of the eddy-induced transport. Using a high-resolution tracer model of the Gulf Stream region, we examine the diffusive and advective properties of lateral eddy-induced transport of dynamically passive tracers, reevaluate the utility of the flux–gradient relation, and propose an alternative approach based on modeling the local eddy forcing by a combination of diffusion and generalized eddy-induced advection. Mesoscale eddies are defined by a scale-based spatial filtering, which leads to the importance of new eddy-induced terms, including eddy-mean covariances in the eddy fluxes. The results show that the biases in representing these terms are noticeably reduced by the new approach. A series of targeted simulations in the high-resolution model further demonstrates that the approach outperforms the flux–gradient model in reproducing the stirring and dispersing effect of eddies. Our study indicates potential to upgrade the traditional flux–gradient relation for representing the eddy-induced tracer transport.
Abstract
Lateral mesoscale eddy-induced tracer transport is traditionally represented in coarse-resolution models by the flux–gradient relation. In its most complete form, the relation assumes the eddy tracer flux as a product of the large-scale tracer concentration gradient and an eddy transport coefficient tensor. However, several recent studies reported that the tensor has significant spatiotemporal complexity and is not uniquely defined, that is, it is sensitive to the tracer distributions and to the presence of nondivergent (“rotational”) components of the eddy flux. These issues could lead to significant biases in the representation of the eddy-induced transport. Using a high-resolution tracer model of the Gulf Stream region, we examine the diffusive and advective properties of lateral eddy-induced transport of dynamically passive tracers, reevaluate the utility of the flux–gradient relation, and propose an alternative approach based on modeling the local eddy forcing by a combination of diffusion and generalized eddy-induced advection. Mesoscale eddies are defined by a scale-based spatial filtering, which leads to the importance of new eddy-induced terms, including eddy-mean covariances in the eddy fluxes. The results show that the biases in representing these terms are noticeably reduced by the new approach. A series of targeted simulations in the high-resolution model further demonstrates that the approach outperforms the flux–gradient model in reproducing the stirring and dispersing effect of eddies. Our study indicates potential to upgrade the traditional flux–gradient relation for representing the eddy-induced tracer transport.
Abstract
We analyze approximately four years of heat-flux measurements at two levels, profiles of air temperature, and multiple measurements of the water temperature collected at a coastal zone site. Our analysis considers underestimation of the sea surface flux resulting from vertical divergence of the heat flux between the surface and the lowest flux level. We examine simple relationships of the heat flux to the wind speed and stratification and the potential influence of fetch and temperature advection. The fetch ranges from about 4 to near 400 km. For a given wind-direction sector, the transfer coefficient varies only slowly with increasing instability but decreases significantly with increasing stability. The intention here is not to recommend a new parameterization but rather to establish relationships that underlie the bulk formula that could lead to assessments of uncertainty and improvement of the bulk formula.
Significance Statement
The behavior of surface heat fluxes in the coastal zone is normally more complex than over the open ocean but has a large impact on human activity. Our study examines extensive flux measurements on a tower in the Baltic Sea that allows partitioning of the fluxes according to wind direction without seriously depleting the data for a given wind-direction sector. Because some of the normal assumptions for the usual parameterization are not met, our study examines relationships behind the parameterization of the surface fluxes.
Abstract
We analyze approximately four years of heat-flux measurements at two levels, profiles of air temperature, and multiple measurements of the water temperature collected at a coastal zone site. Our analysis considers underestimation of the sea surface flux resulting from vertical divergence of the heat flux between the surface and the lowest flux level. We examine simple relationships of the heat flux to the wind speed and stratification and the potential influence of fetch and temperature advection. The fetch ranges from about 4 to near 400 km. For a given wind-direction sector, the transfer coefficient varies only slowly with increasing instability but decreases significantly with increasing stability. The intention here is not to recommend a new parameterization but rather to establish relationships that underlie the bulk formula that could lead to assessments of uncertainty and improvement of the bulk formula.
Significance Statement
The behavior of surface heat fluxes in the coastal zone is normally more complex than over the open ocean but has a large impact on human activity. Our study examines extensive flux measurements on a tower in the Baltic Sea that allows partitioning of the fluxes according to wind direction without seriously depleting the data for a given wind-direction sector. Because some of the normal assumptions for the usual parameterization are not met, our study examines relationships behind the parameterization of the surface fluxes.
Abstract
Internal waves are predominantly generated by winds, tide–topography interactions, and balanced flow–topography interactions. Observations of vertical shear of horizontal velocity (uz
, υz
) from lowered acoustic Doppler current profilers (LADCP) profiles conducted during GO-SHIP hydrographic surveys, as well as vessel-mounted sonars, are used to interpret these signals. Vertical directionality of intermediate-wavenumber [
Abstract
Internal waves are predominantly generated by winds, tide–topography interactions, and balanced flow–topography interactions. Observations of vertical shear of horizontal velocity (uz
, υz
) from lowered acoustic Doppler current profilers (LADCP) profiles conducted during GO-SHIP hydrographic surveys, as well as vessel-mounted sonars, are used to interpret these signals. Vertical directionality of intermediate-wavenumber [
Abstract
Tide-induced near-inertial internal waves (NIWs) are generated by tide–topography interaction and are energized by internal tides through triadic resonant interaction of internal waves. They are located above topography and could be in close contact with wind-induced NIWs when the topography is a tall ridge, like in the Luzon Strait of the northern South China Sea (SCS). A natural question arises as to whether there is significant interaction between wind- and tide-induced NIWs. By using moored velocity observations, a satellite-tracked surface drifter dataset, and idealized numerical simulations, we find that in the presence of tide-induced NIWs, the wind can inject slightly more near-inertial energy (NIE), while in the presence of wind-induced NIWs, significantly more tidal energy is transferred to NIWs. Thus, wind- and tide-induced NIWs can mutually enhance each other, producing more NIE than a linear superposition of that generated by wind and tide forcing alone. Increasing wind intensity and tidal excursion lead to saturation of NIE enhancement, while a taller ridge leads to stronger enhancement. The high mixed layer NIE near Luzon Strait is mostly generated by the wind, while the mutual enhancement between wind- and tide-induced NIWs can further enhance this pattern. The interaction between wind- and tide-induced NIWs leads to an enhancement of 25% more NIE. If tide-induced NIWs are neglected, as is usually the case in the estimation of NIE, the total NIE will be underestimated by almost 50%. This might imply that tide-induced NIWs are important for the energetics of NIWs in Luzon Strait.
Significance Statement
Near-inertial internal waves (NIWs) usually occupy the most kinetic energy of internal waves and contribute significantly to ocean mixing. Near the surface they are usually generated by wind forcing, but near the bottom they can be generated by geostrophic or tidal flow interacting with topography. Above the tall ridge in Luzon Strait, wind- and tide-induced NIWs are in close contact, leading to potential interactions. It is found that these NIWs can mutually enhance each other, with most of the additional near-inertial energy (NIE) coming from the tides. If tide-induced NIWs are neglected, the total NIE will be underestimated by almost 50%. This suggests that tide-induced NIWs are important for the energetics of NIWs in Luzon Strait.
Abstract
Tide-induced near-inertial internal waves (NIWs) are generated by tide–topography interaction and are energized by internal tides through triadic resonant interaction of internal waves. They are located above topography and could be in close contact with wind-induced NIWs when the topography is a tall ridge, like in the Luzon Strait of the northern South China Sea (SCS). A natural question arises as to whether there is significant interaction between wind- and tide-induced NIWs. By using moored velocity observations, a satellite-tracked surface drifter dataset, and idealized numerical simulations, we find that in the presence of tide-induced NIWs, the wind can inject slightly more near-inertial energy (NIE), while in the presence of wind-induced NIWs, significantly more tidal energy is transferred to NIWs. Thus, wind- and tide-induced NIWs can mutually enhance each other, producing more NIE than a linear superposition of that generated by wind and tide forcing alone. Increasing wind intensity and tidal excursion lead to saturation of NIE enhancement, while a taller ridge leads to stronger enhancement. The high mixed layer NIE near Luzon Strait is mostly generated by the wind, while the mutual enhancement between wind- and tide-induced NIWs can further enhance this pattern. The interaction between wind- and tide-induced NIWs leads to an enhancement of 25% more NIE. If tide-induced NIWs are neglected, as is usually the case in the estimation of NIE, the total NIE will be underestimated by almost 50%. This might imply that tide-induced NIWs are important for the energetics of NIWs in Luzon Strait.
Significance Statement
Near-inertial internal waves (NIWs) usually occupy the most kinetic energy of internal waves and contribute significantly to ocean mixing. Near the surface they are usually generated by wind forcing, but near the bottom they can be generated by geostrophic or tidal flow interacting with topography. Above the tall ridge in Luzon Strait, wind- and tide-induced NIWs are in close contact, leading to potential interactions. It is found that these NIWs can mutually enhance each other, with most of the additional near-inertial energy (NIE) coming from the tides. If tide-induced NIWs are neglected, the total NIE will be underestimated by almost 50%. This suggests that tide-induced NIWs are important for the energetics of NIWs in Luzon Strait.
Abstract
Small-scale mixing drives the diabatic upwelling that closes the abyssal ocean overturning circulation. Indirect microstructure measurements of in situ turbulence suggest that mixing is bottom enhanced over rough topography, implying downwelling in the interior and stronger upwelling in a sloping bottom boundary layer. Tracer release experiments (TREs), in which inert tracers are purposefully released and their dispersion is surveyed over time, have been used to independently infer turbulent diffusivities—but typically provide estimates in excess of microstructure ones. In an attempt to reconcile these differences, Ruan and Ferrari derived exact tracer-weighted buoyancy moment diagnostics, which we here apply to quasi-realistic simulations. A tracer’s diapycnal displacement rate is exactly twice the tracer-averaged buoyancy velocity, itself a convolution of an asymmetric upwelling/downwelling dipole. The tracer’s diapycnal spreading rate, however, involves both the expected positive contribution from the tracer-averaged in situ diffusion as well as an additional nonlinear diapycnal distortion term, which is caused by correlations between buoyancy and the buoyancy velocity, and can be of either sign. Distortion is generally positive (stretching) due to bottom-enhanced mixing in the stratified interior but negative (contraction) near the bottom. Our simulations suggest that these two effects coincidentally cancel for the Brazil Basin Tracer Release Experiment, resulting in negligible net distortion. By contrast, near-bottom tracers experience leading-order distortion that varies in time. Errors in tracer moments due to realistically sparse sampling are generally small (<20%), especially compared to the
Abstract
Small-scale mixing drives the diabatic upwelling that closes the abyssal ocean overturning circulation. Indirect microstructure measurements of in situ turbulence suggest that mixing is bottom enhanced over rough topography, implying downwelling in the interior and stronger upwelling in a sloping bottom boundary layer. Tracer release experiments (TREs), in which inert tracers are purposefully released and their dispersion is surveyed over time, have been used to independently infer turbulent diffusivities—but typically provide estimates in excess of microstructure ones. In an attempt to reconcile these differences, Ruan and Ferrari derived exact tracer-weighted buoyancy moment diagnostics, which we here apply to quasi-realistic simulations. A tracer’s diapycnal displacement rate is exactly twice the tracer-averaged buoyancy velocity, itself a convolution of an asymmetric upwelling/downwelling dipole. The tracer’s diapycnal spreading rate, however, involves both the expected positive contribution from the tracer-averaged in situ diffusion as well as an additional nonlinear diapycnal distortion term, which is caused by correlations between buoyancy and the buoyancy velocity, and can be of either sign. Distortion is generally positive (stretching) due to bottom-enhanced mixing in the stratified interior but negative (contraction) near the bottom. Our simulations suggest that these two effects coincidentally cancel for the Brazil Basin Tracer Release Experiment, resulting in negligible net distortion. By contrast, near-bottom tracers experience leading-order distortion that varies in time. Errors in tracer moments due to realistically sparse sampling are generally small (<20%), especially compared to the
Abstract
The abyssal overturning circulation is thought to be primarily driven by small-scale turbulent mixing. Diagnosed water-mass transformations are dominated by rough topography “hotspots,” where the bottom enhancement of mixing causes the diffusive buoyancy flux to diverge, driving widespread downwelling in the interior—only to be overwhelmed by an even stronger upwelling in a thin bottom boundary layer (BBL). These water-mass transformations are significantly underestimated by one-dimensional (1D) sloping boundary layer solutions, suggesting the importance of three-dimensional physics. Here, we use a hierarchy of models to generalize this 1D boundary layer approach to three-dimensional eddying flows over realistically rough topography. When applied to the Mid-Atlantic Ridge in the Brazil Basin, the idealized simulation results are roughly consistent with available observations. Integral buoyancy budgets isolate the physical processes that contribute to realistically strong BBL upwelling. The downward diffusion of buoyancy is primarily balanced by upwelling along the sloping canyon sidewalls and the surrounding abyssal hills. These flows are strengthened by the restratifying effects of submesoscale baroclinic eddies and by the blocking of along-ridge thermal wind within the canyon. Major topographic sills block along-thalweg flows from restratifying the canyon trough, resulting in the continual erosion of the trough’s stratification. We propose simple modifications to the 1D boundary layer model that approximate each of these three-dimensional effects. These results provide local dynamical insights into mixing-driven abyssal overturning, but a complete theory will also require the nonlocal coupling to the basin-scale circulation.
Abstract
The abyssal overturning circulation is thought to be primarily driven by small-scale turbulent mixing. Diagnosed water-mass transformations are dominated by rough topography “hotspots,” where the bottom enhancement of mixing causes the diffusive buoyancy flux to diverge, driving widespread downwelling in the interior—only to be overwhelmed by an even stronger upwelling in a thin bottom boundary layer (BBL). These water-mass transformations are significantly underestimated by one-dimensional (1D) sloping boundary layer solutions, suggesting the importance of three-dimensional physics. Here, we use a hierarchy of models to generalize this 1D boundary layer approach to three-dimensional eddying flows over realistically rough topography. When applied to the Mid-Atlantic Ridge in the Brazil Basin, the idealized simulation results are roughly consistent with available observations. Integral buoyancy budgets isolate the physical processes that contribute to realistically strong BBL upwelling. The downward diffusion of buoyancy is primarily balanced by upwelling along the sloping canyon sidewalls and the surrounding abyssal hills. These flows are strengthened by the restratifying effects of submesoscale baroclinic eddies and by the blocking of along-ridge thermal wind within the canyon. Major topographic sills block along-thalweg flows from restratifying the canyon trough, resulting in the continual erosion of the trough’s stratification. We propose simple modifications to the 1D boundary layer model that approximate each of these three-dimensional effects. These results provide local dynamical insights into mixing-driven abyssal overturning, but a complete theory will also require the nonlocal coupling to the basin-scale circulation.
Abstract
The diabatic transformations of the middepth meridional overturning circulation (MOC) as it exits and reenters the South Atlantic to close the AMOC are studied using a state estimate assimilating data into a dynamically consistent ocean model. Virtual Lagrangian parcels in the lower branch of the MOC are followed in their global tour as they return to the upper branch of the MOC. Three return pathways are identified. The first pathway enters the abyssal Indo-Pacific as Circumpolar Deep Water, directly from the northern Antarctic Circumpolar Current (ACC), and before sampling the Antarctic margin. The second pathway sinks to abyssal densities exclusively in the Southern Ocean, then upwells while circulating within the ACC and eventually enters the Indo-Pacific or Atlantic at mid- to upper depths. The third pathway never reaches densities in the abyssal range. Parcels sinking in the Antarctic Bottom Water range upwell to mid- to upper depths south of 55°S. Parcels in all three pathways experience additional diabatic transformations after upwelling in the Southern Ocean, with more diabatic changes north of about 30°S than elsewhere. Diabatic changes are predominantly in the mixed layer of the tropical and subpolar gyres, enabled by Ekman suction. A simple model of the wind-driven flow illustrates that parcels always reach the surface in the tropical and subpolar gyres, regardless of their initial condition, because of coupling among gyres, the Ekman transport, and its return.
Abstract
The diabatic transformations of the middepth meridional overturning circulation (MOC) as it exits and reenters the South Atlantic to close the AMOC are studied using a state estimate assimilating data into a dynamically consistent ocean model. Virtual Lagrangian parcels in the lower branch of the MOC are followed in their global tour as they return to the upper branch of the MOC. Three return pathways are identified. The first pathway enters the abyssal Indo-Pacific as Circumpolar Deep Water, directly from the northern Antarctic Circumpolar Current (ACC), and before sampling the Antarctic margin. The second pathway sinks to abyssal densities exclusively in the Southern Ocean, then upwells while circulating within the ACC and eventually enters the Indo-Pacific or Atlantic at mid- to upper depths. The third pathway never reaches densities in the abyssal range. Parcels sinking in the Antarctic Bottom Water range upwell to mid- to upper depths south of 55°S. Parcels in all three pathways experience additional diabatic transformations after upwelling in the Southern Ocean, with more diabatic changes north of about 30°S than elsewhere. Diabatic changes are predominantly in the mixed layer of the tropical and subpolar gyres, enabled by Ekman suction. A simple model of the wind-driven flow illustrates that parcels always reach the surface in the tropical and subpolar gyres, regardless of their initial condition, because of coupling among gyres, the Ekman transport, and its return.
