Evaluation of Arctic Sea-ice Cover and Thickness Simulated by MITgcm

A regional Arctic Ocean configuration of the Massachusetts Institute of Technology General Circulation Model(MITgcm)is applied to simulate the Arctic sea ice from 1991 to 2012.The simulations are evaluated by comparing them with observations from different sources.The results show that MITgcm can reproduce the interannual and seasonal variability of the sea-ice extent,but underestimates the trend in sea-ice extent,especially in September.The ice concentration and thickness distributions are comparable to those from the observations,with most deviations within the observational uncertainties and less than 0.5 m,respectively.The simulated sea-ice extents are better correlated with observations in September,with a correlation coefficient of 0.95,than in March,with a correlation coefficient of 0.83.However,the distributions of sea-ice concentration are better simulated in March,with higher pattern correlation coefficients(0.98)than in September.When the model underestimates the atmospheric influence on the sea-ice evolution in March,deviations in the sea-ice concentration arise at the ice edges and are higher than those in September.In contrast,when the model underestimates the oceanic boundaries’influence on the September sea-ice evolution,disagreements in the distribution of the sea-ice concentration and its trend are found over most marginal seas in the Arctic Ocean.The uncertainties of the model,whereby it fails to incorporate the atmospheric information in March and oceanic information in September,contribute to varying model errors with the seasons.     

A numerical world ocean general circulation model

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 modle’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 model are presented.     

Cloud and Radiation Processes Simulated by a Coupled Atmosphere-Ocean Model

Using NCC/IAP T63 coupled atmosphere-ocean general circulation model (AOGCM),two 20-yr integra- tions were processed,and their ability to simulate cloud and radiation was analysed in detail.The results show that the model can simulate the basic distribution of cloud cover,and however,obvious differences still exist compared with ISCCP satellite data and ERA reanalysis data.The simulated cloud cover is less in general,especially the abnormal low values in some regions of ocean.By improving the cloud cover scheme, simulated cloud cover in the eastern Pacific and Atlantic,summer hemisphere's oceans from subtropical to mid-latitude is considerably improved.But in the tropical Indian Ocean and West Pacific the cloud cover difference is still evident,mainly due to the deficiency of high cloud simulation in these regions resulting from deep cumulus convection.In terms of the analysis on radiation and cloud radiative forcing,we find that simulation on long wave radiation is better than short wave radiation.The simulation error of short wave radiation is caused mostly by the simulation difference in short wave radiative forcing,sea ice,and snow cover,and also by not involving aerosol's effect.The simulation error of long wave radiation is mainly resulting from deficiency in simulating cloud cover and underlying surface temperature.Corresponding to improvement of cloud cover,the simulated radiation (especially short wave radiation) in eastern oceans, summer hemisphere's oceans from subtropical to mid-latitude is remarkably improved.This also brings obvious improvement to net radiation in these regions.     

Simulating the effects of the 1991 Mount Pinatubo volcanic eruption using the ARPEGE atmosphere general circulation model

The climate changes that occured following the volcanic eruption of Mount Pinatubo in the Phillippines on 15 June 1991 have been simulated using the ARPEGE atmosphere general circulation model （AGCM）. The model was forced by a reconstructed spatial-time distribution of stratospheric aerosols intended for use in long climate simulations. Four statistical ensembles of the AGCM simulations with and without volcanic aerosols over a period of 5 years following the eruption have been made, and the calculated fields have been compared to available observations. The model is able to reproduce some of the observed features after the eruption, such as the winter warming pattern that was observed over the Northern Hemisphere （NH） during the following winters. This pattern was caused by an enhanced Equator-to-pole temperature gradient in the stratosphere that developed due to aerosol heating of the tropics. This in turn led to a strengthening of the polar vortex, which tends to modulate the planetary wave field in such a way that an anomalously positive Arctic Oscillation pattern is produced in the troposphere and at the surface, favouring warm conditions over the NH. During the summer, the model produced a more uniform cooling over the NH.     

A Coupled General Circulation Model for the Tropical Pacific Ocean and Global Atmosphere

On the basis of Zeng’s theoretical design, a coupled general circulation model (CGCM) is developed with its characteristics different from other CGCMs such as the unified vertical coordinates and subtraction of the standard stratification for both atmosphere and ocean, available energy consideration, and so on. The oceanic component is a free surface tropical Pacific Ocean GCM between 30oN and 30oS with horizontal grid spacing of 1o in latitude and 2o in longitude, and with 14 vertical layers. The atmospheric component it a global GCM with low-resolution of 4o in latitude and 5o in longitude, and two layers or equal man in the vertical between the surface and 200 hPa. The atmospheric GCM includes comprehensive physical processes. The coupled model is subjected to seasonally-varying cycle. Several coupling experiments, ranging from straight forward coupling without flux correction to one with flux correction, and to so-called predictor-corrector monthly coupling (PCMC), are conducted to show the existence and final controlling of the climate drift in the coupled system. After removing the climate drift with the PCMC scheme, the coupled model is integrated for more than twenty years. The results show reasonable simulations of the annual mean and its seasonal cycle of the atmospheric and oceanic circulation. The model also produces the coherent interannual variations of the climate system, manifesting the observed El Ni?o / Southern Oscillation (ENSO).     

CAS-ESM2.0 Model Datasets for the CMIP6 Ocean Model Intercomparison Project Phase 1(OMIP1)

As a member of the Chinese modeling groups,the coupled ocean-ice component of the Chinese Academy of Sciences’Earth System Model,version 2.0(CAS-ESM2.0),is taking part in the Ocean Model Intercomparison Project Phase 1(OMIP1)experiment of phase 6 of the Coupled Model Intercomparison Project(CMIP6).The simulation was conducted,and monthly outputs have been published on the ESGF(Earth System Grid Federation)data server.In this paper,the experimental dataset is introduced,and the preliminary performances of the ocean model in simulating the global ocean temperature,salinity,sea surface temperature,sea surface salinity,sea surface height,sea ice,and Atlantic Meridional Overturning Circulation(AMOC)are evaluated.The results show that the model is at quasi-equilibrium during the integration of 372 years,and performances of the model are reasonable compared with observations.This dataset is ready to be downloaded and used by the community in related research,e.g.,multi-ocean-sea-ice model performance evaluation and interannual variation in oceans driven by prescribed atmospheric forcing.     

A Numerical World Ocean General Circulation Model

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.     

NUMERICAL SIMULATIONS OF OCEANIC ELEMENTS IN THE SCS AND ITS NEIGHBORING SEA REGIONS FROM JANUARY TO AUGUST IN 1998

Studies on oceanic conditions in the South China Sea (SCS) and adjacent waters are helpful for thorough understanding of summer monsoons in East Asia. To have a 3-dimensional picture of how the oceanic currents vary, the oceanic elements in the South China Sea (SCS) and its neighboring sea regions in January～August 1998 have been simulated by using the improved Princeton University Ocean Model (POM) in this paper. The main results are in good agreement with that of ocean investigations and other simulations. The results show that the SCS branch of the Kuroshio Current is an important part in the north SCS from January to August; the SCS warm current is reproduced clearly in all months except in winter; there always exists a large-scale anti-cyclonic vortex on the right of the Kuroshio Current from January to August. In the model domain, the surface currents of the SCS have the closest relations with the monsoon with an apparent seasonal variation. In addition, the developing characteristics of the SST in the SCS and its neighboring sea regions before and after the summer monsoon onset are also well simulated by the improved POM. Those are the foundation for developing a coupled regional ocean-atmospheric model system.     

THE IMPACT OF HORIZONTAL RESOLUTION ON THE INTENSITY AND MICROSTRUCTURE OF SUPER TYPHOON USAGI

Typhoon Usagi (1319) was simulated by using the Advanced Weather Research and Forecasting numerical model (WRF) with different horizontal resolution to understand the impact of horizontal resolution on the intensity and characteristics of typhoon’s microstructures (including dynamic and microphysical structure). The simulated results show that the improvement of horizontal resolution from 5 km to 1 km has little impact on the track which is comparable to real results, but has a significant impact on the intensity and microstructures, and especially, the impact on wind speed at 10 m height, the vertical movement and precipitation intensity is the greatest. When the resolution is increased to 1 km, the intensity and characteristics of typhoon’s microstructures can be simulated better. In lower resolution simulations, some structural characteristics, including more asymmetrical and more outward tilted eyewall, and less water vapor flux on sea surface, work together to weaken typhoon intensity.     

A steady two-dimensional climate model with residual circulation

Based on dynamical energy transport and thermodynamic energy balance in the earth’s atmosphere-ocean system a steady two-dimensional climate model with residual circulation is proposed. In the model, we include some important physical processes with feedbacks such as ice caps-albedo, water vapor-tempera-ture, etc. The simulated steady temperature field is very close to that of the real atmosphere. The numerical experiments show that doubling of the atmospheric carbon dioxide results in temperature increase of 1~2oC at the low latitude surface and 6~8oC at the high latitude surface. It is shown that a 6% decrease in the solar constant is required for the -10oC ice edge to move from its present latitude ~70o to~50o.     

Sensitivity of Arctic warming to sea surface temperature distribution over melted sea-ice region in atmospheric general circulation model experiments

Substantial reduction in Arctic sea ice in recent decades has intensified air-sea interaction over the Arctic Ocean and has altered atmospheric states in the Arctic and surrounding high-latitude regions. This study has found that the atmospheric responses related to Arctic sea-ice melt in the cold season (October–March) depend on sea-ice fraction and are very sensitive to in situ sea surface temperature (SST) from a series of atmospheric general circulation model (AGCM) simulations in which multiple combinations of SSTs and sea-ice concentrations are prescribed in the Arctic Ocean. It has been found that the amplitude of surface warming over the melted sea-ice region is controlled by concurrent in situ SST even if these simulations are forced by the same sea-ice concentration. Much of the sensitivity of surface warming to in situ SST are related with large changes in surface heat fluxes such as the outgoing long-wave flux in early winter (October–December) and the sensible and latent heat fluxes for the entire cold season. Vertical extension of surface warming and moistening is sensitive to these changes as well; the associated condensational heating modulates a static stability in the lower troposphere. This study also indicates that changes in SST fields in AGCM simulations must be implemented with extra care, especially in the melted sea-ice region in the Arctic. The statistical method introduced in this study for adjusting SSTs in conjunction with a given sea-ice change can help to model the atmospheric response to sea-ice loss more accurately.     

Heat budget of the upper Arctic Ocean under a warming climate

The heat budget of the upper Arctic Ocean is examined in an ensemble of coupled climate models under idealised increasing CO₂ scenarios. All of the experiments show a strong amplification of surface air temperatures but a smaller increase in sea surface temperature than the rest of the world as heat is lost to the atmosphere as the sea-ice cover is reduced. We carry out a heat budget analysis of the Arctic Ocean in an ensemble of model runs to understand the changes that occur as the Arctic becomes ice free in summer. We find that as sea-ice retreats heat is lost from the ocean surface to the atmosphere contributing to the amplification of Arctic surface temperatures. Furthermore, heat is mixed upwards into the mixed layer as a result of increased upper ocean mixing and there is increased advection of heat into the Arctic as the ice edge retreats. Heat lost from the upper Arctic Ocean to the atmosphere is therefore replenished by mixing of warmer water from below and by increased advection of warm water from lower latitudes. The ocean is therefore able to contribute more to Arctic amplification.     

Simulation of the interannual variability of the wind-driven Arctic sea-ice cover during 1958–1998

A thermodynamic-dynamic sea-ice model based on a granular material rheology developed by Tremblay and Mysak is used to study the interannual variability of the Arctic sea-ice cover during the 41-year period 1958–98. Monthly wind stress forcing derived from the National Centers for Environmental Prediction (NCEP) Reanalysis data is used to produce the year-to-year variations in the sea-ice circulation and thickness. We focus on analyzing the variability of the sea-ice volume in the Arctic Basin and the subsequent changes in sea-ice export into the Greenland Sea via Fram Strait. The relative contributions of the Fram Strait sea-ice thickness and velocity anomalies to the sea-ice export anomalies are first investigated, and the former is shown to be particularly important during several large export events. The sea-ice export anomalies for these events are next linked to prior sea-ice volume anomalies in the Arctic Basin. The origin and evolution of the sea-ice volume anomalies are then related to the sea-ice circulation and atmospheric forcing patterns in the Arctic. Large sea-ice export anomalies are generally preceded by large volume anomalies formed along the East Siberian coast due to anomalous winds which occur when the Arctic High is centered closer than usual to this coastal area. When the center of this High relocates over the Beaufort Sea and the Icelandic Low extends far into the Arctic Basin, the ice volume anomalies are transported to the Fram Strait region via the Transpolar Drift Stream. Finally, the link between the sea-ice export and the North Atlantic Oscillation (NAO) index is briefly discussed. The overall results from this study show that the Arctic Basin and its ice volume anomalies must be considered in order to fully understand the export through Fram Strait. Received: 27 January 1999 / Accepted: 8 July 1999     

Actual and insolation-weighted Northern Hemisphere snow cover and sea-ice between 1973–2002

Actual and insolation-weighted Northern Hemisphere snow cover and sea ice are binned by latitude bands for the years 1973–2002. Antarctic sea-ice is also analyzed for the years 1980–2002. The use of insolation weighting provides an improved estimate of the radiative feedbacks of snow cover and sea-ice into the atmosphere. One conclusion of our assessment is that while a decrease in both areal and insolation-weighted values have occurred, the data does not show a monotonic decrease of either Arctic sea-ice or Northern Hemisphere snow cover. If Arctic perennial sea-ice is decreasing since the total reduction in areal coverage is relatively small, a large portion of it is being replenished each year such that its radiative feedback to the atmosphere is muted. Antarctic sea-ice areal cover shows no significant long-term trend, while there is a slight decrease in the insolation-weighted values for the period 1980–2002. From the early 1990s to 2001, there was a slight increase in both values. The comparison of general circulation model simulations of changes over the last several decades to observed changes in insolation-weighted sea-ice and snow cover should be a priority research topic.     

Time-dependent response of a zonally averaged ocean–atmosphere–sea ice model to Milankovitch forcing

An ocean–atmosphere–sea ice model is developed to explore the time-dependent response of climate to Milankovitch forcing for the time interval 5–3 Myr BP. The ocean component is a zonally averaged model of the circulation in five basins (Arctic, Atlantic, Indian, Pacific, and Southern Oceans). The atmospheric component is a one-dimensional (latitudinal) energy balance model, and the sea-ice component is a thermodynamic model. Two numerical experiments are conducted. The first experiment does not include sea ice and the Arctic Ocean; the second experiment does. Results from the two experiments are used to investigate (1) the response of annual mean surface air and ocean temperatures to Milankovitch forcing, and (2) the role of sea ice in this response. In both experiments, the response of air temperature is dominated by obliquity cycles at most latitudes. On the other hand, the response of ocean temperature varies with latitude and depth. Deep water formed between 45°N and 65°N in the Atlantic Ocean mainly responds to precession. In contrast, deep water formed south of 60°S responds to obliquity when sea ice is not included. Sea ice acts as a time-integrator of summer insolation changes such that annual mean sea-ice conditions mainly respond to obliquity. Thus, in the presence of sea ice, air temperature changes over the sea ice are amplified, and temperature changes in deep water of southern origin are suppressed since water below sea ice is kept near the freezing point.     

A model-based study of ice and freshwater transport variability along both sides of Greenland

We investigate some aspects of the variability of the Arctic freshwater content during the 1965–2002 period using the DRAKKAR eddy admitting global ocean/sea-ice model (12 km resolution in the Arctic). A comparison with recent mooring sections shows that the model realistically represents the major advective exchanges with the Arctic basin, through Bering, Fram and Davis Straits, and the Barents Sea. This allows the separate contributions of the inflows and outflows across each section to be quantified. In the model, the Arctic freshwater content variability is explained by the sea-ice flux at Fram and the combined variations of ocean freshwater inflow (at Bering) and outflow (at Fram and Davis). At all routes, except trough Fram Strait, the freshwater transport variability is mainly accounted for by the liquid component, with small contributions from the sea-ice flux. The ocean freshwater transport variability through both Davis and Fram is controlled by the variability of the export branch (Baffin Island Current and East Greenland Current, respectively), the variability of the inflow branches playing a minor role. We examine the respective role of velocity and salinity fluctuations in the variability of the ocean freshwater transport. Fram and Davis Straits offer a striking contrast in this regard. Freshwater transport variations across Davis Strait are completely determined by the variations of the total volume flux (0.91 correlation). On the other hand, the freshwater transport through Fram Strait depends both on variations of volume transport and salinity. As a result, there is no significant correlation between the variability of freshwater flux at Fram and Davis, although the volume transports on each side of Greenland are strongly anti-correlated (−0.84). Contrary to Davis Strait, the salinity of water carried by the East Greenland Current through Fram Strait varies strongly due to the ice-ocean flux north of Greenland.     

Impact of a projected future Arctic Sea Ice reduction on extratropical storminess and the NAO

The impact of a reduced Arctic sea ice cover on wintertime extratropical storminess is investigated by conducting atmospheric general circulation model (AGCM) experiments. The AGCM ECHAM5 is forced by the present and a projected future seasonal cycle of Arctic sea ice. In the experiment with projected sea-ice concentrations significant reductions in storminess were found during December and January in both midlatitudes and towards the Arctic. However, a substantially larger reduction in extratropical storminess was found in March, despite a smaller change in surface energy fluxes in March than in the other winter months. The projected decrease in storminess is also related to the negative phase of the North Atlantic Oscillation (NAO). The March response is consistent with a forcing from transient and quasi-stationary eddies associated with negative NAO events. The greater sensitivity to sea-ice anomalies in late winter sets this study apart from earlier ones.     

The role of terrestrial snow cover in the climate system

Snow cover is known to exert a strong influence on climate, but quantifying its impact is difficult. This study investigates the global impact of terrestrial snow cover through a pair of GCM simulations run with prognostic snow cover and with all snow cover on land eliminated (NOSNOWCOVER). In this experiment all snowfall over land was converted into its liquid–water equivalent upon reaching the surface. Compared with the control run, NOSNOWCOVER produces mean-annual surface air temperatures up to 5 K higher over northern North America and Eurasia and 8–10 K greater during winter. The globally averaged warming of 0.8 K is one-third as large as the model’s response to 2 × CO₂ forcing. The pronounced surface heating propagates throughout the troposphere, causing changes in surface and upper-air circulation patterns. Despite the large atmospheric warming, the absence of an insulating snow pack causes soil temperatures in NOSNOWCOVER to fall throughout northern Asia and Canada, including extreme wintertime cooling of over 20 K in Siberia and a 70% increase in permafrost area. The absence of snow melt water also affects extratropical surface hydrology, causing significantly drier upper-layer soils and dramatic changes in the annual cycle of runoff. Removing snow cover also drastically affects extreme weather. Extreme cold-air outbreaks (CAOs)—defined relative to the control climatology—essentially disappear in NOSNOWCOVER. The loss of CAOs appears to stem from both the local effects of eliminating snow cover in mid-latitudes and a remote effect over source regions in the Arctic, where −40°C air masses are no longer able to form.     

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

本文发展了一个可以适用于任意水平正交曲线坐标系的海洋模式动力框架,并将其应用于中国科学院大气物理研究所大气科学和地球流体力学数值模拟国家重点实验室发展的气候系统海洋模式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左右。