Dust sources and deposition during the last glacial maximum and current climate: a comparison of model results with paleodata from ice cores and marine sediments

Mineral dust aerosols in the atmosphere have the potential to affect the global climate by influencing the radiative balance of the atmosphere and the supply of micronutrients to the ocean. Ice and marine sediment cores indicate that dust deposition from the atmosphere was at some locations 2–20 times greater during glacial periods, raising the possibility that mineral aerosols might have contributed to climate change on glacial-interglacial time scales. To address this question, we have used linked terrestrial biosphere, dust source, and atmospheric transport models to simulate the dust cycle in the atmosphere for current and last glacial maximum (LGM) climates. We obtain a 2.5-fold higher dust loading in the entire atmosphere and a twenty-fold higher loading in high latitudes, in LGM relative to present. Comparisons to a compilation of atmospheric dust deposition flux estimates for LGM and present in marine sediment and ice cores show that the simulated flux ratios are broadly in agreement with observations; differences suggest where further improvements in the simple dust model could be made. The simulated increase in high-latitude dustiness depends on the expansion of unvegetated areas, especially in the high latitudes and in central Asia, caused by a combination of increased aridity and low atmospheric [CO2]. The existence of these dust source areas at the LGM is supported by pollen data and loess distribution in the northern continents. These results point to a role for vegetation feedbacks, including climate effects and physiological effects of low [CO2], in modulating the atmospheric distribution of dust.

The consolidation and bond strength of rafted sea ice

The consolidation and bond strength of rafted sea ice were investigated through a series of experiments undertaken in the Ice Physics Laboratory at the UCL. To simulate a section of rafted sea ice, blocks of laboratory grown saline ice were stacked in an insulated tank with spacers between adjacent blocks to allow saline water to flood in. The rate of consolidation was then monitored using a combination of temperature readings recorded in the ice and liquid layer, salinity measurements of the liquid layer, and cores taken at specific times of interest. Two states of consolidation were observed: thermodynamic consolidation where the ice blocks were physically bonded but the bond strength was weak, and mechanical consolidation where the bond had reached full strength. Results showed that the rafted ice had physically bonded in less than a day, however it took many more days (6 to 30 depending on the environmental conditions) for the bond to reach maximum strength. Increasing the thickness of the ice, the salinity of the water and the inter-block gap size all increased the consolidation time. Once consolidated, ice cores were taken and sheared using the asymmetric four-point bending method to measure the strength of the bond between the ice blocks. These were then compared to the shear strength of solid ice blocks simulating level sea ice. Our results show that the shear strength of the bond between the rafted ice blocks is about 30% weaker than that of level ice.

Large inert carbon pool in the terrestrial biosphere at the Last Glacial Maximum

During each of the late Pleistocene glacial–interglacial transitions, atmospheric carbon dioxide concentrations rose by almost 100 ppm. The sources of this carbon are unclear, and efforts to identify them are hampered by uncertainties in the magnitude of carbon reservoirs and fluxes under glacial conditions. Here we use oxygen isotope measurements from air trapped in ice cores and ocean carbon-cycle modelling to estimate terrestrial and oceanic gross primary productivity during the Last Glacial Maximum. We find that the rate of gross terrestrial primary production during the Last Glacial Maximum was about 40±10 Pg C yr−1, half that of the pre-industrial Holocene. Despite the low levels of photosynthesis, we estimate that the late glacial terrestrial biosphere contained only 330 Pg less carbon than pre-industrial levels. We infer that the area covered by carbon-rich but unproductive biomes such as tundra and cold steppes was significantly larger during the Last Glacial Maximum, consistent with palaeoecological data. Our data also indicate the presence of an inert carbon pool of 2,300 Pg C, about 700 Pg larger than the inert carbon locked in permafrost today. We suggest that the disappearance of this carbon pool at the end of the Last Glacial Maximum may have contributed to the deglacial rise in atmospheric carbon dioxide concentrations.

Late Quaternary fire regimes of Australasia

We have compiled 223 sedimentary charcoal records from Australasia in order to examine the temporal and spatial variability of fire regimes during the Late Quaternary. While some of these records cover more than a full glacial cycle, here we focus on the last 70,000 years when the number of individual records in the compilation allows more robust conclusions. On orbital time scales, fire in Australasia predominantly reflects climate, with colder periods characterized by less and warmer intervals by more biomass burning. The composite record for the region also shows considerable millennial-scale variability during the last glacial interval (73.5–14.7 ka). Within the limits of the dating uncertainties of individual records, the variability shown by the composite charcoal record is more similar to the form, number and timing of Dansgaard–Oeschger cycles as observed in Greenland ice cores than to the variability expressed in the Antarctic ice-core record. The composite charcoal record suggests increased biomass burning in the Australasian region during Greenland Interstadials and reduced burning during Greenland Stadials. Millennial-scale variability is characteristic of the composite record of the sub-tropical high pressure belt during the past 21 ka, but the tropics show a somewhat simpler pattern of variability with major peaks in biomass burning around 15 ka and 8 ka. There is no distinct change in fire regime corresponding to the arrival of humans in Australia at 50 ± 10 ka and no correlation between archaeological evidence of increased human activity during the past 40 ka and the history of biomass burning. However, changes in biomass burning in the last 200 years may have been exacerbated or influenced by humans.

Systematic study of the impact of fresh water fluxes on the glacial carbon cycle

During glacial periods, atmospheric CO2 concentration increases and decreases by around 15 ppm. At the same time, the climate changes gradually in Antarctica. Such climate changes can be simulated in models when the AMOC (Atlantic Meridional Oceanic Circulation) is weakened by adding fresh water to the North Atlantic. The impact on the carbon cycle is less straightforward, and previous studies give opposite results. Because the models and the fresh water fluxes were different in these studies, it prevents any direct comparison and hinders finding whether the discrepancies arise from using different models or different fresh water fluxes. In this study we use the CLIMBER-2 coupled climate carbon model to explore the impact of different fresh water fluxes. In both preindustrial and glacial states, the addition of fresh water and the resulting slow-down of the AMOC lead to an uptake of carbon by the ocean and a release by the terrestrial biosphere. The duration, shape and amplitude of the fresh water flux all have an impact on the change of atmospheric CO2 because they modulate the change of the AMOC. The maximum CO2 change linearly depends on the time integral of the AMOC change. The different duration, amplitude, and shape of the fresh water flux cannot explain the opposite evolution of ocean and vegetation carbon inventory in different models. The different CO2 evolution thus depends on the AMOC response to the addition of fresh water and the resulting climatic change, which are both model dependent. In CLIMBER-2, the rise of CO2 recorded in ice cores during abrupt events can be simulated under glacial conditions, especially when the sinking of brines in the Southern Ocean is taken into account. The addition of fresh water in the Southern Hemisphere leads to a decline of CO2, contrary to the addition of fresh water in the Northern Hemisphere.

Global patterns of vegetation response to millennial-scale variability and rapid climate change during the last glacial period

Ninety-four sites worldwide have sufficient resolution and dating to document the impact of millennial-scale climate variability on vegetation and fire regimes during the last glacial period. Although Dansgaard–Oeschger (D–O) cycles all show a basically similar gross structure, they vary in the magnitude and the length of the warm and cool intervals. We illustrate the geographic patterns in the climate-induced changes in vegetation by comparing D–O 6, D–O 8 and D–O 19. There is a strong response to both D–O warming events and subsequent cooling, most marked in the northern extratropics. Pollen records from marine cores from the northern extratropics confirm that there is no lag between the change in climate and the vegetation response, within the limits of the dating resolution (50–100 years). However, the magnitude of the change in vegetation is regionally specific and is not a simple function of either the magnitude or the duration of the change in climate as registered in Greenland ice cores. Fire regimes also show an initial immediate response to climate changes, but during cooling intervals there is a slow recovery of biomass burning after the initial reduction, suggesting a secondary control through the recovery of vegetation productivity. In the extratropics, vegetation changes are largely determined by winter temperatures while in the tropics they are largely determined by changes in plant-available water. Tropical vegetation records show changes corresponding to Heinrich Stadials but the response to D–O warming events is less marked than in the northern extratropics. There are very few high-resolution records from the Southern Hemisphere extratropics, but these records also show both a vegetation and fire response to millennial-scale climate variability. It is not yet possible to determine unequivocally whether terrestrial records reflect the asynchroneity apparent in the ice-core records.

DIRTMAP: the geological record of dust.

Atmospheric dust is an important feedback in the climate system, potentially affecting the radiative balance and chemical composition of the atmosphere and providing nutrients to terrestrial and marine ecosystems. Yet the potential impact of dust on the climate system, both in the anthropogenically disturbed future and the naturally varying past, remains to be quantified. The geologic record of dust provides the opportunity to test earth system models designed to simulate dust. Records of dust can be obtained from ice cores, marine sediments, and terrestrial (loess) deposits. Although rarely unequivocal, these records document a variety of processes (source, transport and deposition) in the dust cycle, stored in each archive as changes in clay mineralogy, isotopes, grain size, and concentration of terrigenous materials. Although the extraction of information from each type of archive is slightly different, the basic controls on these dust indicators are the same. Changes in the dust flux and particle size might be controlled by a combination of (a) source area extent, (b) dust emission efficiency (wind speed) and atmospheric transport, (c) atmospheric residence time of dust, and/or (d) relative contributions of dry settling and rainout of dust. Similarly, changes in mineralogy reflect (a) source area mineralogy and weathering and (b) shifts in atmospheric transport. The combination of these geological data with process-based, forward-modelling schemes in global earth system models provides an excellent means of achieving a comprehensive picture of the global pattern of dust accumulation rates, their controlling mechanisms, and how those mechanisms may vary regionally. The Dust Indicators and Records of Terrestrial and MArine Palaeoenvironments (DIRTMAP) data base has been established to provide a global palaeoenvironmental data set that can be used to validate earth system model simulations of the dust cycle over the past 150,000 years.

Atmospheric electric fields during the Carrington flare

Karen Aplin and Giles Harrison examine international records of the 1859 Carrington flare and consider what they mean for our understanding of space weather today. Space weather is increasingly recognized as a hazard to modern societies, and one way to assess the extent of its possible impact is through analysis of historic space weather events. One such event was the massive solar storm of late August and early September 1859. This is now widely known as the “Carrington flare” or “Carrington event” after the visual solar emissions on 1 September first reported by the Victorian astronomer Richard Carrington from his observatory in Redhill, Surrey (Carrington 1859). The related aurorae and subsequent effects on telegraph networks are well documented (e.g. Clark 2007, Boteler 2006), but use of modern techniques, such as analysis of nitrates produced by solar protons in ice cores to retrospectively assess the nature of the solar flare, has proved problematic (Wolff et al. 2012). This means that there is still very little quantitative information about the flare beyond magnetic observations (e.g. Viljanen et al. 2014).

Our changing Sun

Stellar astronomy tells us much about the long-term evolution of our Sun while forensic evidence (for example, cosmic-ray products in ice cores) gives us indications of its fluctuations over the last millennium. However, such studies do not give us a sufficiently detailed understanding of solar change over the last century to allow us to detect and quantify any role that the Sun might have played in the observed rise in average surface temperatures on Earth. This paper describes recent research that has filled this gap by applying advances in our understanding of the effects and structure of the solar wind to historical data on the Earth's magnetic field.

Mid- to late Holocene Indian Ocean monsoon variability recorded in four speleothems from Socotra Island, Yemen

Four stalagmites covering the last 7.0 ka were sampled on Socotra, an island in the northern Indian Ocean to investigate the evolution of the northeast Indian Ocean Monsoon (IOM) since the mid Holocene. On Socotra, rain is delivered at the start of the southwest IOM in May–June and at the start of the northeast IOM from September to December. The Haggeher Mountains act as a barrier forcing precipitation brought by the northeast winds to fall preferentially on the eastern side of the island, where the studied caves are located. δ18O and δ13C and Mg/Ca and Sr/Ca signals in the stalagmites reflect precipitation amounts brought by the northeast winds. For stalagmite STM6, this amount effect is amplified by kinetic effects during calcite deposition. Combined interpretation of the stalagmites' signals suggest a weakening of the northeast precipitation between 6.0 and 3.8 ka. After 3.8 ka precipitation intensities remain constant with two superimposed drier periods, between 0 and 0.6 ka and from 2.2 to 3.8 ka. No link can be established with Greenland ice cores and with the summer IOM variability. In contrast to the stable northeast rainy season suggested by the records in this study, speleothem records from western Socotra indicate a wettening of the southwest rainy season on Socotra after 4.4 ka. The local wettening of western Socotra could relate to a more southerly path (more over the Indian Ocean) taken by the southwest winds. Stalagmite STM5, sampled at the fringe between both rain areas displays intermediate δ18O values. After 6.2 ka, similar precipitation changes are seen between eastern Socotra and northern Oman indicating that both regions are affected similarly by the monsoon. Different palaeoclimatologic records from the Arabian Peninsula currently located outside the ITCZ migration pathway display an abrupt drying around 6 ka due to their disconnection from the southwest rain influence. Records that are nowadays still receiving rain by the southwest winds, suggest a more gradual drying reflecting the weakening of the southwest monsoon.

The heliospheric Hale cycle over the last 300 years and its implications for a “lost” late 18th century solar cycle

A Hale cycle, one complete magnetic cycle of the Sun, spans two complete Schwabe cycles (also referred to as sunspot and, more generally, solar cycles). The approximately 22-year Hale cycle is seen in magnetic polarities of both sunspots and polar fields, as well as in the intensity of galactic cosmic rays reaching Earth, with odd- and even-numbered solar cycles displaying qualitatively different waveforms. Correct numbering of solar cycles also underpins empirical cycle-to-cycle relations which are used as first-order tests of stellar dynamo models. There has been much debate about whether the unusually long solar cycle 4 (SC4), spanning- 1784–1799, was actually two shorter solar cycles combined as a result of poor data coverage in the original Wolf sunspot number record. Indeed, the group sunspot number does show a small increase around 1794–1799 and there is evidence of an increase in the mean latitude of sunspots at this time, suggesting the existence of a cycle ‘‘4b’’. In this study, we use cosmogenic radionuclide data and associated reconstructions of the heliospheric magnetic field (HMF) to show that the Hale cycle has persisted over the last 300 years and that data prior to 1800 are more consistent with cycle 4 being a single long cycle (the ‘‘no SC4b’’ scenario). We also investigate the effect of cycle 4b on the HMF using an open solar flux (OSF) continuity model, in which the OSF source term is related to sunspot number and the OSF loss term is determined by the heliospheric current sheet tilt, assumed to be a simple function of solar cycle phase. The results are surprising; Without SC4b, the HMF shows two distinct peaks in the 1784–1799 interval, while the addition of SC4b removes the secondary peak, as the OSF loss term acts in opposition to the later rise in sunspot number. The timing and magnitude of the main SC4 HMF peak is also significantly changed by the addition of SC4b. These results are compared with the cosmogenic isotope reconstructions of HMF and historical aurora records. These data marginally favour the existence of SC4b (the ‘‘SC4b’’ scenario), though the result is less certain than that based on the persistence of the Hale cycle. Thus while the current uncertainties in the observations preclude any definitive conclusions, the data favour the ‘‘no SC4b’’ scenario. Future improvements to cosmogenic isotope reconstructions of the HMF, through either improved modelling or additional ice cores from well-separated geographic locations, may enable questions of the existence of SC4b and the phase of Hale cycle prior to the Maunder minimum to be settled conclusively.

Late pleistocene ice export events into the Arctic Ocean from the M'Clure Strait Ice Stream, Canadian Arctic Archipelago

Rapidly-flowing sectors of an ice sheet (ice streams) can play ail important role in abrupt climate change through tile delivery of icebergs and meltwater and tile Subsequent disruption of ocean thermohaline circulation (e.g., the North Atlantic's Heinrich events). Recently, several cores have been raised from the Arctic Ocean which document the existence of massive ice export events during tile Late Pleistocene and whose provenance has been linked to Source regions in the Canadian Arctic Archipelago. In this paper, satellite imagery is used to map glacial geomorphology in the vicinity of Victoria Island, Banks Island and Prince of Wales Island (Canadian Arctic) in order to reconstruct ice flow patterns in the highly complex glacial landscape. A total of 88 discrete flow-sets are mapped and of these, 13 exhibit the characteristic geomorphology of palaeo-ice streams (i.e., parallel patterns of large, highly elongated mega-scale glacial lineations forming a convergent flow pattern with abrupt lateral margins). Previous studies by other workers and cross-cutting relationships indicate that the majority of these ice streams are relatively young and operated during or immediately prior to deglaciation. Our new mapping, however, documents a large (> 700 km long; 110 km wide) and relatively old ice stream imprint centred in M'Clintock Channel and converging into Viscount Melville Sound. A trough mouth fan located on the continental shelf Suggests that it extended along M'Clure Strait and was grounded at tile shelf edge. The location of the M'Clure Strait Ice Stream exactly matches the Source area of 4 (possibly 5) major ice export events recorded in core PS 1230 raised from Fram Strait, the major ice exit for the Arctic Ocean. These ice export events occur at similar to 12.9, similar to 15.6, similar to 22 and 29.8 ka (C-14 yr BP) and we argue that they record vigorous episodes of activity of the M'Clure Strait Ice Stream. The timing of these events is remarkably similar to the North Atlantic's Heinrich events and we take this as evidence that the M'Clure Strait Ice Stream was also activated around the same time. This may hold important implications for tile cause of the North Atlantic's Heinrich events and hints at tile possibility of a pall-ice sheet response. (c) 2005 Elsevier B.V. All rights reserved.

Ice supersaturation in the ECMWF integrated forecast system

A parametrization for ice supersaturation is introduced into the ECMWF Integrated Forecast System (IFS), compatible with the cloud scheme that allows partial cloud coverage. It is based on the simple, but often justifiable, diagnostic assumption that the ice nucleation and subsequent depositional growth time-scales are short compared to the model time step, thus supersaturation is only permitted in the clear-sky portion of the grid cell. Results from model integrations using the new scheme are presented, which is demonstrated to increase upper-tropospheric humidity, decrease high-level cloud cover and, to a much lesser extent, cloud ice amounts, all as expected from simple arguments. Evaluation of the relative distribution of supersaturated humidity amounts shows good agreement with the observed climatology derived from in situ aircraft observations. With the new scheme, the global distribution of frequency of occurrence of supersaturated regions compares well with remotely sensed microwave limb sounder (MLS) data, with the most marked errors of underprediction occurring in regions where the model is known to underpredict deep convection. Finally, it is also demonstrated that the new scheme leads to improved predictions of permanent contrail cloud over southern England, which indirectly implies upper-tropospheric humidity fields are better represented for this region.

Theory and observations of ice particle evolution in cirrus using Doppler radar: evidence for aggregation

Vertically pointing Doppler radar has been used to study the evolution of ice particles as they sediment through a cirrus cloud. The measured Doppler fall speeds, together with radar-derived estimates for the altitude of cloud top, are used to estimate a characteristic fall time tc for the `average' ice particle. The change in radar reflectivity Z is studied as a function of tc, and is found to increase exponentially with fall time. We use the idea of dynamically scaling particle size distributions to show that this behaviour implies exponential growth of the average particle size, and argue that this exponential growth is a signature of ice crystal aggregation.