Gaussian anamorphosis in the analysis step of the EnKF: a joint state-variable/observation approach

The analysis step of the (ensemble) Kalman filter is optimal when (1) the distribution of the background is Gaussian, (2) state variables and observations are related via a linear operator, and (3) the observational error is of additive nature and has Gaussian distribution. When these conditions are largely violated, a pre-processing step known as Gaussian anamorphosis (GA) can be applied. The objective of this procedure is to obtain state variables and observations that better fulfil the Gaussianity conditions in some sense. In this work we analyse GA from a joint perspective, paying attention to the effects of transformations in the joint state variable/observation space. First, we study transformations for state variables and observations that are independent from each other. Then, we introduce a targeted joint transformation with the objective to obtain joint Gaussianity in the transformed space. We focus primarily in the univariate case, and briefly comment on the multivariate one. A key point of this paper is that, when (1)-(3) are violated, using the analysis step of the EnKF will not recover the exact posterior density in spite of any transformations one may perform. These transformations, however, provide approximations of different quality to the Bayesian solution of the problem. Using an example in which the Bayesian posterior can be analytically computed, we assess the quality of the analysis distributions generated after applying the EnKF analysis step in conjunction with different GA options. The value of the targeted joint transformation is particularly clear for the case when the prior is Gaussian, the marginal density for the observations is close to Gaussian, and the likelihood is a Gaussian mixture.

Ensemble downscaling in coupled solar wind-magnetosphere modeling for space weather forecasting

Advanced forecasting of space weather requires simulation of the whole Sun-to-Earth system, which necessitates driving magnetospheric models with the outputs from solar wind models. This presents a fundamental difficulty, as the magnetosphere is sensitive to both large-scale solar wind structures, which can be captured by solar wind models, and small-scale solar wind “noise,” which is far below typical solar wind model resolution and results primarily from stochastic processes. Following similar approaches in terrestrial climate modeling, we propose statistical “downscaling” of solar wind model results prior to their use as input to a magnetospheric model. As magnetospheric response can be highly nonlinear, this is preferable to downscaling the results of magnetospheric modeling. To demonstrate the benefit of this approach, we first approximate solar wind model output by smoothing solar wind observations with an 8 h filter, then add small-scale structure back in through the addition of random noise with the observed spectral characteristics. Here we use a very simple parameterization of noise based upon the observed probability distribution functions of solar wind parameters, but more sophisticated methods will be developed in the future. An ensemble of results from the simple downscaling scheme are tested using a model-independent method and shown to add value to the magnetospheric forecast, both improving the best estimate and quantifying the uncertainty. We suggest a number of features desirable in an operational solar wind downscaling scheme.

The Environmental Data Abstraction Library (EDAL): a modular approach to processing and visualising large environmental data

The Environmental Data Abstraction Library provides a modular data management library for bringing new and diverse datatypes together for visualisation within numerous software packages, including the ncWMS viewing service, which already has very wide international uptake. The structure of EDAL is presented along with examples of its use to compare satellite, model and in situ data types within the same visualisation framework. We emphasize the value of this capability for cross calibration of datasets and evaluation of model products against observations, including preparation for data assimilation.

Predicting the orientation of lipid cubic phase films

Lipid cubic phase films are of increasingly widespread importance, both in the analysis of the cubic phases themselves by techniques including microscopy and X-ray scattering, and in their applications, especially as electrode coatings for electrochemical sensors and for templates for the electrodeposition of nanostructured metal. In this work we demonstrate that the crystallographic orientation adopted by these films is governed by minimization of interfacial energy. This is shown by the agreement between experimental data obtained using grazing-incidence small-angle X-ray scattering (GI-SAXS), and the predicted lowest energy orientation determined using a theoretical approach we have recently developed. GI-SAXS data show a high degree of orientation for films of both the double diamond phase and the gyroid phase, with the [111] and [110] directions respectively perpendicular to the planar substrate. In each case, this matches the lowest energy facet calculated for that particular phase.

Flexible selection of a single treatment incorporating short-term endpoint information in a phase II/III clinical trial

Seamless phase II/III clinical trials in which an experimental treatment is selected at an interim analysis have been the focus of much recent research interest. Many of the methods proposed are based on the group sequential approach. This paper considers designs of this type in which the treatment selection can be based on short-term endpoint information for more patients than have primary endpoint data available. We show that in such a case, the familywise type I error rate may be inflated if previously proposed group sequential methods are used and the treatment selection rule is not specified in advance. A method is proposed to avoid this inflation by considering the treatment selection that maximises the conditional error given the data available at the interim analysis. A simulation study is reported that illustrates the type I error rate inflation and compares the power of the new approach with two other methods: a combination testing approach and a group sequential method that does not use the short-term endpoint data, both of which also strongly control the type I error rate. The new method is also illustrated through application to a study in Alzheimer's disease. © 2015 The Authors. Statistics in Medicine Published by John Wiley & Sons Ltd.

Geometric polarimetry − part I: spinors and wave states

A new formal approach for representation of polarization states of coherent and partially coherent electromagnetic plane waves is presented. Its basis is a purely geometric construction for the normalised complex-analytic coherent wave as a generating line in the sphere of wave directions, and whose Stokes vector is determined by the intersection with the conjugate generating line. The Poincare sphere is now located in physical space, simply a coordination of the wave sphere, its axis aligned with the wave vector. Algebraically, the generators representing coherent states are represented by spinors, and this is made consistent with the spinor-tensor representation of electromagnetic theory by means of an explicit reference spinor we call the phase flag. As a faithful unified geometric representation, the new model provides improved formal tools for resolving many of the geometric difficulties and ambiguities that arise in the traditional formalism.

Aligned platinum nanowire networks from surface-oriented lipid cubic phase templates

Mesoporous metal structures featuring a bicontinuous cubic morphology have a wide range of potential applications and novel opto-electronic properties, often orientation-dependent. We describe the production of nanostructured metal films 1–2 microns thick featuring 3D-periodic ‘single diamond’ morphology that show high out-of-plane alignment, with the (111) plane oriented parallel to the substrate. These are produced by electrodeposition of platinum through a lipid cubic phase (QII) template. Further investigation into the mechanism for the orientation revealed the surprising result that the QII template, which is tens of microns thick, is polydomain with no overall orientation. When thicker platinum films are grown, they also show increased orientational disorder. These results suggest that polydomain QII samples display a region of uniaxial orientation at the lipid/substrate interface up to approximately 2.8 ± 0.3 μm away from the solid surface. Our approach gives previously unavailable information on the arrangement of cubic phases at solid interfaces, which is important for many applications of QII phases. Most significantly, we have produced a previously unreported class of oriented nanomaterial, with potential applications including metamaterials and lithographic masks.

An approximate dynamic programming approach for improving accuracy of lossy data compression by Bloom filters

Bloom filters are a data structure for storing data in a compressed form. They offer excellent space and time efficiency at the cost of some loss of accuracy (so-called lossy compression). This work presents a yes-no Bloom filter, which as a data structure consisting of two parts: the yes-filter which is a standard Bloom filter and the no-filter which is another Bloom filter whose purpose is to represent those objects that were recognised incorrectly by the yes-filter (that is, to recognise the false positives of the yes-filter). By querying the no-filter after an object has been recognised by the yes-filter, we get a chance of rejecting it, which improves the accuracy of data recognition in comparison with the standard Bloom filter of the same total length. A further increase in accuracy is possible if one chooses objects to include in the no-filter so that the no-filter recognises as many as possible false positives but no true positives, thus producing the most accurate yes-no Bloom filter among all yes-no Bloom filters. This paper studies how optimization techniques can be used to maximize the number of false positives recognised by the no-filter, with the constraint being that it should recognise no true positives. To achieve this aim, an Integer Linear Program (ILP) is proposed for the optimal selection of false positives. In practice the problem size is normally large leading to intractable optimal solution. Considering the similarity of the ILP with the Multidimensional Knapsack Problem, an Approximate Dynamic Programming (ADP) model is developed making use of a reduced ILP for the value function approximation. Numerical results show the ADP model works best comparing with a number of heuristics as well as the CPLEX built-in solver (B&B), and this is what can be recommended for use in yes-no Bloom filters. In a wider context of the study of lossy compression algorithms, our researchis an example showing how the arsenal of optimization methods can be applied to improving the accuracy of compressed data.

Phase 1

Exhibition investigating the role of different methods of drawing in the visualization of architectural space

Pasolini and Third World hunger: an approach to Cinema Novo through 'La ricotta'

Among the links between Pier Paolo Pasolini and Brazilian Cinema Novo, one of the most inspiring is the political approach to hunger and consumption. In this text, I analyse this topic to look at how some of the aesthetic ideas in Pasolini’s La ricotta (1963) can also be found in some of the most important films of Cinema Novo. In 'La ricotta' (1963), the irresistible need to eat of a subproletarian interacts and clashes with his responsibilities as an actor in a movie version of the Passion of Christ, so that the film creates a complex network of relations between film shooting, social differences, art, hunger, consumption, time and light, which turns the film set into a space for displaying political relations, differences, exploitation and revolution. The correspondences between these concepts and some aggression techniques of Cinema Novo are numerous and confirm the capacity of Pasolini’s film to project ideas on cinema and politics beyond its particular production context.