Mixing and general circulation dynamics : theory and observations

Submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution February 1987This thesis studies the role of cross-isopycnal mixing in general circulation dynamics, from both the theoretical and observational points of view. The first two chapters discuss some theoretical aspects of cross-isopycnal mixing in the oceans. In chapter one, an integral constraint relating the interior stratification and air-sea heat fluxes is derived, based on the condition that the total mass of water of given density is constant in a steady state ocean. Two simple models are then used to examine the way the numerically small mixing, together with air-sea fluxes, determines the average vertical density stratification of the oceans, and the deep buoyancy driven circulation. In chapter two, a more complete model of a deep flow driven by cross isopycnal diffusion is presented, motivated by the Mediterranean outflow into the North Atlantic. Mixing in this model is responsible for the determination of the detailed structure of the flow and density field, while in the models of the first chapter it was allowed to determine only the average vertical density stratification. In chapter three, a hydrographic data set from the Mediterranean sea is analyzed by inverse methods. The purpose is to examine the importance of mixing when trying to explain tracer distributions in the ocean. The time-mean circulation and the appropriate mixing coefficients are calculated from the hydrographic data. We conclude that the numerically small cross isopycnal mixing processes are crucial to the dynamics, yet difficult to parameterize and measure using available hydrographic data.NSF grants OCE-8521685 and OCE-8017791 supported me during my studies in the joint program

Factors Affecting ENSO’s Period

Accurately capturing the observed mean period of ENSO in general circulation models (GCMs) is often challenging, and it is therefore useful to understand which parameters and processes affect this period. A computationally efficient simulation-based approach is used to extract both the dominant eigenvalues and corresponding eigenvectors of the linearized model from the Zebiak–Cane intermediate-complexity model of ENSO without having to directly construct the linearization. The sensitivity of the period to a variety of parameters is examined, including atmosphere–ocean coupling, atmospheric heating parameterization, thermocline depth zonal profile, western boundary reflection coefficient, atmospheric and ocean wave speeds or Rossby radii of deformation, ocean decay time, and the strength of the annual cycle. In addition to the sensitivity information, the spatial structures of the main fields (SST, thermocline thickness, and more) that are involved in period changes are obtained to aid in the physical interpretation of the sensitivities. There are three main time lags that together compose one-half of a model ENSO period: the Rossby-plus-Kelvin wave propagation time for a wind-caused central Pacific disturbance to propagate to the western ocean and back, SST dynamics that determine the lag between eastern ocean thermocline anomalies and eastern ocean SST anomalies, and the “accumulation” lag of integrating a sufficient delayed wave signal arriving from the western ocean to cancel the eastern ocean anomalies. For any of the parameter changes considered, the eigenvector changes show that the largest contributor to the period change is from changes to the last of these three mechanisms. Physical mechanisms that affect this accumulation delay are discussed, and the case is made that any significant change to ENSO’s period is in turn likely to involve changes to this delay

Applying engineering feedback analysis tools to climate dynamics

The application of feedback analysis tools from engineering control theory to problems in climate dynamics is discussed through two examples. First, the feedback coupling between the thermohaline circulation and wind-driven circulation in the North Atlantic Ocean is analyzed with a relatively simple model, in order to better understand the coupled system dynamics. The simulation behavior is compared with analysis using root locus (in the linear regime) and describing functions (to predict limit cycle amplitude). The second example does not directly involve feedback, but rather uses simulation-based identification of low-order dynamics to understand parameter sensitivity in a model of El Nino/Southern Oscillation dynamics. The eigenvalue and eigenvector sensitivity can be used both to better understand physics and to tune more complex models. Finally, additional applications are discussed where control tools may be relevant to understand existing feedbacks in the climate system, or even to introduce new ones

Two-Way Feedback Interaction between the Thermohaline and Wind-Driven Circulations

The thermohaline circulation (THC) affects the meridional atmospheric temperature gradient and therefore the atmospheric wind and the wind-driven ocean circulation. The wind-driven circulation (WDC), in turn, affects the THC by the advection of salinity anomalies into deep-water formation sites. This paper considers this two-way coupling between the WDC and THC using a simple box-type model and analysis tools from engineering feedback control. The two-way feedback can have a significant effect on the dynamics of the coupled system. For a reasonable choice of parameters, the feedback destabilizes the THC equilibrium for low freshwater forcing. For higher freshwater forcing, the feedback results in a new stable equilibrium instead of the large amplitude oscillation that develops without feedback. It is expected that the analysis approach used here may be broadly applicable to the study of feedback interconnections of other climate systems as well

Methods of testing parameterizations: Vertical ocean mixing

The ocean's velocity field is characterized by an exceptional variety of scales. While the small-scale oceanic turbulence responsible for the vertical mixing in the ocean is of scales a few centimeters and smaller, the oceanic general circulation is characterized by horizontal scales of thousands of kilometers. In oceanic general circulation models that are typically run today, the vertical structure of the ocean is represented by a few tens of discrete grid points. Such models cannot explicitly model the small-scale mixing processes, and must, therefore, find ways to parameterize them in terms of the larger-scale fields. Finding a parameterization that is both reliable and plausible to use in ocean models is not a simple task. Vertical mixing in the ocean is the combined result of many complex processes, and, in fact, mixing is one of the less known and less understood aspects of the oceanic circulation. In present models of the oceanic circulation, the many complex processes responsible for vertical mixing are often parameterized in an oversimplified manner. Yet, finding an adequate parameterization of vertical ocean mixing is crucial to the successful application of ocean models to climate studies. The results of general circulation models for quantities that are of particular interest to climate studies, such as the meridional heat flux carried by the ocean, are quite sensitive to the strength of the vertical mixing. We try to examine the difficulties in choosing an appropriate vertical mixing parameterization, and the methods that are available for validating different parameterizations by comparing model results to oceanographic data. First, some of the physical processes responsible for vertically mixing the ocean are briefly mentioned, and some possible approaches to the parameterization of these processes in oceanographic general circulation models are described in the following section. We then discuss the role of the vertical mixing in the physics of the large-scale ocean circulation, and examine methods of validating mixing parameterizations using large-scale ocean models

Rain driven by receding ice sheets as a cause of past climate change

The Younger Dryas cold period, which interrupted the transition from the last ice age to modern conditions in Greenland, is one of the most dramatic incidents of abrupt climate change reconstructed from paleoclimate proxy records. Changes in the Atlantic Ocean overturning circulation in response to freshwater fluxes from melting ice are frequently invoked to explain this and other past climate changes. Here we propose an alternative mechanism in which the receding glacial ice sheets cause the atmospheric circulation to enter a regime with greater net precipitation in the North Atlantic region. This leads to a significant reduction in ocean overturning circulation, causing an increase in sea ice extent and hence colder temperatures. Positive feedbacks associated with sea ice amplify the cooling. We support the proposed mechanism with the results of a state-of-the-art global climate model. Our results suggest that the atmospheric precipitation response to receding glacial ice sheets could have contributed to the Younger Dryas cooling, as well as to other past climate changes involving the ocean overturning circulation

Spatiotemporal Dynamics of Ice Streams Due to a Triple-Valued Sliding Law

We show that a triple-valued sliding law can be heuristically motivated by the transverse spatial structure of an ice-stream velocity field using a simple one-dimensional model. We then demonstrate that such a sliding law can lead to some interesting stream-like patterns and time-oscillatory solutions. We find a generation of rapid stream-like solutions within a slow ice-sheet flow, separated by narrow internal boundary layers (shear margins), and analyse numerical simulations in two horizontal dimensions over a homogeneous bed and including longitudinal shear stresses. Different qualitative behaviours are obtained by changing a single physical parameter, a mass source magnitude, leading to changes from a slow creeping flow to a relaxation oscillation of the stream pattern, and to steady ice-stream-like solution. We show that the adjustment of the ice-flow shear margins to changes in the driving stress in the one-dimensional approximation is governed by a form of the Ginzburg–Landau equation and use stability analysis to understand this adjustment. In the model analysed here, the width scale of the stream is not set spontaneously by the ice flow dynamics, but rather, it is related to the mass source intensity and spatial distribution.Earth and Planetary Science

Interaction and Variability of Ice Streams under a Triple-Valued Sliding Law and Non-Newtonian Rheology

Ice streams are regions of fast flowing glacier ice that transport a significant portion of the total ice flux from present ice sheets. The flow pattern of ice streams can vary both temporally and spatially. In particular, ice streams can become stagnant and change their path. We study the dynamics of ice streams using an idealized model of an isothermal and power law viscous ice flow that includes horizontal (lateral) shear stresses. The basal sliding law is assumed to be triple-valued. We investigate the spatiotemporal patterns formed because of the flow over a flat bed, fed from an upstream mass source. The ice flows from the mass source region through one or two gaps in a prescribed upstream topographic ridge which restricts the flow, leading to the formation of one or two ice streams. We find a relation between the parameters of the ice rheology and the width of the ice stream shear margins and show how these parameters can affect the minimum width of an ice stream. We also find that complex asymmetric spatiotemporal patterns can result from the interaction of two ice streams sharing a common mass source. The rich spatiotemporal variability is found to mostly be a result of the triple-valued sliding law, but non-Newtonian effects are found to play a significant role in setting a more realistic shear margin width and allowing for relevant time scales of the variability.Earth and Planetary Science