Impact of oceanic processes on the carbon cycle during the last termination

During the last termination (from ~18 000 years ago to ~9000 years ago), the climate significantly warmed and the ice sheets melted. Simultaneously, atmospheric CO2 increased from ~190 ppm to ~260 ppm. Although this CO2 rise plays an important role in the deglacial warming, the reasons for its evolution are difficult to explain. Only box models have been used to run transient simulations of this carbon cycle transition, but by forcing the model with data constrained scenarios of the evolution of temperature, sea level, sea ice, NADW formation, Southern Ocean vertical mixing and biological carbon pump. More complex models (including GCMs) have investigated some of these mechanisms but they have only been used to try and explain LGM versus present day steady-state climates. In this study we use a coupled climate-carbon model of intermediate complexity to explore the role of three oceanic processes in transient simulations: the sinking of brines, stratification-dependent diffusion and iron fertilization. Carbonate compensation is accounted for in these simulations. We show that neither iron fertilization nor the sinking of brines alone can account for the evolution of CO2, and that only the combination of the sinking of brines and interactive diffusion can simultaneously simulate the increase in deep Southern Ocean δ13C. The scenario that agrees best with the data takes into account all mechanisms and favours a rapid cessation of the sinking of brines around 18 000 years ago, when the Antarctic ice sheet extent was at its maximum. In this scenario, we make the hypothesis that sea ice formation was then shifted to the open ocean where the salty water is quickly mixed with fresher water, which prevents deep sinking of salty water and therefore breaks down the deep stratification and releases carbon from the abyss. Based on this scenario, it is possible to simulate both the amplitude and timing of the long-term CO2 increase during the last termination in agreement with ice core data. The atmospheric δ13C appears to be highly sensitive to changes in the terrestrial biosphere, underlining the need to better constrain the vegetation evolution during the termination.

The time-course of morphological processing in a second language

We report findings from psycholinguistic experiments investigating the detailed timing of processing morphologically complex words by proficient adult second (L2) language learners of English in comparison to adult native (L1) speakers of English. The first study employed the masked priming technique to investigate -ed forms with a group of advanced Arabic-speaking learners of English. The results replicate previously found L1/L2 differences in morphological priming, even though in the present experiment an extra temporal delay was offered after the presentation of the prime words. The second study examined the timing of constraints against inflected forms inside derived words in English using the eye-movement monitoring technique and an additional acceptability judgment task with highly advanced Dutch L2 learners of English in comparison to adult L1 English controls. Whilst offline the L2 learners performed native-like, the eye-movement data showed that their online processing was not affected by the morphological constraint against regular plurals inside derived words in the same way as in native speakers. Taken together, these findings indicate that L2 learners are not just slower than native speakers in processing morphologically complex words, but that the L2 comprehension system employs real-time grammatical analysis (in this case, morphological information) less than the L1 system.

Late Quaternary highstands at Lake Chilwa, Malawi: frequency, timing and possible forcing mechanisms in the last 44ka

The extensive shoreline deposits of Lake Chilwa, southern Malawi, a shallow water body today covering 600 km2 of a basin of 7500 km2, are investigated for their record of late Quaternary highstands. OSL dating, applied to 36 samples from five sediment cores from the northern and western marginal sand ridges, reveal a highstand record spanning 44 ka. Using two different grouping methods, highstand phases are identified at 43.7–33.3 ka, 26.2–21.0 ka and 17.9–12.0 ka (total error method) or 38.4–35.5 ka, 24.3–22.3 ka, 16.2–15.1 ka and 13.5–12.7 ka (Finite Mixture Model age components) with two further discrete events recorded at 11.01 ± 0.76 ka and 8.52 ± 0.56 ka. Highstands are comparable to the timing of wet phases from other basins in East and southern Africa, demonstrating wet conditions in the region before the LGM, which was dry, and a wet Lateglacial, which commenced earlier in the southern compared to northern hemisphere in East Africa. We find no evidence that wet phases are insolation driven, but analysis of the dataset and GCM modelling experiments suggest that Heinrich events may be associated with enhanced monsoon activity in East Africa in both timing and as a possible causal mechanism.

A review of time use models of residential electricity demand

Residential electricity demand in most European countries accounts for a major proportion of overall electricity consumption. The timing of residential electricity demand has significant impacts on carbon emissions and system costs. This paper reviews the data and methods used in time use studies in the context of residential electricity demand modelling. It highlights key issues which are likely to become more topical for research on the timing of electricity demand following the roll-out of smart metres.

Agriculture modifies the seasonal decline of breeding success in a tropical wild bird population

1.Habitat conversion for agriculture is a major driver of biodiversity loss, but our understanding of the demographic processes involved remains poor. We typically investigate the impacts of agriculture in isolation even though populations are likely to experience multiple, concurrent changes in the environment (e.g. land and climate change). Drivers of environmental change may interact to affect demography but the mechanisms have yet to be explored fully in wild populations. 2.Here, we investigate the mechanisms linking agricultural land-use with breeding success using long-term data for the formerly Critically Endangered Mauritius kestrel Falco punctatus; a tropical forest specialist that also occupies agricultural habitats. We specifically focused on the relationship between breeding success, agriculture and the timing of breeding because the latter is sensitive to changes in climatic conditions (spring rainfall), and enables us to explore the interactive effects of different (land and climate) drivers of environmental change. 3.Breeding success, measured as egg survival to fledging, declines seasonally in this population, but we found that the rate of this decline became increasingly rapid as the area of agriculture around a nest site increased. If the relationship between breeding success and agriculture was used in isolation to estimate the demographic impact of agriculture it would significantly under-estimate breeding success in dry (early) springs, and over-estimate breeding success in wet (late) springs. 4.Analysis of prey delivered to nests suggests that the relationship between breeding success and agriculture might be due, in part, to spatial variation in the availability of native, arboreal geckos. 5.Synthesis and applications. Agriculture modifies the seasonal decline in breeding success in this population. As springs are becoming wetter in our study area and since the kestrels breed later in wetter springs, the impact of agriculture on breeding success will become worse over time. Our results suggest that forest restoration designed to reduce the detrimental impacts of agriculture on breeding may also help reduce the detrimental effects of breeding late due to wetter springs. Our results therefore highlight the importance of considering the interactive effects of environmental change when managing wild populations.

Use of mutants to dissect the role of ethylene signalling in organ senescence and the regulation of yield in Arabidopsis thaliana

The role of ethylene in regulating organ senescence in Arabidopsis has been investigated by studying the development of mutants that have an attenu- ated capacity to perceive the gas. The onset of leaf senescence and floral organ abscission was delayed in the ethylene-insensitive mutant etr1. The photosynthetic life span of rosette leaves was similarly extended in the gain- of-function mutant ers2, and this mutant also exhibited a delay in the timing of pod dehiscence primarily as a con- sequence of an extension in the final stages of senescence. A detailed analysis of yield revealed that whilst thousand grain weight was increased, by as much as 20 %, in etr1, ein4, and the loss-of-function mutant etr2, only the latter showed a significant increase in total weight of seeds produced per plant. The other studied mutants exhibited a reduction in total seed yield of almost 40 %. These observations are discussed in the context of the possible role of ethylene in regulating organ senescence and their significance in the breeding of crop plants with enhanced phenotypic characteristics.

Development of the Surface Urban Energy and Water balance Scheme (SUEWS) for cold climate cities

The Surface Urban Energy and Water Balance Scheme (SUEWS) is developed to include snow. The processes addressed include accumulation of snow on the different urban surface types: snow albedo and density aging, snow melting and re-freezing of meltwater. Individual model parameters are assessed and independently evaluated using long-term observations in the two cold climate cities of Helsinki and Montreal. Eddy covariance sensible and latent heat fluxes and snow depth observations are available for two sites in Montreal and one in Helsinki. Surface runoff from two catchments (24 and 45 ha) in Helsinki and snow properties (albedo and density) from two sites in Montreal are also analysed. As multiple observation sites with different land-cover characteristics are available in both cities, model development is conducted independent of evaluation. The developed model simulates snowmelt related runoff well (within 19% and 3% for the two catchments in Helsinki when there is snow on the ground), with the springtime peak estimated correctly. However, the observed runoff peaks tend to be smoother than the simulated ones, likely due to the water holding capacity of the catchments and the missing time lag between the catchment and the observation point in the model. For all three sites the model simulates the timing of the snow accumulation and melt events well, but underestimates the total snow depth by 18–20% in Helsinki and 29–33% in Montreal. The model is able to reproduce the diurnal pattern of net radiation and turbulent fluxes of sensible and latent heat during cold snow, melting snow and snow-free periods. The largest model uncertainties are related to the timing of the melting period and the parameterization of the snowmelt. The results show that the enhanced model can simulate correctly the exchange of energy and water in cold climate cities at sites with varying surface cover.

Comparing projections of future changes in runoff from hydrological and biome models in ISI-MIP

Future changes in runoff can have important implications for water resources and flooding. In this study, runoff projections from ISI-MIP (Inter-sectoral Impact Model Inter-comparison Project) simulations forced with HadGEM2-ES bias-corrected climate data under the Representative Concentration Pathway 8.5 have been analysed for differences between impact models. Projections of change from a baseline period (1981-2010) to the future (2070-2099) from 12 impacts models which contributed to the hydrological and biomes sectors of ISI-MIP were studied. The biome models differed from the hydrological models by the inclusion of CO2 impacts and most also included a dynamic vegetation distribution. The biome and hydrological models agreed on the sign of runoff change for most regions of the world. However, in West Africa, the hydrological models projected drying, and the biome models a moistening. The biome models tended to produce larger increases and smaller decreases in regionally averaged runoff than the hydrological models, although there is large inter-model spread. The timing of runoff change was similar, but there were differences in magnitude, particularly at peak runoff. The impact of vegetation distribution change was much smaller than the projected change over time, while elevated CO2 had an effect as large as the magnitude of change over time projected by some models in some regions. The effect of CO2 on runoff was not consistent across the models, with two models showing increases and two decreases. There was also more spread in projections from the runs with elevated CO2 than with constant CO2. The biome models which gave increased runoff from elevated CO2 were also those which differed most from the hydrological models. Spatially, regions with most difference between model types tended to be projected to have most effect from elevated CO2, and seasonal differences were also similar, so elevated CO2 can partly explain the differences between hydrological and biome model runoff change projections. Therefore, this shows that a range of impact models should be considered to give the full range of uncertainty in impacts studies.

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.

Grammatical gender and the notion of default: insights from language acquisition

The aim of this study is to investigate the contrast in the timing of acquisition of grammatical gender attested in Dutch and Greek child learners. Greek children show precocious acquisition of neuter gender in particular, while Dutch children experience a long delay in the acquisition of neuter nouns, which extends to school age. For both Dutch and Greek, neuter has been claimed to be the default gender value on grounds of syntactic distribution in contexts where gender agreement is inert. To reconcile the contrast between the learner and the language facts in Dutch, as well as the contrast in the timing between Greek and Dutch monolingual child learners, we consider two sets of criteria to define the notion of default: one set pertains to the notion of linguistic default and the other to the notion of learner default. We suggest that, whereas Greek neuter is both the linguistic and the learner default value, Dutch neuter is the linguistic but not the learner default, leading to a learnability problem.

Early, late or very late? Timing acquisition and bilingualism

Research on child bilingualism accounts for differences in the course and the outcomes of monolingual and different types of bilingual language acquisition primarily from two perspectives: age of onset of exposure to the language(s) and the role of the input (Genesee, Paradis, & Crago, 2004; Meisel, 2009; Unsworth et al., 2014). Some findings suggest that early successive bilingual children may pattern similarly to simultaneous bilingual children, passing through different trajectories from child L2 learners due to a later age of onset in the latter group. Studies on bilingual development have also shown that input quantity in bilingual acquisition is considerably reduced, i.e., in each of their two languages, bilingual children are likely exposed to much less input than their monolingual peers (Paradis & Genesee, 1996; Unsworth, 2013b). At the same time, simultaneous bilingual children develop and attain competence in the two languages, sometimes without even an attested age delay compared to monolingual children (Paradis, Genesee & Crago, 2011). The implication is that even half of the input suffices for early language development, at least with respect to ‘core’ aspects of language, in whatever way ‘core’ is defined.My aim in this article is to consider how an additional, linguistic variable interacts with age of onset and input in bilingual development, namely, the timing in L1 development of the phenomena examined in bilingual children’s performance. Specifically, I will consider timing differences attested in the monolingual development of features and structures, distinguishing between early, late or ‘very late’ acquired phenomena. I will then argue that this three-way distinction reflects differences in the role of narrow syntax: early phenomena are core, parametric and narrowly syntactic, in contrast to late and very late phenomena, which involve syntax-external or even language-external resources too. I explore the consequences of these timing differences in monolingual development for bilingual development. I will review some findings from early (V2 in Germanic, grammatical gender in Greek), late (passives) and very late (grammatical gender in Dutch) phenomena in the bilingual literature and argue that early phenomena can differentiate between simultaneous and (early) successive bilingualism with an advantage for the former group, while the other two reveal similarly (high or low) performance across bilingual groups, differentiating them from monolinguals. The paper proposes that questions about the role of age of onset and language input in early bilingual development can only be meaningfully addressed when the properties and timing of the phenomena under investigation are taken into account.

Aerosol microphysics simulations of the Mt.~Pinatubo eruption with the UM-UKCA composition-climate model

We use a stratosphere–troposphere composition–climate model with interactive sulfur chemistry and aerosol microphysics, to investigate the effect of the 1991 Mount Pinatubo eruption on stratospheric aerosol properties. Satellite measurements indicate that shortly after the eruption, between 14 and 23 Tg of SO2 (7 to 11.5 Tg of sulfur) was present in the tropical stratosphere. Best estimates of the peak global stratospheric aerosol burden are in the range 19 to 26 Tg, or 3.7 to 6.7 Tg of sulfur assuming a composition of between 59 and 77 % H2SO4. In light of this large uncertainty range, we performed two main simulations with 10 and 20 Tg of SO2 injected into the tropical lower stratosphere. Simulated stratospheric aerosol properties through the 1991 to 1995 period are compared against a range of available satellite and in situ measurements. Stratospheric aerosol optical depth (sAOD) and effective radius from both simulations show good qualitative agreement with the observations, with the timing of peak sAOD and decay timescale matching well with the observations in the tropics and mid-latitudes. However, injecting 20 Tg gives a factor of 2 too high stratospheric aerosol mass burden compared to the satellite data, with consequent strong high biases in simulated sAOD and surface area density, with the 10 Tg injection in much better agreement. Our model cannot explain the large fraction of the injected sulfur that the satellite-derived SO2 and aerosol burdens indicate was removed within the first few months after the eruption. We suggest that either there is an additional alternative loss pathway for the SO2 not included in our model (e.g. via accommodation into ash or ice in the volcanic cloud) or that a larger proportion of the injected sulfur was removed via cross-tropopause transport than in our simulations. We also critically evaluate the simulated evolution of the particle size distribution, comparing in detail to balloon-borne optical particle counter (OPC) measurements from Laramie, Wyoming, USA (41° N). Overall, the model captures remarkably well the complex variations in particle concentration profiles across the different OPC size channels. However, for the 19 to 27 km injection height-range used here, both runs have a modest high bias in the lowermost stratosphere for the finest particles (radii less than 250 nm), and the decay timescale is longer in the model for these particles, with a much later return to background conditions. Also, whereas the 10 Tg run compared best to the satellite measurements, a significant low bias is apparent in the coarser size channels in the volcanically perturbed lower stratosphere. Overall, our results suggest that, with appropriate calibration, aerosol microphysics models are capable of capturing the observed variation in particle size distribution in the stratosphere across both volcanically perturbed and quiescent conditions. Furthermore, additional sensitivity simulations suggest that predictions with the models are robust to uncertainties in sub-grid particle formation and nucleation rates in the stratosphere.

The response of ionospheric convection in the polar cap to substorm activity

We report multi-instrument observations during an isolated substorm on 17 October 1989. The EISCAT radar operated in the SP-UK-POLI mode measuring ionospheric convection at latitudes 71°-78°. SAMNET and the EISCAT Magnetometer Cross provide information on the timing of substorm expansion phase onset and subsequent intensifications, as well as the location of the field aligned and ionospheric currents associated with the substorm current wedge. IMP-8 magnetic field data are also included. Evidence of a substorm growth phase is provided by the equatorward motion of a flow reversal boundary across the EISCAT radar field of view at 2130 MLT, following a southward turning of the interplanetary magnetic field (IMF). We infer that the polar cap expanded as a result of the addition of open magnetic flux to the tail lobes during this interval. The flow reversal boundary, which is a lower limit to the polar cap boundary, reached an invariant latitude equatorward of 71° by the time of the expansion phase onset. A westward electrojet, centred at 65.4°, occurred at the onset of the expansion phase. This electrojet subsequently moved poleward to a maximum of 68.1° at 2000 UT and also widened. During the expansion phase, there is evidence of bursts of plasma flow which are spatially localised at longitudes within the substorm current wedge and which occurred well poleward of the westward electrojet. We conclude that the substorm onset region in the ionosphere, defined by the westward electrojet, mapped to a part of the tail radially earthward of the boundary between open and closed magnetic flux, the “distant” neutral line. Thus the substorm was not initiated at the distant neutral line, although there is evidence that it remained active during the expansion phase. It is not obvious whether the electrojet mapped to a near-Earth neutral line, but at its most poleward, the expanded electrojet does not reach the estimated latitude of the polar cap boundary.

Diet and lifestyle in Bronze Age Northwest Spain: the collective burial of Cova do Santo

A multidisciplinary investigation of the collective burial of Cova do Santo is presented as a novel approach to understand daily life during the Bronze Age in Northwest Iberia. The research is focused on three main aspects: i) taphonomy and patterns of disposal, ii) paleopathology and -demography as indicators of health status and lifestyle, and iii) stable isotope analysis to reconstruct paleodiet and to investigate the timing of the introduction of millet to the Iberian Peninsula. Osteological analyses were performed on 64 bones (61 human and 3 animal); additionally, bone collagen was extracted from 15 samples (13 human and 2 animal) and analyzed for its carbon and nitrogen stable isotopes composition. The radiocarbon age of the human remains is consistent with the Middle Bronze Age (c. 1890 to 1600 cal BC). The recovered remains belonged to a minimum number of 14 individuals with an estimated age at death of forty years or younger. This relatively young age is in contrast to a high prevalence of degenerative joint disease in the group. The isotopic results suggest a very homogeneous diet, which was almost exclusively based on C3 plants and terrestrial animal products. Overall, the data suggest that the studied population belonged to a period prior to the introduction of spring or summer-grown crops such as millets. The collective burial from the cave of Cova de Santo, Galicia, currently represents the largest assemblage of prehistoric human remains from Northwest Spain and the relatively good preservation of the bones offers a unique opportunity to investigate daily life in Northern Iberia during the Bronze Age.

Propagation of the Madden–Julian Oscillation and scale interaction with the diurnal cycle in a high-resolution GCM

The Madden–Julian Oscillation (MJO) is the chief source of tropical intra-seasonal variability, but is simulated poorly by most state-of-the-art GCMs. Common errors include a lack of eastward propagation at the correct frequency and zonal extent, and too small a ratio of eastward- to westward-propagating variability. Here it is shown that HiGEM, a high-resolution GCM, simulates a very realistic MJO with approximately the correct spatial and temporal scale. Many MJO studies in GCMs are limited to diagnostics which average over a latitude band around the equator, allowing an analysis of the MJO’s structure in time and longitude only. In this study a wider range of diagnostics is applied. It is argued that such an approach is necessary for a comprehensive analysis of a model’s MJO. The standard analysis of Wheeler and Hendon (Mon Wea Rev 132(8):1917–1932, 2004; WH04) is applied to produce composites, which show a realistic spatial structure in the MJO envelopes but for the timing of the peak precipitation in the inter-tropical convergence zone, which bifurcates the MJO signal. Further diagnostics are developed to analyse the MJO’s episodic nature and the “MJO inertia” (the tendency to remain in the same WH04 phase from one day to the next). HiGEM favours phases 2, 3, 6 and 7; has too much MJO inertia; and dies out too frequently in phase 3. Recent research has shown that a key feature of the MJO is its interaction with the diurnal cycle over the Maritime Continent. This interaction is present in HiGEM but is unrealistically weak.