The depression of tropical snow lines at the Last Glacial Maximum: what can we learn from climate and model experiments?

Analyses of simulations of the last glacial maximum (LGM) made with 17 atmospheric general circulation models (AGCMs) participating in the Paleoclimate Modelling Intercomparison Project, and a high-resolution (T106) version of one of the models (CCSR1), show that changes in the elevation of tropical snowlines (as estimated by the depression of the maximum altitude of the 0 °C isotherm) are primarily controlled by changes in sea-surface temperatures (SSTs). The correlation between the two variables, averaged for the tropics as a whole, is 95%, and remains >80% even at a regional scale. The reduction of tropical SSTs at the LGM results in a drier atmosphere and hence steeper lapse rates. Changes in atmospheric circulation patterns, particularly the weakening of the Asian monsoon system and related atmospheric humidity changes, amplify the reduction in snowline elevation in the northern tropics. Colder conditions over the tropical oceans combined with a weakened Asian monsoon could produce snowline lowering of up to 1000 m in certain regions, comparable to the changes shown by observations. Nevertheless, such large changes are not typical of all regions of the tropics. Analysis of the higher resolution CCSR1 simulation shows that differences between the free atmospheric and along-slope lapse rate can be large, and may provide an additional factor to explain regional variations in observed snowline changes.

Incorporating toxicokinetics into an individual-based model for more realistic pesticide exposure estimates: A case study of the wood mouse

The potential risk of agricultural pesticides to mammals typically depends on internal concentrations within individuals, and these are determined by the amount ingested and by absorption, distribution, metabolism, and excretion (ADME). Pesticide residues ingested depend, amongst other things, on individual spatial choices which determine how much and when feeding sites and areas of pesticide application overlap, and can be calculated using individual-based models (IBMs). Internal concentrations can be calculated using toxicokinetic (TK) models, which are quantitative representations of ADME processes. Here we provide a population model for the wood mouse (Apodemus sylvaticus) in which TK submodels were incorporated into an IBM representation of individuals making choices about where to feed. This allows us to estimate the contribution of individual spatial choice and TK processes to risk. We compared the risk predicted by four IBMs: (i) “AllExposed-NonTK”: assuming no spatial choice so all mice have 100% exposure, no TK, (ii) “AllExposed-TK”: identical to (i) except that the TK processes are included where individuals vary because they have different temporal patterns of ingestion in the IBM, (iii) “Spatial-NonTK”: individual spatial choice, no TK, and (iv) “Spatial-TK”: individual spatial choice and with TK. The TK parameters for hypothetical pesticides used in this study were selected such that a conventional risk assessment would fail. Exposures were standardised using risk quotients (RQ; exposure divided by LD50 or LC50). We found that for the exposed sub-population including either spatial choice or TK reduced the RQ by 37–85%, and for the total population the reduction was 37–94%. However spatial choice and TK together had little further effect in reducing RQ. The reasons for this are that when the proportion of time spent in treated crop (PT) approaches 1, TK processes dominate and spatial choice has very little effect, and conversely if PT is small spatial choice dominates and TK makes little contribution to exposure reduction. The latter situation means that a short time spent in the pesticide-treated field mimics exposure from a small gavage dose, but TK only makes a substantial difference when the dose was consumed over a longer period. We concluded that a combined TK-IBM is most likely to bring added value to the risk assessment process when the temporal pattern of feeding, time spent in exposed area and TK parameters are at an intermediate level; for instance wood mice in foliar spray scenarios spending more time in crop fields because of better plant cover.

Evaluation of boundary-layer type in a weather forecast model utilising long-term Doppler lidar observations

Many studies evaluating model boundary-layer schemes focus either on near-surface parameters or on short-term observational campaigns. This reflects the observational datasets that are widely available for use in model evaluation. In this paper we show how surface and long-term Doppler lidar observations, combined in a way to match model representation of the boundary layer as closely as possible, can be used to evaluate the skill of boundary-layer forecasts. We use a 2-year observational dataset from a rural site in the UK to evaluate a climatology of boundary layer type forecast by the UK Met Office Unified Model. In addition, we demonstrate the use of a binary skill score (Symmetric Extremal Dependence Index) to investigate the dependence of forecast skill on season, horizontal resolution and forecast leadtime. A clear diurnal and seasonal cycle can be seen in the climatology of both the model and observations, with the main discrepancies being the model overpredicting cumulus capped and decoupled stratocumulus capped boundary-layers and underpredicting well mixed boundary-layers. Using the SEDI skill score the model is most skillful at predicting the surface stability. The skill of the model in predicting cumulus capped and stratocumulus capped stable boundary layer forecasts is low but greater than a 24 hr persistence forecast. In contrast, the prediction of decoupled boundary-layers and boundary-layers with multiple cloud layers is lower than persistence. This process based evaluation approach has the potential to be applied to other boundary-layer parameterisation schemes with similar decision structures.

A numerical model of the ionospheric signatures of time-varying magnetic reconnection: III. Quasi-instantaneous convection responses in the Cowley-Lockwood paradigm

Using a numerical implementation of the Cowley and Lockwood (1992) model of flow excitation in the magnetosphere–ionosphere (MI) system, we show that both an expanding (on a _12-min timescale) and a quasiinstantaneous response in ionospheric convection to the onset of magnetopause reconnection can be accommodated by the Cowley–Lockwood conceptual framework. This model has a key feature of time dependence, necessarily considering the history of the coupled MI system. We show that a residual flow, driven by prior magnetopause reconnection, can produce a quasi-instantaneous global ionospheric convection response; perturbations from an equilibrium state may also be present from tail reconnection, which will superpose constructively to give a similar effect. On the other hand, when the MI system is relatively free of pre-existing flow, we can most clearly see the expanding nature of the response. As the open-closed field line boundary will frequently be in motion from such prior reconnection (both at the dayside magnetopause and in the cross-tail current sheet), it is expected that there will usually be some level of combined response to dayside reconnection.

The greening of Arabia: multiple opportunities for human occupation of the Arabian peninsula during the Late Pleistocene inferred from an ensemble of climate model simulations

Climate models are potentially useful tools for addressing human dispersals and demographic change. The Arabian Peninsula is becoming increasingly significant in the story of human dispersals out of Africa during the Late Pleistocene. Although characterised largely by arid environments today, emerging climate records indicate that the peninsula was wetter many times in the past, suggesting that the region may have been inhabited considerably more than hitherto thought. Explaining the origins and spatial distribution of increased rainfall is challenging because palaeoenvironmental research in the region is in an early developmental stage. We address environmental oscillations by assembling and analysing an ensemble of five global climate models (CCSM3, COSMOS, HadCM3, KCM, and NorESM). We focus on precipitation, as the variable is key for the development of lakes, rivers and savannas. The climate models generated here were compared with published palaeoenvironmental data such as palaeolakes, speleothems and alluvial fan records as a means of validation. All five models showed, to varying degrees, that the Arabia Peninsula was significantly wetter than today during the Last Interglacial (130 ka and 126/125 ka timeslices), and that the main source of increased rainfall was from the North African summer monsoon rather than the Indian Ocean monsoon or from Mediterranean climate patterns. Where available, 104 ka (MIS 5c), 56 ka (early MIS 3) and 21 ka (LGM) timeslices showed rainfall was present but not as extensive as during the Last Interglacial. The results favour the hypothesis that humans potentially moved out of Africa and into Arabia on multiple occasions during pluvial phases of the Late Pleistocene.

Terrestrial biosphere changes over the last 120 kyr and their impact on ocean δ 13C

A new global synthesis and biomization of long (>40 kyr) pollen-data records is presented, and used with simulations from the HadCM3 and FAMOUS climate models to analyse the dynamics of the global terrestrial biosphere and carbon storage over the last glacial–interglacial cycle. Global modelled (BIOME4) biome distributions over time generally agree well with those inferred from pollen data. The two climate models show good agreement in global net primary productivity (NPP). NPP is strongly influenced by atmospheric carbon dioxide (CO2) concentrations through CO2 fertilization. The combined effects of modelled changes in vegetation and (via a simple model) soil carbon result in a global terrestrial carbon storage at the Last Glacial Maximum that is 210–470 Pg C less than in pre-industrial time. Without the contribution from exposed glacial continental shelves the reduction would be larger, 330–960 Pg C. Other intervals of low terrestrial carbon storage include stadial intervals at 108 and 85 ka BP, and between 60 and 65 ka BP during Marine Isotope Stage 4. Terrestrial carbon storage, determined by the balance of global NPP and decomposition, influences the stable carbon isotope composition (δ13C) of seawater because terrestrial organic carbon is depleted in 13C. Using a simple carbon-isotope mass balance equation we find agreement in trends between modelled ocean δ13C based on modelled land carbon storage, and palaeo-archives of ocean δ13C, confirming that terrestrial carbon storage variations may be important drivers of ocean δ13C changes.

Effects of orange juice formulation on prebiotic functionality using an in vitro colonic model sytem

A three-stage continuous fermentative colonic model system was used to monitor in vitro the effect of different orange juice formulations on prebiotic activity. Three different juices with and without Bimuno, a GOS mixture containing galactooligosaccharides (B-GOS) were assessed in terms of their ability to induce a bifidogenic microbiota. The recipe development was based on incorporating 2.75g B-GOS into a 250 ml serving of juice (65°Brix of concentrate juice). Alongside the production of B-GOS juice, a control juice - orange juice without any additional Bimuno and a positive control juice, containing all the components of Bimuno (glucose, galactose and lactose) in the same relative proportions with the exception of B-GOS were developed. Ion Exchange Chromotography analysis was used to test the maintenance of bimuno components after the production process. Data showed that sterilisation had no significant effect on concentration of B-GOS and simple sugars. The three juice formulations were digested under conditions resembling the gastric and small intestinal environments. Main bacterial groups of the faecal microbiota were evaluated throughout the colonic model study using 16S rRNA-based fluorescence in situ hybridization (FISH). Potential effects of supplementation of the juices on microbial metabolism were studied measuring short chain fatty acids (SCFAs) using gas chromatography. Furthermore, B-GOS juices showed positive modulations of the microbiota composition and metabolic activity. In particular, numbers of faecal bifidobacteria and lactobacilli were significantly higher when B-GOS juice was fermented compared to controls. Furthermore, fermentation of B-GOS juice resulted in an increase in Roseburia subcluster and concomitantly increased butyrate production, which is of potential benefit to the host. In conclusion, this study has shown B-GOS within orange juice can have a beneficial effect on the fecal microbiota.

Synergistic inhibition of Haemonchus contortus exsheathment by flavonoid monomers and condensed tannins

This study investigated the separate and combined anthelmintic (AH) effects of different phenolic compounds, including condensed tannins and flavonoids, all of which are known to occur in willow leaves, a potentially valuable dry season feed. A range of contrasting model tannins, which span the whole range of willow tannins, were isolated from tilia flowers, goat willow leaves, black currant leaves and red currant leaves. All together, the tested compounds represented the major tannin types (procyanidins and prodelphinidins) and flavonoid types (flavonols, flavones and flavanones). The larval exsheathment inhibition assay (LEIA) was used to assess their in vitro effects on Haemonchus contortus third stage larvae. Arbutin, vanillic acid, and taxifolin proved to be ineffective whereas naringenin, quercetin and luteolin were highly effective at 250 μM concentrations. Procyanidin (PC) tannins tended to be less active than prodelphinidin tannins (PD). Experiments with combinations of tannins and quercetin or luteolin revealed for the first time the existence of synergistic AH effects between tannins and flavonoid monomers. They also provided evidence that synergistic effects appear to occur at slightly lower concentrations of PC than PD. This suggests that the AH activity of condensed tannins can be significantly enhanced by the addition of quercetin or luteolin. This information may prove useful for plant breeding or selection and for designing optimal feed mixtures.

Tracer Leakage from Modeled Agulhas Rings

In a numerical, isopycnal, ocean model the mixing is investigated with the environment of two idealized Agulhas rings, one that splits, and one that remains coherent. The evolution of a passive tracer , initially contained within the rings, shows that tracer leakage is associated with the for mation of filaments in the early stage of ring evolution. These filaments reach down to the ther mocline. In the deepest layers leakage occurs on a larger scale. Self-advection of the rings is ver y irregular , and it is not possible to compute a Lagrangian boundar y i n order to estimate the transport of leakage from the rings. T o describe the processes that gover n tracer leakage, in a coordinate frame moving with the ring a kinematic separatrix is defined in the streamfunction field for the nondivergent flow . Initially , filaments arise because of the elongation of the ring, which is mainly gover ned by an m 5 2 instability that is collaborating with differential rotation. Because of beta, the symmetr y i s destroyed related to the separatrix associated with a stagnation point in the flow . The filament upstream of the stagnation point grows much faster and is associated with the bulk of tracer leakage. Mixing is enhanced by time dependence of the separatrix. As a result, there are no large differences between the leakage from a coherent ring, where the m 5 2 instability equilibrates, and from a splitting ring, where the m 5 2 instability keeps growing, which confir ms that the amount of leakage is mainly gover ned by the ring’ s initial defor mation combined with unsteady self-advection of the ring and not by the splitting of the ring. The decay of tracer content in the ther mocline shows that in the first months up to 40% of the ring water can be mixed with the environment. In deeper layers the decay of tracer content may reach up to 90%.

The combined influences of westerly phase of the quasi biennial oscillation and 11-year solar maximum conditions on the Northern Hemisphere extratropical winter circulation

The combined influences of the westerly phase of the quasi-biennial oscillation (QBO-W) and solar maximum (Smax) conditions on the Northern Hemisphere extratropical winter circulation are investigated using reanalysis data and Center for Climate System Research/National Institute for Environmental Studies chemistry climate model (CCM) simulations. The composite analysis for the reanalysis data indicates strengthened polar vortex in December followed by weakened polar vortex in February–March for QBO-W during Smax (QBO-W/Smax) conditions. This relationship need not be specific to QBO-W/Smax conditions but may just require strengthened vortex in December, which is more likely under QBO-W/Smax. Both the reanalysis data and CCM simulations suggest that dynamical processes of planetary wave propagation and meridional circulation related to QBO-W around polar vortex in December are similar in character to those related to Smax; furthermore, both processes may work in concert to maintain stronger vortex during QBO-W/Smax. In the reanalysis data, the strengthened polar vortex in December is associated with the development of north–south dipole tropospheric anomaly in the Atlantic sector similar to the North Atlantic oscillation (NAO) during December–January. The structure of the north–south dipole anomaly has zonal wavenumber 1 (WN1) component, where the longitude of anomalous ridge overlaps with that of climatological ridge in the North Atlantic in January. This implies amplification of the WN1 wave and results in the enhancement of the upward WN1 propagation from troposphere into stratosphere in January, leading to the weakened polar vortex in February–March. Although WN2 waves do not play a direct role in forcing the stratospheric vortex evolution, their tropospheric response to QBO-W/Smax conditions appears to be related to the maintenance of the NAO-like anomaly in the high-latitude troposphere in January. These results may provide a possible explanation for the mechanisms underlying the seasonal evolution of wintertime polar vortex anomalies during QBO-W/Smax conditions and the role of troposphere in this evolution.

Simulating the impact on health of internalising the cost of carbon in food prices combined with a tax on sugar-sweetened beverages

Rising greenhouse gas emissions (GHGEs) have implications for health and up to 30 % of emissions globally are thought to arise from agriculture. Synergies exist between diets low in GHGEs and health however some foods have the opposite relationship, such as sugar production being a relatively low source of GHGEs. In order to address this and to further characterise a healthy sustainable diet, we model the effect on UK non-communicable disease mortality and GHGEs of internalising the social cost of carbon into the price of food alongside a 20 % tax on sugar sweetened beverages (SSBs). Developing previously published work, we simulate four tax scenarios: (A) a GHGEs tax of £2.86/tonne of CO2 equivalents (tCO2e)/100 g product on all products with emissions greater than the mean across all food groups (0.36 kgCO2e/100 g); (B) scenario A but with subsidies on foods with emissions lower than 0.36 kgCO2e/100 g such that the effect is revenue neutral; (C) scenario A but with a 20 % sales tax on SSBs; (D) scenario B but with a 20 % sales tax on SSBs. An almost ideal demand system is used to estimate price elasticities and a comparative risk assessment model is used to estimate changes to non-communicable disease mortality. We estimate that scenario A would lead to 300 deaths delayed or averted, 18,900 ktCO2e fewer GHGEs, and £3.0 billion tax revenue; scenario B, 90 deaths delayed or averted and 17,100 ktCO2e fewer GHGEs; scenario C, 1,200 deaths delayed or averted, 18,500 ktCO2e fewer GHGEs, and £3.4 billion revenue; and scenario D, 2,000 deaths delayed or averted and 16,500 ktCO2e fewer GHGEs. Deaths averted are mainly due to increased fibre and reduced fat consumption; a SSB tax reduces SSB and sugar consumption. Incorporating the social cost of carbon into the price of food has the potential to improve health, reduce GHGEs, and raise revenue. The simple addition of a tax on SSBs can mitigate negative health consequences arising from sugar being low in GHGEs. Further conflicts remain, including increased consumption of unhealthy foods such as cakes and nutrients such as salt.

Evaluating and optimizing oral formulations of live bacterial vaccines using a gastro-small intestine model

Gastrointestinal (GI) models that mimic physiological conditions in vitro are important tools for developing and optimizing biopharmaceutical formulations. Oral administration of live attenuated bacterial vaccines (LBV) can safely and effectively promote mucosal immunity but new formulations are required that provide controlled release of optimal numbers of viable bacterial cells, which must survive gastrointestinal transit overcoming various antimicrobial barriers. Here, we use a gastro-small intestine gut model of human GI conditions to study the survival and release kinetics of two oral LBV formulations: the licensed typhoid fever vaccine Vivotif comprising enteric coated capsules; and an experimental formulation of the model vaccine Salmonella Typhimurium SL3261 dried directly onto cast enteric polymer films and laminated to form a polymer film laminate (PFL). Neither formulation released significant numbers of viable cells when tested in the complete gastro-small intestine model. The poor performance in delivering viable cells could be attributed to a combination of acid and bile toxicity plus incomplete release of cells for Vivotif capsules, and to bile toxicity alone for PFL. To achieve effective protection from intestinal bile in addition to effective acid resistance, bile adsorbent resins were incorporated into the PFL to produce a new formulation, termed BR-PFL. Efficient and complete release of 4.4x107 live cells per dose was achieved from BR-PFL at distal intestinal pH, with release kinetics controlled by the composition of the enteric polymer film, and no loss in viability observed in any stage of the GI model. Use of this in vitro GI model thereby allowed rational design of an oral LBV formulation to maximize viable cell release.

Uncertainties in tidally adjusted estimates of sea level rise flooding (bathtub model) for the Greater London

Sea-level rise (SLR) from global warming may have severe consequences for coastal cities, particularly when combined with predicted increases in the strength of tidal surges. Predicting the regional impact of SLR flooding is strongly dependent on the modelling approach and accuracy of topographic data. Here, the areas under risk of sea water flooding for London boroughs were quantified based on the projected SLR scenarios reported in Intergovernmental Panel on Climate Change (IPCC) fifth assessment report (AR5) and UK climatic projections 2009 (UKCP09) using a tidally-adjusted bathtub modelling approach. Medium- to very high-resolution digital elevation models (DEMs) are used to evaluate inundation extents as well as uncertainties. Depending on the SLR scenario and DEMs used, it is estimated that 3%–8% of the area of Greater London could be inundated by 2100. The boroughs with the largest areas at risk of flooding are Newham, Southwark, and Greenwich. The differences in inundation areas estimated from a digital terrain model and a digital surface model are much greater than the root mean square error differences observed between the two data types, which may be attributed to processing levels. Flood models from SRTM data underestimate the inundation extent, so their results may not be reliable for constructing flood risk maps. This analysis provides a broad-scale estimate of the potential consequences of SLR and uncertainties in the DEM-based bathtub type flood inundation modelling for London boroughs.