Triboelectric charging of volcanic ash from the 2011 Grímsvötn eruption

The plume from the 2011 eruption of Grímsvötn was highly electrically charged, as shown by the considerable lightning activity measured by the United Kingdom Met Office’s low-frequency lightning detection network. Previous measurements of volcanic plumes have shown that ash particles are electrically charged up to hundreds of kilometers away from the vent, which indicates that the ash continues to charge in the plume [R. G. Harrison, K. A. Nicoll, Z. Ulanowski, and T. A. Mather, Environ. Res. Lett. 5 024004 (2010); H. Hatakeyama J. Meteorol. Soc. Jpn. 27 372 (1949)]. In this Letter, we study triboelectric charging of different size fractions of a sample of volcanic ash experimentally. Consistently with previous work, we find that the particle size distribution is a determining factor in the charging. Specifically, our laboratory experiments demonstrate that the normalized span of the particle size distribution plays an important role in the magnitude of charging generated. The influence of the normalized span on plume charging suggests that all ash plumes are likely to be charged, with implications for remote sensing and plume lifetime through scavenging effects.

High-resolution provenance of desert dust deposited on Mt. Elbrus, Caucasus in 2009–2012 using snow pit and firn core records

The first record of dust deposition events on Mt. Elbrus, Caucasus Mountains derived from a snow pit and a shallow firn core is presented for the 2009–2012 period. A combination of isotopic analysis, SEVIRI red-greenblue composite imagery, MODIS atmospheric optical depth fields derived using the Deep Blue algorithm, air mass trajectories derived using the HYSPLIT model and analyses of meteorological data enabled identification of dust source regions with high temporal (hours) and spatial (ca. 20–100 km) resolution. Seventeen dust deposition events were detected; fourteen occurred in March–June, one in February and two in October. Four events originated in the Sahara, predominantly in northeastern Libya and eastern Algeria. Thirteen events originated in the Middle East, in the Syrian Desert and northern Mesopotamia, from a mixture of natural and anthropogenic sources. Dust transportation from Sahara was associated with vigorous Saharan depressions, strong surface winds in the source region and mid-tropospheric southwesterly flow with daily winds speeds of 20–30 m s−1 at 700 hPa level. Although these events were less frequent than those originating in the Middle East, they resulted in higher dust concentrations in snow. Dust transportation from the Middle East was associated with weaker depressions forming over the source region, high pressure centred over or extending towards the Caspian Sea and a weaker southerly or southeasterly flow towards the Caucasus Mountains with daily wind speeds of 12–18 m s−1 at 700 hPa level. Higher concentrations of nitrates and ammonium characterised dust from the Middle East deposited on Mt. Elbrus in 2009 indicating contribution of anthropogenic sources. The modal values of particle size distributions ranged between 1.98 µm and 4.16 µm. Most samples were characterised by modal values of 2.0– 2.8 µm with an average of 2.6 µm and there was no signifi- cant difference between dust from the Sahara and the Middle East.

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.

On the barrier properties of the cornea: a microscopy study of the penetration of fluorescently labeled nanoparticles, polymers, and sodium fluorescein

Overcoming the natural defensive barrier functions of the eye remains one of the greatest challenges of ocular drug delivery. Cornea is a chemical and mechanical barrier preventing the passage of any foreign bodies including drugs into the eye, but the factors limiting penetration of permeants and nanoparticulate drug delivery systems through the cornea are still not fully understood. In this study, we investigate these barrier properties of the cornea using thiolated and PEGylated (750 and 5000 Da) nanoparticles, sodium fluorescein, and two linear polymers (dextran and polyethylene glycol). Experiments used intact bovine cornea in addition to bovine cornea de-epithelialized or tissues pretreated with cyclodextrin. It was shown that corneal epithelium is the major barrier for permeation; pretreatment of the cornea with β-cyclodextrin provides higher permeation of low molecular weight compounds, such as sodium fluorescein, but does not enhance penetration of nanoparticles and larger molecules. Studying penetration of thiolated and PEGylated (750 and 5000 Da) nanoparticles into the de-epithelialized ocular tissue revealed that interactions between corneal surface and thiol groups of nanoparticles were more significant determinants of penetration than particle size (for the sizes used here). PEGylation with polyethylene glycol of a higher molecular weight (5000 Da) allows penetration of nanoparticles into the stroma, which proceeds gradually, after an initial 1 h lag phase.

Observing the in situ chiral modification of Ni nanoparticles using scanning transmission X-ray microspectroscopy

Enantioselective heterogeneous hydrogenation of Cdouble bond; length as m-dashO bonds is of great potential importance in the synthesis of chirally pure products for the pharmaceutical and fine chemical industries. One of the most widely studied examples of such a reaction is the hydrogenation of β-ketoesters and β-diketoesters over Ni-based catalysts in the presence of a chiral modifier. Here we use scanning transmission X-ray microscopy combined with near-edge X-ray absorption fine structure spectroscopy (STXM/NEXAFS) to investigate the adsorption of the chiral modifier, namely (R,R)-tartaric acid, onto individual nickel nanoparticles. The C K-edge spectra strongly suggest that tartaric acid deposited onto the nanoparticle surfaces from aqueous solutions undergoes a keto-enol tautomerisation. Furthermore, we are able to interrogate the Ni L2,3-edge resonances of individual metal nanoparticles which, combined with X-ray diffraction (XRD) patterns showed them to consist of a pure nickel phase rather than the more thermodynamically stable bulk nickel oxide. Importantly, there appears to be no “particle size effect” on the adsorption mode of the tartaric acid in the particle size range ~ 90–~ 300 nm.

Structure and magnetic properties of polydisperse ferrofluids: a molecular dynamics study

We study by Langevin molecular dynamics simulations systematically the influence of polydispersity in the particle size, and subsequently in the dipole moment, on the physical properties of ferrofluids. The polydispersity is in a first approximation modeled by a bidisperse system that consists of small and large particles at different ratios of their volume fractions. In the first part of our investigations the total volume fraction of the system is fixed, and the volume fraction phi(L) of the large particles is varied. The initial susceptibility chi and magnetization curve of the systems show a strong dependence on the value of phi(L). With the increase of phi(L), the magnetization M of the system has a much faster increment at weak fields, and thus leads to a larger chi. We performed a cluster analysis that indicates that this is due to the aggregation of the large particles in the systems. The average size of these clusters increases with increasing phi(L). In the second part of our investigations, we fixed the volume fraction of the large particles, and increased the volume fraction phi(S) of the small particles in order to study their influence on the chain formation of the large ones. We found that the average aggregate size formed by large particles decreases when phi(S) is increased, demonstrating a significant effect of the small particles on the structural properties of the system. A topological analysis of the structure reveals that the majority of the small particles remain nonaggregated. Only a small number of them are attracted to the ends of the chains formed by large particles.

Optical properties of Saharan dust aerosol and contribution from the coarse mode as measured during the Fennec 2011 aircraft campaign

New in-situ aircraft measurements of Saharan dust originating from Mali, Mauritania and Algeria taken during the Fennec 2011 aircraft campaign over a remote part of the Sahara Desert are presented. Size distributions extending to 300 μm are shown, representing measurements extending further into the coarse mode than previously published for airborne Saharan dust. A significant coarse mode was present in the size distribution measurements with effective diameter (deff) from 2.3 to 19.4 μm and coarse mode volume median diameter (dvc) from 5.8 to 45.3 μm. The mean size distribution had a larger relative proportion of coarse mode particles than previous aircraft measurements. The largest particles (with deff >12 μm, or dvc >25 μm) were only encountered within 1 km of the ground. Number concentration, mass loading and extinction coefficient showed inverse relationships to dust age since uplift. Dust particle size showed a weak exponential relationship to dust age. Two cases of freshly uplifted dust showed quite different characteristics of size distribution and number concentration. Single Scattering Albed (SSA) values at 550 nm calculated from the measured size distributions revealed high absorption ranging from 0.70 to 0.97 depending on the refractive index. SSA was found to be strongly related to deff. New instrumentation revealed that direct measurements, behind Rosemount inlets, overestimate SSA by up to 0.11 when deff is greater than 2 μm. This is caused by aircraft inlet inefficiencies and sampling losses. Previous measurements of SSA from aircraft measurements may also have been overestimates for this reason. Radiative transfer calculations indicate that the range of SSAs during Fennec 2011 can lead to underestimates in shortwave atmospheric heating rates by 2.0 to 3.0 times if the coarse mode is neglected. This will have an impact on Saharan atmospheric dynamics and circulation,which should be taken into account by numerical weather prediction and climate models.

Impact of atmospheric transport on the evolution of microphysical and optical properties of Saharan dust

Saharan dust affects the climate by altering the radiation balance and by depositing minerals to the Atlantic Ocean. Both are dependent on particle size. We present aircraft measurements comprising 42 profiles of size distribution (0.1–300 µm), representing freshly uplifted dust, regional aged dust, and dust in the Saharan Air Layer (SAL) over the Canary Islands. The mean effective diameter of dust in SAL profiles is 4.5 µm smaller than that in freshly uplifted dust, while the vertical structure changes from a low shallow layer (0–1.5 km) to a well-mixed deep Saharan dust layer (0–5 km). Size distributions show a loss of 60 to 90% of particles larger than 30 µm 12 h after uplift. The single scattering albedo (SSA) increases from 0.92 to 0.94 to 0.95 between fresh, aged, and SAL profiles: this is enough to alter heating rates by 26%. Some fresh dust close to the surface shows SSA as low as 0.85

Atmospheric chemistry and physics in the atmosphere of a developed megacity (London): an overview of the REPARTEE experiment and its conclusions

The REgents PARk and Tower Environmental Experiment (REPARTEE) comprised two campaigns in London in October 2006 and October/November 2007. The experiment design involved measurements at a heavily trafficked roadside site, two urban background sites and an elevated site at 160–190 m above ground on the BT Tower, supplemented in the second campaign by Doppler lidar measurements of atmospheric vertical structure. A wide range of measurements of airborne particle physical metrics and chemical composition were made as well as measurements of a considerable range of gas phase species and the fluxes of both particulate and gas phase substances. Significant findings include (a) demonstration of the evaporation of traffic-generated nanoparticles during both horizontal and vertical atmospheric transport; (b) generation of a large base of information on the fluxes of nanoparticles, accumulation mode particles and specific chemical components of the aerosol and a range of gas phase species, as well as the elucidation of key processes and comparison with emissions inventories; (c) quantification of vertical gradients in selected aerosol and trace gas species which has demonstrated the important role of regional transport in influencing concentrations of sulphate, nitrate and secondary organic compounds within the atmosphere of London; (d) generation of new data on the atmospheric structure and turbulence above London, including the estimation of mixed layer depths; (e) provision of new data on trace gas dispersion in the urban atmosphere through the release of purposeful tracers; (f) the determination of spatial differences in aerosol particle size distributions and their interpretation in terms of sources and physico-chemical transformations; (g) studies of the nocturnal oxidation of nitrogen oxides and of the diurnal behaviour of nitrate aerosol in the urban atmosphere, and (h) new information on the chemical composition and source apportionment of particulate matter size fractions in the atmosphere of London derived both from bulk chemical analysis and aerosol mass spectrometry with two instrument types.

Remarkable dynamics of nanoparticles in the urban atmosphere

Nanoparticles emitted from road traffic are the largest source of respiratory exposure for the general public living in urban areas. It has been suggested that adverse health effects of airborne particles may scale with airborne particle number, which if correct, focuses attention on the nanoparticle (less than 100 nm) size range which dominates the number count in urban areas. Urban measurements of particle size distributions have tended to show a broadly similar pattern dominated by a mode centred on 20–30 nm diameter emitted by diesel engine exhaust. In this paper we report the results of measurements of particle number concentration and size distribution made in a major London park as well as on the BT Tower, 160 m aloft. These measurements taken during the REPARTEE project (Regents Park and BT Tower experiment) show a remarkable shift in particle size distributions with major losses of the smallest particle class as particles are advected away from the traffic source. In the Park, the traffic related mode at 20–30 nm diameter is much reduced with a new mode at <10 nm. Size distribution measurements also revealed higher number concentrations of sub-50 nm particles at the BT Tower during days affected by higher turbulence as determined by Doppler Lidar measurements and are indicative of loss of nanoparticles from air aged during less turbulent conditions. These results are suggestive of nanoparticle loss by evaporation, rather than coagulation processes. The results have major implications for understanding the impacts of traffic-generated particulate matter on human health.

Intercomparison and evaluation of global aerosol microphysical properties among AeroCom models of a range of complexity

Many of the next generation of global climate models will include aerosol schemes which explicitly simulate the microphysical processes that determine the particle size distribution. These models enable aerosol optical properties and cloud condensation nuclei (CCN) concentrations to be determined by fundamental aerosol processes, which should lead to a more physically based simulation of aerosol direct and indirect radiative forcings. This study examines the global variation in particle size distribution simulated by 12 global aerosol microphysics models to quantify model diversity and to identify any common biases against observations. Evaluation against size distribution measurements from a new European network of aerosol supersites shows that the mean model agrees quite well with the observations at many sites on the annual mean, but there are some seasonal biases common to many sites. In particular, at many of these European sites, the accumulation mode number concentration is biased low during winter and Aitken mode concentrations tend to be overestimated in winter and underestimated in summer. At high northern latitudes, the models strongly underpredict Aitken and accumulation particle concentrations compared to the measurements, consistent with previous studies that have highlighted the poor performance of global aerosol models in the Arctic. In the marine boundary layer, the models capture the observed meridional variation in the size distribution, which is dominated by the Aitken mode at high latitudes, with an increasing concentration of accumulation particles with decreasing latitude. Considering vertical profiles, the models reproduce the observed peak in total particle concentrations in the upper troposphere due to new particle formation, although modelled peak concentrations tend to be biased high over Europe. Overall, the multi-model-mean data set simulates the global variation of the particle size distribution with a good degree of skill, suggesting that most of the individual global aerosol microphysics models are performing well, although the large model diversity indicates that some models are in poor agreement with the observations. Further work is required to better constrain size-resolved primary and secondary particle number sources, and an improved understanding of nucleation and growth (e.g. the role of nitrate and secondary organics) will improve the fidelity of simulated particle size distributions.

Aerosol microphysics simulations of the Mt.~Pinatubo eruption with the UM-UKCA composition-climate model

We use a stratosphere–troposphere composition–climate model with interactive sulfur chemistry and aerosol microphysics, to investigate the effect of the 1991 Mount Pinatubo eruption on stratospheric aerosol properties. Satellite measurements indicate that shortly after the eruption, between 14 and 23 Tg of SO2 (7 to 11.5 Tg of sulfur) was present in the tropical stratosphere. Best estimates of the peak global stratospheric aerosol burden are in the range 19 to 26 Tg, or 3.7 to 6.7 Tg of sulfur assuming a composition of between 59 and 77 % H2SO4. In light of this large uncertainty range, we performed two main simulations with 10 and 20 Tg of SO2 injected into the tropical lower stratosphere. Simulated stratospheric aerosol properties through the 1991 to 1995 period are compared against a range of available satellite and in situ measurements. Stratospheric aerosol optical depth (sAOD) and effective radius from both simulations show good qualitative agreement with the observations, with the timing of peak sAOD and decay timescale matching well with the observations in the tropics and mid-latitudes. However, injecting 20 Tg gives a factor of 2 too high stratospheric aerosol mass burden compared to the satellite data, with consequent strong high biases in simulated sAOD and surface area density, with the 10 Tg injection in much better agreement. Our model cannot explain the large fraction of the injected sulfur that the satellite-derived SO2 and aerosol burdens indicate was removed within the first few months after the eruption. We suggest that either there is an additional alternative loss pathway for the SO2 not included in our model (e.g. via accommodation into ash or ice in the volcanic cloud) or that a larger proportion of the injected sulfur was removed via cross-tropopause transport than in our simulations. We also critically evaluate the simulated evolution of the particle size distribution, comparing in detail to balloon-borne optical particle counter (OPC) measurements from Laramie, Wyoming, USA (41° N). Overall, the model captures remarkably well the complex variations in particle concentration profiles across the different OPC size channels. However, for the 19 to 27 km injection height-range used here, both runs have a modest high bias in the lowermost stratosphere for the finest particles (radii less than 250 nm), and the decay timescale is longer in the model for these particles, with a much later return to background conditions. Also, whereas the 10 Tg run compared best to the satellite measurements, a significant low bias is apparent in the coarser size channels in the volcanically perturbed lower stratosphere. Overall, our results suggest that, with appropriate calibration, aerosol microphysics models are capable of capturing the observed variation in particle size distribution in the stratosphere across both volcanically perturbed and quiescent conditions. Furthermore, additional sensitivity simulations suggest that predictions with the models are robust to uncertainties in sub-grid particle formation and nucleation rates in the stratosphere.

A coupled cloud physics–radiation parameterization of the bulk optical properties of cirrus and its impact on the Met Office unified model global atmosphere 5.0 configuration

A new coupled cloud physics–radiation parameterization of the bulk optical properties of ice clouds is presented. The parameterization is consistent with assumptions in the cloud physics scheme regarding particle size distributions (PSDs) and mass–dimensional relationships. The parameterization is based on a weighted ice crystal habit mixture model, and its bulk optical properties are parameterized as simple functions of wavelength and ice water content (IWC). This approach directly couples IWC to the bulk optical properties, negating the need for diagnosed variables, such as the ice crystal effective dimension. The parameterization is implemented into the Met Office Unified Model Global Atmosphere 5.0 (GA5) configuration. The GA5 configuration is used to simulate the annual 20-yr shortwave (SW) and longwave (LW) fluxes at the top of the atmosphere (TOA), as well as the temperature structure of the atmosphere, under various microphysical assumptions. The coupled parameterization is directly compared against the current operational radiation parameterization, while maintaining the same cloud physics assumptions. In this experiment, the impacts of the two parameterizations on the SW and LW radiative effects at TOA are also investigated and compared against observations. The 20-yr simulations are compared against the latest observations of the atmospheric temperature and radiative fluxes at TOA. The comparisons demonstrate that the choice of PSD and the assumed ice crystal shape distribution are as important as each other. Moreover, the consistent radiation parameterization removes a long-standing tropical troposphere cold temperature bias but slightly warms the southern midlatitudes by about 0.5 K.

Microphysical properties of cold frontal rainbands

Observations have been obtained within an intense (precipitation rates > 50 mm h−1 ) narrow cold-frontal rainband (NCFR) embedded within a broader region of stratiform precipitation. In situ data were obtained from an aircraft which flew near a steerable dual-polarisation Doppler radar. The observations were obtained to characterise the microphysical properties of cold frontal clouds, with an emphasis on ice and precipitation formation and development. Primary ice nucleation near cloud top (−55◦ C) appeared to be enhanced by convective features. However, ice multiplication led to the largest ice particle number concentrations being observed at relatively high temperatures (> −10◦ C). The multiplication process (most likely rime splintering) occurs when stratiform precipitation interacts with supercooled water generated in the NCFR. Graupel was notably absent in the data obtained. Ice multiplication processes are known to have a strong impact in glaciating isolated convective clouds, but have rarely been studied within larger organised convective systems such as NCFRs. Secondary ice particles will impact on precipitation formation and cloud dynamics due to their relatively small size and high number density. Further modelling studies are required to quantify the effects of rime splintering on precipitation and dynamics in frontal rainbands. Available parametrizations used to diagnose the particle size distributions do not account for the influence of ice multiplication. This deficiency in parametrizations is likely to be important in some cases for modelling the evolution of cloud systems and the precipitation formation. Ice multiplication has significant impact on artefact removal from in situ particle imaging probes.

Simple measures of climate, soil properties and plant traits predict national scale grassland soil carbon stocks

1. Soil carbon (C) storage is a key ecosystem service. Soil C stocks play a vital role in soil fertility and climate regulation, but the factors that control these stocks at regional and national scales are unknown, particularly when their composition and stability are considered. As a result, their mapping relies on either unreliable proxy measures or laborious direct measurements. 2. Using data from an extensive national survey of English grasslands we show that surface soil (0-7cm) C stocks in size fractions of varying stability can be predicted at both regional and national scales from plant traits and simple measures of soil and climatic conditions. 3. Soil C stocks in the largest pool, of intermediate particle size (50-250 µm), were best explained by mean annual temperature (MAT), soil pH and soil moisture content. The second largest C pool, highly stable physically and biochemically protected particles (0.45-50 µm), was explained by soil pH and the community abundance weighted mean (CWM) leaf nitrogen (N) content, with the highest soil C stocks under N rich vegetation. The C stock in the small active fraction (250-4000 µm) was explained by a wide range of variables: MAT, mean annual precipitation, mean growing season length, soil pH and CWM specific leaf area; stocks were higher under vegetation with thick and/or dense leaves. 4. Testing the models describing these fractions against data from an independent English region indicated moderately strong correlation between predicted and actual values and no systematic bias, with the exception of the active fraction, for which predictions were inaccurate. 5. Synthesis and Applications: Validation indicates that readily available climate, soils and plant survey data can be effective in making local- to landscape-scale (1-100,000 km2) soil C stock predictions. Such predictions are a crucial component of effective management strategies to protect C stocks and enhance soil C sequestration.