Publications

Abstract: Climate change over High Mountain Asia (HMA, including the Tibetan Plateau) is investigated over the period 1979-2014 and in future projections following the four Shared Socioeconomic Pathways: SSP1-2.6, SSP2-4.5, SSP3-7.0 and SSP5-8.5. The skill of 26 Coupled Model Intercomparison Project phase 6 (CMIP6) models is estimated for near-surface air temperature, snow cover extent and total precipitation, and 10 of them are used to describe their projections until 2100. Similarly to previous CMIP models, this new generation of general circulation models (GCMs) shows a mean cold bias over this area reaching -1.9 [ -8.2 to 2.9] degrees C (90 % confidence interval) in comparison with the Climate Research Unit (CRU) observational dataset, associated with a snow cover mean overestimation of 12 % [ -13 % to 43 %], corresponding to a relative bias of 52 % [ -53 % to 183 %] in comparison with the NOAA Climate Data Record (CDR) satellite dataset. The temperature and snow cover model biases are more pronounced in winter. Simulated precipitation rates are overestimated by 1.5 [0.3 to 2 9] mm d(-1), corresponding to a relative bias of 143 % [31 % to 281 %], but this might be an apparent bias caused by the undercatch of solid precipitation in the APHRODI1E (Asian Precipitation-Highly-Resolved Observational Data Integration Towards Evaluation of Water Resources) observational reference. For most models, the cold surface bias is associated with an overestimation of snow cover extent, but this relationship does not hold for all models, suggesting that the processes of the origin of the biases can differ from one model to another. A significant correlation between snow cover bias and surface elevation is found, and to a lesser extent between temperature bias and surface elevation, highlighting the model weaknesses at high elevation. The models with the best performance for temperature are not necessarily the most skillful for the other variables, and there is no clear relationship between model resolution and model skill. This highlights the need for a better understanding of the physical processes driving the climate in this complex topographic area, as well as for further parameterization developments adapted to such areas. A dependency of the simulated past trends on the model biases is found for some variables and seasons; however, some highly biased models fall within the range of observed trends, suggesting that model bias is not a robust criterion to discard models in trend analysis. The HMA median warming simulated over 2081-2100 with respect to 1995-2014 ranges from 1.9 [1.2 to 2.7] degrees C for SSP1-2.6 to 6.5 [4.9 to 9.0] degrees C for SSP5-8.5. This general warming is associated with a relative median snow cover extent decrease from -9.4 % [ – 16.4 % to -5.0 %] to -32.2 % [ -49.1 % to -25.0 %] and a relative median precipitation increase from 8.5 % [4.8 % to 18.2 %] to 24.9 % [14.4 % to 48.1 %] by the end of the century in these respective scenarios. The warming is 11 % higher over HMA than over the other Northern Hemisphere continental surfaces, excluding the Arctic area. Seasonal temperature, snow cover and precipitation changes over HMA show a linear relationship with the global surface air temperature (GSAT), except for summer snow cover which shows a slower decrease at strong levels of GSAT.

Abstract: Here, we define weather regimes in the East Antarctica-Southern Ocean sector based on daily anomalies of 700 hPa geopotential height derived from ERA5 reanalysis during 1979-2018. Most regimes and their preferred transitions depict synoptic-scale disturbances propagating eastwards off the Antarctic coastline. While regime sequences are generally short, their interannual variability is strongly driven by the polarity of the Southern Annular Mode (SAM). Regime occurrences are then intersected with atmospheric rivers (ARs) detected over the same region and period. ARs are equiprobable throughout the year, but clearly concentrate during regimes associated with a strong atmospheric ridges/blockings on the eastern part of the domain, which act to channel meridional advection of heat and moisture from the lower latitudes towards Antarctica. Both regimes and ARs significantly shape climate variability in Antarctica. Regimes favorable to AR occurrences are associated with anomalously warm and humid conditions in coastal Antarctica and, to a lesser extent, the hinterland parts of the Antarctic plateau. These anomalies are strongly enhanced during AR events, with warmer anomalies and dramatically amplified snowfall amounts. Large-scale conditions favoring AR development are finally explored. They show weak dependency to the SAM, but particularly strong atmospheric ridges/blockings over the Southern Ocean appear as the most favorable pattern, in which ARs can be embedded, and to which they contribute.

Abstract: This study presents the global climate model IPSL-CM6A-LR developed at Institut Pierre-Simon Laplace (IPSL) to study natural climate variability and climate response to natural and anthropogenic forcings as part of the sixth phase of the Coupled Model Intercomparison Project (CMIP6). This article describes the different model components, their coupling, and the simulated climate in comparison to previous model versions. We focus here on the representation of the physical climate along with the main characteristics of the global carbon cycle. The model's climatology, as assessed from a range of metrics (related in particular to radiation, temperature, precipitation, and wind), is strongly improved in comparison to previous model versions. Although they are reduced, a number of known biases and shortcomings (e.g., double Intertropical Convergence Zone [ITCZ], frequency of midlatitude wintertime blockings, and El Nino-Southern Oscillation [ENSO] dynamics) persist. The equilibrium climate sensitivity and transient climate response have both increased from the previous climate model IPSL-CM5A-LR used in CMIP5. A large ensemble of more than 30 members for the historical period (1850-2018) and a smaller ensemble for a range of emissions scenarios (until 2100 and 2300) are also presented and discussed.

Abstract: Permafrost is a ubiquitous phenomenon in the Arctic. Its future evolution is likely to control changes in northern high-latitude hydrology and biogeochemistry. Here we evaluate the permafrost dynamics in the global models participating in the Coupled Model Intercomparison Project (present generation – CMIP6; previous generation – CMIP5) along with the sensitivity of permafrost to climate change. Whilst the northern high-latitude air temperatures are relatively well simulated by the climate models, they do introduce a bias into any subsequent model estimate of permafrost. Therefore evaluation metrics are defined in relation to the air temperature. This paper shows that the climate, snow and permafrost physics of the CMIP6 multi-model ensemble is very similar to that of the CMIP5 multi-model ensemble. The main differences are that a small number of models have demonstrably better snow insulation in CMIP6 than in CMIP5 and a small number have a deeper soil profile. These changes lead to a small overall improvement in the representation of the permafrost extent. There is little improvement in the simulation of maximum summer thaw depth between CMIP5 and CMIP6. We suggest that more models should include a better-resolved and deeper soil profile as a first step towards addressing this. We use the annual mean thawed volume of the top 2 m of the soil defined from the model soil profiles for the permafrost region to quantify changes in permafrost dynamics. The CMIP6 models project that the annual mean frozen volume in the top 2 m of the soil could decrease by 10 %-40% degrees C-1 of global mean surface air temperature increase.

Abstract: Climate models and/or their output are usually bias-corrected for climate impact studies. The underlying assumption of these corrections is that climate biases are essentially stationary between historical and future climate states. Under very strong climate change, the validity of this assumption is uncertain, so the practical benefit of bias corrections remains an open question. Here, this issue is addressed in the context of bias correcting the climate models themselves. Employing the ARPEGE, LMDZ and CanAM4 atmospheric models, we undertook experiments in which one centre's atmospheric model takes another centre's coupled model as observations during the historical period, to define the bias correction, and as the reference under future projections of strong climate change, to evaluate its impact. This allows testing of the stationarity assumption directly from the historical through future periods for three different models. These experiments provide evidence for the validity of the new bias-corrected model approach. In particular, temperature, wind and pressure biases are reduced by 40-60% and, with few exceptions, more than 50% of the improvement obtained over the historical period is on average preserved after 100 years of strong climate change. Below 3 degrees C global average surface temperature increase, these corrections globally retain 80% of their benefit. Empirical bias corrections in climate models based on historical data improve future projections of climate change, even in strong change over 100 years, suggest experiments with three climate models.

Abstract: Regional changes in dry-season water availability over recent decades can be attributed to human-induced climate change, according to analyses of global reconstructions. Human-induced climate change impacts the hydrological cycle and thus the availability of water resources. However, previous assessments of observed warming-induced changes in dryness have not excluded natural climate variability and show conflicting results due to uncertainties in our understanding of the response of evapotranspiration. Here we employ data-driven and land-surface models to produce observation-based global reconstructions of water availability from 1902 to 2014, a period during which our planet experienced a global warming of approximately 1 degrees C. Our analysis reveals a spatial pattern of changes in average water availability during the driest month of the year over the past three decades compared with the first half of the twentieth century, with some regions experiencing increased and some decreased water availability. The global pattern is consistent with climate model estimates that account for anthropogenic effects, and it is not expected from natural climate variability, supporting human-induced climate change as the cause. There is regional evidence of drier dry seasons predominantly in extratropical latitudes and including Europe, western North America, northern Asia, southern South America, Australia and eastern Africa. We also find that the intensification of the dry season is generally a consequence of increasing evapotranspiration rather than decreasing precipitation.

Abstract: The Greenland Ice Sheet has been a major contributor to global sea-level rise in recent decades(1,2), and it is expected to continue to be so(3). Although increases in glacier flow(4-6) and surface melting(7-9) have been driven by oceanic(10-12) and atmospheric(13,14) warming, the magnitude and trajectory of the ice sheet's mass imbalance remain uncertain. Here we compare and combine 26 individual satellite measurements of changes in the ice sheet's volume, flow and gravitational potential to produce a reconciled estimate of its mass balance. The ice sheet was close to a state of balance in the 1990s, but annual losses have risen since then, peaking at 345 +/- 66 billion tonnes per year in 2011. In all, Greenland lost 3,902 +/- 342 billion tonnes of ice between 1992 and 2018, causing the mean sea level to rise by 10.8 +/- 0.9 millimetres. Using three regional climate models, we show that the reduced surface mass balance has driven 1,964 +/- 565 billion tonnes (50.3 per cent) of the ice loss owing to increased meltwater runoff. The remaining 1,938 +/- 541 billion tonnes (49.7 per cent) of ice loss was due to increased glacier dynamical imbalance, which rose from 46 +/- 37 billion tonnes per year in the 1990s to 87 +/- 25 billion tonnes per year since then. The total rate of ice loss slowed to 222 +/- 30 billion tonnes per year between 2013 and 2017, on average, as atmospheric circulation favoured cooler conditions(15) and ocean temperatures fell at the terminus of Jakobshavn Isbr AE(16). Cumulative ice losses from Greenland as a whole have been close to the rates predicted by the Intergovernmental Panel on Climate Change for their high-end climate warming scenario(17), which forecast an additional 70 to 130 millimetres of global sea-level rise by 2100 compared with their central estimate.

Abstract: Future sea surface temperature and sea-ice concentration from coupled ocean-atmosphere general circulation models such as those from the CMIP5 experiment are often used as boundary forcings for the downscaling of future climate experiments. Yet, these models show some considerable biases when compared to the observations over present climate. In this paper, existing methods such as an absolute anomaly method and a quantile-quantile method for sea surface temperature (SST) as well as a look-up table and a relative anomaly method for sea-ice concentration (SIC) are presented. For SIC, we also propose a new analogue method. Each method is objectively evaluated with a perfect model test using CMIP5 model experiments and some real-case applications using observations. We find that with respect to other previously existing methods, the analogue method is a substantial improvement for the bias correction of future SIC. Consistency between the constructed SST and SIC fields is an important constraint to consider, as is consistency between the prescribed sea-ice concentration and thickness; we show that the latter can be ensured by using a simple parameterisation of sea-ice thickness as a function of instantaneous and annual minimum SIC.

Abstract: This work presents snapshot simulations of the late 20th and late 21st century Antarctic climate under the RCP8.5 scenario carried out with an empirically bias-corrected global atmospheric general circulation model (AGCM), forced with bias-corrected sea-surface temperatures and sea ice and run with about 100-km resolution over Antarctica. The bias correction substantially improves the simulated mean late 20th century climate. The simulated atmospheric circulation of the bias-corrected model compares very favorably to the best available AMIP (Atmospheric Model Intercomparison Project)-type climate models. The simulated interannual circulation variability is improved by the bias correction. Depending on the metric, a slight improvement or degradation is found in the simulated variability on synoptic timescales. The simulated climate change over the 21st century is broadly similar in the corrected and uncorrected versions of the atmospheric model, and atmospheric circulation patterns are not geographically “pinned” by the applied bias correction. These results suggest that the method presented here can be used for bias-corrected climate projections. Finally, the authors discuss different possible choices in terms of the place of bias corrections and other intermediate steps in the modeling chain leading from global coupled climate simulations to impact assessment. Plain Language Summary Climate models are necessary and irreplaceable tools for climate projections, but despite continuous improvement, they still have biases, and their spatial resolution is too low to provide actionable climate change information at relevant small spatial scales. We present a method combining bias corrections and high-resolution climate modeling that allows improving climate projections at regional scales.

Abstract: The importance of northern peatlands in the global carbon cycle has been recognized, especially for long-term changes. Yet, the complex interactions between climate and peatland hydrology, carbon storage, and area dynamics make it challenging to represent these systems in land surface models. This study describes how peatlands are included as an independent sub-grid hydrological soil unit (HSU) in the ORCHIDEE-MICT land surface model. The peatland soil column in this tile is characterized by multilayered vertical water and carbon transport and peat-specific hydrological properties. The cost-efficient version of TOPMODEL and the scheme of peatland initiation and development from the DYPTOP model are implemented and adjusted to simulate spatial and temporal dynamics of peatland. The model is tested across a range of northern peatland sites and for gridded simulations over the Northern Hemisphere (> 30 degrees N). Simulated northern peatland area (3.9 million km(2)), peat carbon stock (463 Pg C), and peat depth are generally consistent with observed estimates of peatland area (3.4-4.0 million km2), peat carbon (270-540 Pg C), and data compilations of peat core depths. Our results show that both net primary production (NPP) and heterotrophic respiration (HR) of northern peatlands increased over the past century in response to CO2 and climate change. NPP increased more rapidly than HR, and thus net ecosystem production (NEP) exhibited a positive trend, contributing a cumulative carbon storage of 11.13 PgC since 1901, most of it being realized after the 1950s.

Abstract: Permafrost warming and potential soil carbon (SOC) release after thawing may amplify climate change, yet model estimates of present-day and future permafrost extent vary widely, partly due to uncertainties in simulated soil temperature. Here, we derive thermal diffusivity, a key parameter in the soil thermal regime, from depth-specific measurements of monthly soil temperature at about 200 sites in the high latitude regions. We find that, among the tested soil properties including SOC, soil texture, bulk density, and soil moisture, SOC is the dominant factor controlling the variability of diffusivity among sites. Analysis of the CMIP5 model outputs reveals that the parameterization of thermal diffusivity drives the differences in simulated present-day permafrost extent among these models. The strong SOC-thermics coupling is crucial for projecting future permafrost dynamics, since the response of soil temperature and permafrost area to a rising air temperature would be impacted by potential changes in SOC.

Abstract: Because all climate models exhibit biases, their use for assessing future climate change requires implicitly assuming or explicitly postulating that the biases are stationary or vary predictably. This hypothesis, however, has not been, and cannot be, tested directly. This work shows that under very large climate change the bias patterns of key climate variables exhibit a striking degree of stationarity. Using only correlation with a model's preindustrial bias pattern, a model's 4xCO(2) bias pattern is objectively and correctly identified among a large model ensemble in almost all cases. This outcome would be exceedingly improbable if bias patterns were independent of climate state. A similar result is also found for bias patterns in two historical periods. This provides compelling and heretofore missing justification for using such models to quantify climate perturbation patterns and for selecting well-performing models for regional downscaling. Furthermore, it opens the way to extending bias corrections to perturbed states, substantially broadening the range of justified applications of climate models.

Abstract: Widely present in boreal regions, peatlands contain large carbon stocks because of their hydrologic properties and high water content, which makes primary productivity exceed decomposition rates. We have enhanced the global land surface model ORCHIDEE by introducing a hydrological representation of northern peatlands. These peatlands are represented as a new plant functional type (PFT) in the model, with specific hydrological properties for peat soil. In this paper, we focus on the representation of the hydrology of northern peatlands and on the evaluation of the hydrological impact of this implementation. A prescribed map based on the inventory of Yu et al. (2010) defines peatlands as a fraction of a grid cell represented as a PFT comparable to C3 grasses, with adaptations to reproduce shallow roots and higher photosynthesis stress. The treatment of peatland hydrology differs from that of other vegetation types by the fact that runoff from other soil types is partially directed towards the peatlands (instead of directly to the river network). The evaluation of this implementation was carried out at different spatial and temporal scales, from site evaluation to larger scales such as the watershed scale and the scale of all northern latitudes. The simulated net ecosystem exchanges agree with observations from three FLUXNET sites. Water table positions were generally close to observations, with some exceptions in winter. Compared to other soils, the simulated peat soils have a reduced seasonal variability in water storage. The seasonal cycle of the simulated extent of inundated peatlands is compared to flooded area as estimated from satellite observations. The model is able to represent more than 89.5% of the flooded areas located in peatland areas, where the modelled extent of inundated peatlands reaches 0.83 x 10(6) km(2). However, the extent of peatlands in northern latitudes is too small to substantially impact the large-scale terrestrial water storage north of 45 degrees N. Therefore, the inclusion of peatlands has a weak impact on the simulated river discharge rates in boreal regions.

Abstract: High-latitude areas are very sensitive to global warming, which has significant impacts on soil temperatures and associated processes governing permafrost evolution. This study aims to improve first-layer soil temperature retrievals during winter. This key surface state variable is strongly affected by snow's geophysical properties and their associated uncertainties (e.g., thermal conductivity) in land surface climate models. We used infrared MODIS land-surface temperatures (LST) and Advanced Microwave Scanning Radiometer for EOS (AMSR-E) brightness temperatures (Tb) at 10.7 and 18.7 GHz to constrain the Canadian Land Surface Scheme (CLASS), driven by meteorological reanalysis data and coupled with a simple radiative transfer model. The Tb polarization ratio (horizontal/vertical) at 10.7 GHz was selected to improve snowpack density, which is linked to the thermal conductivity representation in the model. Referencing meteorological station soil temperature measurements, we validated the approach at four different sites in the North American tundra over a period of up to 8 years. Results show that the proposed method improves simulations of the soil temperature under snow (Tg) by 64% when using remote sensing (RS) data to constrain the model, compared to model outputs without satellite data information. The root mean square error (RMSE) between measured and simulated Tg under the snow ranges from 1.8 to 3.5 K when using RS data. Improved temporal monitoring of the soil thermal state, along with changes in snow properties, will improve our understanding of the various processes governing soil biological, hydrological, and permafrost evolution.

Abstract: Peatlands store substantial amounts of carbon and are vulnerable to climate change. We present a modified version of the Organising Carbon and Hydrology In Dynamic Ecosystems (ORCHIDEE) land surface model for simulating the hydrology, surface energy, and CO2 fluxes of peatlands on daily to annual timescales. The model includes a separate soil tile in each 0.5 degrees grid cell, defined from a global peatland map and identified with peat-specific soil hydraulic properties. Runoff from non-peat vegetation within a grid cell containing a fraction of peat is routed to this peat soil tile, which maintains shallow water tables. The water table position separates oxic from anoxic decomposition. The model was evaluated against eddy-covariance (EC) observations from 30 northern peatland sites, with the maximum rate of carboxylation (V-cmax) being optimized at each site. Regarding short-term day-to-day variations, the model performance was good for gross primary production (GPP) (r(2) = 0.76; Nash-Sutcliffe modeling efficiency, MEF = 0.76) and ecosystem respiration (ER, r(2) = 0.78, MEF = 0.75), with lesser accuracy for latent heat fluxes (LE, r(2) = 0.42, MEF = 0.14) and and net ecosystem CO2 exchange (NEE, r(2) = 0.38, MEF = 0.26). Seasonal variations in GPP, ER, NEE, and energy fluxes on monthly scales showed moderate to high r(2) values (0.57-0.86). For spatial across-site gradients of annual mean GPP, ER, NEE, and LE, r(2) values of 0.93, 0.89, 0.27, and 0.71 were achieved, respectively. Water table (WT) variation was not well predicted (r(2) < 0.1), likely due to the uncertain water input to the peat from surrounding areas. However, the poor performance of WT simulation did not greatly affect predictions of ER and NEE. We found a significant relationship between optimized V-cmax and latitude (temperature), which better reflects the spatial gradients of annual NEE than using an average V-cmax value.

Abstract: Large herbivores are a major agent in ecosystems, influencing vegetation structure, and carbon and nutrient flows. During the last glacial period, a mammoth steppe ecosystem prevailed in the unglaciated northern lands, supporting a high diversity and density of megafaunal herbivores. The apparent discrepancy between abundant megafauna and the expected low vegetation productivity under a generally harsher climate with a lower CO2 concentration, termed the productivity paradox, requires large-scale quantitative analysis using process-based ecosystem models. However, most of the current global dynamic vegetation models (DGVMs) lack explicit representation of large herbivores. Here we incorporated a grazing module in a DGVM based on physiological and demographic equations for wild large grazers, taking into account feedbacks of large grazers on vegetation. The model was applied globally for present-day and the Last Glacial Maximum (LGM). The present-day results of potential grazer biomass, combined with an empirical land-use map, infer a reduction in wild grazer biomass by 79-93% owing to anthropogenic land replacement of natural grasslands. For the LGM, we find that the larger mean body size of mammalian herbivores than today is the crucial clue to explain the productivity paradox, due to a more efficient exploitation of grass production by grazers with a large body size.

Abstract: It is important that climate models can accurately simulate the terrestrial carbon cycle in the Arctic due to the large and potentially labile carbon stocks found in permafrost-affected environments, which can lead to a positive climate feedback, along with the possibility of future carbon sinks from northward expansion of vegetation under climate warming. Here we evaluate the simulation of tundra carbon stocks and fluxes in three land surface schemes that each form part of major Earth system models (JSBACH, Germany; JULES, UK; ORCHIDEE, France). We use a site-level approach in which comprehensive, high-frequency datasets allow us to disentangle the importance of different processes. The models have improved physical permafrost processes and there is a reasonable correspondence between the simulated and measured physical variables, including soil temperature, soil moisture and snow. We show that if the models simulate the correct leaf area index (LAI), the standard C3 photosynthesis schemes produce the correct order of magnitude of carbon fluxes. Therefore, simulating the correct LAI is one of the first priorities. LAI depends quite strongly on climatic variables alone, as we see by the fact that the dynamic vegetation model can simulate most of the differences in LAI between sites, based almost entirely on climate inputs. However, we also identify an influence from nutrient limitation as the LAI becomes too large at some of the more nutrient-limited sites. We conclude that including moss as well as vascular plants is of primary importance to the carbon budget, as moss contributes a large fraction to the seasonal CO2 flux in nutrient-limited conditions. Moss photosynthetic activity can be strongly influenced by the moisture content of moss, and the carbon uptake can be significantly different from vascular plants with a similar LAI. The soil carbon stocks depend strongly on the rate of input of carbon from the vegetation to the soil, and our analysis suggests that an improved simulation of photosynthesis would also lead to an improved simulation of soil carbon stocks. However, the stocks are also influenced by soil carbon burial (e.g. through cryoturbation) and the rate of heterotrophic respiration, which depends on the soil physical state. More detailed below-ground measurements are needed to fully evaluate biological and physical soil processes. Furthermore, even if these processes are well modelled, the soil carbon profiles cannot resemble peat layers as peat accumulation processes are not represented in the models. Thus, we identify three priority areas for model development: (1) dynamic vegetation including (a) climate and (b) nutrient limitation effects; (2) adding moss as a plant functional type; and an (3) improved vertical profile of soil carbon including peat processes.

Abstract: We review recent literature on atmospheric, surface ocean and sea-ice observations and modeling results in the Antarctic sector and relate the observed climatic trends with the potential changes in the surface mass balance (SMB) of the ice sheet since 1900. Estimates of regional scale SMB distribution and trends remain subject to large uncertainties. Approaches combining and comparing multiple satellite and model-based assessments of ice sheet mass balance aim at reducing these knowledge gaps. During the last decades, significant changes in atmospheric circulation occurred around Antarctica, due to the exceptional positive trend in the Southern Annular Mode and to the climate variability observed in the tropical Pacific at the end of the twentieth century. Even though climate over the East Antarctic Ice-Sheet remained quite stable, a warming and precipitation increase was observed over the West Antarctic Ice-Sheet and over the West Antarctic Peninsula (AP) during the twentieth century. However, the high regional climate variability overwhelms climate changes associated to human drivers of global temperature changes, as reflected by a slight recent decadal cooling trend over the AP. Climate models still fail to accurately reproduce the multi-decadal SMB trends at a regional scale, and progress has to be achieved in reproducing atmospheric circulation changes related to complex ocean/ice/atmosphere interactions. Complex processes are also still insufficiently considered, such as (1) specific polar atmospheric processes (clouds, drifting snow, and stable boundary layer physics), (2) surface firn physics involved in the surface drag variations, or in firn air depletion and albedo feedbacks. Finally, progress in reducing the uncertainties relative to projections of the future SMB of Antarctica will largely depend on climate model capability to correctly consider teleconnections with low and mid-latitudes, and on the ability to correct them for biases, taking into account the coupling between ocean, ice, and atmosphere in high southern latitudes.

Abstract: Realistic projection of future climate-carbon (C) cycle feedbacks requires better understanding and an improved representation of the C cycle in permafrost regions in the current generation of Earth system models. Here we evaluated 10 terrestrial ecosystem models for their estimates of net primary productivity (NPP) and responses to historical climate change in permafrost regions in the Northern Hemisphere. In comparison with the satellite estimate from the Moderate Resolution Imaging Spectroradiometer (MODIS; 246 +/- 6gCm(-2) yr (-1)), most models produced higher NPP (309 +/- 12 g Cm-2 yr(-1)) over the permafrost region during 2000-2009. By comparing the simulated gross primary productivity (GPP) with a flux tower-based database, we found that although mean GPP among the models was only overestimated by 10% over 1982-2009, there was a twofold discrepancy among models (380 to 800 g Cm-2 yr(-1)), which mainly resulted from differences in simulated maximum monthly GPP (GPP(max)). Most models overestimated C use efficiency (CUE) as compared to observations at both regional and site levels. Further analysis shows that model variability of GPP and CUE are nonlinearly correlated to variability in specific leaf area and the maximum rate of carboxylation by the enzyme Rubisco at 25 degrees C (V-cmax_(25)), respectively. Themodels also varied in their sensitivities of NPP, GPP, and CUE to historical changes in climate and atmospheric CO2 concentration. These results indicate that model predictive ability of the C cycle in permafrost regions can be improved by better representation of the processes controlling CUE and GPP(max) as well as their sensitivity to climate change.

Abstract: The atmospheric concentration of CO2 increased from 190 to 280 ppm between the last glacial maximum 21,000 years ago and the pre-industrial era(1,2). This CO2 rise and its timing have been linked to changes in the Earth's orbit, ice sheet configuration and volume, and ocean carbon storage(2,3). The ice-core record of delta(CO2)-C-13 (refs 2,4) in the atmosphere can help to constrain the source of carbon, but previous modelling studies have failed to capture the evolution of delta(CO2)-C-13 over this period(5). Here we show that simulations of the last deglaciation that include a permafrost carbon component can reproduce the ice core records between 21,000 and 10,000 years ago. We suggest that thawing permafrost, due to increasing summer insolation in the northern hemisphere, is the main source of CO2 rise between 17,500 and 15,000 years ago, a period sometimes referred to as the Mystery Interval(6). Together with a fresh water release into the North Atlantic, much of the CO2 variability associated with the Bolling-Allerod/Younger Dryas period similar to 15,000 to similar to 12,000 years ago can also be explained. In simulations of future warming we find that the permafrost carbon feedback increases global mean temperature by 10-40% relative to simulations without this feedback, with the magnitude of the increase dependent on the evolution of anthropogenic carbon emissions.

Abstract: A significant portion of the large amount of carbon (C) currently stored in soils of the permafrost region in the Northern Hemisphere has the potential to be emitted as the greenhouse gases CO2 and CH4 under a warmer climate. In this study we evaluated the variability in the sensitivity of permafrost and C in recent decades among land surface model simulations over the permafrost region between 1960 and 2009. The 15 model simulations all predict a loss of near-surface permafrost (within 3m) area over the region, but there are large differences in the magnitude of the simulated rates of loss among the models (0.2 to 58.8x10(3)km(2)yr(-1)). Sensitivity simulations indicated that changes in air temperature largely explained changes in permafrost area, although interactions among changes in other environmental variables also played a role. All of the models indicate that both vegetation and soil C storage together have increased by 156 to 954TgCyr(-1) between 1960 and 2009 over the permafrost region even though model analyses indicate that warming alone would decrease soil C storage. Increases in gross primary production (GPP) largely explain the simulated increases in vegetation and soil C. The sensitivity of GPP to increases in atmospheric CO2 was the dominant cause of increases in GPP across the models, but comparison of simulated GPP trends across the 1982-2009 period with that of a global GPP data set indicates that all of the models overestimate the trend in GPP. Disturbance also appears to be an important factor affecting C storage, as models that consider disturbance had lower increases in C storage than models that did not consider disturbance. To improve the modeling of C in the permafrost region, there is the need for the modeling community to standardize structural representation of permafrost and carbon dynamics among models that are used to evaluate the permafrost C feedback and for the modeling and observational communities to jointly develop data sets and methodologies to more effectively benchmark models.

Abstract: The Land Surface, Snow and Soil Moisture Model Intercomparison Project (LS3MIP) is designed to provide a comprehensive assessment of land surface, snow and soil moisture feedbacks on climate variability and climate change, and to diagnose systematic biases in the land modules of current Earth system models (ESMs). The solid and liquid water stored at the land surface has a large influence on the regional climate, its variability and predictability, including effects on the energy, water and carbon cycles. Notably, snow and soil moisture affect surface radiation and flux partitioning properties, moisture storage and land surface memory. They both strongly affect atmospheric conditions, in particular surface air temperature and precipitation, but also large-scale circulation patterns. However, models show divergent responses and representations of these feedbacks as well as systematic biases in the underlying processes. LS3MIP will provide the means to quantify the associated uncertainties and better constrain climate change projections, which is of particular interest for highly vulnerable regions (densely populated areas, agricultural regions, the Arctic, semi-arid and other sensitive terrestrial ecosystems). The experiments are subdivided in two components, the first addressing systematic land biases in offline mode (“LMIP”, building upon the 3rd phase of Global Soil Wetness Project; GSWP3) and the second addressing land feedbacks attributed to soil moisture and snow in an integrated framework (“LFMIP”, building upon the GLACE-CMIP blueprint).

Abstract: Thawing of permafrost in a warming climate is governed by a complex interplay of different processes of which only conductive heat transfer is taken into account in most model studies. However, observations in many permafrost landscapes demonstrate that lateral and vertical movement of water can have a pronounced influence on the thaw trajectories, creating distinct landforms, such as thermokarst ponds and lakes, even in areas where permafrost is otherwise thermally stable. Novel process parameterizations are required to include such phenomena in future projections of permafrost thaw and subsequent climatic-triggered feedbacks. In this study, we present a new land-surface scheme designed for permafrost applications, CryoGrid 3, which constitutes a flexible platform to explore new parameterizations for a range of permafrost processes. We document the model physics and employed parameterizations for the basis module CryoGrid 3, and compare model results with in situ observations of surface energy balance, surface temperatures, and ground thermal regime from the Samoylov permafrost observatory in NE Siberia. The comparison suggests that CryoGrid 3 can not only model the evolution of the ground thermal regime in the last decade, but also consistently reproduce the chain of energy transfer processes from the atmosphere to the ground. In addition, we demonstrate a simple 1-D parameterization for thaw processes in permafrost areas rich in ground ice, which can phenomenologically reproduce both formation of thermokarst ponds and subsidence of the ground following thawing of ice-rich subsurface layers. Long-term simulation from 1901 to 2100 driven by reanalysis data and climate model output demonstrate that the hydrological regime can both accelerate and delay permafrost thawing. If meltwater from thawed ice-rich layers can drain, the ground subsides, as well as the formation of a talik, are delayed. If the meltwater pools at the surface, a pond is formed that enhances heat transfer in the ground and leads to the formation of a talik. The model results suggest that the trajectories of future permafrost thaw are strongly influenced by the cryostratigraphy, as determined by the late Quaternary history of a site.

Abstract: We use the BISICLES adaptive mesh ice sheet model to carry out one, two, and three century simulations of the fast-flowing ice streams of the West Antarctic Ice Sheet, deploying sub-kilometer resolution around the grounding line since coarser resolution results in substantial underestimation of the response. Each of the simulations begins with a geometry and velocity close to present-day observations, and evolves according to variation in meteoric ice accumulation rates and oceanic ice shelf melt rates. Future changes in accumulation and melt rates range from no change, through anomalies computed by atmosphere and ocean models driven by the El and A1B emissions scenarios, to spatially uniform melt rate anomalies that remove most of the ice shelves over a few centuries. We find that variation in the resulting ice dynamics is dominated by the choice of initial conditions and ice shelf melt rate and mesh resolution, although ice accumulation affects the net change in volume above flotation to a similar degree. Given sufficient melt rates, we compute grounding line retreat over hundreds of kilometers in every major ice stream, but the ocean models do not predict such melt rates outside of the Amundsen Sea Embayment until after 2100. Within the Amundsen Sea Embayment the largest single source of variability is the onset of sustained retreat in Thwaites Glacier, which can triple the rate of eustatic sea level rise.

Abstract: We present an approach to estimate the feedback from large-scale thawing of permafrost soils using a simplified, data-constrained model that combines three elements: soil carbon (C) maps and profiles to identify the distribution and type of C in permafrost soils; incubation experiments to quantify the rates of C lost after thaw; and models of soil thermal dynamics in response to climate warming. We call the approach the Permafrost Carbon Network Incubation-Panarctic Thermal scaling approach (PInc-PanTher). The approach assumes that C stocks do not decompose at all when frozen, but once thawed follow set decomposition trajectories as a function of soil temperature. The trajectories are determined according to a three-pool decomposition model fitted to incubation data using parameters specific to soil horizon types. We calculate litterfall C inputs required to maintain steady-state C balance for the current climate, and hold those inputs constant. Soil temperatures are taken from the soil thermal modules of ecosystem model simulations forced by a common set of future climate change anomalies under two warming scenarios over the period 2010 to 2100. Under a medium warming scenario (RCP4.5), the approach projects permafrost soil C losses of 12.2-33.4 Pg C; under a high warming scenario (RCP8.5), the approach projects C losses of 27.9-112.6 Pg C. Projected C losses are roughly linearly proportional to global temperature changes across the two scenarios. These results indicate a global sensitivity of frozen soil C to climate change (gamma sensitivity) of -14 to -19 PgC degrees C-1 on a 100 year time scale. For CH4 emissions, our approach assumes a fixed saturated area and that increases in CH4 emissions are related to increased heterotrophic respiration in anoxic soil, yielding CH4 emission increases of 7% and 35% for the RCP4.5 and RCP8.5 scenarios, respectively, which add an additional greenhouse gas forcing of approximately 10-18%. The simplified approach presented here neglects many important processes that may amplify or mitigate C release from permafrost soils, but serves as a data-constrained estimate on the forced, large-scale permafrost C response to warming.

Abstract: A warming climate is altering land-atmosphere exchanges of carbon, with a potential for increased vegetation productivity as well as the mobilization of permafrost soil carbon stores. Here we investigate land-atmosphere carbon dioxide (CO2) cycling through analysis of net ecosystem productivity (NEP) and its component fluxes of gross primary productivity (GPP) and ecosystem respiration (ER) and soil carbon residence time, simulated by a set of land surface models (LSMs) over a region spanning the drainage basin of Northern Eurasia. The retrospective simulations cover the period 1960-2009 at 0.5 degrees resolution, which is a scale common among many global carbon and climate model simulations. Model performance benchmarks were drawn from comparisons against both observed CO2 fluxes derived from site-based eddy covariance measurements as well as regional-scale GPP estimates based on satellite remote-sensing data. The site-based comparisons depict a tendency for overestimates in GPP and ER for several of the models, particularly at the two sites to the south. For several models the spatial pattern in GPP explains less than half the variance in the MODIS MOD17 GPP product. Across the models NEP increases by as little as 0.01 to as much as 0.79 g Cm-2 yr(-2), equivalent to 3 to 340% of the respective model means, over the analysis period. For the multimodel average the increase is 135% of the mean from the first to last 10 years of record (1960-1969 vs. 2000-2009), with a weakening CO2 sink over the latter decades. Vegetation net primary productivity increased by 8 to 30% from the first to last 10 years, contributing to soil carbon storage gains. The range in regional mean NEP among the group is twice the multimodel mean, indicative of the uncertainty in CO2 sink strength. The models simulate that inputs to the soil carbon pool exceeded losses, resulting in a net soil carbon gain amid a decrease in residence time. Our analysis points to improvements in model elements controlling vegetation productivity and soil respiration as being needed for reducing uncertainty in land-atmosphere CO2 exchange. These advances will require collection of new field data on vegetation and soil dynamics, the development of benchmarking data sets from measurements and remote-sensing observations, and investments in future model development and intercomparison studies.

Abstract: Processes that describe the distribution of vegetation and ecosystem succession after disturbance are an important component of dynamic global vegetation models (DGVMs). The vegetation dynamics module (ORC-VD) within the process-based ecosystem model ORCHIDEE (Organizing Carbon and Hydrology in Dynamic Ecosystems) has not been updated and evaluated since many years and is known to produce unrealistic results. This study presents a new parameterization of ORC-VD for mid- to high-latitude regions in the Northern Hemisphere, including processes that influence the existence, mortality and competition between tree functional types. A new set of metrics is also proposed to quantify the performance of ORC-VD, using up to five different data sets of satellite land cover, forest biomass from remote sensing and inventories, a data-driven estimate of gross primary productivity (GPP) and two gridded data sets of soil organic carbon content. The scoring of ORC-VD derived from these metrics integrates uncertainties in the observational data sets. This multi-data set evaluation framework is a generic method that could be applied to the evaluation of other DGVM models. The results of the original ORC-VD published in 2005 for mid- to high-latitudes and of the new parameterization are evaluated against the above-described data sets. Significant improvements were found in the modeling of the distribution of tree functional types north of 40 degrees N. Three additional sensitivity runs were carried out to separate the impact of different processes or drivers on simulated vegetation distribution, including soil freezing which limits net primary production through soil moisture availability in the root zone, elevated CO2 concentration since 1850, and the effects of frequency and severity of extreme cold events during the spin-up phase of the model.

Abstract: A variable-resolution atmospheric general circulation model (AGCM) is used for climate change projections over the Antarctic. The present-day simulation uses prescribed observed sea surface conditions, while a set of five simulations for the end of the twenty-first century (2070-99) under the Special Report on Emissions Scenarios (SRES) A1B scenario uses sea surface condition anomalies from selected coupled ocean atmosphere climate models from phase 3 of the Coupled Model Intercomparison Project (CMIP3). Analysis of the results shows that the prescribed sea surface condition anomalies have a very strong influence on the simulated climate change on the Antarctic continent, largely dominating the direct effect of the prescribed greenhouse gas concentration changes in the AGCM simulations. Complementary simulations with idealized forcings confirm these results. An analysis of circulation changes using self-organizing maps shows that the simulated climate change on regional scales is not principally caused by shifts of the frequencies of the dominant circulation patterns, except for precipitation changes in some coastal regions. The study illustrates that in some respects the use of bias-corrected sea surface boundary conditions in climate projections with a variable-resolution atmospheric general circulation model has some distinct advantages over the use of limited-area atmospheric circulation models directly forced by generally biased coupled climate model output.

Abstract: About 75 % of the Antarctic surface mass gain occurs over areas below 2,000 m asl, which cover 40 % of the grounded ice-sheet. As the topography is complex in many of these regions, surface mass balance modelling is highly dependent on horizontal resolution, and studying the impact of Antarctica on the future rise in sea level requires physical approaches. We have developed a computationally efficient, physical downscaling model for high-resolution (15 km) long-term surface mass balance (SMB) projections. Here, we present results of this model, called SMHiL (surface mass balance high-resolution downscaling), which was forced with the LMDZ4 atmospheric general circulation model to assess Antarctic SMB variability in the twenty first and the twenty second centuries under two different scenarios. The higher resolution of SMHiL better reproduces the geographical patterns of SMB and increase significantly the averaged SMB over the grounded ice-sheet for the end of the twentieth century. A comparison with more than 3200 quality-controlled field data shows that LMDZ4 and SMHiL reproduce the observed values equally well. Nevertheless, field data below 2,000 m asl are too scarce to efficiently show the added value of SMHiL and measuring the SMB in these undocumented areas should be a future scientific priority. Our results suggest that running LMDZ4 at a finer resolution (15 km) may give a future increase in SMB in Antarctica that is about 30 % higher than by using its standard resolution (60 km) due to the higher increase in precipitation in coastal areas at 15 km. However, a part (similar to 15 %) of these discrepancies could be an artefact from SMHiL since it neglects the foehn effect and likely overestimates the precipitation increase. Future changes in the Antarctic SMB at low elevations will result from the competition between higher snow accumulation and runoff. For this reason, developing downscaling models is crucial to represent processes in sufficient detail and correctly model the SMB in coastal areas.

Abstract: Given the growing evidence for megalakes in the geological record, assessing their impact on climate and vegetation is important for the validation of palaeoclimate simulations and therefore the accuracy of model-data comparison in lacustrine environments. Megalake Chad (MLC) occurrences are documented not only for the mid-Holocene but also for the Mio-Pliocene (Schuster et al., 2009). At this time, the surface covered by water would have reached up to similar to 350 000 km(2) (Ghienne et al., 2002; Schuster et al., 2005; Leblanc et al., 2006), making it an important evaporation source, possibly modifying climate and vegetation in the Chad Basin. We investigated the impact of such a giant continental water area in two different climatic backgrounds within the Paleoclimate Model Intercomparison Project phase 3 (PMIP3): the late Pliocene (3.3 to 3 Ma, i. e. the mid-Piacenzian warm period) and the mid-Holocene (6 kyr BP). In all simulations including MLC, precipitation is drastically reduced above the lake surface because deep convection is inhibited by overlying colder air. Meanwhile, convective activity is enhanced around MLC because of the wind increase generated by the flat surface of the megalake, transporting colder and moister air towards the eastern shore of the lake. The effect of MLC on precipitation and temperature is not sufficient to widely impact vegetation patterns. Nevertheless, tropical savanna is present in the Chad Basin in all climatic configurations, even without MLC presence, showing that the climate itself is the driver of favourable environments for sustainable hominid habitats.

Abstract: We present an updated and quality controlled surface mass balance (SMB) database for the Antarctic ice sheet. Importantly, the database includes formatted metadata, such as measurement technique, elevation, time covered, etc, which allows any user to filter out the data. Here, we discard data with limited spatial and temporal representativeness, too small measurement accuracy, or lack of quality control. Applied to the database, this filtering process gives four times more reliable data than when applied to previously available databases. New data with high spatial resolution are now available over long traverses, and at low elevation in some areas. However, the quality control led to a considerable reduction in the spatial density of data in several regions, particularly over West Antarctica. Over interior plateaus, where the SMB is low, the spatial density of measurements remains high. This quality controlled dataset was compared to results from ERA-Interim reanalysis to assess whether field data allow us to reconstruct an accurate description of the main SMB distribution features in Antarctica. We identified large areas where data gaps impede model validation: except for very few areas (e. g., Adelie Land), measurements in the elevation range between 200 m and 1000 m above sea level are not regularly distributed and do not allow a thorough validation of models in such regions with complex topography, where the highest scattering of SMB values is reported. Clearly, increasing the spatial density of field measurements at low elevations, in the Antarctic Peninsula and in West Antarctica is a scientific priority.

Abstract: Approximately 1700 Pg of soil carbon (C) are stored in the northern circumpolar permafrost zone, more than twice as much C than in the atmosphere. The overall amount, rate, and form of C released to the atmosphere in a warmer world will influence the strength of the permafrost C feedback to climate change. We used a survey to quantify variability in the perception of the vulnerability of permafrost C to climate change. Experts were asked to provide quantitative estimates of permafrost change in response to four scenarios of warming. For the highest warming scenario (RCP 8.5), experts hypothesized that C release from permafrost zone soils could be 19-45 Pg C by 2040, 162-288 Pg C by 2100, and 381-616 Pg C by 2300 in CO2 equivalent using 100-year CH4 global warming potential (GWP). These values become 50 % larger using 20-year CH4 GWP, with a third to a half of expected climate forcing coming from CH4 even though CH4 was only 2.3 % of the expected C release. Experts projected that two-thirds of this release could be avoided under the lowest warming scenario (RCP 2.6). These results highlight the potential risk from permafrost thaw and serve to frame a hypothesis about the magnitude of this feedback to climate change. However, the level of emissions proposed here are unlikely to overshadow the impact of fossil fuel burning, which will continue to be the main source of C emissions and climate forcing.

Abstract: The last glacial period has been punctuated by two types of abrupt climatic events, the Dansgaard-Oeschger (DO) and Heinrich (HE) events. These events, recorded in Greenland ice and in marine sediments, involved changes in the Atlantic Meridional Overturning Circulation (AMOC) and led to major changes in the terrestrial biosphere. Here we use the dynamical global vegetation model ORCHIDEE to simulate the response of vegetation to abrupt changes in the AMOC strength. We force ORCHIDEE offline with outputs from the IPSL_CM4 general circulation model, in which the AMOC is forced to change by adding freshwater fluxes in the North Atlantic. We investigate the impact of a collapse and recovery of the AMOC, at different rates, and focus on Western Europe, where many pollen records are available for comparison. The impact of an AMOC collapse on the European mean temperatures and precipitations simulated by the GCM is relatively small but sufficient to drive an important regression of forests and expansion of grasses in ORCHIDEE, in qualitative agreement with pollen data for an HE event. On the contrary, a run with a rapid shift of the AMOC to a hyperactive state of 30 Sv, mimicking the warming phase of a DO event, does not exhibit a strong impact on the European vegetation compared to the glacial control state. For our model, simulating the impact of an HE event thus appears easier than simulating the abrupt transition towards the interstadial phase of a DO. For both a collapse or a recovery of the AMOC, the vegetation starts to respond to climatic changes immediately but reaches equilibrium about 200 yr after the climate equilibrates, suggesting a possible bias in the climatic reconstructions based on pollen records, which assume equilibrium between climate and vegetation. However, our study does not take into account vegetation feedbacks on the atmosphere.

Abstract: Soil freezing is a major feature of boreal regions with substantial impact on climate. The present paper describes the implementation of the thermal and hydrological effects of soil freezing in the land surface model ORCHIDEE, which includes a physical description of continental hydrology. The new soil freezing scheme is evaluated against analytical solutions and in-situ observations at a variety of scales in order to test its numerical robustness, explore its sensitivity to parameterization choices and confront its performance to field measurements at typical application scales. Our soil freezing model exhibits a low sensitivity to the vertical discretization for spatial steps in the range of a few millimetres to a few centimetres. It is however sensitive to the temperature interval around the freezing point where phase change occurs, which should be 1 A degrees C to 2 A degrees C wide. Furthermore, linear and thermodynamical parameterizations of the liquid water content lead to similar results in terms of water redistribution within the soil and thermal evolution under freezing. Our approach does not allow firm discrimination of the performance of one approach over the other. The new soil freezing scheme considerably improves the representation of runoff and river discharge in regions underlain by permafrost or subject to seasonal freezing. A thermodynamical parameterization of the liquid water content appears more appropriate for an integrated description of the hydrological processes at the scale of the vast Siberian basins. The use of a subgrid variability approach and the representation of wetlands could help capture the features of the Arctic hydrological regime with more accuracy. The modeling of the soil thermal regime is generally improved by the representation of soil freezing processes. In particular, the dynamics of the active layer is captured with more accuracy, which is of crucial importance in the prospect of simulations involving the response of frozen carbon stocks to future warming. A realistic simulation of the snow cover and its thermal properties, as well as the representation of an organic horizon with specific thermal and hydrological characteristics, are confirmed to be a pre-requisite for a realistic modeling of the soil thermal dynamics in the Arctic.

Abstract: The high albedo of snow largely determines the climate of polar regions by controlling energy absorption at the surface. In Antarctica, where light-absorbing impurities are few, snow albedo is mostly determined by the size of snow grains(1). Snow metamorphism, the process of grain Coarsening, occurs at a rate that increases with temperature(2,3), so that snow albedo generally decreases as temperature increases. This increases energy absorption at the surface and atmospheric warming ensues, leading to a positive snow-albedo feedback. Here we use passive microwave satellite data and model outputs to show that this feedback is inhibited by small increases in precipitation. This is explained by the fact that grain coarsening in Antarctica is more sensitive to the deposition of small grains on the surface than previously assumed. We deduce that projected future increases in precipitation(4) can increase snow albedo by 0.4% on average during the twenty-first century and hence overcompensate the expected albedo decrease owing to warming (0.3% for 3 degrees C). Albedo-change projections in the Coupled Model Intercomparison Projects 3 and 5 do not reach a consensus on the sign and amplitude of this compensation, showing the need for a finer representation of the impact of precipitation on albedo in Antarctica.

Abstract: Predicting the climate for the future and how it will impact ice sheet evolution requires coupling ice sheet models with climate models. However, before we attempt to develop a realistic coupled setup, we propose, in this study, to first analyse the impact of a model simulated climate on an ice sheet. We undertake this exercise for a set of regional and global climate models. Modelled near surface air temperature and precipitation are provided as upper boundary conditions to the GRISLI (GRenoble Ice Shelf and Land Ice model) hybrid ice sheet model (ISM) in its Greenland configuration. After 20 kyrs of simulation, the resulting ice sheets highlight the differences between the climate models. While modelled ice sheet sizes are generally comparable to the observed one, there are considerable deviations among the ice sheets on regional scales. These deviations can be explained by biases in temperature and precipitation near the coast. This is especially true in the case of global models. But the deviations between the climate models are also due to the differences in the atmospheric general circulation. To account for these differences in the context of coupling ice sheet models with climate models, we conclude that appropriate down-scaling methods will be needed. In some cases, systematic corrections of the climatic variables at the interface may be required to obtain realistic results for the Greenland ice sheet (GIS).

Abstract: This work focuses on the Late Saalian (140 ka) Eurasian ice sheets' surface mass balance (SMB) sensitivity to changes in sea surface temperatures (SST). An Atmospheric General Circulation Model (AGCM), forced with two preexisting Last Glacial Maximum (LGM, 21 ka) SST reconstructions, is used to compute climate at 140 and 21 ka (reference glaciation). Contrary to the LGM, the ablation almost stopped at 140 ka due to the climatic cooling effect from the large ice sheet topography. Late Saalian SST are simulated using an AGCM coupled with a mixed layer ocean. Compared to the LGM, these 140 ka SST show an inter-hemispheric asymmetry caused by the larger ice-albedo feedback, cooling climate. The resulting Late Saalian ice sheet SMB is smaller due to the extensive simulated sea ice reducing the precipitation. In conclusion, SST are important for the stability and growth of the Late Saalian Eurasian ice sheet.

Abstract: Permafrost soils contain enormous amounts of organic carbon, which could act as a positive feedback to global climate change due to enhanced respiration rates with warming. We have used a terrestrial ecosystem model that includes permafrost carbon dynamics, inhibition of respiration in frozen soil layers, vertical mixing of soil carbon from surface to permafrost layers, and CH(4) emissions from flooded areas, and which better matches new circumpolar inventories of soil carbon stocks, to explore the potential for carbon-climate feedbacks at high latitudes. Contrary to model results for the Intergovernmental Panel on Climate Change Fourth Assessment Report (IPCC AR4), when permafrost processes are included, terrestrial ecosystems north of 60 degrees N could shift from being a sink to a source of CO(2) by the end of the 21st century when forced by a Special Report on Emissions Scenarios ( SRES) A2 climate change scenario. Between 1860 and 2100, the model response to combined CO(2) fertilization and climate change changes from a sink of 68 Pg to a 27 + -7 Pg sink to 4 + -18 Pg source, depending on the processes and parameter values used. The integrated change in carbon due to climate change shifts from near zero, which is within the range of previous model estimates, to a climate-induced loss of carbon by ecosystems in the range of 25 + -3 to 85 + -16 Pg C, depending on processes included in the model, with a best estimate of a 62 + -7 Pg C loss. Methane emissions from high-latitude regions are calculated to increase from 34 Tg CH(4)/y to 41-70 TgCH(4)/y, with increases due to CO(2) fertilization, permafrost thaw, and warming-induced increased CH(4) flux densities partially offset by a reduction in wetland extent.

Abstract: The examination of more than 1500 paleohydrological dated records collected between 10 and 28 degrees N during the last 50 years have been used to improve our knowledge and understanding of the Sahara and Sahel vulnerability to the Atlantic monsoon changes in the long-term. We have analyzed the distribution of water bodies (mainly lakes and wetlands) over time and space: the central Saharan massifs played a major role in favoring water supply to the lowlands throughout the whole African Humid Period. In addition, distinct East-West dynamics is recorded with humidity starting – and stopping – several millennia earlier to the east than to the west of the Sahara. A series of time lags are discussed: (1) between the maximum of deep (fresh water) lake formation during the early Holocene and the maximum of water body extensions during the mid-Holocene which highlight the primary role of aquifer water level in lake response to climate change (2) between the hydrological history of the Sahara and the Sahel and the forcings – mainly insolation changes – during the early and mid-Holocene which involves complex interactions between remnant ice sheet in the Northern Hemisphere, open water bodies in the Sahara and Sahel and the Atlantic monsoon system. (C) 2011 Elsevier Ltd. All rights reserved.

Abstract: We compare the present and last interglacial periods as recorded in Antarctic water stable isotope records now available at various temporal resolutions from six East Antarctic ice cores: Vostok, Taylor Dome, EPICA Dome C (EDC), EPICA Dronning Maud Land (EDML), Dome Fuji and the recent TALDICE ice core from Talos Dome. We first review the different modern site characteristics in terms of ice flow, meteorological conditions, precipitation intermittency and moisture origin, as depicted by meteorological data, atmospheric reanalyses and Lagrangian moisture source diagnostics. These different factors can indeed alter the relationships between temperature and water stable isotopes. Using five records with sufficient resolution on the EDC3 age scale, common features are quantified through principal component analyses. Consistent with instrumental records and atmospheric model results, the ice core data depict rather coherent and homogenous patterns in East Antarctica during the last two interglacials. Across the East Antarctic plateau, regional differences, with respect to the common East Antarctic signal, appear to have similar patterns during the current and last interglacials. We identify two abrupt shifts in isotopic records during the glacial inception at TALDICE and EDML, likely caused by regional sea ice expansion. These regional differences are discussed in terms of moisture origin and in terms of past changes in local elevation histories, which are compared to ice sheet model results. Our results suggest that elevation changes may contribute significantly to inter-site differences. These elevation changes may be underestimated by current ice sheet models.

Abstract: Vegetation reconstructions from pollen data for the Last Glacial Maximum (LGM), 21 ky ago, reveal lanscapes radically different from the modern ones, with, in particular, a massive regression of forested areas in both hemispheres. Two main factors have to be taken into account to explain these changes in comparison to today's potential vegetation: a generally cooler and drier climate and a lower level of atmospheric CO(2). In order to assess the relative impact of climate and atmospheric CO(2) changes on the global vegetation, we simulate the potential modern vegetation and the glacial vegetation with the dynamical global vegetation model ORCHIDEE, driven by outputs from the IPSLCM4v1 atmosphere-ocean general circulation model, under modern or glacial CO(2) levels for photosynthesis. ORCHIDEE correctly reproduces the broad features of the glacial vegetation. Our modelling results support the view that the physiological effect of glacial CO(2) is a key factor to explain vegetation changes during glacial times. In our simulations, the low atmospheric CO(2) is the only driver of the tropical forests regression, and explains half of the response of temperate and boreal forests to glacial conditions. Our study shows that the sensitivity to CO(2) changes depends on the background climate over a region, and also depends on the vegetation type, needleleaf trees being much more sensitive than broadleaf trees in our model. This difference of sensitivity leads to a dominance of broadleaf types in the remaining simulated forests, which is not supported by pollen data, but nonetheless suggests a potential impact of CO(2) on the glacial vegetation assemblages. It also modifies the competitivity between the trees and makes the amplitude of the response to CO(2) dependent on the initial vegetation state.

Abstract: This study evaluates the climatic impact of possible future changes in high-latitude inland water surface (IWS) area. We carried out a set of climate-change experiments with an atmospheric general circulation model in which different scenarios of future changes of IWS extent were prescribed. The simulations are based on the SRES-A1B greenhouse gas emission scenario and represent the transient climatic state at the end of the 21st century. Our results indicate that the impact of a reduction in IWS extent depends on the season considered: the total disappearance of IWS would lead to cooling during cold seasons and to warming in summer. In the annual mean, the cooling effect would be dominant. In an experiment in which the future change of prescribed IWS extent is prescribed as a function of the simulated changes of permafrost extent, we find that these changes are self-consistent in the sense that their effects on the simulated temperature and precipitation patterns would not be contradictory to the underlying scenario of changes in IWS extent. In this “best guess” simulation, the projected changes in IWS extent would reduce future near-surface warn-ling over large parts of northern Eurasia by about 20% during the cold season, while the impact in North America and during summer is less clear. As a whole, the direct climatic impact of future IWS changes is likely to be moderate.

Abstract: This paper describes atmospheric general circulation model climate change experiments in which the Arctic sea-ice thickness is either fixed to 3 m or somewhat more realistically parameterized in order to take into account essentially the spatial variability of Arctic sea-ice thickness, which is, to a first approximation, a function of ice type (perennial or seasonal). It is shown that, both at present and at the end of the twenty-first century (under the SRES-A1B greenhouse gas scenario), the impact of a variable sea-ice thickness compared to a uniform value is essentially limited to the cold seasons and the lower troposphere. However, because first-year ice is scarce in the Central Arctic today, but not under SRES-A1B conditions at the end of the twenty-first century, and because the impact of a sea-ice thickness reduction can be masked by changes of the open water fraction, the spatial and temporal patterns of the effect of sea-ice thinning on the atmosphere differ between the two periods considered. As a consequence, not only the climate simulated at a given period, but also the simulated Arctic climate change over the twenty-first century is affected by the way sea-ice thickness is prescribed.

Abstract: This paper presents the major characteristics of the Institut Pierre Simon Laplace (IPSL) coupled ocean-atmosphere general circulation model. The model components and the coupling methodology are described, as well as the main characteristics of the climatology and interannual variability. The model results of the standard version used for IPCC climate projections, and for intercomparison projects like the Paleoclimate Modeling Intercomparison Project (PMIP 2) are compared to those with a higher resolution in the atmosphere. A focus on the North Atlantic and on the tropics is used to address the impact of the atmosphere resolution on processes and feedbacks. In the North Atlantic, the resolution change leads to an improved representation of the storm-tracks and the North Atlantic oscillation. The better representation of the wind structure increases the northward salt transports, the deep-water formation and the Atlantic meridional overturning circulation. In the tropics, the ocean-atmosphere dynamical coupling, or Bjerknes feedback, improves with the resolution. The amplitude of ENSO (El Nio-Southern oscillation) consequently increases, as the damping processes are left unchanged.

Abstract: This work uses an atmospheric general circulation model (AGCM) asynchronously coupled to an equilibrium vegetation model to investigate whether vegetation feedbacks could be one of the reasons why the Late Saalian ice sheet (140 kyrs BP) in Eurasia was substantially larger than the Last Glacial Maximum(LGM, 21 kyrs BP) Eurasian ice sheet. The modeled vegetation changes induce a regional cooling for the Late Saalian while they cause a slight regional warming for LGM. As a result, ablation along the margins of the Late Saalian ice sheet is significantly reduced, leading to an increased surface mass balance, while there are no significant mass balance changes observed from vegetation feedbacks at LGM. Citation: Colleoni, F., G. Krinner, and M. Jakobsson (2009), Sensitivity of the Late Saalian (140 kyrs BP) and LGM (21 kyrs BP) Eurasian ice sheet surface mass balance to vegetation feedbacks, Geophys. Res. Lett., 36, L08704, doi:10.1029/2009GL037200.

Abstract: All climate models participating in the Fourth Assessment Report of the Intergovernmental Panel on Climate Change, as made available by the Program for Climate Model Diagnosis and Intercomparison (PCMDI) as the Coupled Model Intercomparison Project 3 (CMIP3) archive, predict a significant surface warming of Antarctica by the end of the 21st century under a moderate (SRESA1B) greenhouse-gas scenario. All models but one predict a concurrent precipitation increase but with a large scatter of results. The models with finer horizontal resolution tend to predict a larger precipitation increase. Because modeled Antarctic surface mass balance is known to be sensitive to horizontal resolution, extrapolating predictions from the different models with respect to model resolution may provide simple yet better multi-model estimates of Antarctic precipitation change than mere averaging or even more complex approaches. Using such extrapolation, a conservative estimate of the predicted precipitation increase at the end of the 21st century is +30 kg m(-2) a(-1) on the grounded ice sheet, corresponding to a >1 mm a(-1) sea-level rise. About three-quarters of this rise originates from the marginal regions of the Antarctic ice sheet with surface elevation below 2250 m. This is where field programs are most urgently needed to better understand and monitor accumulation at the surface of Antarctica, and to improve and verify prediction models.

Abstract: Here we perform a detailed comparison between climate model results and climate reconstructions in western Europe and the Mediterranean area for the mid-Piacenzian warm interval (ca 3 Myr ago) of the Late Pliocene epoch. This region is particularly well suited for such a comparison as several quantitative climate estimates from local pollen records are available. They show evidence for temperatures significantly warmer than today over the whole area, mean annual precipitation higher in northwestern Europe and equivalent to modern values in its southwestern part. To improve our comparison, we have performed high resolution simulations of the mid-Piacenzian climate using the LMDz atmospheric general circulation model (AGCM) with a stretched grid which allows a finer resolution over Europe. In a first step, we applied the PRISM2 (Pliocene Research, Interpretation, and Synoptic Mapping) boundary conditions except that we used modern terrestrial vegetation. Second, we simulated the vegetation for this period by forcing the ORCHIDEE (Organizing Carbon and Hydrology in Dynamic Ecosystems) dynamic global vegetation model (DGVM) with the climatic outputs from the AGCM. We then supplied this simulated terrestrial vegetation cover as an additional boundary condition in a second AGCM run. This gives us the opportunity to investigate the model's sensitivity to the simulated vegetation changes in a global warming context. Model results and data show a great consistency for mean annual temperatures, indicating increases by up to 4 degrees C in the study area, and some disparities, in particular in the northern Mediterranean sector, as regards winter and summer temperatures. Similar continental mean annual precipitation and moisture patterns are predicted by the model, which broadly underestimates the wetter conditions indicated by the data in northwestern Europe. The biogeophysical effects due to the changes in vegetation simulated by ORCHIDEE are weak, both in terms of the hydrological cycle and of the temperatures, at the regional scale of the European and Mediterranean mid-latitudes. In particular, they do not contribute to improve the model-data comparison. Their main influence concerns seasonal temperatures, with a decrease of the temperatures of the warmest month, and an overall reduction of the intensity of the continental hydrological cycle.

Abstract: We modify the soil component of the ORCHIDEE terrestrial carbon cycle model to include vertically-discretized soil carbon. With this model, we investigate the feedback of considering thermal insulation by soil carbon, which modifies the soil thermal regime by lowering the thermal conductivity and increasing the heat capacity of a carbon-rich soil, on the total carbon stocks the model builds up. In addition, we demonstrate the effects of diffusive vertical mixing of soil organic matter by cryoturbation on the total carbon stocks that the model builds up in mineral soils in equilibrium with a steady climate. We show that including these two effects together leads to up to 30% higher soil carbon stocks in the top meter of permafrost soils, as well as large stocks of carbon below 1m in the upper permafrost soil layers. The vertical profile of partitioning of carbon between different lability pools is also affected, as the slower pools are more deeply mixed; also the time to reach equilibrium lengthens considerably. These effects are largest in the coldest regions such as Eastern Siberia. The inclusion of cryoturbative mixing and insulation by soil carbon leads to better agreement with estimates of high-latitude soil carbon stocks, where substantial amounts of carbon are found in permafrost regions, to depths of three meters; however we do not include peat, Yedoma, or alluvial deposition processes here, so the total carbon stocks are still lower than observed. Citation: Koven, C., P. Friedlingstein, P. Ciais, D. Khvorostyanov, G. Krinner, and C. Tarnocai (2009), On the formation of high-latitude soil carbon stocks: Effects of cryoturbation and insulation by organic matter in a land surface model, Geophys. Res. Lett., 36, L21501, doi:10.1029/2009GL040150.

Abstract: In the companion paper (Part I), we presented a model of permafrost carbon cycle to study the sensitivity of frozen carbon stocks to future climate warming. The mobilization of deep carbon stock of the frozen Pleistocene soil in the case of rapid stepwise increase of atmospheric temperature was considered. In this work, we adapted the model to be used also for floodplain tundra sites and to account for the processes in the soil active layer. The new processes taken into account are litter input and decomposition, plant-mediated transport of methane, and leaching of exudates from plant roots. The SRES-A2 transient climate warming scenario of the IPSL CM4 climate model is used to study the carbon fluxes from the carbon-rich Pleistocene soil with seasonal active-layer carbon cycling on top of it. For a point to the southwest from the western branch of Yedoma Ice Complex, where the climate warming is strong enough to trigger self-sustainable decomposition processes, about 256 kgC m(-2), or 70% of the initial soil carbon stock under present-day climate conditions, are emitted to the atmosphere in about 120 yr, including 20 kgC m(-2) released as methane. The total average flux of CO2 and methane emissions to the atmosphere during this time is of 2.1 kgC m(-2) yr(-1). Within the Yedoma, whose most part of the territory remains relatively cold, the emissions are much smaller: 0.2 kgC m(-2) yr(-1) between 2050 and 2100 for Yakutsk area. In a test case with saturated upper-soil meter, when the runoff is insufficient to evacuate the meltwater, 0.05 kgCH(4) m(-2) yr(-1) on average are emitted as methane during 250 yr starting from 2050. The latter can translate to the upper bound of 1 GtC yr(-1) in CO2 equivalent from the 1 million km(2) area of the Yedoma.

Abstract: This article reports on high-resolution (60 km) atmospheric general circulation model simulations of the Antarctic climate for the periods 1981-2000 and 2081-2100. The analysis focuses on the surface mass balance change, one of the components of the total ice sheet mass balance, and its impact on global eustatic sea level. Contrary to previous simulations, in which the authors directly used sea surface boundary conditions produced by a coupled ocean-atmosphere model for the last decades of both centuries, an anomaly method was applied here in which the present-day simulations use observed sea surface conditions, while the simulations for the end of the twenty-first century use the change in sea surface conditions taken from the coupled simulations superimposed on the present-day observations. It is shown that the use of observed oceanic boundary conditions clearly improves the simulation of the present-day Antarctic climate, compared to model runs using boundary conditions from a coupled climate model. Moreover, although the spatial patterns of the simulated climate change are similar, the two methods yield significantly different estimates of the amplitude of the future climate and surface mass balance change over the Antarctic continent. These differences are of similar magnitude as the intermodel dispersion in the current Intergovernmental Panel on Climate Change (IPCC) exercise: selecting a method for generating boundary conditions for a high-resolution model may be just as important as selecting the climate model itself. Using the anomaly method, the simulated mean surface mass balance change over the grounded ice sheet from 1981-2000 to 2081-2100 is 43-mm water equivalent per year, corresponding to a eustatic sea level decrease of 1.5 mm yr(-1). A further result of this work is that future continental-mean surface mass balance changes are dominated by the coastal regions, and that high-resolution models, which better resolve coastal processes, tend to predict stronger precipitation changes than models with lower spatial resolution.

Abstract: A record of accumulation and ablation from a network of 47 stakes at a coastal blue-ice area in Terre Adelie, Antarctica, is presented and analyzed. The record covers early 2004 to early 2006, from 25 field surveys including some in austral winter. The two years are very different, with a virtually null surface mass balance during the 2004 winter but large accumulation during the 2005 winter. A snow/ice energy- and mass-balance model is used to reproduce the accumulation and ablation record. A parameterization for snow erosion by wind is included. Input meteorology is from the European Centre for Medium-Range Weather Forecasts (ECMWF) analyses and forecasts, corrected using I year of local meteorological observations from an automatic weather station. Model results agree reasonably well with the observations. Wind erosion is the largest contributor to ablation, removing much of the precipitation. Sublimation and, to a lesser extent, melt/runoff together account for >60 cm w.e. of ablation in 2 years, mainly in summer. Although the record is short, it confirms high interannual variability and thus high sensitivity to meteorology and climate. Monitoring and understanding the mass balance of such coastal blue-ice areas may help monitor and detect climate change in the Antarctic coastal regions.

Abstract: [1] On the basis of thousands of surface mass balance (SMB) field measurements over the entire Antarctic ice sheet it is currently estimated that more than 2 Gt of ice accumulate each year at the surface of Antarctica. However, these estimates suffer from large uncertainties. Various problems affect Antarctic SMB measurements, in particular, limited or unwarranted spatial and temporal representativeness, measurement inaccuracy, and lack of quality control. We define quality criteria on the basis of ( 1) an up-to-date review and quality rating of the various SMB measurement methods and ( 2) essential information ( location, dates of measurements, time period covered by the SMB values, and primary data sources) related to each SMB data. We apply these criteria to available SMB values from Queen Mary to Victoria lands (90 degrees – 180 degrees E Antarctic sector) from the early 1950s to present. This results in a new set of observed SMB values for the 1950 – 2005 time period with strong reduction in density and coverage but also expectedly reduced inaccuracies and uncertainties compared to other compilations. The quality-controlled SMB data set also contains new results from recent field campaigns ( International Trans-Antarctic Scientific Expedition (ITASE), Russian Antarctic Expedition (RAE), and Australian National Antarctic Research Expeditions (ANARE) projects) which comply with the defined quality criteria. A comparative evaluation of climate model results against the quality-controlled updated SMB data set and other widely used ones illustrates that such Antarctic SMB studies are significantly affected by the quality of field SMB values used as reference.

Abstract: The Last Glacial Maximum has been one of the first foci of the Paleoclimate Modelling Intercomparison Project (PMIP). During its first phase, the results of 17 atmosphere general circulation models were compared to paleoclimate reconstructions. One of the largest discrepancies in the simulations was the systematic underestimation, by at least 10 degrees C, of the winter cooling over Europe and the Mediterranean region observed in the pollen-based reconstructions. In this paper, we investigate the progress achieved to reduce this inconsistency through a large modelling effort and improved temperature reconstructions. We show that increased model spatial resolution does not significantly increase the simulated LGM winter cooling. Further, neither the inclusion of a vegetation cover compatible with the LGM climate, nor the interactions with the oceans simulated by the atmosphere-ocean general circulation models run in the second phase of PMIP result in a better agreement between models and data. Accounting for changes in interannual variability in the interpretation of the pollen data does not result in a reduction of the reconstructed cooling. The largest recent improvement in the model-data comparison has instead arisen from a new climate reconstruction based on inverse vegetation modelling, which explicitly accounts for the CO2 decrease at LGM and which substantially reduces the LGM winter cooling reconstructed from pollen assemblages. As a result, the simulated and observed LGM winter cooling over Western Europe and the Mediterranean area are now in much better agreement.

Abstract: The recent number isotopic records extracted from Andean ice cores ( South America) has illustrated the key role such archives can play in past climate reconstructions. Nevertheless, interpreting isotopic archives as quantified climate proxies requires an understanding of which climate parameters control the stable isotopic composition of water. Mesoscale modeling sheds new light on the meteorological mechanisms dominant during austral summer. Here we focus on the variability of the South Atlantic Convergence Zone ( SACZ) and its repercussions on upstream regions. The SACZ is a major component of the South American Monsoon System ( SAMS). The present study uses the isotopic signature of the SAMS, as simulated by the stable water isotope enabled regional circulation model REMOiso to answer the question: how does the SAMS affect the isotopic composition of precipitation during the wet season? In order to analyze the internal, purely atmospheric variability mode, the model was forced by climatological sea-surface temperatures. We investigate the concurrent intraseasonal variability of meteorological and isotopic parameters at pentad ( 5 days) interval using empirical orthogonal functions ( EOFs). REMOiso reproduces the main meteorological characteristics of the SAMS consistent with observations as well as previous modeling studies. Furthermore, we demonstrate that delta O-18 integrates both circulation and precipitation variability. This new evidence contributes to the comprehension of the delta O-18 signal in tropical South America, highlighting the internal atmospheric variability, as opposed to external forcing by Pacific and Atlantic sea-surface temperature.

Abstract: The LMDZ4 general circulation model is the atmospheric component of the IPSL-CM4 coupled model which has been used to perform climate change simulations for the 4th IPCC assessment report. The main aspects of the model climatology (forced by observed sea surface temperature) are documented here, as well as the major improvements with respect to the previous versions, which mainly come form the parametrization of tropical convection. A methodology is proposed to help analyse the sensitivity of the tropical Hadley-Walker circulation to the parametrization of cumulus convection and clouds. The tropical circulation is characterized using scalar potentials associated with the horizontal wind and horizontal transport of geopotential (the Laplacian of which is proportional to the total vertical momentum in the atmospheric column). The effect of parametrized physics is analysed in a regime sorted framework using the vertical velocity at 500 hPa as a proxy for large scale vertical motion. Compared to Tiedtke's convection scheme, used in previous versions, the Emanuel's scheme improves the representation of the Hadley-Walker circulation, with a relatively stronger and deeper large scale vertical ascent over tropical continents, and suppresses the marked patterns of concentrated rainfall over oceans. Thanks to the regime sorted analyses, these differences are attributed to intrinsic differences in the vertical distribution of convective heating, and to the lack of self-inhibition by precipitating downdraughts in Tiedtke's parametrization. Both the convection and cloud schemes are shown to control the relative importance of large scale convection over land and ocean, an important point for the behaviour of the coupled model.

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

Abstract: Future climate warming is expected to enhance plant growth in temperate ecosystems and to increase carbon sequestration(1,2). But although severe regional heatwaves may become more frequent in a changing climate(3,4), their impact on terrestrial carbon cycling is unclear. Here we report measurements of ecosystem carbon dioxide fluxes, remotely sensed radiation absorbed by plants, and country- level crop yields taken during the European heatwave in 2003. We use a terrestrial biosphere simulation model(5) to assess continental- scale changes in primary productivity during 2003, and their consequences for the net carbon balance. We estimate a 30 per cent reduction in gross primary productivity over Europe, which resulted in a strong anomalous net source of carbon dioxide ( 0.5 Pg Cyr(-1)) to the atmosphere and reversed the effect of four years of net ecosystem carbon sequestration(6). Our results suggest that productivity reduction in eastern and western Europe can be explained by rainfall deficit and extreme summer heat, respectively. We also find that ecosystem respiration decreased together with gross primary productivity, rather than accelerating with the temperature rise. Model results, corroborated by historical records of crop yields, suggest that such a reduction in Europe's primary productivity is unprecedented during the last century. An increase in future drought events could turn temperate ecosystems into carbon sources, contributing to positive carbon- climate feedbacks already anticipated in the tropics and at high latitudes(1,2).

Abstract: The Heinrich Event 1, the most recent of the glacial North Atlantic large iceberg discharges, is well documented in continental and marine records, but this large perturbation of the climate system has rarely been simulated. Here we propose a preliminary model-data comparison for this period, which we compare to the Last Glacial Maximum state. The pollen record from one specific core from the western Mediterranean Sea (ODP site 976) is analysed both in terms of vegetation distribution and climatic implication. The climate and vegetation of both periods are then simulated and compared to the pollen-based data.

Abstract: Based on an interactive global biogenic emission and dynamic vegetation model, we investigate the evolution of volatile organic compound (VOC) emissions by the terrestrial biosphere in four scenarios: the Last Glacial Maximum (21,000 years BP), the preindustrial (1850s), present-day (1990s) and the future (2100). The combined effects of foliar expansion, climate change and ecosystems redistribution impact strongly on biogenic emissions. Total biogenic VOC emissions increase from 331 TgC/yr at the LGM to 702 TgC/yr at the preindustrial, 725 TgC/yr at present-day and to 1251 TgC/yr under future conditions. If the tropics remain a major source region, a substantial decrease in VOC emissions is calculated over Amazonia for 2100 due to the recession of tropical forests in response to climate change. The Northern Hemisphere becomes a significant source of VOC in the future and globally, emissions increase by 27% for isoprene and 51% for monoterpenes compared to the present.

Abstract: Large proglacial lakes cool regional summer climate because of their large heat capacity, and have been shown to modify precipitation through mesoscale atmospheric feedbacks, as in the case of Lake Agassiz(1). Several large ice-dammed lakes, with a combined area twice that of the Caspian Sea, were formed in northern Eurasia about 90,000 years ago, during the last glacial period when an ice sheet centred over the Barents and Kara seas(2) blocked the large northbound Russian rivers(3). Here we present high-resolution simulations with an atmospheric general circulation model that explicitly simulates the surface mass balance of the ice sheet. We show that the main influence of the Eurasian proglacial lakes was a significant reduction of ice sheet melting at the southern margin of the Barents – Kara ice sheet through strong regional summer cooling over large parts of Russia. In our simulations, the summer melt reduction clearly outweighs lake-induced decreases in moisture and hence snowfall, such as has been reported earlier for Lake Agassiz1. We conclude that the summer cooling mechanism from proglacial lakes accelerated ice sheet growth and delayed ice sheet decay in Eurasia and probably also in North America.

Abstract: We introduced a simple scheme of soil freezing in the LMDz3.3 atmospheric general circulation model (AGCM) to examine the potential effects of this parameterization on simulated future boreal climate change. In this multi-layer soil scheme, soil heat capacity and conductivity are dependent on soil water content, and a parameterization of the thermal and hydrological effects of water phase changes is included. The impact of these new features is evaluated against observations. By comparing present-day and 2 x CO2 AGCM simulations both with and without the parameterization of soil freezing the role of soil freezing in climate change is analysed. Soil freezing does not have significant global impacts, but regional effects on simulated climate and climate change are important. In present-day conditions, hydrological effects due to freezing lead to dryer summers. In 2 x CO2 climate, thermal effects due to freeze/ thaw cycles are more pronounced and contribute to enhance the expected future overall winter warming. Impact of soil freezing on climate sensitivity is not uniform: the annual mean warming is amplified in North America (+ 15%) and Central Siberia (+ 36%) whereas it is reduced in Eastern Siberia ( – 23%). Nevertheless, all boreal lands undergo a strong attenuation of the warming during summertime. In agreement with some previous studies, these results indicate once more that soil freezing effects are significant on regional boreal climate. But this study also demonstrates its importance on regional boreal climate change and thus the necessity to include soil freezing in regional climate change predictions.

Abstract: The role of lakes and wetlands in present-day high latitude climate is quantified using a general circulation model of the atmosphere. The atmospheric model includes a lake module which is presented and validated. Seasonal and spatial wetland distribution can either be prescribed or calculated as a function of the hydrological budget of the wetlands themselves and of continental soil whose runoff feeds them. Simulated wetland extent is discussed both in simulations forced by observed climate and in general circulation model simulations. In off-line simulations, forced by ECMWF reanalyses, the lake model simulates correctly observed lake ice durations, while the wetland extent is somewhat underestimated in the boreal regions. Coupled to the general circulation model, the lake model yields satisfying ice durations, although the climate model biases have impacts on the modeled lake ice conditions. Boreal wetland extents are overestimated in the general circulation model as simulated precipitation is too high. The impact of inundated surfaces on the simulated climate is strongest in summer when these surfaces are ice-free. Wetlands seem to play a more important role than lakes in cooling the boreal regions in summer and in humidifying the atmosphere.

Abstract: For Central Greenland, water isotope analysis indicates a temperature difference of about 10degreesC since the Last Glacial Maximum (LGM). However, borehole thermometry and gas diffusion thermometry indicate that LGM surface temperatures were about 20degreesC colder than today. Two general circulation model studies have shown that changes in the seasonal precipitation timing in Central Greenland might have caused a warm bias in the LGM water isotope proxy temperatures, and that this bias could explain the difference in the estimated paleotemperatures. Here we present an analysis of a number of atmospheric general circulation model simulations mostly done within the framework of the Paleoclimate Modeling Intercomparison Project. The models suggest that the seasonal cycle of precipitation and surface mass balance over Central Greenland at the LGM might have been very different from today. This supports the idea that the accuracy of the water isotope thermometry at the LGM in Greenland might be compromised as a result of a modified surface mass balance seasonality. However, the models disagree on the amplitude and sign of the bias. For Central East Antarctica, a strong seasonality effect on the LGM isotopic signal is not simulated by any of the analyzed models. For the mid-Holocene (6 kyr BP) the models suggest relatively weak isotope paleothermometry biases linked to changes in the surface mass balance seasonality over both ice sheets. (C) 2003 Elsevier B.V. All rights reserved.

Abstract: Because many of the synoptic cyclones south of the 60degreesS parallel originate from 60degreesS and lower latitudes, nudging an atmospheric general circulation model (AGCM) with meteorological analyses at the periphery of the Antarctic region may be expected to exert a strong control on the atmospheric circulation inside the region. Here, the ECMWF reanalyses are used to nudge the atmospheric circulation of the Laboratoire de Meteorologie Dynamique Zoom (LMDZ) stretched-grid AGCM in a 15-yr simulation spanning the 1979-93 period. The horizontal resolution (grid spacing) in the model reaches similar to100 km south of 60degreesS. Nudging is exerted along the 60degreesS parallel, and this is called lateral nudging for the Antarctic region. Nudging is also performed farther north, near 50degrees and 40degreesS, but this is not essential for the results discussed here. Surface pressure and winds in the atmospheric column are nudged without relaxation to maximize control by the meteorological analyses, at the expense of some "noise'' confined to the latitudes where nudging is exerted. The performances of lateral nudging are evaluated with respect to station observations, the free (unnudged) model, the ECMWF reanalyses, and in limited instances with respect to nudging the surface pressure only. It is shown that the free model has limited but persistent surface pressure and geopotential defects in the Antarctic region, which are efficiently corrected by lateral nudging. Also, the laterally nudged simulations confirm, and to some extent correct, a geopotential deficiency of the ECMWF reanalyses over the east Antarctic continent previously identified by others. The monthly mean variability of surface climate at several stations along a coast-to-pole transect is analyzed. A significant fraction of the observed variability of surface pressure and temperature is reproduced. The fraction is often less than in the reanalyses. However, the differences are not large considering that the nudged model is forced at distances of hundreds to thousands of kilometers whereas the reanalyses are forced at much shorter distances, in principle right at each station site by the very station data. The variability of surface wind is significantly less well reproduced than that of pressure and temperature in both the nudged model and the reanalyses. Carefully adjusted polar physics in the LMDZ model seems to compensate for a distant observational constraint in the cases when the nudged model results appear similar or even superior to the reanalyses. Lateral nudging is less computationally intensive than global nudging, and it induces realistic variability and chronology while leaving full expression of the model physics in the region of interest. Laterally nudging an AGCM with meteorological analyses can offer complementary value over the analyses themselves, not only by producing additional atmospheric information not available from the analyses, but also by correcting possible regional defects in the analyses.

Abstract: This paper reports on analyses of the surface pressure, surface wind, and seasonality of temperature and precipitation in the central parts of Greenland and Antarctica as simulated by different general circulation models (GCMs) for the present and for the Last Glacial Maximum (LGM) climates. These parameters, in addition to the mean surface temperature, influence the air content of the ice either directly or through their influence on the pore volume. To correctly interpret the air content of the ice in terms of past ice sheet elevation changes, the variations of these parameters must therefore be known. Most of the simulations discussed here have been carried out within the framework of the Paleoclimate Modeling Intercomparison Project. Moreover, a stretched grid GCM has been used with a high rc solution over the ice sheets. We show that not taking into account changes of surface pressure at constant altitude between the LGM and today leads to an overestimation of past ice sheet elevation up to 150 m, while wind speed changes are too weak to have a significant influence on ice core air content. The results concerning changes of the amplitude or phase of the seasonal variations of precipitation and temperature are somewhat ambiguous. Most, but not; all, of the models suggest an intensification of the seasonal cycle of surface temperatures over central Greenland, and, to a lesser extent, over central East; Antarctica. Neglecting these changes might lead to an underestimation of past elevation by up to 140 m for the Greenland ice sheet, but this number is subject to large uncertainties.

Abstract: The intra-annual variability of Antarctic precipitation from the European Centre for Medium-range Weather Forecasts short-term meteorological forecasts and from climate simulations by the ARPEGE and LMD-Zoom general circulation models is presented and discussed. The spatial resolution of forecasts and simulations is high over the Antarctic region, about 100 km, so that the impact of topography and small-scale atmospheric dynamics are better resolved than with mc,re conventional model grids (about 300 km). All the models and forecasts show that the seasonality of precipitation is spatially very variable. Meridional coast-to-interior contrasts are marked, but equally strong variations are unexpectedly found where more homogeneity might be expected because of the homogeneity of the environment, e.g. on the high Antarctic plateau. Neither the forecasts nor the simulations confirm that precipitation is mostly maximum in winter over much of East Antarctica as suggested by scarce and potentially unreliable observations (Bromwich, 1988). Spring and fall maxima are quite frequent too, though summer maxima are rare. Daily precipitation statistics show more spatial pattern, with increasingly infrequent precipitation as distance from the coast toward the interior of the ice sheet increases. Several aspects of the intra-annual variability of precipitation can be interpreted in terms of atmospheric dynamics, but at both daily and seasonal time-scales the different forecasts and climate simulations often locally and regionally disagree with each other. Discrimination between models and their ability to reproduce the dynamics of Antarctic hydrology and progress on simulating such aspects of the Antarctic climate, is limited by the lack of reliable observation of precipitation variability.

Abstract: A high resolution GCM is used to examine the effect of changes in local surface climate parameters on the ice sheets that can influence the interpretation of the isotopic signal of the ice from deep cores. The model suggests that the 10 degrees C difference between the LGM surface temperature deduced from borehole thermometry and that deduced from the water isotope analysis to a great extent may be due to a modification of the precipitation seasonality in central Greenland. For central East Antarctica, the model tends to suggest a weak opposite bias.