A Gaussian-mixture ensemble transform filter

We generalize the popular ensemble Kalman filter to an ensemble transform filter, in which the prior distribution can take the form of a Gaussian mixture or a Gaussian kernel density estimator. The design of the filter is based on a continuous formulation of the Bayesian filter analysis step. We call the new filter algorithm the ensemble Gaussian-mixture filter (EGMF). The EGMF is implemented for three simple test problems (Brownian dynamics in one dimension, Langevin dynamics in two dimensions and the three-dimensional Lorenz-63 model). It is demonstrated that the EGMF is capable of tracking systems with non-Gaussian uni- and multimodal ensemble distributions. Copyright © 2011 Royal Meteorological Society

Ensemble transform Kalman-Bucy filters

Two recent works have adapted the Kalman–Bucy filter into an ensemble setting. In the first formulation, the ensemble of perturbations is updated by the solution of an ordinary differential equation (ODE) in pseudo-time, while the mean is updated as in the standard Kalman filter. In the second formulation, the full ensemble is updated in the analysis step as the solution of single set of ODEs in pseudo-time. Neither requires matrix inversions except for the frequently diagonal observation error covariance. We analyse the behaviour of the ODEs involved in these formulations. We demonstrate that they stiffen for large magnitudes of the ratio of background error to observational error variance, and that using the integration scheme proposed in both formulations can lead to failure. A numerical integration scheme that is both stable and is not computationally expensive is proposed. We develop transform-based alternatives for these Bucy-type approaches so that the integrations are computed in ensemble space where the variables are weights (of dimension equal to the ensemble size) rather than model variables. Finally, the performance of our ensemble transform Kalman–Bucy implementations is evaluated using three models: the 3-variable Lorenz 1963 model, the 40-variable Lorenz 1996 model, and a medium complexity atmospheric general circulation model known as SPEEDY. The results from all three models are encouraging and warrant further exploration of these assimilation techniques.

Rheology of discrete failure regimes of anisotropic sea ice

A rheological model of sea ice is presented that incorporates the orientational distribution of ice thickness in leads embedded in isotropic floe ice. Sea ice internal stress is determined by coulombic, ridging and tensile failure at orientations where corresponding failure criteria are satisfied at minimum stresses. Because sea ice traction increases in thinner leads and cohesion is finite, such failure line angles are determined by the orientational distribution of sea ice thickness relative to the imposed stresses. In contrast to the isotropic case, sea ice thickness anisotropy results in these failure lines becoming dependent on the stress magnitude. Although generally a given failure criteria type can be satisfied at many directions, only two at most are considered. The strain rate is determined by shearing along slip lines accompanied by dilatancy and closing or opening across orientations affected by ridging or tensile failure. The rheology is illustrated by a yield curve determined by combining coulombic and ridging failure for the case of two pairs of isotropically formed leads of different thicknesses rotated with regard to each other, which models two events of coulombic failure followed by dilatancy and refreezing. The yield curve consists of linear segments describing coulombic and ridging yield as failure switches from one lead to another as the stress grows. Because sliding along slip lines is accompanied by dilatancy, at typical Arctic sea ice deformation rates a one-day-long deformation event produces enough open water that these freshly formed slip lines are preferential places of ridging failure.

Using attribute importance rankings within discrete choice experiments: an application to valuing bread attributes

We present a new Bayesian econometric specification for a hypothetical Discrete Choice Experiment (DCE) incorporating respondent ranking information about attribute importance. Our results indicate that a DCE debriefing question that asks respondents to rank the importance of attributes helps to explain the resulting choices. We also examine how mode of survey delivery (online and mail) impacts model performance, finding that results are not substantively a§ected by the mode of survey delivery. We conclude that the ranking data is a complementary source of information about respondent utility functions within hypothetical DCEs

Computing Fresnel integrals via modified trapezium rules

In this paper we propose methods for computing Fresnel integrals based on truncated trapezium rule approximations to integrals on the real line, these trapezium rules modified to take into account poles of the integrand near the real axis. Our starting point is a method for computation of the error function of complex argument due to Matta and Reichel (J Math Phys 34:298–307, 1956) and Hunter and Regan (Math Comp 26:539–541, 1972). We construct approximations which we prove are exponentially convergent as a function of N , the number of quadrature points, obtaining explicit error bounds which show that accuracies of 10−15 uniformly on the real line are achieved with N=12 , this confirmed by computations. The approximations we obtain are attractive, additionally, in that they maintain small relative errors for small and large argument, are analytic on the real axis (echoing the analyticity of the Fresnel integrals), and are straightforward to implement.

Estimating correlated observation error statistics using an ensemble transform Kalman filter

For certain observing types, such as those that are remotely sensed, the observation errors are correlated and these correlations are state- and time-dependent. In this work, we develop a method for diagnosing and incorporating spatially correlated and time-dependent observation error in an ensemble data assimilation system. The method combines an ensemble transform Kalman filter with a method that uses statistical averages of background and analysis innovations to provide an estimate of the observation error covariance matrix. To evaluate the performance of the method, we perform identical twin experiments using the Lorenz ’96 and Kuramoto-Sivashinsky models. Using our approach, a good approximation to the true observation error covariance can be recovered in cases where the initial estimate of the error covariance is incorrect. Spatial observation error covariances where the length scale of the true covariance changes slowly in time can also be captured. We find that using the estimated correlated observation error in the assimilation improves the analysis.

Multiple, discrete arcs on sunward convecting field lines in the 14-15 MLT region

Ionospheric plasma flow measurements and simultaneous observations of thin (∼0.2° invariant latitude (ILAT)), multiple, longitudinally extended auroral arcs of transient nature within 74°-76° ILAT and 1030-1130 UT (∼14-15 MLT) on January 12, 1989, are reported. The auroral structures appeared within the luminous belt of strong 630.0-nm emissions located predominantly on sunward convecting field lines equatorward of the convection reversal boundary as identified by the European Incoherent Scatter UHF radar. The events occurred during a period of several hours quasi-steady solar wind speed (∼ 700 km s−1) and a radially orientated interplanetary magnetic field (IMF) with a weak northward tilt (IMF Bz>0). These typical dayside auroral features are related to previous studies of auroral activity related to the upward region 1 current in the postnoon sector. The discrete auroral events presented here may result from magnetosheath plasma injections into the low-latitude boundary layer (LLBL) and an associated dynamo mechanism. An alternative explanation invokes kinetic Alfvén waves, triggered either by Kelvin-Helmholtz instability at the inner (or outer) edge of the LLBL or by pressure pulse induced magnetopause surface waves.

Acoustic scattering : high frequency boundary element methods and unified transform methods

We describe some recent advances in the numerical solution of acoustic scattering problems. A major focus of the paper is the efficient solution of high frequency scattering problems via hybrid numerical-asymptotic boundary element methods. We also make connections to the unified transform method due to A. S. Fokas and co-authors, analysing particular instances of this method, proposed by J. A. De-Santo and co-authors, for problems of acoustic scattering by diffraction gratings.

Low complexity equalization of HCM systems with DPFFT demodulation over doubly-selective channels

To mitigate the inter-carrier interference (ICI) of doubly-selective (DS) fading channels, we consider a hybrid carrier modulation (HCM) system employing the discrete partial fast Fourier transform (DPFFT) demodulation and the banded minimum mean square error (MMSE) equalization in this letter. We first provide the discrete form of partial FFT demodulation, then apply the banded MMSE equalization to suppress the residual interference at the receiver. The proposed algorithm has been demonstrated, via numerical simulations, to be its superior over the single carrier modulation (SCM) system and circularly prefixed orthogonal frequency division multiplexing (OFDM) system over a typical DS channel. Moreover, it represents a good trade-off between computational complexity and performance.

Application of two-dimensional wavelet transform in near-shore X-band radar images

Among existing remote sensing applications, land-based X-band radar is an effective technique to monitor the wave fields, and spatial wave information could be obtained from the radar images. Two-dimensional Fourier Transform (2-D FT) is the common algorithm to derive the spectra of radar images. However, the wave field in the nearshore area is highly non-homogeneous due to wave refraction, shoaling, and other coastal mechanisms. When applied in nearshore radar images, 2-D FT would lead to ambiguity of wave characteristics in wave number domain. In this article, we introduce two-dimensional Wavelet Transform (2-D WT) to capture the non-homogeneity of wave fields from nearshore radar images. The results show that wave number spectra by 2-D WT at six parallel space locations in the given image clearly present the shoaling of nearshore waves. Wave number of the peak wave energy is increasing along the inshore direction, and dominant direction of the spectra changes from South South West (SSW) to West South West (WSW). To verify the results of 2-D WT, wave shoaling in radar images is calculated based on dispersion relation. The theoretical calculation results agree with the results of 2-D WT on the whole. The encouraging performance of 2-D WT indicates its strong capability of revealing the non-homogeneity of wave fields in nearshore X-band radar images.

Discontinuous Galerkin methods for the p-biharmonic equation from a discrete variational perspective

We study discontinuous Galerkin approximations of the p-biharmonic equation for p∈(1,∞) from a variational perspective. We propose a discrete variational formulation of the problem based on an appropriate definition of a finite element Hessian and study convergence of the method (without rates) using a semicontinuity argument. We also present numerical experiments aimed at testing the robustness of the method.

A moment problem for random discrete measures

Let X be a locally compact Polish space. A random measure on X is a probability measure on the space of all (nonnegative) Radon measures on X. Denote by K(X) the cone of all Radon measures η on X which are of the form η =

Molecular design of a discrete chain-folding polyimide for controlled inkjet deposition of supramolecular polymers

A supramolecular polymer based upon two complementary polymer components is formed by sequential deposition from solution in THF, using a piezoelectric drop-on-demand inkjet printer. Highly efficient cycloaddition or ‘click’ chemistry afforded a well-defined poly(ethylene glycol) featuring chain-folding diimide end groups, which possesses greatly enhanced solubility in THF relative to earlier materials featuring random diimide sequences. Blending the new polyimide with a complementary poly(ethylene glycol) system bearing pyrene end groups (which bind to the chain-folding diimide units) overcomes the limited solubility encountered previously with chain-folding polyimides in inkjet printing applications. The solution state properties of the resulting polymer blend were assessed via viscometry to confirm the presence of a supramolecular polymer before depositing the two electronically complementary polymers by inkjet printing techniques. The novel materials so produced offer an insight into ways of controlling the properties of printed materials through tuning the structure of the polymer at the (supra)molecular level.

Ectomycorrhizal symbiosis can enhance plant nutrition through improved access to discrete organic nutrient patches of high resource quality

It is known that roots can respond to patches of fertility; however, root proliferation is often too slow to exploit resources fully, and organic nutrient patches may be broken down and leached, immobilized or chemically fixed before they are invaded by the root system. The ability of fungal hyphae to exploit resource patches is far greater than that of roots due to their innate physiological and morphological plasticity, which allows comprehensive exploration and rapid colonization of resource patches in soils. The fungal symbionts of ectomycorrhizal plants excrete significant quantities of enzymes such as chitinases, phosphatases and proteases. These might allow the organic residue to be tapped directly for nutrients such as N and P. Pot experiments conducted with nutrient-stressed ectomycorrhizal and control willow plants showed that when high quality organic nutrient patches were added, they were colonized rapidly by the ectomycorrhizal mycelium. These established willows (0.5 m tall) were colonized by Hebeloma syrjense P. Karst. for 1 year prior to nutrient patch addition. Within days after patch addition, colour changes in the leaves of the mycorrhizal plants (reflecting improved nutrition) were apparent, and after I month the concentration of N and P in the foliage of mycorrhizal plants was significantly greater than that in non-mycorrhizal plants subject to the same nutrient addition. It seems likely that the mycorrhizal plants were able to compete effectively with the wider soil microbiota and tap directly into the high quality organic resource patch via their extra-radical mycelium. We hypothesize that ectomycorrhizal plants may reclaim some of the N and P invested in seed production by direct recycling from failed seeds in the soil. The rapid exploitation of similar discrete, transient, high-quality nutrient patches may have led to underestimations when determining the nutritional benefits of ectomycorrhizal colonization.