Model inter-comparison between statistical and dynamic model assessments of the long-term stability of blanket peat in Great Britain (1940-2099)

We compared output from 3 dynamic process-based models (DMs: ECOSSE, MILLENNIA and the Durham Carbon Model) and 9 bioclimatic envelope models (BCEMs; including BBOG ensemble and PEATSTASH) ranging from simple threshold to semi-process-based models. Model simulations were run at 4 British peatland sites using historical climate data and climate projections under a medium (A1B) emissions scenario from the 11-RCM (regional climate model) ensemble underpinning UKCP09. The models showed that blanket peatlands are vulnerable to projected climate change; however, predictions varied between models as well as between sites. All BCEMs predicted a shift from presence to absence of a climate associated with blanket peat, where the sites with the lowest total annual precipitation were closest to the presence/absence threshold. DMs showed a more variable response. ECOSSE predicted a decline in net C sink and shift to net C source by the end of this century. The Durham Carbon Model predicted a smaller decline in the net C sink strength, but no shift to net C source. MILLENNIA predicted a slight overall increase in the net C sink. In contrast to the BCEM projections, the DMs predicted that the sites with coolest temperatures and greatest total annual precipitation showed the largest change in carbon sinks. In this model inter-comparison, the greatest variation in model output in response to climate change projections was not between the BCEMs and DMs but between the DMs themselves, because of different approaches to modelling soil organic matter pools and decomposition amongst other processes. The difference in the sign of the response has major implications for future climate feedbacks, climate policy and peatland management. Enhanced data collection, in particular monitoring peatland response to current change, would significantly improve model development and projections of future change.

Social networks and positive and negative affect

We followed 100 university students in the UK for one week, instructing them to record all face-to-face, phone and digital contacts during the day as well as their positive and negative affect. We wanted to see how positive and negative affect spread around a social network while taking into account participants’ socio-demographic data, personality, general health and gratitude scores. We focused on the participants’ connections with those in their class; excluding friends and family outside this group. The data was analysed using actor-based models implemented in SIENA. Results show differences between positive and negative affect dynamics in this environment and an influence of personality traits on the average number and rate of communication.

The effects of landscape modifications on the long-term persistence of animal populations

Background: The effects of landscape modifications on the long-term persistence of wild animal populations is of crucial importance to wildlife managers and conservation biologists, but obtaining experimental evidence using real landscapes is usually impossible. To circumvent this problem we used individual-based models (IBMs) of interacting animals in experimental modifications of a real Danish landscape. The models incorporate as much as possible of the behaviour and ecology of four species with contrasting life-history characteristics: skylark (Alauda arvensis), vole (Microtus agrestis), a ground beetle (Bembidion lampros) and a linyphiid spider (Erigone atra). This allows us to quantify the population implications of experimental modifications of landscape configuration and composition. Methodology/Principal Findings: Starting with a real agricultural landscape, we progressively reduced landscape complexity by (i) homogenizing habitat patch shapes, (ii) randomizing the locations of the patches, and (iii) randomizing the size of the patches. The first two steps increased landscape fragmentation. We assessed the effects of these manipulations on the long-term persistence of animal populations by measuring equilibrium population sizes and time to recovery after disturbance. Patch rearrangement and the presence of corridors had a large effect on the population dynamics of species whose local success depends on the surrounding terrain. Landscape modifications that reduced population sizes increased recovery times in the short-dispersing species, making small populations vulnerable to increasing disturbance. The species that were most strongly affected by large disturbances fluctuated little in population sizes in years when no perturbations took place. Significance: Traditional approaches to the management and conservation of populations use either classical methods of population analysis, which fail to adequately account for the spatial configurations of landscapes, or landscape ecology, which accounts for landscape structure but has difficulty predicting the dynamics of populations living in them. Here we show how realistic and replicable individual-based models can bridge the gap between non-spatial population theory and non-dynamic landscape ecology. A major strength of the approach is its ability to identify population vulnerabilities not detected by standard population viability analyses.

Property company performance and real interest rates: a regime-switching approach

Linear models of property market performance may be misspecified if there exist distinct states where the market drivers behave in different ways. This paper examines the applicability of non-linear regime-based models. A Self Exciting Threshold Autoregressive (SETAR) model is applied to property company share data, using the real rate of interest to define regimes. Distinct regimes appear exhibiting markedly different market behaviour. The model both casts doubt on the specification of conventional linear models and offers the possibility of developing effective trading rules for real estate equities.

What people study when they study Twitter: classifying Twitter related academic papers

PURPOSE: Since its introduction in 2006, messages posted to the microblogging system Twitter have provided a rich dataset for researchers, leading to the publication of over a thousand academic papers. This paper aims to identify this published work and to classify it in order to understand Twitter based research. DESIGN/METHODOLOGY/APPROACH: Firstly the papers on Twitter were identified. Secondly, following a review of the literature, a classification of the dimensions of microblogging research was established. Thirdly, papers were qualitatively classified using open coded content analysis, based on the paper’s title and abstract, in order to analyze method, subject, and approach. FINDINGS: The majority of published work relating to Twitter concentrates on aspects of the messages sent and details of the users. A variety of methodological approaches are used across a range of identified domains. RESEARCH LIMITATIONS/IMPLICATIONS: This work reviewed the abstracts of all papers available via database search on the term “Twitter” and this has two major implications: 1) the full papers are not considered and so works may be misclassified if their abstract is not clear, 2) publications not indexed by the databases, such as book chapters, are not included. ORIGINALITY/VALUE: To date there has not been an overarching study to look at the methods and purpose of those using Twitter as a research subject. Our major contribution is to scope out papers published on Twitter until the close of 2011. The classification derived here will provide a framework within which researchers studying Twitter related topics will be able to position and ground their work

Deriving global flood hazard maps of fluvial floods through a physical model cascade

Global flood hazard maps can be used in the assessment of flood risk in a number of different applications, including (re)insurance and large scale flood preparedness. Such global hazard maps can be generated using large scale physically based models of rainfall-runoff and river routing, when used in conjunction with a number of post-processing methods. In this study, the European Centre for Medium Range Weather Forecasts (ECMWF) land surface model is coupled to ERA-Interim reanalysis meteorological forcing data, and resultant runoff is passed to a river routing algorithm which simulates floodplains and flood flow across the global land area. The global hazard map is based on a 30 yr (1979–2010) simulation period. A Gumbel distribution is fitted to the annual maxima flows to derive a number of flood return periods. The return periods are calculated initially for a 25×25 km grid, which is then reprojected onto a 1×1 km grid to derive maps of higher resolution and estimate flooded fractional area for the individual 25×25 km cells. Several global and regional maps of flood return periods ranging from 2 to 500 yr are presented. The results compare reasonably to a benchmark data set of global flood hazard. The developed methodology can be applied to other datasets on a global or regional scale.

Mathematics applied to the climate system: outstanding challenges and recent progress

The societal need for reliable climate predictions and a proper assessment of their uncertainties is pressing. Uncertainties arise not only from initial conditions and forcing scenarios, but also from model formulation. Here, we identify and document three broad classes of problems, each representing what we regard to be an outstanding challenge in the area of mathematics applied to the climate system. First, there is the problem of the development and evaluation of simple physically based models of the global climate. Second, there is the problem of the development and evaluation of the components of complex models such as general circulation models. Third, there is the problem of the development and evaluation of appropriate statistical frameworks. We discuss these problems in turn, emphasizing the recent progress made by the papers presented in this Theme Issue. Many pressing challenges in climate science require closer collaboration between climate scientists, mathematicians and statisticians. We hope the papers contained in this Theme Issue will act as inspiration for such collaborations and for setting future research directions.

Pediastrum species as potential indicators of lake-level change in tropical South America

We explored the potential for using Pediastrum (Meyen), a genus of green alga commonly found in palaeoecological studies, as a proxy for lake-level change in tropical South America. The study site, Laguna La Gaiba (LLG) (17°45′S, 57°40′W), is a broad, shallow lake located along the course of the Paraguay River in the Pantanal, a 135,000-km2 tropical wetland located mostly in western Brazil, but extending into eastern Bolivia. Fourteen surface sediment samples were taken from LLG across a range of lake depths (2-5.2 m) and analyzed for Pediastrum. We found seven species, of which P. musteri (Tell et Mataloni), P. argentiniense (Bourr. et Tell), and P. cf. angulosum (Ehrenb.) ex Menegh. were identified as potential indicators of lake level. Results of the modern dataset were applied to 31 fossil Pediastrum assemblages spanning the early Holocene (12.0 kyr BP) to present to infer past lake level changes qualitatively. Early Holocene (12.0-9.8 kyr BP) assemblages do not show a clear signal, though abundance of P. simplex (Meyen) suggests relatively high lake levels. Absence of P. musteri, characteristic of deep, open water, and abundance of macrophyte-associated taxa indicate lake levels were lowest from 9.8 to 3.0 kyr BP. A shift to wetter conditions began at 4.4 kyr BP, indicated by the appearance of P. musteri, though inferred lake levels did not reach modern values until 1.4 kyr BP. The Pediastrum-inferred mid-Holocene lowstand is consistent with lower precipitation, previously inferred using pollen from this site, and is also in agreement with evidence for widespread drought in the South American tropics during the middle Holocene. An inference for steadily increasing lake level from 4.4 kyr BP to present is consistent with diatom-inferred water level rise at Lake Titicaca, and demonstrates coherence with the broad pattern of increasing monsoon strength from the late Holocene until present in tropical South America.

The AVOID programme’s new simulations of the global benefits of stringent climate change mitigation

Quantitative simulations of the global-scale benefits of climate change mitigation are presented, using a harmonised, self-consistent approach based on a single set of climate change scenarios. The approach draws on a synthesis of output from both physically-based and economics-based models, and incorporates uncertainty analyses. Previous studies have projected global and regional climate change and its impacts over the 21st century but have generally focused on analysis of business-as-usual scenarios, with no explicit mitigation policy included. This study finds that both the economics-based and physically-based models indicate that early, stringent mitigation would avoid a large proportion of the impacts of climate change projected for the 2080s. However, it also shows that not all the impacts can now be avoided, so that adaptation would also therefore be needed to avoid some of the potential damage. Delay in mitigation substantially reduces the percentage of impacts that can be avoided, providing strong new quantitative evidence for the need for stringent and prompt global mitigation action on greenhouse gas emissions, combined with effective adaptation, if large, widespread climate change impacts are to be avoided. Energy technology models suggest that such stringent and prompt mitigation action is technologically feasible, although the estimated costs vary depending on the specific modelling approach and assumptions.

Modelling human visual navigation using multi-view scene reconstruction

It is often assumed that humans generate a 3D reconstruction of the environment, either in egocentric or world-based coordinates, but the steps involved are unknown. Here, we propose two reconstruction-based models, evaluated using data from two tasks in immersive virtual reality. We model the observer’s prediction of landmark location based on standard photogrammetric methods and then combine location predictions to compute likelihood maps of navigation behaviour. In one model, each scene point is treated independently in the reconstruction; in the other, the pertinent variable is the spatial relationship between pairs of points. Participants viewed a simple environment from one location, were transported (virtually) to another part of the scene and were asked to navigate back. Error distributions varied substantially with changes in scene layout; we compared these directly with the likelihood maps to quantify the success of the models. We also measured error distributions when participants manipulated the location of a landmark to match the preceding interval, providing a direct test of the landmark-location stage of the navigation models. Models such as this, which start with scenes and end with a probabilistic prediction of behaviour, are likely to be increasingly useful for understanding 3D vision.

The time-course of morphological constraints: a study of plurals inside derived words

The avoidance of regular but not irregular plurals inside compounds (e.g. *rats eater vs. mice eater) has been one of the most widely studied morphological phenomena in the psycholinguistics literature. To examine whether the constraints that are responsible for this contrast have any general significance beyond compounding, we investigated derived word forms containing regular and irregular plurals in two experiments. Experiment 1 was an offline acceptability judgment task, and experiment 2 measured eye movements during reading derived words containing regular and irregular plurals and uninflected base nouns. The results from both experiments show that the constraint against regular plurals inside compounds generalizes to derived words. We argue that this constraint cannot be reduced to phonological properties, but is instead morphological in nature. The eye-movement data provide detailed information on the time-course of processing derived word forms indicating that early stages of processing are affected by a general constraint that disallows inflected words from feeding derivational processes, and that the more specific constraint against regular plurals comes in at a subsequent later stage of processing. We argue that these results are consistent with stage-based models of language processing.

Ecosystem effects of CO2 concentration: evidence from past climates

Atmospheric CO2 concentration has varied from minima of 170-200 ppm in glacials to maxima of 280-300 ppm in the recent interglacials. Photosynthesis by C-3 plants is highly sensitive to CO2 concentration variations in this range. Physiological consequences of the CO2 changes should therefore be discernible in palaeodata. Several lines of evidence support this expectation. Reduced terrestrial carbon storage during glacials, indicated by the shift in stable isotope composition of dissolved inorganic carbon in the ocean, cannot be explained by climate or sea-level changes. It is however consistent with predictions of current process-based models that propagate known physiological CO2 effects into net primary production at the ecosystem scale. Restricted forest cover during glacial periods, indicated by pollen assemblages dominated by non-arboreal taxa, cannot be reproduced accurately by palaeoclimate models unless CO2 effects on C-3-C-4 plant competition are also modelled. It follows that methods to reconstruct climate from palaeodata should account for CO2 concentration changes. When they do so, they yield results more consistent with palaeoclimate models. In conclusion, the palaeorecord of the Late Quaternary, interpreted with the help of climate and ecosystem models, provides evidence that CO2 effects at the ecosystem scale are neither trivial nor transient.

Effects of noise and by-catch on a Danish harbour porpoise population

Ships and wind turbines generate noise, which can have a negative impact on marine mammal populations by scaring animals away. Effective modelling of how this affects the populations has to take account of the location and timing of disturbances. Here we construct an individual-based model of harbour porpoises in the Inner Danish Waters. Individuals have their own energy budgets constructed using established principles of physiological ecology. Data are lacking on the spatial distribution of food which is instead inferred from knowledge of time-varying porpoise distributions. The model produces plausible patterns of population dynamics and matches well the age distribution of porpoises caught in by-catch. It estimates the effect of existing wind farms as a 10% reduction in population size when food recovers fast (after two days). Proposed new wind farms and ships do not result in further population declines. The population is however sensitive to variations in mortality resulting from by-catch and to the speed at which food recovers after being depleted. If food recovers slowly the effect of wind turbines becomes negligible, whereas ships are estimated to have a significant negative impact on the population. Annual by-catch rates ≥10% lead to monotonously decreasing populations and to extinction, and even the estimated by-catch rate from the adjacent area (approximately 4.1%) has a strong impact on the population. This suggests that conservation efforts should be more focused on reducing by-catch in commercial gillnet fisheries than on limiting the amount of anthropogenic noise. Individual-based models are unique in their ability to take account of the location and timing of disturbances and to show their likely effects on populations. The models also identify deficiencies in the existing database and can be used to set priorities for future field research.

On the cause of a magnetospheric flux transfer event

We present a detailed investigation of a magnetospheric flux transfer event (FTE) seen by the Active Magnetospheric Tracer Explorer (AMPTE) UKS and IRM satellites around 1046 UT on October 28, 1984. This event has been discussed many times previously in the literature and has been cited as support for a variety of theories of FTE formation. We make use of a model developed to reproduce ion precipitations seen in the cusp ionosphere. The analysis confirms that the FTE is well explained as a brief excursion into an open low-latitude boundary layer (LLBL), as predicted by two theories of magnetospheric FTEs: namely, that they are bulges in the open LLBL due to reconnection rate enhancements or that they are indentations of the magnetopause by magnetosheath pressure increases (but in the presence of ongoing steady reconnection). The indentation of the inner edge of the open LLBL that these two models seek to explain is found to be shallow for this event. The ion model reproduces the continuous evolution of the ion distribution function between the sheath-like population at the event center and the surrounding magnetospheric populations; it also provides an explanation of the high-pressure core of the event as comprising field lines that were reconnected considerably earlier than those that are draped over it to give the event boundary layer. The magnetopause transition parameter is used to isolate a field rotation on the boundaries of the core, which is subjected to the tangential stress balance test. The test identifies this to be a convecting structure, which is neither a rotational discontinuity (RD) nor a contact discontinuity, but could possibly be a slow shock. In addition, evidence for ion reflection off a weak RD on the magnetospheric side of this structure is found. The event structure is consistent in many ways with features predicted for the open LLBL by analytic MHD theories and by MHD and hybrid simulations. The de Hoffman-Teller velocity of the structure is significantly different from that of the magnetosheath flow, indicating that it is not an indentation caused by a high-pressure pulse in the sheath but is consistent with the motion of newly opened field lines (different from the sheath flow because of the magnetic tension force) deduced from the best fit to the ion data. However, we cannot here rule out the possibility that the sheath flow pattern has changed in the long interval between the two satellites observing the FTE and subsequently emerging into the magnetosheath; thus this test is not conclusive in this particular case. Analysis of the fitted elapsed time since reconnection shows that the core of the event was reconnected in one pulse and the event boundary layer was reconnected in a subsequent pulse. Between these two pulses is a period of very low (but nonzero) reconnection rate, which lasts about 14 mins. Thus the analysis supports, but does not definitively verify, the concept that the FTE is a partial passage into an open LLBL caused by a traveling bulge in that layer produced by a pulse in reconnection rate.

Speleothems as sensitive recorders of volcanic eruptions - the Bronze Age Minoan eruption recorded in a stalagmite from Turkey

Mounting evidence exists that variations in sulphur content in stalagmites are closely linked to changes in volcanic or anthropogenic atmospheric sulphur. The strong dependency of sulphur on soil pH and ecosystem storage, however, can result in a delay of several years to decades in the registration of volcanic eruptions and anthropogenic emissions by stalagmites. Here we present synchrotron-radiation based trace element analysis performed on a precisely-dated section of a stalagmite from Sofular Cave in Northern Turkey. As this section covers the time interval of the intensively studied Minoan volcanic eruption between 1600 and 1650 BC, we can test whether this vigorous eruption can be traced in a stalagmite. Of all measured trace elements, only bromine shows a clear short-lived peak at 1621±251621±25 BC, whereas sulphur and molybdenum show peaks later at 1617±251617±25 and 1589±251589±25 respectively. We suggest that all trace element peaks are related to the Minoan eruption, whereas the observed phasing of bromine, molybdenum and sulphur is related to differences in their retention rates in the soil above Sofular Cave. For the first time, we can show that bromine appears to be an ideal volcanic tracer in stalagmites, as it is a prominent volatile component in volcanic eruptions, can be easily leached in soils and rapidly transferred from the atmosphere through the soil and bedrock into the cave and stalagmite respectively. Highly resolved oxygen and carbon isotope profiles indicate that the Minoan eruption had no detectable climatic and environmental impact in Northern Turkey.