On the use of current meter data to assess the realism of ocean model simulations

The evaluation of ocean simulations against observed datasets is essential to assess their realism and to guide model development, but often remains qualitative, and ignores certain datasets. This paper presents a three-dimensional, quantitative comparison of a 1/6° Atlantic numerical simulation (CLIPPER) with the WOCE current meter dataset in terms of mean velocity and eddy kinetic energy. Our metrics reveal the good behaviour of CLIPPER open boundary conditions and forcing with respect to full-depth current records. Due to its still moderate resolution, however, the model globally underestimates the observed mean speeds and eddy activity. This discrepancy is barely noticeable at low latitudes but increases toward the poles, probably since the poleward decrease of the Rossby radius exceeds that of the horizontal grid step. At least in this eddy-admitting regime, it is suggested that the numerics of geopotential-coordinate models like ours dissipate mean and eddy momentum at depth and adversely affect current–topography interactions.     

A global ocean mesh to overcome the North Pole singularity

A semi-analytical method is presented for constructing a global orthogonal curvilinear ocean mesh which has no singularity point inside the computational domain since the mesh poles are moved to land points. The method involves defining an analytical set of mesh parallels in the stereographic polar plan, computing the associated set of mesh meridians, and projecting the resulting mesh onto the sphere. The set of mesh parallels proposed here is defined as a series of embedded circles. The resulting mesh presents no loss of continuity in either the mesh lines or the scale factors over the whole ocean domain, as the mesh is not a composite mesh. Thus, the Bering Strait can be opened without specific treatment. The equator is a mesh line, which provides a better numerical solution for equatorial dynamics. The resolution can be easily controlled through the definition of three analytical functions which can increase resolution and/or maintain a low ratio of anisotropy. The mesh has been implemented in the LODYC general circulation ocean model. Results of a semi-diagnostic simulation are shown.     

On the climate response of the low-latitude Pacific Ocean to changes in the global freshwater cycle

Under global warming, the predicted intensification of the global freshwater cycle will modify the net freshwater flux at the ocean surface. Since the freshwater flux maintains ocean salinity structures, changes to the density-driven ocean circulation are likely. A modified ocean circulation could further alter the climate, potentially allowing rapid changes, as seen in the past. The relevant feedback mechanisms and timescales are poorly understood in detail, however, especially at low latitudes where the effects of salinity are relatively subtle. In an attempt to resolve some of these outstanding issues, we present an investigation of the climate response of the low-latitude Pacific region to changes in freshwater forcing. Initiated from the present-day thermohaline structure, a control run of a coupled ocean–atmosphere general circulation model is compared with a perturbation run in which the net freshwater flux is prescribed to be zero over the ocean. Such an extreme experiment helps to elucidate the general adjustment mechanisms and their timescales. The atmospheric greenhouse gas concentrations are held constant, and we restrict our attention to the adjustment of the upper 1,000 m of the Pacific Ocean between 40°N and 40°S, over 100 years. In the perturbation run, changes to the surface buoyancy, near-surface vertical mixing and mixed-layer depth are established within 1 year. Subsequently, relative to the control run, the surface of the low-latitude Pacific Ocean in the perturbation run warms by an average of 0.6°C, and the interior cools by up to 1.1°C, after a few decades. This vertical re-arrangement of the ocean heat content is shown to be achieved by a gradual shutdown of the heat flux due to isopycnal (i.e. along surfaces of constant density) mixing, the vertical component of which is downwards at low latitudes. This heat transfer depends crucially upon the existence of density-compensating temperature and salinity gradients on isopycnal surfaces. The timescale of the thermal changes in the perturbation run is therefore set by the timescale for the decay of isopycnal salinity gradients in response to the eliminated freshwater forcing, which we demonstrate to be around 10–20 years. Such isopycnal heat flux changes may play a role in the response of the low-latitude climate to a future accelerated freshwater cycle. Specifically, the mechanism appears to represent a weak negative sea surface temperature feedback, which we speculate might partially shield from view the anthropogenically-forced global warming signal at low latitudes. Furthermore, since the surface freshwater flux is shown to play a role in determining the ocean’s thermal structure, it follows that evaporation and/or precipitation biases in general circulation models are likely to cause sea surface temperature biases.     

A degradation approach to accelerate simulations to steady-state in a 3-D tracer transport model of the global ocean

 We have developed a new method to accelerate tracer simulations to steady-state in a 3-D global ocean model, run off-line. Using this technique, our simulations for natural ¹⁴C ran 17 times faster when compared to those made with the standard non-accelerated approach. For maximum acceleration we wish to initialize the model with tracer fields that are as close as possible to the final equilibrium solution. Our initial tracer fields were derived by judiciously constructing a much faster, lower-resolution (degraded), off-line model from advective and turbulent fields predicted from the parent on-line model, an ocean general circulation model (OGCM). No on-line version of the degraded model exists; it is based entirely on results from the parent OGCM. Degradation was made horizontally over sets of four adjacent grid-cell squares for each vertical layer of the parent model. However, final resolution did not suffer because as a second step, after allowing the degraded model to reach equilibrium, we used its tracer output to re-initialize the parent model (at the original resolution). After re-initialization, the parent model must then be integrated only to a few hundred years before reaching equilibrium. To validate our degradation-integration technique (DEGINT), we compared ¹⁴C results from runs with and without this approach. Differences are less than 10‰ throughout 98.5% of the ocean volume. Predicted natural ¹⁴C appears reasonable over most of the ocean. In the Atlantic, modeled Δ¹⁴C indicates that as observed, the North Atlantic Deep Water (NADW) fills the deep North Atlantic, and Antartic Intermediate Water (AAIW) infiltrates northward; conversely, simulated Antarctic Bottom Water (AABW) does not penetrate northward beyond the equator as it should. In the Pacific, in surface eastern equatorial waters, the model produces a north–south assymetry similar to that observed; other global ocean models do not, because their resolution is inadequate to resolve equatorial dynamics properly, particularly the intense equatorial undercurrent. The model’s oldest water in the deep Pacific (at −239‰) is close to that observed (−248‰), but is too deep. Surface waters in the Southern Ocean are too rich in natural ¹⁴C due to inadequacies in the OGCM’s thermohaline forcing. Received: 18 March 1997 / Accepted: 27 July 1997     

Impact of westerly wind bursts on the warm pool of the TOGA-COARE domain in an OGCM

 A primitive equation model is used to investigate the warm pool equilibrium of the tropical Pacific ocean. Attention is focused on the upper ocean. The oceanic response is described using an isothermal approach applied to warm waters contained in the TOGA-COARE domain. The heat balance shows that all the terms, atmospheric surface fluxes, advection and diffusion, operate in the heat bugdet with different time scales. Over long periods, diffusive heat fluxes transfer heat received from the atmosphere out of the warm pool trough the top of the main thermocline. Over short periods, the impact of westerly wind bursts modifies this balance: atmospheric heating is reversed, diffusion is enhanced and advective heat transports out of the warm pool operate through zonal and vertical contributions. We were able to relate the two latter processes to zonal jets and Ekman pumping, respectively. Conversely, the meridional contribution always represents a source of heat, mainly due to the tropical wind convergence. The modelling results clearly show that except during strong wind events, entrainment cooling is not an important component of the budget. The inability to remove heat is due to the salt stratification which needs to be first reduced or even destroyed by westerly wind bursts to activate heat entrainment into deeper layers. Finally, we suggest that the near zero estimate for the surface heat flux entering the warm pool may be extended to longer periods including seaosnal to interannual time scale. Received: 16 December 1996 / Accepted: 8 July 1997     

The role of the intra-daily SST variability in the Indian monsoon variability and monsoon-ENSO–IOD relationships in a global coupled model

The impact of diurnal SST coupling and vertical oceanic resolution on the simulation of the Indian Summer Monsoon (ISM) and its relationships with El Ni?o-Southern Oscillation (ENSO) and Indian Ocean Dipole (IOD) events are studied through the analysis of four integrations of a high resolution Coupled General Circulation Model (CGCM), but with different configurations. The only differences between the four integrations are the frequency of coupling between the ocean and atmosphere for the Sea Surface Temperature (SST) parameter (2 vs. 24?h coupling) and/or the vertical oceanic resolution (31 vs. 301 levels) in the CGCM. Although the summer mean tropical climate is reasonably well captured with all the configurations of the CGCM and is not significantly modified by changing the frequency of SST coupling from once to twelve per day, the ISM–ENSO teleconnections are rather poorly simulated in the two simulations in which SST is exchanged only once per day, independently of the vertical oceanic resolution used in the CGCM. Surprisingly, when 2?h SST coupling is implemented in the CGCM, the ISM–ENSO teleconnection is better simulated, particularly, the complex lead-lag relationships between the two phenomena, in which a weak ISM occurs during the developing phase of an El Ni?o event in the Pacific, are closely resembling the observed ones. Evidence is presented to show that these improvements are related to changes in the characteristics of the model’s El Ni?o which has a more realistic evolution in its developing and decaying phases, a stronger amplitude and a shift to lower frequencies when a 2-hourly SST coupling strategy is implemented without any significant changes in the basic state of the CGCM. As a consequence of these improvements in ENSO variability, the lead relationships between Indo-Pacific SSTs and ISM rainfall resemble the observed patterns more closely, the ISM–ENSO teleconnection is strengthened during boreal summer and ISM rainfall power spectrum is in better agreement with observations. On the other hand, the ISM–IOD teleconnection is sensitive to both SST coupling frequency and the vertical oceanic resolution, but increasing the vertical oceanic resolution is degrading the ISM–IOD teleconnection in the CGCM. These results highlight the need of a proper assessment of both temporal scale interactions and coupling strategies in order to improve current CGCMs. These results, which must be confirmed with other CGCMs, have also important implications for dynamical seasonal prediction systems or climate change projections of the monsoon.     

Towards a parametrization of river discharges into ocean general circulation models: a closure through energy conservation

Diagnostic methods are defined in order to compare two numerical simulations of ocean dynamics in a region of freshwater influence. The first one is a river plume simulation based on a high resolution numerical configuration of the POM coastal ocean model in which mixing parametrizations have been previously defined. The second one is a simulation based on the NEMO Global Ocean Model used for climate simulations in its half-a-degree configuration in which a river inflow is represented as precipitation on two coastal grid cells. Both simulations are forced with the same freshwater inflows and wind stresses. The divergence of volumetric fluxes above and below the halocline are compared. Results show that when an upwelling wind blows, the two models display similar behavior although the impact of lack of precision can be observed in the NEMO configuration. When a downwelling wind blows, the NEMO Global Ocean configuration can not reproduce the coastally trapped baroclinic dynamics because its grid resolution is too coarse. To find a parametrization to help represent these dynamics in ocean general circulation models, a method based on energy conservation is investigated. This method shows that it is possible to link the energy fluxes provided by river inflows to the divergence of energy fluxes integrated over the grid cells of ocean general circulation models. A parametrization of the dynamics created by freshwater inflows is deduced from this method. This enabled creation of a box model that proved to have the same behavior as the fluxes previously computed from the high resolution configuration.     

The roles of surface heat flux and ocean heat transport convergence in determining Atlantic Ocean temperature variability

The temperature variability of the Atlantic Ocean is investigated using an eddy-permitting (1/4°) global ocean model (ORCA-025) forced with historical surface meteorological fields from 1958 to 2001. The simulation of volume-averaged temperature and the vertical structure of the zonally averaged temperature trends are compared with those from observations. In regions with a high number of observations, in particular above a depth of 500 m and between 22° N and 65° N, the model simulation and the dataset are in good agreement. The relative contribution of variability in ocean heat transport (OHT) convergence and net surface heat flux to changes in ocean heat content is investigated with a focus on three regions: the subpolar and subtropical gyres and the tropics. The surface heat flux plays a relatively minor role in year-to-year changes in the subpolar and subtropical regions, but in the tropical North Atlantic, its role is of similar significance to the ocean heat transport convergence. The strongest signal during the study period is a cooling of the subpolar gyre between 1970 and 1990, which subsequently reversed as the mid-latitude OHT convergence transitioned from an anomalously weak to an anomalously strong state. We also explore whether model OHT anomalies can be linked to surface flux anomalies through a Hovmöller analysis of the Atlantic sector. At low latitudes, increased ocean heat gain coincides with anomalously strong northward transport, whereas at mid-high latitudes, reduced ocean heat loss is associated with anomalously weak heat transport.     

Physical processes contributing to the water mass transformation of the Indonesian Throughflow

The properties of the waters that move from the Pacific to the Indian Ocean via passages in the Indonesian archipelago are observed to vary with along-flow-path distance. We study an ocean model of the Indonesian Seas with reference to the observed water property distributions and diagnose the mechanisms and magnitude of the water mass transformations using a thermodynamical methodology. This model includes a key parameterization of mixing due to baroclinic tidal dissipation and simulates realistic water property distributions in all of the seas within the archipelago. A combination of air–sea forcing and mixing is found to significantly change the character of the Indonesian Throughflow (ITF). Around 6 Sv (approximately 1/3 the model net ITF transport) of the flow leaves the Indonesian Seas with reduced density. Mixing transforms both the intermediate depth waters (transforming 4.3 Sv to lighter density) and the surface waters (made denser despite the buoyancy input by air–sea exchange, net transformation?=?2 Sv). The intermediate transformation to lighter waters suggests that the Indonesian transformation contributes significantly to the upwelling of cold water in the global conveyor belt. The mixing induced by the wind is not driving the transformation. In contrast, the baroclinic tides have a major role in this transformation. In particular, they are the only source of energy acting on the thermocline and are responsible for creating the homostad thermocline water, a characteristic of the Indonesian outflow water. Furthermore, they cool the sea surface temperature by between 0.6 and 1.5°C, and thus allow the ocean to absorb more heat from the atmosphere. The additional heat imprints its characteristics into the thermocline. The Indonesian Seas cannot only be seen as a region of water mass transformation (in the sense of only transforming water masses in its interior) but also as a region of water mass formation (as it modifies the heat flux and induced more buoyancy flux). This analysis is complemented with a series of companion numerical experiments using different representations of the mixing and advection schemes. All the different schemes diagnose a lack of significant lateral mixing in the transformation.     

The seasonal cycle in coupled ocean-atmosphere general circulation models

We examine the seasonal cycle of near-surface air temperature simulated by 17 coupled ocean-atmosphere general circulation models participating in the Coupled Model Intercomparison Project (CMIP). Nine of the models use ad hoc “flux adjustment” at the ocean surface to bring model simulations close to observations of the present-day climate. We group flux-adjusted and non-flux-adjusted models separately and examine the behavior of each class. When averaged over all of the flux-adjusted model simulations, near-surface air temperature falls within 2?K of observed values over the oceans. The corresponding average over non-flux-adjusted models shows errors up to ～6?K in extensive ocean areas. Flux adjustments are not directly applied over land, and near-surface land temperature errors are substantial in the average over flux-adjusted models, which systematically underestimates (by ～5?K) temperature in areas of elevated terrain. The corresponding average over non-flux-adjusted models forms a similar error pattern (with somewhat increased amplitude) over land. We use the temperature difference between July and January to measure seasonal cycle amplitude. Zonal means of this quantity from the individual flux-adjusted models form a fairly tight cluster (all within ～30% of the mean) centered on the observed values. The non-flux-adjusted models perform nearly as well at most latitudes. In Southern Ocean mid-latitudes, however, the non-flux-adjusted models overestimate the magnitude of January-minus-July temperature differences by ～5?K due to an overestimate of summer (January) near-surface temperature. This error is common to five of the eight non-flux-adjusted models. Also, over Northern Hemisphere mid-latitude land areas, zonal mean differences between July and January temperatures simulated by the non-flux-adjusted models show a greater spread (positive and negative) about observed values than results from the flux-adjusted models. Elsewhere, differences between the two classes of models are less obvious. At no latitude is the zonal mean difference between averages over the two classes of models greater than the standard deviation over models. The ability of coupled GCMs to simulate a reasonable seasonal cycle is a necessary condition for confidence in their prediction of long-term climatic changes (such as global warming), but it is not a sufficient condition unless the seasonal cycle and long-term changes involve similar climatic processes. To test this possible connection, we compare seasonal cycle amplitude with equilibrium warming under doubled atmospheric carbon dioxide for the models in our data base. A small but positive correlation exists between these two quantities. This result is predicted by a simple conceptual model of the climate system, and it is consistent with other modeling experience, which indicates that the seasonal cycle depends only weakly on climate sensitivity.