Assessing the sensitivity of the North Atlantic Ocean circulation to freshwater perturbation in various glacial climate states

A striking characteristic of glacial climate in the North Atlantic region is the recurrence of abrupt shifts between cold stadials and mild interstadials. These shifts have been associated with abrupt changes in Atlantic Meridional Overturning Circulation (AMOC) mode, possibly in response to glacial meltwater perturbations. However, it is poorly understood why they were more clearly expressed during Marine Isotope Stage 3 (MIS3, ~60?C27?ka BP) than during Termination 1 (T1, ~18?C10?ka BP) and especially around the Last Glacial Maximum (LGM, ~23?C19?ka BP). One clue may reside in varying climate forcings, making MIS3 and T1 generally milder than LGM. To investigate this idea, we evaluate in a climate model how ice sheet size, atmospheric greenhouse gas concentration and orbital insolation changes between 56?ka BP (=56k), 21k and 12.5k affect the glacial AMOC response to additional freshwater forcing. We have performed three ensemble simulations with the earth system model LOVECLIM using those forcings. We find that the AMOC mode in the mild glacial climate type (56k and 12.5k), with deep convection in the Labrador Sea and the Nordic Seas, is more sensitive to a constant 0.15?Sv freshwater forcing than in the cold type (21k), with deep convection mainly south of Greenland and Iceland. The initial AMOC weakening in response to freshwater forcing is larger in the mild type due to an early shutdown of Labrador Sea deep convection, which is completely absent in the 21k simulation. This causes a larger fraction of the freshwater anomaly to remain at surface in the mild type compared to the cold type. After 200?years, a weak AMOC is established in both climate types, as further freshening is compensated by an anomalous salt advection from the (sub-)tropical North Atlantic. However, the slightly fresher sea surface in the mild type facilitates further weakening of the AMOC, which occurs when a surface buoyancy threshold (?0.6?kg?m^?3 surface density anomaly to the 56k reference state) is stochastically crossed in the Nordic Seas. While described details are model-specific, our results imply that a more northern location of deep convection sites during milder glacial times may have amplified frequency and amplitude of abrupt climate shifts.     

Mid-Holocene monsoons: a multi-model analysis of the inter-hemispheric differences in the responses to orbital forcing and ocean feedbacks

The response of monsoon circulation in the northern and southern hemisphere to 6?ka orbital forcing has been examined in 17 atmospheric general circulation models and 11 coupled ocean–atmosphere general circulation models. The atmospheric response to increased summer insolation at 6?ka in the northern subtropics strengthens the northern-hemisphere summer monsoons and leads to increased monsoonal precipitation in western North America, northern Africa and China; ocean feedbacks amplify this response and lead to further increase in monsoon precipitation in these three regions. The atmospheric response to reduced summer insolation at 6?ka in the southern subtropics weakens the southern-hemisphere summer monsoons and leads to decreased monsoonal precipitation in northern South America, southern Africa and northern Australia; ocean feedbacks weaken this response so that the decrease in rainfall is smaller than might otherwise be expected. The role of the ocean in monsoonal circulation in other regions is more complex. There is no discernable impact of orbital forcing in the monsoon region of North America in the atmosphere-only simulations but a strong increase in precipitation in the ocean–atmosphere simulations. In contrast, there is a strong atmospheric response to orbital forcing over northern India but ocean feedback reduces the strength of the change in the monsoon although it still remains stronger than today. Although there are differences in magnitude and exact location of regional precipitation changes from model to model, the same basic mechanisms are involved in the oceanic modulation of the response to orbital forcing and this gives rise to a robust ensemble response for each of the monsoon systems. Comparison of simulated and reconstructed changes in regional climate suggest that the coupled ocean–atmosphere simulations produce more realistic changes in the northern-hemisphere monsoons than atmosphere-only simulations, though they underestimate the observed changes in precipitation in all regions. Evaluation of the southern-hemisphere monsoons is limited by lack of quantitative reconstructions, but suggest that model skill in simulating these monsoons is limited.     

Simulation of paleoclimate over East Asia at 6 ka BP and 21 ka BP by a regional climate model

Using a regional climate model with detailed land surface processes (RegCM2), East Asian monsoon climates at 6 ka BP and 21 ka BP are simulated by prescribing vegetation and employing paleovegetation respectively in order to examine land surface effects on East Asian climate system and the potential mechanisms for climate change. The RegCM2 with a 120 × 120 km² resolution has simulated the enlargement of the seasonal cycle of insolation, the temperature rising the whole year, and the reduction of perpetual snow in high latitudes at 6 ka BP. The simulation shows the East Asian summer monsoon strengthening, precipitation and P–E increasing, and the monsoon rain belt shifting westwards and northwards. Effect of paleovegetation included in the modeling reduced surface albedo and caused an increase in the winter temperature, which led to weakening of the winter continental cold anticyclone over China. The results make the seasonal characteristics of simulated temperature changes in better agreement with the geological records, and are an improvement over previous simulations of Paleoclimate Modeling Intercomparison Project (PMIP). The RegCM2 simulated the 21 ka BP climate with lowered temperature throughout the year, and with precipitation reduced in most areas of East Asia (but increased in both the Tibetan Plateau and Central Asia). Low temperature over East Asia led to the strengthening of the East Asian winter monsoon and the shrinking of the summer monsoon. The effect of paleovegetation included in the experiment has enlarged the glacial climate influence in East Asia, which is closer to geological data than the PMIP simulations directly driven by insolation, glaciation and low CO₂ concentration.     

AN AGCM +SSiB MODEL SIMULATION ON CHANGES IN PALAEOMONSOON CLIMATE AT 21 KA BP IN CHINA*

The numerical simulation experiment of climate at Last Glacial Maximum (LGM.21 ka BP) in China is made by using an atmospheric general circulation model (AGCM) coupled with land surface processes (AGCM+SSiB) and earth orbital parameters and boundary forcing conditions at21 ka.The modeled climate features are compared with reconstructed conditions at 21 ka from paleo-lake data and pollen data.The results show that the simulated climate conditions at 21 ka in China are fairly comparable with paleo-climatological data.The climate features at 21 ka in China from the experiment are characterized by a drier in the east and a wetter in the west and in the Tibetan Plateau as well.According to the analysis of distribution of pressure and precipitation,as well as the intensity of atmospheric circulation at 21 ka,monsoon circulation in eastern Asia was significantly weak comparing with the present.In the Tibetan Plateau,the intensity of summer monsoon circulation was strengthened,and winter monsoon was a little stronger than the present.The simulation with given forcing boundary conditions,especially the different vegetation coverage,can reproduce the climate condition at the LGM in China,and therefore provides dynamical mechanisms on the climate changes at 21 ka.     

AN AGCM SSiB MODEL SIMULATION ON CHANGES IN PALAEOMONSOON CLIMATE AT 21 KA BP IN CHINA

The numerical simulation experiment of climate at Last Glacial Maximum (LGM.21 ka BP)in China is made by using an atmospheric general circulation model (AGCM) coupled with landsurface processes (AGCM SSiB) and earth orbital parameters and boundary forcing conditions at21 ka.The modeled climate features are compared with reconstructed conditions at 21 ka frompaleo-lake data and pollen data.The results show that the simulated climate conditions at 21 ka inChina are fairly comparable with paleo-climatological data.The climate features at 21 ka in Chinafrom the experiment are characterized by a drier in the east and a wetter in the west and in theTibetan Plateau as well.According to the analysis of distribution of pressure and precipitation,aswell as the intensity of atmospheric circulation at 21 ka,monsoon circulation in eastern Asia wassignificantly weak comparing with the present.In the Tibetan Plateau,the intensity of summermonsoon circulation was strengthened,and winter monsoon was a little stronger than the present.The simulation with given forcing boundary conditions,especially the different vegetationcoverage,can reproduce the climate condition at the LGM in China,and therefore providesdynamical mechanisms on the climate changes at 21 ka.     

Global monsoons in the mid-Holocene and oceanic feedback

The response of the six major summer monsoon systems (the North American monsoon, the northern Africa monsoon, the Asia monsoon, the northern Australasian monsoon, the South America monsoon and the southern Africa monsoon) to mid-Holocene orbital forcing has been investigated using a coupled ocean–atmosphere general circulation model (FOAM), with the focus on the distinct roles of the direct insolation forcing and oceanic feedback. The simulation result is also found to compare well with the NCAR CSM. The direct effects of the change in insolation produce an enhancement of the Northern Hemisphere monsoons and a reduction of the Southern Hemisphere monsoons. Ocean feedbacks produce a further enhancement of the northern Africa monsoon and the North American monsoon. However, ocean feedbacks appear to weaken the Asia monsoon, although the overall effect (direct insolation forcing plus ocean feedback) remains a strengthened monsoon. The impact of ocean feedbacks on the South American and southern African monsoons is relatively small, and therefore these regions, especially the South America, experienced a reduced monsoon regime compared to present. However, there is a strong ocean feedback on the northern Australian monsoon that negates the direct effects of orbital changes and results in a strengthening of austral summer monsoon precipitation in this region. A new synthesis is made for mid-Holocene paleoenvironmental records and is compared with the model simulations. Overall, model simulations produce changes in regional climates that are generally consistent with paleoenvironmental observations.     

用一个耦合的全球格点大气环流模式-植被模式模拟中全新世气候变化

用一个耦合的全球格点大气环流模式－植被模式模拟中全新世的气候变化，模拟试验中考虑了地球轨道参数的变化，而其他强迫条件均取成现今值。结果表明，耦合的模式能够模拟出较今强的大尺度夏季风，特别是亚洲－非洲季风，而其他季节和区域的变化值一般都比较小。季风环流和季风降水都大幅度地增大了。结果还显示，耦合模式模拟的大尺度季风系统的变化同单纯大气环流模式模拟的结果非常相似，但是，在非洲北部季风区耦合模式模拟的降水和温度变化较单纯大气模式模拟的值要大，而且，耦合模式模拟的冬季降温值要比单纯大气模式模拟的结果小。     

High frequency pulses of East Asian monsoon climate in the last two glaciations: link with the North Atlantic

Three loess sections in the Loess Plateau of China have been studied to characterize the variations of the East Asian monsoon climate in the Late Quaternary period. Paleo-weathering profiles based on two weathering indices reveal a series of spatially correlative intervals of lower weathering intensity (LW), in the last glacial and late penultimate glacial loess, indicating significantly cooler or/and drier conditions. Most of them are identified at or near the transitional boundaries between loess and soil units, suggesting that climatic conditions favorable for LW events tended to occur at or near major climatic boundaries. Twenty-eight radiocarbon and thermoluminescence measurements, combined with Kukla's magnetic susceptibility age model date these events at ∼14, ∼21, ∼29, ∼38, ∼52, ∼71, ∼135, and ∼145 ka BP, with ∼5–10 ka frequency inlaid within the orbitally induced ∼20-ka periodicity. The ages of the first six events in the last glacial period are therefore close to those of the coarse-grained Heinrich layers in the North Atlantic Ocean, which resulted from massive discharges of icebergs. The results indicate that the climate in the area affected by the East Asian monsoon has experienced high-frequency changes more or less synchronous with the Heinrich events. These high-frequency changes have been reported by Porter and An, based on the grain-size time series from the Luochuan loess section. Our data also reveal that similar events also occurred during the penultimate glaciation. The mechanisms linking the Heinrich events and the East Asian monsoon climate may be similar to that driving the Younger Dryas event which has been identified in a large number of records in China and in the surrounding oceans. Received: 16 October 1995 / Accepted: 10 May 1996     

Potential role of winter rainfall in explaining increased moisture in the Mediterranean and Middle East during periods of maximum orbitally-forced insolation seasonality

Precession-related forcing of seasonal insolation changes in the northern hemisphere (NH) alternates between maximum NH seasonality (summer perihelion–increased insolation; winter aphelion–decreased insolation) and minimum NH seasonality (summer aphelion, and winter perihelion). With maximum NH seasonality, climate models simulate stronger NH summer monsoons that bring increased precipitation to North Africa and South and East Asia, in agreement with the in-phase relation of precipitation and NH summer insolation found in many paleoclimatic records. However paleoclimatic records in parts of the Mediterranean, the Middle East, and the interior of Asia also indicate increased moisture at times of maximum NH seasonality, a change not always clearly linked to stronger summer monsoons—either because these regions are at or beyond the boundaries of the present-day monsoon or because the observations allow multiple causal interpretations, or both. This study focuses on the possible role of changes in NH winter climate in explaining these wetter episodes. Using climate model simulations, we show that the ‘NH winter aphelion–decreased NH winter insolation’ orbital configuration is linked to the Mediterranean storm track and increased winter rains in the Mediterranean, the Middle East, and interior Asia. We conclude that wetter periods at precession time scales in these particular regions may have resulted either from increased wintertime storm track precipitation, or from a combination of increased winter and summer rainfall. Given this seasonal ambiguity, both possibilities need to be considered.     

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.     

Multiple factors causing Holocene lake-level change in monsoonal and arid central Asia as identified by model experiments

Lake-level records provide a rich resource of information about past changes in effective moisture, but water-balance fluctuations can be driven by a number of different climate variables and it is often difficult to pinpoint their exact cause. This understanding is essential, however, for reconciling divergent paleo-records or for making predictions about future lake-level variations. This research uses a series of models, the NCAR CCSM3, a lake energy-balance and a lake water-balance model, to examine the reasons for lake-level changes in monsoonal Asia and arid central Asia between the early (8.5 ka), middle (6.0 ka) and late (ca. 1800 AD) Holocene. Our results indicate that the components of the lake water balance responsible for lake-level changes varied by region and through time. High lake levels at 8.5 and 6.0 ka in the monsoon region were caused by the combined effects of low lake evaporation and high precipitation. The low lake evaporation resulted from low winter solar radiation and high summer cloud cover. Precipitation associated with the mid-latitude westerlies increased from the early to middle Holocene and maintained high lake levels throughout most of arid central Asia ca. 6 ka. The modeled evolution of lake level in arid central Asia from the mid to late Holocene was spatially heterogeneous, due to different sensitivities of the northern and southern parts of the region to seasonally-changing insolation, particularly regarding the duration of lake ice cover. The model results do not suggest that precipitation and lake evaporation changes compete with one another in forcing lake-level change, as has been hypothesized.     

Modeling the time-dependent response of the Asian summer monsoon to obliquity forcing in a coupled GCM: a PHASEMAP sensitivity experiment

To study the time-dependent response of the Asian summer monsoon to obliquity forcing, we analyze a 284,000-year long transient simulation produced by a fully coupled global climate model (GCM) using a new phase mapping (PHASEMAP) approach. Here we focus on understanding the phase response of monsoonal circulation to insolation forcing at the Earth-orbital obliquity band (41 Kyr). Our results show that the East Asian summer monsoon (EASM) can be divided into two geographic regions: the North East Asian summer monsoon (NEASM) and the South East Asian summer monsoon (SEASM). The Indian summer monsoon (ISM) and the SEASM are in phase at the obliquity band, strengthened with an increase in obliquity from Obliquity minima (Omin) to Obliquity maxima (Omax). The NEASM is out of phase with the ISM and SEASM, weakened with an increase in obliquity from Omin to Omax. We hypothesize that the inverse phase between the NEASM and the ISM at the obliquity band results from an ISM–NEASM teleconnection linked to the formation mechanism of the Bonin High.     

The role of phase locking in a simple model for glacial dynamics

Glacial–interglacial oscillations are often described by simple conceptual models. Relatively few models, however, are accompanied by analytical solutions, though detailed analytical investigation of climate models often leads to deeper understanding of the climate system. Here we study a simple conceptual model for glacial dynamics, a simplified version of the sea-ice-switch mechanism of Gildor and Tziperman (Paleoceanography 15:605–615, 2000), and provide a detailed analytical treatment for this model. We show that when the model is forced by a simplified insolation forcing it exhibits rich dynamics and passes through a series of bifurcations before being completely phase-locked to the insolation forcing. Our model suggests that even when the glacial cycles are self-sustained, insolation forcing has a major role on the complexity of glacial cycles: (1) it is possible to obtain glacial–interglacial oscillations for a wider parameters range when the amplitude of the insolation forcing is larger; (2) in addition, the ice-volume becomes more periodic; (3) when the period of the ice-volume is minimal the ice-volume is symmetric and for larger period is more asymmetric; (4) the ice-volume can be either periodic, higher order periodic, or quasi-periodic.     

Impact of Earth’s orbit and freshwater fluxes on Holocene climate mean seasonal cycle and ENSO characteristics

We use a state-of-the-art 3-dimensional coupled model to investigate the relative impact of long term variations in the Holocene insolation forcing and of a freshwater release in the North Atlantic. We show that insolation has a greater effect on seasonality and La Ni?a events and is the major driver of sea surface temperature changes. In contrast, the variations in precipitation reflect changes in El Ni?o events. The impact of ice-sheet melting may have offset the impact of insolation on El Ni?o Southern Oscillation variability at the beginning of the Holocene. These simulations provide a coherent framework to refine the interpretation of proxy data and show that changes in seasonality may bias the projection of relationships established between proxy indicators and climate variations in the east Pacific from present day records.     

Simulation of the evolutionary response of global summer monsoons to orbital forcing over the past 280,000 years

We describe the evolutionary response of northern and southern hemisphere summer monsoons to orbital forcing over the past 280,000 years using a fully coupled general circulation ocean-atmosphere model in which the orbital forcing is accelerated by a factor of 100. We find a strong and positive response of northern (southern) summer monsoon precipitation to northern (southern) summer insolation forcing. On average, July (January) precipitation maxima and JJA (DJF) precipitation maxima have high coherence and are approximately in phase with June (December) insolation maxima, implying an average lag between forcing and response of about 30° of phase at the precession period. The average lag increases to over 40° for 4-month precipitation averages, JJAS (DJFM). The phase varies from region to region. The average JJA (DJF) land temperature maxima also lag the June orbital forcing maxima by about 30° of phase, whereas ocean temperature maxima exhibit a lag of about 60° of phase at the precession period. Using generalized measures of the thermal and hydrologic processes that produce monsoons, we find that the summer monsoon precipitation indices for the six regions all fall within the phase limits of the process indices for the respective hemispheres. Selected observational studies from four of the six monsoon regions report approximate in-phase relations of summer monsoon proxies to summer insolation. However other observational studies report substantial phase lags of monsoon proxies and a strong component of forcing associated with glacial-age boundary conditions or other factors. An important next step will be to include glacial-age boundary condition forcing in long, transient paleoclimate simulations, along with orbital forcing.     

Orbital forcing and role of the latitudinal insolation/temperature gradient

Orbital forcing of the climate system is clearly shown in the Earths record of glacial–interglacial cycles, but the mechanism underlying this forcing is poorly understood. Traditional Milankovitch theory suggests that these cycles are driven by changes in high latitude summer insolation, yet this forcing is dominated by precession, and cannot account for the importance of obliquity in the Ice Age record. Here, we investigate an alternative forcing based on the latitudinal insolation gradient (LIG), which is dominated by both obliquity (in summer) and precession (in winter). The insolation gradient acts on the climate system through differential solar heating, which creates the Earths latitudinal temperature gradient (LTG) that drives the atmospheric and ocean circulation. A new pollen-based reconstruction of the LTG during the Holocene is used to demonstrate that the LTG may be much more sensitive to changes in the LIG than previously thought. From this, it is shown how LIG forcing of the LTG may help explain the propagation of orbital signatures throughout the climate system, including the Monsoon, Arctic Oscillation and ocean circulation. These relationships are validated over the last (Eemian) Interglacial, which occurred under a different orbital configuration to the Holocene. We conclude that LIG forcing of the LTG explains many criticisms of classic Milankovitch theory, while being poorly represented in climate models.     

The Seasonal Climate and Low Frequency Oscillation in the Simulated Mid-Holocene Megathermal Climate

l.IntroductionThegoalofPMIP(PaleoclimateModellingIntercomparisonProject)istocomparethestate-ofthe-artclimatemodels'simulationsforthepastclimateonthebasisofthepaleoclimatedata.Oneofthemostinterestingperiodsisthemid-Holoceneduring8.5-3.OkaBPwiththemaximumin7.2-6.OkaBP(beforepresent)inChina(Shietal.,l992).Itwasre-vealedthatthetemPeraturewasl-4'Chigherandtheprecipitationwas4O%-loO%largercomparedtothepresentoverChina(Shietal.,l992;Anetal..l99l;Kongetal.,l99O,l99l).Therehavebeensomesimula…     

Calendar effect on phase study in paleoclimate transient simulation with orbital forcing

Several studies have shown that the use of different calendars in paleoclimate simulations can cause artificial phase shifts on insolation forcing and climatic responses. However, these important calendar corrections are still often neglected. In this paper, the phase shifts at the precession band is quantitatively assessed by converting the model data of the transient GCM climate simulation of Kutzbach et al. (Clim Dyn 30:567?C579, 2008) from the ??fixed-day?? calendar to the ??fixed-angular?? calendar with a new and efficient approach. We find that insolation has a big phase shift in September?COctober?CNovember (SON) when the vernal equinox (VE) is fixed to March 21. At high latitude, the phase bias is up to 60° (about 3650?years). The insolation phase bias in SON in Southern Hemisphere (SH) is especially important because it can influence the timing of the SH summer monsoon response due to the large heat capacity of ocean. The calendar correction has minor effect (±2°) on the phase relationships between forcing and precipitation responses of the six global summer monsoons studied in Kutzbach et al. (2008). After correcting the calendar effect, especial on SH ocean temperature, the new phase wheel results are more similar for both hemispheres. The results suggest that the calendar effect should be corrected before discussing the dynamics between orbital forcing and climatic responses in phase studies of transient simulations.     

The Holocene paleohydrology and climatic history of the northern great sandy desert - fitzroy trough: With special reference to the history of the northwest Australian monsoon

The Holocene history of the northwest Australian monsoon is considered in the light of new stratigraphic findings from the southern margin of the monsoon regime. A discussion of the present-day climatology of northwestern Australia and synoptic controls on precipitation in the areas at the southern margin of the monsoon regime is given as background information for the interpretation of stratigraphic findings from swamp and alluvial sequences. The stratigraphic findings point to little change in the intensity of the northwest Australian monsoon since ca. 6.5 ka. A more tentative inference drawn from the field evidence is that the monsoon activity may well have been reduced in this region during the Early Holocene; during this time insolation levels were also reduced. The possible relationship of the northwest Australian monsoon to insolation levels is explored. It is clear that monsoon behavior cannot be accounted for as a simple linear response of monsoon intensity to variations in insolation. Possible controls are discussed with particular emphasis on the plausibility of changes in the damping of convective motions providing the mechanism. Gill's model of heat-induced circulation is used to provide some evidence in favour of this conjecture.     

Fully coupled climate/dynamical vegetation model simulations over Northern Africa during the mid-Holocene

 The climate and vegetation patterns of the middle Holocene (6000 years ago; 6 ka) over Northern Africa are simulated using a fully-synchronous climate and dynamical vegetation model. The coupled model predicts a northward shift in tropical rainforest and tropical deciduous forest vegetation by about 5 degrees of latitude, and an increase in grassland at the present-day simulated Saharan boundaries. The northward expansion of vegetation over North Africa at 6 ka is initiated by an orbitally-induced amplification of the summer monsoon, and enhanced by feedback effects induced by the vegetation. These combined processes lead to a major reduction in Saharan desert area at 6 ka relative to present-day of about 50%. However, as shown in previous asynchronous modelling studies, the coupled climate/vegetation model does not fully reproduce the vegetation patterns inferred from palaeoenvironmental records, which suggest that steppe vegetation may have existed across most of Northern Africa. Orbital changes produce an intensification of monsoonal precipitation during the peak rainy season (July to September), whilst vegetation feedbacks, in addition to producing further increases in the peak intensity, play an important role in extending the rainy season from May/June through to November. The orbitally induced increases in precipitation are relatively uniform from west to east, in contrast to vegetation feedback-induced increases in precipitation which are concentrated in western North Africa. Annual-average precipitation increases caused by vegetation feedbacks are simulated to be of similar importance to orbital effects in the west, whilst they are relatively unimportant farther to the east. The orbital, vegetation and combined orbital and vegetation-induced changes in climate, from the simulations presented in this study, have been compared with results from previous modelling studies over the appropriate North African domain. Consequently, the important role of vegetation parametrizations in determining the magnitude of vegetation feedbacks has been illustrated. Further modelling studies which include the effects of changes in ocean temperature and changes in soil properties may be needed, along with additional observations, to resolve the discrepancy between model predictions of vegetation and palaeorecords for North Africa. Received: 15 June 1999 / Accepted: 14 December 1999