The role of ocean thermodynamics and dynamics in Asian summer monsoon changes during the mid-Holocene

Studies of climate change 6,000 years before present using atmospheric general circulation models (AGCMs) suggest the enhancement and northward shift of the summer Asian and African monsoons in the Northern Hemisphere. Although enhancement of the African monsoonal precipitation by ocean coupling is a common and robust feature, contradictions exist between analyses of the role of the ocean in the strength of the Asian monsoon. We investigated the role of the ocean in the Asian monsoon and sought to clarify which oceanic mechanisms played an important role using three ocean coupling schemes: MIROC, an atmosphere–ocean coupled general circulation model [C]; an AGCM extracted from MIROC coupled with a mixed-layer ocean model [M]; and the same AGCM, but with prescribed sea surface temperatures [A]. The effect of “ocean dynamics” is quantified through differences between experiments [C] and [M]. The effect of “ocean thermodynamics” is quantified through differences between experiments [M] and [A]. The precipitation change for the African and Asian monsoon area suggested that the ocean thermodynamics played an important role. In particular, the enhancement of the Asian monsoonal precipitation was most vigorous in the AGCM simulations, but mitigated in early summer in ocean coupled cases, which were not significantly different from each other. The ocean feedbacks were not significant for the precipitation change in late summer. On the other hand, in Africa, ocean thermodynamics contributed to the further enhancement of the precipitation from spring to autumn, and the ocean dynamics had a modest impact in enhancing precipitation in late summer.     

Stability dependent parameterisations of vertical mixing in ocean general circulation models

Summary Parameterisations of mixing induced through shear instability, internal wave breaking, and double diffusion are investigated in simulations of ocean climate using a global ocean general circulation model (OGCM). Focus is placed on the sensitivity of the large scale circulation, water mass formation and transport of heat as measures of the model's ability to represent current climate. The model resolution is typical of OGCMs being coupled to atmospheric. GCMs in climate models and the parameterisations investigated are all computationally inexpensive enough to allow for integrations on long time scales. Under the assumption of constant vertical eddy coefficients (the control case), the model climatology displays acceptable values of North Atlantic Deep Water formation, Antarctic Circumpolar Current (ACC) transport, and Indonesian through-flow but an excessively deep and diffuse pycnocline structure with weak stratification in the deep ocean. It is found that various circulation and water mass properties are sensitive to the choice of parameterisation of vertical mixing and that determining a scheme which works satisfactorily over all regions (tropical, mid-latitude, and polar) of the domain is not straightforward. Parameterisations of internal wave breaking or upper ocean shear instability lead to some improvements in the model water mass formation. ACC and poleward heat transport when compared to the control case whereas parameterisations of double diffusive processes did not. Based on these and other results, various recommendations are made for mixing parameterisations in ocean climate models.With 8 Figures     

On the validity of vertical dispersion coefficients derived from a Markov-chain scheme for use in gaussian diffusion models

This paper addresses the following question: how do the σ _z values derived from vertical concentration distributions computed by a Markov-chain diffusion model compare with the Σ values which must be used in a gaussian diffusion model to give the same ground level concentration distribution as the one computed by the sophisticated model?     

Principles for capturing the upstream effects of deep sills in low resolution ocean models

Principles for incorporating the upstream effects of deep sills into numerical ocean circulation models using nonlinear analytical hydraulic models are discussed within the context of reduced gravity flow. A method is developed allowing the upstream influence of a numerically unresolvable deep sill or width contraction to be reproduced. The method consists of placing an artificial boundary in the numerical model's overflowing layer at some distance upstream of the actual sill or width contraction of the deep strait. Given the model state at time t, the dependent flow variables are then predicted at the artificial boundary at time t + Δt by using the method of characteristics in combination with quasi-steady hydraulic laws. The calculation requires the use of Riemann invariants and examples are given for a simple nonrotating flow and for rotating channel flow with uniform potential vorticity. The computation is considerably simplified by linearizing the relevant equations in the vicinity of the artificial boundary, resulting in a linear wave reflection problem. The reflection coefficients for the two cases are calculated and these can be used directly to numerically satisfy the boundary condition in a straightforward way.     

On the effects of vertical eddy viscosity on rossby and eady waves in the ocean

The effects of vertical eddy viscosity on simple mesoscale waves in the ocean are studied. The decay of Rossby waves is investigated by one-dimensional depth-dependent linear stability problems which are derived for the interior non-viscous or viscous quasigeostrophic flow using parameterizations of the top and bottom boundary layers corresponding to Ekman suction, no-stress and bottom-stress boundary conditions.The non-slip condition at the bottom yielding an O(E_v^1/2)-Ekman layer causes very short damping times for the 0th Rossby mode. This suggests that this boundary condition is not suitable for mesoscale wave studies, because a Rossby wave fit for the MODE eddy can be done satisfactorily without any damping. Reasonable results for damping times of Rossby waves are obtained by prescribing the bottom stress, resulting from the constant-stress layer at the bottom, and the free-slip condition at the surface. The growth rates of Eady waves are reexamined using this bottom-stress condition.Vertical viscosity in the interior of the ocean, e.g. internal wave induced viscosity, may have a significant influence on the dynamics of the mesoscale motions, comparable to that of the boundary layers in some cases. The results are compatible with the sparse observations available.     

On the influence of frictional parameterization in wind-driven ocean circulation models

In a series of numerical experiments the wind-driven ocean circulation is studied in an idealized, rectangular model ocean, which is forced by steady zonal winds and damped by lateral and/or bottom friction. The problem as described by the barotropic vorticity equation is characterized by a Rossby number (R) and horizontal and/or vertical Ekman numbers (E_L, E_B) only.With free-slip conditions at the boundaries steady solutions for all chosen values of R are obtained, provided the diffusivity is sufficiently large. For both the forms of frictional parameterization a northern boundary current emerges with an eastward penetration scale depending on R. The recirculation pattern in the oceanically relevant ‘intermediate’ range of R is strongly affected by the type of friction. If lateral diffusion dominates bottom friction, a strong recirculating sub-gyre emerges in the northwestern corner of the basin. Its shape resembles the vertically integrated transport fields in recent eddy resolving model (EGCM) studies. The maximum transport is increased to values several times larger than the Sverdrup transport. The increase in transport is coupled with a development of closed contours of potential vorticity, enabling a nearly free inertial flow.This behaviour provides a sharp contrast to the bottom friction case (Veronis) where inertial recirculation only takes place with values of R so large that the eastward jet reaches the eastern boundary. It is shown that the linear friction law puts a strong constraint on the flow by preventing an intense recirculation in a small part of the basin.A reduction of the diffusivity (E_L) in the lateral friction case leads to quasi-steady solutions. The interaction with eddies becomes an integral part of the time mean energetics but does not influence the recirculation character of the flow.The main conclusion of the study is that the horizontal structure of the EGCM-transport fields can be explained in terms of a steady barotropic model where lateral friction represents the dominant dissipation mechanism.     

On the incompatibility of ocean and atmosphere models and the need for flux adjustments

The surface heat and freshwater fluxes from equilibrium ocean (OGCM) and atmospheric (AGCM) general circulation model climates are examined in order to determine the minimum flux adjustment required to prevent climate drift upon coupling. This is accomplished by integrating an OGCM with specified surface fluxes. It is shown that a dramatic climate drift of the coupled system is inevitable unless ocean meridional heat and freshwater (salt) transports are used as constraints for tuning the AGCM present-day climatology. It is further shown that the magnitude of the mismatch between OGCM and AGCM fluxes is not as important for climate drift as the difference in OGCM and implied AGCM meridional heat and freshwater (salt) transports. Hence a minimum flux adjustment is proposed, which is zonally-uniform in each basin and of small magnitude compared to present flux adjustments. This minimum flux adjustment acts only to correct the AGCM implied oceanic meridional transports of heat and freshwater (salt). A slight extension is also proposed to overcome the drift in the surface waters when the minimum flux adjustment is used. Finally, it is suggested that the flux adjustments which arise from current methods used to determine them are all very similar, leading to adjustment fields which are significantly larger than both AGCM and climatological fields over large regions.     

On vertical interpolation and truncation in connexion with use of sigma system models

Abstract

The Mackenzie Shelf in the Canadian Beaufort Sea receives large amounts of freshwater runoff in winter and, yet, it also produces ventilating water masses by brine rejection from growing ice. We examine physical and chemical data to see how these contradictory processes can occur juxtaposed on the shelf. Measurements of salinity and δ¹⁸O both from ice cores and the water column are used to infer the separation into two convective regimes due to the under‐ice topography of the system of large pressure ridges that forms at the boundary between landfast ice and pack ice. Outside this ridge system the ice cover is subject to frequent openings due to offshore ice motion. The inner regime is thus dominated by the impoundment of Mackenzie River water, whereas the outer regime is subject to brine enhancement. This paper compares freezing processes and system evolution for these two regimes in winter.     

Global ocean circulation and equator-pole heat transport as a function of ocean GCM resolution

To determine whether resolution of smaller scales is necessary to simulate large-scale ocean climate correctly, I examine results from a global ocean GCM run with different horizontal grid spacings. The horizontal grid spacings span a range from coarse resolutions traditionally used in climate modeling to nearly the highest resolution attained with today's computers. The experiments include four cases employing 4°, 2°, 1° and 1/2° spacing in latitude and longitude, which were run with minimal differences among them, i.e., in a controlled experiment. Two additional cases, 1/2° spacing with a more scale-selective sub-gridscale mixing of heat and momentum, and approximate 1/2° spacing, are also included. The 1/2° run resolves most of the observed mesoscale eddy energy in the ocean. Artificial constraints on the model tend to minimize differences among the different resolution cases. Nevertheless, the simulations show significant changes as resolution increases. These changes generally but not always bring the model into better agreement with observations. Differences are typically more noticeable when comparing the 4° and 2° runs than when comparing the 2° and 1° runs or the 1° and 1/2° runs. A reasonable conclusion to draw for current studies with coupled ocean-atmosphere GCMs is that the ocean grid spacing could be set to about 1° to accrue the benefits of enhanced resolution without paying an excessively steep price in computer-time cost. The model's poleward heat transport at 1/2° grid spacing peaks at about 1 × 10¹⁵ W in the Northern Hemisphere and 0.5 × 10¹⁵ W in the Southern Hemisphere. These values are significantly below observations, a problem typical of ocean GCMs even when they are less constrained than in the present study. This present problem is alleviated somewhat in the 1/2° run. In this case, however, the eddies resolved by the model generally act to counter rather than to reinforce the heat transport of the mean flow. Improved heat transport may result less from enhanced resolution than from other changes made in this version of the model, such as more accurate wind forcing.     

On the realization of a system of reference in four dimensions for ocean dynamics

The instantaneous position of the sea surface from satellites can be related in principle to an absolute frame of reference in four dimensions, consisting of three space and one time coordinate. A basis therefore exists for referring variations in sea-surface topography, irrespective of Earth-space location, to such a system of reference when used in dynamic equations which apply at the air/sea interface. Such data can be evaluated in both deep oceans and on continental shelves without recourse to subjective judgements on the characteristics of ocean dynamics. In practice, the realization of such a system of reference is only as good as the tracking data used with geodetic concepts.A major problem in the determination of sea-surface topography as a global field is likely to be the establishment of tracking coverage with adequate precision in high latitudes. Surface gravity data could be used to improve orbits in regions where no direct tracking is available but only at the 1–2 m level, due to the data being sampled in relation to the sea surface and not the geoid. It will not be possible to determine sea-surface topography with wavelengths less than about 10³ km in the absence of measurements of higher derivatives of the geopotential with adequate precision in the regions of interest.Asvariations in the sea-surface topography which are not of tidal origin have smaller magnitudes than the sea-surface topography itself, it follows that a proportionate improvement in the tracking resolution will have to precede the unambiguous detection of such variations without restriction of wavelength. Four-dimensional concepts are also essential for the resolution of the tidal signal from satellite altimetry data.The application of principles underlying four-dimensional geodesy to ocean dynamics in the context of a cohesive satellite altimetry program therefore constitutes the necessary basis for synoptic monitoring of the air/sea interface.     

Acoustic thermometry data compared with two ocean models: the importance of Rossby waves and ENSO in modifying the ocean interior

An interpretation is made of interannual changes in acoustic travel time between Oahu and seven receivers at distances of 3000–4000 km. Measurements were made in late 1983, and over two 5-month intervals between 1987 and 1989. Previous publications demonstrated that these changes stem from variations in temperature. Two hydrodynamic ocean models are used to identify plausible oceanic features that could cause these variations. They are from the Naval Research Laboratory and the Florida State University at (1/8)° and (1/6)° resolution, respectively, and are forced with different interannual wind sets for more than a decade. Modelled El Niño's and La Niña's generate poleward travelling Kelvin waves on the eastern boundary of the Pacific. These excite Rossby waves that propagate westward at mid-latitudes. Rossby waves are the dominant model features which affect the modelled acoustic travel times, and hence section-averaged temperatures in the eastern North Pacific. These waves yield travel times whose standard deviations and rates of changes are similar to the measurements. In the observations, some sections separated by less than 500 km exhibit trends in heat content with opposite signs. Similar variability can be explained with modelled Rossby waves. Model wavelengths less than 500 km, eddies, and seasonal cycles induced by seasonal winds yield travel times that are two orders of magnitude too small to account for the data.     

On stochastic stability of regional ocean models to finite-amplitude perturbations of initial conditions

We consider error propagation near an unstable equilibrium state (classified as an unstable focus) for spatially uncorrelated and correlated finite-amplitude initial perturbations using short- (up to several weeks) and intermediate (up to 2 months) range forecast ensembles produced by a barotropic regional ocean model. An ensemble of initial perturbations is generated by the Latin Hypercube design strategy, and its optimal size is estimated through the Kullback–Liebler distance (the relative entropy). Although the ocean model is simple, the prediction error (PE) demonstrates non-trivial behavior similar to that existing in 3D ocean circulation models. In particular, in the limit of zero horizontal viscosity, the PE at first decays with time for all scales due to dissipation caused by non-linear bottom friction, and then grows faster than (quasi)-exponentially. Statistics of a prediction time scale (the irreversible predictability time (IPT)) quickly depart from Gaussian (the linear predictability regime) and becomes Weibullian (the non-linear predictability regime) as amplitude of initial perturbations grows. A transition from linear to non-linear predictability is clearly detected by the specific behavior of IPT variance. A new analytical formula for the model predictability horizon is introduced and applied to estimate the limit of predictability for the ocean model.     

On the bogusing of tropical cyclones in numerical models: The influence of vertical structure

Summary In this study, idealised conditions are used to study the influence of vertical structure of the bogus vortex on its motion in numerical models by comparing the resultant forecast tracks. Two vortices were used: one has a cyclonic circulation throughout the troposphere and the other has an upper tropospheric anticyclone. Both vortices have the same structure in the middle and lower troposphere. The two vortices were inserted into four different environmental flows on a beta-plane: (a) a resting atmosphere; (b) a uniform flow; (c) a horozontal shear flow and (d) a vertical shear flow. The results show that the forecast tracks are very sensitive to the vertical structure of the bogus vortex, especially when the environmental flow is very weak, or is westerly and has a cyclonic horizontal shear. However, this sensitivity is reduced in moderate vertical shear. This motion sensitivity is found to arise from the vertical coupling mechanism by which the upper-and lower-level circulations interact with each other when a horizontal displacement occurs between them.The vertical structure of the bogus vortex can also affect the intensity of the model cyclone, depending on the configuration of the environmental flow. In general, the bogus vortex without an upper-level anticyclone will intensify quicker and will develop more intense than the one with an upper-level anticyclone. The vertical coupling mechanism can result in different asymmetric rainfall pattern in cyclone core region depending on the vertical structure of the bogus vortex. The asymmetric divergent flow associated with these convective asymmetries may in turn further influence the vortex motion. It is suggested that care needs to be taken in determining the vertical structure of the bogus vortex in numerical models.With 14 Figures     

On the cool skin of the ocean

Previous data relating sea-surface temperature to heat flux across the air-sea interface were reanalyzed with a common formula for the wind-stress coefficient. An expression is proposed for the nondimensional thickness of the thermal sublayer: the expression increases with wind velocity at light winds and has a value of 7 when the wind velocity reaches 7 m s^–1.     

Effectiveness of the Bjerknes stability index in representing ocean dynamics

The El Niño-Southern Oscillation (ENSO) is a naturally occurring coupled phenomenon originating in the tropical Pacific Ocean that relies on ocean–atmosphere feedbacks. The Bjerknes stability index (BJ index), derived from the mixed-layer heat budget, aims to quantify the ENSO feedback process in order to explore the linear stability properties of ENSO. More recently, the BJ index has been used for model intercomparisons, particularly for the CMIP3 and CMIP5 models. This study investigates the effectiveness of the BJ index in representing the key ENSO ocean feedbacks—namely the thermocline, zonal advective, and Ekman feedbacks—by evaluating the amplitudes and phases of the BJ index terms against the corresponding heat budget terms from which they were derived. The output from Australian Community Climate and Earth System Simulator Ocean Model (a global ocean/sea ice flux-forced model) is used to calculate the heat budget in the equatorial Pacific. Through the model evaluation process, the robustness of the BJ index terms are tested. We find that the BJ index overestimates the relative importance of the thermocline feedback to the zonal advective feedback when compared with the corresponding terms from the heat budget equation. The assumption of linearity between variables in the BJ index formulation is the primary reason for these differences. Our results imply that a model intercomparison relying on the BJ index to explain ENSO behavior is not necessarily an accurate quantification of dynamical differences between models that are inherently nonlinear. For these reasons, the BJ index may not fully explain underpinning changes in ENSO under global warming scenarios.     

Effects of atmospheric dynamics and ocean resolution on bi-stability of the thermohaline circulation examined using the Grid ENabled Integrated Earth system modelling (GENIE) framework

We have used the Grid ENabled Integrated Earth system modelling (GENIE) framework to undertake a systematic search for bi-stability of the ocean thermohaline circulation (THC) for different surface grids and resolutions of 3-D ocean (GOLDSTEIN) under a 3-D dynamical atmosphere model (IGCM). A total of 407,000 years were simulated over a three month period using Grid computing. We find bi-stability of the THC despite significant, quasi-periodic variability in its strength driven by variability in the dynamical atmosphere. The position and width of the hysteresis loop depends on the choice of surface grid (longitude-latitude or equal area), but is less sensitive to changes in ocean resolution. For the same ocean resolution, the region of bi-stability is broader with the IGCM than with a simple energy-moisture balance atmosphere model (EMBM). Feedbacks involving both ocean and atmospheric dynamics are found to promote THC bi-stability. THC switch-off leads to increased import of freshwater at the southern boundary of the Atlantic associated with meridional overturning circulation. This is counteracted by decreased freshwater import associated with gyre and diffusive transports. However, these are localised such that the density gradient between North and South is reduced tending to maintain the THC off state. THC switch-off can also generate net atmospheric freshwater input to the Atlantic that tends to maintain the off state. The ocean feedbacks are present in all resolutions, across most of the bi-stable region, whereas the atmosphere feedback is strongest in the longitude–latitude grid and around the transition where the THC off state is disappearing. Here the net oceanic freshwater import due to the overturning mode weakens, promoting THC switch-on, but the atmosphere counteracts this by increasing net freshwater input. This increases the extent of THC bi-stability in this version of the model. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.