Synoptic-scale perturbations in AGCM simulations of the present and Last Glacial Maximum climates

 The conditions of development of mid-latitude depressions (synoptic eddies) in the winter Northern Hemisphere mid-latitudes at the Last Glacial Maximum (LGM, 21 000 years ago) are very different from the present ones: this period is characterised by a general cooling of the extra-tropics, with massive ice sheets over the Northern Hemisphere continents and sea-ice extending very far south over the North Atlantic. The present work uses regression analysis to study the characteristics of the synoptic eddies in present-day and LGM climate simulations by the Atmospheric General Circulation Model (AGCM) of the UK Universities' Global Atmospheric Programme (UGAMP). In the LGM experiment, the structure of the Pacific eddies is similar to the present-day (PD) situation, but they are weaker. On the other hand, the Atlantic eddies show an increased zonal wavelength and a much shallower structure in the temperature and vertical wind perturbations. To understand the changes of these characteristics from present-day to LGM, we compare them to those computed for the most unstable modes of the corresponding mean flows, determined using a dry primitive equation model. A normal-mode stability analysis is carried both on zonally symmetric and asymmetric flows for each of the Northern Hemisphere storm-tracks. The changes in the most unstable normal modes found by both these analyses give a good account of changes in the structure of the perturbations as retrieved from the AGCM, suggesting that changes in the mean state (especially the temperature gradient) is the main driver of these changes. However in the case of the present-day Atlantic storm-track, the growth rate of these modes is found to be very low compared to the other cases. A complementary analysis evaluates the importance of non-modal growth, in the form of downstream development of perturbations, for each of the storm-tracks. This type of growth is found to be especially important in the case of the present-day Atlantic storm-track. Received: 29 September 1999 / Accepted: 17 November 1999     

Simulations of the Last Glacial Maximum climates using a general circulation model: prescribed versus computed sea surface temperatures

 The climate during the Last Glacial Maximum (LGM) has been simulated using the UK Universities Global Atmospheric Modelling Programme (UGAMP) general circulation model (GCM) with both prescribed sea surface temperatures (SSTs) based on the CLIMAP reconstruction and computed SSTs with a simple thermodynamic slab ocean. Consistent with the Paleoclimate Modelling Intercomparison Project (PMIP), the other boundary conditions include the large changes in ice-sheet topography and geography, a lower sea level, a lower concentration of CO₂ in the atmosphere, and a slightly different insolation pattern at the top of the atmosphere. The results are analysed in terms of changes in atmospheric circulation. Emphasis is given to the changes in surface temperatures, planetary waves, storm tracks and the associated changes in distribution of precipitation. The model responds in a similar manner to the changes in boundary conditions to previous studies in global mean statistics, but differs in its treatment of regional climates. Results also suggest that both the land ice sheets and sea ice introduce significant changes in planetary waves and transient eddy activity, which in turn affect regional climates. The computed SST simulations predict less sea ice and cooler tropical temperatures than those based on CLIMAP SSTs. It is unclear as to whether this is a model and/or a data problem, but the resulting changes in land temperatures and precipitation can be large. Snow mass budget analysis suggests that there is net ice loss along the southern edges of the Laurentide and Fennoscandian ice sheets and net ice gain over other parts of the two ice sheets. The net accumulation is mainly due to the decrease in ablation in the cold climate rather than to the changes in snowfall. The characteristics of the Greenland ice-sheet mass balance in the LGM simulations is also quite different from those in the present-day (PD) simulations. The ablation in the LGM simulations is negligible while it is a very important process in the ice mass budget in the PD simulations. Received: 10 January 1997 / Accepted: 11 December 1997     

Tropical paleoclimates at the Last Glacial Maximum: comparison of Paleoclimate Modeling Intercomparison Project (PMIP) simulations and paleodata

 Seventeen simulations of the Last Glacial Maximum (LGM) climate have been performed using atmospheric general circulation models (AGCM) in the framework of the Paleoclimate Modeling Intercomparison Project (PMIP). These simulations use the boundary conditions for CO₂, insolation and ice-sheets; surface temperatures (SSTs) are either (a) prescribed using CLIMAP data set (eight models) or (b) computed by coupling the AGCM with a slab ocean (nine models). The present-day (PD) tropical climate is correctly depicted by all the models, except the coarser resolution models, and the simulated geographical distribution of annual mean temperature is in good agreement with climatology. Tropical cooling at the LGM is less than at middle and high latitudes, but greatly exceeds the PD temperature variability. The LGM simulations with prescribed SSTs underestimate the observed temperature changes except over equatorial Africa where the models produce a temperature decrease consistent with the data. Our results confirm previous analyses showing that CLIMAP (1981) SSTs only produce a weak terrestrial cooling. When SSTs are computed, the models depict a cooling over the Pacific and Indian oceans in contrast with CLIMAP and most models produce cooler temperatures over land. Moreover four of the nine simulations, produce a cooling in good agreement with terrestrial data. Two of these model results over ocean are consistent with new SST reconstructions whereas two models simulate a homogeneous cooling. Finally, the LGM aridity inferred for most of the tropics from the data, is globally reproduced by the models with a strong underestimation for models using computed SSTs. Received: 9 September 1998 / Accepted: 18 March 1999     

Changes in fire regimes since the Last Glacial Maximum: an assessment based on a global synthesis and analysis of charcoal data

Fire activity has varied globally and continuously since the last glacial maximum (LGM) in response to long-term changes in global climate and shorter-term regional changes in climate, vegetation, and human land use. We have synthesized sedimentary charcoal records of biomass burning since the LGM and present global maps showing changes in fire activity for time slices during the past 21,000 years (as differences in charcoal accumulation values compared to pre-industrial). There is strong broad-scale coherence in fire activity after the LGM, but spatial heterogeneity in the signals increases thereafter. In North America, Europe and southern South America, charcoal records indicate less-than-present fire activity during the deglacial period, from 21,000 to ∼11,000 cal yr BP. In contrast, the tropical latitudes of South America and Africa show greater-than-present fire activity from ∼19,000 to ∼17,000 cal yr BP and most sites from Indochina and Australia show greater-than-present fire activity from 16,000 to ∼13,000 cal yr BP. Many sites indicate greater-than-present or near-present activity during the Holocene with the exception of eastern North America and eastern Asia from 8,000 to ∼3,000 cal yr BP, Indonesia and Australia from 11,000 to 4,000 cal yr BP, and southern South America from 6,000 to 3,000 cal yr BP where fire activity was less than present. Regional coherence in the patterns of change in fire activity was evident throughout the post-glacial period. These complex patterns can largely be explained in terms of large-scale climate controls modulated by local changes in vegetation and fuel load. The readers are requested to refer to the section “List of contributors” for the complete list of author affiliation details.     

The Last Glacial Maximum climate over Europe and western Siberia: a PMIP comparison between models and data

 Under the framework of the Palaeoclimate Modelling Intercomparison Project (PMIP), 17 climate models, 16 of which are atmospheric general circulation models, have been run to simulate the climate of the Last Glacial Maximum (21 000 years ago) using the same set of boundary conditions. Parallel to these numerical experiments, new, consistent, data bases have been developed on a continental scale. The present work compares the range of the model responses to the large perturbation corresponding to the conditions of the Last Glacial Maximum with consistently derived climate reconstructions from pollen records over Europe and western Siberia. It accounts for the differences in the model results due to the models themselves and directly compares this “error bar” due to the models to the uncertainties in the climate reconstructions from the pollen records. Overall the Last Glacial Maximum climate simulated by the models over western Europe is warmer, especially in winter, and wetter than the one depicted by the reconstructions. This is the region where the reconstructed increase in temperature, precipitation and moisture index from the Last Glacial Maximum to the present conditions is largest. The same disagreement, but of smaller amplitude, is found over Central Europe and the eastern Mediterranean Basin, while models and data are in broad agreement over western Siberia. The numerous modelling results allow a study of the link between the changes in atmospheric circulation and those in temperature, and an interpretation of the discrepancies in precipitation in terms of those in temperature. Received: 1 February 2000 / Accepted: 9 May 2000     

Comparing modeled and observed changes in mineral dust transport and deposition to Antarctica between the Last Glacial Maximum and current climates

Mineral dust aerosols represent an active component of the Earth’s climate system, by interacting with radiation directly, and by modifying clouds and biogeochemistry. Mineral dust from polar ice cores over the last million years can be used as paleoclimate proxy, and provide unique information about climate variability, as changes in dust deposition at the core sites can be due to changes in sources, transport and/or deposition locally. Here we present results from a study based on climate model simulations using the Community Climate System Model. The focus of this work is to analyze simulated differences in the dust concentration, size distribution and sources in current climate conditions and during the Last Glacial Maximum at specific ice core locations in Antarctica, and compare with available paleodata. Model results suggest that South America is the most important source for dust deposited in Antarctica in current climate, but Australia is also a major contributor and there is spatial variability in the relative importance of the major dust sources. During the Last Glacial Maximum the dominant source in the model was South America, because of the increased activity of glaciogenic dust sources in Southern Patagonia-Tierra del Fuego and the Southernmost Pampas regions, as well as an increase in transport efficiency southward. Dust emitted from the Southern Hemisphere dust source areas usually follow zonal patterns, but southward flow towards Antarctica is located in specific areas characterized by southward displacement of air masses. Observations and model results consistently suggest a spatially variable shift in dust particle sizes. This is due to a combination of relatively reduced en route wet removal favouring a generalized shift towards smaller particles, and on the other hand to an enhanced relative contribution of dry coarse particle deposition in the Last Glacial Maximum.     

Impact of vegetation changes on the dynamics of the atmosphere at the Last Glacial Maximum

Much work is under way to identify and quantify the feedbacks between vegetation and climate. Palaeoclimate modelling may provide a mean to address this problem by comparing simulations with proxy data. We have performed a series of four simulations of the Last Glacial Maximum (LGM, 21,000 years ago) using the climate model HadSM3, to test the sensitivity of climate to various changes in vegetation: a global change (according to a previously discussed simulation of the LGM with HadSM3 coupled to the dynamical vegetation model TRIFFID); a change only north of 35°N; a change only south of 35°N; and a variation in stomatal opening induced by the reduction in atmospheric CO₂ concentration. We focus mainly on the response of temperature, precipitation, and atmosphere dynamics. The response of continental temperature and precipitation mainly results from regional interactions with vegetation. In Eurasia, particularly Siberia and Tibet, the response of the biosphere substantially enhances the glacial cooling through a positive feedback loop between vegetation, temperature, and snow-cover. In central Africa, the decrease in tree fraction reduces the amount of precipitation. Stomatal opening is not seen to play a quantifiable role. The atmosphere dynamics, and more specifically the Asian summer monsoon system, are significantly altered by remote changes in vegetation: the cooling in Siberia and Tibet act in concert to shift the summer subtropical front southwards, weaken the easterly tropical jet and the momentum transport associated with it. By virtue of momentum conservation, these changes in the mid-troposphere circulation are associated with a slowing of the Asian summer monsoon surface flow. The pattern of moisture convergence is slightly altered, with moist convection weakening in the western tropical Pacific and strengthening north of Australia.     

关于末次盛冰期青藏高原大范围冰盖存在可能性的再研究

通过4组数值敏感性试验,利用BIOME3生态模式对末次盛冰期中国大陆植被分布状况进行了数值模拟,再次研究了该时期青藏高原大范围冰盖存在的可能性问题.结果表明,末次盛冰期大陆植被发生了大规模的变化,总体呈退化趋势;在降温5℃、降水减少10%、二氧化碳下降145×10-6和地球轨道参数轻微变化的情况下,青藏高原约1/2的地区已经被冰覆盖;进一步的降温和干燥则会加大高原冰的覆盖面积和大陆其余地区植被的退化.     

Last Glacial Maximum climate of the former Soviet Union and Mongolia reconstructed from pollen and plant macrofossil data

 An improved concept of the best analogues method was used to reconstruct the Last Glacial Maximum (LGM) climate from a set of botanical records from the former Soviet Union and Mongolia. Terrestrial pollen and macrofossil taxa were grouped into broad classes – plant functional types (PFTs), defined by the ecological and climatic parameters used in the BIOME1 model. PFT scores were then calibrated in terms of modern climate using 1245 surface pollen spectra from Eurasia and North America. In contrast to individual taxa, which exhibit great variability and may not be present in the palaeoassemblages, even in suitable climates, PFTs are more characteristic of the vegetation types. The modified method thus allows climate reconstruction at time intervals with partial direct analogues of modern vegetation (e.g. the LGM). At 18 kBP, mean temperatures were 20–29 °C colder than today in winter and 5–11 °C colder in summer in European Russia and Ukraine. Sites from western Georgia show negative, but moderate temperature anomalies compared to today: 8–11 °C in January and 5–7 °C in July. LGM winters were 7–15 °C colder and summers were 1–7 °C colder in Siberia and Mongolia. Annual precipitation sums were 50–750 mm lower than today across northern Eurasia, suggesting a weakening of the Atlantic and Pacific influences. Reconstructed drought index shows much drier LGM conditions in northern and mid-latitude Russia, but similar to or slightly wetter than today around the Black Sea and in Mongolia, suggesting compensation of precipitation losses by lower-than-present evaporation. Received: 11 May 1998 / Accepted: 25 September 1998     

GCM simulations of the Last Glacial Maximum surface climate of Greenland and Antarctica

 The LMDz variable grid GCM was used to simulate the Last Glacial Maximum (LGM, 21 ky Bp.) climate of Greenland and Antarctica at a spatial resolution of about 100 km.The high spatial resolution allows to investigate the spatial variability of surface climate change signals, and thus to address the question whether the sparse ice core data can be viewed as representative for the regional scale climate change. This study addresses primarily surface climate parameters because these can be checked against the, limited, ice core record. The changes are generally stronger for Greenland than for Antarctica, as the imposed changes of the forcing boundary conditions (e.g., sea surface temperatures) are more important in the vicinity of Greenland. Over Greenland, and to a limited extent also in Antarctica, the climate shows stronger changes in winter than in summer. The model suggests that the linear relationship between the surface temperature and inversion strength is modified during the LGM. The temperature dependency of the moisture holding capacity of the atmosphere alone cannot explain the strong reduction in snowfall over central Greenland; atmospheric circulation changes also play a crucial role. Changes in the high frequency variability of snowfall, atmospheric pressure and temperature are investigated and possible consequences for the interpretation of ice core records are discussed. Using an objective cyclone tracking scheme, the importance of changes of the atmospheric dynamics off the coasts of the ice sheets, especially for the high frequency variability of surface climate parameters, is illustrated. The importance of the choice of the LGM ice sheet topography is illustrated for Greenland, where two different topographies have been used, yielding results that differ quite strongly in certain nontrivial respects. This means that the paleo-topography is a significant source of uncertainty for the modelled paleoclimate. The sensitivity of the Greenland LGM climate to the prescribed sea surface conditions is examined by using two different LGM North Atlantic data sets. Received: 23 October 1997 / Accepted: 17 March 1998     

A coupled climate model simulation of the Last Glacial Maximum,Part 2: approach to equilibrium

The climate of the last glacial maximum (LGM) is simulated with a coupled climate model. The simulated climate undergoes a rapid adjustment during the first several decades after imposition of LGM boundary conditions, as described in Part 1, and then evolves toward equilibrium over 900 model years. The climate simulated by the coupled model at this period is compared with observationally-based LGM reconstructions and with LGM results obtained with an atmosphere-mixed layer (slab) ocean version of the model in order to investigate the role of ocean dynamics in the LGM climate. Global mean surface air temperature and sea surface temperature (SST) decrease by about 10 °C and 5.6 °C in the coupled model which includes ocean dynamics, compared to decreases of 6.3 and 3.8 °C in slab ocean case. The coupled model simulates a cooling of about 6.5 °C over the tropics, which is larger than that of the CLIMAP reconstruction (1.7 °C) and larger than that of the slab ocean simulation (3.3 °C), but which is in reasonable agreement with some recent proxy estimates. The ocean dynamics of the coupled model captures features found in the CLIMAP reconstructions such as a relative maximum of ocean cooling over the tropical Pacific associated with a mean La Niña-like response and lead to a more realistic SST pattern than in the slab model case. The reduction in global mean precipitation simulated in the coupled model is larger (15%) than that simulated with the slab ocean model (~10%) in conjunction with the enhanced cooling. Some regions, such as the USA and the Mediterranean region, experience increased precipitation in accord with proxy paleoclimate evidence. The overall much drier climate over the ocean leads to higher sea surface salinity (SSS) in most ocean basins except for the North Atlantic where SSS is considerably lower due to an increase in the supply of fresh water from the Mississippi and Amazon rivers and presumably a decrease in salt transport by the weakened North Atlantic overturning circulation. The North Atlantic overturning stream function weakens to less than half of the control run value. The overturning is limited to a shallower depth (less than 1000 m) and its outflow is confined to the Northern Hemisphere. In the Southern Ocean, convection is much stronger than in the control run leading to a stronger overturning stream function associated with enhanced Antarctic Bottom Water formation. As a result, Southern Ocean water masses fill the entire deep ocean. The Antarctic Circumpolar Current (ACC) transport through the Drake Passage increases by about 25%. The ACC transport, despite weaker zonal winds, is enhanced due to changes in bottom pressure torque. The weakening of the overturning circulation in the North Atlantic and the accompanying 30% decrease in the poleward ocean heat transport contrasts with the strengthening of the overturning circulation in the Southern Ocean and a 40% increase in heat transport. As a result, sea ice coverage and thickness are affected in opposite senses in the two hemispheres. The LGM climate simulated by the coupled model is in reasonable agreement with paleoclimate proxy evidence. The dynamical response of the ocean in the coupled model plays an important role in determining the simulated, and undoubtedly, the actual, LGM climate.     

Insolation, Moisture Balance and Climate Change on the South American Altiplano Since the Last Glacial Maximum

Sediment cores from Lake Titicaca contain proxy records of past lake level and hydrologic change on the South American Altiplano. Large downcore shifts in the isotopic composition of organic carbon, C/N, wt.%C_org, %CaCO₃, and % biogenicsilica illustrate the dynamic changes in lake level that occurred during the past 20,000 years. The first cores taken from water depths greater than 50 meters in the northern subbasin of the lake are used to develop and extend the paleolake-level record back to the Last Glacial Maximum (LGM). Quantitative estimates of lake level are developed using transfer functions based on the ¹³C of modern lacustrine organic sources and the ¹³C of modern sedimented organic matter from core-tops. Lake level was slightly higher than modern during much of the post-LGM (20,000–13,500 yr BP) and lake water was freshunder the associated outflow conditions. The Pleistocene/Holocene transition (13,500–7,500 yr BP) was a period of gradual regression, punctuated by minor trangressions. Following a brief highstand at about 7250 yr BP, lake level dropped rapidly to 85 m below the modern level, reaching maximum lowstand conditions by 6250 yr BP. Lake level increased rapidly between 5000yr BP and 4000 yr BP, and less rapidly between 4000 yr BP and 1500 yr BP.Lake level remained relatively high throughout the latest Holocene with only minor fluctuations (<12 meters). Orbitally induced changes in solar insolation, coupled with long-term changes in El Niño-Southern Oscillation variability, are the most likely driving forces behind millennial-scale shifts in lake level that reflect regional-scale changes in the moisture balance of the Atlantic-Amazon-Altiplano hydrologic system.     

Glacial meltwater cooling of the Gulf of Mexico: GCM implications for Holocene and present-day climates

An atmospheric general circulation model, the NCAR CCM, has been used to investigate the possible effects that reduced Gulf of Mexico sea surface temperatures (SST) could have on regional and hemispheric climates. ¹⁸O records and terrestrial evidence indicate at least two major glacial meltwater discharges into the Gulf of Mexico subsequent to the last glacial maximum. It is probable that these discharges reduced Gulf of Mexico SST. We have conducted three numerical experiments, with imposed gulf-wide SST coolings of 3°C, 6°C, and 12°C, and find in all three experiments significant reductions in the North Atlantic storm-track intensity, along with a strong decrease in transient eddy water vapor transport out of the Gulf of Mexico. Surface pressures are higher over the North Atlantic, indicating a reduction of the climatological Icelandic low. The region is generally cooler and drier, with a reduction in precipitation that agrees well with evidence from Greenland ice cores. Other statistically significant changes occur across the Northern Hemisphere, but vary between the three experiments. In particular, warmer, wetter conditions are found over Europe for both the 6°C and 12°C SST reductions, but cooler conditions are found for the 3°C reduction. This indicates a dependence, in both the sign and magnitude of the model response, on the magnitude of the imposed SST anomaly. The results suggest that the present-day North Atlantic storm track is dependent on warm Gulf of Mexico SST for much of its intensity. They also suggest that meltwater-induced coolings may help account, in part, for some of the climatic oscillations that occurred during the last glacial/interglacial transition.     

Climate change impacts on ecosystems and the terrestrial carbon sink: a new assessment

Climate output from the UK Hadley Centre's HadCM2 and HadCM3 experiments for the period 1860 to 2100, with IS92a greenhouse gas forcing, together with predicted patterns of N deposition and increasing CO₂, were input (offline) to the dynamic vegetation model, Hybrid v4.1 (Friend et al., 1997; Friend and White, 1999). This model represents biogeochemical, biophysical and biogeographical processes, coupling the carbon, nitrogen and water cycles on a sub-daily timestep, simulating potential vegetation and transient changes in annual growth and competition between eight generalized plant types in response to climate.Global vegetation carbon was predicted to rise from about 600 to 800 PgC (or to 650 PgC for HadCM3) while the soil carbon pool of about 1100 PgC decreased by about 8%. By the 2080s, climate change caused a partial loss of Amazonian rainforest, C₄ grasslands and temperate forest in areas of southern Europe and eastern USA, but an expansion in the boreal forest area. These changes were accompanied by a decrease in net primary productivity (NPP) of vegetation in many tropical areas, southern Europe and eastern USA (in response to warming and a decrease in rainfall), but an increase in NPP of boreal forests. Global NPP increased from 45 to 50 PgC y⁻¹ in the 1990s to about 65 PgC y⁻¹ in the 2080s (about 58 PgC y⁻¹ for HadCM3). Global net ecosystem productivity (NEP) increased from about 1.3 PgC y⁻¹ in the 1990s to about 3.6 PgC y⁻¹ in the 2030s and then declined to zero by 2100 owing to a loss of carbon from declining forests in the tropics and at warm temperate latitudes — despite strengthening of the carbon sink at northern high latitudes. HadCM3 gave a more erratic temporal evolution of NEP than HadCM2, with a dramatic collapse in NEP in the 2050s.