A comparison of the Vostok ice deuterium record and series from Southern Ocean core MD 88-770 over the last two glacial-interglacial cycles

Taking advantage of the fact that the Vostok deuterium (δD) record now covers almost two entire climatic cycles, we have applied the orbital tuning approach to derive an age-depth relation for the Vostok ice core, which is consistent with the SPECMAP marine time scale. A second age-depth relation for Vostok was obtained by correlating the ice isotope content with estimates of sea surface temperature from Southern Ocean core MD 88-770. Both methods lead to a close correspondence between Vostok and MD 88-770 time series. However, the coherence between the correlated δD and insolation is much lower than between the orbitally tuned δD and insolation. This reflects the lower accuracy of the correlation method with respect to direct orbital tuning. We compared the ice and marine records, set in a common temporal framework, in the time and frequency domains. Our results indicate that changes in the Antarctic air temperature quite clearly lead variations in global ice volume in the obliquity and precession frequency bands. Moreover, the average phase we estimated between the filtered δD and insolation signals at precessional frequencies indicates that variations in the southern high latitude surface temperature could be induced by changes in insolation taking place during a large period of the summer in northern low latitudes or winter in southern low latitudes. The relatively large lag found between Vostok δD variations and obliquity-driven changes in insolation suggests that variations in the local radiative balance are not the only mechanism responsible for the variability in surface temperature at those frequencies. Finally, in contrast to the cross-spectral analysis method used in previous studies, the method we use here to estimate the phases can reveal errors in cross-correlations with orbitally tuned chronologies. Received: 11 April 1995 / Accepted: 19 July 1995     

The impact of cold North Atlantic sea surface temperatures on climate: implications for the Younger Dryas cooling (11–10 k)

The sensitivity of global climate to colder North Atlantic sea surface temperatures is in vestigated with the use of the GISS general circulation model. North Atlantic ocean temperatures 18,000 B.P., resembling those prevalent during the Younger Dryas, were incorporated into the model of the present climate and also into an experiment using orbital parameters and land ice characteristic of 11,000 B.P. The results show that with both 11,000 B.P. and present conditions the colder ocean temperatures produce cooling over western and central Europe, in good agreement with Younger Dryas paleoclimatic evidence. Cooling also occurs over extreme eastern North America, although the precise magnitude and location depends upon the specification of ocean temperature change in the western Atlantic. Despite the presence of increased land ice and colder ocean temperatures, the Younger Dryas summer air temperatures at Northern Hemisphere midlatitudes in the model are warmer than those of today due to changes in the orbital parameters, chiefly precession, and atmospheric subsidence at the perimeter of the ice sheets.     

Phase Modulation Effect of the Rubincam Insolation Variations on Climate Change

Summary For astronomical seasons, Rubincam insolation deviations at latitude 65° N varied from 218.50 Wm⁻² to 225.75 Wm⁻² (3%). The periodicity of the insolation cycles varied from 36.7 Kyr to 44.7 Kyr (20%) due to phase shift. Phase shift of insolation variations can induce asymmetry of the insolation cycles, permitting rapid melting and prolonged glaciation of ice sheets to occur. For instance, an abnormal decrease of the insolation frequency during the longer period of glacial interval would prolong glaciation into deep ice age. In this study, we apply Rubincam’s insolation equations to investigate the phase shift effect of insolation variations on climate change. Using complex transforms of the changing insolation, we have detected a phase modulation signal in the insolation variations. As a result, an especially new and interesting series of the phase-related insolation pulsation is established. The phase modulated insolation is then introduced as a forcing function into energy balance climate models. Results of model computations shed new insights into the spectrum of the paleoclimatic proxy-data. It is shown that phase modulation of the insolation may provide an appropriate and complete external forcing mechanism to which the climate system would respond. The 100 Kyr cycle of the frequency modulation of the Rubincam’s insolation variations does seem adequate to change the climate. Received July 16, 1997 Revised May 18, 1998     

Past and future polar amplification of climate change: climate model intercomparisons and ice-core constraints

Climate model simulations available from the PMIP1, PMIP2 and CMIP (IPCC-AR4) intercomparison projects for past and future climate change simulations are examined in terms of polar temperature changes in comparison to global temperature changes and with respect to pre-industrial reference simulations. For the mid-Holocene (MH, 6,000 years ago), the models are forced by changes in the Earth’s orbital parameters. The MH PMIP1 atmosphere-only simulations conducted with sea surface temperatures fixed to modern conditions show no MH consistent response for the poles, whereas the new PMIP2 coupled atmosphere–ocean climate models systematically simulate a significant MH warming both for Greenland (but smaller than ice-core based estimates) and Antarctica (consistent with the range of ice-core based range). In both PMIP1 and PMIP2, the MH annual mean changes in global temperature are negligible, consistent with the MH orbital forcing. The simulated last glacial maximum (LGM, 21,000 years ago) to pre-industrial change in global mean temperature ranges between 3 and 7°C in PMIP1 and PMIP2 model runs, similar to the range of temperature change expected from a quadrupling of atmospheric CO₂ concentrations in the CMIP simulations. Both LGM and future climate simulations are associated with a polar amplification of climate change. The range of glacial polar amplification in Greenland is strongly dependent on the ice sheet elevation changes prescribed to the climate models. All PMIP2 simulations systematically underestimate the reconstructed glacial–interglacial Greenland temperature change, while some of the simulations do capture the reconstructed glacial–interglacial Antarctic temperature change. Uncertainties in the prescribed central ice cap elevation cannot account for the temperature change underestimation by climate models. The variety of climate model sensitivities enables the exploration of the relative changes in polar temperature with respect to changes in global temperatures. Simulated changes of polar temperatures are strongly related to changes in simulated global temperatures for both future and LGM climates, confirming that ice-core-based reconstructions provide quantitative insights on global climate changes. An erratum to this article can be found at     

Present And Past Ice-Sheet Mass Balance Simulations For Greenland And The Tibetan Plateau

Net annual mass balance was evaluated for Greenland and the Tibetan Plateau using the meteorological forcings from the NCEP reanalysis and two GCMs (FOAM1.0 and CSM1.4) for modern climate and for different time periods extending back to the beginning of the Holocene (11,000 years ago) for the climate models. The ice-sheet budget calculations, using the degree day methodology, were performed on a finer grid than the model output by interpolating monthly precipitation and surface temperature and correcting the latter to account for the GCMs smoothed topography. The computed net mass balance for Greenland in the present day is positive and it ranges between 290–300 mm water equivalent (w.e.)/year for the two models, values close to the NCEP estimate of 250 mm/year. The past climate simulations show that the Greenland mass balance has become slightly more positive since the beginning of the Holocene. The Tibetan Plateaus present-day area average net mass balance is negative and ranges between –1200 and –2000 mm w.e. /year for the two models, values bracketing the NCEP estimate of 1700 mm/year, although the balance is positive over small regions of the plateau consistent with the existence of small ice caps and glaciers. The calculated past mass balance shows an increasingly less negative value for FOAM from 11,000 years ago towards the present and expansion of the positive mass balance areas, mainly due to decreased snow ablation as the summertime insolation decreases with the changes in orbital forcing; in CSM the opposite trend occurs but changes are smaller and less systematic. The result from FOAM shows that the likelihood of ice sheets developing on the Tibetan Plateau may have increased since 11000 years ago, which is consistent with some glacial records.     

Orbital forcing of Arctic climate: mechanisms of climate response and implications for continental glaciation

Progress in understanding how terrestrial ice volume is linked to Earths orbital configuration has been impeded by the cost of simulating climate system processes relevant to glaciation over orbital time scales (10³–10⁵ years). A compromise is usually made to represent the climate system by models that are averaged over one or more spatial dimensions or by three-dimensional models that are limited to simulating particular snapshots in time. We take advantage of the short equilibration time (10 years) of a climate model consisting of a three-dimensional atmosphere coupled to a simple slab ocean to derive the equilibrium climate response to accelerated variations in Earths orbital configuration over the past 165,000 years. Prominent decreases in ice melt and increases in snowfall are simulated during three time intervals near 26, 73, and 117 thousand years ago (ka) when aphelion was in late spring and obliquity was low. There were also significant decreases in ice melt and increases in snowfall near 97 and 142 ka when eccentricity was relatively large, aphelion was in late spring, and obliquity was high or near its long term mean. These glaciation-friendly time intervals correspond to prominent and secondary phases of terrestrial ice growth seen within the marine ¹⁸O record. Both dynamical and thermal effects contribute to the increases in snowfall during these periods, through increases in storm activity and the fraction of precipitation falling as snow. The majority of the mid- to high latitude response to orbital forcing is organized by the properties of sea ice, through its influence on radiative feedbacks that nearly double the size of the orbital forcing as well as its influence on the seasonal evolution of the latitudinal temperature gradient.     

Transient simulation of the last glacial inception. Part I: glacial inception as a bifurcation in the climate system

We study the mechanisms of glacial inception by using the Earth system model of intermediate complexity, CLIMBER-2, which encompasses dynamic modules of the atmosphere, ocean, biosphere and ice sheets. Ice-sheet dynamics are described by the three-dimensional polythermal ice-sheet model SICOPOLIS. We have performed transient experiments starting at the Eemiam interglacial, at 126 ky BP (126,000 years before present). The model runs for 26 kyr with time-dependent orbital and CO₂ forcings. The model simulates a rapid expansion of the area covered by inland ice in the Northern Hemisphere, predominantly over Northern America, starting at about 117 kyr BP. During the next 7 kyr, the ice volume grows gradually in the model at a rate which corresponds to a change in sea level of 10 m per millennium. We have shown that the simulated glacial inception represents a bifurcation transition in the climate system from an interglacial to a glacial state caused by the strong snow-albedo feedback. This transition occurs when summer insolation at high latitudes of the Northern Hemisphere drops below a threshold value, which is only slightly lower than modern summer insolation. By performing long-term equilibrium runs, we find that for the present-day orbital parameters at least two different equilibrium states of the climate system exist—the glacial and the interglacial; however, for the low summer insolation corresponding to 115 kyr BP, we find only one, glacial, equilibrium state, while for the high summer insolation corresponding to 126 kyr BP only an interglacial state exists in the model.

---

Evaluation of a high-resolution regional climate simulation over Greenland

A simulation of the 1991 summer has been performed over south Greenland with a coupled atmosphere–snow regional climate model (RCM) forced by the ECMWF re-analysis. The simulation is evaluated with in-situ coastal and ice-sheet atmospheric and glaciological observations. Modelled air temperature, specific humidity, wind speed and radiative fluxes are in good agreement with the available observations, although uncertainties in the radiative transfer scheme need further investigation to improve the model’s performance. In the sub-surface snow-ice model, surface albedo is calculated from the simulated snow grain shape and size, snow depth, meltwater accumulation, cloudiness and ice albedo. The use of snow metamorphism processes allows a realistic modelling of the temporal variations in the surface albedo during both melting periods and accumulation events. Concerning the surface albedo, the main finding is that an accurate albedo simulation during the melting season strongly depends on a proper initialization of the surface conditions which mainly result from winter accumulation processes. Furthermore, in a sensitivity experiment with a constant 0.8 albedo over the whole ice sheet, the average amount of melt decreased by more than 60%, which highlights the importance of a correctly simulated surface albedo. The use of this coupled atmosphere–snow RCM offers new perspectives in the study of the Greenland surface mass balance due to the represented feedback between the surface climate and the surface albedo, which is the most sensitive parameter in energy-balance-based ablation calculations.     

Seasonally asymmetric transition of the Asian monsoon in response to ice age boundary conditions

Modulation of a monsoon under glacial forcing is examined using an atmosphere?Cocean coupled general circulation model (AOGCM) following the specifications established by Paleoclimate Modelling Intercomparison Project phase 2 (PMIP2) to understand the air?Csea?Cland interaction under different climate forcing. Several sensitivity experiments are performed in response to individual changes in the continental ice sheet, orbital parameters, and sea surface temperature (SST) in the Last Glacial Maximum (LGM: 21?ka) to evaluate the driving mechanisms for the anomalous seasonal evolution of the monsoon. Comparison of the model results in the LGM with the pre-industrial (PI) simulation shows that the Arabian Sea and Bay of Bengal are characterized by enhancement of pre-monsoon convection despite a drop in the SST encompassing the globe, while the rainfall is considerably suppressed in the subsequent monsoon period. In the LGM winter relative to the PI, anomalies in the meridional temperature gradient (MTG) between the Asian continents minus the tropical oceans become positive and are consistent with the intensified pre-monsoon circulation. The enhanced MTG anomalies can be explained by a decrease in the condensation heating relevant to the suppressed tropical convection as well as positive insolation anomalies in the higher latitude, showing an opposing view to a warmer future climate. It is also evident that a latitudinal gradient in the SST across the equator plays an important role in the enhancement of pre-monsoon rainfall. As for the summer, the sensitivity experiments imply that two ice sheets over the northern hemisphere cools the air temperature over the Asian continent, which is consistent with the reduction of MTG involved in the attenuated monsoon. The surplus pre-monsoon convection causes a decrease in the SST through increased heat loss from the ocean surface; in other words, negative ocean feedback is also responsible for the subsequent weakening of summer convection.     

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     

Climate Change in Northern Africa: The Past is Not the Future

By using a climate system model of intermediate complexity, we have simulated long-term natural climate changes occurring over the last 9000 years. The paleo-simulations in which the model is driven by orbital forcing only, i.e., by changes in insolation caused by changes in the Earth's orbit, are compared with sensitivity simulations in which various scenarios of increasing atmospheric CO₂ concentration are prescribed. Focussing on climate and vegetation change in northern Africa, we recapture the strong greening of the Sahara in the early and mid-Holocene (some 9000–6000 years ago), and we show that some expansion of grasslandinto the Sahara is theoretically possible, if the atmospheric CO₂ concentration increases well above pre-industrial values and if vegetation growth is not disturbed. Depending on the rate of CO₂ increase, vegetation migration into the Sahara can be rapid, up to 1/10th of the Saharan area per decade, but could not exceed a coverage of 45%. In ourmodel, vegetation expansion into today's Sahara is triggered by an increase in summer precipitation which is amplified by a positive feedback between vegetation and precipitation. This is valid for simulations with orbital forcing and greenhouse-gas forcing. However, we argue that the mid-Holocene climate optimum some 9000 to 6000 years ago with its marked reduction of deserts in northern Africa is not a direct analogue for future greenhouse-gas induced climate change, as previously hypothesized. Not only does the global pattern of climate change differ between the mid-Holocene model experiments and the greenhouse-gas sensitivity experiments, but the relative role of mechanisms which lead to a reduction of the Sahara also changes. Moreover, the amplitude of simulated vegetation cover changes in northern Africa is less than is estimated for mid-Holocene climate.     

The effect of snow on Antarctic sea ice simulations in a coupled atmosphere-sea ice model

 The effect of a snow cover on sea ice accretion and ablation is estimated based on the ‘zero-layer’ version sea ice model of Semtner, and is examined using a coupled atmosphere-sea ice model including feedbacks and ice dynamics effects. When snow is disregarded in the coupled model the averaged Antarctic sea ice becomes thicker. When only half of the snowfall predicted by the atmospheric model is allowed to land on the ice surface sea ice gets thicker in most of the Weddell and Ross Seas but thinner in East Antarctic in winter, with the average slightly thicker. When twice as much snowfall as predicted by the atmospheric model is assumed to land on the ice surface sea ice also gets much thicker due to the large increase of snow-ice formation. These results indicate the importance of the correct simulation of the snow cover over sea ice and snow-ice formation in the Antarctic. Our results also illustrate the complex feedback effects of the snow cover in global climate models. In this study we have also tested the use of a mean value of 0.16 Wm^-1 K^-1 instead of 0.31 for the thermal conductivity of snow in the coupled model, based on the most recent observations in the eastern Antarctic and Bellingshausen and Amundsen Seas, and have found that the sea ice distribution changes greatly, with the ice becoming much thinner by about 0.2 m in the Antarctic and about 0.4 m in the Arctic on average. This implies that the magnitude of the thermal conductivity of snow is of considerable importance for the simulation of the sea ice distribution. An appropriate value of the thermal conductivity of snow is as crucial as the depth of the snow layer and the snowfall rate in a sea ice model. The coupled climate models require accurate values of the effective thermal conductivity of snow from observations for validating the simulated sea ice distribution under the present climate conditions. Received: 20 November 1997/Accepted: 27 July 1998     

Transient simulation of the last glacial inception. Part II: sensitivity and feedback analysis

The sensitivity of the last glacial-inception (around 115 kyr BP, 115,000 years before present) to different feedback mechanisms has been analysed by using the Earth system model of intermediate complexity CLIMBER-2. CLIMBER-2 includes dynamic modules of the atmosphere, ocean, terrestrial biosphere and inland ice, the last of which was added recently by utilising the three-dimensonal polythermal ice-sheet model SICOPOLIS. We performed a set of transient experiments starting at the middle of the Eemiam interglacial and ran the model for 26,000 years with time-dependent orbital forcing and observed changes in atmospheric CO₂ concentration (CO₂ forcing). The role of vegetation and ocean feedback, CO₂ forcing, mineral dust, thermohaline circulation and orbital insolation were closely investigated. In our model, glacial inception, as a bifurcation in the climate system, appears in nearly all sensitivity runs including a run with constant atmospheric CO₂ concentration of 280 ppmv, a typical interglacial value, and simulations with prescribed present-day sea-surface temperatures or vegetation cover—although the rate of the growth of ice-sheets growth is smaller than in the case of the fully interactive model. Only if we run the fully interactive model with constant present-day insolation and apply present-day CO₂ forcing does no glacial inception appear at all. This implies that, within our model, the orbital forcing alone is sufficient to trigger the interglacial–glacial transition, while vegetation, ocean and atmospheric CO₂ concentration only provide additional, although important, positive feedbacks. In addition, we found that possible reorganisations of the thermohaline circulation influence the distribution of inland ice.     

Individual contribution of insolation and CO2 to the interglacial climates of the past 800,000?years

The individual contributions of insolation and greenhouse gases (GHG) to the interglacial climates of the past 800,000?years are quantified through simulations with a model of intermediate complexity LOVECLIM and using the factor separation technique. The interglacials are compared in terms of their forcings and responses of surface air temperature, vegetation and sea ice. The results show that the relative magnitude of the simulated interglacials is in reasonable agreement with proxy data. GHG plays a dominant role on the variations of the annual mean temperature of both the Globe and the southern high latitudes, whereas, insolation plays a dominant role on the variations of tree fraction, precipitation and of the northern high latitude temperature and sea ice. The Mid-Brunhes Event (MBE) appears to be significant only in GHG and climate variables dominated by it. The results also show that the relative importance of GHG and insolation on the warmth intensity varies from one interglacial to another. For the warmest (MIS-9 and MIS-5) and coolest (MIS-17 and MIS-13) interglacials, GHG and insolation reinforce each other. MIS-11 (MIS-15) is a warm (cool) interglacial due to its high (low) GHG concentration, its insolation contributing to a cooling (warming). MIS-7, although with high GHG concentrations, can not be classified as a warm interglacial due to it large insolation-induced cooling. Related to these two forcings, MIS-19 appears to be the best analogue for MIS-1. In the response to insolation, the annual mean temperatures averaged over the globe and over southern high latitudes are highly linearly correlated with obliquity. However, precession becomes important in the temperature of the northern high latitudes and controls the tree fraction globally. Over the polar oceans, the response during the local winters, although the available energy is small, is larger than during the local summers due to the summer remnant effect. The sensitivity to double CO₂ is the highest for the coolest interglacial.     

Climate changes and the Caspian Sea Level fluctuations in the Holocene and historical time

The influence of climate changes in Eurasia on the Caspian Sea level fluctuations in the Holocene is studied on the basis of historical data. The correspondence of the Late Khvalyn transgression to the time of the thawing of the European (Scandinavian) ice cap is revealed as well as that of the climatic optimum of the Holocene to the Boreal period (10-8 thousand years ago). In the next New-Caspian epoch, the Caspian Sea regression occurred. In the historical time, the maxima of cold snaps and Caspian Sea level are registered every five centuries.     

Investigation of summer synoptic climate controls on the mass balance of Meighen Ice Cap

Abstract

The Meighen Ice Cap synoptic climate classification system, developed from the study of six years of summer meteorological and glaciological observations, appears to account for significant variations in the energy‐ and mass‐balance climates of the ice cap. In relating the summer frequency of the three synoptic types to fourteen years of mass‐balance measurements, it was found that variations in surface conditions, solar angle and type of precipitation could be accounted for by the relative sequence of synoptic types. Further it was shown that the types could be represented by the position of the dominant 500‐mb cold Low influencing Meighen Island, thus providing a link between the mass balance and the general circulation.

Dominance of the winter pattern of a 500‐mb Low in the Hudson Bay –Baffin Island region throughout the summer season is capable of maintaining Meighen Ice Cap at its present size. A shift of the 500‐mb Low from the winter position directly to the Beaufort Sea or adjoining Polar Ocean area is capable of increasing the size of the ice cap. On the other hand, a shift of the 500‐mb vortex to the Asiatic side of the Polar Ocean before taking up position in the Beaufort Sea – Polar Ocean area produces negative mass‐balance conditions. When the 500‐mb Low remains on the Asiatic side of the Polar Ocean during most of the summer season the slow accumulation of two decades of Polar Ocean years is destroyed.     

Climate stability and the development of agricultural societies

Although Modern Man had developed long before the migration from Africa began ～ 55,000 years ago no agricultural societies developed until about ～ 10,000 years ago. In the next 5,000 years agricultures developed independently in at least six regions of the world. It is virtually certain that it was not a chance occurrence that so many new agricultures appeared in the same 5,000 years. What inhibited agriculture world wide for 44,000 years and what changed ～ 10,000 years ago? Here we suggest that a major factor influencing the development of agricultural societies was climate stability. From the experience of four cultures we estimate that the development of agriculture needed ～ 2,000 years of climate free from significant climate variations on time scales of a few centuries. Using the Empirical Mode Decomposition technique specifically designed to exhibit the time history of the amplitude of variations in non-stationary time series such as climate proxy records, we find that between 50,000 years ago and the termination of the Younger Dryas ～ 11,600 years ago there was probably no time span as long as 2,000 years that was free of relatively large century scale variations. Furthermore variations on these time scales appear to have been relatively small since the Younger Dryas (YD) ended, supporting our proposition concerning the importance of climate stability in the history of human culture.     

Late Pliocene to Pleistocene sensitivity of the Greenland Ice Sheet in response to external forcing and internal feedbacks

The timing and nature of ice sheet variations on Greenland over the last ～5 million years remain largely uncertain. Here, we use a coupled climate-vegetation-ice sheet model to determine the climatic sensitivity of Greenland to combined sets of external forcings and internal feedbacks operating on glacial-interglacial timescales. In particular, we assess the role of atmospheric pCO₂, orbital forcing, and vegetation dynamics in modifying thresholds for the onset of glaciation in late Pliocene and Pleistocene. The response of circum-Arctic vegetation to declining levels of pCO₂ (from 400 to 200 ppmv) and decreasing summer insolation includes a shift from boreal forest to tundra biomes, with implications for the surface energy balance. The expansion of tundra amplifies summer surface cooling and heat loss from the ground, leading to an expanded summer snow cover over Greenland. Atmospheric and land surface fields respond to forcing most prominently in late spring-summer and are more sensitive at lower Pleistocene-like levels of pCO₂. We find cold boreal summer orbits produce favorable conditions for ice sheet growth, however simulated ice sheet extents are highly dependent on both background pCO₂ levels and land-surface characteristics. As a result, late Pliocene ice sheet configurations on Greenland differ considerably from late Pleistocene, with smaller ice caps on high elevations of southern and eastern Greenland, even when orbital forcing is favorable for ice sheet growth.     

前寒武纪气候演化中的三个重要科学问题

自地球形成至寒武纪将近40亿年（距今46亿～5.4亿年,通常称为前寒武纪）的气候演变是一个具有特殊难度和挑战性的研究领域,同时也是基础和前沿的研究领域。文章选择了前寒武纪气候演化中的三个重要科学问题进行综述：大气演化、两次全球性的冰川期以及暗弱太阳问题。关于大气演化,本文首先描述了大气成分的演化历史,然后简述了影响大气成分演化的三个基本过程：大气逃逸、两次大气氧含量突然增加、碳酸盐-硅酸盐循环及其对气候系统的负反馈作用。两次全球性的冰川期分别发生在古元古代（距今24亿～21亿年）和新元古代（距今8亿～5.8亿年）,文章简述了其成因以及相关的气候模拟结果。暗弱太阳问题是地球历史气候演化的一个经典问题,论文简要地综述了一些最新的研究成果和观点。     

The global response to Younger Dryas boundary conditions in an AGCM simulation

 Geological evidence points to a global Younger Dryas (YD) climatic oscillation during the last glacial/ present interglacial transition phase. A convincing mechanism to explain this global YD climatic oscillation is not yet available. Nevertheless, a profound understanding of the mechanism behind the YD climate would lead to a better understanding of climate variability. Therefore, the Hamburg atmospheric circulation model was used to perform four numerical experiments on the YD climate. The objective of this study is to improve the understanding of different forcings influencing climate during the last glacial/interglacial transition and to investigate to what extent the model response agrees with global geological evidence of YD climate change. The following boundary conditions were altered: sea surface conditions, ice sheets, insolation and atmospheric CO₂ concentration. Sea surface temperatures based on foraminiferal assemblages proved to produce insufficient winter cooling in the N Atlantic Ocean in two experiments. It is proposed that this discrepancy is caused by uncertainties in the reconstruction method of sea surface temperatures. Therefore, a model-derived set of Atlantic surface ocean conditions was prescribed in a subsequent simulation. However, the latter set represented an Atlantic Ocean without a thermohaline circulation, which is not in agreement with evidence from ocean cores. The global response to the boundary conditions was analysed using three variables, namely surface temperature, zonal wind speed and precipitation. The statistical significance of the changes was tested with a two-tailed t-test. Moreover, the significant responses to cooled oceans were compared with geological evidence of a YD oscillation. This comparison revealed a good match in Europe, Greenland, Atlantic Canada and the N Pacific region, explaining the YD oscillation in these regions as a response to cooled N Atlantic and N Pacific Oceans. However, the results leave the YD climate in other regions completely unexplained. This reflects either an insufficient set of boundary conditions or the important role played by feedbacks within the coupled atmosphere-ocean-ice system. These feedbacks are poorly represented in the used atmospheric model, since ice sheets and the ocean surface conditions have to be prescribed. Received: 30 July 1996 / Accepted: 12 February 1997