Dust fluxes and iron fertilization in Holocene and Last Glacial Maximum climates

Mineral dust aerosols play a major role in present and past climates. To date, we rely on climate models for estimates of dust fluxes to calculate the impact of airborne micronutrients on biogeochemical cycles. Here we provide a new global dust flux data set for Holocene and Last Glacial Maximum (LGM) conditions based on observational data. A comparison with dust flux simulations highlights regional differences between observations and models. By forcing a biogeochemical model with our new data set and using this model's results to guide a millennial-scale Earth System Model simulation, we calculate the impact of enhanced glacial oceanic iron deposition on the LGM-Holocene carbon cycle. On centennial timescales, the higher LGM dust deposition results in a weak reduction of <10?ppm in atmospheric CO2 due to enhanced efficiency of the biological pump. This is followed by a further ~10?ppm reduction over millennial timescales due to greater carbon burial and carbonate compensation.

Advances in understanding mineral dust and boundary layer processes over the Sahara from Fennec aircraft observations

The Fennec climate program aims to improve understanding of the Saharan climate system through a synergy of observations and modelling. We present a description of the Fennec airborne observations during 2011 and 2012 over the remote Sahara (Mauritania and Mali) and the advances in the understanding of mineral dust and boundary layer processes they have provided. Aircraft instrumentation aboard the UK FAAM BAe146 and French SAFIRE Falcon 20 is described, with specific focus on instrumentation specially developed and relevant to Saharan meteorology and dust. Flight locations, aims and associated meteorology are described. Examples and applications of aircraft measurements from the Fennec flights are presented, highlighting new scientific results delivered using a synergy of different instruments and aircraft. These include: (1) the first airborne measurement of dust particles sized up to 300 microns and associated dust fluxes in the Saharan atmospheric boundary layer (SABL), (2) dust uplift from the breakdown of the nocturnal low-level jet before becoming visible in SEVIRI satellite imagery, (3) vertical profiles of the unique vertical structure of turbulent fluxes in the SABL, (4) in-situ observations of processes in SABL clouds showing dust acting as CCN and IN at −15 °C, (5) dual-aircraft observations of the SABL dynamics, thermodynamics and composition in the Saharan heat low region (SHL), (6) airborne observations of a dust storm associated with a cold-pool (haboob) issued from deep convection over the Atlas, (7) the first airborne chemical composition measurements of dust in the SHL region with differing composition, sources (determined using Lagrangian backward trajectory calculations) and absorption properties between 2011 and 2012, (8) coincident ozone and dust surface area measurements suggest coarser particles provide a route for ozone depletion, (9) discrepancies between airborne coarse mode size distributions and AERONET sunphotometer retrievals under light dust loadings. These results provide insights into boundary layer and dust processes in the SHL region – a region of substantial global climatic importance.

Quantifying particle size and turbulent scale dependence of dust uplift in the Sahara using aircraft measurements

The first size-resolved airborne measurements of dust fluxes and the first dust flux measurements from the central Sahara are presented and compared with a parameterization by Kok (2011a). High-frequency measurements of dust size distribution were obtained from 0.16 to 300 µm diameter, and eddy covariance fluxes were derived. This is more than an order of magnitude larger size range than previous flux estimates. Links to surface emission are provided by analysis of particle drift velocities. Number flux is described by a −2 power law between 1 and 144 µm diameter, significantly larger than the 12 µm upper limit suggested by Kok (2011a). For small particles, the deviation from a power law varies with terrain type and the large size cutoff is correlated with atmospheric vertical turbulent kinetic energy, suggesting control by vertical transport rather than emission processes. The measured mass flux mode is in the range 30–100 µm. The turbulent scales important for dust flux are from 0.1 km to 1–10 km. The upper scale increases during the morning as boundary layer depth and eddy size increase. All locations where large dust fluxes were measured had large topographical variations. These features are often linked with highly erodible surface features, such as wadis or dunes. We also hypothesize that upslope flow and flow separation over such features enhance the dust flux by transporting large particles out of the saltation layer. The tendency to locate surface flux measurements in open, flat terrain means these favored dust sources have been neglected in previous studies.

Tracking eolian dust with helium and thorium: Impacts of grain size and provenance

Reconstructions of the deposition rate of windblown mineral dust in ocean sediments offer an important means of tracking past climate changes and of assessing the radiative and biogeochemical impacts of dust in past climates. Dust flux estimates in ocean sediments have commonly been based on the operationally defined lithogenic fraction of sediment samples. More recently, dust fluxes have been estimated from measurements of helium and thorium, as rare isotopes of these elements (He-3 and Th-230) allow estimates of sediment flux, and the dominant isotopes (He-4 and Th-232) are uniquely associated with the lithogenic fraction of marine sediments. In order to improve the fidelity of dust flux reconstructions based on He and Th, we present a survey of He and Th concentrations in sediments from dust source areas in East Asia, Australia and South America. Our data show systematic relationships between He and Th concentrations and grain size, with He concentrations decreasing and Th concentrations increasing with decreasing grain size. We find consistent He and Th concentrations in the fine fraction (<5 μm) of samples from East Asia, Australia and Central South America (Puna-Central West Argentina), with Th concentrations averaging 14 μg/g and He concentrations averaging 2 μcc STP/g. We recommend use of these values for estimating dust fluxes in sediments where dust is dominantly fine-grained, and suggest that previous studies may have systematically overestimated Th-based dust fluxes by 30%. Source areas in Patagonia appear to have lower He and Th contents than other regions, as fine fraction concentrations average 0.8 μcc STP/g and 9 μg/g for 4He and 232Th, respectively. The impact of grain size on lithogenic He and Th concentrations should be taken into account in sediments proximal to dust sources where dust grain size may vary considerably. Our data also have important implications for the hosts of He in long-traveled dust and for the 3He/4He ratio used for terrigenous He in studies of extraterrestrial He in sediments and ice. We also investigate the use of He/Th ratios as a provenance tracer. Our results suggest differences in fine fraction He/Th ratios between East Asia, Australia, central South America and Patagonia, with ratios showing a positive relationship with the geological age of source rocks. He/Th ratios may thus provide useful provenance information, for example allowing separation of Patagonian sources from Puna-Central West Argentina or Australian dust sources. He/Th ratios in open-ocean marine sediments are similar to ratios in the fine fraction of upwind dust source areas. He/Th ratios in mid-latitude South Atlantic sediments suggest that dust in this region primarily derives from the Puna-Central West Argentina region (23–32°S) rather than Patagonia (>38°S). In the equatorial Pacific, He/Th ratios are much lower than in extratropical Pacific sediments or potential source areas measured as a part of this study (East Asia, South America, Australia) for reasons that are at present unclear, complicating their use as provenance tracers in this region.

Reconstruction of dust flux of ODP Hole 198-1209C

Eolian dust in pelagic deep sea sediments can be used to reconstruct ancient wind patterns and paleoenvironmental response to climate change. Traditional methods to determine dust accumulation involve isolating the non-dissolvable aluminosilicate minerals from deep sea sediments through a series of chemical leaches, but cannot differentiate between minerals from eolian, authigenic and volcanogenic sources. Other geochemical proxies, such as sedimentary 232Th and crustal 4He content, have been used to construct high-resolution records of atmospheric dust fluxes to the deep sea during the Quaternary. Here we use sedimentary Th content as a proxy for terrigenous material (eolian dust) in ~58 Myr-old sediments from the Shatsky Rise (ODP Site 1209) and compare our results with previous dust estimates generated using the traditional chemical extraction method and sedimentary 4He(crustal) concentrations. We find excellent agreement between Th-based dust estimates and those generated using the traditional method. In addition our results show a correlation between sedimentary Th and 4He(crustal) content, which suggests a source older than present day Asian loess supplied dust to the central subtropical Pacific Ocean during the early Paleogene.

Chemical concentrations and fluxes from EPICA ice cores EDML and EDC

Continuous sea salt and mineral dust aerosol records have been studied on the two EPICA (European Project for Ice Coring in Antarctica) deep ice cores. The joint use of these records from opposite sides of the East Antarctic plateau allows for an estimate of changes in dust transport and emission intensity as well as for the identification of regional differences in the sea salt aerosol source. The mineral dust flux records at both sites show a strong coherency over the last 150 kyr related to dust emission changes in the glacial Patagonian dust source with three times higher dust fluxes in the Atlantic compared to the Indian Ocean sector of the Southern Ocean (SO). Using a simple conceptual transport model this indicates that transport can explain only 40% of the atmospheric dust concentration changes in Antarctica, while factor 5-10 changes occurred. Accordingly, the main cause for the strong glacial dust flux changes in Antarctica must lie in environmental changes in Patagonia. Dust emissions, hence environmental conditions in Patagonia, were very similar during the last two glacials and interglacials, respectively, despite 2-4 °C warmer temperatures recorded in Antarctica during the penultimate interglacial than today. 2-3 times higher sea salt fluxes found in both ice cores in the glacial compared to the Holocene are difficult to reconcile with a largely unchanged transport intensity and the distant open ocean source. The substantial glacial enhancements in sea salt aerosol fluxes can be readily explained assuming sea ice formation as the main sea salt aerosol source with a significantly larger expansion of (summer) sea ice in the Weddell Sea than in the Indian Ocean sector. During the penultimate interglacial, our sea salt records point to a 50% reduction of winter sea ice coverage compared to the Holocene both in the Indian and Atlantic Ocean sector of the SO. However, from 20 to 80 ka before present sea salt fluxes show only very subdued millennial changes despite pronounced temperature fluctuations, likely due to the large distance of the sea ice salt source to our drill sites.

Map of North African dust emission

Changes in the emission, transport and deposition of aeolian dust have profound effects on regional climate, so that characterizing the lifecycle of dust in observations and improving the representation of dust in global climate models is necessary. A fundamental aspect of characterizing the dust cycle is quantifying surface dust fluxes, yet no spatially explicit estimates of this flux exist for the World's major source regions. Here we present a novel technique for creating a map of the annual mean emitted dust flux for North Africa based on retrievals of dust storm frequency from the Meteosat Second Generation Spinning Enhanced Visible and InfraRed Imager (SEVIRI) and the relationship between dust storm frequency and emitted mass flux derived from the output of five models that simulate dust. Our results suggest that 64 (±16)% of all dust emitted from North Africa is from the Bodélé depression, and that 13 (±3)% of the North African dust flux is from a depression lying in the lee of the Aïr and Hoggar Mountains, making this area the second most important region of emission within North Africa.

Reconstructions of eolian dust accumulation in northwest African margin sediments

Reconstructions of eolian dust accumulation in northwest African margin sediments provide important continuous records of past changes in atmospheric circulation and aridity in the region. Existing records indicate dramatic changes in North African dust emissions over the last 20 ka, but the limited spatial extent of these records and the lack of high-resolution flux data do not allow us to determine whether changes in dust deposition occurred with similar timing, magnitude and abruptness throughout northwest Africa. Here we present new records from a meridional transect of cores stretching from 31°N to 19°N along the northwest African margin. By combining grain size endmember modeling with 230Th-normalized fluxes for the first time, we are able to document spatial and temporal changes in dust deposition under the North African dust plume throughout the last 20 ka. Our results provide quantitative estimates of the magnitude of dust flux changes associated with Heinrich Stadial 1, the Younger Dryas, and the African Humid Period (AHP; ~11.7-5 ka), offering robust targets for model-based estimates of the climatic and biogeochemical impacts of past changes in North African dust emissions. Our data suggest that dust fluxes between 8 and 6 ka were a factor of ~5 lower than average fluxes during the last 2 ka. Using a simple model to estimate the effects of bioturbation on dust input signals, we find that our data are consistent with abrupt, synchronous changes in dust fluxes in all cores at the beginning and end of the AHP. The mean ages of these transitions are 11.8±0.2 ka (1Sigma) and 4.9±0.2 ka, respectively.

Grain-size analysis and geochemistry of sediment core GeoB13602-1

Grain-size, terrigenous element and rock magnetic remanence data of Quaternary marine sediments retrieved at the NW African continental margin off Gambia (gravity core GeoB 13602-1, 13°32.71' N, 17°50.96'W) were jointly analyzed by end-member (EM) unmixing methods to distinguish and budget past terrigenous fluxes. We compare and cross-validate the identified single-parameter EM systems and develop a numerical strategy to calculate associated multi-parameter EM properties. One aeolian and two fluvial EMs were found. The aeolian EM is much coarser than the fluvial EMs and is associated with a lower goethite/hematite ratio, a higher relative concentration of magnetite and lower Al/Si and Fe/K ratios. Accumulation rates and grain sizes of the fluvial sediment appear to be primarily constrained by shore distance (i.e., sea-level fluctuations) and to a lesser extent by changes in hinterland precipitation. High dust fluxes occurred during the Last Glacial Maximum (LGM) and during Heinrich Stadials (HS) while the fluvial input remained unchanged. Our approach reveals that the LGM dust fluxes were ~7 times higher than today's. However, by far the highest dust accumulation occurred during HS 1 (~300 g m**-2 yr** -1), when dust fluxes were ~80 fold higher than today. Such numbers have not yet been reported for NW Africa, and emphasize strikingly different environmental conditions during HSs. They suggest that deflation rate and areal extent of HSs dust sources were much larger due to retreating vegetation covers. Beyond its regional and temporal scope, this study develops new, in principle, generally applicable strategies for multi-method end-member interpretation, validation and flux budgeting calibration.