Last Glacial Maximum Simulations Research Articles

Last glacial maximum permafrost in China from CMIP5 simulations

The distributions of frozen ground and active layer thickness (ALT) during the last glacial maximum (LGM) and pre-industrial periods in China were investigated using Coupled Model Intercomparison Project Phase 5 (CMIP5) model experiments. Compared to the pre-industrial period, the LGM climate was ~5°C colder and featured significantly higher freezing indices on the Tibetan Plateau and in Northeast China. Frozen ground expanded widely in the LGM. The extents of permafrost and seasonally frozen ground in China were 4.11×106km2 and 4.97×106km2, respectively, which are 2.42×106km2 larger and 1.45×106km2 smaller, respectively, than the pre-industrial levels. Moreover, the colder climate and longer duration also resulted in LGM ALT values that were 1.3m less than the pre-industrial values in the permafrost areas common to both periods. Altitudinal permafrost was present mainly on the Tibetan Plateau and adjacent mountains in West China between 28°N and 41°30′N and covered an area of ~2.63×106km2. Latitudinal permafrost was present mainly in Northeast China and occupied an area of 1.48×106km2. The southern limit of latitudinal permafrost was located ~10° of latitude farther south during the LGM than during the pre-industrial period. The LGM simulation results agree reasonably well with previous reconstructions, with the exception of an underestimation in the permafrost extent. Although relatively high-level disagreement exists between the models in terms of the exact locations of the southern limits, the ensemble average is still able to represent the large-scale spatial pattern of frozen ground remarkably well.

How well do simulated last glacial maximum tropical temperatures constrain equilibrium climate sensitivity?

AbstractPrevious work demonstrated a significant correlation between tropical surface air temperature and equilibrium climate sensitivity (ECS) in PMIP (Paleoclimate Modelling Intercomparison Project) phase 2 model simulations of the last glacial maximum (LGM). This implies that reconstructed LGM cooling in this region could provide information about the climate system ECS value. We analyze results from new simulations of the LGM performed as part of Coupled Model Intercomparison Project (CMIP5) and PMIP phase 3. These results show no consistent relationship between the LGM tropical cooling and ECS. A radiative forcing and feedback analysis shows that a number of factors are responsible for this decoupling, some of which are related to vegetation and aerosol feedbacks. While several of the processes identified are LGM specific and do not impact on elevated CO2 simulations, this analysis demonstrates one area where the newer CMIP5 models behave in a qualitatively different manner compared with the older ensemble. The results imply that so‐called Earth System components such as vegetation and aerosols can have a significant impact on the climate response in LGM simulations, and this should be taken into account in future analyses.

Range and niche shifts in response to past climate change in the desert horned lizard (Phrynosoma platyrhinos).

During climate change, species are often assumed to shift their geographic distributions (geographic ranges) in order to track environmental conditions - niches - to which they are adapted. Recent work, however, suggests that the niches do not always remain conserved during climate change but shift instead, allowing populations to persist in place or expand into new areas. We assessed the extent of range and niche shifts in response to the warming climate after the Last Glacial Maximum (LGM) in the desert horned lizard (Phrynosoma platyrhinos), a species occupying the western deserts of North America. We used a phylogeographic approach with mitochondrial DNA sequences to approximate the species range during the LGM by identifying populations that exhibit a genetic signal of population stability versus those that exhibit a signal of a recent (likely post-LGM) geographic expansion. We then compared the climatic niche that the species occupies today with the niche it occupied during the LGM using two models of simulated LGM climate. The genetic analyses indicated that P. platyrhinos persisted within the southern Mojave and Sonoran deserts throughout the latest glacial period and expanded from these deserts northwards, into the western and eastern Great Basin, after the LGM. The climatic niche comparisons revealed that P. platyrhinos expanded its climatic niche after the LGM towards novel, warmer and drier climates that allowed it to persist within the southern deserts. Simultaneously, the species shifted its climatic niche towards greater temperature and precipitation fluctuations after the LGM. We concluded that climatic changes at the end of the LGM promoted both range and niche shifts in this lizard. The mechanism that allowed the species to shift its niche remains unknown, but phenotypic plasticity likely contributes to the species ability to adjust to climate change.

Coupled ice sheet–climate modeling under glacial and pre-industrial boundary conditions

Abstract. In the standard Paleoclimate Modelling Intercomparison Project (PMIP) experiments, the Last Glacial Maximum (LGM) is modeled in quasi-equilibrium with atmosphere–ocean–vegetation general circulation models (AOVGCMs) with prescribed ice sheets. This can lead to inconsistencies between the modeled climate and ice sheets. One way to avoid this problem would be to model the ice sheets explicitly. Here, we present the first results from coupled ice sheet–climate simulations for the pre-industrial times and the LGM. Our setup consists of the AOVGCM ECHAM5/MPIOM/LPJ bidirectionally coupled with the Parallel Ice Sheet Model (PISM) covering the Northern Hemisphere. The results of the pre-industrial and LGM simulations agree reasonably well with reconstructions and observations. This shows that the model system adequately represents large, non-linear climate perturbations. A large part of the drainage of the ice sheets occurs in ice streams. Most modeled ice stream systems show recurring surges as internal oscillations. The Hudson Strait Ice Stream surges with an ice volume equivalent to about 5 m sea level and a recurrence interval of about 7000 yr. This is in agreement with basic expectations for Heinrich events. Under LGM boundary conditions, different ice sheet configurations imply different locations of deep water formation.

Assessing the impact of Laurentide Ice Sheet topography on glacial climate

Abstract. Simulations of past climates require altered boundary conditions to account for known shifts in the Earth system. For the Last Glacial Maximum (LGM) and subsequent deglaciation, the existence of large Northern Hemisphere ice sheets caused profound changes in surface topography and albedo. While ice-sheet extent is fairly well known, numerous conflicting reconstructions of ice-sheet topography suggest that precision in this boundary condition is lacking. Here we use a high-resolution and oxygen-isotope-enabled fully coupled global circulation model (GCM) (GISS ModelE2-R), along with two different reconstructions of the Laurentide Ice Sheet (LIS) that provide maximum and minimum estimates of LIS elevation, to assess the range of climate variability in response to uncertainty in this boundary condition. We present this comparison at two equilibrium time slices: the LGM, when differences in ice-sheet topography are maximized, and 14 ka, when differences in maximum ice-sheet height are smaller but still exist. Overall, we find significant differences in the climate response to LIS topography, with the larger LIS resulting in enhanced Atlantic Meridional Overturning Circulation and warmer surface air temperatures, particularly over northeastern Asia and the North Pacific. These up- and downstream effects are associated with differences in the development of planetary waves in the upper atmosphere, with the larger LIS resulting in a weaker trough over northeastern Asia that leads to the warmer temperatures and decreased albedo from snow and sea-ice cover. Differences between the 14 ka simulations are similar in spatial extent but smaller in magnitude, suggesting that climate is responding primarily to the larger difference in maximum LIS elevation in the LGM simulations. These results suggest that such uncertainty in ice-sheet boundary conditions alone may significantly impact the results of paleoclimate simulations and their ability to successfully simulate past climates, with implications for estimating climate sensitivity to greenhouse gas forcing utilizing past climate states.

New Estimates of Permafrost Evolution during the Last 21 k Years in Eurasia using Numerical Modelling

ABSTRACTThe evolution of past permafrost since the Last Glacial Maximum (LGM) is simulated using the Vrije Universiteit Amsterdam Permafrost (VAMPER) model. This method is different from a proxy‐based approach which translates reconstructed air temperatures to estimate past permafrost extent and thickness. First, a sensitivity analysis was performed to assess the behaviour of the model. Then five case studies within Eurasia were performed using mean annual ground surface temperatures derived from an Earth system model as the surface forcing. In Central and West Siberia, the simulated LGM permafrost thicknesses of 730–940 m and 365–445 m, respectively, agree well with previous estimates. The LGM and present‐day estimates for South Russia (9–15 m) are underestimated, which is likely due to a highly simplified land‐atmosphere coupling. In West and Central Europe, however, the VAMPER model was not able to produce permafrost during LGM conditions, which is due to previously recognised biases of the Earth system model. A supplementary simulation was then performed, resulting in an LGM permafrost thickness estimate of 260–320 m. Average thawing rates are on the order of 1 to 3 cm/yr except for Central Siberia, where permafrost thawed at rates of 0.3 to 0.4 cm/yr. Overall results of these simulations provide a basis for future improvement in modelling the permafrost‐climate relationship over millennia. Copyright © 2013 John Wiley & Sons, Ltd.

Asian monsoon climate during the Last Glacial Maximum: palaeo‐data–model comparisons

The Last Glacial Maximum (LGM) (23–19 ka BP) in the Asian monsoon region is generally described as cool and dry, due to a strong winter monsoon. More recently, however, palaeo‐data and climate model simulations have argued for a more variable LGM Asian monsoon climate with distinct regional differences. We compiled, evaluated, and partly re‐assessed proxy records for the Asian monsoon region in terms of wet/dry climatic conditions based on precipitation and effective moisture, and of sea surface temperatures. The comparison of the palaeo‐data set to LGM simulations by the Climate Community System Model version 3 (CCSM3) shows fairly good agreement: a dry LGM climate in the western and northern part due to a strengthened winter monsoon and/or strengthened westerly winds and wetter conditions in equatorial areas, due to a stronger summer monsoon. Data–model discrepancies are seen in some areas and are ascribed to the fairly coarse resolution of CCSM3 and/or to uncertainties in the reconstructions. Differences are also observed between the reconstructed and simulated northern boundaries of the Intertropical Convergence Zone (ITCZ). The reconstructions estimate a more southern position over southern India and the Bay of Bengal, whereas CCSM3 simulates a more northern position. In Indochina, the opposite is the case. The palaeo‐data indicate that climatic conditions changed around 20–19 ka BP, with some regions receiving higher precipitation and some experiencing drier conditions, which would imply a distinct shift in summer monsoon intensity. This shift was probably triggered by the late LGM sea‐level rise, which led to changes in atmosphere–ocean interactions in the Indian Ocean. The overall good correspondence between reconstructions and CCSM3 suggests that CCSM3 simulates LGM climate conditions over subtropical and tropical areas fairly well. The few high‐resolution qualitative and quantitative palaeo‐records available for the large Asian monsoon region make reconstructions however still uncertain.

Influence of Last Glacial Maximum boundary conditions on the global water isotope distribution in an atmospheric general circulation model

Abstract. To understand the validity of δ18O proxy records as indicators of past temperature change, a series of experiments was conducted using an atmospheric general circulation model fitted with water isotope tracers (Community Atmosphere Model version 3.0, IsoCAM). A pre-industrial simulation was performed as the control experiment, as well as a simulation with all the boundary conditions set to Last Glacial Maximum (LGM) values. Results from the pre-industrial and LGM simulations were compared to experiments in which the influence of individual boundary conditions (greenhouse gases, ice sheet albedo and topography, sea surface temperature (SST), and orbital parameters) were changed each at a time to assess their individual impact. The experiments were designed in order to analyze the spatial variations of the oxygen isotopic composition of precipitation (δ18Oprecip) in response to individual climate factors. The change in topography (due to the change in land ice cover) played a significant role in reducing the surface temperature and δ18Oprecip over North America. Exposed shelf areas and the ice sheet albedo reduced the Northern Hemisphere surface temperature and δ18Oprecip further. A global mean cooling of 4.1 °C was simulated with combined LGM boundary conditions compared to the control simulation, which was in agreement with previous experiments using the fully coupled Community Climate System Model (CCSM3). Large reductions in δ18Oprecip over the LGM ice sheets were strongly linked to the temperature decrease over them. The SST and ice sheet topography changes were responsible for most of the changes in the climate and hence the δ18Oprecip distribution among the simulations.

Sensitivity to Glacial Forcing in the CCSM4

AbstractResults are presented from the Community Climate System Model, version 4 (CCSM4), simulation of the Last Glacial Maximum (LGM) from phase 5 of the Coupled Model Intercomparison Project (CMIP5) at the standard 1° resolution, the same resolution as the majority of the CCSM4 CMIP5 long-term simulations for the historical and future projection scenarios. The forcings and boundary conditions for this simulation follow the protocols of the Paleoclimate Modeling Intercomparison Project, version 3 (PMIP3). Two additional CCSM4 CO2 sensitivity simulations, in which the concentrations are abruptly changed at the start of the simulation to the low 185 ppm LGM concentrations (LGMCO2) and to a quadrupling of the preindustrial concentration (4×CO2), are also analyzed. For the full LGM simulation, the estimated equilibrium cooling of the global mean annual surface temperature is 5.5°C with an estimated radiative forcing of −6.2 W m−2. The radiative forcing includes the effects of the reduced LGM greenhouse gases, ice sheets, continental distribution with sea level lowered by approximately 120 m from the present, and orbital parameters, but not changes to atmospheric aerosols or vegetation biogeography. The LGM simulation has an equilibrium climate sensitivity (ECS) of 3.1(±0.3)°C, comparable to the CCSM4 4×CO2 result. The LGMCO2 simulation shows a greater ECS of 4.2°C. Other responses found at the LGM in CCSM4 include a global precipitation rate decrease at a rate of ~2% °C−1, similar to climate change simulations in the Intergovernmental Panel on Climate Change Fourth Assessment Report (IPCC AR4); a strengthening of the Atlantic meridional overturning circulation (AMOC) with a shoaling of North Atlantic Deep Water and a filling of the deep basin up to sill depth with Antarctic Bottom Water; and an enhanced seasonal cycle accompanied by reduced ENSO variability in the eastern Pacific Ocean’s SSTs.

Sensitivity of Southern Hemisphere circulation to LGM and 4 × CO2 climates

This paper investigates the effect of Last Glacial Maximum (LGM) versus high CO2 world boundary condition on the Southern Hemisphere atmospheric circulation, in particular on the strength and latitudinal position of the near surface Southern Westerly Winds (SWW). PMIP2 and PMIP3 experiments, as well as the “abrupt 4 × CO2” simulations from CMIP5, were analyzed. Robust findings include poleward expansion of the Mean Meridional Circulation (MMC) and intensified and poleward‐shifted SWW in the 4 × CO2 simulations (consistent with recent observations and 21st century climate change projections); and for the LGM, stronger and southward shifted northern hemisphere MMC, and weakened southern Hadley cell. However, six of the eight LGM simulations show a decrease in the SWW, the other two models simulate the opposite. A critical difference between the models is strong coupling between sea‐ice extent, surface temperature gradients, SWW, and Ferrel cell in the two models with stronger and poleward‐shifted SWW.

Southern Hemisphere westerly wind changes during the Last Glacial Maximum: model-data comparison

The Southern Hemisphere (SH) westerly winds are thought to be critical to global ocean circulation, productivity, and carbon storage. For example, an equatorward shift in the winds, though its affect on the Southern Ocean circulation, has been suggested as the leading cause for the reduction in atmospheric CO2 during the Last Glacial period. Despite the importance of the winds, it is currently not clear, from observations or model results, how they behave during the Last Glacial. Here, an atmospheric modelling study is performed to help determine likely changes in the SH westerly winds during the Last Glacial Maximum (LGM). Using LGM boundary conditions, the maximum in SH westerlies is strengthened by ∼+1 m s−1 and moved southward by ∼2° at the 850 hPa pressure level. Boundary layer stabilisation effects over equatorward extended LGM sea-ice can lead to a small apparent equatorward shift in the wind band at the surface. Further sensitivity analysis with individual boundary condition changes indicate that changes in sea surface temperatures are the strongest factor behind the wind change. The HadAM3 atmospheric simulations, along with published PMIP2 coupled climate model simulations, are then assessed against the newly synthesised database of moisture observations for the LGM. Although the moisture data is the most commonly cited evidence in support of a large equatorward shift in the SH winds during the LGM, none of the models that produce realistic LGM precipitation changes show such a large equatorward shift. In fact, the model which best simulates the moisture proxy data is the HadAM3 LGM simulation which shows a small poleward wind shift. While we cannot prove here that a large equatorward shift would not be able to reproduce the moisture data as well, we show that the moisture proxies do not provide an observational evidence base for it.

Can the Last Glacial Maximum constrain climate sensitivity?

We investigate the relationship between the Last Glacial Maximum (LGM) and climate sensitivity across the PMIP2 multi‐model ensemble of GCMs, and find a correlation between tropical temperature and climate sensitivity which is statistically significant and physically plausible. We use this relationship, together with the LGM temperature reconstruction of Annan and Hargreaves (2012), to generate estimates for the equilibrium climate sensitivity. We estimate the equilibrium climate sensitivity to be about 2.5°C with a high probability of being under 4°C, though these results are subject to several important caveats. The forthcoming PMIP3/CMIP5 models were not considered in this analysis, as very few LGM simulations are currently available from these models. We propose that these models will provide a useful validation of the correlation presented here.

Mid-Holocene and Last Glacial Maximum climate simulations with the IPSL model—part I: comparing IPSL_CM5A to IPSL_CM4

The climates of the mid-Holocene (MH), 6,000 years ago, and of the Last Glacial Maximum (LGM), 21,000 years ago, have extensively been simulated, in particular in the framework of the Palaeoclimate Modelling Intercomparion Project. These periods are well documented by paleo-records, which can be used for evaluating model results for climates different from the present one. Here, we present new simulations of the MH and the LGM climates obtained with the IPSL_CM5A model and compare them to our previous results obtained with the IPSL_CM4 model. Compared to IPSL_CM4, IPSL_CM5A includes two new features: the interactive representation of the plant phenology and marine biogeochemistry. But one of the most important differences between these models is the latitudinal resolution and vertical domain of their atmospheric component, which have been improved in IPSL_CM5A and results in a better representation of the mid-latitude jet-streams. The Asian monsoon’s representation is also substantially improved. The global average mean annual temperature simulated for the pre-industrial (PI) period is colder in IPSL_CM5A than in IPSL_CM4 but their climate sensitivity to a CO2 doubling is similar. Here we show that these differences in the simulated PI climate have an impact on the simulated MH and LGM climatic anomalies. The larger cooling response to LGM boundary conditions in IPSL_CM5A appears to be mainly due to differences between the PMIP3 and PMIP2 boundary conditions, as shown by a short wave radiative forcing/feedback analysis based on a simplified perturbation method. It is found that the sensitivity computed from the LGM climate is lower than that computed from 2 × CO2 simulations, confirming previous studies based on different models. For the MH, the Asian monsoon, stronger in the IPSL_CM5A PI simulation, is also more sensitive to the insolation changes. The African monsoon is also further amplified in IPSL_CM5A due to the impact of the interactive phenology. Finally the changes in variability for both models and for MH and LGM are presented taking the example of the El-Nino Southern Oscillation (ENSO), which is very different in the PI simulations. ENSO variability is damped in both model versions at the MH, whereas inconsistent responses are found between the two versions for the LGM. Part 2 of this paper examines whether these differences between IPSL_CM4 and IPSL_CM5A can be distinguished when comparing those results to palaeo-climatic reconstructions and investigates new approaches for model-data comparisons made possible by the inclusion of new components in IPSL_CM5A.

Changes in the intermediate water mass formation rates in the global ocean for the Last Glacial Maximum, mid‐Holocene and pre‐industrial climates

The paleoclimate version of the National Center for Atmospheric Research Community Climate System Model version 3 (NCAR‐CCSM3) is used to analyze changes in the water formation rates in the Atlantic, Pacific, and Indian Oceans for the Last Glacial Maximum (LGM), mid‐Holocene (MH) and pre‐industrial (PI) control climate. During the MH, CCSM3 exhibits a north‐south asymmetric response of intermediate water subduction changes in the Atlantic Ocean, with a reduction of 2 Sv in the North Atlantic and an increase of 2 Sv in the South Atlantic relative to PI. During the LGM, there is increased formation of intermediate water and a more stagnant deep ocean in the North Pacific. The production of North Atlantic Deep Water (NADW) is significantly weakened. The NADW is replaced in large extent by enhanced Antarctic Intermediate Water (AAIW), Glacial North Atlantic Intermediate Water (GNAIW), and also by an intensified of Antarctic Bottom Water (AABW), with the latter being a response to the enhanced salinity and ice formation around Antarctica. Most of the LGM intermediate/mode water is formed at 27.4 <σθ < 29.0 kg/m3, while for the MH and PI most of the subduction transport occurs at 26.5 < σθ < 27.4 kg/m3. The simulated LGM Southern Hemisphere winds are more intense by 0.2–0.4 dyne/cm2. Consequently, increased Ekman transport drives the production of intermediate water (low salinity) at a larger rate and at higher densities when compared to the other climatic periods.

Tropical Cyclone Genesis Factors in Simulations of the Last Glacial Maximum

Abstract Large-scale environmental factors that favor tropical cyclogenesis are calculated and examined in simulations of the Last Glacial Maximum (LGM) from the Paleoclimate Modelling Intercomparison Project Phase 2 (PMIP2). Despite universally colder conditions at the LGM, values of tropical cyclone potential intensity, which both serves as an upper bound on thermodynamically achievable intensity and indicates regions supportive of the deep convection required, are broadly similar in magnitude to those in preindustrial era control simulation. Some regions, including large areas of the central and western North Pacific, feature higher potential intensities at the LGM than they do in the control runs, while other regions including much of the Atlantic and Indian Oceans are lower. Changes in potential intensity are strongly correlated with the degree of surface cooling during the LGM. Additionally, two thermodynamic parameters—one that measures midtropospheric entropy deficits relevant for tropical cyclogenesis and another related to the time required for genesis—are broadly more favorable in the LGM simulation than in the preindustrial era control. A genesis potential index yields higher values for the LGM in much of the western Pacific, a feature common to nearly all of the individual models examined.

The thermal threshold of the Atlantic meridional overturning circulation and its control by wind stress forcing during glacial climate

Paleo proxy data suggest that the Atlantic meridional overturning circulation (AMOC) was shallower and weaker at the Last Glacial Maximum (LGM) than at present. In this study, we have identified the existence of a thermal threshold of the AMOC which may explain why many coupled climate models fail to simulate the weaker AMOC during the LGM. By using results obtained from a coupled climate model and conducting sensitivity simulations with an ocean general circulation model, we found that the sudden transition from the present‐day AMOC to the weaker glacial AMOC occurs when we gradually change the degree of surface cooling from present‐day to glacial conditions. This result is related to response of deep convection in the northern North Atlantic Ocean; moderate cooling enhances deep convection whereas sufficient cooling results in total covering of sea ice there and suppression of deep convection. The findings indicate the existence of a thermal threshold controlling the AMOC, where the present‐day‐type AMOC suddenly shifts to the weaker glacial AMOC once the surface cooling exceeds this threshold. We also demonstrate that wind stress forcing plays a critical role in controlling the value of the thermal threshold. Our study suggests that slight differences in the degree of surface cooling or wind stress forcing for LGM simulations could lead to the very different response of the AMOC during the LGM as reported in previous LGM simulations.

Tropical climate and vegetation changes during Heinrich Event 1: a model-data comparison

Abstract. Abrupt climate changes from 18 to 15 thousand years before present (kyr BP) associated with Heinrich Event 1 (HE1) had a strong impact on vegetation patterns not only at high latitudes of the Northern Hemisphere, but also in the tropical regions around the Atlantic Ocean. To gain a better understanding of the linkage between high and low latitudes, we used the University of Victoria (UVic) Earth System-Climate Model (ESCM) with dynamical vegetation and land surface components to simulate four scenarios of climate-vegetation interaction: the pre-industrial era, the Last Glacial Maximum (LGM), and a Heinrich-like event with two different climate backgrounds (interglacial and glacial). We calculated mega-biomes from the plant-functional types (PFTs) generated by the model to allow for a direct comparison between model results and palynological vegetation reconstructions. Our calculated mega-biomes for the pre-industrial period and the LGM corresponded well with biome reconstructions of the modern and LGM time slices, respectively, except that our pre-industrial simulation predicted the dominance of grassland in southern Europe and our LGM simulation resulted in more forest cover in tropical and sub-tropical South America. The HE1-like simulation with a glacial climate background produced sea-surface temperature patterns and enhanced inter-hemispheric thermal gradients in accordance with the "bipolar seesaw" hypothesis. We found that the cooling of the Northern Hemisphere caused a southward shift of those PFTs that are indicative of an increased desertification and a retreat of broadleaf forests in West Africa and northern South America. The mega-biomes from our HE1 simulation agreed well with paleovegetation data from tropical Africa and northern South America. Thus, according to our model-data comparison, the reconstructed vegetation changes for the tropical regions around the Atlantic Ocean were physically consistent with the remote effects of a Heinrich event under a glacial climate background.

Present and LGM permafrost from climate simulations: contribution of statistical downscaling

Abstract. We quantify the agreement between permafrost distributions from PMIP2 (Paleoclimate Modeling Intercomparison Project) climate models and permafrost data. We evaluate the ability of several climate models to represent permafrost and assess the variability between their results. Studying a heterogeneous variable such as permafrost implies conducting analysis at a smaller spatial scale compared with climate models resolution. Our approach consists of applying statistical downscaling methods (SDMs) on large- or regional-scale atmospheric variables provided by climate models, leading to local-scale permafrost modelling. Among the SDMs, we first choose a transfer function approach based on Generalized Additive Models (GAMs) to produce high-resolution climatology of air temperature at the surface. Then we define permafrost distribution over Eurasia by air temperature conditions. In a first validation step on present climate (CTRL period), this method shows some limitations with non-systematic improvements in comparison with the large-scale fields. So, we develop an alternative method of statistical downscaling based on a Multinomial Logistic GAM (ML-GAM), which directly predicts the occurrence probabilities of local-scale permafrost. The obtained permafrost distributions appear in a better agreement with CTRL data. In average for the nine PMIP2 models, we measure a global agreement with CTRL permafrost data that is better when using ML-GAM than when applying the GAM method with air temperature conditions. In both cases, the provided local information reduces the variability between climate models results. This also confirms that a simple relationship between permafrost and the air temperature only is not always sufficient to represent local-scale permafrost. Finally, we apply each method on a very different climate, the Last Glacial Maximum (LGM) time period, in order to quantify the ability of climate models to represent LGM permafrost. The prediction of the SDMs (GAM and ML-GAM) is not significantly in better agreement with LGM permafrost data than large-scale fields. At the LGM, both methods do not reduce the variability between climate models results. We show that LGM permafrost distribution from climate models strongly depends on large-scale air temperature at the surface. LGM simulations from climate models lead to larger differences with LGM data than in the CTRL period. These differences reduce the contribution of downscaling.

Contrasting wetland CH4 emission responses to simulated glacial atmospheric CO2 in temperate bogs and fens

Wetlands were the largest source of atmospheric methane (CH(4) ) during the Last Glacial Maximum (LGM), but the sensitivity of this source to exceptionally low atmospheric CO(2) concentration ([CO(2) ]) at the time has not been examined experimentally. We tested the hypothesis that LGM atmospheric [CO(2) ] reduced CH(4) emissions as a consequence of decreased photosynthate allocation to the rhizosphere. We exposed minerotrophic fen and ombrotrophic bog peatland mesocosms to simulated LGM (c. 200 ppm) or ambient (c. 400 ppm) [CO(2) ] over 21 months (n = 8 per treatment) and measured gaseous CH(4) flux, pore water dissolved CH(4) and volatile fatty acid (VFA; an indicator of plant carbon supply to the rhizosphere) concentrations. Cumulative CH(4) flux from fen mesocosms was suppressed by 29% (P < 0.05) and rhizosphere pore water [CH(4) ] by c. 50% (P < 0.01) in the LGM [CO(2) ], variables that remained unaffected in bog mesocosms. VFA analysis indicated that changes in plant root exudates were not the driving mechanism behind these results. Our data suggest that the LGM [CO(2) ] suppression of wetland CH(4) emissions is contingent on trophic status. The heterogeneous response may be attributable to differences in species assemblage that influence the dominant CH(4) production pathway, rhizosphere supplemented photosynthesis and CH(4) oxidation.

Links between Rossby Wave Breaking and the North Atlantic Oscillation–Arctic Oscillation in Present-Day and Last Glacial Maximum Climate Simulations

Abstract Upper-tropospheric Rossby wave–breaking processes are examined in coupled ocean–atmosphere simulations of the Last Glacial Maximum (LGM) and of the modern era. LGM statistics of the Northern Hemisphere in winter, computed from the Paleoclimate Modeling Intercomparison Project Phase II (PMIP2) dataset, are compared with those from preindustrial simulations and from the 40-yr ECMWF Re-Analysis (ERA-40). Particular attention is given to the role of wave-breaking events in the North Atlantic Oscillation (NAO) for each simulation. Anticyclonic (AWB) and cyclonic (CWB) wave-breaking events during LGM are shown to be less and more frequent, respectively, than in the preindustrial climate, especially in the Pacific. This is consistent with the slight equatorward shift of the eddy-driven jets in the LGM runs. The most remarkable feature of the simulated LGM climate is that it presents much weaker latitudinal fluctuations of the eddy-driven jets. This is accompanied by less dispersion in the wave-breaking events. A physical interpretation is provided in terms of the fluctuations of the low-level baroclinicity at the entrance of the storm tracks. The NAO in the preindustrial simulations and in ERA-40 is characterized by strong latitudinal fluctuations of the Atlantic eddy-driven jet as well as by significant changes in the nature of the wave breaking. During the positive phase, the eddy-driven jet moves to the north with more AWB events than usual and is well separated from the subtropical African jet. The negative phase exhibits a more equatorward Atlantic jet and more CWB events. In contrast, the LGM NAO is less well marked by the latitudinal vacillation of the Atlantic jet and for some models this property disappears entirely. The LGM NAO corresponds more to acceleration–deceleration or extension–retraction of the Atlantic jet. The hemispheric point of view of the Arctic Oscillation exhibits similar changes.