Fundamental framework and experiments of the third generation of IAP / LASG world ocean general circulation model

A new generation of the IAP / LASG world ocean general circulation model is designed and presented based on the previous 20-layer model, with enhanced spatial resolutions and improved parameterizations. The model uses a triangular-truncated spectral horizontal grid system with its zonal wave number of 63 (T63) to match its atmospheric counterpart of a T63 spectral atmosphere general circulation model in a planned coupled ocean-atmosphere system. There are 30 layers in vertical direction, of which 20 layers are located above 1000 m for better depicting the permanent thermocline. As previous ocean models developed in IAP / LASG, a free surface (rather than “rigid-lid” approximation) is included in this model. Compared with the 20-layer model, some more detailed physical parameterizations are considered, including the along / cross isopycnal mixing scheme adapted from the Gent-MacWilliams scheme. The model is spun up from a motionless state. Initial conditions for temperature and salinity are taken from the three-dimensional distributions of Levitus’ annual mean observation. A preliminary analysis of the first 1000-year integration of a control experiment shows some encouraging improvements compared with the twenty-layer model, particularly in the simulations of permanent thermocline, thermohaline circu?lation, meridional heat transport, etc. resulted mainly from using the isopycnal mixing scheme. However, the use of isopycnal mixing scheme does not significantly improve the simulated equatorial thermocline. A series of numerical experiments show that the most important contribution to the improvement of equatori?al thermocline and the associated equatorial under current comes from reducing horizontal viscosity in the equatorial regions. It is found that reducing the horizontal viscosity in the equatorial Atlantic Ocean may slightly weaken the overturning rate of North Atlantic Deep Water.     

Centennial oscillations in an ocean-ice coupled model

Irregular centennial oscillations, with a spectral peak at 106 years, were obtained from an ocean-ice coupled model for the North Atlantic with realistic coastline and bottom topography. The model’s thermohaline circulation is forced by mixed boundary conditions, i.e., a Haney-type relaxation condition for temperature, but an equivalent virtual salt flux condition for salinity. All forcing fields are taken from the observed monthly mean climatological wind stress and buoyancy fluxes. The oscillations appeared in the form of a surface-intensified tripole in both the sea surface temperature and salinity fields located in the vicinity of the Labrador Sea. The oscillations involve a delicate interplay between heat and fresh water advection by meridional overturning circulation, horizontal gyres, vertical convection, and the seasonal cycle. The oscillations are primarily control?led by the salinity component of the circulation; however, sea ice plays a minor role in driving the oscillations observed in the model. On the other hand, a regular seasonal cycle in the forcing fields is an important ingredient for the centennial oscillations.     

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...     

A two-dimensional zonally averaged ocean carbon cycle model

An ocean carbon cycle model driven by a constant flow field produced by a two-dimensional thermohaline circu-lation model is developed. Assuming that the biogenic carbon in the oceans is in a dynamic equilibrium, the inorganic carbon cycle is investigated. Before the oceanic uptake of CO2 is carried out, the investigation of 14C distributions in the oceans, including natural and bomb-produced 14C, is conducted by using different values of the exchange coefficient of CO2 for different flow fields (different vertical diffusivities) to test the performance of the model. The suitable values of the exchange coefficient and vertical diffusivities are chosen for the carbon cycle model. Under the forcing of given preindustrial atmospheric CO2 concentration of 280 ppmv, the carbon cycle model is integrated for seven thousand years to reach a steady state. For the human perturbation, two methods including the prescribed at-mospheric pCO2 and prescribed industrial emissions are used in this work. The results from the prescribed atmospher-ic pCO2 show that the oceans take up 36% of carbon dioxide released by human activities for the period of 1980-1989, while the results from the prescribed industrial emission rates show that the oceans take up 34% of car-bon dioxide emitted by industrial sources for the same period. By using the simple method of subtracting industrial emission rate from the total atmosphere+ocean accumulating rate, it can be deduced that before industrial revolution a non-industrial source exists, while after 1940 an extra sink is needed, and that a total non-industrial source of 45 GtC is obtained for the period of 1790-1990.     

A Two-Dimensional Zonally Averaged Ocean Carbon Cycle Model

ＡＴｗｏ－ＤｉｍｅｎｓｉｏｎａｌＺｏｎａｌｙＡｖｅｒａｇｅｄＯｃｅａｎＣａｒｂｏｎＣｙｃｌｅＭｏｄｅｌＸｕＹｏｎｇｆｕ（徐永福）ａｎｄＷａｎｇＭｉｎｇｘｉｎｇ（王明星）ＬＡＰＣ,ＩｎｓｔｉｔｕｔｅｏｆＡｔｍｏｓｐｈｅｒｉｃＰｈｙｓｉｃｓ,Ｃｈｉｎｅｓ...     

Error evolution in the dynamics of an ocean general circulation model

The problem of error propagation is considered for spatially uncorrelated errors of the barotropic stream function in an oceanic general circulation model (OGCM). Such errors typically occur when altimetric data from satellites are assimilated into ocean models. It is shown that the error decays at first due to the dissipation of the smallest scales in the error field. The error then grows exponentially before it saturates at the value corresponding to the difference between independent realizations. A simple analytic formula for the error behavior is derived; it matches the numerical results documented for the present primitive-equation ocean model, and other models in the literature.     

Response of IAP/ LASG GOALS model to the coupling of air-Sea fresh water exchange

l.IntroductionWaterPervadesallofthephysicalanddynamicstructureswithintheEarthclimatesys-temthroughamyriadofhydrologicalprocesses.Thehydrologicalcycleinfluencesclimateinavarietyofways,andthereexistsrobustwaterfluxexchangebetweenoceanandatmosphere(Zhouetal.,l999).Theexchangeoffreshwaterbetweenatmosphereandoceanthroughtheair-seainterfaceisoneoftheleastelementsunderstood,butnowisconsideredasoneofthemostimportantelementsoftheclimatesystem,esPeciallyforoceancirculationchangesondecadaltomillennialt…     

Economic impacts on key Barents Sea fisheries arising from changes in the strength of the Atlantic thermohaline circulation

A bioeconomic model of key fisheries of the Barents Sea is run with scenarios generated by an earth system model of intermediate complexity to assess how the Barents Sea fisheries of cod (Gadus morhua) and capelin (Mallotus villosus) are affected by changes in the Atlantic thermohaline circulation (THC) arising from anthropogenic climate change. Changes in hydrographic conditions have an impact on recruitment success and survival rates, which constitute a lasting effect on the stocks. The economic development of the fisheries is determined for the 21st century, considering a purely stock size based and a coupled stock size-hydrography based harvesting strategy. Results show that a substantial weakening of the THC leads to impaired cod stock development, causing the associated fishery to become unprofitable in the long run. Simultaneous improvements in capelin stock development help the capelin fishery, but are insufficient to offset the losses incurred by the cod fishery. A comparison of harvest strategies reveals that in times of high variability in stock development, coupled stock size-hydrography based management leads to more stable economic results of these fisheries than the stock size based fishing strategy.     

A study on the physical processes of the formation of the ENSO cycle

The physical processes involved in the formation of the ENSO cycle,as well as the possible roles of the Hadley circulation （HC）,Walker circulation （WC）,and the propagating waves of the Southern Oscillation/Northern Oscillation （SO/NO） in its formation,were studied using composite and regression methods.The analysis showed that the convection and heat release triggered by ENSO in the central-eastern equatorial Pacific are the primary drivers for the 3-5 year cycle of the HC,WC and the meridional/zonal circulation.The HC plays a key role in the influence of ENSO on the circulation outside the tropics through angular momentum transportation.Meanwhile,the feedback effects of the anomalous circulation in the mid-high latitudes on ENSO are accomplished by the propagating waves of SO/NO associated with the evolutions of HC and WC.These propagating waves are the main agents of the connections among the meridional/zonal circulation outside the tropics,the Asian/Australian monsoon,the anomalous easterly/westerly winds over the tropical Pacific,and ENSO events.It was found that the 3-5 year cycle of the meridional/zonal circulation forced by ENSO is quite different from the several-week variation of the circulation index triggered by the inner dynamic processes of the atmosphere.The former occurs at the global scale with a definite flow pattern,while the latter occurs only in a wide area without a definite flow pattern.Finally,a physical model for the formation of the ENSO cycle composed of two fundamental processes at the basin and global scale,respectively,is proposed.     

Impact of remote forcing,model resolution and bathymetry on predictions of currents on the shelf

Impacts of remote forcing, model resolution and bathymetry on current predictions at two moorings located on the shelf of the Monterey Bay area are investigated. We consider three Monterey Bay model configurations which differ in resolution and bathymetry representation, and we specify open boundary conditions for these three configurations from two larger scale models, which have different accuracy in the representation of the remote forcing (in the form of poleward propagating along the coast coastally-trapped Kelvin type waves).Comparisons of correlations between observed and model currents as well as visual comparisons show that the most critical element in reproducing currents on the shelf is accurate representation of the remote forcing. Our results also show that accurate representation of bathymetry is the second most critical factor in reproducing observed currents.     

Acoustic thermometry data compared with two ocean models: the importance of Rossby waves and ENSO in modifying the ocean interior

An interpretation is made of interannual changes in acoustic travel time between Oahu and seven receivers at distances of 3000–4000 km. Measurements were made in late 1983, and over two 5-month intervals between 1987 and 1989. Previous publications demonstrated that these changes stem from variations in temperature. Two hydrodynamic ocean models are used to identify plausible oceanic features that could cause these variations. They are from the Naval Research Laboratory and the Florida State University at (1/8)° and (1/6)° resolution, respectively, and are forced with different interannual wind sets for more than a decade. Modelled El Niño's and La Niña's generate poleward travelling Kelvin waves on the eastern boundary of the Pacific. These excite Rossby waves that propagate westward at mid-latitudes. Rossby waves are the dominant model features which affect the modelled acoustic travel times, and hence section-averaged temperatures in the eastern North Pacific. These waves yield travel times whose standard deviations and rates of changes are similar to the measurements. In the observations, some sections separated by less than 500 km exhibit trends in heat content with opposite signs. Similar variability can be explained with modelled Rossby waves. Model wavelengths less than 500 km, eddies, and seasonal cycles induced by seasonal winds yield travel times that are two orders of magnitude too small to account for the data.     

变形的经向环流变化方程及其在诊断阻塞高压形成中的应用

本文推导出可以用来讨论大尺度扰动强迫对平均经圈环流作用的变形平均经圈环流变化方程，并利用此方程对北半球冬季4次阻塞高压的形成进行了诊断分析。分析结果表明:由于波动对经圈环流的强迫作用，使高压脊的西部南风加强，而东部北风加强，从而有利于阻塞高压的形成。文中还分别分析了瞬变波及其低频瞬变波部分在阻塞高压形成中对经圈环流的强迫作用，并讨论了低频瞬变波强迫在太平洋和大西洋地区阻塞高压形成的地域差别。     

Simulation of the Atlantic meridional overturning circulation in an atmosphere-ocean global coupled model. Part II: weakening in a climate change experiment: a feedback mechanism

Most state-of-the art global coupled models simulate a weakening of the Atlantic meridional overturning circulation (MOC) in climate change scenarios but the mechanisms leading to this weakening are still being debated. The third version of the CNRM (Centre National de Recherches Météorologiques) global atmosphere-ocean-sea ice coupled model (CNRM-CM3) was used to conduct climate change experiments for the Intergovernmental Panel on Climate Change Fourth Assessment Report (IPCC AR4). The analysis of the A1B scenario experiment shows that global warming leads to a slowdown of North Atlantic deep ocean convection and thermohaline circulation south of Iceland. This slowdown is triggered by a freshening of the Arctic Ocean and an increase in freshwater outflow through Fram Strait. Sea ice melting in the Barents Sea induces a local amplification of the surface warming, which enhances the cyclonic atmospheric circulation around Spitzberg. This anti-clockwise circulation forces an increase in Fram Strait outflow and a simultaneous increase in ocean transport of warm waters toward the Barents Sea, favouring further sea ice melting and surface warming in the Barents Sea. Additionally, the retreat of sea ice allows more deep water formation north of Iceland and the thermohaline circulation strengthens there. The transport of warm and saline waters toward the Barents Sea is further enhanced, which constitutes a second positive feedback.     

Lyapunov exponents of a simple stochastic model of the thermally and wind-driven ocean circulation

A reformulation of the simple model of the thermally and wind-driven ocean circulation introduced by Maas [Tellus 46A (1994) 671] is considered. Under a realistic range of forcing parameters, this model displays multiple attractors, corresponding to thermally direct and indirect circulations. The fixed point associated with the thermally direct circulation is unstable for a broad range of parameters, leading to limit cycles and chaotic behaviour. It is demonstrated that if weather variability is parameterised as stochastic perturbations to the mechanical and buoyancy fluxes, then the leading Lyapunov exponent of the circulation can become positive for sufficiently strong fluctuations in parameter ranges where it is deterministically zero. If the fluctuations are sufficiently small that the stochastic trajectories are not too far from the deterministic attractor, it is demonstrated that the sign of the leading Lyapunov exponent can have a substantial effect on the predictability of the system.     

Design and Numerical Simulation of an Arctic Ocean Circulation and Thermodynamic Sea－Ice Model

In this paper, the first version of a new Arctic Ocean circulation and thermodynamic sea-ice model is presented by the authors based on the framework of a twenty-layer World Oceanic general circulation model developed by Zhang et al. in 1994, The model’s domain covers the Arctic Ocean and Greenland-Norwegian Seas with the horizon-tal resolution of 200 km × 200 km on a stereographic projection plane. In vertical, the model uses the Eta-coordinate (Sigma modified to have quasi-horizontal coordinate surfaces) and has ten unevenly-spaced layers to cover the deep-est water column of 3000 m. Two 150-year integrations of coupling the ocean circulation model with the sea-ice model have been performed with seasonally cyclic surface boundary conditions. The only difference between the two experiments is in the model’s geography. Some preliminary analyses of the experimental results have been done fo-cused on the following aspects: (1) surface layer temperature, salinity and current; (2) the "Atlantic Layer"; (3) sea-ice cover and its seasonal variation. In comparison with the available observational data, these results are accept-able with reasonable accuracy.     

热带印度洋环流动力与季风相互作用研究进展

围绕中国科学院战略性先导科技专项"泛第三极环境变化与绿色丝绸之路建设"之子课题"热带印度洋环流动力与季风相互作用及其影响",从热带印度洋上层海洋环流、盐度变异对印度洋典型海气耦合过程的影响、印度洋海气相互作用及其对泛第三极水汽输送的影响方面回顾了热带印度洋环流动力与季风相互作用的研究进展.针对国内外研究发展的现状,提出热带印度洋上层经向、纬向、垂向流系间三维联动机制,海盆尺度热盐再分配对局地海-气模态变异的响应和反馈机制,热带印度洋典型海气耦合模态对泛第三极地区气候变化的影响机理等关键科学问题亟待解决.开展该子课题研究的最终目标是：全面认识和理解热带印度洋上层环流体系,加深印度洋海洋环流动力与海气相互作用及其对泛第三极经向水汽输送作用的理解,提高泛第三极地区气候预测水平,提升丝绸之路海上观测航道监测保障能力,从而为"一带一路"倡议和"21世纪海上丝绸之路"建设服务.     

Geostrophic Adjustment Problems in a Polar Basin

The geostrophic adjustment of a homogeneous fluid in a circular basin with idealized topography is addressed using a numerical ocean circulation model and analytical process models. When the basin is rotating uniformly, the adjustment takes place via excitation of boundary propagating waves and when topography is present, via topographic Rossby waves. In the numerically derived solution, the waves are damped because of bottom friction, and a quasi-steady geostrophically balanced state emerges that subsequently spins-down on a long time scale. On the f-plane, numerical quasi-steady state solutions are attained well before the system's mechanical energy is entirely dissipated by friction. It is demonstrated that the adjusted states emerging in a circular basin with a step escarpment or a top hat ridge, centred on a line of symmetry, are equivalent to that in a uniform depth semicircular basin, for a given initial condition. These quasi-steady solutions agree well with linear analytical solutions for the latter case in the inviscid limit.

On the polar plane, the high latitude equivalent to the β-plane, no quasi-steady adjusted state emerges from the adjustment process. At intermediate time scales, after the fast Poincaré and Kelvin waves are damped by friction, the solutions take the form of steady-state adjusted solutions on the f-plane. At longer time scales, planetary waves control the flow evolution. An interesting property of planetary waves on a polar plane is a nearly zero eastward group velocity for the waves with a radial mode higher than two and the resulting formation of eddy-like small-scale barotropic structures that remain trapped near the western side of topographic features.

RÉSUMÉ?Traduit par la rédaction] Nous étudions l'ajustement géostrophique d'un fluide homogène dans un bassin circulaire ayant une topographie idéalisée à l'aide d'un modèle numérique de circulation océanique et de modèles analytiques de processus. Quand le bassin est en rotation uniforme, l'ajustement se fait par l'excitation d'ondes de propagation aux limites, et en présence de topographie, par des ondes de Rossby topographiques. Dans la solution numériquement dérivée, les ondes sont amorties à cause frottement contre le fond, et un état quasi-stable géostrophiquement équilibré s’établit pour ensuite décélérer sur une longue période de temps. Sur le plan f, les solutions numériques d’états quasi-stables sont atteintes bien avant que l’énergie du système mécanique soit entièrement dissipée par le frottement. Nous démontrons que les états ajustés apparaissant dans un bassin circulaire avec un accore en forme de marche ou une crête en forme de merlon, centrés sur une ligne de symétrie, sont équivalents à ceux d'un bassin semi-circulaire de profondeur uniforme pour une condition initiale donnée. Les solutions quasi-stables s'accordent bien avec les solutions analytiques linéaires pour le dernier cas dans la limite de la non-viscosité.

Sur le plan polaire, la haute latitude équivalente au plan β, le processus d'ajustement n'aboutit à aucun état ajusté quasi-stable. Aux échelles de temps intermédiaires, après l'amortissement des ondes rapides de Poincaré et Kelvin par le frottement, les solutions prennent la forme de solutions ajustées d’états stables sur le plan f. Aux échelles de temps plus longues, les ondes planétaires déterminent l’évolution de l’écoulement. Une propriété intéressante des ondes planétaires sur un plan polaire est une vitesse de groupe presque nulle vers l'est pour les ondes ayant un mode radial plus grand que deux et la formation résultante de structures barotropiques de petite échelle de type remous qui demeurent emprisonnées près du bord ouest des caractéristiques topographiques.     

Pathways of intraseasonal variability in the Indonesian Throughflow region

The recent INSTANT measurements in the Indonesian archipelago revealed a broad spectrum of time scales that influence Indonesian Throughflow (ITF) variability, from intraseasonal (20–90 days) to interannual. The different time scales are visible in all transport and property fluxes and are the result of remote forcing by both the Pacific and Indian Ocean winds, and local forcing generated within the regional Indonesian seas. This study focuses on the time-dependent three-dimensional intraseasonal variability (ISV) in the ITF region, in particular at the locations of the INSTANT moorings at the Straits of Lombok, Ombai and Timor. Observations from the INSTANT program in combination with output from the Bluelink ocean reanalysis provide a comprehensive picture about the propagation of ISV in the ITF region. The reanalysis assimilates remotely sensed and in situ ocean observations into an ocean general circulation model to create a hindcast of ocean conditions. Data from the reanalysis and observations from the INSTANT program reveal that deep-reaching subsurface ISV in the eastern Indian Ocean and ITF is closely linked with equatorial wind stress anomalies in the central Indian Ocean. Having traveled more than 5000 km in about 14 days, the associated Kelvin waves can be detected as far east as the Banda Sea. ISV near the Straits of Ombai and Timor is also significantly influenced by local wind forcing from within the ITF region. At the INSTANT mooring sites the ocean reanalysis agrees reasonably well with the observations. Intraseasonal amplitudes are about ±1.0 °C and ±0.5 m/s for potential temperature and velocity anomalies. Associated phases of ISV are very similar in observations and the reanalysis. Where differences exist they can be traced back to likely deficits in the reanalysis, namely the lack of tidal dissipation, insufficient spatial resolution of fine-scale bathymetry in the model in narrow straits or errors in surface forcing.     

The Gulf Stream pathway and the impacts of the eddy-driven abyssal circulation and the Deep Western Boundary Current

A hydrodynamic model of the subtropical Atlantic basin and the Intra-Americas Sea (9–47°N) is used to investigate the dynamics of Gulf Stream separation from the western boundary at Cape Hatteras and its mean pathway to the Grand Banks. The model has five isopycnal Lagrangian layers in the vertical and allows realistic boundary geometry, bathymetry, wind forcing, and a meridional overturning circulation (MOC), the latter specified via ports in the northern and southern boundaries. The northward upper ocean branch of the MOC (14 Sv) was always included but the southward Deep Western Boundary Current (DWBC) was excluded in some simulations, allowing investigation of the impacts of the DWBC and the eddy-driven mean abyssal circulation on Gulf Stream separation from the western boundary. The result is resolution dependent with the DWBC playing a crucial role in Gulf Stream separation at 1/16° resolution but with the eddy-driven abyssal circulation alone sufficient to obtain accurate separation at 1/32° resolution and a realistic pathway from Cape Hatteras to the Grand Banks with minimal DWBC impact except southeast of the Grand Banks. The separation from the western boundary is particularly sensitive to the strength of the eddy-driven abyssal circulation. Farther to the east, between 68°W and the Grand Banks, all of the 1/16° and 1/32° simulations with realistic topography (with or without a DWBC) gave similar generally realistic mean pathways with clear impacts of the topographically constrained eddy-driven abyssal circulation versus very unrealistic Gulf Stream pathways between Cape Hatteras and the Grand Banks from otherwise identical simulations run with a flat bottom, in reduced-gravity mode, or with 1/8° resolution and realistic topography. The model is realistic enough to allow detailed model-data comparisons and a detailed investigation of Gulf Stream dynamics. The corresponding linear solution with a Sverdrup interior and Munk viscous western boundary layers, including one from the northward branch of the MOC, yielded two unrealistic Gulf Stream pathways, a broad eastward pathway centered at the latitude of Cape Hatteras and a second wind plus MOC-driven pathway hugging the western boundary to the north. Thus, a high resolution model capable of simulating an inertial jet is required to obtain a single nonlinear Gulf Stream pathway as it separates from the coast. None of the simulations were sufficiently inertial to overcome the linear solution need for a boundary current north of Cape Hatteras without assistance from pathway advection by the abyssal circulation, even though the core speeds of the simulated currents were consistent with observations near separation. In the 1/16° simulation with no DWBC and a 1/32° simulation with high bottom friction and no DWBC the model Gulf Stream overshot the observed separation latitude. With abyssal current assistance the simulated (and the observed) mean Gulf Stream pathway between separation from the western boundary and 70°W agreed closely with a constant absolute vorticity (CAV) trajectory influenced by the angle of the coastline prior to separation. The key abyssal current crosses under the Gulf Stream at 68.5–69°W and advects the Gulf Stream pathway southward to the terminus of an escarpment in the continental slope. There the abyssal current crosses to deeper depths to conserve potential vorticity while passing under the downward-sloping thermocline of the stream and then immediately retroflects eastward onto the abyssal plain, preventing further southward pathway advection. Thus specific topographic features and feedback from the impact of the Gulf Stream on the abyssal current pathway determined the latitude of the stream at 68.5–69°W, a latitude verified by observations. The associated abyssal current was also verified by observations.