Evaluation of modelled regional water balance using lake status data: a comparison of 6ka simulations with the NCAR CCM

Variations in lake area and depth reflect climatically induced changes in the water balance of overflowing as well as closed lakes. A new global data base of lake status has been assembled, and is used to compare two simulations for 6 ka (6000 yr ago) made with successive R15 versions of the NCAR Community Climate Model (CCM). Simulated water balance was expressed as anomalies of annual precipitation minus evaporation (P-E); observed water balance as anomalies of lake status. Comparisons were made visually, by comparing regional averages, and by a statistic that compares the signs of simulated P-E anomalies (smoothly interpolated to the lake sites) with the status anomalies. Both CCM0 and CCM1 showed enhanced Northern-Hemisphere monsoons at 6 ka. Both underestimated the effect, but CCM1 fitted the spatial patterns better. In the northern mid- and high-latitudes the two versions differed more, and fitted the data less satisfactorily. CCM1 performed better than CCM0 in North America and central Eurasia, but not in Europe. Both models (especially CCM0) simulated excessive aridity in interior Eurasia. The models were systematically wrong in the southern mid-latitudes. Problems may have been caused by inadequate treatment of changes in sea-surface conditions in both models. Palaeolake status data will continue to provide a benchmark for the evaluation of modelling improvements.

Holocene alluvial environments at Barking, Lower Thames Valley (London, UK)

Deposit modelling based on archived borehole logs supplemented by a small number of dedicated boreholes is used to reconstruct the main boundary surfaces and the thickness of the main sediment units within the succession of Holocene alluvial deposits underlying the floodplain in the Barking Reach of the Lower Thames Valley. The basis of the modelling exercise is discussed and the models are used to assess the significance of floodplain relief in determining patterns of sedimentation. This evidence is combined with the results of biostratigraphical and geochronological investigations to reconstruct the environmental conditions associated with each successive stage of floodplain aggradation. The two main factors affecting the history and spatial pattern of Holocene sedimentation are shown to be the regional behaviour of relative sea level and the pattern of relief on the surface of the sub-alluvial, Late Devensian Shepperton Gravel. As is generally the case in the Lower Thames Valley, three main stratigraphic units are recognised, the Lower Alluvium, a peat bed broadly equivalent to the Tilbury III peat of Devoy (1979) and an Upper Alluvium. There is no evidence to suggest that the floodplain was substantially re-shaped by erosion during the Holocene. Instead, the relief inherited from the Shepperton Gravel surface was gradually buried either by the accumulation of peat or by deposition of fine-grained sediment from suspension in standing or slow-moving water. The palaeoenvironmental record from Barking confirms important details of the Holocene record observed elsewhere in the Lower Thames Valley, including the presence of Taxus in the valley-floor fen carr woodland between about 5000 and 4000 cal BP, and the subsequent growth of Ulmus on the peat surface.

A refined parametric model for short term load forecasting

Abstract We present a refined parametric model for forecasting electricity demand which performed particularly well in the recent Global Energy Forecasting Competition (GEFCom 2012). We begin by motivating and presenting a simple parametric model, treating the electricity demand as a function of the temperature and day of the data. We then set out a series of refinements of the model, explaining the rationale for each, and using the competition scores to demonstrate that each successive refinement step increases the accuracy of the model’s predictions. These refinements include combining models from multiple weather stations, removing outliers from the historical data, and special treatments of public holidays.

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.

Cusp ion steps, field-aligned currents and poleward moving auroral forms

We predict the field-aligned currents around cusp ion steps produced by pulsed reconnection between the geomagnetic field and an interplanetary magnetic field (IMF) with a B-Y component that is large in magnitude. For B-Y > 0, patches of newly opened flux move westward and eastward in the Northern and Southern Hemispheres, respectively, under the influence of the magnetic curvature force. These flow directions are reversed for B-Y < 0. The speed of this longitudinal motion initially grows with elapsed time since reconnection, but then decays as the newly opened field lines straighten. We predict sheets of field-aligned current on the boundaries between the patches produced by successive reconnection pulses, associated with the difference in the speeds of their longitudinal motion. For low elapsed times since reconnection, near the equatorward edge of the cusp region where the field lines are accelerating, the field-aligned current sheets will be downward or upward in both hemispheres for positive or negative IMF B-Y, respectively. At larger elapsed times since reconnection, as events slow and evolve from the cusp into the mantle region, these field-aligned current directions will be reversed. Observations by the Polar spacecraft on August 26,1998, show the predicted upward current sheets at steps seen in the mantle region for IMF B-Y > 0. Mapped into the ionosphere, the steps coincide with poleward moving events seen by the CUTLASS HF radar. The mapped location of the largest step also coincides with a poleward moving arc seen by the UVI imager on Polar. We show that the arc is consistent with a region of upward field-aligned current that has become unstable, such that a potential drop of about 1 kV formed below the spacecraft. The importance of these observations is that they confirm that the poleward moving events, as seen by the HF radar and the UV imager, are due to pulsed magnetopause reconnection. Milan et al. [2000] noted that the great longitudinal extent of these events means that the required reconnection pulses would have contributed almost all the voltage placed across the magnetosphere at this time. The observations also show that auroral arcs can form on open field lines in response to the pulsed application of voltage at the magnetopause.

Nonlinear regional warming with increasing CO₂ concentration

When considering adaptation measures and global climate mitigation goals, stakeholders need regional-scale climate projections, including the range of plausible warming rates. To assist these stakeholders, it is important to understand whether some locations may see disproportionately high or low warming from additional forcing above targets such as 2 K (ref. 1). There is a need to narrow uncertainty2 in this nonlinear warming, which requires understanding how climate changes as forcings increase from medium to high levels. However, quantifying and understanding regional nonlinear processes is challenging. Here we show that regional-scale warming can be strongly superlinear to successive CO2 doublings, using five different climate models. Ensemble-mean warming is superlinear over most land locations. Further, the inter-model spread tends to be amplified at higher forcing levels, as nonlinearities grow—especially when considering changes per kelvin of global warming. Regional nonlinearities in surface warming arise from nonlinearities in global-mean radiative balance, the Atlantic meridional overturning circulation, surface snow/ice cover and evapotranspiration. For robust adaptation and mitigation advice, therefore, potentially avoidable climate change (the difference between business-as-usual and mitigation scenarios) and unavoidable climate change (change under strong mitigation scenarios) may need different analysis methods.

Precise analysis of array usage in scientific programs

The automatic transformation of sequential programs for efficient execution on parallel computers involves a number of analyses and restructurings of the input. Some of these analyses are based on computing array sections, a compact description of a range of array elements. Array sections describe the set of array elements that are either read or written by program statements. These sections can be compactly represented using shape descriptors such as regular sections, simple sections, or generalized convex regions. However, binary operations such as Union performed on these representations do not satisfy a straightforward closure property, e.g., if the operands to Union are convex, the result may be nonconvex. Approximations are resorted to in order to satisfy this closure property. These approximations introduce imprecision in the analyses and, furthermore, the imprecisions resulting from successive operations have a cumulative effect. Delayed merging is a technique suggested and used in some of the existing analyses to minimize the effects of approximation. However, this technique does not guarantee an exact solution in a general setting. This article presents a generalized technique to precisely compute Union which can overcome these imprecisions.

Rainfall estimation by rain gauge-radar combination: a concurrent multiplicative-additive approach

A procedure (concurrent multiplicative-additive objective analysis scheme [CMA-OAS]) is proposed for operational rainfall estimation using rain gauges and radar data. On the basis of a concurrent multiplicative-additive (CMA) decomposition of the spatially nonuniform radar bias, within-storm variability of rainfall and fractional coverage of rainfall are taken into account. Thus both spatially nonuniform radar bias, given that rainfall is detected, and bias in radar detection of rainfall are handled. The interpolation procedure of CMA-OAS is built on Barnes' objective analysis scheme (OAS), whose purpose is to estimate a filtered spatial field of the variable of interest through a successive correction of residuals resulting from a Gaussian kernel smoother applied on spatial samples. The CMA-OAS, first, poses an optimization problem at each gauge-radar support point to obtain both a local multiplicative-additive radar bias decomposition and a regionalization parameter. Second, local biases and regionalization parameters are integrated into an OAS to estimate the multisensor rainfall at the ground level. The procedure is suited to relatively sparse rain gauge networks. To show the procedure, six storms are analyzed at hourly steps over 10,663 km2. Results generally indicated an improved quality with respect to other methods evaluated: a standard mean-field bias adjustment, a spatially variable adjustment with multiplicative factors, and ordinary cokriging.

The ecological impact of conquest and colonization on a medieval frontier landscape: combined palynological and geochemical analysis of lake sediments from Radzyń Chełminski, northern Poland

Slavic and German colonization of the southern Baltic between the 8th and 15th centuries A.D. is well-documented archaeologically and historically. Despite the large number of pollen profiles from Poland, few palaeoecological studies have examined the ecological impact of a process that was central to the expansion of European, Christian, societies. This study aims to redress this balance through multiproxy analysis of lake sediments from Radzyń Chełminski, Northern Poland, using pollen, element geochemistry (Inductively Coupled-Optical Emission Spectroscopy [ICP-OES]), organic content, and magnetic susceptibility. The close association between lake and medieval settlements presents the ideal opportunity to reconstruct past vegetation and land-use dynamics within a well-documented archaeological, historical, and cultural context. Three broad phases of increasing landscape impact are visible in the pollen and geochemical data dating from the 8th/9th, 10th/11th, and 13th centuries, reflecting successive phases of Slavic and German colonization. This involved the progressive clearance of oak-hornbeam dominated woodland and the development of an increasingly open agricultural landscape. Although the castles and towns of the Teutonic Order remain the most visible signs of medieval colonization, the palynological and geochemical data demonstrate that the major phase of woodland impact occurred during the preceding phase of Slavic expansion; Germans colonists were entering a landscape already significantly altered.

A model analysis of climate and CO2 controls on tree growth in a semi-arid woodland

We used a light-use efficiency model of photosynthesis coupled with a dynamic carbon allocation and tree-growth model to simulate annual growth of the gymnosperm Callitris columellaris in the semi-arid Great Western Woodlands, Western Australia, over the past 100 years. Parameter values were derived from independent observations except for sapwood specific respiration rate, fine-root turnover time, fine-root specific respiration rate and the ratio of fine-root mass to foliage area, which were estimated by Bayesian optimization. The model reproduced the general pattern of interannual variability in radial growth (tree-ring width), including the response to the shift in precipitation regimes that occurred in the 1960s. Simulated and observed responses to climate were consistent. Both showed a significant positive response of tree-ring width to total photosynthetically active radiation received and to the ratio of modeled actual to equilibrium evapotranspiration, and a significant negative response to vapour pressure deficit. However, the simulations showed an enhancement of radial growth in response to increasing atmospheric CO2 concentration (ppm) ([CO2]) during recent decades that is not present in the observations. The discrepancy disappeared when the model was recalibrated on successive 30-year windows. Then the ratio of fine-root mass to foliage area increases by 14% (from 0.127 to 0.144 kg C m-2) as [CO2] increased while the other three estimated parameters remained constant. The absence of a signal of increasing [CO2] has been noted in many tree-ring records, despite the enhancement of photosynthetic rates and water-use efficiency resulting from increasing [CO2]. Our simulations suggest that this behaviour could be explained as a consequence of a shift towards below-ground carbon allocation.

National experiences with insurgencies and counterinsurgencies: are there “National styles”?

There are clearly identifiable patterns in the way in which insurgents operate in certain geographic areas and cultures, and in which states deal with them. These have been stressed by writing about national "ways of war", strategic culture, or national styles. Nevertheless, there have been important ruptures and changes in some of these, so that for Britain, for example, three successive phases corresponding to patterns can be identified. For France, two pronounced rivaling traditions coincided over two centuries. Algerians by contrast changed fundamentally in their fighting style with the different political ideologies they were following. Palestinian insurgency against Israel is also marked by change in approach. Russia and China possibly show the longest continuity in their handling of insurgencies.

Orbital and millenial-scale environmental changes between 64 and 25 ka BP recorded in Black Sea sediments

High-resolution pollen and dinoflagellate cyst records from sediment core M72/5-25-GC1 were used to reconstruct vegetation dynamics in northern Anatolia and surface conditions of the Black Sea between 64 and 20 ka BP. During this period, the dominance of Artemisia in the pollen record indicates a steppe landscape and arid climate conditions. However, the concomitant presence of temperate arboreal pollen suggests the existence of glacial refugia in northern Anatolia. Long-term glacial vegetation dynamics reveal two major arid phases ~64–55 and 40–32 ka BP, and two major humid phases ~54–45 and 28–20 ka BP, correlating with higher and lower summer insolation, respectively. Dansgaard–Oeschger (D–O) cycles are clearly indicated by the 25-GC1 pollen record. Greenland interstadials are characterized by a marked increase in temperate tree pollen, indicating a spread of forests due to warm/wet conditions in northern Anatolia, whereas Greenland stadials reveal cold and arid conditions as indicated by spread of xerophytic biomes. There is evidence for a phase lag of ~500 to 1500 yr between initial warming and forest expansion, possibly due to successive changes in atmospheric circulation in the North Atlantic sector. The dominance of Pyxidinopsis psilata and Spiniferites cruciformis in the dinocyst record indicates brackish Black Sea conditions during the entire glacial period. The decrease of marine indicators (marine dinocysts, acritarchs) at ~54 ka BP and increase of freshwater algae (Pediastrum, Botryococcus) from 32 to 25 ka BP reveals freshening of the Black Sea surface water. This freshening is possibly related to humid phases in the region, to connection between Caspian Sea and Black Sea, to seasonal freshening by floating ice, and/or to closer position of river mouths due to low sea level. In the southern Black Sea, Greenland interstadials are clearly indicated by high dinocyst concentrations and calcium carbonate content, as a result of an increase in primary productivity. Heinrich events show a similar impact on the environment in the northern Anatolia/Black Sea region as Greenland stadials.

Transmathematical basis of infinitely scalable pipeline machines

We describe infinitely scalable pipeline machines with perfect parallelism, in the sense that every instruction of an inline program is executed, on successive data, on every clock tick. Programs with shared data effectively execute in less than a clock tick. We show that pipeline machines are faster than single or multi-core, von Neumann machines for sufficiently many program runs of a sufficiently time consuming program. Our pipeline machines exploit the totality of transreal arithmetic and the known waiting time of statically compiled programs to deliver the interesting property that they need no hardware or software exception handling.

Hydrocarbon biodegradation and soil microbial community response to repeated oil exposure

Bioremediation strategies continue to be developed to mitigate the environmental impact of petroleum hydrocarbon contamination. This study investigated the ability of soil microbiota, adapted by prior exposure, to biodegrade petroleum. Soils from Barrow Is. (W. Australia), a class A nature reserve and home to Australia’s largest onshore oil field, were exposed to Barrow production oil (50 ml/kg soil) and incubated (25 °C) for successive phases of 61 and 100 days. Controls in which oil was not added at Phase I or II were concurrently studied and all treatments were amended with the same levels of additional nutrient and water to promote microbial activity. Prior exposure resulted in accelerated biodegradation of most, but not all, hydrocarbon constituents in the production oil. Molecular biodegradation parameters measured using gas chromatography–mass spectrometry (GC–MS) showed that several aromatic constituents were degraded more slowly with increased oil history. The unique structural response of the soil microbial community was reflected by the response of different phospholipid fatty acid (PLFA) sub-classes (e.g. branched saturated fatty acids of odd or even carbon number) measured using a ratio termed Barrow PLFA ratio (B-PLFAr). The corresponding values of a previously proposed hydrocarbon degrading alteration index showed a negative correlation with hydrocarbon exposure, highlighting the site specificity of PLFA-based ratios and microbial community dynamics. B-PLFAr values increased with each Phase I and II addition of production oil. The different hydrocarbon biodegradation rates and responses of PLFA subclasses to the Barrow production oil probably relate to the relative bioavailability of production oil hydrocarbons. These different effects suggest preferred structural and functional microbial responses to anticipated contaminants may potentially be engineered by controlled pre-exposure to the same or closely related substrates. The bioremediation of soils freshly contaminated with petroleum could benefit from the addition of exhaustively bioremediated soils rich in biota primed for the impacting hydrocarbons.