A nonparametric ensemble transform method for Bayesian inference

Many applications, such as intermittent data assimilation, lead to a recursive application of Bayesian inference within a Monte Carlo context. Popular data assimilation algorithms include sequential Monte Carlo methods and ensemble Kalman filters (EnKFs). These methods differ in the way Bayesian inference is implemented. Sequential Monte Carlo methods rely on importance sampling combined with a resampling step, while EnKFs utilize a linear transformation of Monte Carlo samples based on the classic Kalman filter. While EnKFs have proven to be quite robust even for small ensemble sizes, they are not consistent since their derivation relies on a linear regression ansatz. In this paper, we propose another transform method, which does not rely on any a priori assumptions on the underlying prior and posterior distributions. The new method is based on solving an optimal transportation problem for discrete random variables. © 2013, Society for Industrial and Applied Mathematics

Settlement patterns in the late Mesolithic of western Scotland: the implications of Bayesian analysis of radiocarbon dates and inter-site technological comparisons

The character of settlement patterns within the late Mesolithic communities of north-west Europe is a topic of substantial debate. An important case study concerns the five shell middens on the island of Oronsay, Inner Hebrides, western Scotland. Two conflicting interpretations have been proposed: the evidence from seasonality indicators and stable isotope analysis of human bones has been used to support a model of year-round settlement on this small island; alternatively, the middens have been interpreted as resulting from short-term intermittent visits to Oronsay within a regionally mobile settlement pattern. We contribute to this debate by describing Storakaig, a newly discovered site on the nearby island of Islay, undertaking a Bayesian chronological analysis and providing evidence for technological continuity between Oronsay and sites elsewhere in the region. While this new evidence remains open to alternative interpretation, we suggest that it makes regional mobility rather than year-round settlement on Oronsay a more viable interpretation for the Oronsay middens. Our analysis also confirms the likely overlap of the late Mesolithic with the earliest Neolithic within western Scotland.

Change in convergence and accommodation after two weeks of eye exercises in typical young adults

Abstract: Introduction Although eye exercises appear to help heterophoria, convergence insufficiency and intermittent strabismus, true treatment effects can be confounded by placebo, practice and encouragement factors. This study assessed objective changes in vergence and accommodation responses in typical naïve young adults after two weeks of exercises compared to control conditions to assess the extent of treatment effects occur above other factors. Methods 156 asymptomatic young adults were randomly assigned to 6 exercise groups or 2 no-treatment groups. Treatment targeted i) accommodation, ii)vergence, iii) both, iv) convergence>accommodation, v)accommodation>convergence, or vi) a placebo. All were re-tested under identical conditions, except for the second control group who were additionally encouraged during testing. Objective accommodation and vergence were assessed to a range of targets moving in depth containing combinations of blur, disparity and proximity/looming cues. Results Response gain improved more for less naturalistic targets where more improvement was possible. Convergence exercises improved vergence for near across all targets (P=.035). Mean accommodation changed similarly,but non-significantly. No other treatment group differed significantly from the non-encouraged control group, while encouraging effort produced significantly increased vergence (P=.004) and accommodation (P=.005) gains in the other control group. Conclusions True treatment effects were small, only significantly better after vergence exercises to a non-accommodative target, and were rarely related to response they were designed to improve. Exercising accommodation without convergence made no difference to accommodation to cues containing detail. Additional effort improved objective responses the most, so should be controlled carefully in research, and considered when auditing treatment.

A nonlinear oscillator describing storm track variability

We construct a two-variable model which describes the interaction between local baroclinicity and eddy heat flux in order to understand aspects of the variance in storm tracks. It is a heuristic model for diabatically forced baroclinic instability close to baroclinic neutrality. The two-variable model has the structure of a nonlinear oscillator. It exhibits some realistic properties of observed storm track variability, most notably the intermittent nature of eddy activity. This suggests that apparent threshold behaviour can be more accurately and succinctly described by a simple nonlinearity. An analogy is drawn with triggering of convective events.

Financial implications for companies running standby generators in support of a smart grid

Integrating renewable energy into built environments requires additional attention to the balancing of supply and demand due to their intermittent nature. Demand Side Response (DSR) has the potential to make money for organisations as well as support the System Operator as the generation mix changes. There is an opportunity to increase the use of existing technologies in order to manage demand. Company-owned standby generators are a rarely used resource; their maintenance schedule often accounts for a majority of their running hours. DSR encompasses a range of technologies and organisations; Sustainability First (2012) suggest that the System Operator (SO), energy supply companies, Distribution Network Operators (DNOs), Aggregators and Customers all stand to benefit from DSR. It is therefore important to consider impact of DSR measures to each of these stakeholders. This paper assesses the financial implications of organisations using existing standby generation equipment for DSR in order to avoid peak electricity charges. It concludes that under the current GB electricity pricing structure, there are several regions where running diesel generators at peak times is financially beneficial to organisations. Issues such as fuel costs, Carbon Reduction Commitment (CRC) charges, maintenance costs and electricity prices are discussed.

Ecophysiology of seed dormancy and the control of germination in early spring-flowering Galanthus nivalis and Narcissus pseudonarcissus (Amaryllidaceae)

Seed dormancy induction and alleviation in the winter-flowering moist temperate woodland species Galanthus nivalis and Narcissus pseudonarcissus are complex and poorly understood. Temperature, light and desiccation were investigated to elucidate their role in the germination ecophysiology of these species. Outdoor and laboratory experiments simulating different seasonal temperatures, seasonal durations, and temperature fluctuations; the presence of light during different seasons; and intermittent drying (during the summer period) over several ‘years’ investigated the importance of these factors in germination. Warm summer-like temperatures (20°C) were necessary for germination at subsequent cooler autumn-like temperatures (greatest at 15°C in G. nivalis and 10°C in N. pseudonarcissus). As the warm temperature duration increased so did germination at subsequent cooler temperatures; further germination occurred in subsequent ‘years’ at cooler temperatures following a second, and also third, warm period. Germination was significantly greater in darkness, particularly in G. nivalis. Dormancy increased with seed maturation period in G. nivalis, because seeds extracted from green capsules germinated more readily than those from yellow. Desiccation increased dormancy in an increasing proportion of N. pseudonarcissus seeds the later they were dried in ‘summer’. Seed viability was only slightly reduced by desiccation in N. pseudonarcissus but was poor and variable in G. nivalis. Shoot formation occurred both at the temperature at which germination was greatest and also if 5°C cooler. In summary, continuous hydration of seeds of both species during warm summer-like temperatures results in the gradual release of seed dormancy; thereafter, darkness and cooler temperatures promote germination. Cold temperatures, increased seed maturity (G. nivalis), and desiccation (N. pseudonarcissus) increase dormancy while light inhibits germination.

Motion of the dayside polar cap boundary during substorm cycles: II. Generation of poleward-moving events and polar cap patches by pulses in the magnetopause reconnection rate

Using data from the EISCAT (European Incoherent Scatter) VHF and CUTLASS (Co-operative UK Twin- Located Auroral Sounding System) HF radars, we study the formation of ionospheric polar cap patches and their relationship to the magnetopause reconnection pulses identified in the companion paper by Lockwood et al. (2005). It is shown that the poleward-moving, high-concentration plasma patches observed in the ionosphere by EISCAT on 23 November 1999, as reported by Davies et al. (2002), were often associated with corresponding reconnection rate pulses. However, not all such pulses generated a patch and only within a limited MLT range (11:00–12:00 MLT) did a patch result from a reconnection pulse. Three proposed mechanisms for the production of patches, and of the concentration minima that separate them, are analysed and evaluated: (1) concentration enhancement within the patches by cusp/cleft precipitation; (2) plasma depletion in the minima between the patches by fast plasma flows; and (3) intermittent injection of photoionisation-enhanced plasma into the polar cap. We devise a test to distinguish between the effects of these mechanisms. Some of the events repeat too frequently to apply the test. Others have sufficiently long repeat periods and mechanism (3) is shown to be the only explanation of three of the longer-lived patches seen on this day. However, effect (2) also appears to contribute to some events. We conclude that plasma concentration gradients on the edges of the larger patches arise mainly from local time variations in the subauroral plasma, via the mechanism proposed by Lockwood et al. (2000).

Auroral bright spot sequence near 1400 MLT: Coordinated optical and ion drift observations

Optical observations of a dayside auroral brightening sequence, by means of all-sky TV cameras and meridian scanning photometers, have been combined with EISCAT ion drift observations within the same invariant latitude-MLT sector. The observations were made during a January 1989 campaign by utilizing the high F region ion densities during the maximum phase of the solar cycle. The characteristic intermittent optical events, covering ∼300 km in east-west extent, move eastward (antisunward) along the poleward boundary of the persistent background aurora at velocities of ∼1.5 km s−1 and are associated with ion flows which swing from eastward to westward, with a subsequent return to eastward, during the interval of a few minutes when there is enhanced auroral emission within the radar field of view. The breakup of discrete auroral forms occurs at the reversal (negative potential) that forms between eastward plasma flow, maximizing near the persistent arc poleward boundary, and strong transient westward flow to the south. The reported events, covering a 35 min interval around 1400 MLT, are embedded within a longer period of similar auroral activity between 0830 (1200 MLT) and 1300 UT (1600 MLT). These observations are discussed in relation to recent models of boundary layer plasma dynamics and the associated magnetosphere-ionosphere coupling. The ionospheric events may correspond to large-scale wave like motions of the low-latitude boundary layer (LLBL)/plasma sheet (PS) boundary. On the basis of this interpretation the observed spot size, speed and repetition period (∼10 min) give a wavelength (the distance between spots) of ∼900 km in the present case. The events can also be explained as ionospheric signatures of newly opened flux tubes associated with reconnection bursts at the magnetopause near 1400 MLT. We also discuss these data in relation to random, patchy reconnection (as has recently been invoked to explain the presence of the sheathlike plasma on closed field lines in the LLBL). In view of the lack of IMF data, and the existing uncertainty on the location of the open-closed field line boundary relative to the optical events, an unambiguous discrimination between the different alternatives is not easily obtained.

Midday auroral breakup events and related energy and momentum transfer from the magnetosheath

Combined observations by meridian-scanning photometers, all-sky auroral TV camera and the EISCAT radar permitted a detailed analysis of the temporal and spatial development of the midday auroral breakup phenomenon and the related ionospheric ion flow pattern within the 71°–75° invariant latitude radar field of view. The radar data revealed dominating northward and westward ion drifts, of magnitudes close to the corresponding velocities of the discrete, transient auroral forms, during the two different events reported here, characterized by IMF |BY/BZ| < 1 and > 2, respectively (IMF BZ between −8 and −3 nT and BY > 0). The spatial scales of the discrete optical events were ∼50 km in latitude by ∼500 km in longitude, and their lifetimes were less than 10 min. Electric potential enhancements with peak values in the 30–50 kV range are inferred along the discrete arc in the IMF |BY/BZ| < 1 case from the optical data and across the latitudinal extent of the radar field of view in the |BY/BZ| > 2 case. Joule heat dissipation rates in the maximum phase of the discrete structures of ∼ 100 ergs cm−2 s−1 (0.1 W m−2) are estimated from the photometer intensities and the ion drift data. These observations combined with the additional characteristics of the events, documented here and in several recent studies (i.e., their quasi-periodic nature, their motion pattern relative to the persistent cusp or cleft auroral arc, the strong relationship with the interplanetary magnetic field and the associated ion drift/E field events and ground magnetic signatures), are considered to be strong evidence in favour of a transient, intermittent reconnection process at the dayside magnetopause and associated energy and momentum transfer to the ionosphere in the polar cusp and cleft regions. The filamentary spatial structure and the spectral characteristics of the optical signature indicate associated localized ˜1-kV potential drops between the magnetopause and the ionosphere during the most intense auroral events. The duration of the events compares well with the predicted characteristic times of momentum transfer to the ionosphere associated with the flux transfer event-related current tubes. It is suggested that, after this 2–10 min interval, the sheath particles can no longer reach the ionosphere down the open flux tube, due to the subsequent super-Alfvénic flow along the magnetopause, conductivities are lower and much less momentum is extracted from the solar wind by the ionosphere. The recurrence time (3–15 min) and the local time distribution (∼0900–1500 MLT) of the dayside auroral breakup events, combined with the above information, indicate the important roles of transient magnetopause reconnection and the polar cusp and cleft regions in the transfer of momentum and energy between the solar wind and the magnetosphere.

Kraichnan–Leith–Batchelor similarity theory and two-dimensional inverse cascades

We study the scaling properties and Kraichnan–Leith–Batchelor (KLB) theory of forced inverse cascades in generalized two-dimensional (2D) fluids (α-turbulence models) simulated at resolution 8192x8192. We consider α=1 (surface quasigeostrophic flow), α=2 (2D Euler flow) and α=3. The forcing scale is well resolved, a direct cascade is present and there is no large-scale dissipation. Coherent vortices spanning a range of sizes, most larger than the forcing scale, are present for both α=1 and α=2. The active scalar field for α=3 contains comparatively few and small vortices. The energy spectral slopes in the inverse cascade are steeper than the KLB prediction −(7−α)/3 in all three systems. Since we stop the simulations well before the cascades have reached the domain scale, vortex formation and spectral steepening are not due to condensation effects; nor are they caused by large-scale dissipation, which is absent. One- and two-point p.d.f.s, hyperflatness factors and structure functions indicate that the inverse cascades are intermittent and non-Gaussian over much of the inertial range for α=1 and α=2, while the α=3 inverse cascade is much closer to Gaussian and non-intermittent. For α=3 the steep spectrum is close to that associated with enstrophy equipartition. Continuous wavelet analysis shows approximate KLB scaling ℰ(k)∝k−2 (α=1) and ℰ(k)∝k−5/3 (α=2) in the interstitial regions between the coherent vortices. Our results demonstrate that coherent vortex formation (α=1 and α=2) and non-realizability (α=3) cause 2D inverse cascades to deviate from the KLB predictions, but that the flow between the vortices exhibits KLB scaling and non-intermittent statistics for α=1 and α=2.

Economic crises in England, 1270-1520: a statistical approach

The chapter examines how far medieval economic crises can be identified by analysing the residuals from a simultaneous equation model of the medieval English economy. High inflation, falls in gross domestic product and large intermittent changes in wage rates are all considered as potential indicators of crisis. Potential causal factors include bad harvests, wars and political instability. The chapter suggests that crises arose when a combination of different problems overwhelmed the capacity of government to address them. It may therefore be a mistake to look for a single cause of any crisis. The coincidence of separate problems is a more plausible explanation of many crises.

Effect of fat replacement by inulin on textural and structural properties of short dough biscuits

The aim of this study was to evaluate the effects of inulin as fat replacer on short dough biscuits and their corresponding doughs. A control formulation, with no replacement, and four formulations in which 10, 20, 30, and 40 % of shortening was replaced by inulin were studied. In the dough, shortening was observed surrounding flour components. At higher fat replacement levels, flour was more available for hydration leading to significant (P<0.05) harder doughs: from 2.76 (0.12)N in 10 % fat-replaced biscuits to 5.81 (1.56)N in 30 % fat-replaced ones. Biscuit structure was more continuous than dough structure. A continuous fat layer coated the matrix surface, where starch granules were embedded. In general, weight loss during baking and water activity decreased significantly (P<0.05) as fat replacement increased. Biscuit dimensions and aeration decreased when fat replacement increased, e.g., width gain was +1.20 mm in 10 fat-replaced biscuits and only +0.32 mm in 40 % fat-replaced ones. Panelist found biscuits with 20 % of fat replacement slightly harder than control biscuits. It can be concluded that shortening may be partially replaced, up to 20 %, with inulin. These low fat biscuits are similar than the control biscuits, and they can have additional health benefits derived from inulin presence.

Developments in the quality of treated greywater supplies for buildings, and associated user perception and acceptance

A manufactured aeration and nanofiltration MBR greywater system was tested during continuous operation at the University of Reading, to demonstrate reliability in delivery of high quality treated greywater. Its treatment performance was evaluated against British Standard criteria [BSI (Greywater Systems—Part 1 Code of Practice: BS8525-1:2010. BS Press, 2010); (Greywater Systems—Part 2 Domestic Greywater Treatment, Requirements and Methods: BS 8525-2:2011. BS Press, 2011)]. The low carbon greywater recycling technology produced excellent analytical results as well as consistency in performance. User acceptance of such reliably treated greywater was then evaluated through user perception studies. The results inform the potential supply of treated greywater to student accommodation. Out of 135 questionnaire replies, 95% demonstrated a lack of aversion in one or more attributes, to using treated, recycled greywater.

Eddies and variability in the Mozambique Channel

Between 1995 and 2000, on average 4 eddies per year are observed from satellite altimetry to propagate southward through the Mozambique Channel, into the upstream Agulhas region. Further south, these eddies have been found to control the timing and frequenc yof Agulhas ring shedding. Within the Mozambique Channel, anomalous SSH amplitudes rise to 30 cm ; in agreement with in situ measured velocities. Comparison of an observed velocit ysection with GCM model results shows that the Mozambique Channel eddies in these models are too surface intensified. Also, the number of eddies formed in the models is in disagreement with our observational analysis. Moored current meter measurements observing the passage of three eddies in 2000 are extended to a 5-year time series b yreferencing the anomalous surface currents estimated from altimeter data to a s ynoptic LADCP velocit y measurement. The results show intermittent edd ypassage at the mooring location. A statistical analysis of SSH observations in different parts of the Mozambique Channel shows a southward decrease of the dominant frequency of the variability, going from 7 per year in the extension of the South Equatorial Current north of Madagascar to 4 per year south of Madagascar. The observations suggest that frequency reduction is related to the Rossb ywaves coming in from the east

Fathoming solar effects in Neptune’s atmosphere

Long-duration observations of Neptune’s brightness in two visible wavelengths provide a disk-averaged estimate of its atmospheric aerosol. Brightness variations were previously associated with the 11-year solar cycle, through solar-modulated mechanisms linked with either ultra-violet (UV) or galactic cosmic ray (GCR) effects on atmospheric particles. Here we use a recently extended brightness dataset (1972-2014), with physically realistic modelling to show that rather than alternatives, UV and GCR are likely to be modulating Neptune’s atmosphere in combination. The importance of GCR is further supported by the response of Neptune's atmosphere to an intermittent 1.5 to 1.9 year periodicity, which occurred preferentially in GCR (not UV) during the mid-1980s. This periodicity was detected both at Earth, and in GCR measured by Voyager 2, then near Neptune. A similar coincident variability in Neptune’s brightness suggests nucleation onto GCR ions. Both GCR and UV mechanisms may occur more rapidly than the subsequent atmospheric particle transport.