Impact of vegetation and preferential source areas on global dust aerosol: results from a model study

[1] We present a model of the dust cycle that successfully predicts dust emissions as determined by land surface properties, monthly vegetation and snow cover, and 6-hourly surface wind speeds for the years 1982–1993. The model takes account of the role of dry lake beds as preferential source areas for dust emission. The occurrence of these preferential sources is determined by a water routing and storage model. The dust source scheme also explicitly takes into account the role of vegetation type as well as monthly vegetation cover. Dust transport is computed using assimilated winds for the years 1987–1990. Deposition of dust occurs through dry and wet deposition, where subcloud scavenging is calculated using assimilated precipitation fields. Comparison of simulated patterns of atmospheric dust loading with the Total Ozone Mapping Spectrometer satellite absorbing aerosol index shows that the model produces realistic results from daily to interannual timescales. The magnitude of dust deposition agrees well with sediment flux data from marine sites. Emission of submicron dust from preferential source areas are required for the computation of a realistic dust optical thickness. Sensitivity studies show that Asian dust source strengths are particularly sensitive to the seasonality of vegetation cover.

Observation impact in data assimilation: the effect of non-Gaussian observation error.

Data assimilation methods which avoid the assumption of Gaussian error statistics are being developed for geoscience applications. We investigate how the relaxation of the Gaussian assumption affects the impact observations have within the assimilation process. The effect of non-Gaussian observation error (described by the likelihood) is compared to previously published work studying the effect of a non-Gaussian prior. The observation impact is measured in three ways: the sensitivity of the analysis to the observations, the mutual information, and the relative entropy. These three measures have all been studied in the case of Gaussian data assimilation and, in this case, have a known analytical form. It is shown that the analysis sensitivity can also be derived analytically when at least one of the prior or likelihood is Gaussian. This derivation shows an interesting asymmetry in the relationship between analysis sensitivity and analysis error covariance when the two different sources of non-Gaussian structure are considered (likelihood vs. prior). This is illustrated for a simple scalar case and used to infer the effect of the non-Gaussian structure on mutual information and relative entropy, which are more natural choices of metric in non-Gaussian data assimilation. It is concluded that approximating non-Gaussian error distributions as Gaussian can give significantly erroneous estimates of observation impact. The degree of the error depends not only on the nature of the non-Gaussian structure, but also on the metric used to measure the observation impact and the source of the non-Gaussian structure.

Effect of uncertainty in surface mass balance--elevation feedback on projections of the future sea level contribution of the Greenland ice sheet

We apply a new parameterisation of the Greenland ice sheet (GrIS) feedback between surface mass balance (SMB: the sum of surface accumulation and surface ablation) and surface elevation in the MAR regional climate model (Edwards et al., 2014) to projections of future climate change using five ice sheet models (ISMs). The MAR (Modèle Atmosphérique Régional: Fettweis, 2007) climate projections are for 2000–2199, forced by the ECHAM5 and HadCM3 global climate models (GCMs) under the SRES A1B emissions scenario. The additional sea level contribution due to the SMB– elevation feedback averaged over five ISM projections for ECHAM5 and three for HadCM3 is 4.3% (best estimate; 95% credibility interval 1.8–6.9 %) at 2100, and 9.6% (best estimate; 95% credibility interval 3.6–16.0 %) at 2200. In all results the elevation feedback is significantly positive, amplifying the GrIS sea level contribution relative to the MAR projections in which the ice sheet topography is fixed: the lower bounds of our 95% credibility intervals (CIs) for sea level contributions are larger than the “no feedback” case for all ISMs and GCMs. Our method is novel in sea level projections because we propagate three types of modelling uncertainty – GCM and ISM structural uncertainties, and elevation feedback parameterisation uncertainty – along the causal chain, from SRES scenario to sea level, within a coherent experimental design and statistical framework. The relative contributions to uncertainty depend on the timescale of interest. At 2100, the GCM uncertainty is largest, but by 2200 both the ISM and parameterisation uncertainties are larger. We also perform a perturbed parameter ensemble with one ISM to estimate the shape of the projected sea level probability distribution; our results indicate that the probability density is slightly skewed towards higher sea level contributions.

Regional and monthly and clear-sky aerosol direct radiative effect (and forcing) derived from the GlobAEROSOL-AATSR satellite aerosol product

Using the GlobAEROSOL-AATSR dataset, estimates of the instantaneous, clear-sky, direct aerosol radiative effect and radiative forcing have been produced for the year 2006. Aerosol Robotic Network sun-photometer measurements have been used to characterise the random and systematic error in the GlobAEROSOL product for 22 regions covering the globe. Representative aerosol properties for each region were derived from the results of a wide range of literature sources and, along with the de-biased GlobAEROSOL AODs, were used to drive an offline version of the Met Office unified model radiation scheme. In addition to the mean AOD, best-estimate run of the radiation scheme, a range of additional calculations were done to propagate uncertainty estimates in the AOD, optical properties, surface albedo and errors due to the temporal and spatial averaging of the AOD fields. This analysis produced monthly, regional estimates of the clear-sky aerosol radiative effect and its uncertainty, which were combined to produce annual, global mean values of (−6.7±3.9)Wm−2 at the top of atmosphere (TOA) and (−12±6)Wm−2 at the surface. These results were then used to give estimates of regional, clear-sky aerosol direct radiative forcing, using modelled pre-industrial AOD fields for the year 1750 calculated for the AEROCOM PRE experiment. However, as it was not possible to quantify the uncertainty in the pre-industrial aerosol loading, these figures can only be taken as indicative and their uncertainties as lower bounds on the likely errors. Although the uncertainty on aerosol radiative effect presented here is considerably larger than most previous estimates, the explicit inclusion of the major sources of error in the calculations suggest that they are closer to the true constraint on this figure from similar methodologies, and point to the need for more, improved estimates of both global aerosol loading and aerosol optical properties.

Impacts of increasing the aerosol complexity in the Met Office global numerical weather prediction model.

The inclusion of the direct and indirect radiative effects of aerosols in high-resolution global numerical weather prediction (NWP) models is being increasingly recognised as important for the improved accuracy of short-range weather forecasts. In this study the impacts of increasing the aerosol complexity in the global NWP configuration of the Met Office Unified Model (MetUM) are investigated. A hierarchy of aerosol representations are evaluated including three-dimensional monthly mean speciated aerosol climatologies, fully prognostic aerosols modelled using the CLASSIC aerosol scheme and finally, initialised aerosols using assimilated aerosol fields from the GEMS project. The prognostic aerosol schemes are better able to predict the temporal and spatial variation of atmospheric aerosol optical depth, which is particularly important in cases of large sporadic aerosol events such as large dust storms or forest fires. Including the direct effect of aerosols improves model biases in outgoing long-wave radiation over West Africa due to a better representation of dust. However, uncertainties in dust optical properties propagate to its direct effect and the subsequent model response. Inclusion of the indirect aerosol effects improves surface radiation biases at the North Slope of Alaska ARM site due to lower cloud amounts in high-latitude clean-air regions. This leads to improved temperature and height forecasts in this region. Impacts on the global mean model precipitation and large-scale circulation fields were found to be generally small in the short-range forecasts. However, the indirect aerosol effect leads to a strengthening of the low-level monsoon flow over the Arabian Sea and Bay of Bengal and an increase in precipitation over Southeast Asia. Regional impacts on the African Easterly Jet (AEJ) are also presented with the large dust loading in the aerosol climatology enhancing of the heat low over West Africa and weakening the AEJ. This study highlights the importance of including a more realistic treatment of aerosol–cloud interactions in global NWP models and the potential for improved global environmental prediction systems through the incorporation of more complex aerosol schemes.

Aerosol direct radiative effect of smoke over clouds over the southeast Atlantic Ocean from 2006 to 2009

The aerosol direct radiative effect (DRE) of African smoke was analyzed in cloud scenes over the southeast Atlantic Ocean, using Scanning Imaging Absorption Spectrometer for Atmospheric Chartography (SCIAMACHY) satellite observations and Hadley Centre Global Environmental Model version 2 (HadGEM2) climate model simulations. The observed mean DRE was about 30–35 W m−2 in August and September 2006–2009. In some years, short episodes of high-aerosol DRE can be observed, due to high-aerosol loadings, while in other years the loadings are lower but more prolonged. Climate models that use evenly distributed monthly averaged emission fields will not reproduce these high-aerosol loadings. Furthermore, the simulated monthly mean aerosol DRE in HadGEM2 is only about 6 W m−2 in August. The difference with SCIAMACHY mean observations can be partly explained by an underestimation of the aerosol absorption Ångström exponent in the ultraviolet. However, the subsequent increase of aerosol DRE simulation by about 20% is not enough to explain the observed discrepancy between simulations and observations.

The effect of alkyl substitution on the extraction properties of BTPhen ligands for the partitioning of trivalent minor actinides and lanthanides in spent nuclear fuels

Bis-triazinylphenanthroline ligands (BTPhens), which contain additional alkyl (n-butyl and sec-butyl) groups attached to the triazine rings, have been synthesized, and the effects of this alkyl substitution on their extraction properties with Ln(III) and An(III) cations in simulated nuclear waste solutions have been studied. The speciation of n-butyl-substituted ligand (C4- BTPhen) with some trivalent lanthanide nitrates was elucidated by 1 H-NMR spectroscopic titrations. These experiments have shown that the dominant species in solution were the 1:2 complexes [Ln(III)(BTPhen)2], even at higher Ln(III) concentrations, and the relative stability of 2:1 to 1:1 BTPhen-Ln(III) complexes varied with different lanthanides. As expected, sec-butylsubstituted ligand (sec-C4 BTPhen) showed higher solubility than C4-BTPhen in certain diluents. A greater separation factor (SFAm/Eu = ca. 210) was observed for sec-C4-BTPhen compared to C4-BTPhen (SFAm/Eu = ca. 125) in 1-octanol at 4 M HNO3 solutions. The greater separation factor may be due to the higher solubility of the 2:1 complex for sec-C4-BTPhen at the interface than the 1:1 complex of C4-BTPhen.

Quantifying sources of inter-model diversity in the cloud albedo effect

There is large diversity in simulated aerosol forcing among models that participated in the fifth Coupled Model Intercomparison Project (CMIP5), particularly related to aerosol interactions with clouds. Here we use the reported model data and fitted aerosol-cloud relations to separate the main sources of inter-model diversity in the magnitude of the cloud albedo effect. There is large diversity in the global load and spatial distribution of sulfate aerosol, as well as in global-mean cloud-top effective radius. The use of different parameterizations of aerosol-cloud interactions makes the largest contribution to diversity in modeled radiative forcing (up to -39%, +48% about the mean estimate). Uncertainty in pre-industrial sulfate load also makes a substantial contribution (-15%, +61% about the mean estimate), with smaller contributions from inter-model differences in the historical change in sulfate load and in mean cloud fraction.

Interactive effect of cowpea variety, dose and exposure time on bruchid tolerance to botanical pesticides

Callosobruchus maculatus has for years remained a serious menace in cowpea in Sub-Sahara Africa. The objective of this study was to investigate the effect of genotypic cowpea (Vigna unguiculata (L.) Walp) varieties, time and dose on C. maculatus exposed to powders of Piper guineense and Eugenia aromatica. Irrespective of duration and botanicals, bruchid reared on KDV showed the highest tolerance to both plant materials; while their counterparts from IAR48V were the most susceptible. Median lethal time (LT50) also varied according to the plant materials; with the highest in KDV reared bruchid [P. guineense: KDV (18.31), IAR48V (9.27), IFBV (13.17); E. aromatica: KDV (76.01), IAR48V (5.59), IFBV (6.49)]. There was a significant impact of cowpea variety (V), exposure time (T) and dose (D) on the tolerance of C. maculatus to both plant materials. The effect of all two-way (VxT, VxD, DxT) and three way interactions (V×T×D) on the tolerance of C. maculatus to both plant materials was also significant. Varietal effect was more pronounced in bruchids exposed to E. aromatica; while exposure time was more pronounced in bruchids exposed to P. guineense.

Effect of high temperature and water stress on groundnuts under field conditions

Groundnuts cultivated in the semiarid tropics are often exposed to water stress (mid-season and end season) and high temperature (> 34 °C) during the critical stages of flowering and pod development. This study evaluated the effects of both water stress and high temperature under field conditions at ICRISAT, India. Treatments included two irrigations (full irrigation, 100 % of crop evapotranspiration; and water stress, 40 % of crop evapotranspiration), four temperature treatments from a combination of two sowing dates and heat tunnels with mean temperatures from sowing to maturity of 26.3° (T1), 27.3° (T2), 29.0° (T3) and 29.7 °C (T4) and two genotypes TMV2 and ICGS 11. The heat tunnels were capable of raising the day temperature by > 10 °C compared to ambient. During the 20-day high-temperature treatment at flowering, mean temperatures were 33.8° (T1), 41.6° (T2), 38.7° (T3) and 43.5°C (T4). The effects of water stress and high temperature were additive and temporary for both vegetative and pod yield, and disappeared as soon as high-temperature stress was removed. Water use efficiency was significantly affected by the main effects of temperature and cultivar and not by water stress treatments. Genotypic differences for tolerance to high temperature can be attributed to differences in flowering pattern, flower number, peg-set and harvest index. It can be inferred from this study that genotypes that are tolerant to water stress are also tolerant to high temperature under field conditions. In addition, genotypes with an ability to establish greater biomass and with a significantly greater partitioning of biomass to pod yield would be suitable for sustaining higher yields in semiarid tropics with high temperature and water stress.

Efficient incorporation of channel cross-section geometry uncertainty into regional and global scale flood inundation models

This paper investigates the challenge of representing structural differences in river channel cross-section geometry for regional to global scale river hydraulic models and the effect this can have on simulations of wave dynamics. Classically, channel geometry is defined using data, yet at larger scales the necessary information and model structures do not exist to take this approach. We therefore propose a fundamentally different approach where the structural uncertainty in channel geometry is represented using a simple parameterization, which could then be estimated through calibration or data assimilation. This paper first outlines the development of a computationally efficient numerical scheme to represent generalised channel shapes using a single parameter, which is then validated using a simple straight channel test case and shown to predict wetted perimeter to within 2% for the channels tested. An application to the River Severn, UK is also presented, along with an analysis of model sensitivity to channel shape, depth and friction. The channel shape parameter was shown to improve model simulations of river level, particularly for more physically plausible channel roughness and depth parameter ranges. Calibrating channel Manning’s coefficient in a rectangular channel provided similar water level simulation accuracy in terms of Nash-Sutcliffe efficiency to a model where friction and shape or depth were calibrated. However, the calibrated Manning coefficient in the rectangular channel model was ~2/3 greater than the likely physically realistic value for this reach and this erroneously slowed wave propagation times through the reach by several hours. Therefore, for large scale models applied in data sparse areas, calibrating channel depth and/or shape may be preferable to assuming a rectangular geometry and calibrating friction alone.

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.

Future changes in the Western North Pacific tropical cyclone activity projected by a multidecadal simulation with a 16-km global atmospheric GCM

How tropical cyclone (TC) activity in the northwestern Pacific might change in a future climate is assessed using multidecadal Atmospheric Model Intercomparison Project (AMIP)-style and time-slice simulations with the ECMWF Integrated Forecast System (IFS) at 16-km and 125-km global resolution. Both models reproduce many aspects of the present-day TC climatology and variability well, although the 16-km IFS is far more skillful in simulating the full intensity distribution and genesis locations, including their changes in response to El Niño–Southern Oscillation. Both IFS models project a small change in TC frequency at the end of the twenty-first century related to distinct shifts in genesis locations. In the 16-km IFS, this shift is southward and is likely driven by the southeastward penetration of the monsoon trough/subtropical high circulation system and the southward shift in activity of the synoptic-scale tropical disturbances in response to the strengthening of deep convective activity over the central equatorial Pacific in a future climate. The 16-km IFS also projects about a 50% increase in the power dissipation index, mainly due to significant increases in the frequency of the more intense storms, which is comparable to the natural variability in the model. Based on composite analysis of large samples of supertyphoons, both the development rate and the peak intensities of these storms increase in a future climate, which is consistent with their tendency to develop more to the south, within an environment that is thermodynamically more favorable for faster development and higher intensities. Coherent changes in the vertical structure of supertyphoon composites show system-scale amplification of the primary and secondary circulations with signs of contraction, a deeper warm core, and an upward shift in the outflow layer and the frequency of the most intense updrafts. Considering the large differences in the projections of TC intensity change between the 16-km and 125-km IFS, this study further emphasizes the need for high-resolution modeling in assessing potential changes in TC activity.

Metrics for linking emissions of gases and aerosols to global precipitation changes

Recent advances in understanding have made it possible to relate global precipitation changes directly to emissions of particular gases and aerosols that influence climate. Using these advances, new indices are developed here called the Global Precipitation-change Potential for pulse (GPP_P) and sustained (GPP_S) emissions, which measure the precipitation change per unit mass of emissions. The GPP can be used as a metric to compare the effects of different emissions. This is akin to the global warming potential (GWP) and the global temperature-change potential (GTP) which are used to place emissions on a common scale. Hence the GPP provides an additional perspective of the relative or absolute effects of emissions. It is however recognised that precipitation changes are predicted to be highly variable in size and sign between different regions and this limits the usefulness of a purely global metric. The GPP_P and GPP_S formulation consists of two terms, one dependent on the surface temperature change and the other dependent on the atmospheric component of the radiative forcing. For some forcing agents, and notably for CO2, these two terms oppose each other – as the forcing and temperature perturbations have different timescales, even the sign of the absolute GPP_P and GPP_S varies with time, and the opposing terms can make values sensitive to uncertainties in input parameters. This makes the choice of CO2 as a reference gas problematic, especially for the GPP_S at time horizons less than about 60 years. In addition, few studies have presented results for the surface/atmosphere partitioning of different forcings, leading to more uncertainty in quantifying the GPP than the GWP or GTP. Values of the GPP_P and GPP_S for five long- and short-lived forcing agents (CO2, CH4, N2O, sulphate and black carbon – BC) are presented, using illustrative values of required parameters. The resulting precipitation changes are given as the change at a specific time horizon (and hence they are end-point metrics) but it is noted that the GPPS can also be interpreted as the time-integrated effect of a pulse emission. Using CO2 as a references gas, the GPP_P and GPP_S for the non-CO2 species are larger than the corresponding GTP values. For BC emissions, the atmospheric forcing is sufficiently strong that the GPP_S is opposite in sign to the GTP_S. The sensitivity of these values to a number of input parameters is explored. The GPP can also be used to evaluate the contribution of different emissions to precipitation change during or after a period of emissions. As an illustration, the precipitation changes resulting from emissions in 2008 (using the GPP_P) and emissions sustained at 2008 levels (using the GPP_S) are presented. These indicate that for periods of 20 years (after the 2008 emissions) and 50 years (for sustained emissions at 2008 levels) methane is the dominant driver of positive precipitation changes due to those emissions. For sustained emissions, the sum of the effect of the five species included here does not become positive until after 50 years, by which time the global surface temperature increase exceeds 1 K.

The atomic structure of low-index surfaces of the intermetallic compound InPd

The intermetallic compound InPd (CsCl type of crystal structure with a broad compositional range) is considered as a candidate catalyst for the steam reforming of methanol. Single crystals of this phase have been grown to study the structure of its three low-index surfaces under ultra-high vacuum conditions, using low energy electron diffraction (LEED), X-ray photoemission spectroscopy (XPS), and scanning tunneling microscopy (STM). During surface preparation, preferential sputtering leads to a depletion of In within the top few layers for all three surfaces. The near-surface regions remain slightly Pd-rich until annealing to ∼580 K. A transition occurs between 580 and 660 K where In segregates towards the surface and the near-surface regions become slightly In-rich above ∼660 K. This transition is accompanied by a sharpening of LEED patterns and formation of flat step-terrace morphology, as observed by STM. Several superstructures have been identified for the different surfaces associated with this process. Annealing to higher temperatures (≥750 K) leads to faceting via thermal etching as shown for the (110) surface, with a bulk In composition close to the In-rich limit of the existence domain of the cubic phase. The Pd-rich InPd(111) is found to be consistent with a Pd-terminated bulk truncation model as shown by dynamical LEED analysis while, after annealing at higher temperature, the In-rich InPd(111) is consistent with an In-terminated bulk truncation, in agreement with density functional theory (DFT) calculations of the relative surface energies. More complex surface structures are observed for the (100) surface. Additionally, individual grains of a polycrystalline sample are characterized by micro-spot XPS and LEED as well as low-energy electron microscopy. Results from both individual grains and “global” measurements are interpreted based on comparison to our single crystals findings, DFT calculations and previous literature.