Palaeoclimate Modelling Intercomparison Project Research Articles

No changes in overall AMOC strength in interglacial PMIP4 time slices

Abstract. The Atlantic Meridional Overturning Circulation (AMOC) is a key mechanism of poleward heat transport and an important part of the global climate system. How it responded to past changes in forcing, such as those experienced during Quaternary interglacials, is an intriguing and open question. Previous modelling studies suggest an enhanced AMOC in the mid-Holocene compared to the preindustrial period. In earlier simulations from the Palaeoclimate Modelling Intercomparison Project (PMIP), this arose from feedbacks between sea ice and AMOC changes, which were dependent on resolution. Here we present an initial analysis of recently available PMIP4 simulations for three experiments representing different interglacial conditions – one 127 000 years ago within the Last Interglacial (127 ka, called lig127k), one in the middle of the Holocene (midHolocene, 6 ka), and a preindustrial control simulation (piControl, 1850 CE). Both lig127k and midHolocene have altered orbital configurations compared to piControl. The ensemble mean of the PMIP4 models shows the strength of the AMOC does not markedly change between the midHolocene and piControl experiments or between the lig127k and piControl experiments. Therefore, it appears orbital forcing itself does not alter the overall AMOC. We further investigate the coherency of the forced response in AMOC across the two interglacials, along with the strength of the signal, using eight PMIP4 models which performed both interglacial experiments. Only two models show a stronger change with the stronger forcing, but those models disagree on the direction of the change. We propose that the strong signals in these two models are caused by a combination of forcing and the internal variability. After investigating the AMOC changes in the interglacials, we further explored the impact of AMOC on the climate system, especially on the changes in the simulated surface temperature and precipitation. After identifying the AMOC's fingerprint on the surface temperature and rainfall, we demonstrate that only a small percentage of the simulated surface climate changes could be attributed to the AMOC. Proxy records of sedimentary Pa/Th ratio during the two interglacial periods both show a similar AMOC strength compared to the preindustrial, which fits nicely with the simulated results. Although the overall AMOC strength shows minimal changes, future work is required to explore whether this occurs through compensating variations in the different components of AMOC (such as Iceland–Scotland overflow water). This line of evidence cautions against interpreting reconstructions of past interglacial climate as being driven by AMOC, outside of abrupt events.

Analysing the PMIP4-CMIP6 collection: a workflow and tool (pmip_p2fvar_analyzer v1)

Abstract. Experiment outputs are now available from the Coupled Model Intercomparison Project's sixth phase (CMIP6) and the past climate experiments defined in the Palaeoclimate Modelling Intercomparison Project's fourth phase (PMIP4). All of this output is freely available from the Earth System Grid Federation (ESGF). Yet there is overhead in analysing this resource that may prove complicated or prohibitive. Here we document the steps taken by ourselves to produce ensemble analyses covering past and future simulations. We outline the strategy used to curate, adjust the monthly calendar aggregation and process the information downloaded from the ESGF. The results of these steps were used to perform analysis for several of the initial publications arising from PMIP4. We provide post-processed fields for each simulation, such as climatologies and common measures of variability. Example scripts used to visualise and analyse these fields are provided for several important case studies.

Response of Amazonian forests to mid-Holocene drought: A model-data comparison.

There is a major concern for the fate of Amazonia over the coming century in the face of anthropogenic climate change. A key area of uncertainty is the scale of rainforest dieback to be expected under a future, drier climate. In this study, we use the middle Holocene (ca. 6000years before present) as an approximate analogue for a drier future, given that palaeoclimate data show much of Amazonia was significantly drier than present at this time. Here, we use an ensemble of climate and vegetation models to explore the sensitivity of Amazonian biomes to mid-Holocene climate change. For this, we employ three dynamic vegetation models (JULES, IBIS, and SDGVM) forced by the bias-corrected mid-Holocene climate simulations from seven models that participated in the Palaeoclimate Modelling Intercomparison Project 3 (PMIP3). These model outputs are compared with a multi-proxy palaeoecological dataset to gain a better understanding of where in Amazonia we have most confidence in the mid-Holocene vegetation simulations. A robust feature of all simulations and palaeodata is that the central Amazonian rainforest biome is unaffected by mid-Holocene drought. Greater divergence in mid-Holocene simulations exists in ecotonal eastern and southern Amazonia. Vegetation models driven with climate models that simulate a drier mid-Holocene (100-150mm per year decrease) better capture the observed (palaeodata) tropical forest dieback in these areas. Based on the relationship between simulated rainfall decrease and vegetation change, we find indications that in southern Amazonia the rate of tropical forest dieback was ~125,000km2 per 100mm rainfall decrease in the mid-Holocene. This provides a baseline sensitivity of tropical forests to drought for this region (without human-driven changes to greenhouse gases, fire, and deforestation). We highlight the need for more palaeoecological and palaeoclimate data across lowland Amazonia to constrain model responses.

A Last Glacial Maximum forcing dataset for ocean modelling

Abstract. Model simulations of the Last Glacial Maximum (LGM; ∼ 21 000 years before present) can aid the interpretation of proxy records, can help to gain an improved mechanistic understanding of the LGM climate system, and are valuable for the evaluation of model performance in a different climate state. Ocean-ice only model configurations forced by prescribed atmospheric data (referred to as “forced ocean models”) drastically reduce the computational cost of palaeoclimate modelling compared to fully coupled model frameworks. While feedbacks between the atmosphere and ocean and sea-ice compartments of the Earth system are not present in such model configurations, many scientific questions can be addressed with models of this type. Our dataset supports simulations of the LGM in a forced ocean model set-up while still taking advantage of the complexity of fully coupled model set-ups. The data presented here are derived from fully coupled palaeoclimate simulations of the Palaeoclimate Modelling Intercomparison Project phase 3 (PMIP3). The data are publicly accessible at the National Infrastructure for Research Data (NIRD) Research Data Archive at https://doi.org/10.11582/2020.00052 (Morée and Schwinger, 2020). They consist of 2-D anomaly forcing fields suitable for use in ocean models that employ a bulk forcing approach and are optimized for use with CORE forcing fields. The data include specific humidity, downwelling long-wave and short-wave radiation, precipitation, wind (v and u components), temperature, and sea surface salinity (SSS). All fields are provided as climatological mean anomalies between LGM and pre-industrial (PI) simulations. These anomaly data can therefore be added to any pre-industrial ocean forcing dataset in order to obtain forcing fields representative of LGM conditions as simulated by PMIP3 models. Furthermore, the dataset can be easily updated to reflect results from upcoming and future palaeo-model intercomparison activities.

Development and testing scenarios for implementing land use and land cover changes during the Holocene in Earth system model experiments

Abstract. Anthropogenic changes in land use and land cover (LULC) during the pre-industrial Holocene could have affected regional and global climate. Existing scenarios of LULC changes during the Holocene are based on relatively simple assumptions and highly uncertain estimates of population changes through time. Archaeological and palaeoenvironmental reconstructions have the potential to refine these assumptions and estimates. The Past Global Changes (PAGES) LandCover6k initiative is working towards improved reconstructions of LULC globally. In this paper, we document the types of archaeological data that are being collated and how they will be used to improve LULC reconstructions. Given the large methodological uncertainties involved, both in reconstructing LULC from the archaeological data and in implementing these reconstructions into global scenarios of LULC, we propose a protocol to evaluate the revised scenarios using independent pollen-based reconstructions of land cover and climate. Further evaluation of the revised scenarios involves carbon cycle model simulations to determine whether the LULC reconstructions are consistent with constraints provided by ice core records of CO2 evolution and modern-day LULC. Finally, the protocol outlines how the improved LULC reconstructions will be used in palaeoclimate simulations in the Palaeoclimate Modelling Intercomparison Project to quantify the magnitude of anthropogenic impacts on climate through time and ultimately to improve the realism of Holocene climate simulations.

A Method for Generating Coherent Spatially Explicit Maps of Seasonal Paleoclimates From Site‐Based Reconstructions

AbstractWe describe a method for reconstructing spatially explicit maps of seasonal paleoclimate variables from site‐based reconstructions. Using a 3‐D‐Variational technique, the method finds the best statistically unbiased, and spatially continuous, estimate of the paleoclimate anomalies through combining the site‐based reconstructions and a prior estimate of the paleoclimate state. By assuming a set of correlations in the error of the prior, the resulting climate is smoothed both from month to month and from grid cell to grid cell. The amount of smoothing can be controlled through the choice of two length‐scale values. The method is applied to a set of reconstructions of the climate of the Last Glacial Maximum (ca. 21,000 years ago) for southern Europe derived from pollen data with a prior derived from results from the third phase of the Palaeoclimate Modelling Intercomparison Project. We demonstrate how to choose suitable values for the smoothing length scales for the data sets used in the reconstruction.

Ice-sheet configuration in the CMIP5/PMIP3 Last Glacial Maximum experiments

Abstract. We describe the creation of a data set describing changes related to the presence of ice sheets, including ice-sheet extent and height, ice-shelf extent, and the distribution and elevation of ice-free land at the Last Glacial Maximum (LGM), which were used in LGM experiments conducted as part of the fifth phase of the Coupled Modelling Intercomparison Project (CMIP5) and the third phase of the Palaeoclimate Modelling Intercomparison Project (PMIP3). The CMIP5/PMIP3 data sets were created from reconstructions made by three different groups, which were all obtained using a model-inversion approach but differ in the assumptions used in the modelling and in the type of data used as constraints. The ice-sheet extent in the Northern Hemisphere (NH) does not vary substantially between the three individual data sources. The difference in the topography of the NH ice sheets is also moderate, and smaller than the differences between these reconstructions (and the resultant composite reconstruction) and ice-sheet reconstructions used in previous generations of PMIP. Only two of the individual reconstructions provide information for Antarctica. The discrepancy between these two reconstructions is larger than the difference for the NH ice sheets, although still less than the difference between the composite reconstruction and previous PMIP ice-sheet reconstructions. Although largely confined to the ice-covered regions, differences between the climate response to the individual LGM reconstructions extend over the North Atlantic Ocean and Northern Hemisphere continents, partly through atmospheric stationary waves. Differences between the climate response to the CMIP5/PMIP3 composite and any individual ice-sheet reconstruction are smaller than those between the CMIP5/PMIP3 composite and the ice sheet used in the last phase of PMIP (PMIP2).

ENSO teleconnections with Australian rainfall in coupled model simulations of the last millennium

El Nino-Southern Oscillation is the major source of interannual rainfall variability in the Australian region, with the strongest influence over eastern Australia. The strength of this regional ENSO–rainfall teleconnection varies in the observational record. Climate model simulations of the “last millennium” (850–1850 C.E.) can be used to quantify the natural variability of the relationship between ENSO and Australian rainfall on decadal and longer time scales, providing a baseline for evaluating future projections. In this study, historical and last millennium (LM) simulations from six models were obtained from the Coupled Model Intercomparison Project Phase 5 and Palaeoclimate Modelling Intercomparison Project Phase 3. All models reproduce the observed negative correlation between September to February (SONDJF) eastern Australian rainfall and the NINO3.4 index, with varying skill. In the LM simulations, all models produce decadal-scale cooling over eastern Australia in response to volcanic forcing, as well as a long-term cooling trend. Rainfall variability over the same region is not strongly driven by external forcing, with each model simulating rainfall anomalies of different phase and magnitude. SONDJF eastern Australian rainfall is strongly correlated with ENSO in the LM simulations for all models, although some models simulate periods when the teleconnection weakens substantially for several decades. Changes in ENSO variance play a role in modulating the teleconnection strength, but are not the only factor. The long-term average spatial pattern of the ENSO–Australian rainfall teleconnection is similar in the LM and historical simulations, although the spatial pattern varies over time in the LM simulations.

Climate versus carbon dioxide controls on biomass burning: a model analysis of the glacial–interglacial contrast

Abstract. Climate controls fire regimes through its influence on the amount and types of fuel present and their dryness. CO2 concentration constrains primary production by limiting photosynthetic activity in plants. However, although fuel accumulation depends on biomass production, and hence on CO2 concentration, the quantitative relationship between atmospheric CO2 concentration and biomass burning is not well understood. Here a fire-enabled dynamic global vegetation model (the Land surface Processes and eXchanges model, LPX) is used to attribute glacial–interglacial changes in biomass burning to an increase in CO2, which would be expected to increase primary production and therefore fuel loads even in the absence of climate change, vs. climate change effects. Four general circulation models provided last glacial maximum (LGM) climate anomalies – that is, differences from the pre-industrial (PI) control climate – from the Palaeoclimate Modelling Intercomparison Project Phase~2, allowing the construction of four scenarios for LGM climate. Modelled carbon fluxes from biomass burning were corrected for the model's observed prediction biases in contemporary regional average values for biomes. With LGM climate and low CO2 (185 ppm) effects included, the modelled global flux at the LGM was in the range of 1.0–1.4 Pg C year-1, about a third less than that modelled for PI time. LGM climate with pre-industrial CO2 (280 ppm) yielded unrealistic results, with global biomass burning fluxes similar to or even greater than in the pre-industrial climate. It is inferred that a substantial part of the increase in biomass burning after the LGM must be attributed to the effect of increasing CO2 concentration on primary production and fuel load. Today, by analogy, both rising CO2 and global warming must be considered as risk factors for increasing biomass burning. Both effects need to be included in models to project future fire risks.

Recent and deep pasts in paleoclimate model intercomparison project

Second Workshop of the Palaeoclimate Modelling Intercomparison Project Phase 3;Crewe, United Kingdom, 6–11 May 2012 The Palaeoclimate Modelling Intercomparison Project (PMIP) is an international initiative endorsed by the World Climate Research Programme, the International Geosphere‐Biosphere Programme, and the Past Global Changes Project. Its objective is to provide a coherent framework to assess climate models and investigate climate dynamics on the basis of paleoclimate evidence. To this end, PMIP provides protocols for simulations of a set of past epochs with climate models, centralizes model output and facilitates access to it, and organizes syntheses of paleoclimate observations for model evaluation.

Evaluation of climate models using palaeoclimatic data

There is large uncertainty about the magnitude of warming and how rainfall patterns will change in response to any given scenario of future changes in atmospheric composition and land use. The models used for future climate projections were developed and calibrated using climate observations from the past 40 years. The geologic record of environmental responses to climate changes provides a unique opportunity to test model performance outside this limited climate range. Evaluation of model simulations against palaeodata shows that models reproduce the direction and large-scale patterns of past changes in climate, but tend to underestimate the magnitude of regional changes. As part of the effort to reduce model-related uncertainty and produce more reliable estimates of twenty-first century climate, the Palaeoclimate Modelling Intercomparison Project is systematically applying palaeoevaluation techniques to simulations of the past run with the models used to make future projections. This evaluation will provide assessments of model performance, including whether a model is sufficiently sensitive to changes in atmospheric composition, as well as providing estimates of the strength of biosphere and other feedbacks that could amplify the model response to these changes and modify the characteristics of climate variability.

Pliocene Model Intercomparison Project (PlioMIP): experimental design and boundary conditions (Experiment 2)

Abstract. The Palaeoclimate Modelling Intercomparison Project has expanded to include a model intercomparison for the mid-Pliocene warm period (3.29 to 2.97 million yr ago). This project is referred to as PlioMIP (the Pliocene Model Intercomparison Project). Two experiments have been agreed upon and together compose the initial phase of PlioMIP. The first (Experiment 1) is being performed with atmosphere-only climate models. The second (Experiment 2) utilises fully coupled ocean-atmosphere climate models. Following on from the publication of the experimental design and boundary conditions for Experiment 1 in Geoscientific Model Development, this paper provides the necessary description of differences and/or additions to the experimental design for Experiment 2.

Global vegetation and terrestrial carbon cycle changes after the last ice age

• In current models, the ecophysiological effects of CO₂ create both woody thickening and terrestrial carbon uptake, as observed now, and forest cover and terrestrial carbon storage increases that took place after the last glacial maximum (LGM). Here, we aimed to assess the realism of modelled vegetation and carbon storage changes between LGM and the pre-industrial Holocene (PIH). • We applied Land Processes and eXchanges (LPX), a dynamic global vegetation model (DGVM), with lowered CO₂ and LGM climate anomalies from the Palaeoclimate Modelling Intercomparison Project (PMIP II), and compared the model results with palaeodata. • Modelled global gross primary production was reduced by 27-36% and carbon storage by 550-694 Pg C compared with PIH. Comparable reductions have been estimated from stable isotopes. The modelled areal reduction of forests is broadly consistent with pollen records. Despite reduced productivity and biomass, tropical forests accounted for a greater proportion of modelled land carbon storage at LGM (28-32%) than at PIH (25%). • The agreement between palaeodata and model results for LGM is consistent with the hypothesis that the ecophysiological effects of CO₂ influence tree-grass competition and vegetation productivity, and suggests that these effects are also at work today.

Pliocene Model Intercomparison Project (PlioMIP): experimental design and boundary conditions (Experiment 1)

Abstract. In 2008 the temporal focus of the Palaeoclimate Modelling Intercomparison Project was expanded to include a model intercomparison for the mid-Pliocene warm period (3.29–2.97 million years ago). This project is referred to as PlioMIP (Pliocene Model Intercomparison Project). Two experiments have been agreed upon and comprise phase 1 of PlioMIP. The first (Experiment 1) will be performed with atmosphere-only climate models. The second (Experiment 2) will utilise fully coupled ocean-atmosphere climate models. The aim of this paper is to provide a detailed model intercomparison project description which documents the experimental design in a more detailed way than has previously been done in the literature. Specifically, this paper describes the experimental design and boundary conditions that will be utilised for Experiment 1 of PlioMIP.

The Southern Westerlies during the last glacial maximum in PMIP2 simulations

The Southern Hemisphere westerly winds are an important component of the climate system at hemispheric and global scales. Variations in their intensity and latitudinal position through an ice-age cycle have been proposed as important drivers of global climate change due to their influence on deep-ocean circulation and changes in atmospheric CO2. The position, intensity, and associated climatology of the southern westerlies during the last glacial maximum (LGM), however, is still poorly understood from empirical and modelling standpoints. Here we analyse the behaviour of the southern westerlies during the LGM using four coupled ocean-atmosphere simulations carried out by the Palaeoclimate Modelling Intercomparison Project Phase 2 (PMIP2). We analysed the atmospheric circulation by direct inspection of the winds and by using a cyclone tracking software to indicate storm tracks. The models suggest that changes were most significant during winter and over the Pacific ocean. For this season and region, three out four models indicate decreased wind intensities at the near surface as well as in the upper troposphere. Although the LGM atmosphere is colder and the equator to pole surface temperature gradient generally increases, the tropospheric temperature gradients actually decrease, explaining the weaker circulation. We evaluated the atmospheric influence on the Southern Ocean by examining the effect of wind stress on the Ekman pumping. Again, three of the models indicate decreased upwelling in a latitudinal band over the Southern Ocean. All models indicate a drier LGM than at present with a clear decrease in precipitation south of 40°S over the oceans. We identify important differences in precipitation anomalies over the land masses at regional scale, including a drier climate over New Zealand and wetter over NW Patagonia.

Evaluation of coupled ocean–atmosphere simulations of the mid-Holocene using palaeovegetation data from the northern hemisphere extratropics

We have used the BIOME4 biogeography–biochemistry model and comparison with palaeovegetation data to evaluate the response of six ocean–atmosphere general circulation models to mid-Holocene changes in orbital forcing in the mid- to high-latitudes of the northern hemisphere. All the models produce: (a) a northward shift of the northern limit of boreal forest, in response to simulated summer warming in high-latitudes. The northward shift is markedly asymmetric, with larger shifts in Eurasia than in North America; (b) an expansion of xerophytic vegetation in mid-continental North America and Eurasia, in response to increased temperatures during the growing season; (c) a northward expansion of temperate forests in eastern North America, in response to simulated winter warming. The northward shift of the northern limit of boreal forest and the northward expansion of temperate forests in North America are supported by palaeovegetation data. The expansion of xerophytic vegetation in mid-continental North America is consistent with palaeodata, although the extent may be over-estimated. The simulated expansion of xerophytic vegetation in Eurasia is not supported by the data. Analysis of an asynchronous coupling of one model to an equilibrium-vegetation model suggests vegetation feedback exacerbates this mid-continental drying and produces conditions more unlike the observations. Not all features of the simulations are robust: some models produce winter warming over Europe while others produce winter cooling. As a result, some models show a northward shift of temperate forests (consistent with, though less marked than, the expansion shown by data) and others produce a reduction in temperate forests. Elucidation of the cause of such differences is a focus of the current phase of the Palaeoclimate Modelling Intercomparison Project.

Predictive modelling of tree species distributions on the Iberian Peninsula during the Last Glacial Maximum and Mid‐Holocene

This paper reports a bioclimatic envelope model study of the potential distribution of 19 tree species in the Iberian Peninsula during the Last Glacial Maximum (LGM; 21 000 yr BP) and the Mid‐Holocene (6000 yr BP). Current patterns of tree species richness and distributions are believed to have been strongly influenced by the climate during these periods. The modelling employed novel machine learning techniques, and its accuracy was evaluated using a threshold‐independent method. Two atmospheric general circulation models, UGAMP and ECHAM3 (generated by the Palaeoclimate Modelling Intercomparison Project, PMIP), were used to provide climate scenarios under which the distributions of the 19 tree species were modelled. The results obtained for these scenarios were assessed by agreement measure analysis; they differed significantly for the LGM, but were more similar for the Mid‐Holocene.The results for the LGM support the inferred importance of pines in the Iberian Peninsula at this time, and the presence of evergreen Quercus in the south. Important differences in the altitude at which the modelled species grew were also predicted. During the LGM, some normally higher mountain species potentially became re‐established in the foothills of the Pyrenees. The warm Mid‐Holocene climate is clearly reflected in the predicted expansion of broad‐leaved forests during this period, including the colonization of the northern part of the Iberian Peninsula by evergreen Quercus species.

East European river runoff and Black Sea and Caspian Sea level changes as simulated within the Paleoclimate modeling intercomparison project

During the last late-Pleistocene and post-glacial epochs, the Caspian Sea fluctuated between regressive and transgressive stages. The Black Sea experienced fluctuations too, but their amplitude was smaller because they were mainly controlled by the World Ocean due to water exchange through the Bosporus Strait. In the water-budget equation, the change in storage over a time period is balanced by the sum of the inflows and outflows that occur during the time period. They can be calculated from atmospheric general circulation models (GCMs). A modeling initiative, the Palaeoclimate Modelling Intercomparison Project (PMIP), could help to resolve this question. The PMIP has focused on two slices of the past: the mid-Holocene (6ka calendar years Before Present (BP)) and the last cold event of the Late Quaternary (21ka calendar years BP). The PMIP GCM simulations indicate that changes in level are primarily influenced by changes of river runoff. Ocean precipitation and evaporation together play a minor role. At 21ka BP the estimated level lowering was ∼50 and ∼200m for the Caspian Sea and Black Sea, respectively. This lends credit to the idea of the connection between deep regressive states of Caspian Sea and Black Sea and mature stages of Late Quaternary glacial/cooling/drying planetary events.

Conservatism of ecological niche characteristics in North American plant species over the Pleistocene‐to‐Recent transition

AbstractAim To provide a test of the conservatism of a species’ niche over the last 20,000 years by tracking the distribution of eight pollen taxa relative to climate type as they migrated across eastern North America following the Last Glacial Maximum (LGM).Location North America.Methods We drew taxon occurrence data from the North American pollen records in the Global Pollen Database, representing eight pollen types – all taxa for which ≥5 distinct geographic occurrences were available in both the present day and at the LGM (21,000 years ago ± 3000 years). These data were incorporated into ecological niche models based on present‐day and LGM climatological summaries available from the Palaeoclimate Modelling Intercomparison Project to produce predicted potential geographic distributions for each species at present and at the LGM. The output for each time period was projected onto the ‘other’ time period, and tested using independent known occurrence information from that period.Results The result of our analyses was that all species tested showed general conservatism in ecological characteristics over the climate changes associated with the Pleistocene‐to‐Recent transition.Main conclusions This analysis constitutes a further demonstration of general and pervasive conservatism in ecological niche characteristics over moderate periods of time despite profound changes in climate and environmental conditions. As such, our results reinforce the application of ecological niche modelling techniques to the reconstruction of Pleistocene biodiversity distribution patterns, and to project the future potential distribution range of species in the face of global‐scale climatic changes.

Pollen‐based reconstructions of biome distributions for Australia, Southeast Asia and the Pacific (SEAPAC region) at 0, 6000 and 18,00014C yr BP

AbstractAim This paper documents reconstructions of the vegetation patterns in Australia, Southeast Asia and the Pacific (SEAPAC region) in the mid‐Holocene and at the last glacial maximum (LGM).Methods Vegetation patterns were reconstructed from pollen data using an objective biomization scheme based on plant functional types. The biomization scheme was first tested using 535 modern pollen samples from 377 sites, and then applied unchanged to fossil pollen samples dating to 6000 ± 500 or 18,000 ± 100014C yrbp.Results 1. Tests using surface pollen sample sites showed that the biomization scheme is capable of reproducing the modern broad‐scale patterns of vegetation distribution. The north–south gradient in temperature, reflected in transitions from cool evergreen needleleaf forest in the extreme south through temperate rain forest or wet sclerophyll forest (WSFW) and into tropical forests, is well reconstructed. The transitions from xerophytic through sclerophyll woodlands and open forests to closed‐canopy forests, which reflect the gradient in plant available moisture from the continental interior towards the coast, are reconstructed with less geographical precision but nevertheless the broad‐scale pattern emerges. 2. Differences between the modern and mid‐Holocene vegetation patterns in mainland Australia are comparatively small and reflect changes in moisture availability rather than temperature. In south‐eastern Australia some sites show a shift towards more moisture‐stressed vegetation in the mid‐Holocene with xerophytic woods/scrub and temperate sclerophyll woodland and shrubland at sites characterized today by WSFW or warm‐temperate rain forest (WTRF). However, sites in the Snowy Mountains, on the Southern Tablelands and east of the Great Dividing Range have more moisture‐demanding vegetation in the mid‐Holocene than today. South‐western Australia was slightly drier than today. The single site in north‐western Australia also shows conditions drier than today in the mid‐Holocene. Changes in the tropics are also comparatively small, but the presence of WTRF and tropical deciduous broadleaf forest and woodland in the mid‐Holocene, in sites occupied today by cool‐temperate rain forest, indicate warmer conditions. 3. Expansion of xerophytic vegetation in the south and tropical deciduous broadleaf forest and woodland in the north indicate drier conditions across mainland Australia at the LGM. None of these changes are informative about the degree of cooling. However the evidence from the tropics, showing lowering of the treeline and forest belts, indicates that conditions were between 1 and 9 °C (depending on elevation) colder. The encroachment of tropical deciduous broadleaf forest and woodland into lowland evergreen broadleaf forest implies greater aridity.Main conclusions This study provides the first continental‐scale reconstruction of mid‐Holocene and LGM vegetation patterns from Australia, Southeast Asia and the Pacific (SEAPAC region) using an objective biomization scheme. These data will provide a benchmark for evaluation of palaeoclimate simulations within the framework of the Palaeoclimate Modelling Intercomparison Project.