Impact of Resolution and Optimized ECCO Forcing on Simulations of the Tropical Pacific

2008 ◽  
Vol 25 (1) ◽  
pp. 131-147 ◽  
Author(s):  
I. Hoteit ◽  
B. Cornuelle ◽  
V. Thierry ◽  
D. Stammer

Abstract The sensitivity of the dynamics of a tropical Pacific Massachusetts Institute of Technology (MIT) general circulation model (MITgcm) to the surface forcing fields and to the horizontal resolution is analyzed. During runs covering the period 1992–2002, two different sets of surface forcing boundary conditions are used, obtained 1) from the National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) reanalysis project and 2) from the Estimating the Circulation and Climate of the Ocean (ECCO) assimilation consortium. The “ECCO forcing” is the “NCEP forcing” adjusted by a state estimation procedure using the MITgcm with a 1° × 1° global grid and the adjoint method assimilating a multivariate global ocean dataset. The skill of the model is evaluated against ocean observations available in situ and from satellites. The model domain is limited to the tropical Pacific, with open boundaries located along 26°S, 26°N, and in the Indonesian throughflow. To account for large-scale changes of the ocean circulation, the model is nested in the global time-varying ocean state provided by the ECCO consortium on a 1° grid. Increasing the spatial resolution to 1/3° and using the ECCO forcing fields significantly improves many aspects of the circulation but produces overly strong currents in the western model domain. Increasing the resolution to 1/6° does not yield further improvements of model results. Using the ECCO heat and freshwater fluxes in place of NCEP products leads to improved time-mean model skill (i.e., reduced biases) over most of the model domain, underlining the important role of adjusted heat and freshwater fluxes for improving model representations of the tropical Pacific. Combinations of ECCO and NCEP wind forcing fields can improve certain aspects of the model solutions, but neither ECCO nor NCEP winds show clear overall superiority.

2007 ◽  
Vol 15 (2) ◽  
pp. 107-115 ◽  
Author(s):  
Chris Hill ◽  
Dimitris Menemenlis ◽  
Bob Ciotti ◽  
Chris Henze

Up to 1920 processors of a cluster of distributed shared memory machines at the NASA Ames Research Center are being used to simulate ocean circulation globally at horizontal resolutions of 1/4, 1/8, and 1/16-degree with the Massachusetts Institute of Technology General Circulation Model, a finite volume code that can scale to large numbers of processors. The study aims to understand physical processes responsible for skill improvements as resolution is increased and to gain insight into what resolution is sufficient for particular purposes. This paper focuses on the computational aspects of reaching the technical objective of efficiently performing these global eddy-resolving ocean simulations. At 1/16-degree resolution the model grid contains 1.2 billion cells. At this resolution it is possible to simulate approximately one month of ocean dynamics in about 17 hours of wallclock time with a model timestep of two minutes on a cluster of four 512-way NUMA Altix systems. The Altix systems' large main memory and I/O subsystems allow computation and disk storage of rich sets of diagnostics during each integration, supporting the scientific objective to develop a better understanding of global ocean circulation model solution convergence as model resolution is increased.


2007 ◽  
Vol 37 (2) ◽  
pp. 277-295 ◽  
Author(s):  
A. J. Meijers ◽  
N. L. Bindoff ◽  
J. L. Roberts

Abstract The large-scale volume, heat, and freshwater ocean transports in the Southern Hemisphere are investigated using time-averaged output from a seasonless, high-resolution general circulation model. The ocean circulation is realistic, and property transports are comparable to observations. The Antarctic Circumpolar Current (ACC) carries 144 Sv (Sv ≡ 106 m3 s−1) of water eastward across Drake Passage, increasing to 155 Sv south of Australia because of the Indonesian Throughflow (ITF). There is a clear Indo-Pacific gyre around Australia exchanging −10 Sv, 0.9 PW of heat, and 0.2 Sv of freshwater through the ITF, and there is a 9-Sv leakage from the Tasman Sea to the Indian Ocean. The transport of heat and freshwater by eddies is localized to the upper 1000 m of the water column and specific regions, such as western boundary currents, confluences, and the subantarctic front (SAF). Eddy transport of heat and freshwater is negligible in gyre interiors and south of the SAF but is vital across the northern edge of the ACC, in particular at the Agulhas Retroflection where eddies accomplish almost 100% of the net ocean heat and 60% of the southward freshwater transport. The eddy transport is almost zero across the latitude of Drake Passage while in a quasi-Lagrangian frame eddy transports are significant across the ACC but surprisingly are still smaller than the mean transport of heat. Mean and eddy property transport divergences are found to be strongly compensating in areas of high eddy activity. This is caused by increased baroclinic instability in strong mean flows, which induces an opposing eddy transport. This relationship is observed to be stronger in the case of horizontal heat transport than in corresponding horizontal freshwater transports.


1997 ◽  
Vol 25 ◽  
pp. 111-115 ◽  
Author(s):  
Achim Stössel

This paper investigates the long-term impact of sea ice on global climate using a global sea-ice–ocean general circulation model (OGCM). The sea-ice component involves state-of-the-art dynamics; the ocean component consists of a 3.5° × 3.5° × 11 layer primitive-equation model. Depending on the physical description of sea ice, significant changes are detected in the convective activity, in the hydrographic properties and in the thermohaline circulation of the ocean model. Most of these changes originate in the Southern Ocean, emphasizing the crucial role of sea ice in this marginally stably stratified region of the world's oceans. Specifically, if the effect of brine release is neglected, the deep layers of the Southern Ocean warm up considerably; this is associated with a weakening of the Southern Hemisphere overturning cell. The removal of the commonly used “salinity enhancement” leads to a similar effect. The deep-ocean salinity is almost unaffected in both experiments. Introducing explicit new-ice thickness growth in partially ice-covered gridcells leads to a substantial increase in convective activity, especially in the Southern Ocean, with a concomitant significant cooling and salinification of the deep ocean. Possible mechanisms for the resulting interactions between sea-ice processes and deep-ocean characteristics are suggested.


2007 ◽  
Vol 20 (4) ◽  
pp. 765-771 ◽  
Author(s):  
Markus Jochum ◽  
Clara Deser ◽  
Adam Phillips

Abstract Atmospheric general circulation model experiments are conducted to quantify the contribution of internal oceanic variability in the form of tropical instability waves (TIWs) to interannual wind and rainfall variability in the tropical Pacific. It is found that in the tropical Pacific, along the equator, and near 25°N and 25°S, TIWs force a significant increase in wind and rainfall variability from interseasonal to interannual time scales. Because of the stochastic nature of TIWs, this means that climate models that do not take them into account will underestimate the strength and number of extreme events and may overestimate forecast capability.


2019 ◽  
Vol 12 (8) ◽  
pp. 3745-3758 ◽  
Author(s):  
François Massonnet ◽  
Antoine Barthélemy ◽  
Koffi Worou ◽  
Thierry Fichefet ◽  
Martin Vancoppenolle ◽  
...  

Abstract. The ice thickness distribution (ITD) is one of the core constituents of modern sea ice models. The ITD accounts for the unresolved spatial variability of sea ice thickness within each model grid cell. While there is a general consensus on the added physical realism brought by the ITD, how to discretize it remains an open question. Here, we use the ocean–sea ice general circulation model, Nucleus for European Modelling of the Ocean (NEMO) version 3.6 and Louvain-la-Neuve sea Ice Model (LIM) version 3 (NEMO3.6-LIM3), forced by atmospheric reanalyses to test how the ITD discretization (number of ice thickness categories, positions of the category boundaries) impacts the simulated mean Arctic and Antarctic sea ice states. We find that winter ice volumes in both hemispheres increase with the number of categories and attribute that increase to a net enhancement of basal ice growth rates. The range of simulated mean winter volumes in the various experiments amounts to ∼30 % and ∼10 % of the reference values (run with five categories) in the Arctic and Antarctic, respectively. This suggests that the way the ITD is discretized has a significant influence on the model mean state, all other things being equal. We also find that the existence of a thick category with lower bounds at ∼4 and ∼2 m for the Arctic and Antarctic, respectively, is a prerequisite for allowing the storage of deformed ice and therefore for fostering thermodynamic growth in thinner categories. Our analysis finally suggests that increasing the resolution of the ITD without changing the lower limit of the upper category results in small but not negligible variations of ice volume and extent. Our study proposes for the first time a bi-polar process-based explanation of the origin of mean sea ice state changes when the ITD discretization is modified. The sensitivity experiments conducted in this study, based on one model, emphasize that the choice of category positions, especially of thickest categories, has a primary influence on the simulated mean sea ice states while the number of categories and resolution have only a secondary influence. It is also found that the current default discretization of the NEMO3.6-LIM3 model is sufficient for large-scale present-day climate applications. In all cases, the role of the ITD discretization on the simulated mean sea ice state has to be appreciated relative to other influences (parameter uncertainty, forcing uncertainty, internal climate variability).


1998 ◽  
Vol 11 (8) ◽  
pp. 1997-2015 ◽  
Author(s):  
Bing Ye ◽  
Anthony D. Del Genio ◽  
Kenneth K-W. Lo

Abstract Observed variations of convective available potential energy (CAPE) in the current climate provide one useful test of the performance of cumulus parameterizations used in general circulation models (GCMs). It is found that frequency distributions of tropical Pacific CAPE, as well as the dependence of CAPE on surface wet-bulb potential temperature (Θw) simulated by the Goddard Institute for Space Studies’s GCM, agree well with that observed during the Australian Monsoon Experiment period. CAPE variability in the current climate greatly overestimates climatic changes in basinwide CAPE in the tropical Pacific in response to a 2°C increase in sea surface temperature (SST) in the GCM because of the different physics involved. In the current climate, CAPE variations in space and time are dominated by regional changes in boundary layer temperature and moisture, which in turn are controlled by SST patterns and large-scale motions. Geographical thermodynamic structure variations in the middle and upper troposphere are smaller because of the canceling effects of adiabatic cooling and subsidence warming in the rising and sinking branches of the Walker and Hadley circulations. In a forced equilibrium global climate change, temperature change is fairly well constrained by the change in the moist adiabatic lapse rate and thus the upper troposphere warms to a greater extent than the surface. For this reason, climate change in CAPE is better predicted by assuming that relative humidity remains constant and that the temperature changes according to the moist adiabatic lapse rate change of a parcel with 80% relative humidity lifted from the surface. The moist adiabatic assumption is not symmetrically applicable to a warmer and colder climate: In a warmer regime moist convection determines the tropical temperature structure, but when the climate becomes colder the effect of moist convection diminishes and the large-scale dynamics and radiative processes become relatively important. Although a prediction based on the change in moist adiabat matches the GCM simulation of climate change averaged over the tropical Pacific basin, it does not match the simulation regionally because small changes in the general circulation change the local boundary layer relative humidity by 1%–2%. Thus, the prediction of regional climate change in CAPE is also dependent on subtle changes in the dynamics.


2021 ◽  
Vol 14 (5) ◽  
pp. 2781-2799
Author(s):  
Pengfei Wang ◽  
Jinrong Jiang ◽  
Pengfei Lin ◽  
Mengrong Ding ◽  
Junlin Wei ◽  
...  

Abstract. A high-resolution (1/20∘) global ocean general circulation model with graphics processing unit (GPU) code implementations is developed based on the LASG/IAP Climate System Ocean Model version 3 (LICOM3) under a heterogeneous-compute interface for portability (HIP) framework. The dynamic core and physics package of LICOM3 are both ported to the GPU, and three-dimensional parallelization (also partitioned in the vertical direction) is applied. The HIP version of LICOM3 (LICOM3-HIP) is 42 times faster than the same number of CPU cores when 384 AMD GPUs and CPU cores are used. LICOM3-HIP has excellent scalability; it can still obtain a speedup of more than 4 on 9216 GPUs compared to 384 GPUs. In this phase, we successfully performed a test of 1/20∘ LICOM3-HIP using 6550 nodes and 26 200 GPUs, and on a large scale, the model's speed was increased to approximately 2.72 simulated years per day (SYPD). By putting almost all the computation processes inside GPUs, the time cost of data transfer between CPUs and GPUs was reduced, resulting in high performance. Simultaneously, a 14-year spin-up integration following phase 2 of the Ocean Model Intercomparison Project (OMIP-2) protocol of surface forcing was performed, and preliminary results were evaluated. We found that the model results had little difference from the CPU version. Further comparison with observations and lower-resolution LICOM3 results suggests that the 1/20∘ LICOM3-HIP can reproduce the observations and produce many smaller-scale activities, such as submesoscale eddies and frontal-scale structures.


1997 ◽  
Vol 25 ◽  
pp. 116-120 ◽  
Author(s):  
S. Legutke ◽  
E. Maier-Reimkr ◽  
A. Stössel ◽  
A. Hellbach

A global ocean general circulation model has been coupled with a dynamic thermodynamic sea-ice model. This model has been spun-up in a 1000 year integration using daily atmosphere model data. Main water masses and currents are reproduced as well as the seasonal characteristics of the ice cover of the Northern and Southern Hemispheres. Model results for the Southern Ocean, however, show the ice cover as too thin, and there are large permanent polynyas in the Weddell and Ross Seas. These polynyas are due to a large upward oceanic heat flux caused by haline rejection during the freezing of sea ice. Sensitivity studies were performed to test several ways of treating the sea-surface salinity and the rejected brine. The impact on the ice cover, water-mass characteristics, and ocean circulation are described.


Sign in / Sign up

Export Citation Format

Share Document