Heat,momentum and moisture budgets of the katabatic layer over the melting zone of the west Greenland ice sheet in summer

The budgets of momentum, heat and moisture of the atmospheric boundary layer overlying the melting zone of the west Greenland ice sheet during an 8-day period in summer are calculated. To do so, the governing budget equations are derived and presented in terms of vertically averaged quantities. Moreover, stationarity is assumed in the present study. Measurements collected during the GIMEX-91 experiment are used to calculate the contribution of the different terms in the equations to the budget.During summer, a well developed katabatic wind system is present over the melting zone of the Greenland ice sheet. The budgets show that advection in the katabatic layer is small for momentum, heat and humidity, when the horizontal length scale of the integration area is sufficiently large (>50 km). This indicates that in principle one-dimensional atmospheric models can be used to study the boundary layer over the melting zone of the Greenland ice sheet. The background stratification plays a crucial role in the heat and moisture budget. Vertical divergence of longwave radiation provides one-third and the turbulent flux of sensible heat the rest of the cooling of the boundary layer. Moisture is added to the boundary layer by evaporation which is a significant term in the moisture budget. Negative buoyancy (katabatic forcing) dominates the momentum budget in the downslope direction. Coriolis forcing is important, stressing the large spatial scale of the katabatic winds on the Greenland ice sheet.     

The small ice cap instability in seasonal energy balance models

Results from a two-dimensional energy balance model with a realistic land-ocean distribution show that the small ice cap instability exists in the Southern Hemisphere, but not in the Northern Hemisphere. A series of experiments with a one-dimensional energy balance model with idealized geography are used to study the roles of the seasonal cycle and the land-ocean distribution. The results indicate that the seasonal cycle and land-ocean distribution can influence the strength of the albedo feedback, which is responsible for the small ice cap instability, through two factors: the temperature gradient and the amplitude of the seasonal cycle. The land-ocean distribution in the Southern Hemisphere favors the small ice cap instability, while the land-ocean distribution in the Northern Hemisphere does not. Because of the longitudinal variations of land-ocean distribution in the Northern Hemisphere, the behavior of ice lines in the Northern Hemisphere cannot be simulated and explained by the model with zonally symmetric land-ocean distribution. Model results suggest that the small ice cap instability may be a possible mechanism for the formation of the Antarctic icesheet. The model results cast doubt, however, on the role of the small ice cap instability in Northern Hemisphere glaciations. Offprint requests to: J Huang     

Sensitivity of a Greenland ice sheet model to ice flow and ablation parameters: consequences for the evolution through the last climatic cycle

Sensitivity experiments are conducted to test the influence of poorly known model parameters on the simulation of the Greenland ice sheet by means of a three dimensional numerical model including the mechanical and thermal processes within the ice. Two types of experiments are performed: steady-state climatic conditions and simulations over the last climatic cycle with a climatic forcing derived from the GRIP record. The experiments show that the maximum altitude of the ice sheet depends on the ice flow parameters (deformation and sliding law coefficients, geothermal flux) and that it is low when the ice flow is fast. On the other hand, the maximum altitude is not sensitive to the ablation strength and consequently during the climatic cycle it is driven by changes in accumulation rate. The ice sheet extension shows the opposite sensitivity: it is barely affected by ice flow velocity and the ice covered area is smaller for large ablation coefficients. For colder climates, when there is no ablation, the ice sheet extension depends on the sea level. An interesting result is that the variations with time of the altitude at the ice divide (Summit) do not depend on the parameters we tested. The present modelled ice sheets resulting from the climatic cycle experiments are compared with the present measured ice sheet in order to find the set of parameters that gives the best fit between modelled and measured geometry. It seems that, compared to the parameter set most commonly used, higher ablation rate coefficents must be used. Received: 19 September 1995 / Accepted: 30 May 1996     

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 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 effect of land surface changes on Eemian climate

Transient experiments for the Eemian (128–113 ky BP) were performed with a complex, coupled earth system model, including atmosphere, ocean, terrestrial biosphere and marine biogeochemistry. In order to investigate the effect of land surface parameters (background albedo, vegetation and tree fraction and roughness length) on the simulated changes during the Eemian, simulations with interactive coupling between climate and vegetation were compared with additional experiments in which these feedbacks were suppressed. The experiments show that the influence of land surface on climate is mainly caused by changes in the albedo. For the northern hemisphere high latitudes, land surface albedo is changed partially due to the direct albedo effect of the conversion of grasses into forest, but the indirect effect of forests on snow albedo appears to be the major factor influencing the total absorption of solar radiation. The Western Sahara region experiences large changes in land surface albedo due to the appearance of vegetation between 128 and 120 ky BP. These local land surface albedo changes can be as much as 20%, thereby affecting the local as well as the global energy balance. On a global scale, latent heat loss over land increases more than 10% for 126 ky BP compared to present-day.     

An investigation of the small ice cap instability in the Southern Hemisphere with a coupled atmosphere-sea ice-ocean-terrestrial ice model

 A simple climate model has been developed to investigate the existence of the small ice cap instability in the Southern Hemisphere. The model consists of four coupled components: an atmospheric energy balance model, a thermodynamic snow-sea ice model, an oceanic mixed layer model and a terrestrial ice model. Results from a series of experiments involving different degrees of coupling in the model show that the instability appears only in those cases when an explicit representation of the Antarctic ice sheet is not included in the model. In order to determine which physical processes in the ice sheet model lead to a stabilization of the system we have conducted several sensitivity experiments in each of which a given ice sheet process has been removed from the control formulation of the model. Results from these experiments suggest that the feedback between the elevation of the ice sheet and the snow accumulation-ice ablation balance is responsible for the disappearance of the small ice cap instability in our simulation. In the model, the mass balance of the ice sheet depends on the air temperature at sea level corrected for altitude and it is, therefore, a function of surface elevation. This altitude-mass balance feedback effectively decouples the location of the ice edge from any specific sea level isotherm, thus decreasing the model sensitivity to the albedo-temperature feedback, which is responsible for the appearance of the instability. It is also shown that the elevation-radiative cooling feedback tends to stabilize the ice sheet, although its effect does not seem to be strong enough to remove the instability. Another interesting result is that for those simulations which include the terrestrial ice model with elevation-dependent surface mass balance, hysteresis is exhibited, where for a given level of external forcing, two stable solutions with different, non-zero ice-sheet volume and area and different air and ocean temperature fields occur. However, no unstable transition between the two solutions is ever observed. Our results suggest that the small ice cap instability mechanism could be unsuitable for explaining the inception of glaciation in Antarctica. Received: 14 April 1997 / Accepted: 22 October 1997     

Investigating the variability of the Arctic sea ice thickness in response to a stochastic thermodynamic atmospheric forcing

The natural low frequency variability of the sea-ice thickness in the Arctic is investigated based on a 10 000 years simulation with a one-dimensional thermodynamic sea-ice model forced by random perturbations of the air surface temperature and solar radiation. The simulation results suggest that atmospheric random perturbations are integrated by the sea-ice. Moreover those perturbations occurring at the onset of ice melting force the largest ice thickness anomalies, which are successively amplified in summer by the albedo feedback and damped in winter by the feedback of the heat conduction through the ice. They also result in a global shift of the melting season which, in the mean annual cycle, leads to earlier melting as compared to the mean climatological cycle. The power spectrum of the ice anomalies suggests that the thickness of the perennial ice should vary preferentially on a time scale of approximately 20 years. The shape of the spectrum is consistent with that of a first order Markov process in which the characteristic time scale of the ice fluctuations would be the relaxation time scale associated with the linear feedback. The equivalent Markov model is constructed by linearizing the ice growth rate anomaly equations and allows us to derive an analytical expression of the feedback and of the forcing of the anomalies. The characteristic time scale depends explicitly on those model parameters involved in the atmosphere-ice interaction but also on the mean seasonal characteristics of the forcing and of the ice thickness. Received: 18 August 1999 / Accepted: 10 May 2000     

Recent precipitation trends over the polar ice sheets

Summary Meteorological and glaciological analyses are integrated to examine the precipitation trends during the last three decades over the ice sheets covering Antarctica and Greenland. For Antarctica, the best data source is provided by glaciologically-measured trends of snow accumulation, and for limited sectors of East Antarctica consistency with precipitation amounts calculated from the atmospheric water balance equation is obtained. For Greenland, precipitation rates parameterized from atmospheric analyses yield the only comprehensive depiction. The precipitation rate over Antarctica appears to have increased by about 5% over a time period spanning the accumulation means for the 1955–65 to 1965–75 periods, while over Greenland it has decreased by about 15% since 1983 with a secondary increase over the southern part of the ice sheet starting in 1977. At the end of the 10-year overlapping period, the global sea-level impact of the precipitation changes over Antarctica dominates that for Greenland and yields a net ice-sheet precipitation contribution of roughly 0.02 mm yr^–1. These changes are likely due to marked variations in the cyclonic forcing affecting the ice sheets, but are only weakly reflected in the temperature regime, consistent with the episodic nature of cyclonic precipitation. These conclusions are not founded on high quality data bases. The importance of such changes for understanding global sea-level variations argues for a modest research effort to collect simultaneous meteorological and glaciological observations in order to describe and understand the current precipitation variations over both ice sheets. Some suggestions are offered for steps that could be taken.With 8 Figures     

Energy balance of the intertidal zone of Western Hudson bay II: Ice‐dominated periods and seasonal patterns

Abstract

This study treats the energy balance during fast‐ice and floating‐ice conditions and examines overall seasonal patterns. The rate of ablation of the fast ice was controlled equally by net radiation and air temperature. The ratio of net/solar radiation increased 2.5 times during the ablation period owing to the decrease in ice albedo. Air temperature in the ablation zone was up to 8°C colder than that over the adjacent snow‐free terrestrial surface and remained near 0°Cfor the full ablation period. The sensible heat flux was small and downward (negative), whereas the evaporative heat flux was small and positive. Thus, the energy used in melting the ice was approximately equal to that provided by the net radiation. Above‐freezing air temperatures decreased the albedo through surface melting thus increasing net radiation. This combination of higher temperature and large net radiation was associated with offshore winds and resulted in large ablation relative to periods with colder onshore winds.

The floating‐ice period is one of great variability owing to changing ice conditions, variable current behaviour, tidal cycles and changing wind direction. The intertidal zone acts as a major heat sink, both early and late in the floating‐ice period. The turbulent heat fluxes were small and were either positive or negative. Nearly all of the energy from net radiation was used in melting ice and in warming tidal water during high tide and in warming the residual tidal ponds and in melting stranded ice rafts during low tide.

The overall study period, from May to September, included most of the season of positive radiation balance and above‐freezing temperatures. Winds were dominantly onshore in the first half of the period and equally onshore and offshore in the second half. Wind frequencies resembled longer term averages for other stations on James Bay and Hudson Bay. The ratio of net to solar radiation was at a maximum during the ice‐free period in August, whereas for adjacent terrestrial surfaces, it was largest at the summer solstice. Land‐sea breezes first developed in mid‐July and were influential in making offshore winds the dominant nocturnal regime. As a result, offshore winds were associated with small magnitudes of net radiation. Onshore winds were more than 5°C colder than those blowing offshore and their vapour pressure deficits were three times smaller. Convective heat fluxes were small for onshore winds and very small and usually negative for offshore winds. For all wind directions throughout the period, most of the available radiant energy was used to melt ice and to heat the sea water. This is a pattern similar to that of the ice‐covered or open sea and dissimilar to that of the adjacent terrestrial environment. It implies that the main energy‐balance transitions, during onshore airflow, occur at the high‐tide line.     

Effects of a melted greenland ice sheet on climate, vegetation, and the cryosphere

This paper investigates the possible implications for the earth-system of a melting of the Greenland ice-sheet. Such a melting is a possible result of increased high latitude temperatures due to increasing anthropogenic greenhouse gas emissions. Using an atmosphere-ocean general circulation model (AOGCM), we investigate the effects of the removal of the ice sheet on atmospheric temperatures, circulation, and precipitation. We find that locally over Greenland, there is a warming associated directly with the altitude change in winter, and the altitude and albedo change in summer. Outside of Greenland, the largest signal is a cooling over the Barents sea in winter. We attribute this cooling to a decrease in poleward heat transport in the region due to changes to the time mean circulation and eddies, and interaction with sea-ice. The simulated climate is used to force a vegetation model and an ice-sheet model. We find that the Greenland climate in the absence of an ice sheet supports the growth of trees in southern Greenland, and grass in central Greenland. We find that the ice sheet is likely to regrow following a melting of the Greenland ice sheet, the subsequent rebound of its bedrock, and a return to present day atmospheric CO₂ concentrations. This regrowth is due to the high altitude bedrock in eastern Greenland which allows the growth of glaciers which develop into an ice sheet.     

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.     

Microphysical and radiative effects of ice clouds on diurnal variations of tropical convective and stratiform rainfall

Microphysical and radiative effects of ice clouds on diurnal variations of tropical convective and stratiform rainfall are examined with the equilibrium simulation data from three experiments conducted with a two-dimensional cloud resolving model with imposed temporally and zonally invariant winds and sea surface temperature and zero mean vertical velocity. The experiment without ice radiative effects is compared with the control experiment with ice microphysics (both the ice radiative and microphysical effects) to study effects of ice radiative effects on diurnal rainfall variations whereas it is compared with the experiment without ice microphysics to examine ice microphysical effects on the diurnal rainfall variations. The ice radiative processes mainly affect diurnal cycle of convective rainfall whereas the ice microphysical processes have important impacts on the diurnal cycles of both convective and stratiform rainfall. Turning off the ice radiative effects generally enhances convective rainfall during the morning and evening and suppresses convective rainfall in the afternoon whereas turning off the ice microphysical effects generally suppresses convective and stratiform rainfall during the morning and enhances convective and stratiform rainfall in the afternoon and evening. The ice radiative and microphysical effects on the diurnal cycle of surface rainfall are mainly associated with that of vapor condensation and deposition, which is controlled by air temperature through saturation specific humidity. The ice effects on the diurnal cycle of local temperature tendency are largely explained by that of latent heating since the diurnal cycle of radiation is insensitive to the ice effects.     

The effects of the Arctic haze as determined from airborne radiometric measurements during AGASP II

The interaction of the Aretic winter aerosol (Arctic haze) with solar radiation produces changes in the radiation field that result in the enhancement of scattering and absorption processes which alter the energy balance and solar energy distribution in the Arctic atmosphere-surface system. During the second Arctic Gas and Aerosols Sampling Project (AGASP II) field experiment, we measured radiation parameters using the NOAA WP-3D research aircraft as a platform. State-of-the-art instrumentation was used to measure in situ the absorption of solar radiation by the Arctic atmosphere during severe haze events. Simultaneously with the absorption measurements, we determined optical depths, and total, direct, and scattered radiation fields. All optical measurements were made at spectral bands centered at 412, 500, 675, and 778 nm and with a bandpass of 10 nm. With this selection of spectral regions we concentrated on the measurement of the radiative effects of the aerosol excluding most of the contributions by the gaseous components of the atmosphere. An additional measurement performed during these experiments was the determination of total solar spectrum fluxes. The experimentally determined parameters were used to define an aerosol model that was employed to deduce the absorption by the aerosols over the full solar spectrum and to calculate atmospheric heating rate profiles. The analyses summarized above allowed us to deduce the magnitude of the change in some important parameters. For example, we found changes in instantaneous heating rate of up to about 0.6 K/day. Besides the increased absorption (30 to 40%) and scattering of radiation by the atmosphere, the haze reduces the surface absorption of solar energy by 6 to 10% and the effective planetary albedo over ice surfaces by 3 to 6%. The vertical distribution of the absorbing aerosol is inferred from the flux measurements. Values for the specific absorption of carbon are found to be around 6 m²/g for externally mixed aerosol and about 11.7 m²/g for internally mixed aerosol. A complete study of the radiative effects of the Arctic haze should include infrared measurements and calculations as well as physics of the ice, snow, and water surfaces.     

Response of the Antarctic ice sheet to future greenhouse warming

Possible future changes in land ice volume are mentioned frequently as an important aspect of the greenhouse problem. This paper deals with the response of the Antarctic ice sheet and presents a tentative projection of changes in global sea level for the next few hundred years, due to changes in its surface mass balance. We imposed a temperature scenario, in which surface air temperature rises to 4.2° C in the year 2100 AD and is kept constant afterwards. As GCM studies seem to indicate a higher temperature increase in polar latitudes, the response to a more extreme scenario (warming doubled) has also been investigated. The mass balance model, driven by these temperature perturbations, consists of two parts: the accumulation rate is derived from present observed values and is consequently perturbed in proportion to the saturated vapour pressure at the temperature above the inversion layer. The ablation model is based on the degree-day method. It accounts for the daily temperature cycle, uses a different degree-day factor for snow and ice melting and treats refreezing of melt water in a simple way. According to this mass balance model, the amount of accumulation over the entire ice sheet is presently 24.06 × 10¹¹ m³ of ice, and no runoff takes place. A 1°C uniform warming is then calculated to increase the overall mass balance by an amount of 1.43 × 10¹¹ m³ of ice, corresponding to a lowering of global sea level with 0.36 mm/yr. A temperature increase of 5.3°C is needed for the increase in ablation to become more important than the increase in accumulation and the temperature would have to rise by as much as 11.4°C to produce a zero surface mass balance. Imposing the Bellagio-scenario and accumulating changes in mass balance forward in time (static response) would then lower global sea level by 9 cm by 2100 AD. In a subsequent run with a high-resolution 3-D thermomechanic model of the ice sheet, it turns out that the dynamic response of the ice sheet (as compared to the direct effect of the changes in surface mass balance) becomes significant after 100 years or so. Ice-discharge across the grounding-line increases, and eventually leads to grounding-line retreat. This is particularly evident in the extreme case scenario and is important along the Antarctic Peninsula and the overdeepened outlet glaciers along the East Antarctic coast. Grounding-line retreat in the Ross and Ronne-Filchner ice shelves, on the other hand, is small or absent.     

二维能量平衡模式中极冰对气候的影响——解的稳定性分析(一)

本文从理论上分析了二维能量平衡气候模式解的稳定性。结果表明,解的上分支,即冰界随太阳常数增加而北移的解是稳定的;而解的下分支,即冰界随太阳常数增加而南移的解是不稳定的。由于现在地球上的气候是处在解的上分支,所以可以认为,现在气候状态是稳定的。     

Applying the energy- and water balance model for incorporation of the cryospheric component into a climate model. Part II. Modeled mass balance on the green land ice sheet surface

The Greenland ice sheet is a very important potential source of fresh water inflow to the World Ocean under warming climate conditions. Apparently, it was the same during the Last Interglacial 130-115 thousand years ago. In order to quantify input of the Greenland ice sheet to the rise of the global mean sea level in the past or in the future, we include a surface mass balance model block into the Earth System Model. The computational algorithm is based on the calculation of energy balance on the ice sheet surface. The key tuning parameter of the model is the daily amplitude of air surface temperature. It defines the area and the rate of snow or ice melting. The range of possible values of this parameter is determined during a series of numerical experiments. High sensitivity of meltwater runoff volume to surface air temperature amplitude is revealed.     

The Pangean ice age: studies with a coupled climate-ice sheet model

 Application of an ice sheet model developed for the Pleistocene to the extensive Carboniferous glaciation on Gondwana yields an ice sheet which has several features consistent with observations. While complete deglaciation is not achieved without CO₂ changes, the Milankovich-induced fluctuations in ice sheet volume are comparable to Pleistocene glacial/ interglacial signals. This result is shown to hold for a large fraction of physically reasonable parameter space. The model also exhibits multiple equilibria and sharp bifurcations, as infinitesimal changes in the solar constant or precipitation can lead to a qualitatively different climate. The success of the model in predicting ice location in an environment quite different from the Pleistocene provides additional support for the robustness of the basic model physics and suggests that the model can be applied with some confidence to other pre-Pleistocene glaciations. Received: 30 June 1998 / Accepted: 5 January 1999     

Potential effect of nuclear war smokefall on sea ice

A large nuclear war could produce massive quantities of smoke from burning cities and industries. A portion of this smoke would fall out on Arctic sea ice, thus lowering its albedo and potentially increasing the solar energy absorbed by the ice and the snow that covers it. We use a one-dimensional thermodynamic sea ice model to examine the effect of smokefall on the seasonal variation of sea ice. In particular, we test the sensitivity of the model results to the time of year, duration, and latitude of smokefall.Sea ice thickness variations and the period of summer ice-free conditions are sensitive to the season of smokefall. The largest sea ice perturbations are generated by smokefall in spring. In this case the period of ice-free conditions during the summer can increase by 2 – 3.5 months between 67.5° N and 82.5° N. In any given season, the annual cycle of sea ice is not very sensitive to the duration of smokefall. The equilibrium annual cycle of sea ice variation is restored within a few years of smokefall when the smoke is flushed out of the ice/snow system.Since the sea ice model used here is not a comprehensive global climate model, it is difficult to predict the mid-latitude climate effects of the massive, but temporary, Arctic sea ice changes. However, our results suggest that future global climate model simulations of the effects of nuclear war smoke include interactive sea ice calculations.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

A two-dimensional energy balance climate model including radiation and ice caps-albedo feedback

A simplified two-dimensional energy balance climate model including the solar and infrared radiation transports, the turbulent exchanges of heat in vertical and horizontal directions and the ice caps-albedo feedback is developed The solutions show that if the atmosphere is considered as a grey body and the grey coefficient depends upon the distributions of absorption medium and cloudiness, both horizontal and vertical distributions of temperature are identical to the observation.On the other hand, comparing the models that the atmosphere is considered as a grey body with ones that the infrared radiation is parameterized as a linear function of temperature, as was considered by Budyko, Sellers(1969), then the results show that even though both of them can obtain the earth's surface temperature in agreement with the observation, the sensitiv ity of the climate to the changes of solar constant is very different. In the former case,the requirement for the ice edge to move southward from the normal 72°N to 50°N(i.e. where the glacial climate would take place) is that the solar constant should decrease by 13% to 16%. However, in the latter case, the climate is highly sensitive to the changes of solar radiation. In this case, the requirement of solar radiation occurring in the glacial climate should decrease by 2% to 6%. According to the investigations mentioned above we must be careful when the parameterizations of the radiation and other processes are conducted in a climate model, otherwise the reliability of the results is suspicious.