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.

Effect of high pressure homogenisation and heat treatment on physical properties and stability of almond and hazelnut milks

The effect of high pressure homogenisation (HPH) and heat treatments on physicochemical properties and physical stability of almond and hazelnut milks was studied. Vegetable milks were obtained and homogenised by applying 62, 103 and 172 MPa (MF1, MF2 and MF3, respectively). Untreated and MF3 samples were also submitted to two different heat treatments (85 °C/30 min (LH) or 121 °C/15 min (HH)). Physical and structural properties of the products were greatly affected by heat treatments and HPH. In almond milk, homogenised samples showed a significant reduction in particle size, which turned from bimodal and polydisperse to monodisperse distributions. Particle surface charge, clarity and Whiteness Index were increased and physical stability of samples was improved, without affecting either viscosity or protein stability. Hazelnut beverages showed similar trends, but HPH notably increased their viscosity while change their rheological behaviour, which suggested changes in protein conformation. HH treatments caused an increment of particle size due to the formation oil droplet-protein body clusters, associated with protein denaturation. Samples submitted to the combined treatment MF3 and LH showed the greatest stability.

Enzyme assisted extraction of chitin from shrimp shells (Litopenaeus vannamei)

BACKGROUND: Chemical chitin extraction generates large amounts of wastes and increases partial deacetylation of the product. Therefore, the use of biological methods for chitin extraction is an interesting alternative. The effects of process conditions on enzyme assisted extraction of chitin from the shrimp shells in a systematic way were the focal points of this study. RESULTS: Demineralisation conditions of 25C, 20 min, shells-lactic acid ratio of 1:1.1 w/w; and shells-acetic acid ratio of 1:1.2 w/w, the maximum demineralisation values were 98.64 and 97.57% for lactic and acetic acids, respectively. A total protein removal efficiency of 91.10% by protease from Streptomyces griseus with enzyme-substrate ratio 55 U/g, pH 7.0 and incubation time 3 h is obtained when the particle size range is 50-25 μm, which was identified as the most critical factor. The X-ray diffraction and 13C NMR spectroscopy analysis showed that the lower percent crystallinity and higher degree of acetylation of chitin from enzyme assisted extraction may exhibit better solubility properties and less depolymerisation in comparison with chitin from the chemical extraction. CONCLUSION: The present work investigates the effects of individual factors on process yields, and it has shown that, if the particle size is properly controlled a reaction time of 3 h is more than enough for deproteination by protease. Physicochemical analysis indicated that the enzyme assisted production of chitin seems appropriate to extract chitin, possibly retaining its native structure.

Effect of air-drying pre-treatment on the characterization of forest soil carbon pools

Archived soils could represent a valuable resource for the spatio-temporal inventory of soil carbon stability. However, archived soils are usually air-dried before storage and the impact of a drying pretreatment on physically and chemically-defined C fractions has not yet been fully assessed. Through the comparison of field-moist and corresponding air-dried (at 25oC for 2 weeks) forest soil samples, we examined the effect of air-drying on: a) the quantity and the quality of cold- (CWEC) and hot-water (HWEC) extractable C and b) the concentration of C in physically isolated fractions (free- and intra-aggregate light and organo-mineral). Soil samples were collected from the organic (O) and mineral (A and B) horizons of three different forest soils from southeastern England: (i) Cambisol under Pine (Pinus nigra); (ii) Cambisol under Beech (Fagus sylvatica) and (iii) Gleysol under oak (Quercus robur). CWEC concentrations for dry samples were up to 2 times greater than for corresponding field moist samples and had significantly (p < 0.001) higher phenolic content. However, the effect of drying pretreatment on HWEC, its phenolic content was not significant (p > 0.05) for most samples. Dried soils had significantly (p < 0.001) higher concentrations of free light-C while having lower concentrations of intra-aggregate-C when compared to moist samples (p < 0.001). However, fine silt and clay fractions were not significantly affected by the drying pretreatment (p=0.789). Therefore, based on the results obtained from gleysol and cambisol forest soils studied here, C contents in hot-water extractions and fine particle size physical fractions (< 25µm) seem to be robust measurements for evaluating C fractions in dried stored forest soils. Further soil types should be tested to evaluate the wider generality of these findings.

The partial oxidation of methane over Pd/Al2O3 catalyst nanoparticles studied in-situ by near ambient-pressure x-ray photoelectron spectroscopy

Near ambient-pressure X-ray photoelectron spectroscopy (NAP-XPS) is used to study the chemical state of methane oxidation catalysts in-situ. Al2O3{supported Pd catalysts are prepared with different particle sizes ranging from 4 nm to 10 nm. These catalysts were exposed to conditions similar to those used in the partial oxidation of methane (POM) to syn-gas and simultaneously monitored by NAP-XPS and mass spectrometry. NAP-XPS data show changes in the oxidation state of the palladium as the temperature in- creases, from metallic Pd0 to PdO, and back to Pd0. Mass spectrometry shows an increase in CO production whilst the Pd is in the oxide phase, and the metal is reduced back under presence of newly formed H2. A particle size effect is observed, such that CH4 conversion starts at lower temperatures with larger sized particles from 6 nm to 10 nm. We find that all nanoparticles begin CH4 conversion at lower temperatures than polycrystalline Pd foil.

Effects of coastal managed retreat on mercury biogeochemistry

We investigated the impact of managed retreat on mercury (Hg) biogeochemistry at a site subject to diffuse contamination with Hg. We collected sediment cores from an area of land behind a dyke one year before and one year after it was intentionally breached. These sediments were compared to those of an adjacent mudflat and a salt marsh. The concentration of total mercury (THg) in the sediment doubled after the dyke was breached due to the deposition of fresh sediment that had a smaller particle size, and higher pH. The concentration of methylmercury (MeHg) was 27% lower in the sediments after the dyke was breached. We conclude that coastal flooding during managed retreat of coastal flood defences at this site has not increased the risk of Hg methylation or bioavailability during the first year. As the sediment becomes vegetated, increased activity of Hg-methylating bacteria may accelerate Hg-methylation rate.

Systematic and random errors between collocated satellite ice water path observations

There remains large disagreement between ice-water path (IWP) in observational data sets, largely because the sensors observe different parts of the ice particle size distribution. A detailed comparison of retrieved IWP from satellite observations in the Tropics (!30 " latitude) in 2007 was made using collocated measurements. The radio detection and ranging(radar)/light detection and ranging (lidar) (DARDAR) IWP data set, based on combined radar/lidar measurements, is used as a reference because it provides arguably the best estimate of the total column IWP. For each data set, usable IWP dynamic ranges are inferred from this comparison. IWP retrievals based on solar reflectance measurements, in the moderate resolution imaging spectroradiometer (MODIS), advanced very high resolution radiometer–based Climate Monitoring Satellite Applications Facility (CMSAF), and Pathfinder Atmospheres-Extended (PATMOS-x) datasets, were found to be correlated with DARDAR over a large IWP range (~20–7000 g m -2 ). The random errors of the collocated data sets have a close to lognormal distribution, and the combined random error of MODIS and DARDAR is less than a factor of 2, which also sets the upper limit for MODIS alone. In the same way, the upper limit for the random error of all considered data sets is determined. Data sets based on passive microwave measurements, microwave surface and precipitation products system (MSPPS), microwave integrated retrieval system (MiRS), and collocated microwave only (CMO), are largely correlated with DARDAR for IWP values larger than approximately 700 g m -2 . The combined uncertainty between these data sets and DARDAR in this range is slightly less MODIS-DARDAR, but the systematic bias is nearly an order of magnitude.

Applications of electrified dust and dust devil electrodynamics to Martian atmospheric electricity

Atmospheric transport and suspension of dust frequently brings electrification, which may be substantial. Electric fields of 10 kVm-1 to 100 kVm-1 have been observed at the surface beneath suspended dust in the terrestrial atmosphere, and some electrification has been observed to persist in dust at levels to 5 km, as well as in volcanic plumes. The interaction between individual particles which causes the electrification is incompletely understood, and multiple processes are thought to be acting. A variation in particle charge with particle size, and the effect of gravitational separation explains to, some extent, the charge structures observed in terrestrial dust storms. More extensive flow-based modelling demonstrates that bulk electric fields in excess of 10 kV m-1 can be obtained rapidly (in less than 10 s) from rotating dust systems (dust devils) and that terrestrial breakdown fields can be obtained. Modelled profiles of electrical conductivity in the Martian atmosphere suggest the possibility of dust electrification, and dust devils have been suggested as a mechanism of charge separation able to maintain current flow between one region of the atmosphere and another, through a global circuit. Fundamental new understanding of Martian atmospheric electricity will result from the ExoMars mission, which carries the DREAMS (Dust characterization, Risk Assessment, and Environment Analyser on the Martian Surface)-MicroARES (Atmospheric Radiation and Electricity Sensor) instrumentation to Mars in 2016 for the first in situ measurements.

Edible oil from Tiger nut (Cyperus esculentus L.): mechanical pressing and aqueous enzymatic extraction methods

The tiger nut tuber of the Cyperus esculentus L. plant is an unusual storage system with similar amounts of starch and lipid. The extraction of its oil employing both mechanical pressing and aqueous enzymatic extraction (AEE) methods was investigated and an examination of the resulting products was carried out. The effects of particle size and moisture content of the tuber on the yield of tiger nut oil with pressing were initially studied. Smaller particles were found to enhance oil yields while a range of moisture content was observed to favour higher oil yields. When samples were first subjected to high pressures up to 700 MPa before pressing at 38 MPa there was no increase in the oil yields. Ground samples incubated with a mixture of α- Amylase, Alcalase, and Viscozyme (a mixture of cell wall degrading enzyme) as a pre-treatment, increased oil yield by pressing and 90% of oil was recovered as a result. When aqueous enzymatic extraction was carried out on ground samples, the use of α- Amylase, Alcalase, and Celluclast independently improved extraction oil yields compared to oil extraction without enzymes by 34.5, 23.4 and 14.7% respectively. A mixture of the three enzymes further augmented the oil yield and different operational factors were individually studied for their effects on the process. These include time, total mixed enzyme concentration, linear agitation speed, and solid-liquid ratio. The largest oil yields were obtained with a solid-liquid ratio of 1:6, mixed enzyme concentration of 1% (w/w) and 6 h incubation time although the longer time allowed for the formation of an emulsion. Using stationary samples during incubation surprisingly gave the highest oil yields, and this was observed to be as a result of gravity separation occurring during agitation. Furthermore, the use of high pressure processing up to 300 MPa as a pre-treatment enhanced oil yields but additional pressure increments had a detrimental effect. The quality of oils recovered from both mechanical and aqueous enzymatic extraction based on the percentage free fatty acid (% FFA) and peroxide values (PV) all reflected the good stabilities of the oils with the highest % FFA of 1.8 and PV of 1.7. The fatty acid profiles of all oils also remained unchanged. The level of tocopherols in oils were enhanced with both enzyme aided pressing (EAP) and high pressure processing before AEE. Analysis on the residual meals revealed DP 3 and DP 4 oligosaccharides present in EAP samples but these would require further assessment on their identity and quality.

Whey protein mouth drying influenced by thermal denaturation

Whey proteins are becoming an increasingly popular functional food ingredient. There are, however, sensory properties associated with whey protein beverages that may hinder the consumption of quantities sufficient to gain the desired nutritional benefits. One such property is mouth drying. The influence of protein structure on the mouthfeel properties of milk proteins has been previously reported. This paper investigates the effect of thermal denaturation of whey proteins on physicochemical properties (viscosity, particle size, zeta-potential, pH), and relates this to the observed sensory properties measured by qualitative descriptive analysis and sequential profiling. Mouthcoating, drying and chalky attributes built up over repeated consumption, with higher intensities for samples subjected to longer heating times (p < 0.05). Viscosity, pH, and zeta-potential were found to be similar for all samples, however particle size increased with longer heating times. As the pH of all samples was close to neutral, this implies that neither the precipitation of whey proteins at low pH, nor their acidity, as reported in previous literature, can be the drying mechanisms in this case. The increase in mouth drying with increased heating time suggests that protein denaturation is a contributing factor and a possible mucoadhesive mechanism is discussed.

The impact of two coupled cirrus microphysics-radiation parameterizations on the temperature and specific humidity biases in the tropical tropopause layer in a climate model

The impact of two different coupled cirrus microphysics-radiation parameterizations on the zonally averaged temperature and humidity biases in the tropical tropopause layer (TTL) of a Met Office climate model configuration is assessed. One parameterization is based on a linear coupling between a model prognostic variable, the ice mass mixing ratio, qi, and the integral optical properties. The second is based on the integral optical properties being parameterized as functions of qi and temperature, Tc, where the mass coefficients (i.e. scattering and extinction) are parameterized as nonlinear functions of the ratio between qi and Tc. The cirrus microphysics parameterization is based on a moment estimation parameterization of the particle size distribution (PSD), which relates the mass moment (i.e. second moment if mass is proportional to size raised to the power of 2 ) of the PSD to all other PSD moments through the magnitude of the second moment and Tc. This same microphysics PSD parameterization is applied to calculate the integral optical properties used in both radiation parameterizations and, thus, ensures PSD and mass consistency between the cirrus microphysics and radiation schemes. In this paper, the temperature-non-dependent and temperature-dependent parameterizations are shown to increase and decrease the zonally averaged temperature biases in the TTL by about 1 K, respectively. The temperature-dependent radiation parameterization is further demonstrated to have a positive impact on the specific humidity biases in the TTL, as well as decreasing the shortwave and longwave biases in the cloudy radiative effect. The temperature-dependent radiation parameterization is shown to be more consistent with TTL and global radiation observations.