Evaluation of modern and mid-Holocene seasonal precipitation of the Mediterranean and northern Africa in the CMIP5 simulations

We analyse the spatial expression of seasonal climates of the Mediterranean and northern Africa in pre-industrial (piControl) and mid-Holocene (midHolocene, 6 yr BP) simulations from the fifth phase of the Coupled Model Intercomparison Project (CMIP5). Modern observations show four distinct precipitation regimes characterized by differences in the seasonal distribution and total amount of precipitation: an equatorial band characterized by a double peak in rainfall, the monsoon zone characterized by summer rainfall, the desert characterized by low seasonality and total precipitation, and the Mediterranean zone characterized by summer drought. Most models correctly simulate the position of the Mediterranean and the equatorial climates in the piControl simulations, but overestimate the extent of monsoon influence and underestimate the extent of desert. However, most models fail to reproduce the amount of precipitation in each zone. Model biases in the simulated magnitude of precipitation are unrelated to whether the models reproduce the correct spatial patterns of each regime. In the midHolocene, the models simulate a reduction in winter rainfall in the equatorial zone, and a northward expansion of the monsoon with a significant increase in summer and autumn rainfall. Precipitation is slightly increased in the desert, mainly in summer and autumn, with northward expansion of the monsoon. Changes in the Mediterranean are small, although there is an increase in spring precipitation consistent with palaeo-observations of increased growing-season rainfall. Comparison with reconstructions shows most models underestimate the mid-Holocene changes in annual precipitation, except in the equatorial zone. Biases in the piControl have only a limited influence on midHolocene anomalies in ocean–atmosphere models; carbon-cycle models show no relationship between piControl bias and midHolocene anomalies. Biases in the prediction of the midHolocene monsoon expansion are unrelated to how well the models simulate changes in Mediterranean climate.

A new data set of soil mineralogy for dust-cycle modeling

The mineralogy of airborne dust affects the impact of dust particles on direct and indirect radiative forcing, on atmospheric chemistry and on biogeochemical cycling. It is determined partly by the mineralogy of the dust-source regions and partly by size-dependent fractionation during erosion and transport. Here we present a data set that characterizes the clay and silt-sized fractions of global soil units in terms of the abundance of 12 minerals that are important for dust–climate interactions: quartz, feldspars, illite, smectite, kaolinite, chlorite, vermiculite, mica, calcite, gypsum, hematite and goethite. The basic mineralogical information is derived from the literature, and is then expanded following explicit rules, in order to characterize as many soil units as possible. We present three alternative realizations of the mineralogical maps, taking the uncertainties in the mineralogical data into account. We examine the implications of the new database for calculations of the single scattering albedo of airborne dust and thus for dust radiative forcing.

Attribution of climate extreme events

There is a tremendous desire to attribute causes to weather and climate events that is often challenging from a physical standpoint. Headlines attributing an event solely to either human-induced climate change or natural variability can be misleading when both are invariably in play. The conventional attribution framework struggles with dynamically driven extremes because of the small signal-to-noise ratios and often uncertain nature of the forced changes. Here, we suggest that a different framing is desirable, which asks why such extremes unfold the way they do. Specifically, we suggest that it is more useful to regard the extreme circulation regime or weather event as being largely unaffected by climate change, and question whether known changes in the climate system's thermodynamic state affected the impact of the particular event. Some examples briefly illustrated include 'snowmaggedon' in February 2010, superstorm Sandy in October 2012 and supertyphoon Haiyan in November 2013, and, in more detail, the Boulder floods of September 2013, all of which were influenced by high sea surface temperatures that had a discernible human component.

Root elongation rate is correlated with the length of the bare root apex of maize and lupin roots despite contrasting responses of root growth to compact and dry soils

Background We investigated interacting effects of matric potential and soil strength on root elongation of maize and lupin, and relations between root elongation rates and the length of bare (hairless) root apex. Methods Root elongation rates and the length of bare root apexwere determined formaize and lupin seedlings in sandy loam soil of various matric potentials (−0.01 to −1.6 MPa) and bulk densities (0.9 to 1.5 Mg m−3). Results Root elongation rates slowed with both decreasing matric potential and increasing penetrometer resistance. Root elongation of maize slowed to 10 % of the unimpeded rate when penetrometer resistance increased to 2 MPa, whereas lupin elongated at about 40 % of the unimpeded rate. Maize root elongation rate was more sensitive to changes in matric potential in loosely packed soil (penetrometer resistances <1 MPa) than lupin. Despite these differing responses, root elongation rate of both species was linearly correlated with length of the bare root apex (r2 0.69 to 0.97). Conclusion Maize root elongation was more sensitive to changes in matric potential and mechanical impedance than lupin. Robust linear relationships between elongation rate and length of bare apex suggest good potential for estimating root elongation rates for excavated roots.

Contrasting responses to drought stress in herbaceous perennial legumes

Background and aims Medicago sativa L. is widely grown in southern Australia, but is poorly adapted to dry, hot summers. This study aimed to identify perennial herbaceous legumes with greater resistance to drought stress and explore their adaptive strategies. Methods Ten herbaceous perennial legume species/accessions were grown in deep pots in a sandy, low-phosphorus field soil in a glasshouse. Drought stress was imposed by ceasing to water. A companion M. sativa plant in each pot minimised differences in leaf area and water consumption among species. Plants were harvested when stomatal conductance of stressed plants decreased to around 10% of well watered plants. Results A range of responses to drought stress were identified, including: reduced shoot growth; leaf curling; thicker pubescence on leaves and stems; an increased root:shoot ratio; an increase, decrease or no change in root distribution with depth; reductions in specific leaf area or leaf water potential; and osmotic adjustment. The suite of changes differed substantially among species and, less so, among accessions. Conclusions The inter- and intra-specific variability of responses to drought-stress in the plants examined suggests a wide range of strategies are available in perennial legumes to cope with drying conditions, and these could be harnessed in breeding/selection programs.

Moisture can be the dominant environmental parameter governing cadaver decomposition in soil

Forensic taphonomy involves the use of decomposition to estimate postmortem interval (PMI) or locate clandestine graves. Yet, cadaver decomposition remains poorly understood, particularly following burial in soil. Presently, we do not know how most edaphic and environmental parameters, including soil moisture, influence the breakdown of cadavers following burial and alter the processes that are used to estimate PMI and locate clandestine graves. To address this, we buried juvenile rat (Rattus rattus) cadavers (∼18 g wet weight) in three contrasting soils from tropical savanna ecosystems located in Pallarenda (sand), Wambiana (medium clay), or Yabulu (loamy sand), Queensland, Australia. These soils were sieved (2 mm), weighed (500 g dry weight), calibrated to a matric potential of -0.01 megapascals (MPa), -0.05 MPa, or -0.3 MPa (wettest to driest) and incubated at 22 °C. Measurements of cadaver decomposition included cadaver mass loss, carbon dioxide-carbon (CO2-C) evolution, microbial biomass carbon (MBC), protease activity, phosphodiesterase activity, ninhydrin-reactive nitrogen (NRN) and soil pH. Cadaver burial resulted in a significant increase in CO2-C evolution, MBC, enzyme activities, NRN and soil pH. Cadaver decomposition in loamy sand and sandy soil was greater at lower matric potentials (wetter soil). However, optimal matric potential for cadaver decomposition in medium clay was exceeded, which resulted in a slower rate of cadaver decomposition in the wettest soil. Slower cadaver decomposition was also observed at high matric potential (-0.3 MPa). Furthermore, wet sandy soil was associated with greater cadaver decomposition than wet fine-textured soil. We conclude that gravesoil moisture content can modify the relationship between temperature and cadaver decomposition and that soil microorganisms can play a significant role in cadaver breakdown. We also conclude that soil NRN is a more reliable indicator of gravesoil than soil pH.

Sequential hydrocarbon biodegradation of a soil from arid coastal Australia oil-treated under laboratory controlled conditions

A general consistency in the sequential order of petroleum hydrocarbon reduction in previous biodegradation studies has led to the proposal of several molecularly based biodegradation scales. Few studies have investigated the biodegradation susceptibility of petroleum hydrocarbon products in soil media, however, and metabolic preferences can change with habitat type. A laboratory based study comprising gas chromatography–mass spectrometry (GC–MS) analysis of extracts of oil-treated soil samples incubated for up to 161 days was conducted to investigate the biodegradation of crude oil exposed to sandy soils of Barrow Island, home to both a Class ‘‘A” nature reserve and Australia’s largest on-shore oil field. Biodegradation trends of the hydrocarbon-treated soils were largely consistent with previous reports but some unusual behaviour was recognised both between and within hydrocarbon classes. For example, the n-alkanes persisted at trace levels from day 86 to 161 following the removal of typically more stable dimethyl naphthalenes and methyl phenanthrenes. The relative susceptibility to biodegradation of different di- tri- and tetramethylnaphthalene isomers also showed several features distinct from previous reports. The unique biodegradation behaviour of Barrow Is. soil likely reflects difference in microbial functioning with physiochemical variation in the environment. Correlation of molecular parameters, reduction rates of selected alkyl naphthalene isomers and CO2 respiration values with a delayed (61 d) oil-treated soil identified a slowing of biodegradation with microcosm incubation; a reduced function or population of incubated soil flora might also influence the biodegradation patterns observed.

Responses of leaf traits to climatic gradients: adaptive variation versus compositional shifts

Dynamic global vegetation models (DGVMs) typically rely on plant functional types (PFTs), which are assigned distinct environmental tolerances and replace one another progressively along environmental gradients. Fixed values of traits are assigned to each PFT; modelled trait variation along gradients is thus driven by PFT replacement. But empirical studies have revealed "universal" scaling relationships (quantitative trait variations with climate that are similar within and between species, PFTs and communities); and continuous, adaptive trait variation has been proposed to replace PFTs as the basis for next-generation DGVMs. Here we analyse quantitative leaf-trait variation on long temperature and moisture gradients in China with a view to understanding the relative importance of PFT replacement vs. continuous adaptive variation within PFTs. Leaf area (LA), specific leaf area (SLA), leaf dry matter content (LDMC) and nitrogen content of dry matter were measured on all species at 80 sites ranging from temperate to tropical climates and from dense forests to deserts. Chlorophyll fluorescence traits and carbon, phosphorus and potassium contents were measured at 47 sites. Generalized linear models were used to relate log-transformed trait values to growing-season temperature and moisture indices, with or without PFT identity as a predictor, and to test for differences in trait responses among PFTs. Continuous trait variation was found to be ubiquitous. Responses to moisture availability were generally similar within and between PFTs, but biophysical traits (LA, SLA and LDMC) of forbs and grasses responded differently from woody plants. SLA and LDMC responses to temperature were dominated by the prevalence of evergreen PFTs with thick, dense leaves at the warm end of the gradient. Nutrient (N, P and K) responses to climate gradients were generally similar within all PFTs. Area-based nutrients generally declined with moisture; Narea and Karea declined with temperature, but Parea increased with temperature. Although the adaptive nature of many of these trait-climate relationships is understood qualitatively, a key challenge for modelling is to predict them quantitatively. Models must take into account that community-level responses to climatic gradients can be influenced by shifts in PFT composition, such as the replacement of deciduous by evergreen trees, which may run either parallel or counter to trait variation within PFTs. The importance of PFT shifts varies among traits, being important for biophysical traits but less so for physiological and chemical traits. Finally, models should take account of the diversity of trait values that is found in all sites and PFTs, representing the "pool" of variation that is locally available for the natural adaptation of ecosystem function to environmental change.

Soil phosphorus dynamics and phytoavailability from sewage sludge at different stages in a treatment stream

The effect of different stages of sewage sludge treatment on phosphorus (P) dynamics in amended soils was determined using samples of undigested liquid (UL), anaerobically digested liquid (AD) and dewatered anaerobically digested (DC) sludge. Sludges were taken from three points in the same treatment stream and applied to a sandy loam soil in field-based mesocosms at 4, 8 and 16t ha−1 dry solids. Mesocosms were sown with perennial ryegrass (Lolium perenne cv. Melle), and the sward was harvested after 35 and 70 days to determine yield and foliar P concentration. Soils were also sampled during this period to measure P transformations and the activities of acid phosphomonoesterase and phosphodiesterase. Data show that the AD amended soils had the greatest plant-available and foliar P content up to the second harvest, but the UL amended soils had the greatest enzyme activity. Characterisation of control and 16t ha−1 soils and sludge using solution 31P nuclear magnetic resonance (NMR) spectroscopy after NaOH–EDTA extraction revealed that P was predominantly in the inorganic pool in all three sludge samples, with the highest proportion (of the total extracted P) as inorganic P in the anaerobically digested liquid sludge. After sludge incorporation, P was immobilised to organic species. The majority of organic P was in monoester-P forms, while the remainder of organic P (diester P and phosphonate P) was more susceptible to transformations through time and showed variation with sludge type. These results show that application of sewage sludge at rates as low as 4t ha−1 can have a significant nutritional benefit to ryegrass over an initial 35-day growth and subsequent 35-day re-growth periods. Differences in P transformation, and hence nutritional benefit, between sludge types were evident throughout the experiment. Thus, differences in sludge treatment process alter the edaphic mineralisation characteristics of biosolids derived from the same source material.

Nitrogen dynamics under Lolium perenne after a single application of different sewage sludge types from the same treatment stream

Three sludge types from the same treatment stream (undigested liquid, anaerobically digested liquid and dewatered, anaerobically digested cake) were used in a field based tub study. Amendments (4, 8, and 16 Mg dry solid (ds)ha(-1)) were incorporated into the upper 15 cm of a sandy loam soil prior to sowing with rye-grass (Lolium perenne L.). Nitrogen transformations in the soil were determined for the 80 d period following incorporation. Nitrogen uptake and crop yield were measured in the cut sward 35 and 70 d after sowing. The study showed that application of sewage sludge at rates as low as 4 Mgha(-1) can have a nutritional benefit to rye-grass over the two harvests. Differences in N transformation, and hence crop nutritional benefit, between sludge types were evident throughout the experiment. In particular, the dewatering process changed the mineral N characteristics of the anaerobically digested sludge, which, when not dewatered, outperformed the other sludges in terms of yield and mineralisation rate at both harvests. The dewatered sludge produced the lowest yield of rye-grass. The undigested liquid sludge had the lowest foliar N and soil NO(3)-N concentrations, possibly immobilised as the large oxidisable C component of this sludge was metabolised by the microbial biomass. Correlation data support the concept of preferential uptake of NH(4)-N over NO(3)-N in Lolium perenne. Results are discussed in the context of managing sludge type and application for a plant nutrient source and NO(3)-N release.

A laboratory incubation method for determining the rate of microbiological degradation of skeletal muscle tissue in soil

A controlled laboratory experiment is described, in principle and practice, which can be used for the of determination the rate of tissue decomposition in soil. By way of example, an experiment was conducted to determine the effect of temperature (12°C, 22°C) on the aerobic decomposition of skeletal muscle tissue (Organic Texel × Suffolk lamb (Ovis aries)) in a sandy loam soil. Measurements of decomposition processes included muscle tissue mass loss, microbial CO2 respiration, and muscle tissue carbon (C) and nitrogen (N). Muscle tissue mass loss at 22°C always was greater than at 12°C (p < 0.001). Microbial respiration was greater in samples incubated at 22°C for the initial 21 days of burial (p < 0.01). All buried muscle tissue samples demonstrated changes in C and N content at the end of the experiment. A significant correlation (p < 0.001) was demonstrated between the loss of muscle tissue-derived C (C1) and microbially-respired C (Cm) demonstrating CO2 respiration may be used to predict mass loss and hence biodegradation. In this experiment Q10 (12°C - 22°C) = 2.0. This method is recommended as a useful tool in determining the effect of environmental variables on the rate of decomposition of various tissues and associated materials.

Climate-driven expansion of blanket bogs in Britain during the Holocene

Blanket bog occupies approximately 6 % of the area of the UK today. The Holocene expansion of this hyperoceanic biome has previously been explained as a consequence of Neolithic forest clearance. However, the present distribution of blanket bog in Great Britain can be predicted accurately with a simple model (PeatStash) based on summer temperature and moisture index thresholds, and the same model correctly predicts the highly disjunct distribution of blanket bog worldwide. This finding suggests that climate, rather than land-use history, controls blanket-bog distribution in the UK and everywhere else. We set out to test this hypothesis for blanket bogs in the UK using bioclimate envelope modelling compared with a database of peat initiation age estimates. We used both pollen-based reconstructions and climate model simulations of climate changes between the mid-Holocene (6000 yr BP, 6 ka) and modern climate to drive PeatStash and predict areas of blanket bog. We compiled data on the timing of blanket-bog initiation, based on 228 age determinations at sites where peat directly overlies mineral soil. The model predicts large areas of northern Britain would have had blanket bog by 6000 yr BP, and the area suitable for peat growth extended to the south after this time. A similar pattern is shown by the basal peat ages and new blanket bog appeared over a larger area during the late Holocene, the greatest expansion being in Ireland, Wales and southwest England, as the model predicts. The expansion was driven by a summer cooling of about 2 °C, shown by both pollen-based reconstructions and climate models. The data show early Holocene (pre-Neolithic) blanket-bog initiation at over half of the sites in the core areas of Scotland, and northern England. The temporal patterns and concurrence of the bioclimate model predictions and initiation data suggest that climate change provides a parsimonious explanation for the early Holocene distribution and later expansion of blanket bogs in the UK, and it is not necessary to invoke anthropogenic activity as a driver of this major landscape change.

Ice-sheet configuration in the CMIP5/PMIP3 Last Glacial Maximum experiments

We describe the creation of a data set describing changes related to the presence of ice sheets, including ice-sheet extent and height, ice-shelf extent, and the distribution and elevation of ice-free land at the Last Glacial Maximum (LGM), which were used in LGM experiments conducted as part of the fifth phase of the Coupled Modelling Intercomparison Project (CMIP5) and the third phase of the Palaeoclimate Modelling Intercomparison Project (PMIP3). The CMIP5/PMIP3 data sets were created from reconstructions made by three different groups, which were all obtained using a model-inversion approach but differ in the assumptions used in the modelling and in the type of data used as constraints. The ice-sheet extent in the Northern Hemisphere (NH) does not vary substantially between the three individual data sources. The difference in the topography of the NH ice sheets is also moderate, and smaller than the differences between these reconstructions (and the resultant composite reconstruction) and ice-sheet reconstructions used in previous generations of PMIP. Only two of the individual reconstructions provide information for Antarctica. The discrepancy between these two reconstructions is larger than the difference for the NH ice sheets, although still less than the difference between the composite reconstruction and previous PMIP ice-sheet reconstructions. Although largely confined to the ice-covered regions, differences between the climate response to the individual LGM reconstructions extend over the North Atlantic Ocean and Northern Hemisphere continents, partly through atmospheric stationary waves. Differences between the climate response to the CMIP5/PMIP3 composite and any individual ice-sheet reconstruction are smaller than those between the CMIP5/PMIP3 composite and the ice sheet used in the last phase of PMIP (PMIP2).

Evaluation of the realism of climate reconstruction using the Coexistence Approach with modern pollen samples from the Tibetan Plateau

We apply the Coexistence Approach (CoA) to reconstruct mean annual precipitation (MAP), mean annual temperature (MAT), mean temperature of thewarmestmonth (MTWA) and mean temperature of the coldest month (MTCO) at 44 pollen sites on the Qinghai–Tibetan Plateau. The modern climate ranges of the taxa are obtained (1) from county-level presence/absence data and (2) from data on the optimum and range of each taxon from Lu et al. (2011). The CoA based on the optimumand range data yields better predictions of observed climate parameters at the pollen sites than that based on the county-level data. The presence of arboreal pollen, most of which is derived fromoutside the region, distorts the reconstructions. More reliable reconstructions are obtained using only the non-arboreal component of the pollen assemblages. The root mean-squared error (RMSE) of the MAP reconstructions are smaller than the RMSE of MAT, MTWA and MTCO, suggesting that precipitation gradients are the most important control of vegetation distribution on the Qinghai–Tibetan Plateau. Our results show that CoA could be used to reconstruct past climates in this region, although in areas characterized by open vegetation the most reliable estimates will be obtained by excluding possible arboreal contaminants.

Evaluation of CMIP5 palaeo-simulations to improve climate projections

Structural differences among models account for much of the uncertainty in projected climate changes, at least until the mid-twenty-first century. Recent observations encompass too limited a range of climate variability to provide a robust test of the ability to simulate climate changes. Past climate changes provide a unique opportunity for out-of-sample evaluation of model performance. Palaeo-evaluation has shown that the large-scale changes seen in twenty-first-century projections, including enhanced land–sea temperature contrast, latitudinal amplification, changes in temperature seasonality and scaling of precipitation with temperature, are likely to be realistic. Although models generally simulate changes in large-scale circulation sufficiently well to shift regional climates in the right direction, they often do not predict the correct magnitude of these changes. Differences in performance are only weakly related to modern-day biases or climate sensitivity, and more sophisticated models are not better at simulating climate changes. Although models correctly capture the broad patterns of climate change, improvements are required to produce reliable regional projections.