R. C. Pacanowski

The Pliocene Paradox (Mechanisms for a Permanent El Niño)

During the early Pliocene, 5 to 3 million years ago, globally averaged temperatures were substantially higher than they are today, even though the external factors that determine climate were essentially the same. In the tropics, El Niño was continual (or “permanent”) rather than intermittent. The appearance of northern continental glaciers, and of cold surface waters in oceanic upwelling zones in low latitudes (both coastal and equatorial), signaled the termination of those warm climate conditions and the end of permanent El Niño. This led to the amplification of obliquity (but not precession) cycles in equatorial sea surface temperatures and in global ice volume, with the former leading the latter by several thousand years. A possible explanation is that the gradual shoaling of the oceanic thermocline reached a threshold around 3 million years ago, when the winds started bringing cold waters to the surface in low latitudes. This introduced feedbacks involving ocean-atmosphere interactions that, along with ice-albedo feedbacks, amplified obliquity cycles. A future melting of glaciers, changes in the hydrological cycle, and a deepening of the thermocline could restore the warm conditions of the early Pliocene.

Simulated Climate and Climate Change in the GFDL CM2.5 High-Resolution Coupled Climate Model

AbstractThe authors present results for simulated climate and climate change from a newly developed high-resolution global climate model [Geophysical Fluid Dynamics Laboratory Climate Model version 2.5 (GFDL CM2.5)]. The GFDL CM2.5 has an atmospheric resolution of approximately 50 km in the horizontal, with 32 vertical levels. The horizontal resolution in the ocean ranges from 28 km in the tropics to 8 km at high latitudes, with 50 vertical levels. This resolution allows the explicit simulation of some mesoscale eddies in the ocean, particularly at lower latitudes.Analyses are presented based on the output of a 280-yr control simulation; also presented are results based on a 140-yr simulation in which atmospheric CO2 increases at 1% yr−1 until doubling after 70 yr.Results are compared to GFDL CM2.1, which has somewhat similar physics but a coarser resolution. The simulated climate in CM2.5 shows marked improvement over many regions, especially the tropics, including a reduction in the double ITCZ and an i...

GFDL's CM2 Global Coupled Climate Models. Part II: The Baseline Ocean Simulation

The current generation of coupled climate models run at the Geophysical Fluid Dynamics Laboratory (GFDL) as part of the Climate Change Science Program contains ocean components that differ in almost every respect from those contained in previous generations of GFDL climate models. This paper summarizes the new physical features of the models and examines the simulations that they produce. Of the two new coupled climate model versions 2.1 (CM2.1) and 2.0 (CM2.0), the CM2.1 model represents a major improvement over CM2.0 in most of the major oceanic features examined, with strikingly lower drifts in hydrographic fields such as temperature and salinity, more realistic ventilation of the deep ocean, and currents that are closer to their observed values. Regional analysis of the differences between the models highlights the importance of wind stress in determining the circulation, particularly in the Southern Ocean. At present, major errors in both models are associated with Northern Hemisphere Mode Waters and outflows from overflows, particularly the Mediterranean Sea and Red Sea.

Spurious Diapycnal Mixing Associated with Advection in a z-Coordinate Ocean Model

Abstract This paper discusses spurious diapycnal mixing associated with the transport of density in a z-coordinate ocean model. A general method, based on the work of Winters and collaborators, is employed for empirically diagnosing an effective diapycnal diffusivity corresponding to any numerical transport process. This method is then used to quantify the spurious mixing engendered by various numerical representations of advection. Both coarse and fine resolution examples are provided that illustrate the importance of adequately resolving the admitted scales of motion in order to maintain a small amount of mixing consistent with that measured within the ocean’s pycnocline. Such resolution depends on details of the advection scheme, momentum and tracer dissipation, and grid resolution. Vertical transport processes, such as convective adjustment, act as yet another means to increase the spurious mixing introduced by dispersive errors from numerical advective fluxes.

A Global Ocean-Atmosphere Climate Model. Part II. The Oceanic Circulation

Abstract A numerical experiment has been carried out with a joint model of the ocean and atmosphere. The 12-level model of the world ocean predicts the fields of horizontal velocity, temperature and salinity. It includes the effects of bottom topography, and a simplified model of polar pack ice. The numerical experiment allows the joint ocean-atmosphere model to seek an equilibrium over the equivalent of 270 years in the ocean time scale. The initial state of the ocean is uniform stratification and complete rest. Although the final temperature distribution is more zonal than it should be, the major western boundary currents and the equatorial undercurrent are successfully predicted. The calculated salinity field has the correct observed range, and correctly indicates that the Atlantic is saltier than the Pacific. It also predicts that the surface waters of the North Pacific are less saline than the surface waters of the South Pacific in accordance with observations. The pack ice model predicts heavy ice i...

Isoneutral Diffusion in a z-Coordinate Ocean Model

This paper considers the requirements that must be satisfied in order to provide a stable and physically based isoneutral tracer diffusion scheme in a z-coordinate ocean model. Two properties are emphasized: 1) downgradient orientation of the diffusive fluxes along the neutral directions and 2) zero isoneutral diffusive flux of locally referenced potential density. It is shown that the Cox diffusion scheme does not respect either of these properties, which provides an explanation for the necessity to add a nontrivial background horizontal diffusion to that scheme. A new isoneutral diffusion scheme is proposed that aims to satisfy the stated properties and is found to require no horizontal background diffusion.

Transient Response in a Z-Level Ocean Model That Resolves Topography with Partial Cells

Ocean simulations are in part determined by topographic waves with speeds and spatial scales dependent on bottom slope. By their very nature, discrete z-level ocean models have problems accurately representing bottom topography when slopes are less than the grid cell aspect ratio Dz/Dx. In such regions, the dispersion relation for topographic waves is inaccurate. However, bottom topography can be accurately represented in discrete zlevel models by allowing bottom-most grid cells to be partially filled with land. Consequently, gently sloping bottom topography is resolved on the scale of horizontal grid resolution and the dispersion relation for topographic waves is accurately approximated. In contrast to the standard approach using full cells, partial cells imply that all grid points within a vertical level are not necessarily at the same depth and problems arise with pressure gradient errors and the spurious diapycnal diffusion. However, both problems have been effectively dealt with. Differences in flow fields between simulations with full cells and partial cells can be significant, and simulations with partial cells are more robust than with full cells. Partial cells provide a superior representation of topographic waves when compared to the standard method employing full cells.

Tracer Conservation with an Explicit Free Surface Method for z-Coordinate Ocean Models

This paper details a free surface method using an explicit time stepping scheme for use in z-coordinate ocean models. One key property that makes the method especially suitable for climate simulations is its very stable numerical time stepping scheme, which allows for the use of a long density time step, as commonly employed with coarse-resolution rigid-lid models. Additionally, the effects of the undulating free surface height are directly incorporated into the baroclinic momentum and tracer equations. The novel issues related to local and global tracer conservation when allowing for the top cell to undulate are the focus of this work. The method presented here is quasi-conservative locally and globally of tracer when the baroclinic and tracer time steps are equal. Important issues relevant for using this method in regional as well as large-scale climate models are discussed and illustrated, and examples of scaling achieved on parallel computers provided.

The Thermal Structure of the Upper Ocean

Abstract The salient feature of the oceanic thermal structure is a remarkably shallow thermocline, especially in the Tropics and subtropics. What factors determine its depth? Theories for the deep thermohaline circulation provide an answer that depends on oceanic diffusivity, but they deny the surface winds an explicit role. Theories for the shallow ventilated thermocline take into account the influence of the wind explicitly, but only if the thermal structure in the absence of any winds, the thermal structure along the eastern boundary, is given. To complete and marry the existing theories for the oceanic thermal structure, this paper invokes the constraint of a balanced heat budget for the ocean. The oceanic heat gain occurs primarily in the upwelling zones of the Tropics and subtropics and depends strongly on oceanic conditions, specifically the depth of the thermocline. The heat gain is large when the thermocline is shallow but is small when the thermocline is deep. The constraint of a balanced heat b...

The Effect of Salinity on the Wind-Driven Circulation and the Thermal Structure of the Upper Ocean

Studies of the effect of a freshening of the surface waters in high latitudes on the oceanic circulation have thus far focused almost entirely on the deep thermohaline circulation and its poleward heat transport. Here it is demonstrated, by means of an idealized general circulation model, that a similar freshening can also affect the shallow, wind-driven circulation of the ventilated thermocline and its heat transport from regions of gain (mainly in the upwelling zones of low latitudes) to regions of loss in higher latitudes. A freshening that decreases the surface density gradient between low and high latitudes reduces this poleward heat transport, thus forcing the ocean to gain less heat in order to maintain a balanced heat budget. The result is a deepening of the equatorial thermocline. (The deeper the thermocline in equatorial upwelling zones is, the less heat the ocean gains.) For a sufficiently strong freshwater forcing, the poleward heat transport all but vanishes, and permanently warm conditions prevail in the Tropics. The approach to warm oceanic conditions is shown to introduce a bifurcation mechanism for the north‐south asymmetry of the thermal and salinity structure of the upper ocean.