Print version ISSN 0102-261X
Rev. Bras. Geof. vol.27 no.2 São Paulo Apr./June 2009
Volume and heat transports in the world oceans from an ocean general circulation model
Luiz Paulo de Freitas AssadI; Audalio Rebelo Torres JuniorII; Wilton Zumpichiatti ArrudaIII; Affonso da Silveira Mascarenhas JuniorIV; Luiz LandauV
IUniversidade Federal do Rio de Janeiro, COPPE, Centro de Tecnologia, Bloco I, sala 214, Laboratório de Métodos Computacionais em Engenharia, Av. Athos da Silveira Ramos, 149, Cidade Universitária, Ilha do Fundão, 21941-909 Rio de Janeiro, RJ, Brazil. Phone: +55 (21) 2562-8419 - E-mail: firstname.lastname@example.org
IIUniversidade Federal do Rio de Janeiro, CCMN, Instituto de Geociências, Departamento de Meteorologia, Laboratório de Modelagem de Processos Marinhos e Atmosféricos, Av. Athos da Silveira Ramos, s/n, Cidade Universitária, Ilha do Fundão, 21949-900 Rio de Janeiro, RJ, Brazil. Phone: +55 (21) 2562-9470 ext. 21- E-mail: email@example.com
IIIUniversidade Federal do Rio de Janeiro, CCMN, Instituto de Matemática, Av. Athos da Silveira Ramos, 149, Bl. C, CT, Cidade Universitária, Ilha do Fundão, 21949-900 Rio de Janeiro, RJ, Brazil. Phone: +55 (21) 2562-7508 ext. 222 - E-mail: firstname.lastname@example.org
IVCentro Internacional para la Investigación del Fenómeno El Niño (CIIFEN), Escobedo #1204 y 9 de Octubre, P.O. Box 09014237, Guayaquil, Ecuador. Phone: +59 (34) 251-4770; Fax: +59 (34) 251-4771 - E-mail: email@example.com
VUniversidade Federal do Rio de Janeiro, COPPE, Centro de Tecnologia, Bloco I, sala 214, Laboratório de Métodos Computacionais em Engenharia, Av. Athos da Silveira Ramos, 149, Cidade Universitária, Ilha do Fundão, 21941-909 Rio de Janeiro, RJ, Brazil. Phone: +55 (21) 2562-8415 - E-mail: firstname.lastname@example.org
O monitoramento dos transportes de volume e calor nos oceanos é de fundamental importância para o estudo do sistema climático, sua variabilidade e possíveis mudanças. A aplicação de Modelos de Circulação Global dos Oceanos em estudos climáticos necessita que os dados dinâmicos e termodinâmicos gerados pelos mesmos se encontrem em equilíbrio. Esse trabalho apresenta alguns resultados obtidos a partir da aplicação do Modular Ocean Model , em sua versão 4.0, inicializado com campos de temperatura, salinidade, velocidade e elevação da superfície livre em equilíbrio, oriundos do experimento de assimilação de dados oceânicos conduzido pelo Geophysical Fluid Dynamics Laboratory (GFDL). O uso desse conjunto de dados objetivou a diminuição do tempo de "spin-up" do modelo oceânico utilizado. O modelo foi integrado por sete anos, ao fim dos quais foram estimados transportes de volume e calor em diferentes seções ao redor do globo. Tais estimativas demonstraram boa concordância com valores estimados em trabalhos anteriores. Os resultados demonstraram a ocorrência de um ciclo sazonal oceânico bem definido a partir do segundo ano de integração do modelo, assim como importantes aspectos da circulação oceânica global como o "conveyor belt", os giros subtropicais oceânicos e o escoamento associado à Corrente Circumpolar Antártica. O experimento realizado demonstrou a importância desse modelo como fonte de dados a ser utilizado como condições iniciais e de contorno em modelos regionais de circulação oceânica, assim como para longas integrações a serem aplicadas em estudos de variabilidade e mudanças climáticas.
Palavras-chave: clima oceânico, modelagem computacional, transporte de volume, transporte de calor.
Monitoring the volume and heat transports around the world oceans is of fundamental importance in the study of the climate system, its variability, and possible changes. The application of an Oceanic General Circulation Model for climatic studies needs that its dynamic and thermodynamic fields are in equilibrium. The time spent by the model to reach this equilibrium is called spin-up time. This work presents some results obtained from the application of the Modular Ocean Model version 4.0 initialized with temperature, salinity, velocity and sea surface height data already in equilibrium from an ocean data assimilation experiment conducted by Geophysical Fluid Dynamics Laboratory (GFDL). The use of this dataset as initial condition aimed to diminish the spin-up time of the oceanic model. The model was integrated for seven years. Volume and heat transports in different sections around the world oceans showed to be in good agreement with the literature. The results showed a well defined seasonal cycle starting at the second integration year, also some important dynamic and thermodynamic aspects of the Global Ocean Circulation, as the great conveyor belt, are well reproduced. The obtained results could constitute an important data source to be used as initial and boundary conditions in regional ocean model experiments as well as for long integration runs in order to study oceanic climate variability.
Keywords: ocean climate, computational modeling, volume transport, heat transport.
Monitoring the volume and heat transports around the worldoceans is of fundamental importance to understand the climate system, its variability and possible changes. In the last years, most of the climate studies focused on the following methods: global ocean numerical modeling (Maltrud & McClean, 2003; Stammer et al., 2003) and hydrographic inverse methods(Schmitz, 1996; Ganachaud & Wunsch, 2000).
According to Peixoto & Oort (1992) the "climate" state could be understood as a set of average quantities that characterizes the structure and the behavior of the atmosphere, hydrosphere and cryosphere for a particular period. The application of an Ocean Global Circulation Model (OGCM) to make climatic forecasts requires the model to be in dynamic and thermodynamic equilibrium.
To start the integration of an OGCM it is necessary to specify initial conditions. The initial fields are basically represented by the following variables: temperature, salinity, velocity field and sea surface height. In general, models are started from a state of rest with climatological distribution of temperature and salinity fields imposed. On this manner the model physics will generate dynamical and thermodynamical fields in equilibrium state.The time needed by the model to reach this equilibrium state is called spin-up time and could range from days for the upper layers to hundred of years for the deep layers of the ocean. Long spin-up times are related with the low velocity propagation of some internal ocean modes. Semtner & Chervin (1992) use a high resolution OGCM to study the general circulation of the oceans. With horizontal resolution of 0.5º× 0.5º and 20 levels in the vertical, the model took 32.5 years to reach its dynamic equilibrium starting from a state of rest with imposed climatological temperature and salinity fields from Levitus (1982). A way to initialize an OGCM, that dramatically reduces the spin-up time, is the use of all initial fields already in thermodynamic and dynamic equilibrium and not only the mass structure. These data can be obtained from another numerical experiment that has already reached a "climate" state. Using temperature and salinity distributions from Navy's Modular Ocean Data Assimilation System 1/8º climatological product (Fox et al., 2002 apud Maltrud & McClean, 2003), Maltrud & McClean (2003) initialized a 1/10º resolution Parallel Ocean Program general circulation model. The modelreached quasi-equilibrium for the mesoscale upper ocean processes after 15 years integration. Using a regional model for the South Atlantic, Matano & Philander (1993) estimate a spin-up time of 3 years. The model was initialized not only with a climatologic density field but also with the vertical structure of the transport for the Antarctic Circumpolar Current (ACC) andAgulhas Current (AC). This work presents the application of this kind of strategy to initialize a coarse resolution OGCM and also some results from 7 years of model integration (Assad, 2006).
A brief description of the model and some characteristics of the applied methodology are given in the next section where also it will be described the grid setup, model parameterizations, forcing, initial and boundary conditions applied in the simulation. The ocean volume and heat transports annual mean values and time series analysis of some dynamic and thermodynamic prognostic transports are presented and discussed.
In this work we use the Modular Ocean Model 4.0 (MOM 4.0) developed by GFDL under the administration of National Oceanic and Atmospheric Administration (NOAA). This z-coordinate model has been largely adopted by the scientific community for global climatic research investigations. It is mathematically formulated with a Navier-Stokes set of equations under the Boussinesq and hydrostatic approximations. The system of continuous equations are basically completed by an equation of state of sea-water, the continuity equation for an incompressible fluid, and conservation equations for temperature and salinity. A complete description of this model can be found in Pacanowsky& Griffies (1999).
The prognostic variables are space distributed following the staggered Arakawa B-grid system. The MOM 4.0 code also allows the construction of spherical global grids with numerical poles not necessarily coincident with the geographical ones. This method minimizes problems related with the convergence of meridians over the ocean that happens mainly on North Pole.The grid constructed on this work uses the tripolar grid method by Murray (1996). It adopts normal spherical coordinates south of 65º N while North of this latitude two numerical poles situated on land regions of Asia and North America are used.
The longitudinal resolution is 1º worldwide while the latitudinal resolution increases from 1º to 1/3º within the 10ºN-10ºS equatorial band. On the vertical, 50 levels are implemented. In order to accommodate a higher resolution near the ocean surface, the first 22 levels are located in the top 220 m depth. The thickness of each level varies from 10 m to a maximum of 366.6 m near the sea bottom (5500 m). Regions shallower than 40 m are not considered. This grid follows the oceanic component of the GFDL's climate model presented on the 4th IPCC Assessment of Global Climate Changes (Griffies et al., 2005). The bathymetry used is distributed by the Southampton Oceanography Centre. It is a composite of different products: remote sensing data from Smith & Sandwell (1997) between 72ºS and 72ºN; the International Bathymetric Chart of the Oceans (Jakobssen et al., 2000; Griffies et al., 2005) North of 72ºN; and ETOPO 5 South of 72ºS.
In this section, we present the parameterizations related with the model turbulence, such as vertical mixing coefficients, neutral physics, and horizontal friction. We also describe some aspects related with the representation of the penetration of solar shortwave radiation into the model.
For the tracer diffusion, it is used a modified formulation of canonical Bryan & Lewis (1979). In the upper ocean, we assume a diffusion coefficient of 0.1 × 10-4 m2s-1 in the tropics and 0.3 × 10-4 m2s-1 elsewhere. For the deeper ocean a uniform value of 1.2 × 10-4 m2s-1 in is used. The neutral physics parameterization was applied by using a depth-independent constant neutral diffusivity value of 600 m2s-1. In MOM 4.0, regions where the neutral slope steepens, such as near the upper ocean boundary layer and in convective regions, neutral diffusivity is exponentially converted to horizontal diffusion. The conversion starts when the neutral slope becomes steeper than 1/500.
The horizontal friction parameterization was applied with an anisotropic viscosity scheme (Large et al., 2001) within the equatorial band from 20ºS to 20ºN. Outside this band, the viscosity reverts to the traditional isotropic scheme with a grid size-dependent and vertically constant background viscosity added to a horizontal shear-dependent Smagorinsky viscosity (Smagorinsky, 1963). For the vertical viscosity, it was chosen a depth independent background viscosity of 10-4 m2s-1. In addition to the background vertical diffusivity and viscosity, we adopt the K-profile parameterization of diapycnal mixing (Large et al., 1994).This parameterization prescribes added levels of tracer and velocity mixing in regions where mixing is under-represented, such as, in surface boundary layer (Griffies et al., 2005).
The absorption of shortwave solar radiation in the upperlayers of the ocean varies significantly in space and time. In regions with low chlorophyll concentrations this penetration could reach 20 to 30 m depth. The high vertical resolution (10 meters in the first 220 meters) in this work makes necessary the correctrepresentation of this process. So, we use a climatologic seasonal chlorophyll field from SeaWIFS sensor together with theoptical model by Morel & Antoine (1994). This model allows a sunlight penetration of 150 m. Below this level, the shortwave penetration is ignored.
In this work we use, as initial condition, the data from the Ocean Data Assimilation for Seasonal to Interannual Prediction experiment (ODASI) conducted by GFDL (Sun et al., 2007). The ODASI also uses MOM 4.0 and generated 40 years of monthly data between 1963 and 2003. Another important characteristic of this dataset is the high vertical spatial resolution which is similar to the one used in this work, especially for the upper layers of the ocean. The month of January 1985 was chosen to be our initial condition since it does not present strong climate anomalies,such as El Niño.
The sea surface boundary conditions are taken from the climatological dataset of Ocean Model Intercomparison Project (OMIP) (Röeske, 2001). This version of OMIP dataset was produced by ECMWF (European Centre for Medium-Range Weather Forecasts) under ERA-15 project. This project generated 15 years of validated data from 1979 to 1993 by applying data assimilationtechniques on numerical experiments. The variables we use as surface forcing are meridional and zonal wind stress components, net shortwave and long-wave radiation, sensible heat flux, specific humidity flux, frozen and liquid precipitation. The OMIP outputs consist of daily means of an annual cycle except for the specific humidity and sensible heat flux which were represented by the monthly means of an annual cycle.
The boundary fields described above are cyclically imposed to the model until it reaches an equilibrium state. The model was integrated for seven years with an external mode time step of 80s and a baroclinic time-step of 4800s. Using 8 processors of an ALTIX 350 computer, the model was able to complete one simulated year in about 1.5 day, generating about 2 GB output per model year.
RESULTS AND DISCUSSIONS
The integrated global ocean kinetic energy shows a well defined cyclic behavior starting at second integration year (Fig. 1). The most energetic periods occur during the austral winter and the less energetic periods during austral summer. This result indicates that the ocean kinetic energy in the annual period is dominated by the wind stress annual energy cycle (Oort et al., 1994).Also it is possible to identify a small time growing tendency ofthe kinetic energy curve. This trend is linearly estimated as 1.82 × 1011 Watts/day, what represents only 0.0003% of the difference between the maximum and minimum global values of ocean kinetic energy. This fact can be attributed to the difference between the OMIP wind stress filed used on this work and the ODASI wind stress forcing field that was used to produce the initial conditions. According to Wunsch (1998) the largest contribution to the global wind-stress field comes from thezonal wind-stress which responds by 96% of the global work done by the wind. The main regions of zonal wind energy input are the Southern Ocean, Kuroshio region, and Gulf Stream/NorthAtlantic Current regions.
Figure 2 shows the zonal mean of the zonal component of the wind stress for OMIP and ODASI wind products and the difference between them. We can observe that the zonal wind component in OMIP is more intense than in ODASI on most latitudes, especially for the high latitude bands of both hemispheres where the most energetic wind-stress field is found. The difference between OMIP and ODASI wind-stress globally averagedin space and time is 6.22 × 1011 Watts.day-1 that is of the same order as the global kinetic energy grow rate (1.82 × 1011 Watts.day-1). It is important to note that the first integration year is not considered on the kinetic energy growth rate estimate toavoid noises associated with the initial warm-up phase.
The cyclic behavior of the globally integrated kinetic energy was not restricted to the upper layers but also can be seen indeeper ocean layers where the kinetic energy peaks are much less intense. On Figure 3, we plot the ocean kinetic energy integrated for each level between 1000 m and 5500 m for the whole South Atlantic basin.
Volume transport analysis
To analyze the model response to the imposed forcing theworld ocean was divided in nine sections (three meridional and six zonal), where volume and heat transports were estimated (Fig. 4). In order to allow a comparison with previous works, we chose a set of sections following closely the ones used by Ganachaud & Wunsch (2000) and Maltrud & McClean (2003). Section transports were integrated in three different vertical layers: 0-1500 m (upper layer), 1500-3500 m (deep layer) and 3500-5500 m (bottom layer) according to the classification by Maltrud& McClean (2003).
Figure 5 displays the annual mean (for the last six years of integration) of the vertically integrated volume transport (intensity and direction) and associated variance across each monitored section (Fig. 4) in three different levels 0-1250 m (upper layer), 1250-3500 m (deep layer) and 3500-5500 m (bottom layer) for sections 1-6, and in the whole water column for sections 7-9.
The zonal sections situated in the Atlantic Ocean (sections 1 and 2) present similar behavior with a northward volume transport for the upper layer and a southward transport for the deep layers. These results are in accordance with the general ocean circulation in Atlantic sector where heat is exported to the Northin the upper layers and to the South by the deep water movement. The obtained values agree with literature estimates. Applying an inverse calculation method on WOCE hydrographic data, Ganachaud & Wunsch (2000) obtain 16 ± 3 Sv (1 Sv = 106 m3.s-1) and 16 ± 2 Sv northward volume transports in the upper layer for sections 1 and 2 respectively. Using an eddy resolving OGCM, Maltrud & McClean (2003) find average volume transports of 15 Sv and 16 Sv in the upper layer of the same sections. For the deep layer we get southward transports of17.64 ± 1.06 Sv and 18.13 ± 0.13 Sv for sections 1 and 2 respectively. Schmitz (1996) find southward transports of 18 Svfor both sections while Maltrud & McClean (2003) find 17 Sv and 19 Sv. In section 1, at bottom layer, we find southward volume transport. This result does not agree with estimates obtained by Schmitz (1996), Ganachaud & Wunsch (2000) andMaltrud & McClean (2003) that found northward transport. This discrepancy is due to the positioning of North Atlantic Deep Water flow at this section whose lower bound is at deeper levels (deeper than 3500 m). Also, the northward flow of Antarctic BottomWater, with upper salinity and temperature values of 34.7 and 0.3ºC in this region (Pickard & Emery, 1990), occupies depths deeper than 3700 m and it is confined to the western side of the basin. If we consider the upper limit of the bottom layer at 3800 m, we estimate a northward volume transport of 1.17 ± 0.41 Sv in this layer, what is comparable with the literature.
In section 2 it is also observed southward transport of 2.87 ± 0.6 Sv. In this section Ganachaud & Wunsch (2000) foundalso southward transport of 4 ± 2 Sv.
The grid used in the present experiment has a free communication between North Atlantic and Arctic Oceans and, consequently, ensures the existence of non-zero resultant volumetransports. In Figure 6, two zonal and one meridional section situated between Europe and Groveland show the existence of inflows at high latitudes of the North Atlantic. The values have the same order of the ones obtained by Stammer et al. (2003)on close sections. These transports could explain non zero net volume transports through the monitored Atlantic Ocean areas.
In section 3, at the upper layer, our model gives a southward transport of 15.15 ± 9.09 Sv. This value agrees with Schmitz(1996) and Maltrud & McClean (2003) that find 18 Sv and 16 Sv, respectively. The northward deep layer transport presents much less intensity (1.13 Sv) and stronger variability (± 3.11) than the bottom layer values. Ganachaud & Wunsch (2000) also find northward volume transport of 3 ± 5 Sv in this layeron a close section.
The zonal Pacific sections (4 and 5) present northward upper layer volume transports of 11.07 ± 30.59 Sv and 6.9 ±16.78 Sv. Section 6 presents southward volume transports for all the monitored layers.
The Indonesian throughflow (Fig. 1, section 10) presents 13.36 ± 1.33 Sv westward volume transport for the wholewater column. Ganachaud & Wunsch (2000), Maltrud & McClean (2003), and Stammer et al. (2003) find for this section:15 ± 5 Sv, 12 Sv, and 11.5 ± 5 Sv, respectively (Fig. 5).
Figure 5 also shows that the most intense transports are associated with the ACC zonal flux. Due to ACC barotropic nature on its meridional sections (Fig. 4; sections 7, 8, and 9) we calculate the volume transports integrated on the whole water. Significant transport variances are found on these sections probably related with strong OMIP zonal wind stress variance on the ACC zonal band.
The time series plots of the vertically integrated eastward volume transport of the ACC at Drake Passage (Fig. 4, section 7) has an annual cyclic behavior with maximum values on austral winter and minimum values on austral summer (Fig. 7).The magnitude of this transport ranges from 137 Sv to 148 Sv. The mean zonal volume transport across the Drake Passage is 143.1 ± 11.1 Sv, which is close to Whitworth & Peterson (1985) estimate of 124.7 ± 9.9 Sv based on multi-years bottom pressure and current meter moorings records. These authors also find the same annual cycle obtained. Based on inverse methods applied to WOCE hydrographic data, Rintoul (1991) find 130 ± 13 Sv, while Ganachaud & Wunsch (2000) find 140 ± 6 Sv.
In order to make a volume balance, the world ocean was divided in six different areas (Fig. 8).
On Table 1 we present the average balance of volume inside each monitored area on Figure 8, what help us to evaluate the model ability to produce results that converge to an ocean climate state (Table 1). These values should be analyzed together with the boundary values related with the hydrological cycle used in the model. In this experiment the hydrological cycle was represented by the following variables (in m.s-1): Precipitation minus evaporation, and runoff. The balance of volume (in m3.s-1) in each region is divided by the respective area in order to allow a comparison with hydrological variables that are expressed in m.s-1. As seen before, the annual mean values related with the balance of volume in each monitored area are insignificant when compared with the hydrological cycle variables and would represent an insignificant sea surface elevation growth.
In Figure 9 it is shown a time series plot of the volume transport balance for each monitored area. It is possible to observe a well defined cyclical behavior in all monitored areas.
Heat transport analysis
Before we start, it is important to emphasize that what we call "heat transport" in this section is, in fact, the heat advection (advective heat transport).
In Figure 10, we plot the mean vertically integrated (from surface to bottom) heat transport and respective variance across each monitored section (Fig. 4). The Atlantic and Pacific zonal sections (Fig. 4; sections 1, 2, 4, 5, and 6) present northward heat transports, while the Indian Ocean zonal section (Section 3) presents a southward heat transport that agrees with other estimates (Ganachaud & Wunsch, 2000; Stammer et al., 2003).
Section 1 presents a northward heat transport of 0.60 ± 0.02 PW (1 PW = 1015 Watts), what reinforces the peculiar behavior of the South Atlantic basin as exporter of heat to theNorth Atlantic. The same behavior is observed across section 2, but with an enhanced northward heat transport that could be associated with the strong upper layer (0-1500m) northward heat and volume transports associated with the Gulf Stream. Northward resultant heat transports were also observed in the works of Ganachaud & Wunsch (2000) and Stammer et al. (2003) atclose sections.
The most intense heat transports are associated with the ACC zonal flux through its meridional sections (Fig. 4; sections 7, 8, and 9). The time series plot for the vertically integrated heat transport perpendicular to each monitored section also reveals a well defined annual cycle behavior (not shown). For the Drake Passage section (Fig. 4; section 7), maximum values are found during austral summer and minimum during austral winter (Fig. 11). The magnitude of this transport ranges from 1.48 to 1.58 PW with a mean value of 1.52 ± 0.0008 PW. These values agree with previous estimate of 1.3 PW by Ganachaud & Wunsch (2000).
Another significant westward heat transport is observed in Indonesia meridional section (Fig. 4; section 10). In this section it is estimated an annual mean value of 1.07 PW ± 0.08 PW. Stammer et al. (2003) find an westward heat transport of 1.12 PW ± 0.52 PW from an ocean general circulation model and Ganachaud & Wunsch (2000) find 1.4 PW westward heat transport bythe application of inverse methods.
Table 2 shows the gain (positive values) or loss (negative values) of heat by each monitored ocean area shown in Figure 8. It reveals a net gain of heat by southern hemisphere ocean surface and a net loss of heat by the northern hemisphere ocean surface.
Figure 12 reveals the presence of an intense northward heat transport between equator and 40ºN and a southward heat transport between the equator and 40ºS what indicates that the conducted experiment was able to reproduce the so-called conveyor belt circulation. Using a geostrophic method on data obtained during the WOCE project, Ganachaud & Wunsch (2000) estimate the mean value of -0.6 ± 0.3 PW for the global meridional heat transport across 30ºS and 2.2 ± 0.6 PW on 12ºN. On the same latitudes we find -0.2 ± 0.01 PW and 1.75 ± 0.03 PW, respectively.
The heat transport analysis also reveals a well defined cycle for this property. Again, it is also possible to see a small growth tendency that could be related with the imposition of the new heat flux boundary condition to the model. The ODASI uses initial conditions generated by a numerical experiment where sea surface temperature data assimilation techniques were applied.As presented before, the heat flux boundary conditions used in this experiment are the climatological OMIP data. In Table 2, it is presented the resultant sea surface heat flux integrated in each monitored area. Positive values are found on SouthernHemisphere, what demonstrate the importance of these oceansas heat reservoirs of the Earth's climate system. In contrast, Northern Hemisphere ocean present negative surface heat fluxes. This fact can be related with the ratio of land/ocean area that is almost three times bigger on Northern Hemisphere thanon the Southern Hemisphere (Peixoto & Oort, 1992). Another important point to emphasize is the high variance values associated with each sea surface heat flux estimate, especially on equatorial regions. This fact could be associated with surface velocities that are strong influenced by wind field fluctuations that are very significant on this region.
The estimated annual mean volume and advective heat transports across the sections in Figure 4 agree with existing in literature. It can be observed some significant variances associated with meridional transports in low latitude regions, especially, in upper levels due to strong wind spatial and temporal variations in these regions. The ACC volume transport at Drake Passage reveals a well defined seasonal cycle starting at the second integration year and followed the same seasonal behavior found in the scientific literature. Also, high ACC volume transport variances are found, what seems to be related with barotropic fluctuations of the OMIP zonal wind-stress variance on that latitude band. Some important physical aspects of the ocean global circulation, as the so-called conveyor-belt meridional heat circulation, could also be observed on model results. The northward advective heat transport found in the 32º South Atlantic section is in accordancewith the scientific literature.
The initialization strategy used on the present work was successful in decreasing the model spin-up time. This work shows that the chosen model configuration and initialization strategy is able to reproduce the main balances in the world oceans, what makes this model suitable for longer integration periods for studying low frequency climate events. The obtained results could constitute an important data source to be used as initial and boundary conditions in regional ocean model experiments, as well as, for long integration runs in order to study oceanic climate variability.
The authors would like to thank the ocean and climate group of GFDL for the strong help in configuring the ocean model and turning available the needed data to initialize the experiment. The third author would also like to thank the National Council forScientific and Technological Development (CNPq) for the support under VARICONF project (Grant #476472/2006-7).
ASSAD LPF. 2006. Influência do Campo de Vento anômalo tipo ENSO na Dinâmica do Atlântico Sul. Doctorate Thesis, Universidade Federal do Rio de Janeiro, Rio de Janeiro, 222 p. [ Links ]
BRYAN K & LEWIS LJ. 1979. A water mass model of the world ocean. Journal of Geophysical Research, 84: 2503-2517. [ Links ]
GANACHAUD A & WUNSCH C. 2000. Improved estimates of global circulation, heat transport and mixing from hydrographic data. Nature, 408: 453-457. [ Links ]
GRIFFIES SM, GNANADESIKAN A, DIXON KW, DUNNE JP, GERDES R, HARRISON MJ, ROSATI A, RUSSELL JL, SAMUELS BL, SPELMAN MJ, WINTON M & ZHANG R. 2005. Formulation of an ocean model for global climate simulations. Ocean Science, 1: 45-79. [ Links ]
JAKOBSSEN M, CHERKIS N, WOODWARD J, MacNAB R & COAKLEY B. 2000. New grid of Arctic bathymetry aids scientists and mapmakers, EOS Trans. Am. Geophys. Union, p. 81, 89, 93, 96. [ Links ]
LARGE WG, McWILLIAMS JC & DONEY SC. 1994. Oceanic vertical mixing: A review and a model with non-local boundary layer parameterization. Rev. Geophys., 32: 363-403. [ Links ]
LARGE WG, DANABASOGLU G, McWILLIAMS JC, GENT PR & BRYAN FO. 2001. Equatorial circulation of a global ocean climate model with anisotropic horizontal viscosity. Journal of Physical Oceanography, 31: 518-536. [ Links ]
LEVITUS S. 1982. Climatological Atlas of the World Ocean. - NOAA Prof. Pap. No. 13. U.S. Government Printing Office, Washington, DC, 173 pp. [ Links ]
MALTRUD ME & McCLEAN JL. 2003. An eddy resolving global 1/10º ocean simulation. Ocean Modelling, 8: 31-54. [ Links ]
MATANO RP & PHILANDER SGH. 1993. Heat and Mass Balances of the South Atlantic Ocean calculated from a Numerical Model. Journal ofGeophysical Research, 98: 977-984. [ Links ]
MOREL A & ANTOINE D. 1994. Heating rate within the upper ocean in relation to its bio-optical state. Journal of Physical Oceanography, 24: 1652-1665. [ Links ]
MURRAY RJ. 1996. Explicit generation of orthogonal grids for ocean models. Journal of Computational Physics, 126: 251-273. [ Links ]
OORT AH, ANDERSON LA & PEIXOTO JP. 1994. Estimates of the energy cycle of the oceans. Journal of Geophysical Research, 99: 7665-7688. [ Links ]
PACANOWSKY RC & GRIFFIES SM. 1999. The MOM3 Manual. Geophysical Fluid Dynamics Laboratory/NOAA, Princeton, USA, p. 680. [ Links ]
PEIXOTO JP & OORT AH. 1992. Physics of Climate. American Institute of Physics. 520 pp. [ Links ]
PICKARD GL & EMERY WJ. 1990. Descriptive Physical Oceanography: An Introduction. Pergamon Press. 320 pp. [ Links ]
RINTOUL SR. 1991. South Atlantic Interbasin Exchange. Journal ofGeophysical Research, 96: 2675-2692. [ Links ]
RÖESKE F. 2001. An Atlas of Surface Fluxes based on the ECMWF Re-Analysis - a Climatological Dataset to force Global Ocean General Circulation Models. Max-Planck Institut für Meteorologie, Hamburg. Report no. 323. ISSN 0937-1060. [ Links ]
SCHMITZ WJ. 1996. On The World Ocean Circulation: Volume II. Technical Report. Woods Hole Oceanographic Institution. WHOI-96-08. 241 p. [ Links ]
SEMTNER AJ & CHERVIN RM. 1992. Ocean General Circulation from Global Eddy-Resolving Model. Journal of Geophysical Research, 97: 5493-5550. [ Links ]
SMAGORINSKY J. 1963. General Circulation experiments with the primitive equations I. The basic experiment. Monthly Weather Review,91(3): 99-164. [ Links ]
SMITH WHF & SANDWELL DT. 1997. Global sea floor topography from satellite altimetry and ship depth soundings. Science, 277: 1956-1961. [ Links ]
STAMMER D, WUNSCH C, GIERING R, ECKERT C, HEIMBACH P, MAROTZKE J, ADCROFT A, HILL CN & MARSHALL J. 2003. Volume, heat, and freshwater transports of the global ocean circulation 1993-2000, estimated from a general circulation model constrained by World Ocean circulation Experiment (WOCE) data. Journal of Geophysical Research, 108: pp. 7-1-7-23. [ Links ]
SUN C, RIENECKER MM, ROSATI A, HARRISON M, WITTENBERG A, KEPPENNE CL, JACOB JP & KOVACH RM. 2007. Comparison and Sensitivity of ODASI Ocean Analyses in the Tropical Pacific. Monthly Weather Review, 135: 2242-2264. [ Links ]
WHITWORTH T III & PETERSON RG. 1985. The volume transport of the Antarctic Circumpolar Current from bottom pressure measurements.Journal of Physical Oceanography, 15: 810-816. [ Links ]
WUNSCH C. 1998. The work done by the wind on the oceanic general circulation. Journal of Physical Oceanography, 28: 2332-2340. [ Links ]
Recebido em 29 setembro, 2008 / Aceito em 9 junho, 2009
Received on September 29, 2008 / Accepted on June 9, 2009
NOTES ABOUT THE AUTHORS
Luiz Paulo de Freitas Assad. Bachelor in Oceanography (UERJ/1996), Master in Physical Oceanography (USP/2000) and Doctor in Civil Engineering in the area of computational modeling applied to the Environmental Engineering (UFRJ/2006). At present, he is researcher at COPPE/UFRJ. His area of interest is Physical Oceanography, with emphasis in ocean computational modeling, acting mainly in the following subjects: global and regional ocean computational modeling, ocean climate variability, and analysis of oceanographic data.
Audalio Rebelo Torres Junior. Bachelor in Meteorology (UFRJ/1982), specialized in Signal Processing Neural Networks (Petrobras/1992), Master in Oceanic Engineering (UFRJ/1995), and Doctor in Oceanic Engineering (UFRJ/2005). At present, he is professor at the Meteorology Department/UFRJ. The Physical Oceanography is his area of interest, working mainly in the following subjects: El Niño, ocean circulation, climate, numerical modeling, teleconnections processes, and oceanography of the South Atlantic.
Wilton Zumpichiatti Arruda. Bachelor in Computer Science (UFRJ/1990), Master and Doctor in Mathematics (UFRJ/1991, 1996), and Ph.D. in Physical Oceanography at Florida State University (2002). He is professor at Mathematical Institute/UFRJ since 1995. His areas of interest are: ocean circulation, ocean modeling, climate changes, and remote sensing among others.
Affonso da Silveira Mascarenhas Junior. Bachelor in Physics (USP/1967), Master in Meteorology (INPE/1976), Master in Physical Oceanography at the Massachusetts Institute of Technology (1979), and Ph.D. in Physical Oceanography at the Massachusetts Institute of Technology (1981). Professor at USP (1972-1993), Universidad Autónoma de Baja California (1993-2008), and nowadays he is director of the International Centre of Investigation on the El Niño Phenomenon. His areas of interest are: observational and synoptical physical oceanography, sea meteorology, and remote sensing.
Luiz Landau. Bachelor in Civil Engineering (PUC-RJ/1973), Master in Civil Engineering (UFRJ/1976), and Doctor of Civil Engineering (UFRJ/1983). At present, he is professor at UFRJ. His area of interest is Civil Engineering, with emphasis in structures, working mainly in the following subjects: finite elements methods, computational methods, oil producing systems, virtual reality, and biomarkers.