Sensitivity of the thermohaline circulation and climate to ocean exchanges in a simple coupled model

 A new simple, coupled climate model is presented and used to investigate the sensitivity of the thermohaline circulation and climate to ocean vertical and horizontal exchange. As formulated, the model highlights the role of thin, ocean surface layers in the communication between the atmosphere and the subsurface ocean. Model vertical exchange is considered to be an analogue to small-scale, diapycnal mixing and convection (when present) in the ocean. Model horizontal exchange is considered to be an analogue to the effects of the wind-driven circulation. For small vertical exchange in the ocean, the model exhibits only one steady-state solution: a relatively cold, mid-high-latitude climate associated with a weak, salinity-driven circulation (“off ” mode). For large vertical and horizontal exchange in the ocean, the model also exhibits only one steady-state solution: a relatively warm, mid-high-latitude climate associated with a strong, thermally-driven circulation (“on” mode). For sufficiently weak horizontal exchange but large enough vertical exchange, both modes are possible stable, steady-state solutions. When model parameters are calibrated to fit tracer distributions of the modern ocean-atmosphere system, only the “on” mode is possible in this standard case. This suggests that the wind-driven circulation in consort with diapycnal mixing suppresses the “off ” mode in the modern ocean-atmosphere system. Since both diapycnal mixing and the wind-driven circulation would be expected to increase in a cold climate with greater meridional temperature gradients and enhanced winds, vertical and horizontal exchange in the ocean are probably associated with strong negative feedbacks which tend to stabilize climate. These results point to the need to resolve ocean wind-driven circulation and to greatly improve the treatment of ocean diapycnal mixing in more complete models of the climate system. Received: 16 November 1999 / Accepted: 19 June 2000     

Influence of bottom topography roughness on the jet and inertial recirculation of a mid-latitude gyre

Several numerical experiments are conducted to examine the influence of mesoscale, bottom topography roughness on the inertial circulation of a wind-driven, mid-latitude ocean gyre. The ocean model is based on the quasi-geostrophic formulation, and is eddy-resolving as it features high vertical and horizontal resolutions (six layers and a 10 km grid). An antisymmetrical double-gyre wind stress curl forces the baroclinic modes and generates a strong surface jet. In the case of a flat bottom, inertia and inverse energy cascade force the barotropic mode, and the resulting circulation features strong, barotropic, inertial gyres. The sea-floor roughness inhibits the inertial circulation in the deep layers; the barotropic component of the flow is then forced by eddy-topography interactions, and its energy concentrates at the scales of the topography. As a result, the baroclinicity of the flow is intesified: the barotropic mode is reduced with regard to the baroclinic modes, and the bottom flow (constrained by the mesoscale sea-floor roughness) is decoupled from the surface flow (forced by the gyre-scale wind). Rectified, mesoscale bottom circulation induces an interfacial form stress at the thermocline, which enhances horizontal shear instability and opposes the eastward penetration of the jet. The mean jet is consequently shortened, but the instantaneous jet remains very turbulent, with meanders of large meridional extent. The sea-floor roughness modifies the energy pathways, and the eddies have an even more important role in the establishment of the mean circulation: below the thermocline, rectification processes are dominant, and eddies transfer energy toward permanent mesoscale circulations strongly correlated with topography, whereas above the thermocline mean flow and eddy generation are influenced by the mean bottom circulation through interfacial stress. The topography modifies the vorticity of the barotropic and highest baroclinic modes. Vorticity accumulates at the small topographic scales, and the vorticity content of the highest modes, which is very weak in the flat-bottom case, increases significantly. Few changes occur in surface-intensified modes. In the deep layers of the model, the inverse correlation between relative vorticity and topography at small scales ensures the homogenization of the potential vorticity, which mainly retains the largest scales of the bottom flow and the scale of β.     

The impact of a downslope water-transport parametrization in a global ocean general circulation model

The impact of a downslope water-transport parametrization on the circulation and water mass characteristics of a global depth-level ocean general circulation model is investigated. The spreading of dense water from the formation regions into the deep ocean is known to be poorly represented in depth-level models with no bottom boundary layer resolved or attached. The new scheme is simple and intends to parametrize the effects of various oceanographic processes (rather than the processes themselves) that help dense water to descend topographic slopes by which the formation regions are separated from the world ocean. The new scheme significantly improves the large scale properties of the North Atlantic Deep Water. Changes in the North Atlantic circulation, however, are rather small. In the Southern Ocean, the exchange between the dense water formation regions on the continental shelves and the deep ocean is strengthened at the expense of deep water mass formation by open ocean convection. In all three ocean basins, the density of the deep and bottom water is higher with the new parametrization, which brings the simulations closer to observations in the Atlantic and Indian Oceans. In the Pacific Ocean, however, where the density has already been well reproduced without the downslope transport, it becomes slightly too high. The results are in agreement with those from other model studies.     

Impact of Surface Exchange Coefficients and Sea Surface Cooling on Tropical Cyclone Simulation

The surface flux exchange associated with the exchange coefficients and upper ocean conditions is essential to the development of tropical cyclones (TCs). Using the Weather Research and Forecasting (WRF) model, the present study has investigated the impact of exchange coefficients and ocean coupling during Super Typhoon Saomai (2006). Firstly, two experiments with different formula of roughness are conducted. The experiment with the Donelan formula for drag coefficient (Cd) and ramped formula for enthalpy coefficient (Ck) can simulate stronger intensity compared to other experiments due to the increased surface wind and enthalpy fluxes. That is because the new formulas allows for a smaller Cd and larger Ck in the high wind regime than the former formulas did. Moreover, two coupled simulations between WRF and a one-dimensional ocean model are conducted to examine the feedback of sea surface cooling to the TC. In the experiments with a horizontal uniform mixed layer depth of 70 m, the sea surface cooling is too weak to change the evolution of TC. While in the experiment with an input mixed layer calculated using the Hybrid Coordinate Ocean Model (HYCOM) data, the significant sea surface cooling induces obvious impact on TC intensity and structure. Under the negative feedback of sea surface cooling, the sensible and latent heat fluxes decreases, especially in the right part of Saomai (2006). The negative feedback with coupled ocean model plays a vital role in simulating the intensity and structure of TC.     

On the bottom relief effect on wind-driven currents in a homogeneous ocean

The three-dimensional model of stationary wind-driven currents in a homogeneous ocean of a variable depth is investigated. The model is linear but includes horizontal and vertical turbulent mixings. Two cases of the behaviour of the isolines of the function $\text{?/H}$ are considered, namely: (1) all isolines $\text{?/H}$ start at one part of the coastline and end in another part of it, and (2) a certain isoline $\text{?/H}$ exists which is tangential to the coastline. Here ? is the Coriolis parameter, and H is the depth of the ocean. The first case is the simplest one; it arises in particular if H = constant and the coasts are meridional. The second case is marked by the boundary current separation from the coast. The paper deals with the boundary layers which arise at the surface, bottom, side boundary and inside the ocean.     

Learning about the ocean carbon cycle from observational constraints and model simulations of multiple tracers

A key question in studies of the potential for reducing uncertainty in climate change projections is how additional observations may be used to constrain models. We examine the case of ocean carbon cycle models. The reliability of ocean models in projecting oceanic CO₂ uptake is fundamentally dependent on their skills in simulating ocean circulation and air–sea gas exchange. In this study we demonstrate how a model simulation of multiple tracers and utilization of a variety of observational data help us to obtain additional information about the parameterization of ocean circulation and air–sea gas exchange, relative to approaches that use only a single tracer. The benefit of using multiple tracers is based on the fact that individual tracer holds unique information with regard to ocean mixing, circulation, and air–sea gas exchange. In a previous modeling study, we have shown that the simulation of radiocarbon enables us to identify the importance of parameterizing sub-grid scale ocean mixing processes in terms of diffusive mixing along constant density surface (isopycnal mixing) and the inclusion of the effect of mesoscale eddies. In this study we show that the simulation of phosphate, a major macronutrient in the ocean, helps us to detect a weak isopycnal mixing in the upper ocean that does not show up in the radiocarbon simulation. We also show that the simulation of chlorofluorocarbons (CFCs) reveals excessive upwelling in the Southern Ocean, which is also not apparent in radiocarbon simulations. Furthermore, the updated ocean inventory data of man-made radiocarbon produced by nuclear tests (bomb ¹⁴C) enable us to recalibrate the rate of air–sea gas exchange. The progressive modifications made in the model based on the simulation of additional tracers and utilization of updated observational data overall improve the model’s ability to simulate ocean circulation and air–sea gas exchange, particularly in the Southern Ocean, and has great consequence for projected CO₂ uptake. Simulated global ocean uptake of anthropogenic CO₂ from pre-industrial time to the present day by both previous and updated models are within the range of observational-based estimates, but with substantial regional difference, especially in the Southern Ocean. By year 2100, the updated model estimated CO₂ uptake are 531 and 133 PgC (1PgC?=?10¹⁵ gram carbon) for the global and Southern Ocean respectively, whereas the previous version model estimated values are 540 and 190 PgC.     

A mechanism of deep-ocean mixing due to near-inertial waves generated by flow over bottom topography

This article describes a remarkable effect of the horizontal component of the Coriolis vector, which is traditionally neglected in geophysical fluid wave theory, for the generation of inertio-gravity waves caused by flow over topography. Earlier an initial-value approach was used to calculate the evolutions of internal waves generated by atmospheric forcing in the linear Boussinesq model as done by many investigators for the study of near-inertial waves in the ocean. In this study, however, we focused on the mechanism of generating near-inertial waves by a bottom forcing with various periods. It is shown that the non-traditional Coriolis effects can significantly enhance the generation of near-inertial waves. Moreover, the intensity of generated near-inertial waves increases as the period of bottom forcing approaches to the local inertial period, suggesting the importance of tidal flows for generating near-inertial waves in deep ocean near critical latitudes. It is inferred that such near-inertial waves may contribute to a mechanism of deep-ocean mixing.     

The role of lateral ocean physics in the upper ocean thermal balance of a coupled ocean-atmosphere GCM

 The sensitivity of the upper ocean thermal balance of an ocean-atmosphere coupled GCM to lateral ocean physics is assessed. Three 40-year simulations are performed using horizontal mixing, isopycnal mixing, and isopycnal mixing plus eddy induced advection. The thermal adjustment of the coupled system is quite different between the simulations, confirming the major role of ocean mixing on the heat balance of climate. The initial adjustment phase of the upper ocean (SST) is used to diagnose the physical mechanisms involved in each parametrisation. When the lateral ocean physics is modified, significant changes of SST are seen, mainly in the southern ocean. A heat budget of the annual mixed layer (defined as the “bowl”) shows that these changes are due to a modified heat transfer between the bowl and the ocean interior. This modified heat intake of the ocean interior is directly due to the modified lateral ocean physics. In isopycnal diffusion, this heat exchange, especially marked at mid-latitudes, is both due to an increased effective surface of diffusion and to the sign of the isopycnal gradients of temperature at the base of the bowl. As this gradient is proportional to the isopycnal gradient of salinity, this confirms the strong role of salinity in the thermal balance of the coupled system. The eddy induced advection also leads to increased exchanges between the bowl and the ocean interior. This is both due to the shape of the bowl and again to the existence of a salinity structure. The lateral ocean physics is shown to be a significant contributor to the exchanges between the diabatic and the adiabatic parts of the ocean. Received: 24 January 2000 / Accepted: 11 September 2000     

The axial angular momentum balance of a global ocean general circulation model

In this study we examine the axial angular momentum balance of a non-eddy-resolving global ocean general circulation model, from the perspective of the geographical and seasonal variability of angular momentum and from the perspective of the torques acting on the ocean through its surfaces. Our purpose is to provide an estimate of the magnitude of the seasonal storage of angular momentum in the ocean and hence the oceanic excitation of variability in length of day, and to elucidate the role of the ocean in transferring angular momentum between the atmosphere and the Earth's crust. We provide an assessment of the reliability of the model results by examining the sensitivity of the angular momentum and torque distributions to several model parameters.Although the Southern Ocean region containing the Antarctic Circumpolar Current (ACC) makes the largest contribution to both the annual mean oceanic angular momentum and its seasonal variability, inclusion of the rest of the world ocean reduces both of these quantities to about two-thirds of the value of the Southern Ocean alone. The annual, global mean angular momentum is found to be insensitive to most model choices except for the isopycnal diffusivity. The seasonal variability, on the other hand, is insensitive to the isopycnal diffusivity, but sensitive to the smoothness of the representation of topography and moderately sensitive to horizontal and vertical friction parameterizations. The torque balance at all latitudes, including within the Antarctic circumpolar belt, is between wind stress and bottom pressure torques. Horizontal friction torques are small but non-negligible. Bottom friction and storage of angular momentum are negligible in angular momentum budgets on seasonal time scales. Two commonly used wind stress climatologies, one based on historical marine meteorological observations and the other based on operational weather analyses, differ in the sign of the globally integrated wind stress torque.     

Steering of upper ocean currents and fronts by the topographically constrained abyssal circulation

A two-layer theory is used to investigate (1) the steering of upper ocean current pathways by topographically constrained abyssal currents that do not impinge on the bottom topography and (2) its application to upper ocean – topographic coupling via flow instabilities where topographically constrained eddy-driven deep mean flows in turn steer the mean pathways of upper ocean currents and associated fronts. In earlier studies the two-layer theory was applied to ocean models with low vertical resolution (2–6 layers). Here we investigate its relevance to complex ocean general circulation models (OGCMs) with high vertical resolution that are designed to simulate a wide range of ocean processes. The theory can be easily applied to models ranging from idealized to complex OGCMs, provided it is valid for the application. It can also be used in understanding some persistent features seen in observed ocean frontal pathways (over deep water) derived from satellite imagery and other data. To facilitate its application, a more thorough explanation of the theory is presented that emphasizes its range of validity. Three regions of the world ocean are used to investigate its application to eddy-resolving ocean models with high vertical resolution, including one where an assumption of the two-layer theory is violated. Results from the OGCMs with high vertical resolution are compared to those from models with low vertical resolution and to observations. In the Kuroshio region upper ocean – topographic coupling via flow instabilities and a modest seamount complex are used to explain the observed northward mean meander east of Japan where the Kuroshio separates from the coast. The Japan/East Sea (JES) is used to demonstrate the impact of upper ocean – topographic coupling in a relatively weak flow regime. East of South Island, New Zealand, the Southland Current is an observed western boundary current that flows in a direction counter to the demands of Sverdrup flow and counter to the direction simulated in nonlinear global flat bottom and reduced gravity models. A model with high vertical resolution (and topography extending through any number of layers) and a model with low vertical resolution (and vertically compressed but otherwise realistic topography confined to the lowest layer) both simulate a Southland Current in the observed direction with dynamics depending on the configuration of the regional seafloor. However, the dynamics of these simulations are very different because the Campbell Plateau and Chatham Rise east and southeast of New Zealand are rare features of the world ocean where the topography intrudes into the stratified water column over a relatively broad area but lies deeper than the nominal 200 m depth of the continental shelf break, violating a limitation of the two-layer theory. Observations confirm the results from the high vertical resolution model. Overall, the model simulations show increasingly widespread upper ocean – topographic coupling via flow instabilities as the horizontal resolution of the ocean models is increased, but fine resolution of mesoscale variability and the associated flow instabilities are required to obtain sufficient coupling. As a result, this type of coupling is critical in distinguishing between eddy-resolving and eddy-permitting ocean models in regions where it occurs.     

An upper ocean general circulation model for climate studies: global simulation with seasonal cycle

The global ocean circulation with a seasonal cycle has been simulated with a two-and-a-half layer upper-ocean model. This model was developed for the purpose of coupling to an atmospheric general circulation model for climate studies on decadal time scales. The horizontal resolution is 4° latitude by 5° longitude and is thus not eddy-resolving. Effects of bottom topography are neglected. In the vertical, the model resolves the oceanic mixed layer and the thermocline. A thermodynamic sea-ice model is coupled to the mixed layer. The model is forced at the surface with seasonally varying (a) observed wind stress, (b) heat fluxes, as defined by an atmospheric equilibrium temperature, and (c) Newtonian-type surface salt fluxes. The second layer is coupled to the underlying deep ocean through Newtonian-type diffusive heat and salt fluxes, convective overturning, and mass entrainment in the upwelling regions of the subpolar gyres. The overall global distributions of mixed layer temperature, salinity and thickness are favorably reproduced. Inherent limitations due to coarse horizontal resolution result in large mixed-layer temperature errors near continental boundaries and in weak current systems. Sea ice distributions agree well with observations except in the interiors of the Ross and Weddell Seas. A realistic time rate of change of heat storage is simulated. There is also realistic heat transport from low to high latitudes.     

Time-dependent topographic meandering of a baroclinic current

The effects of baroclinic instability of a broad ocean current, flowing in an ocean basin with a plane sloping bottom, on the path of the current are studied. The set of equations governing this path and its variation with depth are the vorticity equation and the heat equation. It is assumed that the vertical and horizontal temperature contrasts are comparable as suggested by observations of the Gulf Stream. When quasi-geostrophy is assumed in addition, this implies that the leading contribution to the heat equation does not contain the vertical advection of the basic stratification. This corresponds to the long-wave approximation of the usual baroclinic-instability problem. The heat equation determines the vertical variation of the path and when this is combined with the vorticity equation, the equation governing the path at one level is obtained. The path equation requires a specification of the direction and curvature at the inlet and these conditions are taken to be time-dependent. When these conditions contain frequencies for which the current is unstable, meanders in the path of the current increase in amplitude downstream of the inlet. When the path at the inlet changes suddenly from one parallel to the isobaths to one making a small angle with them, the region of instability in which the amplitude of the meanders increases, is confined to a restricted segment of the path, at soms distance from the inlet. This region becomes advected with the basic current, and its extent increases with time. The amplitude of the meanders in this region increases while their wavelength decreases in time because the shorter waves are unstabler. The increase in amplitude and decrease in wavelength in a restricted segment of the path could lead to eddy formation in a finite-amplitude model and may therefore suggest a mechanism for eddy formation in the Gulf Stream.     

Climate Change on Newfoundland and Labrador Shelves: Results From a Regional Downscaled Ocean and Sea-Ice Model Under an A1B Forcing Scenario 2011–2069

Climate change may affect ocean and ice conditions in coastal oceans and thus have significant impacts on coastal infrastructure, marine navigation, and marine ecosystems. In this study a three-dimensional ice–ocean model is developed to examine likely changes of ocean and ice conditions over the Newfoundland and Labrador Shelves in response to climate change. The model is configured with a horizontal grid of approximately 7?km and a vertical grid of 46 levels and is run from 1979 to 2069. The projection period is 2011 to 2069 under a median emission scenario A1B used by the Intergovernmental Panel on Climate Change. For the projection period, the surface atmospheric forcing fields used are from the Canadian Regional Climate Model over the North Atlantic. The open boundary conditions come from the Canadian Global Climate Model, Version 3 (CGCM3), adjusted for the 1981–2010 mean of the Simple Ocean Data Assimilation model output. The simulated fields over the 1981–2010 period have patterns consistent with observations. Over the Newfoundland and Labrador Shelves during the projection period, the model shows general trends of warming, freshening, and decreasing ice. From 2011 to 2069, the model projects that under A1B sea surface temperature will increase by 1.4°C; bottom temperature will increase by 1.6°C; sea surface salinity will decrease by 0.7; bottom salinity will decrease by 0.3; and sea-ice extent will decrease by 70%. The sea level will rise by 0.11?m at the St. John's tide-gauge station because of oceanographic change, and the freshwater transport of the Labrador Current will double as a result of freshening. The regional ice–ocean model reproduces more realistic present climate conditions and projects considerably different future climate conditions than CGCM3.     

A numerical world ocean general circulation model

This paper describes a numerical model of the world ocean based on the fully primitive equations. A “Standard” ocean state is introduced into the equations of the model and the perturbed thermodynamic variables are used in the modle’s calculations. Both a free upper surface and a bottom topography are included in the model and a sigma coordinate is used to normalize the model’s vertical component. The model has four unevenly-spaced layers and 4 × 5 horizontal resolution based on C-grid system. The finite-difference scheme of the model is designed to conserve the gross available energy in order to avoid fictitious energy generation or decay.The model has been tested in response to the annual mean surface wind stress, sea level air pressure and sea level air temperature as a preliminary step to its further improvement and its coupling with a global atmospheric general circulation model. Some of results, including currents, temperature and sea surface elevation simulated by the model are presented.     

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.     

The formation and maintenance of the North Water Polynya

Abstract

This paper investigates the formation and maintenance of the North Water Polynya, Baffin Bay in winter using a multi‐category sea‐ice model coupled with the Princeton ocean model. Monthly climatological atmospheric data from the National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) reanalysis provides the forcing. An objectively‐analysed climatology provides the initial ocean temperature and salinity. Wind stress drives the ice in a cyclonic gyre around northern Baffin Bay. Localized regions of thin ice form where wind drives ice away from coastlines or fast ice. The regions of thin ice are characterized by enhanced ice growth, exceeding 1.2 m mo^?1. In the regions of thin ice, surface ocean heat flux is also enhanced and is between 30–60 W m^?2. Surface heat flux is, in part, attributable to convective mixing and entrainment driven by ice growth. The surface heat flux reflects advection of the warm West Greenland Current. Heat and salt balances show that horizontal advective exchange counterbalances surface fluxes of heat and salt.     

A study on the nonlinear stability of fronts in the ocean on a sloping continental shelf

1.IntroductionArnol'd(1965,1969)variationalprincipleandapriorestimatemethodisessentiallyageneralizationofLyapunovstabilitymethodforfinite--dimensionaldynamicalsystemsininfinite--dimensionalones,andhestudiedthenonlinearstabilityof2--dimensionalincompressibleidealfluidmotionbyuseofthismethod,andestablishedtwotheoremswhichareArnol,d'sfirsttheoremandArnol'd'ssecondtheorem.Eversincethe1980's,manyscientistshavebeenworkingonthissubject,Holmetal.(1985);MclntyreandShepherd(1987);Zeng(1989);Muetal.(1…     

Estimation of the rate of water renewal in the ocean conveyor

The imitating numerical model of the turbulent exchange between the ocean conveyor belt and the surrounding waters is proposed. The water particle exchange through the belt surface of the given size during its movement along the deep conveyor branch is considered. The length of the turbulent jump is parameterized with the aid of coefficients of vertical and horizontal turbulent diffusion and its direction is simulated using the Monte Carlo method. By means of the computation of trajectories of 10⁶ separate particles, their spatial distribution is determined for different time intervals after the deep water formation. It is demonstrated that due to the lateral exchange, the water in the given part of the belt is renewed by 60% during 25 years and by 75% during 75 years. A comparison with the computations based on the usual method is carried out for the vertical exchange: it is revealed that in both cases the considerable renewal of water in the conveyor belt occurs during the time period of about 10⁴ years.     

Modeling the tropical Pacific Ocean using a regional coupled climate model

A high-resolution tropical Pacific general circulation model (GCM) coupled to a global atmospheric GCM is described in this paper. The atmosphere component is the 5°×4°global general circulation model of the Institute of Atmospheric Physics (IAP) with 9 levels in the vertical direction. The ocean component with a horizontal resolution of 0.5°, is based on a low-resolution model (2°×1°in longitude-latitude).Simulations of the ocean component are first compared with its previous version. Results show that the enhanced ocean horizontal resolution allows an improved ocean state to be simulated; this involves (1) an apparent decrease in errors in the tropical Pacific cold tongue region, which exists in many ocean models,(2) more realistic large-scale flows, and (3) an improved ability to simulate the interannual variability and a reduced root mean square error (RMSE) in a long time integration. In coupling these component models, a monthly "linear-regression" method is employed to correct the model's exchanged flux between the sea and the atmosphere. A 100-year integration conducted with the coupled GCM (CGCM) shows the effectiveness of such a method in reducing climate drift. Results from years 70 to 100 are described.The model produces a reasonably realistic annual cycle of equatorial SST. The large SSTA is confined to the eastern equatorial Pacific with little propagation. Irregular warm and cold events alternate with a broad spectrum of periods between 24 and 50 months, which is very realistic. But the simulated variability is weaker than the observed and is also asymmetric in the sense of the amplitude of the warm and cold events.     

Combined wavelet and principal component analysis (WEof) of a scale-oriented model of coastal ocean gravity waves

Coastal ocean numerical modeling is basically the representation of the dynamics of the coastal ocean in a chosen range of length scales and over an associated frequency band, including the modeling of both coherent processes and associated transient processes. The ocean dynamical features can be individually identified by combining wavelet analysis for time and frequency localization and principal component analysis to “decorrelate” physically consistent structures. In the present paper, the so-called WEof analysis is applied for the extraction of external gravity waves and internal gravity wave lower modes in a simple case of a flat bottom, constant Brunt-Väisälä ocean. It is shown that, with some well known restrictive assumptions, WEof analysis is an efficient candidate for the recognition of frequency localized dynamical processes.