A model of late Pleistocene ice sheet growth with realistic geography and simplified cryodynamics and geodynamics

A global two-dimensional one-level seasonal energy-balance model is asynchronously coupled to vertically integrated ice-flow models (which depend both on latitude and longitude) to study the response of the atmosphere-ocean-cryosphere-lithosphere system to solar forcing for the last ice age cycle of the late Pleistocene. The model simulates the position of the North American and European ice sheet complexes at the last glacial maximum satisfactorily. Both the geographic distributions of the ice volumes delivered by the model and their masses are a reasonable approximation to those inferred on the basis of relative sea level data (Tushingham and Peltier 1990). The sensitivity of the coupled model over the last glacial-interglacial cycle to solar forcing is nevertheless low, which suggests that further physical mechanisms will have to be added to the model (such as explicit basal sliding and ice shelves which would respond to sea-level variations and therefore permit marine incursions), if it is to adequately simulate the terminations that control the 10⁵ year ice age cycle. One should also incorporate long-term variations of the greenhouse gases (Manabe et al. 1985b).This paper was presented at the International Conference on Modelling of Global Climate Change and Variability, held in Hamburg 11–15 September 1989 under the auspices of the Meteorological Institute of the University of Hamburg and the Max Planck Institute for Meteorology. Guest Editor for these papers is Dr. L. Dümenil     

The effects of aggressive mitigation on steric sea level rise and sea ice changes

With an increasing political focus on limiting global warming to less than 2 °C above pre-industrial levels it is vital to understand the consequences of these targets on key parts of the climate system. Here, we focus on changes in sea level and sea ice, comparing twenty-first century projections with increased greenhouse gas concentrations (using the mid-range IPCC A1B emissions scenario) with those under a mitigation scenario with large reductions in emissions (the E1 scenario). At the end of the twenty-first century, the global mean steric sea level rise is reduced by about a third in the mitigation scenario compared with the A1B scenario. Changes in surface air temperature are found to be poorly correlated with steric sea level changes. While the projected decreases in sea ice extent during the first half of the twenty-first century are independent of the season or scenario, especially in the Arctic, the seasonal cycle of sea ice extent is amplified. By the end of the century the Arctic becomes sea ice free in September in the A1B scenario in most models. In the mitigation scenario the ice does not disappear in the majority of models, but is reduced by 42 % of the present September extent. Results for Antarctic sea ice changes reveal large initial biases in the models and a significant correlation between projected changes and the initial extent. This latter result highlights the necessity for further refinements in Antarctic sea ice modelling for more reliable projections of future sea ice.     

中国地区夏季云粒子尺寸的时空分布特征

利用CloudSat卫星资料,分析了2006～2008年中国地区夏季月平均云粒子有效半径的垂直和区域变化特征。结果显示,水云粒子有效半径在对流层低层达到最大,并随高度增加而减小。30°N纬度带的水云相对以南及以北纬度带的粒子有效半径偏大。6月水云粒子有效半径较大,应与梅雨季节有密切联系。对于中国北部和中部,水云粒子有效半径在西部较东部偏大,而在南部地区,东西部差异不明显。不同纬度带上的冰云粒子有效半径相类似,在冰云下边界最大,随高度增加而减小。水云和冰云的云粒子尺度的年际变化不明显。对上述特征的成因分析表明,高原地形以及东亚夏季风对月平均云粒子有效半径具有明显影响。所揭示的云粒子有效半径特征为天气和气候模式改进、人工影响天气及云—辐射—气候相互作用等研究提供了重要的基础信息。     

Greenland Ice Sheet Contribution to Future Global Sea Level Rise based on CMIP5 Models简

Sea level rise （SLR） is one of the major socioeconomic risks associated with global warming. Mass losses from the Greenland ice sheet （GrIS） will be partially responsible for future SLR, although there are large uncertainties in modeled climate and ice sheet behavior. We used the ice sheet model SICOPOLIS （Simulation COde for POLythermal Ice Sheets） driven by climate projections from 20 models in the fifth phase of the Coupled Model Intercomparison Project （CMIP5） to estimate the GrlS contribution to global SLR. Based on the outputs of the 20 models, it is estimated that the GrIS will contribute 0-16 （0-27） cm to global SLR by 2100 under the Representative Concentration Pathways （RCP） 4.5 （RCP 8.5） scenarios. The projected SLR increases further to 7-22 （7-33） cm with 2~basal sliding included. In response to the results of the multimodel ensemble mean, the ice sheet model projects a global SLR of 3 cm and 7 cm （10 cm and 13 cm with 2~basal sliding） under the RCP 4.5 and RCP 8.5 scenarios, respectively. In addition, our results suggest that the uncertainty in future sea level projection caused by the large spread in climate projections could be reduced with model-evaluation and the selective use of model outputs.     

Atmospheric Precursors of and Response to Anomalous Arctic Sea Ice in CMIP5 Models

This study examines pre-industrial control simulations from CMIP5 climate models in an effort to better understand the complex relationships between Arctic sea ice and the stratosphere, and between Arctic sea ice and cold winter temperatures over Eurasia. We present normalized regressions of Arctic sea-ice area against several atmospheric variables at extended lead and lag times. Statistically significant regressions are found at leads and lags, suggesting both atmospheric precursors of, and responses to, low sea ice; but generally, the regressions are stronger when the atmosphere leads sea ice, including a weaker polar stratospheric vortex indicated by positive polar cap height anomalies. Significant positive midlatitude eddy heat flux anomalies are also found to precede low sea ice. We argue that low sea ice and raised polar cap height are both a response to this enhanced midlatitude eddy heat flux. The so-called "warm Arctic, cold continents" anomaly pattern is present one to two months before low sea ice, but is absent in the months following low sea ice, suggesting that the Eurasian cooling and low sea ice are driven by similar processes. Lastly, our results suggest a dependence on the geographic region of low sea ice, with low Barents–Kara Sea ice correlated with a weakened polar stratospheric vortex, whilst low Sea of Okhotsk ice is correlated with a strengthened polar vortex. Overall, the results support a notion that the sea ice, polar stratospheric vortex and Eurasian surface temperatures collectively respond to large-scale changes in tropospheric circulation.     

Cirrus cloud horizontal and vertical inhomogeneity effects in a GCM

Summary A set of the inhomogeneity factor for high-level clouds derived from the ISCCP D1 dataset averaged over a five-year period has been incorporated in the UCLA atmospheric GCM to investigate the effect of cirrus cloud inhomogeneity on climate simulation. The inclusion of this inhomogeneous factor improves the global mean planetary albedo by about 4% simulated from the model. It also produces changes in solar fluxes and OLRs associated with changes in cloud fields, revealing that the cloud inhomogeneity not only affects cloud albedo directly, but also modifies cloud and radiation fields. The corresponding difference in the geographic distribution of precipitation is as large as 7 mm day⁻¹. Using the climatology cloud inhomogeneity factor also produces a warmer troposphere related to changes in the cloudiness and the corresponding radiative heating, which, to some extent, corrects the cold bias in the UCLA AGCM. The region around 14 km, however, is cooler associated with increase in the reflected solar flux that leads to a warmer region above. An interactive parameterization for mean effective ice crystal size based on ice water content and temperature has also been developed and incorporated in the UCLA AGCM. The inclusion of the new parameterization produces substantial differences in the zonal mean temperature and the geographic distribution of precipitation, radiative fluxes, and cloud cover with respect to the control run. The vertical distribution of ice crystal size appears to be an important factor controlling the radiative heating rate and the consequence of circulation patterns, and hence must be included in the cloud-radiation parameterization in climate models to account for realistic cloud processes in the atmosphere.     

Evaluation of Arctic Sea Ice Drift and its Relationship with Near-surface Wind and Ocean Current in Nine CMIP6 Models from China

The simulated Arctic sea ice drift and its relationship with the near-surface wind and surface ocean current during 1979-2014 in nine models from China that participated in the sixth phase of the Coupled Model Intercomparison Project(CMIP6)are examined by comparison with observational and reanalysis datasets.Most of the models reasonably represent the Beaufort Gyre(BG)and Transpolar Drift Stream(TDS)in the spatial patterns of their long-term mean sea ice drift,while the detailed location,extent,and strength of the BG and TDS vary among the models.About two-thirds of the models agree with the observation/reanalysis in the sense that the sea ice drift pattern is consistent with the near-surface wind pattern.About the same proportion of models shows that the sea ice drift pattern is consistent with the surface ocean current pattern.In the observation/reanalysis,however,the sea ice drift pattern does not match well with the surface ocean current pattern.All nine models missed the observational widespread sea ice drift speed acceleration across the Arctic.For the Arctic basin-wide spatial average,five of the nine models overestimate the Arctic long-term(1979-2014)mean sea ice drift speed in all months.Only FGOALS-g3 captures a significant sea ice drift speed increase from 1979 to 2014 both in spring and autumn.The increases are weaker than those in the observation.This evaluation helps assess the performance of the Arctic sea ice drift simulations in these CMIP6 models from China.     

渤海海冰数值模式的研究进展

根据渤海海冰冰情,从热力和动力两个方面综述了渤海海冰数值模式的研究进展。从热力参数方案和气候特点上,展望新一代渤海中小尺度短期海冰热力模式;从海冰动力模式中3种常用数值方法的主要特点和实际应用情况上,论述SPH(Smoothed Particle Hydrodynamics)方法在渤海海冰短期预报的发展前景。在此基础上,结合海洋气象的发展,讨论海冰短期预报的关键技术。将海冰数值预报模式与传统的中长期海冰统计预报有机地结合起来将是渤海海冰模式未来的发展方向。     

A simulation of sea‐ice motion and distribution off Newfoundland during LIMEX,March 1987

Abstract

A coupled ice‐ocean dynamical model is applied to the simulation of sea‐ice motion and distribution off Newfoundland during the Labrador Ice Margin Experiment (LIMEX), March 1987. In the model, the ice is coupled to a barotropic ocean through an Ekman layer that deepens with increasing wind speed. A 6‐hourly gridded wind dataset was used as input to drive the ice and the ocean. The results show that ice velocities with ice‐ocean coupling are appreciably higher than those without coupling because of the generation of wind‐driven coastal currents. This suggests that coupled ice‐ocean dynamics should always be considered in short‐term sea‐ice models. The model gives reasonable agreement with the observed ice edge except in the southern boundary where ice‐melt has a strong influence on the ice‐edge position. Ocean currents, sea level and ice velocities computed from the model are in qualitative agreement with limited current‐meter, tide‐gauge, and ice drifter trajectory data.     

A regional model for studies of atmosphere-ice-ocean interaction in the western Arctic

Summary Asa step in the development of a fully coupled regional model of the atmosphere-ice-ocean system, atmospheric and sea ice models have been adapted to a western Arctic domain centered on the Bering Strait. Lateral boundary conditions derived from operational analyses drive the models through simulations on grids having horizontal resolutions of 21 km and 7 km. Sensitivities to the presence of sea ice are large after only 48 hours, by which time the surface temperatures in the Bering and Chukchi Seas are 10–15°C higher without sea ice than with sea ice. The temperatures, in turn, modify the fields of sea level pressure, surface wind and precipitation. By influencing the surface wind stress through the static static stability, the surface state feeds back to the surface momentum exchange, ice/ocean transport, and the rate of formation of new ice. The results also show a resolution-dependence of the surface winds, precipitation rates and new ice formation rates, particularly in areas in which the coastal configuration and topography are spatially complex. The experiments will be augmented by the implementation of an ocean model on the same grids.With 12 Figures     

On the radiative properties of ice clouds: Light scattering,remote sensing,and radiation parameterization

Presented is a review of the radiative properties of ice clouds from three perspectives: light scattering simulations, remote sensing applications, and broadband radiation parameterizations appropriate for numerical models. On the subject of light scattering simulations, several classical computational approaches are reviewed, including the conventional geometric-optics method and its improved forms, the finite-difference time domain technique, the pseudo-spectral time domain technique, the discrete dipole approximation method, and the T-matrix method, with specific applications to the computation of the singlescattering properties of individual ice crystals. The strengths and weaknesses associated with each approach are discussed.With reference to remote sensing, operational retrieval algorithms are reviewed for retrieving cloud optical depth and effective particle size based on solar or thermal infrared(IR) bands. To illustrate the performance of the current solar- and IR-based retrievals, two case studies are presented based on spaceborne observations. The need for a more realistic ice cloud optical model to obtain spectrally consistent retrievals is demonstrated. Furthermore, to complement ice cloud property studies based on passive radiometric measurements, the advantage of incorporating lidar and/or polarimetric measurements is discussed.The performance of ice cloud models based on the use of different ice habits to represent ice particles is illustrated by comparing model results with satellite observations. A summary is provided of a number of parameterization schemes for ice cloud radiative properties that were developed for application to broadband radiative transfer submodels within general circulation models(GCMs). The availability of the single-scattering properties of complex ice habits has led to more accurate radiation parameterizations. In conclusion, the importance of using nonspherical ice particle models in GCM simulations for climate studies is proven.     

Twentieth century simulation of the southern hemisphere climate in coupled models. Part II: sea ice conditions and variability

We examine the representation of the mean state and interannual variability of Antarctic sea ice in six simulations of the twentieth century from coupled models participating in the Intergovernmental Panel on Climate Change fourth assessment report. The simulations exhibit a largely seasonal southern hemisphere ice cover, as observed. There is a considerable scatter in the monthly simulated climatological ice extent among different models, but no consistent bias when compared to observations. The scatter in maximum winter ice extent among different models is correlated to the strength of the climatological zonal winds suggesting that wind forced ice transport is responsible for much of this scatter. Observations show that the leading mode of southern hemisphere ice variability exhibits a dipole structure with anomalies of one sign in the Atlantic sector associated with anomalies of the opposite sign in the Pacific sector. The observed ice anomalies also exhibit eastward propagation with the Antarctic circumpolar current, as part of the documented Antarctic circumpolar wave phenomenon. Many of the models do simulate dipole-like behavior in sea ice anomalies as the leading mode of ice variability, but there is a large discrepancy in the eastward propagation of these anomalies among the different models. Consistent with observations, the simulated Antarctic dipole-like variations in the ice cover are led by sea-level pressure anomalies in the Amundsen/ Bellingshausen Sea. These are associated, to different degrees in different models, with both the southern annular mode and the El Nino-Southern Oscillation (ENSO). There are indications that the magnitude of the influence of ENSO on the southern hemisphere ice cover is related to the strength of ENSO events simulated by the different models.     

The sea ice mass budget of the Arctic and its future change as simulated by coupled climate models

Arctic sea ice mass budgets for the twentieth century and projected changes through the twenty-first century are assessed from 14 coupled global climate models. Large inter-model scatter in contemporary mass budgets is strongly related to variations in absorbed solar radiation, due in large part to differences in the surface albedo simulation. Over the twenty-first century, all models simulate a decrease in ice volume resulting from increased annual net melt (melt minus growth), partially compensated by reduced transport to lower latitudes. Despite this general agreement, the models vary considerably regarding the magnitude of ice volume loss and the relative roles of changing melt and growth in driving it. Projected changes in sea ice mass budgets depend in part on the initial (mid twentieth century) ice conditions; models with thicker initial ice generally exhibit larger volume losses. Pointing to the importance of evolving surface albedo and cloud properties, inter-model scatter in changing net ice melt is significantly related to changes in downwelling longwave and absorbed shortwave radiation. These factors, along with the simulated mean and spatial distribution of ice thickness, contribute to a large inter-model scatter in the projected onset of seasonally ice-free conditions.     

A continuous simulation of global ice volume over the past 1 million years with 3-D ice-sheet models

Sea-level records show large glacial-interglacial changes over the past million years, which on these time scales are related to changes of ice volume on land. During the Pleistocene, sea-level changes induced by ice volume are largely caused by the waxing and waning of the large ice sheets in the Northern Hemisphere. However, the individual contributions of ice in the Northern and Southern Hemisphere are poorly constrained. In this study, for the first time a fully coupled system of four 3-D ice-sheet models is used, simulating glaciations on Eurasia, North America, Greenland and Antarctica. The ice-sheet models use a combination of the shallow ice and shelf approximations to determine sheet, shelf and sliding velocities. The framework consists of an inverse forward modelling approach to derive a self-consistent record of temperature and ice volume from deep-sea benthic δ¹⁸O data over the past 1 million years, a proxy for ice volume and temperature. It is shown that for both eustatic sea level and sea water δ¹⁸O changes, the Eurasian and North American ice sheets are responsible for the largest part of the variability. The combined contribution of the Antarctic and Greenland ice sheets is about 10 % for sea level and about 20 % for sea water δ¹⁸O during glacial maxima. However, changes in interglacials are mainly caused by melt of the Greenland and Antarctic ice sheets, with an average time lag of 4 kyr between melt and temperature. Furthermore, we have tested the separate response to changes in temperature and sea level for each ice sheet, indicating that ice volume can be significantly influenced by changes in eustatic sea level alone. Hence, showing the importance of a simultaneous simulation of all four ice sheets. This paper describes the first complete simulation of global ice-volume variations over the late Pleistocene with the possibility to model changes above and below present-day ice volume, constrained by observations of benthic δ¹⁸O proxy data.     

中国地区夏季6～8月云水含量的垂直分布特征

基于观测资料的夏季云水含量时空分布情况对于数值天气预报、气候预测以及人工影响天气试验都十分重要。本文利用CloudSat卫星资料, 分析了2006～2008年中国地区夏季月平均云水含量的垂直和区域变化特征。结果显示, 青藏高原地形以及东亚夏季风对月平均云含水量分布具有明显影响。中国中部纬度上对流层中层的月平均液态水含量比南部及北部的量值大。各月平均云液水含量垂直廓线存在两个不同高度上的峰值区, 原因可能主要是受大尺度参数的控制, 以及受到青藏高原和东亚季风环流的影响。平均冰水含量纬向垂直分布的高值区主要在对流层中上部。本文中所揭示的云水含量特征为天气和气候模式改进、人工影响天气及云—辐射相互作用提供了重要的基础信息。     

Snowpack model estimates of the mass balance of the Greenland ice sheet and its changes over the twentyfirst century

The performance of a snow cover model in capturing the ablation on the Greenland ice sheet is evaluated. This model allows an explicit calculation of the formation of melt water, of the fraction of melt water which re-freezes, and of runoff in the ablation region. The input climate variables to the snowpack model come from two climate models. While the higher resolution general circulation model (ECHAM 4), is closest to observations in its estimate of accumulation, it fails to give accurate results in its predictions of runoff, primarily in the southern half of the ice sheet. The two-dimensional low-resolution climate model (MIT 2D LO) produces estimates of runoff from the Greenland ice sheet within the range of uncertainty of the Inter governmental Panel on Climate Change (IPCC1) 1995 estimates. Both models reproduce some of the characteristics of the extent of the wet snow zone observed with satellite remote sensing; the MIT model is closer to observations in terms of areal extent and intensity of the melting in the southern half of the ice-sheet in July and August while the ECHAM model reproduces melting in the northern half of the ice sheet well. Changes in runoff from Greenland and Antarctica are often cited as one of the major concerns linked to anthropogenic changes in climate. Because it is based on physical principles and relies on the surface energy balance as input, the snow cover model can respond to the current climatic forcing as well as to future changes in climate on the century time scale without the limitations inherent in empirical parametrizations. For a reference climate scenario similar to the IPCC's IS92a, the model projects that the Greenland ice sheet does not contribute significantly to changes in the level of the ocean over the twenty-first century. Increases in accumulation over the central portion of the ice sheet offset most of the increase in melting and runoff, which takes place along the margins of the ice sheet. The range of uncertainty in the predictions of sea-level rise is estimated by repeating the calculation with the MIT model for seven climate change scenarios. The range is –0.5 to 1.7 cm.     

冰云短波辐射特性参数化

该文对卷云和高层云分别构造了15种冰晶尺度分布，3种云高和云厚，总共30种冰水含量、90种冰水程长模式。对尺度小于与大于30 μm的冰晶粒子的单次散射特性分别用表面积等效冰球Mie理论和射线光学理论进行计算。采用delta-Eddington法计算了在4种不同地面反射率下冰云的多重散射辐射传输特性，提出了冰云辐射特性的参数化公式，指出参数化公式可以用于气候模式或大气环流模式。     

青藏高原那曲地区一次深对流云垂直结构的多源卫星和地基雷达观测对比分析

为了加强对青藏高原深对流云垂直结构的深入认识,利用TRMM、CloudSat和Aqua多源卫星观测资料及地基垂直指向雷达（C波段调频连续波雷达和KA波段毫米波云雷达）资料,对第三次青藏高原大气科学试验期间2014年7月9日13-16时（北京时）发生在那曲气象站附近的深厚强对流云和那曲气象站以西100 km左右的深厚弱对流云的垂直结构特征进行了分析,得到的结果如下:（1）深厚强对流云和深厚弱对流云的水平尺度均较小（10-20 km）,垂直发展高度较高（15-16 km,均指海拔高度）;深厚强对流云在0℃层以下雷达反射率因子递增非常快,表明对流云内固态降水粒子下落至0℃层以下后融化过程有很重要的作用;在对流减弱阶段有明显的0℃层亮带出现,亮带位于5.5 km左右（距地1 km）;（2）对比TRMM测雨雷达和C波段调频连续波雷达观测到的雷达反射率因子,发现TRMM测雨雷达在11 km以下存在高估;（3）深对流云主要为冰相云,云内10 km以上主要是丰富小冰粒子,而10 km以下是较少的大冰晶粒子;深厚强对流云和深厚弱对流云的微物理过程都主要包括混合相过程和冰化过程,混合相过程分为两种:一种是-25℃（深厚强对流云）或-29℃（深厚弱对流云）高度以下以凇附增长为主,另一种是该高度以上主要以冰晶聚合、凝华增长为主,该过程冰晶粒子有效半径增长较快。这些空基和地基的观测证据进一步揭示了青藏高原深对流云的垂直结构特征,为模式模拟青藏高原深对流云的检验提供了依据。      

Climate scenarios of sea level rise for the northeast Atlantic Ocean: a study including the effects of ocean dynamics and gravity changes induced by ice melt

Here we present a set of regional climate scenarios of sea level rise for the northeast Atlantic Ocean. In this study, the latest observations and results obtained with state-of-the-art climate models are combined. In addition, regional effects due to ocean dynamics and changes in the Earth’s gravity field induced by melting of land-based ice masses have been taken into account. The climate scenarios are constructed for the target years 2050 and 2100, for both a moderate and a large rise in global mean atmospheric temperature (2 °C and 4 °C in 2100 respectively). The climate scenarios contain contributions from changes in ocean density (global thermal expansion and local steric changes related to changing ocean dynamics) and changes in ocean mass (melting of mountain glaciers and ice caps, changes in the Greenland and Antarctic ice sheets, and (minor) terrestrial water-storage contributions). All major components depend on the global temperature rise achieved in the target periods considered. The resulting set of climate scenarios represents our best estimate of twenty-first century sea level rise in the northeast Atlantic Ocean, given the current understanding of the various contributions. For 2100, they yield a local rise of 30 to 55 cm and 40 to 80 cm for the moderate and large rise in global mean atmospheric temperature, respectively.     

On modelling the seasonal thermodynamic cycle of sea ice in studies of climatic change

Recent work in modelling climatic changes due to increased atmospheric CO₂ has shown the maximum change to occur in the polar regions as a result of seasonal reductions in sea ice coverage. Typically, sea ice thermodynamics is modelled in a very simple way, whereby the storage of both sensible and latent heat within the ice is ignored, and the effects of snow cover on conductivity and on surface albedo and of oceanic heat flux on bottom ablation may also be neglected. This paper considers whether omission of these processes is justified within the context of quantitatively determining regional climatic changes. A related question, whether omission of ice dynamics can be justified, is not considered.Relatively complete one-dimensional models of sea-ice thermodynamics have previously been developed and tested for a variety of environmental conditions by Maykut and Untersteiner (1969, 1971) and by Semtner (1976). A simpler model which neglects the storage of sensible and latent heat is described in the Appendix to Semtner (1976). In that model, the errors in annual-mean ice thickness which would arise from neglect of heat storage can be compensated by increases in albedo and in conductivity. Here we examine the seasonal cycle of ice thickness predicted by such a model and find significant errors in phase (one month lead) and in amplitude (50% overestimate). The amplitude errors are enhanced as snowfall and oceanic heat flux diminish (or are neglected). This suggests that substantial errors may occur in climate simulations which use very simple formulations of sea ice thermodynamics, whereby early and excessive melting exaggerates the seasonal disappearance of sea ice.To illustrate the above point, two models are configured to examine the local response of Arctic sea ice to a quadrupling of atmospheric CO₂. The first model neglects a number of physical processes and mimics the behavior of sea ice found in Manabe and Stouffer (1980), both for present and enhanced levels of CO₂. The more complete second model gives a better simulation of Arctic ice for the present level of CO₂ and shows a reduced response to CO₂ quadrupling relative to that in Manabe and Stouffer (1980). In particular, the change in surface temperature is cut by a factor of two. In view of this result, a more complete treatment of sea ice thermodynamics would seem warranted in further studies of climate change. Only a minor computational increase is required.A portion of this study is supported by the U.S. Department of Energy as a part of its Carbon Dioxide Research Program.The National Center for Atmospheric Research is sponsored by the National Science Foundation.