Abstract

A model of the convectively driven mixed layer is used to interpret BOMEX Phase III data. It is concluded that the compensating downdrafts and the moisture transports due to trade cumulus must be included explicitly in any theory of the dynamics of the tropical mixed layer.

Abstract

The fundamental modes of oscillation of a coupled atmosphere–ocean basin system in the presence of a spatially varying oceanic basic state are investigated by formulating and solving an eigenvalue problem, thereby extending the work of Hirst. The model reduces essentially to the linearized Zebiak and Cane model as discussed by Battisti and Hirst. With conventionally chosen basic states, the unstable eigenmode closely resembles the El Niño–Southern Oscillation (ENSO) cycle in these models.

It is shown that the unstable low-frequency eigenfunction consists primarily of a Kelvin mode and a gravest equatorial Rossby mode, and the oscillation can be understood in particularly simple term essentially those proposed by Suarez and Schopf and others. The oscillatory nature of the ENSO cycle can be explained by a transition mechanism resulting from the interaction of these two equatorial (but not necessarily propagating) modes. A growing unstable positive wind anomaly in the central Pacific produces a growing eastward-propagating downwelling Kelvin mode and a growing westward-propagating upwelling equatorial Rossby mode. The down-welling Kelvin mode propagates eastward and enhances the growing warm phase of the ENSO. On the other hand, the upwelling Rossby mode propagates westward and produces an upwelling Kelvin mode via rejection at the western boundary. This growing Kelvin mode propagates to the central and eastern Pacific where it then grows without propagation, cools the warm anomaly, eventually changes the phase of the warm event to cold, and therefore switches the sign of the air–sea coupled instability in the eastern Pacific. The regular ENSO cycle is the repeated application of this mechanism.

The nature of the propagation of the ENSO anomalies is shown to be sensitive to the meridional profile of the upwelling velocity near the equator. The sea surface temperature (SST) anomaly changes synchronously (i.e., without propagation) in the eastern Pacific only if the entrainment velocity is tightly confined meridionally to the equator, while it begins to propagate eastward if the entrainment velocity expands in the meridional direction, all other parameters held constant.

In examining the parameter dependence of the unstable modes, it was found that two nonoscillatory solutions appear as a transition from the oscillatory solution as the air–sea coupling parameter and the Rayleigh friction parameter of the ocean are increased.

Abstract

Arakawa's recent parameterization of the effects of a cumulus ensemble on the large-scale environment is applied to the problem of conditional instability of the second kind (CISK). In particular, Charney's linear, two-level, line-symmetry CISK model of the ITCZ is re-examined using a simplified non-entraining cloud version of the Arakawa scheme. It is found that the growth rate is maximum, in fact infinite, at some reasonable mesoscale rather than at cumulus scale as is characteristic of Charney's solution. A more accurate semi-analytic model of CISK is considered and it is found that a separable, line-symmetric CISK solution is always possible under very general conditions. In both the two-level and semi-analytic models of CISK, it is proved that a necessary condition for the existence of a growing solution is that the mass flux into the clouds exceeds the Ekman pumping out of the boundary layer, or equivalently, that the air between the clouds must subside and therefore heat the environment by adiabatic compression.

Abstract

The effect of nonlinearities on a previously investigated coupled atmosphere–ocean basin mode is examined. The nonlinearity in the thermodynamic equation for sea surface temperature arises mainly from the dependence of subsurface temperature on the thermocline depth anomaly in the parameterization of entrainment into the mixed layer. This nonlinearity ultimately suppresses the linear growth of the unstable mode and equilibrates it at a finite amplitude. Because this nonlinearity acts differently for warm and cold states, the warm states are enhanced at finite amplitude. It is found that multiple equilibrium states appear as the coupling coefficient increases and as the reflection coefficient of the oceanic Rossby mode at the western boundary decreases. The finite-amplitude warm equilibrium state turns out to be stable, but the finite-amplitude cold state is unstable. The explicit inclusion of the dependence of the coupling strength on the warm and cold sea surface temperature anomalies modulates the sinusoidal-like oscillation and increases the period, but aperiodic solutions could not be obtained.

Abstract

Several simple numerical experiments are conducted, using both single- and double-hemisphere ocean basins under symmetric steady forcing to study de ocean's thermohaline circulation. It is shown that a stable steady state obtained under a restoring surface boundary condition on salinity becomes unstable upon a switch to a flux boundary condition. The polar halocline catastrope of F. Bryan occurs. It is shown that further integration of this collapsed state ultimately yields a steady, stable one-cell circulation with the approach being essentially chaotic but with significant energy at decadal period. The two-hemisphere ocean passes through many stages in which violent overturning occurs O(80 × 10¹ m³ a⁻¹). These flushes occurs in both hemispheres and are of one-cell structure. The time period between them Bushes varies from seveal hundred to about one thousand years.

A single 12-vertical-level hemispheric basin, spun up from an initial state of rest under mixed boundary conditions (restoring boundary condition on temperature and flux boundary condition on salinity), never reaches a study gate. Three characteristic stages are observed in the integration: a stage where the system oscillates with decadal time scale, a stage when the system undergoes a violent overturning flush, and a Quiescent stage in which either deep water is forming or the themohaline circulation is in a collapsed state. These three characteristic stage are also present in 33 level single- and double-hemisphere runs. The decadal time wide is associated primarily with the advection of positive salinity anomalies into the region of deep-water formation from the midocean region between the subtropical and subpolar gyres. Upon increasing the resolution to 33 levels a steady is reached. The resulting steady state is fundamentally different from the one obtained under the same resolution and restoring boundary conditions in that it is more energetic and has much warmer basin mean temperature. These differences are due to a change in the location of deep-water formation.

The dependence of the results on the type a convection scheme used, vertical resolution and time-stepping procedure (synchronous or asynchronous integration) is also studied in order to separate physical processes from those that might be numerical artifacts. Sufficient vertical resolution is shown to be important in obtaining realistic models of the thermohaline circulation. It is shown that a steady state, which is stable under asynchronous integration and mixed boundary conditions may become unstable upon a switch to synchronous integration. It is also shown that the steady state obtained under restoring boundary conditions only changes slightly upon a switch to synchronous integration. Under mixed boundary conditions the steady state is shown to be very sensitive to the choice of surface tracer time step even while integrating asynchronously. Upon a Switch in this time step a polar halocline catastrophe way be induced.

The implications of the present study for future ocean climate modles are discussed.

Abstract

A series of numerical experiments is conducted with a three-dimensional ocean general circulation model and a two-dimensional counterpart both designed for efficient integration over diffusive (millennial) time scales. With strong steady salinity fluxes (salting at low latitudes and freshening at high), basin mean temperature and several other diagnostics show a series of self-sustaining oscillations. The oscillations termed deep decoupling oscillations, exhibit halocline catastrophes at regular intervals, followed by warming deep decoupled phases (when the deep overturning is weak), cooling flushes, and in the lower range of salinity forcing, a coupled phase when the deep ocean advective/diffusive heat balance is almost, but not quite, met. It is suggested that oscillations arise when a steady overturning circulation encounters a contradiction: the poleward salt and heat transport needed to maintain convection in the polar ocean requires more overturning than is consistent with the reduced thermocline depth that results. This hypothesis is supported by the sensitivity to variations in the vertical diffusivity: increased vertical diffusivity stabilizes oscillating solutions into steady, thermally direct circulations.

Although deep decoupling oscillations appear in both two- and three-dimensional models, they occur over a much broader range of forcing in the three-dimensional model. This is shown to be due to heat and salt transports by the horizontal plane (gyre) motions in the three-dimensional model that intensify in the upper polar ocean in response to the formation of a halocline and eventually destabilize it. Increasing the wind stress in the three-dimensional model and the horizontal diffusivity in the two-dimensional model stabilizes oscillating solutions. The amplitude, shape, and period of the oscillations are also sensitive to the strength of the salinity forcing.

Another kind of oscillation, termed a loop oscillation, with a smaller amplitude and an overturning time scale, is found in some of the more weakly forced experiments with both models. These oscillations are shown to be a result of the advection of salinity anomalies by the deep overturning, affecting its strength in a manner that leads to their further amplification by feedback from the salinity flux boundary condition. A simple thermohaline loop model demonstrates the essential advective mechanism for this kind of oscillation.

Abstract

A two-dimensional (in a vertical and meridional plane) model for steady equatorial undercurrents is described. Compared to the primitive equation model, the zonal pressure gradient and associated zonal temperature gradients (both vary vertically) are prescribed in this model, and all other terms involving zonal variations are ignored. With zonal pressure gradients resembling actual ocean gradients, model undercurrents agree well with observations as far as the main features are concerned. In particular, the model simulates a stronger undercurrent in the Pacific than in the Atlantic, suggesting that a weaker zonal wind stress, a shallower thermocline, a more surface-confined zonal pressure gradient, and an associated larger magnitude of near-surface zonal temperature gradient around 30°W in the Atlantic than around 150°W in the Pacific, which is related to the longitudinal structure of the zonal wind stress and longitudinal basin extent, are the cause of this difference. An argument based on geostrophy and heat balance is also given.

The model is used to examine the dynamic nature and heat balance of steady equatorial undercurrents for a symmetric circulation about the equator. With a full, nonlinear heat balance, an undercurrent is generated in both linear and nonlinear dynamic balances, but the dynamical features are different in the two cases. In the nonlinear dynamic case, vertical-momentum transports play a key role; in the linear dynamic case, though the eastward zonal pressure gradient provides a necessary forcing, the existence of the undercurrent also relies on the meridional diffusive momentum transport near the surface, which is positive instead of negative. For a doubling of zonal wind stress and a fixed vertical profile of zonal pressure gradient, the speed of the undercurrent core increases by about 25% in the nonlinear case but remains unchanged in the linear case; surface temperature increases by about 1.3 K in the nonlinear case and decreases by 3 K in the linear case.

Within the undercurrent core, the dominant momentum balance is between the zonal pressure gradient and meridional diffusive friction, and the heat balance is between zonal and vertical advections. It is proposed that the position of the undercurrent core relative to the thermocline reflects different advective heat balances: the undercurrent core is above (or below) the thermocline if the net heat advection balance tends to heat (or cool). The fact that the undercurrent core is more or less in the thermocline suggests that three-dimensional advective heat transports almost cancel each other.

Abstract

No abstract available.

Abstract

In order to test certain aspects of the ENSO mechanism proposed by Suarez and Schopf and by Battisti and Hirst, we force a shallow water ocean with 28 years of observed (FSU) winds and decompose both the atmospheric forcing and ocean response into equatorial modes.

The proposed mechanism was verified in the following sense. For each warm phase of the cycle, a downwelling Kelvin mode directly forced by the weakening of the trade winds in the central Pacific was found. The same wind anomaly forces an upwelling gravest Rossby mode which propagates freely westward and reflects as a freely propagating Kelvin mode at the western boundary. This Kelvin mode returns to the scene of the original warming and acts as a retarded forcing (of opposite sign) eventually switching the phase of the oscillation from warm to cold. The cold phase then proceeds by the same mechanism but with modes of opposite signs, to eventually switch to the warm phase.

The contribution of the higher Rossby modes to this process is estimated. It is found that almost all the Kelvin response in the western Pacific results from the reflection of the gravest symmetric equatorial Rossby mode so that the ENSO cycle is defined by the interaction of only two ocean modes, the Kelvin and lowest Rossby mode. These modes are sometimes forced directly by the wind as part of an intrinsic coupled atmosphere–ocean mode and sometimes propagate freely where the wind forcing is negligible. The distinction between these two manifestations of the ocean modes is of the greatest importance and is stressed throughout the paper.

The case of the warm event of 1976 is an interesting example of the failure of the switching mechanism, and is discussed in the conclusion.

Abstract

In centered difference models of ocean circulation, two grid-point computational modes can be excited if grid Reynolds and Peclet numbers are greater than two. The Bryan-Cox General Circulation Model (GCM) is used to show the dramatic effect that this instability has on the equatorial thermohaline circulation. In many recent numerical calculations researchers have used 12 vertical levels. It is shown that this resolution produces an artificial cell at the equator when typical values of the vertical diffusivity and viscosity parameters are used. This artifical cell rotates counter to the primary cell driven by deep water formation at high latitudes, is driven by downwelling at the eastern boundary near the equator and is 40% the strength of the primary cell for the parameters used in the present study. When the vertical resolution is increased the cell vanishes. It is suggested therefore that higher vertical resolution should be used in Bryan-Cox GCM deep-ocean modeling studies when current values of the vertical diffusivity and viscosity parameters are used.