Historical and idealized climate model experiments: an intercomparison of Earth system models of intermediate complexity

Both historical and idealized climate model experiments are performed with a variety of Earth system models of intermediate complexity (EMICs) as part of a community contribution to the Intergovernmental Panel on Climate Change Fifth Assessment Report. Historical simulations start at 850 CE and continue through to 2005. The standard simulations include changes in forcing from solar luminosity, Earth's orbital configuration, CO2, additional greenhouse gases, land use, and sulphate and volcanic aerosols. In spite of very different modelled pre-industrial global surface air temperatures, overall 20th century trends in surface air temperature and carbon uptake are reasonably well simulated when compared to observed trends. Land carbon fluxes show much more variation between models than ocean carbon fluxes, and recent land fluxes appear to be slightly underestimated. It is possible that recent modelled climate trends or climate–carbon feedbacks are overestimated resulting in too much land carbon loss or that carbon uptake due to CO2 and/or nitrogen fertilization is underestimated. Several one thousand year long, idealized, 2 × and 4 × CO2 experiments are used to quantify standard model characteristics, including transient and equilibrium climate sensitivities, and climate–carbon feedbacks. The values from EMICs generally fall within the range given by general circulation models. Seven additional historical simulations, each including a single specified forcing, are used to assess the contributions of different climate forcings to the overall climate and carbon cycle response. The response of surface air temperature is the linear sum of the individual forcings, while the carbon cycle response shows a non-linear interaction between land-use change and CO2 forcings for some models. Finally, the preindustrial portions of the last millennium simulations are used to assess historical model carbon-climate feedbacks. Given the specified forcing, there is a tendency for the EMICs to underestimate the drop in surface air temperature and CO2 between the Medieval Climate Anomaly and the Little Ice Age estimated from palaeoclimate reconstructions. This in turn could be a result of unforced variability within the climate system, uncertainty in the reconstructions of temperature and CO2, errors in the reconstructions of forcing used to drive the models, or the incomplete representation of certain processes within the models. Given the forcing datasets used in this study, the models calculate significant land-use emissions over the pre-industrial period. This implies that land-use emissions might need to be taken into account, when making estimates of climate–carbon feedbacks from palaeoclimate reconstructions.

Amplified inception of European Little Ice Age by sea ice-ocean-atmosphere feedbacks

The inception of the Little Ice Age (~1400–1700 AD) is believed to have been driven by an interplay of external forcing and climate system internal variability. While the hemispheric signal seems to have been dominated by solar irradiance and volcanic eruptions, the understanding of mechanisms shaping the climate on a continental scale is less robust. In an ensemble of transient model simulations and a new type of sensitivity experiments with artificial sea ice growth, the authors identify a sea ice–ocean–atmosphere feedback mechanism that amplifies the Little Ice Age cooling in the North Atlantic–European region and produces the temperature pattern suggested by paleoclimatic reconstructions. Initiated by increasing negative forcing, the Arctic sea ice substantially expands at the beginning of the Little Ice Age. The excess of sea ice is exported to the subpolar North Atlantic, where it melts, thereby weakening convection of the ocean. Consequently, northward ocean heat transport is reduced, reinforcing the expansion of the sea ice and the cooling of the Northern Hemisphere. In the Nordic Seas, sea surface height anomalies cause the oceanic recirculation to strengthen at the expense of the warm Barents Sea inflow, thereby further reinforcing sea ice growth. The absent ocean–atmosphere heat flux in the Barents Sea results in an amplified cooling over Northern Europe. The positive nature of this feedback mechanism enables sea ice to remain in an expanded state for decades up to a century, favoring sustained cold periods over Europe such as the Little Ice Age. Support for the feedback mechanism comes from recent proxy reconstructions around the Nordic Seas.

A centrifugal ice microtome for measurements of atmospheric CO₂ on air trapped in polar ice cores

For atmospheric CO2 reconstructions using ice cores, the technique to release the trapped air from the ice samples is essential for the precision and accuracy of the measurements. We present here a new dry extraction technique in combination with a new gas analytical system that together show significant improvements with respect to current systems. Ice samples (3–15 g) are pulverised using a novel centrifugal ice microtome (CIM) by shaving the ice in a cooled vacuum chamber (−27 °C) in which no friction occurs due to the use of magnetic bearings. Both, the shaving principle of the CIM and the use of magnetic bearings have not been applied so far in this field. Shaving the ice samples produces finer ice powder and releases a minimum of 90% of the trapped air compared to 50%–70% when needle crushing is employed. In addition, the friction-free motion with an optimized design to reduce contaminations of the inner surfaces of the device result in a reduced system offset of about 2.0 ppmv compared to 4.9 ppmv. The gas analytical part shows a higher precision than the corresponding part of our previous system by a factor of two, and all processes except the loading and cleaning of the CIM now run automatically. Compared to our previous system, the complete system shows a 3 times better measurement reproducibility of about 1.1 ppmv (1 σ) which is similar to the best reproducibility of other systems applied in this field. With this high reproducibility, no replicate measurements are required anymore for most future measurement campaigns resulting in a possible output of 12–20 measurements per day compared to a maximum of 6 with other systems.

Continuous Flow Analysis of labile iron in ice-cores

The important active and passive role of mineral dust aerosol in the climate and the global carbon cycle over the last glacial/interglacial cycles has been recognized. However, little data on the most important aeolian dust-derived biological micronutrient, iron (Fe), has so far been available from ice-cores from Greenland or Antarctica. Furthermore, Fe deposition reconstructions derived from the palaeoproxies particulate dust and calcium differ significantly from the Fe flux data available. The ability to measure high temporal resolution Fe data in polar ice-cores is crucial for the study of the timing and magnitude of relationships between geochemical events and biological responses in the open ocean. This work adapts an existing flow injection analysis (FIA) methodology for low-level trace Fe determinations with an existing glaciochemical analysis system, continuous flow analysis (CFA) of ice-cores. Fe-induced oxidation of N,N′-dimethyl-p-pheylenediamine (DPD) is used to quantify the biologically more important and easily leachable Fe fraction released in a controlled digestion step at pH ∼1.0. The developed method was successfully applied to the determination of labile Fe in ice-core samples collected from the Antarctic Byrd ice-core and the Greenland Ice-Core Project (GRIP) ice-core.

Ice-Bed Coupling Beneath and Beyond Ice Streams: Byrd Glacier, Antarctica

Ice sheet thickness is determined mainly by the strength of ice-bed coupling that controls holistic transitions from slow sheet flow to fast streamflow to buttressing shelf flow. Byrd Glacier has the largest ice drainage system in Antarctica and is the fastest ice stream entering Ross Ice Shelf. In 2004 two large subglacial lakes at the head of Byrd Glacier suddenly drained and increased the terminal ice velocity of Byrd Glacier from 820 m yr(-1) to 900 m yr(-1). This resulted in partial ice-bed recoupling above the lakes and partial decoupling along Byrd Glacier. An attempt to quantify this behavior is made using flowband and flowline models in which the controlling variable for ice height above the bed is the floating fraction phi of ice along the flowband and flowline. Changes in phi before and after drainage are obtained from available data, but more reliable data in the map plane are required before Byrd Glacier can be modeled adequately. A holistic sliding velocity is derived that depends on phi, with contributions from ice shearing over coupled beds and ice stretching over uncoupled beds, as is done in state-of-the-art sliding theories.

Local Rates of Ice-Sheet Thickness Change in Greenland

The rate of ice-sheet thickness change is calculated for 10 sites in Greenland by comparing measured values of ice vertical velocity and snow-accumulation rate. Vertical velocities are derived from repeat surveys of markers using precision global positioning system techniques, and accumulation rates are determined from stratigraphic analysis of firn cores. The results apply to time-scales covered by the firn-core records, which in most cases are a few decades. A spectrum of thickness-change rates is obtained, ranging from substantial thinning to slow thickening. The sites where ice-sheet thinning is indicated are located near the ice-sheet margin or in outlet glacier catchments. Interior and high-elevation sites are predominantly in balance or thickening slowly. Uncertainties in the rates of thickness change are dominated by errors in the determination of accumulation rates. The results of this work are broadly comparable with regional estimates of mass balance obtained from the analysis of catchment input vs discharge.

Climate Variability in West Antarctica Derived from Annual Accumulation-Rate Records from ITASE Firn/Ice Cores

Thirteen annually resolved accumulation-rate records covering the last similar to 200 years from the Pine Island-Thwaites and Ross drainage systems and the South Pole are used to examine climate variability over West Antarctica. Accumulation is controlled spatially by the topography of the ice sheet, and temporally by changes in moisture transport and cyclonic activity. A comparison of mean accumulation since 1970 at each site to the long-term mean indicates an increase in accumulation for sites located in the western sector of the Pine Island-Thwaites drainage system. Accumulation is negatively associated with the Southern Oscillation Index (Sol) for sites near the ice divide, and periods of sustained negative Sol (1940-42, 1991-95) correspond to above-mean accumulation at most sites. Correlations of the accumulation-rate records with sea-level pressure (SLP) and the SOI suggest that accumulation near the ice divide and in the Ross drainage system may be associated with the midlatitudes. The post-1970 increase in accumulation coupled with strong SLP-accumulation-rate correlations near the coast suggests recent intensification of cyclonic activity in the Pine Island-Thwaites drainage system.

An Ice-Core Proxy for Antarctic Circumpolar Zonal Wind Intensity

Using US National Centers for Environmental Prediction/US National Center for Atmospheric Research re-analysis data, we investigate the relationships between crustal ion (nssCa(2+)) concentrations from three West Antarctic ice cores, namely, Siple Dome (SD), ITASE00-1 (IT001) and ITASE01-5 (IT015), and primary components of the climate system, namely, air pressure/geopotential height, zonal (u) and meridional (v) wind strength. Linear correlation analyses between nssCa(2+) concentrations and both air-pressure and wind fields for the period of overlap between records indicate that the SD nssCa(2+) variation is positively correlated with spring circumpolar zonal wind, while IT001 nssCa(2+) has a positive correlation with circumpolar zonal wind throughout the year (r > 0.3, p < 0.01). Intensified Southern Westerlies circulation is conducive to transport of more crustal aerosols to both sites. Further correlation analyses between nssCa(2+) concentrations from SD and IT001 and atmospheric circulation suggest that the high inland plateau (represented by core IT001) is largely influenced by transport from the upper troposphere. IT015 nssCa(2+) is negatively correlated with westerly wind in October and November, suggesting that stronger westerly circulation may weaken the transport of crustal species to IT015. Correlations of nssCa(2+) from the three ice cores with the Antarctic Oscillation index are consistent with results developed from the wind-field investigation. In addition, calibration between nssCa(2+) concentration and the multivariate El Nino-Southern Oscillation (ENSO) index shows that crustal species transport to IT001 is enhanced during strong ENSO events.

Possible Effects of Anthropogenically-Increased CO2 on the Dynamics of Climate: Implications for Ice-Age Cycles

A dynamical model, developed to account for the observed major variations of global ice mass and atmospheric CO2 during the late Cenozoic, is used to provide a quantitative demonstration of the possibility that the anthropogenically-forced increase of atmospheric CO2, if maintained over a long period of time (perhaps by tectonic forcing), could displace the climatic system from an unstable regime of oscillating ice ages into a more stable regime representative of the pre-Pleistocene. This stable regime is characterized by orbitally-forced oscillations that are of much weaker amplitude than prevailed during the Pleistocene.

State of the Antarctic and Southern Ocean Climate System

This paper reviews developments in our understanding of the state of the Antarctic and Southern Ocean climate and its relation to the global climate system over the last few millennia. Climate over this and earlier periods has not been stable, as evidenced by the occurrence of abrupt changes in atmospheric circulation and temperature recorded in Antarctic ice core proxies for past climate. Two of the most prominent abrupt climate change events are characterized by intensification of the circumpolar westerlies (also known as the Southern Annular Mode) between similar to 6000 and 5000 years ago and since 1200-1000 years ago. Following the last of these is a period of major trans-Antarctic reorganization of atmospheric circulation and temperature between A. D. 1700 and 1850. The two earlier Antarctic abrupt climate change events appear linked to but predate by several centuries even more abrupt climate change in the North Atlantic, and the end of the more recent event is coincident with reorganization of atmospheric circulation in the North Pacific. Improved understanding of such events and of the associations between abrupt climate change events recorded in both hemispheres is critical to predicting the impact and timing of future abrupt climate change events potentially forced by anthropogenic changes in greenhouse gases and aerosols. Special attention is given to the climate of the past 200 years, which was recorded by a network of recently available shallow firn cores, and to that of the past 50 years, which was monitored by the continuous instrumental record. Significant regional climate changes have taken place in the Antarctic during the past 50 years. Atmospheric temperatures have increased markedly over the Antarctic Peninsula, linked to nearby ocean warming and intensification of the circumpolar westerlies. Glaciers are retreating on the peninsula, in Patagonia, on the sub-Antarctic islands, and in West Antarctica adjacent to the peninsula. The penetration of marine air masses has become more pronounced over parts of West Antarctica. Above the surface, the Antarctic troposphere has warmed during winter while the stratosphere has cooled year-round. The upper kilometer of the circumpolar Southern Ocean has warmed, Antarctic Bottom Water across a wide sector off East Antarctica has freshened, and the densest bottom water in the Weddell Sea has warmed. In contrast to these regional climate changes, over most of Antarctica, near-surface temperature and snowfall have not increased significantly during at least the past 50 years, and proxy data suggest that the atmospheric circulation over the interior has remained in a similar state for at least the past 200 years. Furthermore, the total sea ice cover around Antarctica has exhibited no significant overall change since reliable satellite monitoring began in the late 1970s, despite large but compensating regional changes. The inhomogeneity of Antarctic climate in space and time implies that recent Antarctic climate changes are due on the one hand to a combination of strong multidecadal variability and anthropogenic effects and, as demonstrated by the paleoclimate record, on the other hand to multidecadal to millennial scale and longer natural variability forced through changes in orbital insolation, greenhouse gases, solar variability, ice dynamics, and aerosols. Model projections suggest that over the 21st century the Antarctic interior will warm by 3.4 degrees +/- 1 degrees C, and sea ice extent will decrease by similar to 30%. Ice sheet models are not yet adequate enough to answer pressing questins about the effect of projected warming on mass balance and sea level. Considering the potentially major impacts of a warming climate on Antarctica, vigorous efforts are needed to better understand all aspects of the highly coupled Antarctic climate system as well as its influence on the Earth's climate and oceans.

Thermal Conditions at the Bed of the Laurentide Ice Sheet in Maine During Deglaciation: Implications for Esker Formation

The University of Maine Ice Sheet Model was used to study basal conditions during retreat of the Laurentide ice sheet in Maine. Within 150 km of the margin, basal melt rates average similar to 5 mm a(-1) during retreat. They decline over the next 100km, so areas of frozen bed develop in northern Maine during retreat. By integrating the melt rate over the drainage area typically subtended by an esker, we obtained a discharge at the margin of similar to 1.2 m(3) s(-1). While such a discharge could have moved the material in the Katahdin esker, it was likely too low to build the esker in the time available. Additional water from the glacier surface was required. Temperature gradients in the basal ice increase rapidly with distance from the margin. By conducting upward into the ice all of the additional viscous heat produced by any perturbation that increases the depth of flow in a flat conduit in a distributed drainage system, these gradients inhibit the formation of sharply arched conduits in which an esker can form. This may explain why eskers commonly seem to form near the margin and are typically segmented, with later segments overlapping onto earlier ones.

Solar Forcing Recorded by Aerosol Concentrations in Coastal Antarctic Glacier Ice, McMurdo Dry Valleys

Ice-core chemistry data from Victoria Lower Glacier, Antarctica, suggest, at least for the last 50 years, a direct influence of solar activity variations on the McMurdo Dry Valleys (MDV) climate system via controls on air-mass input from two competing environments: the East Antarctic ice sheet and the Ross Sea. During periods of increased solar activity, when total solar irradiance is relatively high, the MDV climate system appears to be dominated by air masses originating from the Ross Sea, leading to higher aerosol deposition. During reduced solar activity, the Antarctic interior seems to be the dominant air-mass source, leading to lower aerosol concentration in the ice-core record. We propose that the sensitivity of the MDV to variations in solar irradiance is caused by strong albedo differences between the ice-free MDV and the ice sheet.

Spatial and Temporal Melt Variability at Helheim Glacier, East Greenland, and Its Effect on Ice Dynamics

Understanding the behavior of large outlet glaciers draining the Greenland Ice Sheet is critical for assessing the impact of climate change on sea level rise. The flow of marine-terminating outlet glaciers is partly governed by calving-related processes taking place at the terminus but is also influenced by the drainage of surface runoff to the bed through moulins, cracks, and other pathways. To investigate the extent of the latter effect, we develop a distributed surface-energy-balance model for Helheim Glacier, East Greenland, to calculate surface melt and thereby estimate runoff. The model is driven by data from an automatic weather station operated on the glacier during the summers of 2007 and 2008, and calibrated with independent measurements of ablation. Modeled melt varies over the deployment period by as much as 68% relative to the mean, with melt rates approximately 77% higher on the lower reaches of the glacier trunk than on the upper glacier. We compare melt variations during the summer season to estimates of surface velocity derived from global positioning system surveys. Near the front of the glacier, there is a significant correlation (on >95% levels) between variations in runoff (estimated from surface melt) and variations in velocity, with a 1 day delay in velocity relative to melt. Although the velocity changes are small compared to accelerations previously observed following some calving events, our findings suggest that the flow speed of Helheim Glacier is sensitive to changes in runoff. The response is most significant in the heavily crevassed, fast-moving region near the calving front. The delay in the peak of the cross-correlation function implies a transit time of 12-36 h for surface runoff to reach the bed.

Controls on the Basal Water Pressure in Subglacial Channels Near the Margin of the Greenland Ice Sheet

Assuming a channelized drainage system in steady state, we investigate the influence of enhanced surface melting on the water pressure in subglacial channels, compared to that of changes in conduit geometry, ice rheology and catchment variations. The analysis is carried out for a specific part of the western Greenland ice-sheet margin between 66 degrees N and 66 degrees 30' N using new high-resolution digital elevation models of the subglacial topography and the ice-sheet surface, based on an airborne ice-penetrating radar survey in 2003 and satellite repeat-track interferometric synthetic aperture radar analysis of European Remote-sensing Satellite 1 and 2 (ERS-1/-2) imagery, respectively. The water pressure is calculated up-glacier along a likely subglacial channel at distances of 1, 5 and 9 km from the outlet at the ice margin, using a modified version of Rothlisberger's equation. Our results show that for the margin of the western Greenland ice sheet, the water pressure in subglacial channels is not sensitive to realistic variations in catchment size and mean surface water input compared to small changes in conduit geometry and ice rheology.

Fire in ice: two millennia of boreal forest fire history from the Greenland NEEM ice core

Biomass burning is a major source of greenhouse gases and influences regional to global climate. Pre-industrial fire-history records from black carbon, charcoal and other proxies provide baseline estimates of biomass burning at local to global scales spanning millennia, and are thus useful to examine the role of fire in the carbon cycle and climate system. Here we use the specific biomarker levoglucosan together with black carbon and ammonium concentrations from the North Greenland Eemian (NEEM) ice cores (77.49° N, 51.2° W; 2480 m a.s.l) over the past 2000 years to infer changes in boreal fire activity. Increases in boreal fire activity over the periods 1000–1300 CE and decreases during 700–900 CE coincide with high-latitude NH temperature changes. Levoglucosan concentrations in the NEEM ice cores peak between 1500 and 1700 CE, and most levoglucosan spikes coincide with the most extensive central and northern Asian droughts of the past millennium. Many of these multi-annual droughts are caused by Asian monsoon failures, thus suggesting a connection between low- and high-latitude climate processes. North America is a primary source of biomass burning aerosols due to its relative proximity to the Greenland Ice Cap. During major fire events, however, isotopic analyses of dust, back trajectories and links with levoglucosan peaks and regional drought reconstructions suggest that Siberia is also an important source of pyrogenic aerosols to Greenland.