Search Results
You are looking at 1  10 of 41 items for :
 Author or Editor: Paola Cessi x
 Journal of Physical Oceanography x
 Refine by Access: All Content x
Abstract
A parameterization for eddy buoyancy fluxes for use in coarsegrid models is developed and tested against eddyresolving simulations. The development is based on the assumption that the eddies are adiabatic (except near the surface) and the observation that the flux of buoyancy is affected by barotropic, depthindependent eddies. Like the previous parameterizations of Gent and McWilliams (GM) and Visbeck et al. (VMHS), the horizontal flux of a tracer is proportional to the local largescale horizontal gradient of the tracer through a transfer coefficient assumed to be given by the product of a typical eddy velocity scale and a typical mixing length. The proposed parameterization differs from GM and VMHS in the selection of the eddy velocity scale, which is based on the kinetic energy balance of baroclinic eddies. The three parameterizations are compared to eddyresolving computations in a variety of forcing configurations and for several sets of parameters. The VMHS and the energy balance parameterizations perform best in the tests considered here.
Abstract
A parameterization for eddy buoyancy fluxes for use in coarsegrid models is developed and tested against eddyresolving simulations. The development is based on the assumption that the eddies are adiabatic (except near the surface) and the observation that the flux of buoyancy is affected by barotropic, depthindependent eddies. Like the previous parameterizations of Gent and McWilliams (GM) and Visbeck et al. (VMHS), the horizontal flux of a tracer is proportional to the local largescale horizontal gradient of the tracer through a transfer coefficient assumed to be given by the product of a typical eddy velocity scale and a typical mixing length. The proposed parameterization differs from GM and VMHS in the selection of the eddy velocity scale, which is based on the kinetic energy balance of baroclinic eddies. The three parameterizations are compared to eddyresolving computations in a variety of forcing configurations and for several sets of parameters. The VMHS and the energy balance parameterizations perform best in the tests considered here.
Abstract
The current paradigm for the meridional overturning cell and the associated middepth stratification is that the wind stress in the subpolar region of the Southern Ocean drives a northward Ekman flow, which, together with the global diapycnal mixing across the lower boundary of the middepth waters, feeds the upper branch of the interhemispheric overturning. The resulting mass transport proceeds to the Northern Hemisphere of the North Atlantic, where it sinks, to be eventually returned to the Southern Ocean at depth. Seemingly, the wind stress in the Atlantic basin plays no role. This asymmetry occurs because the Ekman transport in the Atlantic Ocean is assumed to return geostrophically at depths much shallower than those occupied by the interhemispheric overturning. However, this vertical separation fails in the North Atlantic subpolar gyre region. Using a conceptual model and an ocean general circulation model in an idealized geometry, we show that the westerly wind stress in the northern part of the Atlantic provides two opposing effects. Mechanically, the return of the Ekman transport in the North Atlantic opposes sinking in this region, reducing the total overturning and deepening the middepth stratification; thermodynamically, the subpolar gyre advects salt poleward, promoting Northern Hemisphere sinking. Depending on which mechanism prevails, increased westerly winds in the Northern Hemisphere can reduce or augment the overturning.
Abstract
The current paradigm for the meridional overturning cell and the associated middepth stratification is that the wind stress in the subpolar region of the Southern Ocean drives a northward Ekman flow, which, together with the global diapycnal mixing across the lower boundary of the middepth waters, feeds the upper branch of the interhemispheric overturning. The resulting mass transport proceeds to the Northern Hemisphere of the North Atlantic, where it sinks, to be eventually returned to the Southern Ocean at depth. Seemingly, the wind stress in the Atlantic basin plays no role. This asymmetry occurs because the Ekman transport in the Atlantic Ocean is assumed to return geostrophically at depths much shallower than those occupied by the interhemispheric overturning. However, this vertical separation fails in the North Atlantic subpolar gyre region. Using a conceptual model and an ocean general circulation model in an idealized geometry, we show that the westerly wind stress in the northern part of the Atlantic provides two opposing effects. Mechanically, the return of the Ekman transport in the North Atlantic opposes sinking in this region, reducing the total overturning and deepening the middepth stratification; thermodynamically, the subpolar gyre advects salt poleward, promoting Northern Hemisphere sinking. Depending on which mechanism prevails, increased westerly winds in the Northern Hemisphere can reduce or augment the overturning.
Abstract
It is well established that the mean transport through Bering Strait is balanced by a sea level difference between the North Pacific and the Arctic Ocean, but no mechanism has been proposed to explain this sea level difference. It is argued that the sea level difference across Bering Strait, which geostrophically balances the northward throughflow, is associated with the sea level difference between the North Pacific and the North Atlantic/Arctic. In turn, the latter difference is caused by deeper middepth isopycnals in the IndoPacific than in the Atlantic, especially in the northern high latitudes because there is deep water formation in the Atlantic, but not in the Pacific. Because the depth of the middepth isopycnals is associated with the dynamics of the upper branch of the meridional overturning circulation (MOC), a model is formulated that quantitatively relates the sea level difference between the North Pacific and the Arctic/North Atlantic with the wind stress in the Antarctic Circumpolar region, since this forcing powers the MOC, and with the outcropping isopycnals shared between the Northern Hemisphere and the Antarctic circumpolar region, since this controls the location of deep water formation. This implies that if the sinking associated with the MOC were to occur in the North Pacific, rather than the North Atlantic, then the Bering Strait flow would reverse. These predictions, formalized in a theoretical box model, are confirmed by a series of numerical experiments in a simplified geometry of the World Ocean, forced by steady surface wind stress, temperature, and freshwater flux.
Abstract
It is well established that the mean transport through Bering Strait is balanced by a sea level difference between the North Pacific and the Arctic Ocean, but no mechanism has been proposed to explain this sea level difference. It is argued that the sea level difference across Bering Strait, which geostrophically balances the northward throughflow, is associated with the sea level difference between the North Pacific and the North Atlantic/Arctic. In turn, the latter difference is caused by deeper middepth isopycnals in the IndoPacific than in the Atlantic, especially in the northern high latitudes because there is deep water formation in the Atlantic, but not in the Pacific. Because the depth of the middepth isopycnals is associated with the dynamics of the upper branch of the meridional overturning circulation (MOC), a model is formulated that quantitatively relates the sea level difference between the North Pacific and the Arctic/North Atlantic with the wind stress in the Antarctic Circumpolar region, since this forcing powers the MOC, and with the outcropping isopycnals shared between the Northern Hemisphere and the Antarctic circumpolar region, since this controls the location of deep water formation. This implies that if the sinking associated with the MOC were to occur in the North Pacific, rather than the North Atlantic, then the Bering Strait flow would reverse. These predictions, formalized in a theoretical box model, are confirmed by a series of numerical experiments in a simplified geometry of the World Ocean, forced by steady surface wind stress, temperature, and freshwater flux.
Abstract
An inertial gyre with characteristics very similar to the recirculation observed in eddyresolving general circulation models is obtained with a simple, analytically tractable, twolayer model. The recirculating gyre is contained in a box of simple geometry, which isolates it from the sverdrup interior. The gyre is forced by prescribing anomalous values of potential vorticity northward in the subtropical gyre or can be though of as a rough parametrization of diabatic forcing. In both cases the forcing is confined to the water above the thermocline, which is represented by the upper layer and is transmitted to the abyssal ocean through interfacial friction.
The condition for the abyssal water to be set in motion, is derived and for oceanic values the recirculation goes all the way to the bottom. When this occurs the center of the gyre is dominated by a barotropic flow, while the baroclinic flow is confined to the edges of the gyre. The width and strength of the gyre can be easily calculated in the limit of long, narrow gyres. The meridional scale of the gyre is directly proportional to the vorticity anomaly injected at the northern boundary, and the barotropic part of the transport is proportional to the cube of the abyssal gyre width, in close analogy with the results found by Cessi, Ierley and Young in a one layer model.
Abstract
An inertial gyre with characteristics very similar to the recirculation observed in eddyresolving general circulation models is obtained with a simple, analytically tractable, twolayer model. The recirculating gyre is contained in a box of simple geometry, which isolates it from the sverdrup interior. The gyre is forced by prescribing anomalous values of potential vorticity northward in the subtropical gyre or can be though of as a rough parametrization of diabatic forcing. In both cases the forcing is confined to the water above the thermocline, which is represented by the upper layer and is transmitted to the abyssal ocean through interfacial friction.
The condition for the abyssal water to be set in motion, is derived and for oceanic values the recirculation goes all the way to the bottom. When this occurs the center of the gyre is dominated by a barotropic flow, while the baroclinic flow is confined to the edges of the gyre. The width and strength of the gyre can be easily calculated in the limit of long, narrow gyres. The meridional scale of the gyre is directly proportional to the vorticity anomaly injected at the northern boundary, and the barotropic part of the transport is proportional to the cube of the abyssal gyre width, in close analogy with the results found by Cessi, Ierley and Young in a one layer model.
Abstract
A simple mechanism for the generation of flow in the eastern region of the oceans is presented. It is shown that even a small amount of upwellingfavorable wind stress parallel to the boundary makes the subducted geostrophic contours veer into the eastern boundary. This sweep toward the wall ventilates the top few hundred meters of the eastern subtropical gyre. The latitudinal extent of the eastern ventilation is calculated. Along a meridional boundary the subsurface flow is eastward and appears to provide the mass source for the poleward undercurrents observed at the eastern boundary of the ocean.
Abstract
A simple mechanism for the generation of flow in the eastern region of the oceans is presented. It is shown that even a small amount of upwellingfavorable wind stress parallel to the boundary makes the subducted geostrophic contours veer into the eastern boundary. This sweep toward the wall ventilates the top few hundred meters of the eastern subtropical gyre. The latitudinal extent of the eastern ventilation is calculated. Along a meridional boundary the subsurface flow is eastward and appears to provide the mass source for the poleward undercurrents observed at the eastern boundary of the ocean.
Abstract
In an earlier paper it was shown that the presence of bottom topography can substantially alter the Sverdrup prediction for the largescale wind driven flow. In their model the abyssal water is set in motion by the eddy field, parameterized as lateral diffusion of potential vorticity. If the topography has a structure in the eastwest direction, then the solution found in the inviscid limit by Cessi and Pedlosky predicts the occurrence of strong jets in the interior of the model ocean. In this note I used a numerical model to test whether the jets predicted by the analytic solutions survive when inertia and diffusion are included explicitly.
In the inviscid limit, according to the structure of the topography these internal jets can occur in both vertically homogeneous and stratified models. Specifically, if the topographic slope changes sign, then one kind of jets is observed both in baroclinic and barotropic models. In the present work it is shown that this phenomenon is observed with moderate amounts of diffusion and is not disturbed by the occurrence of recirculating inertial gyres within the basin.
If the topographic slope is constant, then another kind of internal jets is theoretically predicted in the inviscid limit, and it occurs in stratified models only. Numerical calculations did not exhibit this kind of internal jets in the presence of inertia and weak diffusion, and the reasons for this failure are rationalized.
Abstract
In an earlier paper it was shown that the presence of bottom topography can substantially alter the Sverdrup prediction for the largescale wind driven flow. In their model the abyssal water is set in motion by the eddy field, parameterized as lateral diffusion of potential vorticity. If the topography has a structure in the eastwest direction, then the solution found in the inviscid limit by Cessi and Pedlosky predicts the occurrence of strong jets in the interior of the model ocean. In this note I used a numerical model to test whether the jets predicted by the analytic solutions survive when inertia and diffusion are included explicitly.
In the inviscid limit, according to the structure of the topography these internal jets can occur in both vertically homogeneous and stratified models. Specifically, if the topographic slope changes sign, then one kind of jets is observed both in baroclinic and barotropic models. In the present work it is shown that this phenomenon is observed with moderate amounts of diffusion and is not disturbed by the occurrence of recirculating inertial gyres within the basin.
If the topographic slope is constant, then another kind of internal jets is theoretically predicted in the inviscid limit, and it occurs in stratified models only. Numerical calculations did not exhibit this kind of internal jets in the presence of inertia and weak diffusion, and the reasons for this failure are rationalized.
Abstract
Welander's flip–flop model exhibits oscillations when forced by stochastic white noise (with zero mean) even in the region of parameters where the deterministic system has a globally stable fixed point. Perturbations away from the attracting solutions decay exponentially in time, without any oscillation. Thus, the oscillation that appears when the system is stochastically forced is not related to an eigenfrequency of the linearized system.
The characteristics of this noiseinduced oscillation are contrasted with those obtained when the damped harmonic oscillator is forced stochastically. In the case of a stochastically forced damped harmonic oscillator the spectral peak coincides with the frequency of the oscillator, and the amplitude of the oscillations is proportional to the strength of the noise. In the stochastically forced flipflop model the amplitude of the oscillations is independent of the strength of the noise and the spectral peak moves to lower frequencies as the amplitude of the noise is reduced. Moreover, for noise below a critical threshold, no spectral peak is obtained.
The flip–flop model shares four characteristics with the thermohaline oscillations observed in OGCMs:

The freshwater flux determines whether the system oscillates or settles into a steady state. The period of the oscillations is very sensitive to the freshwater flux and becomes arbitrarily long near the transition from steady to periodic behavior.

The oscillations are of finite amplitude even just past the threshold value of the freshwater flux that separates periodic behavior from a steady equilibrium.

One extremum of the oscillation excursion is close to the value of the steady state that exists below the threshold for transition.

When the deterministic system reaches a steady state, oscillations can be excited by adding a stochastic component to the freshwater flux. The period of the resulting oscillations decreases with increasing noise amplitude, while the amplitude of the oscillations is insensitive to the amplitude of the noise.
Abstract
Welander's flip–flop model exhibits oscillations when forced by stochastic white noise (with zero mean) even in the region of parameters where the deterministic system has a globally stable fixed point. Perturbations away from the attracting solutions decay exponentially in time, without any oscillation. Thus, the oscillation that appears when the system is stochastically forced is not related to an eigenfrequency of the linearized system.
The characteristics of this noiseinduced oscillation are contrasted with those obtained when the damped harmonic oscillator is forced stochastically. In the case of a stochastically forced damped harmonic oscillator the spectral peak coincides with the frequency of the oscillator, and the amplitude of the oscillations is proportional to the strength of the noise. In the stochastically forced flipflop model the amplitude of the oscillations is independent of the strength of the noise and the spectral peak moves to lower frequencies as the amplitude of the noise is reduced. Moreover, for noise below a critical threshold, no spectral peak is obtained.
The flip–flop model shares four characteristics with the thermohaline oscillations observed in OGCMs:

The freshwater flux determines whether the system oscillates or settles into a steady state. The period of the oscillations is very sensitive to the freshwater flux and becomes arbitrarily long near the transition from steady to periodic behavior.

The oscillations are of finite amplitude even just past the threshold value of the freshwater flux that separates periodic behavior from a steady equilibrium.

One extremum of the oscillation excursion is close to the value of the steady state that exists below the threshold for transition.

When the deterministic system reaches a steady state, oscillations can be excited by adding a stochastic component to the freshwater flux. The period of the resulting oscillations decreases with increasing noise amplitude, while the amplitude of the oscillations is insensitive to the amplitude of the noise.
Abstract
A modified Stommel twobox model is considered as a minimal representation of the buoyancydriven ocean circulation. In the limit of fast temperature relaxation only the salinity evolves in time while the temperature is clamped to the prescribed ambient value. The box model has no intrinsic variability: just two linearly stable and one unstable equilibria. A finite perturbation is needed to shift the system from one stable equilibrium to the other. The minimum amplitude and duration in time of the perturbation are calculated.
A stochastic component of the freshwater flux forcing is then added to model the effect of changes in the global hydrological cycle due to the “weather.” The stochastic forcing is a source of extrinsic time dependence. The salinity gradient obeys an equation analogous to the trajectory of a viscous particle in a doublewelled potential, subject to Brownian agitation. If the amplitude of the stochastic driving is above a certain threshold, then there is a finite probability of switching from one stable equilibrium to the other. The threshold variance and the average residence time in each equilibrium are calculated. For timescales on the order of the average residence time or longer, the box model behaves like a random telegraph process.
The stochastic driving also induces a “rattling” around each steady equilibrium whose frequency is proportional to the curvature of the potential well at each equilibrium. The probability of being in each well can be calculated and, within each equilibrium, the box model behaves like an OrnsteinUhlenbeck process.
Finally the spectrum of the salinity gradient is calculated analytically using standard approximations in stochastic processes. The approximate analytical results are in excellent agreement with those obtained by direct computation.
Abstract
A modified Stommel twobox model is considered as a minimal representation of the buoyancydriven ocean circulation. In the limit of fast temperature relaxation only the salinity evolves in time while the temperature is clamped to the prescribed ambient value. The box model has no intrinsic variability: just two linearly stable and one unstable equilibria. A finite perturbation is needed to shift the system from one stable equilibrium to the other. The minimum amplitude and duration in time of the perturbation are calculated.
A stochastic component of the freshwater flux forcing is then added to model the effect of changes in the global hydrological cycle due to the “weather.” The stochastic forcing is a source of extrinsic time dependence. The salinity gradient obeys an equation analogous to the trajectory of a viscous particle in a doublewelled potential, subject to Brownian agitation. If the amplitude of the stochastic driving is above a certain threshold, then there is a finite probability of switching from one stable equilibrium to the other. The threshold variance and the average residence time in each equilibrium are calculated. For timescales on the order of the average residence time or longer, the box model behaves like a random telegraph process.
The stochastic driving also induces a “rattling” around each steady equilibrium whose frequency is proportional to the curvature of the potential well at each equilibrium. The probability of being in each well can be calculated and, within each equilibrium, the box model behaves like an OrnsteinUhlenbeck process.
Finally the spectrum of the salinity gradient is calculated analytically using standard approximations in stochastic processes. The approximate analytical results are in excellent agreement with those obtained by direct computation.
Abstract
The role of the relative geometry of mechanical forcing (wind stress) and buoyancy forcing (prescribed surface temperature) in the maintenance of the main thermocline is explored. In particular, the role of the wind stress curl in enhancing or suppressing the generation of baroclinic eddies is studied in simplified domains. The dependence of key quantities, such as the depth of the thermocline and the maximum heat transport, on the external parameters such as diapycnal mixing and dissipation rate is examined. Qualitatively different regimes are found depending on the relative phase of the wind stress and surface buoyancy distribution. The most efficient arrangement for eddy generation has Ekman pumping (suction) in conjunction with high (low) surface buoyancy. This corresponds to the situation found in the midlatitudes, where the surface Ekman flow carries heat toward the warmer region (i.e., upgradient of the surface temperature). In this case, strong eddy fluxes are generated in order to counteract the upgradient heat transport by the Ekman cell. The result is a thermocline whose depth is independent of the diapycnal diffusivity. However, the competition between these opposing heat fluxes leads to a weak net heat transport, proportional to the diffusivity responsible for the diabatic forcing. This arrangement of wind stress provides a large source of available potential energy on which eddies can grow, so the mechanical energy balance for the eddies is consistent with a substantial eddy heat flux. When the same surface temperature distribution is paired with the opposite wind stress curl, the mean flow produces a sink, rather than a source, of available potential energy and eddies are suppressed. With this arrangement, typical of low latitudes and the subpolar regions, the Ekman overturning cell carries heat downgradient of the surface temperature. Thus, the net heat transport is almost entirely due to the Ekman flow and is independent of the diapycnal diffusivity. At the same time the thermocline is a thin, diffusive boundary layer. Quantitative scalings for the thermocline depth and the poleward heat transport in these two limiting cases are contrasted and successfully compared with eddyresolving computations.
Abstract
The role of the relative geometry of mechanical forcing (wind stress) and buoyancy forcing (prescribed surface temperature) in the maintenance of the main thermocline is explored. In particular, the role of the wind stress curl in enhancing or suppressing the generation of baroclinic eddies is studied in simplified domains. The dependence of key quantities, such as the depth of the thermocline and the maximum heat transport, on the external parameters such as diapycnal mixing and dissipation rate is examined. Qualitatively different regimes are found depending on the relative phase of the wind stress and surface buoyancy distribution. The most efficient arrangement for eddy generation has Ekman pumping (suction) in conjunction with high (low) surface buoyancy. This corresponds to the situation found in the midlatitudes, where the surface Ekman flow carries heat toward the warmer region (i.e., upgradient of the surface temperature). In this case, strong eddy fluxes are generated in order to counteract the upgradient heat transport by the Ekman cell. The result is a thermocline whose depth is independent of the diapycnal diffusivity. However, the competition between these opposing heat fluxes leads to a weak net heat transport, proportional to the diffusivity responsible for the diabatic forcing. This arrangement of wind stress provides a large source of available potential energy on which eddies can grow, so the mechanical energy balance for the eddies is consistent with a substantial eddy heat flux. When the same surface temperature distribution is paired with the opposite wind stress curl, the mean flow produces a sink, rather than a source, of available potential energy and eddies are suppressed. With this arrangement, typical of low latitudes and the subpolar regions, the Ekman overturning cell carries heat downgradient of the surface temperature. Thus, the net heat transport is almost entirely due to the Ekman flow and is independent of the diapycnal diffusivity. At the same time the thermocline is a thin, diffusive boundary layer. Quantitative scalings for the thermocline depth and the poleward heat transport in these two limiting cases are contrasted and successfully compared with eddyresolving computations.
Abstract
The objective of this study is to investigate the timedependent circulation in a closed basin where the steady circulation is included and long Rossby wave speeds are consistent with observations. Specifically, the largescale baroclinic eigenmodes of a twolayer rectangular basin forced by surface wind stress in the limit of small dissipation are examined. Lowfrequency modes with small decay rates independent of friction result when the constraint of mass conservation is enforced. The magnitude of the wind stress is found to be critical to the eigenspectrum. For all forcing magnitudes, including forcings with closed geostrophic contours, oscillatory modes with decay rates independent of friction emerge. For forcings with closed geostrophic contours, two important classes of eigenmodes with comparable decay rates emerge: purely decaying modes confined to the region of closed contours, and basinscale oscillatory modes. The purely decaying modes also exist without the constraint of total mass conservation but their decay rates depend on the magnitude of friction to leading order.
Abstract
The objective of this study is to investigate the timedependent circulation in a closed basin where the steady circulation is included and long Rossby wave speeds are consistent with observations. Specifically, the largescale baroclinic eigenmodes of a twolayer rectangular basin forced by surface wind stress in the limit of small dissipation are examined. Lowfrequency modes with small decay rates independent of friction result when the constraint of mass conservation is enforced. The magnitude of the wind stress is found to be critical to the eigenspectrum. For all forcing magnitudes, including forcings with closed geostrophic contours, oscillatory modes with decay rates independent of friction emerge. For forcings with closed geostrophic contours, two important classes of eigenmodes with comparable decay rates emerge: purely decaying modes confined to the region of closed contours, and basinscale oscillatory modes. The purely decaying modes also exist without the constraint of total mass conservation but their decay rates depend on the magnitude of friction to leading order.