Fundamental Framework and Experiments of the Third Generation of IAP/LASG World Ocean General Circulation Model

1.IntroductionThefirstbaroclinicoceanicgeneralcirculationmodel(OGCM)developedattheStateKeyLaboratoryofNumericalModelingforAtmosphericSciencesandGeophysicalFluidDynamics(LASG),InstituteofAtmosphericPhysics(IAP)isafour--layermodelwithitshorizontalresol...     

Sensitivity of tropical Atlantic climate to mixing in a coupled ocean–atmosphere model

The sensitivity of tropical Atlantic climate to upper ocean mixing is investigated using an ocean-only model and a coupled ocean–atmosphere model. The upper ocean thermal structure and associated atmospheric circulation prove to be strongly related to the strength of upper ocean mixing. Using the heat balance in the mixed layer it is shown that an excessively cold equatorial cold tongue can be attributed to entrainment flux at the base of the oceanic mixed layer, that is too large. Enhanced entrainment efficiency acts to deepen the mixed layer and causes strong reduction in the upper ocean divergence in the central equatorial Atlantic. As a result, the simulated sea surface temperature, thermocline structure, and upwelling velocities are close to the observed estimates. In the coupled model, the seasonal migration of the Intertropical Convergence Zone (ITCZ) reduces when the entrainment efficiency in the oceanic mixed layer is enhanced. The precipitation rates decrease in the equatorial region and increase along 10°N, resulting in a more realistic Atlantic Marine ITCZ. The reduced meridional surface temperature gradient in the eastern tropical Atlantic prohibits the development of convective precipitation in the southeastern part of the tropical Atlantic. Also, the simulation of tropical Atlantic variability as expressed in the meridional gradient mode and the eastern cold tongue mode improves when the entrainment efficiency is enhanced.     

North African climate variability. Part 1: Tropical thermocline coupling

Summary Tropical ocean thermocline variability is studied using gridded data assimilated by an ocean model in the period 1950–2000. The dominant patterns and variability are identified using EOF analysis applied to E–W depth slices of sea temperatures averaged over the tropics. After removing the annual cycle, an east–west ‘see-saw’ with an interannual to decadal rhythm is the leading mode in each of the tropical basins. In the case of the leading mode in the Pacific, the thermocline oscillation forms a dipole structure, but in the (east) Atlantic and (southwest) Indian Ocean there is a single center of action. The interaction of the ocean thermocline and atmospheric Walker circulations is studied through cross-modulus analysis of wavelet-filtered EOF time scores. Our study demonstrates how tropical ocean thermocline variability contributes to zonal circulation anomalies in the atmosphere. The equatorial Pacific thermocline oscillation explains 62 and 53% of the variability of the Pacific and Atlantic zonal overturning circulations, the latter driving convective polarity between North Africa and South America. The Pacific sea-saw leads the Atlantic zonal circulation by a few months.     

Simulation of ~(14)C in IAP/LASG L30T63 Ocean Model

1. Introduction14C is a radioactive isotope of carbon, whose halftime is about 5730 yr. Under natural circumstance,14C comes into being in the stratosphere because of thenitrogen explosion of cosmic radial. 14C is mixed inthe atmosphere, absorbed by oceans, and then trans-ported into deep ocean. 14C radioactively reduces asits age increases. The reduction process of 14C canuncover the evolvement process of sea water's age andrenewal time.14 C content is usually described with one in athousan…     

一个具有高分辨率海洋分量的海气耦合模式

发展了中国科学院大气物理研究所（IAP）大气科学与地球流体力学数值模拟国家重点实验室(LASG)全球海-陆-气耦合系统模式的一个新版本（GOALS-5），其中海洋分量是基于30层高分辨率的海洋环流模式。对该耦合模式成功地进行了30年的积分，基本上克服了气候飘移。与GOALS较早的几个版本模拟的SST相比，北半球夏季中高纬地区的误差显著改善，可能主要是由于采用了GM90的等密度混合方案使经向热输送增强的结果。对赤道中东太平洋冷暖事件有较好的模拟能力，并且对SST变率的模拟在位置和强度上比旧版本有明显改善。     

An investigation of tropical Atlantic bias in a high-resolution coupled regional climate model

Coupled atmosphere–ocean general circulation models (AOGCMs) commonly fail to simulate the eastern equatorial Atlantic boreal summer cold tongue and produce a westerly equatorial trade wind bias. This tropical Atlantic bias problem is investigated with a high-resolution (27-km atmosphere represented by the Weather Research and Forecasting Model, 9-km ocean represented by the Regional Ocean Modeling System) coupled regional climate model. Uncoupled atmospheric simulations test climate sensitivity to cumulus, land-surface, planetary boundary layer, microphysics, and radiation parameterizations and reveal that the radiation scheme has a pronounced impact in the tropical Atlantic. The CAM radiation simulates a dry precipitation (up to ?90%) and cold land-surface temperature (up to ?8?K) bias over the Amazon related to an over-representation of low-level clouds and almost basin-wide westerly trade wind bias. The Rapid Radiative Transfer Model and Goddard radiation simulates doubled Amazon and Congo Basin precipitation rates and a weak eastern Atlantic trade wind bias. Season-long high-resolution coupled regional model experiments indicate that the initiation of the warm eastern equatorial Atlantic sea surface temperature (SST) bias is more sensitive to the local rather than basin-wide trade wind bias and to a wet Congo Basin instead of dry Amazon—which differs from AOGCM simulations. Comparisons between coupled and uncoupled simulations suggest a regional Bjerknes feedback confined to the eastern equatorial Atlantic amplifies the initial SST, wind, and deepened thermocline bias, while barrier layer feedbacks are relatively unimportant. The SST bias in some CRCM simulations resembles the typical AOGCM bias indicating that increasing resolution is unlikely a simple solution to this problem.     

任意正交曲线坐标系下的海洋模式动力框架的发展与评估

本文发展了一个可以适用于任意水平正交曲线坐标系的海洋模式动力框架,并将其应用于中国科学院大气物理研究所大气科学和地球流体力学数值模拟国家重点实验室发展的气候系统海洋模式LICOM2.0（LASG/IAP Climate system Ocean Model,version2.0）。在经纬网格坐标系下,新的动力框架与LICOM2.0原有的动力框架模拟结果完全一致。基于新的动力框架,海洋模式可采用能够准确描述北冰洋地形的三极网格,克服了LICOM2.0经纬网格版本必须将北极点处理为孤岛的缺陷,从而显著改进了模式对于北冰洋环流和北大西洋经圈翻转流函数（AMOC）的模拟能力。此外,引进三极网格还可以避免模式网格距随纬度增加而急剧减小带来的计算不稳定,在LICOM2.0的三极网格版本中,模式不需要采用任何空间滤波方案仍然能够保证计算的稳定性,从而与LICOM2.0的经纬网格版本相比,极大地提高了模式的并行效率,这一点在当水平分辨率提高到0.1度时表现得尤为明显,海洋模式的并行加速比可以从经纬网格版本的5.8左右提高到三极网格版本的15.0左右。     

Estimates of Anthropogenic CO_2 Uptake in a Global Ocean Model

A global ocean general circulation model (L30T63) is employed to study the uptake and distribution of anthropogenic CO2 in the ocean. A subgrid-scale mixing scheme called GM90 is used in the model. There are two main GM90 parameters including isopycnal diffusivity and skew (thickness) diffusivity. Sensitivities of the ocean circulation and the redistribution of dissolved anthropogenic CO2 to these two parameters are examined. Two runs estimate the global oceanic anthropogenic CO2 uptake to be 1.64 and 1.73 ...     

Seasonal variability of the equatorial undercurrent termination and associated salinity maximum in the Gulf of Guinea

The termination of the Equatorial Undercurrent (EUC) in the eastern equatorial Atlantic during boreal summer and fall, and the fate of the associated saline water masses, are analyzed from in situ hydrological and currents data collected during 19 hydrographic cruises between 2000 and 2007, complemented by observations from Argo profiling floats and PIRATA moorings, and from a numerical simulation of the Tropical Atlantic Ocean for the period 1993–2007. An intense variability of the circulation and hydrological properties is evidenced from observations in the upper thermocline (24.5–26.2 isopycnal layer) between June and November. During early boreal summer, saline water masses are transported eastward in the upper thermocline to the African coast within the EUC, and recirculate westward on both sides of the EUC. In mid-boreal summer, the EUC weakens in the upper thermocline and the equatorial salinity maximum disappears due to intense mixing with the surface waters during the upwelling season. The extra-equatorial salinity maxima are also partially eroded during the boreal summer, with a slight poleward migration of the southern hemisphere maximum until late boreal summer. The upper EUC reappears in September, feeding again the eastern equatorial Atlantic with saline waters until boreal spring. During December–January, numerical results suggest a second seasonal weakening of the EUC in the Gulf of Guinea, with a partial erosion of the associated equatorial salinity maximum.     

The influence of numerical advection schemes on the results of ocean general circulation models

The dependence of results from coarse-resolution models of the North Atlantic circulation on the numerical advection algorithm is studied. In particular, the sensitivity of parameters relevant for climate simulations as e.g., meridional transport of mass and heat and main thermocline thickness is investigated. Three algorithms were considered: (a) a central difference scheme with constant values for horizontal and vertical diffusion, (b) an upstream scheme with no explicit diffusion, and (c) a flux-corrected transport (FCT) scheme with constant and strictly isopycnal diffusion. The temporal evolution of the three models on time scales of centuries is markedly different, the upstream scheme resulting in much shorter adjustment time whereas the central difference scheme is slower and controlled by vertical diffusion rather than advection. In the steady state, the main thermocline structure is much less diffusive in the FCT calculation which also has much lower heat transport. Both horizontal circulation and overturning in the meridional-vertical plane are strongest in the upstream-model. The results are discussed in terms of the effective vertical (diapycnal) mixing in the different models. A significant increase in vertical resolution would be required to eliminate the high sensitivity due to the numerical algorithms, and allow physically motivated mixing formulations to become effective.This paper was presented at the International Conference on Modelling of Global Climate Change and Variability, held in Hamburg 11–15 September 1989 under the auspices of the Meteorological Institute of the University of Hamburg and the Max Planck Institute for Meteorology. Guest Editor for these papers is Dr. L. Dümenil     

Estimates of anthropogenic CO<Subscript>2</Subscript> uptake in a global ocean model

A global ocean general circulation model (L30T63) is employed to study the uptake and distribution of anthropogenic CO₂ in the ocean. A subgrid-scale mixing scheme called GM90 is used in the model. There are two main GM90 parameters including isopycnal diffusivity and skew (thickness) diffusivity. Sensitivities of the ocean circulation and the redistribution of dissolved anthropogenic CO₂ to these two parameters are examined. Two runs estimate the global oceanic anthropogenic CO₂ uptake to be 1.64 and 1.73 Pg C yr^-1 for the 1990s, and that the global ocean contained 86.8 and 92.7 Pg C of anthropogenic CO₂ at the end of 1994, respectively. Both the total inventory and uptake from our model are smaller than the data-based estimates. In this presentation, the vertical distributions of anthropogenic CO₂ at three meridional sections are discussed and compared with the available data-based estimates. The inventory in the individual basins is also calculated. Use of large isopycnal diffusivity can generally improve the simulated results, including the exchange flux, the vertical distribution patterns, inventory, storage, etc. In terms of comparison of the vertical distributions and column inventory, we find that the total inventory in the Pacific Ocean obtained from our model is in good agreement with the data-based estimate, but a large difference exists in the Atlantic Ocean, particularly in the South Atlantic. The main reasons are weak vertical mixing and that our model generates small exchange fluxes of anthropogenic CO₂ in the Southern Ocean. Improvement in the simulation of the vertical transport and sea ice in the Southern Ocean is important in future work.     

Impact of vertical mixing induced by small vertical scale structures above and within the equatorial thermocline on the tropical Pacific in a CGCM

Oceanic vertical mixing is known to influence the state of the equatorial ocean which affects the climate system, including the amplitude of El Niño/Southern Oscillation (ENSO). Recent measurements of ocean currents at high vertical resolution capture numerous small vertical scale structures (SVSs) within and above the equatorial thermocline that contribute significantly to vertical mixing but which are not sufficiently resolved by coarse resolution ocean models. We investigate the impact of the vertical mixing induced by the SVSs on the mean state and interannual variability in the tropical Pacific by using a coupled general circulation model. The vertical mixing induced by the SVSs is represented as an elevated vertical diffusivity from the surface down to the 20 °C isotherm depth, a proxy for the depth of the thermocline. We investigate different forms for the elevated mixing. It is found that the SVS-induced mixing strongly affect the mean state of the ocean leading to a warming of sea surface temperature (SST) and associated deepening and sharpening of the thermocline in the eastern equatorial Pacific. We find that the SST warming induced by the elevated mixing is further strengthened through the Bjerknes feedback and SST-shortwave flux feedback. We also find a reduction in the number of large amplitude ENSO events and in certain cases an increase in the skewness of ENSO.     

潮汐混合对大西洋经圈翻转环流(AMOC)模拟影响的数值模拟研究

海洋中的潮汐混合对大西洋经圈翻转环流AMOC(Atlantic Meridional Overturning Circulation)模拟的影响是海洋环流模式研究的热点问题之一。本文采用IAP/LASG发展的气候系统海洋模式LICOM(LASG/IAP Climate system Ocean Model)及与海冰耦合模式进行了有无潮汐混合方案的试验,重点探讨了潮汐混合对AMOC强度模拟的影响。结果显示,引入潮汐混合后模拟的AMOC强度极大值比对照试验增加约1倍,更接近RAPID(Rapid Climate Change Programme)观测。而且,潮汐混合试验中模拟的AMOC上层环流深度(3200 m)比对照试验加深1000 m左右,同样更接近RAPID观测。海洋底部的垂直混合增强,使海洋层结变得更加不稳定,加强了北大西洋高纬地区,特别是拉布拉多海等地区的深对流,这是AMOC加强的直接原因。同时,潮汐混合试验中上层海洋环流也加强,增加了中低纬副热带高盐海水向高纬输送,使表层增密,海洋层结更加不稳定,也可以进一步增强AMOC。     

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     

A Coupling Experiment of an Atmosphere and an Ocean Model with a Monthly Anomaly Exchange Scheme

A nine-layer spectral atmospheric general circulation model is coupled to a twenty-layer global oceanic general circulation model with the “prediction-correction” monthly anomaly exchange scheme which has been proposed at the Institute of Atmospheric Physics (IAP). A forty-year integration of the coupled model shows that the CGCM is fairly successful in keeping a reasonable pattern of the modelled SST although most of the Pacific become warmer than those given by the uncoupled ocean model. The model tends to reach a more realistic state than the uncoupled one in terms of downward surface heat flux into ocean particularly in the equatorial Pacific region. Also, the model is capable to simulate interannual variability of sea surface temperature in tropical region.     

Transports and budgets of volume,heat, and salt from a global eddy-resolving ocean model

The results from an integration of a global ocean circulation model have been condensed into an analysis of the volume, heat, and salt transports among the major ocean basins. Transports are also broken down between the model's Ekman, thermocline, and deep layers. Overall, the model does well. Horizontal exchanges of mass, heat, and salt between ocean basins have reasonable values; and the volume of North Atlantic Deep Water (NADW) transport is in general agreement with what limited observations exist. On a global basis the zonally integrated meridional heat transport is poleward at all latitudes except for the latitude band 30°S to 45°S. This anomalous transport is most likely a signature of the model's inability to form Antarctic Intermediate (AAIW) and Antarctic bottom water (AABW) properly. Eddy heat transport is strong at the equator where its convergence heats the equatorial Pacific about twice as much as it heats the equatorial Atlantic. The greater heating in the Pacific suggests that mesoscale eddies may be a vital mechanism for warming and maintaining an upwelling portion of the global conveyor-belt circulation. The model's fresh water transport compares well with observations. However, in the Atlantic there is an excessive southward transport of fresh water due to the absence of the Mediterranean outflow and weak northward flow of AAIW. Eddies in the mid-latitudes act to redistribute heat and salt down the mean gradients. Residual fluxes calculated from a sum of the computed advective (including eddies), forced, and stored fluxes of heat and salt represent transport mostly due to vertical sub-grid scale mixing processes. Perhaps the model's greatest weakness is the lack of strong AAIW and AABW circulation cells. Accurate thermohaline forcing in the North Atlantic (based on numerous hydrographic observations) helps the model adequately produce NADW. In contrast, the southern ocean is an area of sparse observation. Better thermohaline observations in this area may be needed if models such as this are to produce the deep convection that will achieve more accurate simulations of the global 3-dimensional circulation.     

A fast version of LASG/IAP climate system model and its 1000-year control integration

A fast version of the State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geo- physical Fluid Dynamics(LASG)/Institute of Atmospheric Physics(IAP)climate system model is briefly documented.The fast coupled model employs a low resolution version of the atmospheric component Grid Atmospheric Model of IAP/LASG(GAMIL),with the other parts of the model,namely an oceanic component LASG/IAP Climate Ocean Model(LICOM),land component Common Land Model(CLM),and sea ice component from National Center for Atmospheric Research Community Climate System Model (NCAR CCSM2),as the same as in the standard version of LASG/IAP Flexible Global Ocean Atmosphere Land System model(FGOALS g).The parameterizations of physical and dynamical processes of the at- mospheric component in the fast version are identical to the standard version,although some parameter values are different.However,by virtue of reduced horizontal resolution and increased time-step of the most time-consuming atmospheric component,it runs faster by a factor of 3 and can serve as a useful tool for long- term and large-ensemble integrations.A 1000-year control simulation of the present-day climate has been completed without flux adjustments.The final 600 years of this simulation has virtually no trends in global mean sea surface temperatures and is recommended for internal variability studies.Several aspects of the control simulation’s mean climate and variability are evaluated against the observational or reanalysis data. The strengths and weaknesses of the control simulation are evaluated.The mean atmospheric circulation is well simulated,except in high latitudes.The Asian-Australian monsoonal meridional cell shows realistic features,however,an artificial rainfall center is located to the eastern periphery of the Tibetan Plateau persists throughout the year.The mean bias of SST resembles that of the standard version,appearing as a"double ITCZ"(Inter-Tropical Convergence Zone)associated with a westward extension of the equatorial eastern Pacific cold tongue.The sea ice extent is acceptable but has a higher concentration.The strength of Atlantic meridional overturning is 27.5 Sv.Evidence from the 600-year simulation suggests a modulation of internal variability on ENSO frequency,since both regular and irregular oscillations of ENSO are found during the different time periods of the long-term simulation.     

Numerical Experiment for the Impact of the Ozone Hole over Antarctica on the Global Climate

A numerical experiment has been carried out with IAP (Institute of Atmospheric Physics) 9-layer general circulation model to investigate the influence of the Antarctic Ozone Hole on the global climate. The results show that the changes of total amount of ozone over higher latitude and polar region of the Southern Hemisphere affect not only the climate in the Southern Hemisphere, but also that in the Northern Hemisphere significantly. In the next spring, although the total amount of ozone over Antarctica has returned to the normal value, the influences of Ozone Hole still exist. Suppported by LASG and the National Key Project of Fundamental Research “Climate Dynamics and Climate Prediction Theory.≓     

The Annual Modes of Tropical Precipitation Simulated by the LASG/IAP Coupled Ocean-Atmosphere Model FGOALS_s1.1

This paper evaluates the performance of a coupled general circulation model FGOALS_s1.1 developed by LASG/IAP in simulating the annual modes of tropical precipitation.To understand the impacts of air-sea coupling on the annual modes,the result of an off-line simulation of the atmospheric component of FGOALS_s1.1,i.e.,LASG/IAP atmospheric general circulation model SAMIL,is also analyzed.FGOALS_s1.1 can reasonably reproduce major characteristics of the annual mean precipitation.Nonetheless,the coupled model shows overestimation of precipitation over the equatorial Pacific and tropical South Pacific,and underestimation of precipitation over the northern equatorial Pacific.The monsoon mode simulated by FGOALS_s1.1 shows an equatorial anti-symmetric structure,which is consistent with the observation.The bias of the coupled model in simulating monsoon mode resembles that of SAMIL,especially over the subtropics.The main deficiency of FGOALS_s1.1 is its failure in simulating the spring-fall asymmetric mode.This is attributed to the false phase of sea surface temperature anomaly (SSTA) annual cycleover the equatorial central-castern Pacific and Indian Ocean,which leads to the bias of the Walker circulation over the equatorial Pacific and the anti-Walker circulation over the Indian Ocean in boreal spring and fall.In addition,the domains of the western North Pacific monsoon and Indian monsoon simulated by the coupled model are smaller than the observation.The study suggests that the bias of the fully coupled oceanatmosphere model can only be partly attributed to the bias of the atmospheric component.The performance of FGOALS-s1.1 in simulating the annual cycle of equatorial SST deserves further improvement.     

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.