Due to the decrease in grid size associated with the convergence of meridians toward the poles inspherical coordinates, the time steps in many global climate models with finite-difference method are restrictedto be un...Due to the decrease in grid size associated with the convergence of meridians toward the poles inspherical coordinates, the time steps in many global climate models with finite-difference method are restrictedto be unpleasantly small. To overcome the problem, a reduced grid is introduced to LASG/IAP world oceangeneral circulation models. The reduced grid is implemented successfully in the coarser resolutions versionmodel L30T63 at first. Then, it is carried out in the improved version model LICOM with finer resolutions. Inthe experiment with model L30T63, under time step unchanged though, execution time per single model run isshortened significantly owing to the decrease of grid number and filtering execution in high latitudes. Resultsfrom additional experiments with L30T63 show that the time step of integration can be quadrupled at most inreduced grid with refinement ratio 3. In the experiment with model LICOM and with the model’s original timestep unchanged, the model covered area is extended to the whole globe from its original case with the grid pointof North Pole considered as an isolated island and the results of experiment are shown to be acceptable.展开更多
Atlantic meridional overturning circulation(AMOC)plays an important role in transporting heat meridionally in the Earth’s climate system and is also a key metrical tool to verify oceanic general circulation models.Tw...Atlantic meridional overturning circulation(AMOC)plays an important role in transporting heat meridionally in the Earth’s climate system and is also a key metrical tool to verify oceanic general circulation models.Two OMIP(Ocean Model Intercomparison Project phase 1 and 2)simulations with LICOM3(version 3 of the LASG/IAP Climate System Ocean Model)developed at the State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics(LASG),Institute of Atmospheric Physics(IAP),are compared in this study.Both simulations well reproduce the fundamental characteristics of the AMOC,but the OMIP1 simulation shows a significantly stronger AMOC than the OMIP2 simulation.Because the LICOM3 configurations are identical between these two experiments,any differences between them must be attributed to the surface forcing data.Further analysis suggests that sea surface salinity(SSS)differences should be mainly responsible for the enhanced AMOC in the OMIP1 simulation,but sea surface temperature(SST)also play an unignorable role in modulating AMOC.In the North Atlantic,where deep convection occurs,the SSS in OMIP1 is more saline than that in OMIP1.We find that in the major region of deep convection,the change of SSS has more significant effect on density than the change of SST.As a result,the SSS was more saline than that in OMIP2,leading to stronger deep convection and subsequently intensify the AMOC.We conduct a series of numerical experiments with LICOM3,and the results confirmed that the changes in SSS have more significant effect on the strength of AMOC than the changes in SST.展开更多
An eddy-permitting, quasi-global oceanic general circulation model, LICOM (LASG/IAP (State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics, Institute of Atmospheric Physi...An eddy-permitting, quasi-global oceanic general circulation model, LICOM (LASG/IAP (State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics, Institute of Atmospheric Physics) Climate System Ocean Model), with a uniform grid of 0.5? × 0.5? is established. Forced by wind stresses from Hellerman and Rosenstain (1983), a 40-yr integration is conducted with sea surface temperature and salinity being restored to the Levitus 94 datasets. The evaluation of the annual mean climatology of the LICOM control run shows that the large-scale circulation can be well reproduced. A comparison between the LICOM control run and a parallel integration of L30T63, which has the same framework but a coarse resolution, is also made to con?rm the impact of resolution on the model performance. On account of the reduction of horizontal viscosity with the enhancement of the horizontal resolution, LICOM improves the simulation with respect to not only the intensity of the large scale circulations, but also the magnitude and structure of the Equatorial Undercurrent and South Equatorial Current. Taking advantage of the ?ne grid size, the pathway of the Indonesian Through?ow (ITF) is better represented in LICOM than in L30T63. The transport of ITF in LICOM is more convergent in the upper layer. As a consequence, the Indian Ocean tends to get warmer in LICOM. The poleward heat transports for both the global and individual basins are also signi?cantly improved in LICOM. A decomposed analysis indicates that the transport due to the barotropic gyre, which primarily stands for the barotropic e?ect of the western boundary currents, plays a crucial role in making the di?erence.展开更多
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 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 circulation, 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 equatorial 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.展开更多
Different versions of a new nine-layer general circulation model which is rhomboidally truncated at zonal wavenumber 15(L9R15)are introduced in this paper.On using the observed global monthly sea surfaCe temperature(S...Different versions of a new nine-layer general circulation model which is rhomboidally truncated at zonal wavenumber 15(L9R15)are introduced in this paper.On using the observed global monthly sea surfaCe temperature(SST)and sea ice(SI)data from 1979 to 1988 offered by the internahonal Atmospheric Model Iute-comparison Program(AMIP),these different model versions were integrated for the ten-year AMIP period. Results show that the model iscapable of simulating the basic states of the atmosphere and its interannual variability,and in performing reasonablesensitivity experiments.展开更多
The dynamical framework of the nine-level version of the IAP AGCM is presented in this paper. The emphasis of the model's description is put on the following two aspects:(1) A model's standard atmosphere, whic...The dynamical framework of the nine-level version of the IAP AGCM is presented in this paper. The emphasis of the model's description is put on the following two aspects:(1) A model's standard atmosphere, which is a satisfactory approximation to the observed troposphere and lower stratosphere standard atmosphere, is introduced into the equations of the model to permit a more accurate calculation of the vertical transport terms, especially near the tropopause; (2) The vertical levels of the model are carefully selected to guarantee a smooth dependence of layer thickness upon pressure in order to reduce the truncation error involved in the unequal interval vertical finite-differencing. For testing the model, two kinds of linear baroclinic Rossby-Haurwitz waves, one of which has a dynamically stable vertical structure and the other has a relatively unstable one, are constructed to provide initial conditions for numerical experiments. The two waves have been integrated for more than 300 days and 100 days respectively by using the model and both of them are propagating westward with almost identical phase-speed during the time period of the integrations. No obvious change of the wave patterns is found at the levels in the model's troposphere. The amplitudes of both two waves at the uppermost level, however, exhibit rather significant oscillation with time, of which the periods are exactly 20 days and 25 days espectively.The explanation of this interesting phenomena is still under investigation.展开更多
This paper describes a numerical model of the world ocean based on the fully primitive equations. A 'Standard' ocean state is introduced into the equations of the model and the perturbed thermodynamic variable...This paper describes a numerical model of the world ocean based on the fully primitive equations. A 'Standard' ocean state is introduced into the equations of the model and the perturbed thermodynamic variables are used in the modlc's calculations. Both a free upper surface and a bottom topography are included in the model and a sigma coordinate is used to normalize the model's vertical component. The model has four unevenly-spaced layers and 4 × 5 horizontal resolution based on C-grid system. The finite-difference scheme of the model is designed to conserve the gross available energy in order to avoid fictitious energy generation or decay.The model has been tested in response to the annual mean surface wind stress, sea level air pressure and sea level air temperature as a preliminary step to its further improvement and its coupling with a global atmospheric general circulation model. Some of results, including currents, temperature and sea surface elevation simulated by the mode! arc presented.展开更多
The direct radiative forcing （DRF） of sulfate aerosols depends highly on the atmospheric sulfate loading and the meteorology, both of which undergo strong regional and seasonal variations. Because the optical proper...The direct radiative forcing （DRF） of sulfate aerosols depends highly on the atmospheric sulfate loading and the meteorology, both of which undergo strong regional and seasonal variations. Because the optical properties of sulfate aerosols are also sensitive to atmospheric relative humidity, in this study we first examine the scheme for optical properties that considers hydroscopic growth. Next, we investigate the seasonal and regional distributions of sulfate DRF using the sulfate loading simulated from NCAR CAM-Chem together with the meteorology modeled from a spectral atmospheric general circulation model （AGCM） developed by LASG-IAP. The global annual-mean sulfate loading of 3.44 mg m-2 is calculated to yield the DRF of -1.03 and -0.57 W m-2 for clear-sky and all-sky conditions, respectively. However, much larger values occur on regional bases. For example, the maximum all-sky sulfate DRF over Europe, East Asia, and North America can be up to -4.0 W m-2. The strongest all-sky sulfate DRF occurs in the Northern Hemispheric July, with a hemispheric average of -1.26 W m-2. The study results also indicate that the regional DRF are strongly affected by cloud and relative humidity, which vary considerably among the regions during different seasons. This certainly raises the issue that the biases in model-sinmlated regional meteorology can introduce biases into the sulfate DRF. Hence, the model processes associated with atmospheric humidity and cloud physics should be modified in great depth to improve the simulations of the LASG-IAP AGCM and to reduce the uncertainty of sulfate direct effects on global and regional climate in these simulations.展开更多
This paper investigates the processes behind the double ITCZ phenomenon, a common problem in Coupled ocean-atmosphere General Circulation Models (CGCMs), using a CGCM—FGCM-0 (Flexible General Circulat...This paper investigates the processes behind the double ITCZ phenomenon, a common problem in Coupled ocean-atmosphere General Circulation Models (CGCMs), using a CGCM—FGCM-0 (Flexible General Circulation Model, version 0). The double ITCZ mode develops rapidly during the ?rst two years of the integration and becomes a perennial phenomenon afterwards in the model. By way of Singular Value Decomposition (SVD) for SST, sea surface pressure, and sea surface wind, some air-sea interactions are analyzed. These interactions prompt the anomalous signals that appear at the beginning of the coupling to develop rapidly. There are two possible reasons, proved by sensitivity experiments: (1) the overestimated east-west gradient of SST in the equatorial Paci?c in the ocean spin-up process, and (2) the underestimated amount of low-level stratus over the Peruvian coast in CCM3 (the Community Climate Model, Version Three). The overestimated east-west gradient of SST brings the anomalous equatorial easterly. The anomalous easterly, a?ected by the Coriolis force in the Southern Hemisphere, turns into an anomalous westerly in a broad area south of the equator and is enhanced by atmospheric anomalous circulation due to the underestimated amount of low-level stratus over the Peruvian coast simulated by CCM3. The anomalous westerly leads to anomalous warm advection that makes the SST warm in the southeast Paci?c. The double ITCZ phenomenon in the CGCM is a result of a series of nonlocal and nonlinear adjustment processes in the coupled system, which can be traced to the uncoupled models, oceanic component, and atmospheric component. The zonal gradient of the equatorial SST is too large in the ocean component and the amount of low-level stratus over the Peruvian coast is too low in the atmosphere component.展开更多
Based on the theoretical spectral model of inertial internal wave breaking （fine structure） proposed previ- ously, in which the effects of the horizontal Coriolis frequency component f-tilde on a potential isopycnal...Based on the theoretical spectral model of inertial internal wave breaking （fine structure） proposed previ- ously, in which the effects of the horizontal Coriolis frequency component f-tilde on a potential isopycnal are taken into account, a parameterization scheme of vertical mixing in the stably stratified interior be- low the surface mixed layer in the ocean general circulation model （OGCM） is put forward preliminarily in this paper. Besides turbulence, the impact of sub-mesoscale oceanic processes （including inertial internal wave breaking product） on oceanic interior mixing is emphasized. We suggest that adding the inertial inter- hal wave breaking mixing scheme （F-scheme for short） put forward in this paper to the turbulence mixing scheme of Canuto et al. （T-scheme for short） in the OGCM, except the region from 15°S to 15°N. The numeri- cal results ofF-scheme by usingWOA09 data and an OGCM （LICOM, LASG/IAP climate system ocean model） over the global ocean are given. A notable improvement in the simulation of salinity and temperature over the global ocean is attained by using T-scheme adding F-scheme, especially in the mid- and high-latitude regions in the simulation of the intermediate water and deep water. We conjecture that the inertial internal wave breaking mixing and inertial forcing of wind might be one of important mechanisms maintaining the ventilation process. The modeling strength of the Atlantic meridional overturning circulation （AMOC） by using T-scheme adding F-scheme may be more reasonable than that by using T-scheme alone, though the physical processes need to be further studied, and the overflow parameterization needs to be incorporated. A shortcoming in F-scheme is that in this paper the error of simulated salinity and temperature by using T-scheme adding F-scheme is larger than that by using T-scheme alone in the subsurface layer.展开更多
In an earlier study, the Atmospheric Models Intercomparison Program (AMIP) simulations of African climate using the nine-layer gridpoint atmospheric general circulation model were found to be closely related to the ob...In an earlier study, the Atmospheric Models Intercomparison Program (AMIP) simulations of African climate using the nine-layer gridpoint atmospheric general circulation model were found to be closely related to the observed European Centre for Medium Range Weather Forecast (ECMWF) temperature data at 500 and 850 hPa. This paper presents the analysis of the simulation of African climate using the Global Ocean-Atmosphere-Land System Model (IAP/LASG GOALS) and the nine-layer spectral general circulation model rhomboidally truncated at zonal wave number 15 (L9R15) developed at the Institute of Atmospheric Physics (IAP), Chinese Academy of Sciences, Beijing. Both model simulations were not significantly different from the National Center for Environmental Prediction (NCEP) Reanalysis monthly mean data for 1980-1995 in the case of surface air temperature, sea level pressure and precipitation, with the GOALS reproducing the seasonal mean climate over Africa better. The implications of the encouraging results in developing a local area model for Nigeria have been discussed. The great role of topography in the developing of general circulation models for numerical modelling of weather and climate has been stressed.展开更多
Satellite observations of SSTs have revealed the existence of unstable waves in the equatorial eastern Pacific and Atlantic oceans. These waves have a 20-40-day periodicity with westward phase speeds of 0.4-0.6 m s^-1...Satellite observations of SSTs have revealed the existence of unstable waves in the equatorial eastern Pacific and Atlantic oceans. These waves have a 20-40-day periodicity with westward phase speeds of 0.4-0.6 m s^-1 and wavelengths of 1000-2000 km during boreal summer and fall. They are generally called tropical instability waves （TIWs）. This study investigates TIWs simulated by a high-resolution coupled atmosphere-ocean general circulation model （AOGCM）. The horizontal resolution of the model is 120 km in the atmosphere, and 30 km longitude by 20 km latitude in the ocean. Model simulations show good agreement with the observed main features associated with TIWs. The results of energetics analysis reveal that barotropic energy conversion is responsible for providing the main energy source for TIWs by extracting energy from the meridional shear of the climatological-mean equatorial currents in the mixed layer. This deeper and northward-extended wave activity appears to gain its energy through baroclinic conversion via buoyancy work, which further contributes to the asymmetric distribution of TIWs. It is estimated that the strong cooling effect induced by equatorial upwelling is partially （-30%-40%） offset by the equatorward heat flux due to TIWs in the eastern tropical Pacific during the seasons when TIWs are active. The atmospheric mixed layer just above the sea surface responds to the waves with enhanced or reduced vertical mixing. Furthermore, the changes in turbulent mixing feed back to sea surface evaporation, favoring the westward propagation of TIWs. The atmosphere to the south of the Equator also responds to TIWs in a similar way, although TIWs are much weaker south of the Equator.展开更多
A three-dimensional density field associated with mesoscaie unstable waves generated by the 3-D, primitive-equation model (Wang and Ikeda, 1996) is provided to the quasi-geostrophic pressure tendency and ω-equations,...A three-dimensional density field associated with mesoscaie unstable waves generated by the 3-D, primitive-equation model (Wang and Ikeda, 1996) is provided to the quasi-geostrophic pressure tendency and ω-equations, and to the (ageostrophic) Q-vector equation. Diagnostic analyses, analogous to the approaches in meteorology: ω-equation and Q-vector method, are for the first time developed to examine the mesoscaie dynamical processes and mechanisms of the unstable waves propagating in the mid-latitude ocean. The weaknesses and strengths of these two diagnostic approaches are evaluated and compared to the model results. The Q-vector method is then recommended to diagnose the vertical motion associated with the mesoscaie dynamics from a hydrographic CTD (conductivity-temperature-depth) array, while the quasi-geostrophic equations produce some small-scale features (errors) in the diagnosed fields.展开更多
The eight main tidal constituents have been implemented in the global ocean general circulation model with approximate 1° horizontal resolution.Compared with the observation data,the patterns of the tidal amplitu...The eight main tidal constituents have been implemented in the global ocean general circulation model with approximate 1° horizontal resolution.Compared with the observation data,the patterns of the tidal amplitudes and phases had been simulated fairly well.The responses of mean circulation,temperature and salinity are further investigated in the global sense.When implementing the tidal forcing,wind-driven circulations are reduced,especially those in coastal regions.It is also found that the upper cell transport of the Atlantic meridional overturning circulation（AMOC） reduces significantly,while its deep cell transport is slightly enhanced from 9×106m3/s to 10×106 m3/s.The changes of circulations are all related to the increase of a bottom friction and a vertical viscosity due to the tidal forcing.The temperature and salinity of the model are also significantly affected by the tidal forcing through the enhanced bottom friction,mixing and the changes in mean circulation.The largest changes occur in the coastal regions,where the water is cooled and freshened.In the open ocean,the changes are divided into three layers：cooled and freshened on the surface and below 3 000 m,and warmed and salted in the middle in the open ocean.In the upper two layers,the changes are mainly caused by the enhanced mixing,as warm and salty water sinks and cold and fresh water rises;whereas in the deep layer,the enhancement of the deep overturning circulation accounts for the cold and fresh changes in the deep ocean.展开更多
In this study, we develop a variable-grid global ocean general circulation model （OGCM） with a fine grid （1/6）° covering the area from 20°S-50°N and from 99°-150°E, and use the model to in...In this study, we develop a variable-grid global ocean general circulation model （OGCM） with a fine grid （1/6）° covering the area from 20°S-50°N and from 99°-150°E, and use the model to investigate the isopycnal surface circulation in the South China Sea （SCS）. The simulated results show four layer structures in vertical： the surface and subsurface circulation of the SCS are characterized by the monsoon driven circulation, with basin-scaled cyclonic gyre in winter and anti-cyclonic gyre in summer. The intermediate layer circulation is opposite to the upper layer, showing anti-cyclonic gyre in winter but cyclonic gyre in summer. The circulation in the deep layer is much weaker in spring and summer, with the maximum velocity speed below 0.6 cm/s. In fall and winter, the SCS deep layer circulation shows strong east boundary current along the west coast of Philippine with the velocity speed at 1.5 m/s, which flows southward in fall and northward in winter. The results have also revealed a fourlayer vertical structure of water exchange through the Luzon Strait. The dynamics of the intermediate and deep circulation are attributed to the monsoon driving and the Luzon Strait transport forcing.展开更多
An informal review is presented of recent developments in numerical simulation of the global atmospheric circulation with very fine numerical resolution models, The focus is on results obtained recently with versions ...An informal review is presented of recent developments in numerical simulation of the global atmospheric circulation with very fine numerical resolution models, The focus is on results obtained recently with versions of the GFDL SKYHI model and the Atmospheric Model for the Earth Simulator （AFES） global atmospheric models. These models have been run with effective horizontal grid resolution of -10-40 km and fine vertical resolution. The results presented demonstrate the utility of such models for the study of a diverse range of phenomena, Specifically the models are shown to simulate the development of tropical cyclones with peak winds and minimum central pressures comparable to those of the most intense hurricanes actually observed, More fundamentally, the spectrum of energy content in the mesoscale in the flow can be reproduced by these models down to near the smallest explicitly-resolved horizontal scales, in the middle atmosphere it is shown that increasing horizontal resolution can lead to significantly improved overall simulation of the global-scale circulation, The application of the models to two specific problems requiring very fine resolution global will be discussed, The spatial and temporal variability of the vertical eddy flux of zonal momentum associated with gravity waves near the tropopause is evaluated in the very fine resolution AFES model, This is a subject of great importance for understanding and modelling the flow in the middle atmosphere, Then the simulation of the small scale variations of the semidiurnal surface pressure oscillation is analyzed, and the signature of significant topographic modulation of the semidiurnal atmospheric tide is identified.展开更多
This paper describes the large scale aspects of the seasonal surface heat budget and discusses its main forcing mechanisms in the tropical Western Pacific Ocean. The high - resolution general circulation model (Semtne...This paper describes the large scale aspects of the seasonal surface heat budget and discusses its main forcing mechanisms in the tropical Western Pacific Ocean. The high - resolution general circulation model (Semtner & Chervin, 1992) used in this study reproduced well the observed upper-layer thermal structure and circulation. It is shown that at least on the average of the study region (20 °N -20°N, west boundary-160 °E) the semiannual variation is a dominant signal for all heat budget components and is presumably due to the sun's passing across the equator twice a year; but that the components have substantial differences in amplitude. The local Ekman divergence in the region does not change significantly through the year. As a result, the change in surface heat content is roughly half due to ocean ?atmosphere heat exchange and half due to heat advection by remotely forced verti-cal motion. Horizontal currents do not play a significant role directly by advection, because the wat-er which enters the region is not very much different in temperature from the water which leaves it.展开更多
In this study the structure and seasonal variations of deep mean circulation in the East/lapan Sea （E/S） were numerically simulated using a mid-resolution ocean general circulation model with two different parameter...In this study the structure and seasonal variations of deep mean circulation in the East/lapan Sea （E/S） were numerically simulated using a mid-resolution ocean general circulation model with two different parameterizations for the eddy-topography interaction （ETI）. The strong deep mean circulations observed in the EIS are well reproduced when using the ETI parameterizations. The seasonal variability in the EIS deep layer is shown by using ETI parameterization based on the potential vorticity approach, while it is not shown in the statistical dynamical parameterization. The driving mechanism of the strong deep mean currents in the E/S are discussed by investigating the effects of model grids and parameterizations. The deep mean circulation is more closely related to the baroclinic process and potential vorticity than it is to the wind driven circulation.展开更多
基金National Natural Science Foundation of China (40233031)
文摘Due to the decrease in grid size associated with the convergence of meridians toward the poles inspherical coordinates, the time steps in many global climate models with finite-difference method are restrictedto be unpleasantly small. To overcome the problem, a reduced grid is introduced to LASG/IAP world oceangeneral circulation models. The reduced grid is implemented successfully in the coarser resolutions versionmodel L30T63 at first. Then, it is carried out in the improved version model LICOM with finer resolutions. Inthe experiment with model L30T63, under time step unchanged though, execution time per single model run isshortened significantly owing to the decrease of grid number and filtering execution in high latitudes. Resultsfrom additional experiments with L30T63 show that the time step of integration can be quadrupled at most inreduced grid with refinement ratio 3. In the experiment with model LICOM and with the model’s original timestep unchanged, the model covered area is extended to the whole globe from its original case with the grid pointof North Pole considered as an isolated island and the results of experiment are shown to be acceptable.
基金Supported by the Strategic Priority Research Program of the Chinese Academy of Sciences(No.XDA19060102)the National Natural Science Foundation of China(Nos.91958201,42130608)the Strategic Priority Research Program of the Chinese Academy of Sciences(No.XDB42000000)。
文摘Atlantic meridional overturning circulation(AMOC)plays an important role in transporting heat meridionally in the Earth’s climate system and is also a key metrical tool to verify oceanic general circulation models.Two OMIP(Ocean Model Intercomparison Project phase 1 and 2)simulations with LICOM3(version 3 of the LASG/IAP Climate System Ocean Model)developed at the State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics(LASG),Institute of Atmospheric Physics(IAP),are compared in this study.Both simulations well reproduce the fundamental characteristics of the AMOC,but the OMIP1 simulation shows a significantly stronger AMOC than the OMIP2 simulation.Because the LICOM3 configurations are identical between these two experiments,any differences between them must be attributed to the surface forcing data.Further analysis suggests that sea surface salinity(SSS)differences should be mainly responsible for the enhanced AMOC in the OMIP1 simulation,but sea surface temperature(SST)also play an unignorable role in modulating AMOC.In the North Atlantic,where deep convection occurs,the SSS in OMIP1 is more saline than that in OMIP1.We find that in the major region of deep convection,the change of SSS has more significant effect on density than the change of SST.As a result,the SSS was more saline than that in OMIP2,leading to stronger deep convection and subsequently intensify the AMOC.We conduct a series of numerical experiments with LICOM3,and the results confirmed that the changes in SSS have more significant effect on the strength of AMOC than the changes in SST.
基金the Chinese Academy of Sciences "Innovation Program" (Grant No. ZKCX2-SW-210) the National Key Program for Developing Basic Sciences (Grant No.ZKCX2-SW-210)+1 种基金 the National Key Program for Developing Basic Sciences (Grant G 1999043808 , G2000078502) the National Natural Science Foundation of China under Grant Nos. 40233031 , 40231004.
文摘An eddy-permitting, quasi-global oceanic general circulation model, LICOM (LASG/IAP (State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics, Institute of Atmospheric Physics) Climate System Ocean Model), with a uniform grid of 0.5? × 0.5? is established. Forced by wind stresses from Hellerman and Rosenstain (1983), a 40-yr integration is conducted with sea surface temperature and salinity being restored to the Levitus 94 datasets. The evaluation of the annual mean climatology of the LICOM control run shows that the large-scale circulation can be well reproduced. A comparison between the LICOM control run and a parallel integration of L30T63, which has the same framework but a coarse resolution, is also made to con?rm the impact of resolution on the model performance. On account of the reduction of horizontal viscosity with the enhancement of the horizontal resolution, LICOM improves the simulation with respect to not only the intensity of the large scale circulations, but also the magnitude and structure of the Equatorial Undercurrent and South Equatorial Current. Taking advantage of the ?ne grid size, the pathway of the Indonesian Through?ow (ITF) is better represented in LICOM than in L30T63. The transport of ITF in LICOM is more convergent in the upper layer. As a consequence, the Indian Ocean tends to get warmer in LICOM. The poleward heat transports for both the global and individual basins are also signi?cantly improved in LICOM. A decomposed analysis indicates that the transport due to the barotropic gyre, which primarily stands for the barotropic e?ect of the western boundary currents, plays a crucial role in making the di?erence.
文摘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 circulation, 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 equatorial 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.
文摘Different versions of a new nine-layer general circulation model which is rhomboidally truncated at zonal wavenumber 15(L9R15)are introduced in this paper.On using the observed global monthly sea surfaCe temperature(SST)and sea ice(SI)data from 1979 to 1988 offered by the internahonal Atmospheric Model Iute-comparison Program(AMIP),these different model versions were integrated for the ten-year AMIP period. Results show that the model iscapable of simulating the basic states of the atmosphere and its interannual variability,and in performing reasonablesensitivity experiments.
文摘The dynamical framework of the nine-level version of the IAP AGCM is presented in this paper. The emphasis of the model's description is put on the following two aspects:(1) A model's standard atmosphere, which is a satisfactory approximation to the observed troposphere and lower stratosphere standard atmosphere, is introduced into the equations of the model to permit a more accurate calculation of the vertical transport terms, especially near the tropopause; (2) The vertical levels of the model are carefully selected to guarantee a smooth dependence of layer thickness upon pressure in order to reduce the truncation error involved in the unequal interval vertical finite-differencing. For testing the model, two kinds of linear baroclinic Rossby-Haurwitz waves, one of which has a dynamically stable vertical structure and the other has a relatively unstable one, are constructed to provide initial conditions for numerical experiments. The two waves have been integrated for more than 300 days and 100 days respectively by using the model and both of them are propagating westward with almost identical phase-speed during the time period of the integrations. No obvious change of the wave patterns is found at the levels in the model's troposphere. The amplitudes of both two waves at the uppermost level, however, exhibit rather significant oscillation with time, of which the periods are exactly 20 days and 25 days espectively.The explanation of this interesting phenomena is still under investigation.
文摘This paper describes a numerical model of the world ocean based on the fully primitive equations. A 'Standard' ocean state is introduced into the equations of the model and the perturbed thermodynamic variables are used in the modlc's calculations. Both a free upper surface and a bottom topography are included in the model and a sigma coordinate is used to normalize the model's vertical component. The model has four unevenly-spaced layers and 4 × 5 horizontal resolution based on C-grid system. The finite-difference scheme of the model is designed to conserve the gross available energy in order to avoid fictitious energy generation or decay.The model has been tested in response to the annual mean surface wind stress, sea level air pressure and sea level air temperature as a preliminary step to its further improvement and its coupling with a global atmospheric general circulation model. Some of results, including currents, temperature and sea surface elevation simulated by the mode! arc presented.
基金supported jointly by the grant from National Basic Research Program of China(Grant No.2012CB955303)and from the Office of Biological and Environmental Sciences,US Department of Energy
文摘The direct radiative forcing （DRF） of sulfate aerosols depends highly on the atmospheric sulfate loading and the meteorology, both of which undergo strong regional and seasonal variations. Because the optical properties of sulfate aerosols are also sensitive to atmospheric relative humidity, in this study we first examine the scheme for optical properties that considers hydroscopic growth. Next, we investigate the seasonal and regional distributions of sulfate DRF using the sulfate loading simulated from NCAR CAM-Chem together with the meteorology modeled from a spectral atmospheric general circulation model （AGCM） developed by LASG-IAP. The global annual-mean sulfate loading of 3.44 mg m-2 is calculated to yield the DRF of -1.03 and -0.57 W m-2 for clear-sky and all-sky conditions, respectively. However, much larger values occur on regional bases. For example, the maximum all-sky sulfate DRF over Europe, East Asia, and North America can be up to -4.0 W m-2. The strongest all-sky sulfate DRF occurs in the Northern Hemispheric July, with a hemispheric average of -1.26 W m-2. The study results also indicate that the regional DRF are strongly affected by cloud and relative humidity, which vary considerably among the regions during different seasons. This certainly raises the issue that the biases in model-sinmlated regional meteorology can introduce biases into the sulfate DRF. Hence, the model processes associated with atmospheric humidity and cloud physics should be modified in great depth to improve the simulations of the LASG-IAP AGCM and to reduce the uncertainty of sulfate direct effects on global and regional climate in these simulations.
基金the National Natural Science Foundation of China under Grant Nos.40221503,40231004, 40233031.
文摘This paper investigates the processes behind the double ITCZ phenomenon, a common problem in Coupled ocean-atmosphere General Circulation Models (CGCMs), using a CGCM—FGCM-0 (Flexible General Circulation Model, version 0). The double ITCZ mode develops rapidly during the ?rst two years of the integration and becomes a perennial phenomenon afterwards in the model. By way of Singular Value Decomposition (SVD) for SST, sea surface pressure, and sea surface wind, some air-sea interactions are analyzed. These interactions prompt the anomalous signals that appear at the beginning of the coupling to develop rapidly. There are two possible reasons, proved by sensitivity experiments: (1) the overestimated east-west gradient of SST in the equatorial Paci?c in the ocean spin-up process, and (2) the underestimated amount of low-level stratus over the Peruvian coast in CCM3 (the Community Climate Model, Version Three). The overestimated east-west gradient of SST brings the anomalous equatorial easterly. The anomalous easterly, a?ected by the Coriolis force in the Southern Hemisphere, turns into an anomalous westerly in a broad area south of the equator and is enhanced by atmospheric anomalous circulation due to the underestimated amount of low-level stratus over the Peruvian coast simulated by CCM3. The anomalous westerly leads to anomalous warm advection that makes the SST warm in the southeast Paci?c. The double ITCZ phenomenon in the CGCM is a result of a series of nonlocal and nonlinear adjustment processes in the coupled system, which can be traced to the uncoupled models, oceanic component, and atmospheric component. The zonal gradient of the equatorial SST is too large in the ocean component and the amount of low-level stratus over the Peruvian coast is too low in the atmosphere component.
基金The National Natural Science Foundation of China under contract No.41275084the Key Program of National Natural Science Foundation of China under contract No.41030855
文摘Based on the theoretical spectral model of inertial internal wave breaking （fine structure） proposed previ- ously, in which the effects of the horizontal Coriolis frequency component f-tilde on a potential isopycnal are taken into account, a parameterization scheme of vertical mixing in the stably stratified interior be- low the surface mixed layer in the ocean general circulation model （OGCM） is put forward preliminarily in this paper. Besides turbulence, the impact of sub-mesoscale oceanic processes （including inertial internal wave breaking product） on oceanic interior mixing is emphasized. We suggest that adding the inertial inter- hal wave breaking mixing scheme （F-scheme for short） put forward in this paper to the turbulence mixing scheme of Canuto et al. （T-scheme for short） in the OGCM, except the region from 15°S to 15°N. The numeri- cal results ofF-scheme by usingWOA09 data and an OGCM （LICOM, LASG/IAP climate system ocean model） over the global ocean are given. A notable improvement in the simulation of salinity and temperature over the global ocean is attained by using T-scheme adding F-scheme, especially in the mid- and high-latitude regions in the simulation of the intermediate water and deep water. We conjecture that the inertial internal wave breaking mixing and inertial forcing of wind might be one of important mechanisms maintaining the ventilation process. The modeling strength of the Atlantic meridional overturning circulation （AMOC） by using T-scheme adding F-scheme may be more reasonable than that by using T-scheme alone, though the physical processes need to be further studied, and the overflow parameterization needs to be incorporated. A shortcoming in F-scheme is that in this paper the error of simulated salinity and temperature by using T-scheme adding F-scheme is larger than that by using T-scheme alone in the subsurface layer.
文摘In an earlier study, the Atmospheric Models Intercomparison Program (AMIP) simulations of African climate using the nine-layer gridpoint atmospheric general circulation model were found to be closely related to the observed European Centre for Medium Range Weather Forecast (ECMWF) temperature data at 500 and 850 hPa. This paper presents the analysis of the simulation of African climate using the Global Ocean-Atmosphere-Land System Model (IAP/LASG GOALS) and the nine-layer spectral general circulation model rhomboidally truncated at zonal wave number 15 (L9R15) developed at the Institute of Atmospheric Physics (IAP), Chinese Academy of Sciences, Beijing. Both model simulations were not significantly different from the National Center for Environmental Prediction (NCEP) Reanalysis monthly mean data for 1980-1995 in the case of surface air temperature, sea level pressure and precipitation, with the GOALS reproducing the seasonal mean climate over Africa better. The implications of the encouraging results in developing a local area model for Nigeria have been discussed. The great role of topography in the developing of general circulation models for numerical modelling of weather and climate has been stressed.
基金supported by the Postdoctoral Fellow ship given by the Japan Society for the Promotion of Sciencesupported by the Kyousei and Kakushin Projects of the ministry of Education, Culture,Sports, Science, and Technology of Japan, the Core Research for Evolutional Science and Technology of the Japan Science and Technology Agencythe National Basic Research Program of China (Grant No. 2006CB403606)
文摘Satellite observations of SSTs have revealed the existence of unstable waves in the equatorial eastern Pacific and Atlantic oceans. These waves have a 20-40-day periodicity with westward phase speeds of 0.4-0.6 m s^-1 and wavelengths of 1000-2000 km during boreal summer and fall. They are generally called tropical instability waves （TIWs）. This study investigates TIWs simulated by a high-resolution coupled atmosphere-ocean general circulation model （AOGCM）. The horizontal resolution of the model is 120 km in the atmosphere, and 30 km longitude by 20 km latitude in the ocean. Model simulations show good agreement with the observed main features associated with TIWs. The results of energetics analysis reveal that barotropic energy conversion is responsible for providing the main energy source for TIWs by extracting energy from the meridional shear of the climatological-mean equatorial currents in the mixed layer. This deeper and northward-extended wave activity appears to gain its energy through baroclinic conversion via buoyancy work, which further contributes to the asymmetric distribution of TIWs. It is estimated that the strong cooling effect induced by equatorial upwelling is partially （-30%-40%） offset by the equatorward heat flux due to TIWs in the eastern tropical Pacific during the seasons when TIWs are active. The atmospheric mixed layer just above the sea surface responds to the waves with enhanced or reduced vertical mixing. Furthermore, the changes in turbulent mixing feed back to sea surface evaporation, favoring the westward propagation of TIWs. The atmosphere to the south of the Equator also responds to TIWs in a similar way, although TIWs are much weaker south of the Equator.
文摘A three-dimensional density field associated with mesoscaie unstable waves generated by the 3-D, primitive-equation model (Wang and Ikeda, 1996) is provided to the quasi-geostrophic pressure tendency and ω-equations, and to the (ageostrophic) Q-vector equation. Diagnostic analyses, analogous to the approaches in meteorology: ω-equation and Q-vector method, are for the first time developed to examine the mesoscaie dynamical processes and mechanisms of the unstable waves propagating in the mid-latitude ocean. The weaknesses and strengths of these two diagnostic approaches are evaluated and compared to the model results. The Q-vector method is then recommended to diagnose the vertical motion associated with the mesoscaie dynamics from a hydrographic CTD (conductivity-temperature-depth) array, while the quasi-geostrophic equations produce some small-scale features (errors) in the diagnosed fields.
基金The National Key Program for Developing Basic Sciences of China under contract No.2013CB956204the National Natural Science Foundation of China under contract Nos 41275084 and 41576025the Strategic Priority Research of the Chinese Academy of Science under contract Nos XDA01020304 and DA05110302
文摘The eight main tidal constituents have been implemented in the global ocean general circulation model with approximate 1° horizontal resolution.Compared with the observation data,the patterns of the tidal amplitudes and phases had been simulated fairly well.The responses of mean circulation,temperature and salinity are further investigated in the global sense.When implementing the tidal forcing,wind-driven circulations are reduced,especially those in coastal regions.It is also found that the upper cell transport of the Atlantic meridional overturning circulation（AMOC） reduces significantly,while its deep cell transport is slightly enhanced from 9×106m3/s to 10×106 m3/s.The changes of circulations are all related to the increase of a bottom friction and a vertical viscosity due to the tidal forcing.The temperature and salinity of the model are also significantly affected by the tidal forcing through the enhanced bottom friction,mixing and the changes in mean circulation.The largest changes occur in the coastal regions,where the water is cooled and freshened.In the open ocean,the changes are divided into three layers：cooled and freshened on the surface and below 3 000 m,and warmed and salted in the middle in the open ocean.In the upper two layers,the changes are mainly caused by the enhanced mixing,as warm and salty water sinks and cold and fresh water rises;whereas in the deep layer,the enhancement of the deep overturning circulation accounts for the cold and fresh changes in the deep ocean.
基金The National High Technology Research and Development Program(863 Program)of China under contract No.2013AA09A506the National Natural Science Foundation of China-Shandong Joint Fund for Marine Science Research Centers under contract No.U1406404+1 种基金the National Basic Research Program(973 Program)of China under contract No.2011CB956000the National Natural Science Foundation of China under contract No.40476016
文摘In this study, we develop a variable-grid global ocean general circulation model （OGCM） with a fine grid （1/6）° covering the area from 20°S-50°N and from 99°-150°E, and use the model to investigate the isopycnal surface circulation in the South China Sea （SCS）. The simulated results show four layer structures in vertical： the surface and subsurface circulation of the SCS are characterized by the monsoon driven circulation, with basin-scaled cyclonic gyre in winter and anti-cyclonic gyre in summer. The intermediate layer circulation is opposite to the upper layer, showing anti-cyclonic gyre in winter but cyclonic gyre in summer. The circulation in the deep layer is much weaker in spring and summer, with the maximum velocity speed below 0.6 cm/s. In fall and winter, the SCS deep layer circulation shows strong east boundary current along the west coast of Philippine with the velocity speed at 1.5 m/s, which flows southward in fall and northward in winter. The results have also revealed a fourlayer vertical structure of water exchange through the Luzon Strait. The dynamics of the intermediate and deep circulation are attributed to the monsoon driving and the Luzon Strait transport forcing.
基金National Science Foundation Grant ATM- 02-19120the Japan Agency for Marine-Earth ScienceTechnology （JAMSTEC） through its sponsorship of the International Pacific Research Center.
文摘An informal review is presented of recent developments in numerical simulation of the global atmospheric circulation with very fine numerical resolution models, The focus is on results obtained recently with versions of the GFDL SKYHI model and the Atmospheric Model for the Earth Simulator （AFES） global atmospheric models. These models have been run with effective horizontal grid resolution of -10-40 km and fine vertical resolution. The results presented demonstrate the utility of such models for the study of a diverse range of phenomena, Specifically the models are shown to simulate the development of tropical cyclones with peak winds and minimum central pressures comparable to those of the most intense hurricanes actually observed, More fundamentally, the spectrum of energy content in the mesoscale in the flow can be reproduced by these models down to near the smallest explicitly-resolved horizontal scales, in the middle atmosphere it is shown that increasing horizontal resolution can lead to significantly improved overall simulation of the global-scale circulation, The application of the models to two specific problems requiring very fine resolution global will be discussed, The spatial and temporal variability of the vertical eddy flux of zonal momentum associated with gravity waves near the tropopause is evaluated in the very fine resolution AFES model, This is a subject of great importance for understanding and modelling the flow in the middle atmosphere, Then the simulation of the small scale variations of the semidiurnal surface pressure oscillation is analyzed, and the signature of significant topographic modulation of the semidiurnal atmospheric tide is identified.
文摘This paper describes the large scale aspects of the seasonal surface heat budget and discusses its main forcing mechanisms in the tropical Western Pacific Ocean. The high - resolution general circulation model (Semtner & Chervin, 1992) used in this study reproduced well the observed upper-layer thermal structure and circulation. It is shown that at least on the average of the study region (20 °N -20°N, west boundary-160 °E) the semiannual variation is a dominant signal for all heat budget components and is presumably due to the sun's passing across the equator twice a year; but that the components have substantial differences in amplitude. The local Ekman divergence in the region does not change significantly through the year. As a result, the change in surface heat content is roughly half due to ocean ?atmosphere heat exchange and half due to heat advection by remotely forced verti-cal motion. Horizontal currents do not play a significant role directly by advection, because the wat-er which enters the region is not very much different in temperature from the water which leaves it.
基金The Research Program on Climate Change Adaptation(RECCA)of the Ministry of Education,Culture,Sports,Science and Technology(MEXT)of Japan
文摘In this study the structure and seasonal variations of deep mean circulation in the East/lapan Sea （E/S） were numerically simulated using a mid-resolution ocean general circulation model with two different parameterizations for the eddy-topography interaction （ETI）. The strong deep mean circulations observed in the EIS are well reproduced when using the ETI parameterizations. The seasonal variability in the EIS deep layer is shown by using ETI parameterization based on the potential vorticity approach, while it is not shown in the statistical dynamical parameterization. The driving mechanism of the strong deep mean currents in the E/S are discussed by investigating the effects of model grids and parameterizations. The deep mean circulation is more closely related to the baroclinic process and potential vorticity than it is to the wind driven circulation.
