Polar MM5 Simulations of the Winter Climate of the Laurentide Ice Sheet at the LGM

Optimized regional climate simulations are conducted using the Polar MM5, a version of the fifth-generation Pennsylvania State University-NCAR Mesoscale Model (MM5), with a 60-km horizontal resolution domain over North America during the Last Glacial Maximum (LGM, 21 000 calendar years ago), when much of the continent was covered by the Laurentide Ice Sheet (LIS). The objective is to describe the LGM annual cycle at high spatial resolution with an emphasis on the winter atmospheric circulation. Output from a tailored NCAR Community Climate Model version 3 (CCM3) simulation of the LGM climate is used to provide the initial and lateral boundary conditions for Polar MM5. LGM boundary conditions include continental ice sheets, appropriate orbital forcing, reduced CO2 concentration, paleovegetation, modified sea surface temperatures, and lowered sea level. Polar MM5 produces a substantially different atmospheric response to the LGM boundary conditions than CCM3 and other recent GCM simulations. In particular, from November to April the upper-level flow is split around a blocking anticyclone over the LIS, with a northern branch over the Canadian Arctic and a southern branch impacting southern North America. The split flow pattern is most pronounced in January and transitions into a single, consolidated jet stream that migrates northward over the LIS during summer. Sensitivity experiments indicate that the winter split flow in Polar MM5 is primarily due to mechanical forcing by LIS, although model physics and resolution also contribute to the simulated flow configuration. Polar MM5 LGM results are generally consistent with proxy climate estimates in the western United States, Alaska, and the Canadian Arctic and may help resolve some long-standing discrepancies between proxy data and previous simulations of the LGM climate.

Denitrification in the Hypolimnion of Permanently Ice-Covered Lake Bonney, Antarctica

The distribution of denitrification was investigated in the hypolimnion of the east and west lobes of permanently ice-covered Lake Bonney, Taylor Valley, Antarctica. Anomalously high concentrations of dissolved inorganic nitrogen (DIN; nitrate, nitrite, ammonium and nitrous oxide) in the oxygen-depleted hypolimnion of the east lobe of the Lake implied that denitrification is or was active in the west, but not in the east lobe. While previous investigations reported no detectable denitrification in the east lobe, we measured active denitrification in samples from both the east and west lobes. In the west lobe, measured denitrification rates exhibited a maximum at the depth of the chemocline and denitrification was not detectable in either the oxic surface waters or in the deep water where nitrate was absent. In the east lobe, denitrification was detected below the chemocline, at the depths where ammonium, nitrate, nitrite and nitrous oxide are all present at anomalously high levels, Trace metal availability was manipulated in incubation experiments in order to determine whether trace metal toxicity in the east lobe could explain the difference in nitrogen cycling between the 2 lobes. There were no consistent stimulatory effects of metal chelators or nutrient addition on the rate of denitrification in either lobe, so that the mechanisms underlying the unusual N cycle of the east lobe remain unknown. We conclude that all the ingredients necessary to allow denitrification to occur are present in the east lobe. However, even though denitrification could be detected under certain conditions in incubations, denitrification is inhibited under the in situ conditions of the lake.

Spatial and temporal patterns of greenness on the Yamal Peninsula, Russia: interactions of ecological and social factors affecting the Arctic normalized difference vegetation index

The causes of a greening trend detected in the Arctic using the normalized difference vegetation index (NDVI) are still poorly understood. Changes in NDVI are a result of multiple ecological and social factors that affect tundra net primary productivity. Here we use a 25 year time series of AVHRR-derived NDVI data (AVHRR: advanced very high resolution radiometer), climate analysis, a global geographic information database and ground-based studies to examine the spatial and temporal patterns of vegetation greenness on the Yamal Peninsula, Russia. We assess the effects of climate change, gas-field development, reindeer grazing and permafrost degradation. In contrast to the case for Arctic North America, there has not been a significant trend in summer temperature or NDVI, and much of the pattern of NDVI in this region is due to disturbances. There has been a 37% change in early-summer coastal sea-ice concentration, a 4% increase in summer land temperatures and a 7% change in the average time-integrated NDVI over the length of the satellite observations. Gas-field infrastructure is not currently extensive enough to affect regional NDVI patterns. The effect of reindeer is difficult to quantitatively assess because of the lack of control areas where reindeer are excluded. Many of the greenest landscapes on the Yamal are associated with landslides and drainage networks that have resulted from ongoing rapid permafrost degradation. A warming climate and enhanced winter snow are likely to exacerbate positive feedbacks between climate and permafrost thawing. We present a diagram that summarizes the social and ecological factors that influence Arctic NDVI. The NDVI should be viewed as a powerful monitoring tool that integrates the cumulative effect of a multitude of factors affecting Arctic land-cover change.

Long-term perspective on wildfires in the western USA

Understanding the causes and consequences of wildfires in forests of the western United States requires integrated information about fire, climate changes, and human activity on multiple temporal scales. We use sedimentary charcoal accumulation rates to construct long-term variations in fire during the past 3,000 y in the American West and compare this record to independent fire-history data from historical records and fire scars. There has been a slight decline in burning over the past 3,000 y, with the lowest levels attained during the 20th century and during the Little Ice Age (LIA, ca. 1400-1700 CE Common Era]). Prominent peaks in forest fires occurred during the Medieval Climate Anomaly (ca. 950-1250 CE) and during the 1800s. Analysis of climate reconstructions beginning from 500 CE and population data show that temperature and drought predict changes in biomass burning up to the late 1800s CE. Since the late 1800s, human activities and the ecological effects of recent high fire activity caused a large, abrupt decline in burning similar to the LIA fire decline. Consequently, there is now a forest ``fire deficit'' in the western United States attributable to the combined effects of human activities, ecological, and climate changes. Large fires in the late 20th and 21st century fires have begun to address the fire deficit, but it is continuing to grow.

Improving accuracy and precision of ice core δD(CH4) analyses using methane pre-pyrolysis and hydrogen post-pyrolysis trapping and subsequent chromatographic separation

Firn and polar ice cores offer the only direct palaeoatmospheric archive. Analyses of past greenhouse gas concentrations and their isotopic compositions in air bubbles in the ice can help to constrain changes in global biogeochemical cycles in the past. For the analysis of the hydrogen isotopic composition of methane (δD(CH4) or δ2H(CH4)) 0.5 to 1.5 kg of ice was hitherto used. Here we present a method to improve precision and reduce the sample amount for δD(CH4) measurements in (ice core) air. Pre-concentrated methane is focused in front of a high temperature oven (pre-pyrolysis trapping), and molecular hydrogen formed by pyrolysis is trapped afterwards (post-pyrolysis trapping), both on a carbon-PLOT capillary at −196 °C. Argon, oxygen, nitrogen, carbon monoxide, unpyrolysed methane and krypton are trapped together with H2 and must be separated using a second short, cooled chromatographic column to ensure accurate results. Pre- and post-pyrolysis trapping largely removes the isotopic fractionation induced during chromatographic separation and results in a narrow peak in the mass spectrometer. Air standards can be measured with a precision better than 1‰. For polar ice samples from glacial periods, we estimate a precision of 2.3‰ for 350 g of ice (or roughly 30 mL – at standard temperature and pressure (STP) – of air) with 350 ppb of methane. This corresponds to recent tropospheric air samples (about 1900 ppb CH4) of about 6 mL (STP) or about 500 pmol of pure CH4.

Influence of microphysical schemes on atmospheric water in the Weather Research and Forecasting model

This study examines how different microphysical parameterization schemes influence orographically induced precipitation and the distributions of hydrometeors and water vapour for midlatitude summer conditions in the Weather Research and Forecasting (WRF) model. A high-resolution two-dimensional idealized simulation is used to assess the differences between the schemes in which a moist air flow is interacting with a bell-shaped 2 km high mountain. Periodic lateral boundary conditions are chosen to recirculate atmospheric water in the domain. It is found that the 13 selected microphysical schemes conserve the water in the model domain. The gain or loss of water is less than 0.81% over a simulation time interval of 61 days. The differences of the microphysical schemes in terms of the distributions of water vapour, hydrometeors and accumulated precipitation are presented and discussed. The Kessler scheme, the only scheme without ice-phase processes, shows final values of cloud liquid water 14 times greater than the other schemes. The differences among the other schemes are not as extreme, but still they differ up to 79% in water vapour, up to 10 times in hydrometeors and up to 64% in accumulated precipitation at the end of the simulation. The microphysical schemes also differ in the surface evaporation rate. The WRF single-moment 3-class scheme has the highest surface evaporation rate compensated by the highest precipitation rate. The different distributions of hydrometeors and water vapour of the microphysical schemes induce differences up to 49 W m−2 in the downwelling shortwave radiation and up to 33 W m−2 in the downwelling longwave radiation.

Temperature reconstruction from 10 to 120 kyr b2k from the NGRIP ice core

In order to reconstruct the temperature of the North Greenland Ice Core Project (NGRIP) site, new measurements of δ15N have been performed covering the time period from the beginning of the Holocene to Dansgaard–Oeschger (DO) event 8. Together with previously measured and mostly published δ15N data, we present for the first time a NGRIP temperature reconstruction for the whole last glacial period from 10 to 120 kyr b2k (thousand years before 2000 AD) including every DO event based on δ15N isotope measurements combined with a firn densification and heat diffusion model. The detected temperature rises at the onset of DO events range from 5 °C (DO 25) up to 16.5 °C (DO 11) with an uncertainty of ±3 °C. To bring measured and modelled data into agreement, we had to reduce the accumulation rate given by the NGRIP ss09sea06bm timescale in some periods by 30 to 35%, especially during the last glacial maximum. A comparison between reconstructed temperature and δ18Oice data confirms that the isotopic composition of the stadial was strongly influenced by seasonality. We evidence an anticorrelation between the variations of the δ18Oice sensitivity to temperature (referred to as α) and obliquity in agreement with a simple Rayleigh distillation model. Finally, we suggest that α might be influenced by the Northern Hemisphere ice sheet volume.

Temperature and precipitation signal in two Alpine ice cores over the period 1961-2001

Water stable isotope ratios and net snow accumulation in ice cores are commonly interpreted as temperature or precipitation proxies. However, only in a few cases has a direct calibration with instrumental data been attempted. In this study we took advantage of the dense network of observations in the European Alpine region to rigorously test the relationship of the annual and seasonal resolved proxy data from two highly resolved ice cores with local temperature and precipitation. We focused on the time period 1961–2001 with the highest amount and quality of meteorological data and the minimal uncertainty in ice core dating (±1 year). The two ice cores were retrieved from the Fiescherhorn glacier (northern Alps, 3900 m a.s.l.), and Grenzgletscher (southern Alps, 4200 m a.s.l.). A parallel core from the Fiescherhorn glacier allowed assessing the reproducibility of the ice core proxy data. Due to the orographic barrier, the two flanks of the Alpine chain are affected by distinct patterns of precipitation. The different location of the two glaciers therefore offers a unique opportunity to test whether such a specific setting is reflected in the proxy data. On a seasonal scale a high fraction of δ18O variability was explained by the seasonal cycle of temperature (~60% for the ice cores, ~70% for the nearby stations of the Global Network of Isotopes in Precipitation – GNIP). When the seasonality is removed, the correlations decrease for all sites, indicating that factors other than temperature such as changing moisture sources and/or precipitation regimes affect the isotopic signal on this timescale. Post-depositional phenomena may additionally modify the ice core data. On an annual scale, the δ18O/temperature relationship was significant at the Fiescherhorn, whereas for Grenzgletscher this was the case only when weighting the temperature with precipitation. In both cases the fraction of interannual temperature variability explained was ~20%, comparable to the values obtained from the GNIP stations data. Consistently with previous studies, we found an altitude effect for the δ18O of −0.17‰/100 m for an extended elevation range combining data of the two ice core sites and four GNIP stations. Significant correlations between net accumulation and precipitation were observed for Grenzgletscher during the entire period of investigation, whereas for Fiescherhorn this was the case only for the less recent period (1961–1977). Local phenomena, probably related to wind, seem to partly disturb the Fiescherhorn accumulation record. Spatial correlation analysis shows the two glaciers to be influenced by different precipitation regimes, with the Grenzgletscher reflecting the characteristic precipitation regime south of the Alps and the Fiescherhorn accumulation showing a pattern more closely linked to northern Alpine stations.

A new sensitive method for the quantification of glyoxal and methylglyoxal in snow and ice by stir bar sorptive extraction and liquid desorption-HPLC-ESI-MS

In this study, the development of a new sensitive method for the analysis of alpha-dicarbonyls glyoxal (G) and methylglyoxal (MG) in environmental ice and snow is presented. Stir bar sorptive extraction with in situ derivatization and liquid desorption (SBSE-LD) was used for sample extraction, enrichment, and derivatization. Measurements were carried out using high-performance liquid chromatography coupled to electrospray ionization tandem mass spectrometry (HPLC-ESI-MS/MS). As part of the method development, SBSE-LD parameters such as extraction time, derivatization reagent, desorption time and solvent, and the effect of NaCl addition on the SBSE efficiency as well as measurement parameters of HPLC-ESI-MS/MS were evaluated. Calibration was performed in the range of 1–60 ng/mL using spiked ultrapure water samples, thus incorporating the complete SBSE and derivatization process. 4-Fluorobenzaldehyde was applied as internal standard. Inter-batch precision was <12 % RSD. Recoveries were determined by means of spiked snow samples and were 78.9 ± 5.6 % for G and 82.7 ± 7.5 % for MG, respectively. Instrumental detection limits of 0.242 and 0.213 ng/mL for G and MG were achieved using the multiple reaction monitoring mode. Relative detection limits referred to a sample volume of 15 mL were 0.016 ng/mL for G and 0.014 ng/mL for MG. The optimized method was applied for the analysis of snow samples from Mount Hohenpeissenberg (close to the Meteorological Observatory Hohenpeissenberg, Germany) and samples from an ice core from Upper Grenzgletscher (Monte Rosa massif, Switzerland). Resulting concentrations were 0.085–16.3 ng/mL for G and 0.126–3.6 ng/mL for MG. Concentrations of G and MG in snow were 1–2 orders of magnitude higher than in ice core samples. The described method represents a simple, green, and sensitive analytical approach to measure G and MG in aqueous environmental samples.

Eemian interglacial reconstructed from a Greenland folded ice core

Efforts to extract a Greenland ice core with a complete record of the Eemian interglacial (130,000 to 115,000 years ago) have until now been unsuccessful. The response of the Greenland ice sheet to the warmer-than-present climate of the Eemian has thus remained unclear. Here we present the new North Greenland Eemian Ice Drilling ('NEEM') ice core and show only a modest ice-sheet response to the strong warming in the early Eemian. We reconstructed the Eemian record from folded ice using globally homogeneous parameters known from dated Greenland and Antarctic ice-core records. On the basis of water stable isotopes, NEEM surface temperatures after the onset of the Eemian (126,000 years ago) peaked at 8±4 degrees Celsius above the mean of the past millennium, followed by a gradual cooling that was probably driven by the decreasing summer insolation. Between 128,000 and 122,000 years ago, the thickness of the northwest Greenland ice sheet decreased by 400±250 metres, reaching surface elevations 122,000 years ago of 130±300 metres lower than the present. Extensive surface melt occurred at the NEEM site during the Eemian, a phenomenon witnessed when melt layers formed again at NEEM during the exceptional heat of July 2012. With additional warming, surface melt might become more common in the future.

Terrestrial and submarine evidence for the extent and timing of the Last Glacial Maximum and the onset of deglaciation on the maritime-Antarctic and sub-Antarctic islands

This paper is the maritime and sub–Antarctic contribution to the Scientific Committee for Antarctic Research (SCAR) Past Antarctic Ice Sheet Dynamics (PAIS) community Antarctic Ice Sheet reconstruction. The overarching aim for all sectors of Antarctica was to reconstruct the Last Glacial Maximum (LGM) ice sheet extent and thickness, and map the subsequent deglaciation in a series of 5000 year time slices. However, our review of the literature found surprisingly few high quality chronological constraints on changing glacier extents on these timescales in the maritime and sub–Antarctic sector. Therefore, in this paper we focus on an assessment of the terrestrial and offshore evidence for the LGM ice extent, establishing minimum ages for the onset of deglaciation, and separating evidence of deglaciation from LGM limits from those associated with later Holocene glacier fluctuations. Evidence included geomorphological descriptions of glacial landscapes, radiocarbon dated basal peat and lake sediment deposits, cosmogenic isotope ages of glacial features and molecular biological data. We propose a classification of the glacial history of the maritime and sub–Antarctic islands based on this assembled evidence. These include: (Type I) islands which accumulated little or no LGM ice; (Type II) islands with a limited LGM ice extent but evidence of extensive earlier continental shelf glaciations; (Type III) seamounts and volcanoes unlikely to have accumulated significant LGM ice cover; (Type IV) islands on shallow shelves with both terrestrial and submarine evidence of LGM (and/or earlier) ice expansion; (Type V) Islands north of the Antarctic Polar Front with terrestrial evidence of LGM ice expansion; and (Type VI) islands with no data. Finally, we review the climatological and geomorphological settings that separate the glaciological history of the islands within this classification scheme.

Peak Last Glacial weathering intensity on the North American continent recorded by the authigenic Hf isotope composition of North Atlantic deep-sea sediments

We have retrieved radiogenic hafnium (Hf) isotope compositions (ɛHf) from authigenic Fe–Mn oxyhydroxides of deep northwest Atlantic sediments deposited over the past 26 ka to investigate the oceanic evidence of changes in dissolved weathering inputs from NE America during the last deglaciation. The extraction of seawater-derived Hf isotopic compositions from Fe–Mn oxyhydroxides is not a standard procedure. Comparisons between the Al/Hf ratios and Hf isotopic compositions of the chemically extracted authigenic phase on the one hand, and those of the corresponding detrital fractions on the other, provide evidence that the composition of past seawater has been reliably obtained for most sampled depths with our leaching procedures. This is endorsed most strongly by data for a sediment core from 4250 m water depth at the deeper Blake Ridge, for which consistent replicates were produced throughout. The Hf isotopic composition of the most recent sample in this core also closely matches that of nearby present day central North Atlantic seawater. Comparison with previously published seawater Nd and Pb isotope compositions obtained on the same cores shows that both Hf and Pb were released incongruently during incipient chemical weathering, but responded differently to the deglacial retreat of the Laurentide Ice Sheet. Hafnium was released more congruently during peak glacial conditions of the Last Glacial Maximum (LGM) and changed to typical incongruent interglacial ɛHf signatures either during or shortly after the LGM. This indicates that some zircon-derived Hf was released to seawater during the LGM. Conversely, there is no clear evidence for an increase in the influence of weathering of Lu-rich mineral phases during deglaciation, possibly since relatively unradiogenic Hf contributions from feldspar weathering were superimposed. While the authigenic Pb isotope signal in the same marine sediment samples traced peak chemical weathering rates on continental North America during the transition to the Holocene a similar incongruent excursion is notably absent in the Hf isotope record. The early change towards more radiogenic ɛHf in relation to the LGM may provide direct evidence for the transition from a cold-based to a warm-based Laurentide Ice Sheet on the Atlantic sector of North America.

Isotopic constraints on marine and terrestrial N2O emissions during the last deglaciation

Nitrous oxide (N2O) is an important greenhouse gas and ozone-depleting substance that has anthropogenic as well as natural marine and terrestrial sources. The tropospheric N2O concentrations have varied substantially in the past in concert with changing climate on glacial–interglacial and millennial timescales. It is not well understood, however, how N2O emissions from marine and terrestrial sources change in response to varying environmental conditions. The distinct isotopic compositions of marine and terrestrial N2O sources can help disentangle the relative changes in marine and terrestrial N2O emissions during past climate variations. Here we present N2O concentration and isotopic data for the last deglaciation, from 16,000 to 10,000 years before present, retrieved from air bubbles trapped in polar ice at Taylor Glacier, Antarctica. With the help of our data and a box model of the N2O cycle, we find a 30 per cent increase in total N2O emissions from the late glacial to the interglacial, with terrestrial and marine emissions contributing equally to the overall increase and generally evolving in parallel over the last deglaciation, even though there is no a priori connection between the drivers of the two sources. However, we find that terrestrial emissions dominated on centennial timescales, consistent with a state-of-the-art dynamic global vegetation and land surface process model that suggests that during the last deglaciation emission changes were strongly influenced by temperature and precipitation patterns over land surfaces. The results improve our understanding of the drivers of natural N2O emissions and are consistent with the idea that natural N2O emissions will probably increase in response to anthropogenic warming.

A community-based geological reconstruction of Antarctic Ice Sheet deglaciation since the Last Glacial Maximum

A robust understanding of Antarctic Ice Sheet deglacial history since the Last Glacial Maximum is important in order to constrain ice sheet and glacial-isostatic adjustment models, and to explore the forcing mechanisms responsible for ice sheet retreat. Such understanding can be derived from a broad range of geological and glaciological datasets and recent decades have seen an upsurge in such data gathering around the continent and Sub-Antarctic islands. Here, we report a new synthesis of those datasets, based on an accompanying series of reviews of the geological data, organised by sector. We present a series of timeslice maps for 20 ka, 15 ka, 10 ka and 5 ka, including grounding line position and ice sheet thickness changes, along with a clear assessment of levels of confidence. The reconstruction shows that the Antarctic Ice sheet did not everywhere reach the continental shelf edge at its maximum, that initial retreat was asynchronous, and that the spatial pattern of deglaciation was highly variable, particularly on the inner shelf. The deglacial reconstruction is consistent with a moderate overall excess ice volume and with a relatively small Antarctic contribution to meltwater pulse 1a. We discuss key areas of uncertainty both around the continent and by time interval, and we highlight potential priorities for future work. The synthesis is intended to be a resource for the modelling and glacial geological community.

Evaluation of aeolian dust records obtained from Polar Ice Cores

When an ice core sample is analysed for its aeolian dust content, it is melted and the particles detected are suspended in water. Consequently, dust measurement techniques employed in the ice core community differ from those used for in-situ studies of airborne dust. Methods commonly used to classify insolubles suspended in a liquid are either based on the particles’ interaction with light or on the detection of resistive pulses by means of Coulter counting. Data sets obtained with Coulter counters are widely accepted as references and other techniques are judged against their ability to reproduce these. Unfortunately, optically acquired ice core dust records were found to differ. By analyzing two sections of the NEEM dust record, two different evaluation procedures are discussed before a third protocol is proposed. It is found that relative changes in the archived dust load can be reproduced, while the simultaneous attainment of absolute concentrations or changes in the grain size frequency histograms in high resolution remains difficult.