Simulated last glacial maximum D14Catm and the deep glacial ocean carbon reservoir

∆14Catm has been estimated as 420 ± 80‰ (IntCal09) during the Last Glacial Maximum (LGM) compared to preindustrial times (0‰), but mechanisms explaining this difference are not yet resolved. ∆14Catm is a function of both cosmogenic production in the high atmosphere and of carbon cycling and partitioning in the Earth system. 10Be-based reconstructions show a contribution of the cosmogenic production term of only 200 ± 200‰ in the LGM. The remaining 220‰ have thus to be explained by changes in the carbon cycle. Recently, Bouttes et al. (2010, 2011) proposed to explain most of the difference in pCO2atm and δ13C between glacial and interglacial times as a result of brine-induced ocean stratification in the Southern Ocean. This mechanism involves the formation of very saline water masses that contribute to high carbon storage in the deep ocean. During glacial times, the sinking of brines is enhanced and more carbon is stored in the deep ocean, lowering pCO2atm. Moreover, the sinking of brines induces increased stratification in the Southern Ocean, which keeps the deep ocean well isolated from the surface. Such an isolated ocean reservoir would be characterized by a low ∆14C signature. Evidence of such 14C-depleted deep waters during the LGM has recently been found in the Southern Ocean (Skinner et al. 2010). The degassing of this carbon with low ∆14C would then reduce ∆14Catm throughout the deglaciation. We have further developed the CLIMBER-2 model to include a cosmogenic production of 14C as well as an interactive atmospheric 14C reservoir. We investigate the role of both the sinking of brine and cosmogenic production, alongside iron fertilization mechanisms, to explain changes in ∆14Catm during the last deglaciation. In our simulations, not only is the sinking of brine mechanism consistent with past ∆14C data, but it also explains most of the differences in pCO2atm and ∆14Catm between the LGM and preindustrial times. Finally, this study represents the first time to our knowledge that a model experiment explains glacial-interglacial differences in pCO2atm, δ13C, and ∆14C together with a coherent LGM climate.

The AeroCom evaluation and intercomparison of organic aerosol in global models

This paper evaluates the current status of global modeling of the organic aerosol (OA) in the troposphere and analyzes the differences between models as well as between models and observations. Thirty-one global chemistry transport models (CTMs) and general circulation models (GCMs) have participated in this intercomparison, in the framework of AeroCom phase II. The simulation of OA varies greatly between models in terms of the magnitude of primary emissions, secondary OA (SOA) formation, the number of OA species used (2 to 62), the complexity of OA parameterizations (gas-particle partitioning, chemical aging, multiphase chemistry, aerosol microphysics), and the OA physical, chemical and optical properties. The diversity of the global OA simulation results has increased since earlier AeroCom experiments, mainly due to the increasing complexity of the SOA parameterization in models, and the implementation of new, highly uncertain, OA sources. Diversity of over one order of magnitude exists in the modeled vertical distribution of OA concentrations that deserves a dedicated future study. Furthermore, although the OA / OC ratio depends on OA sources and atmospheric processing, and is important for model evaluation against OA and OC observations, it is resolved only by a few global models. The median global primary OA (POA) source strength is 56 Tg a−1 (range 34–144 Tg a−1) and the median SOA source strength (natural and anthropogenic) is 19 Tg a−1 (range 13–121 Tg a−1). Among the models that take into account the semi-volatile SOA nature, the median source is calculated to be 51 Tg a−1 (range 16–121 Tg a−1), much larger than the median value of the models that calculate SOA in a more simplistic way (19 Tg a−1; range 13–20 Tg a−1, with one model at 37 Tg a−1). The median atmospheric burden of OA is 1.4 Tg (24 models in the range of 0.6–2.0 Tg and 4 between 2.0 and 3.8 Tg), with a median OA lifetime of 5.4 days (range 3.8–9.6 days). In models that reported both OA and sulfate burdens, the median value of the OA/sulfate burden ratio is calculated to be 0.77; 13 models calculate a ratio lower than 1, and 9 models higher than 1. For 26 models that reported OA deposition fluxes, the median wet removal is 70 Tg a−1 (range 28–209 Tg a−1), which is on average 85% of the total OA deposition. Fine aerosol organic carbon (OC) and OA observations from continuous monitoring networks and individual field campaigns have been used for model evaluation. At urban locations, the model–observation comparison indicates missing knowledge on anthropogenic OA sources, both strength and seasonality. The combined model–measurements analysis suggests the existence of increased OA levels during summer due to biogenic SOA formation over large areas of the USA that can be of the same order of magnitude as the POA, even at urban locations, and contribute to the measured urban seasonal pattern. Global models are able to simulate the high secondary character of OA observed in the atmosphere as a result of SOA formation and POA aging, although the amount of OA present in the atmosphere remains largely underestimated, with a mean normalized bias (MNB) equal to −0.62 (−0.51) based on the comparison against OC (OA) urban data of all models at the surface, −0.15 (+0.51) when compared with remote measurements, and −0.30 for marine locations with OC data. The mean temporal correlations across all stations are low when compared with OC (OA) measurements: 0.47 (0.52) for urban stations, 0.39 (0.37) for remote stations, and 0.25 for marine stations with OC data. The combination of high (negative) MNB and higher correlation at urban stations when compared with the low MNB and lower correlation at remote sites suggests that knowledge about the processes that govern aerosol processing, transport and removal, on top of their sources, is important at the remote stations. There is no clear change in model skill with increasing model complexity with regard to OC or OA mass concentration. However, the complexity is needed in models in order to distinguish between anthropogenic and natural OA as needed for climate mitigation, and to calculate the impact of OA on climate accurately.

Sensitivity of a coupled climate model to canopy interception capacity

The canopy interception capacity is a small but key part of the surface hydrology, which affects the amount of water intercepted by vegetation and therefore the partitioning of evaporation and transpiration. However, little research with climate models has been done to understand the effects of a range of possible canopy interception capacity parameter values. This is in part due to the assumption that it does not significantly affect climate. Near global evapotranspiration products now make evaluation of canopy interception capacity parameterisations possible. We use a range of canopy water interception capacity values from the literature to investigate the effect on climate within the climate model HadCM3. We find that the global mean temperature is affected by up to -0.64 K globally and -1.9 K regionally. These temperature impacts are predominantly due to changes in the evaporative fraction and top of atmosphere albedo. In the tropics, the variations in evapotranspiration affect precipitation, significantly enhancing rainfall. Comparing the model output to measurements, we find that the default canopy interception capacity parameterisation overestimates canopy interception loss (i.e. canopy evaporation) and underestimates transpiration. Overall, decreasing canopy interception capacity improves the evapotranspiration partitioning in HadCM3, though the measurement literature more strongly supports an increase. The high sensitivity of climate to the parameterisation of canopy interception capacity is partially due to the high number of light rain-days in the climate model that means that interception is overestimated. This work highlights the hitherto underestimated importance of canopy interception capacity in climate model hydroclimatology and the need to acknowledge the role of precipitation representation limitations in determining parameterisations.

Infrared and millimetre-wave scintillometry in the suburban environment – Part 1: Structure parameters

Scintillometry, a form of ground-based remote sensing, provides the capability to estimate surface heat fluxes over scales of a few hundred metres to kilometres. Measurements are spatial averages, making this technique particularly valuable over areas with moderate heterogeneity such as mixed agricultural or urban environments. In this study, we present the structure parameters of temperature and humidity, which can be related to the sensible and latent heat fluxes through similarity theory, for a suburban area in the UK. The fluxes are provided in the second paper of this two-part series. A millimetre-wave scintillometer was combined with an infrared scintillometer along a 5.5 km path over northern Swindon. The pairing of these two wavelengths offers sensitivity to both temperature and humidity fluctuations, and the correlation between wavelengths is also used to retrieve the path-averaged temperature–humidity correlation. Comparison is made with structure parameters calculated from an eddy covariance station located close to the centre of the scintillometer path. The performance of the measurement techniques under different conditions is discussed. Similar behaviour is seen between the two data sets at sub-daily timescales. For the two summer-to-winter periods presented here, similar evolution is displayed across the seasons. A higher vegetation fraction within the scintillometer source area is consistent with the lower Bowen ratio observed (midday Bowen ratio < 1) compared with more built-up areas around the eddy covariance station. The energy partitioning is further explored in the companion paper.

Infrared and millimetre-wave scintillometry in the suburban environment – Part 2: Large-area sensible and latent heat fluxes

A millimetre-wave scintillometer was paired with an infrared scintillometer, enabling estimation of large-area evapotranspiration across northern Swindon, a suburban area in the UK. Both sensible and latent heat fluxes can be obtained using this "two-wavelength" technique, as it is able to provide both temperature and humidity structure parameters, offering a major advantage over conventional single-wavelength scintillometry. The first paper of this two-part series presented the measurement theory and structure parameters. In this second paper, heat fluxes are obtained and analysed. These fluxes, estimated using two-wavelength scintillometry over an urban area, are the first of their kind. Source area modelling suggests the scintillometric fluxes are representative of 5–10 km2. For comparison, local-scale (0.05–0.5 km2) fluxes were measured by an eddy covariance station. Similar responses to seasonal changes are evident at the different scales but the energy partitioning varies between source areas. The response to moisture availability is explored using data from 2 consecutive years with contrasting rainfall patterns (2011–2012). This extensive data set offers insight into urban surface-atmosphere interactions and demonstrates the potential for two-wavelength scintillometry to deliver fluxes over mixed land cover, typically representative of an area 1–2 orders of magnitude greater than for eddy covariance measurements. Fluxes at this scale are extremely valuable for hydro-meteorological model evaluation and assessment of satellite data products

JULES-crop: a parametrisation of crops in the Joint UK Land Environment Simulator

Studies of climate change impacts on the terrestrial biosphere have been completed without recognition of the integrated nature of the biosphere. Improved assessment of the impacts of climate change on food and water security requires the development and use of models not only representing each component but also their interactions. To meet this requirement the Joint UK Land Environment Simulator (JULES) land surface model has been modified to include a generic parametrisation of annual crops. The new model, JULES-crop, is described and evaluation at global and site levels for the four globally important crops; wheat, soybean, maize and rice. JULES-crop demonstrates skill in simulating the inter-annual variations of yield for maize and soybean at the global and country levels, and for wheat for major spring wheat producing countries. The impact of the new parametrisation, compared to the standard configuration, on the simulation of surface heat fluxes is largely an alteration of the partitioning between latent and sensible heat fluxes during the later part of the growing season. Further evaluation at the site level shows the model captures the seasonality of leaf area index, gross primary production and canopy height better than in the standard JULES. However, this does not lead to an improvement in the simulation of sensible and latent heat fluxes. The performance of JULES-crop from both an Earth system and crop yield model perspective is encouraging. However, more effort is needed to develop the parametrisation of the model for specific applications. Key future model developments identified include the introduction of processes such as irrigation and nitrogen limitation which will enable better representation of the spatial variability in yield.

Detection of coherent airstreams using cluster analysis: application to an extratropical cyclone

Flow in geophysical fluids is commonly summarized by coherent streams, for example conveyor belt flows in extratropical cyclones or jet streaks in the upper troposphere. Typically, parcel trajectories are calculated from the flow field and subjective thresholds are used to distinguish coherent streams of interest. This methodology contribution develops a more objective approach to distinguish coherent airstreams within extratropical cyclones. Agglomerative clustering is applied to trajectories along with a method to identify the optimal number of cluster classes. The methodology is applied to trajectories associated with the low-level jets of a well-studied extratropical cyclone. For computational efficiency, a constraint that trajectories must pass through these jet regions is applied prior to clustering; the partitioning into different airstreams is then performed by the agglomerative clustering. It is demonstrated that the methodology can identify the salient flow structures of cyclones: the warm and cold conveyor belts. A test focusing on the airstreams terminating at the tip of the bent-back front further demonstrates the success of the method in that it can distinguish fine-scale flow structure such as descending sting jet airstreams.

Glacial-interglacial temperature change in the tropical West Pacific: a comparison of stalagmite-based paleo-thermometers

In the tropics, geochemical records from stalagmites have so far mainly been used to qualitatively reconstruct changes in precipitation, but several new methods to reconstruct past temperatures from stalagmite material have emerged recently: i) liquid–vapor homogenization of fluid inclusion water ii) noble gas concentrations in fluid inclusion water, iii) the partitioning of oxygen isotopes between fluid inclusion water and calcite, and iv) the abundance of the 13C18O16O (‘clumped’) isotopologue in calcite. We present, for the first time, a direct comparison of these four paleo-thermometers by applying them to a fossil stalagmite covering nearly two glacial–interglacial cycles (Marine Isotope Stages (MIS) 12–9) and to two modern stalagmites, all from northern Borneo. The temperature estimates from the different methods agree in most cases within errors for both the old and recent samples; reconstructed formation temperatures of the recent samples match within 2-sigma errors with measured cave temperatures. However, slight but systematic deviations are observed between noble gas and liquid–vapor homogenization temperatures. Whereas the temperature sensitivity of fluid inclusion δ18O and clumped isotopes is currently debated, we find that the calibration of Tremaine et al. (2011) for fluid inclusion δ18O and a synthetic calcite-based clumped isotope calibration (Ziegler et al., in prep.) yield temperature estimates consistent with the other methods. All methods (with the potential exception of clumped isotopes) show excellent agreement on the amplitude of glacial–interglacial temperature change, indicating temperature shifts of 4–5 °C. This amplitude is similar to the amplitude of Mg/Ca-based regional sea surface temperature records, when correcting for sea level driven changes in cave elevation. Our reconstruction of tropical temperature evolution over the time period from 440 to 320 thousand years ago (ka) adds support to the view that climate sensitivity to varying greenhouse forcing is substantial also in the deep tropics.

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.

Persistence in diving American mink

Background American mink forage on land and in water, with aquatic prey often constituting a large proportion of their diet. Their long, thin body shape and relatively poor insulation make them vulnerable to heat loss, particularly in water, yet some individuals dive over 100 times a day. At the level of individual dives, previous research found no difference in dive depth or duration, or the total number of dives per day between seasons, but mink did appear to make more dives per active hour in winter than in summer. There was also no difference in the depth or duration of individual dives between the sexes, but there was some evidence that females made more dives per day than males. However, because individual mink dives tend to be extremely short in duration, persistence (quantified as the number of consecutive dives performed) may be a more appropriate metric with which to compare diving behaviour under different scenarios. Results Mink performed up to 28 consecutive dives, and dived continually for up to 36 min. Periods of more loosely aggregated diving (termed ‘aquatic activity sessions’) comprised up to 80 dives, carried out over up to 162.8 min. Contrary to our predictions, persistence was inversely proportional to body weight, with small animals more persistent than large ones, and (for females, but not for males) increased with decreasing temperature. For both sexes, persistence was greater during the day than during the night. Conclusions The observed body weight effect may point to inter-sexual niche partitioning, since in mink the smallest animals are females and the largest are males. The results may equally point to individual specialism’s, since persistence was also highly variable among individuals. Given the energetic costs involved, the extreme persistence of some animals observed in winter suggests that the costs of occasional prolonged activity in cold water are outweighed by the energetic gains. Analysing dive persistence can provide information on an animal’s physical capabilities for performing multiple dives and may reveal how such behaviour is affected by different conditions. Further development of monitoring and biologging methodology to allow quantification of hunting success, and thus the rewards obtained under alternative scenarios, would be insightful.

Anger, preoccupied attachment and domain disorganisation in borderline personality disorder

Emotional dysregulation and attachment insecurity have been reported in borderline personality disorder (BPD). Domain disorganization, evidenced in poor regulation of emotions and behaviors in relation to the demands of different social domains, may be a distinguishing feature of BPD. Understanding the interplay between these factors may be critical for identifying interacting processes in BPD and potential subtypes of BPD. Therefore, we examined the joint and interactive effects of anger, preoccupied attachment, and domain disorganization on BPD traits in a clinical sample of 128 psychiatric patients. The results suggest that these factors contribute to BPD both independently and in interaction, even when controlling for other personality disorder traits and Axis I symptoms. In regression analyses, the interaction between anger and domain disorganization predicted BPD traits. In recursive partitioning analyses, two possible paths to BPD were identified: high anger combined with high domain disorganization and low anger combined with preoccupied attachment. These results may suggest possible subtypes of BPD or possible mechanisms by which BPD traits are established and maintained.

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.

Carbon trading for phosphorus gain: the balance between rhizosphere carboxylates and arbuscular mycorrhizal symbiosis in plant phosphorus acquisition

Two key plant adaptations for phosphorus (P) acquisition are carboxylate exudation into the rhizosphere and mycorrhizal symbioses. These target different soil P resources, presumably with different plant carbon costs. We examined the effect of inoculation with arbuscular mycorrhizal fungi (AMF) on amount of rhizosphere carboxylates and plant P uptake for 10 species of low-P adapted Kennedia grown for 23 weeks in low-P sand. Inoculation decreased carboxylates in some species (up to 50%), decreased plant dry weight (21%) and increased plant P content (23%). There was a positive logarithmic relationship between plant P content and the amount of rhizosphere citric acid for inoculated and uninoculated plants. Causality was indicated by experiments using sand where little citric acid was lost from the soil solution over 2 h and citric acid at low concentrations desorbed P into the soil solution. Senesced leaf P concentration was often low and P-resorption efficiencies reached >90%. In conclusion, we propose that mycorrhizally mediated resource partitioning occurred because inoculation reduced rhizosphere carboxylates, but increased plant P uptake. Hence, presumably, the proportion of plant P acquired from strongly sorbed sources decreased with inoculation, while the proportion from labile inorganic P increased. Implications for plant fitness under field conditions now require investigation.

Variability and trends in England and Wales precipitation

The England and Wales precipitation (EWP) dataset is a homogeneous time series of daily accumulations from 1931 to 2014, composed from rain gauge observations spanning the region. The daily regional-average precipitation statistics are shown to be well described by a Weibull distribution, which is used to define extremes in terms of percentiles. Computed trends in annual and seasonal precipitation are sensitive to the period chosen, due to large variability on interannual and decadal timescales. Atmospheric circulation patterns associated with seasonal precipitation variability are identified. These patterns project onto known leading modes of variability, all of which involve displacements of the jet stream and storm-track over the eastern Atlantic. The intensity of daily precipitation for each calendar season is investigated by partitioning all observations into eight intensity categories contributing equally to the total precipitation in the dataset. Contrary to previous results based on shorter periods, no significant trends of the most intense categories are found between 1931 and 2014. The regional-average precipitation is found to share statistical properties common to the majority of individual stations across England and Wales used in previous studies. Statistics of the EWP data are examined for multi-day accumulations up to 10 days, which are more relevant for river flooding. Four recent years (2000, 2007, 2008 and 2012) have a greater number of extreme events in the 3-and 5-day accumulations than any previous year in the record. It is the duration of precipitation events in these years that is remarkable, rather than the magnitude of the daily accumulations.

Effects of aphid infestation on Cd and Zn concentration in wheat

The effect of infestation by the aphid Metopolophium dirhodum on the concentration and mass partitioning of Cd and Zn was studied in wheat plants. Results show that infestation did not affect the concentration of either metal in the roots or shoots of wheat, but elevated concentrations in the ears. This appeared to result from the concentration of metal in the smaller ear mass of infested plants. Infestation did not significantly affect the mass partitioning of either metal in any part of the plant, some 10% of both metals being allocated to the roots. However, the two metals contrasted in their partitioning in the aerial parts of the plant, with ca. 60% of Cd mass partitioned in the shoots and ca. 67% of Zn mass partitioned in the ears. The possible effects of infestation on the transfer of Cd and Zn from the soil to cereal aphids are discussed.