Gaussian quadrature 4D-VAR

A new incremental four-dimensional variational (4D-Var) data assimilation algorithm is introduced. The algorithm does not require the computationally expensive integrations with the nonlinear model in the outer loops. Nonlinearity is accounted for by modifying the linearization trajectory of the observation operator based on integrations with the tangent linear (TL) model. This allows us to update the linearization trajectory of the observation operator in the inner loops at negligible computational cost. As a result the distinction between inner and outer loops is no longer necessary. The key idea on which the proposed 4D-Var method is based is that by using Gaussian quadrature it is possible to get an exact correspondence between the nonlinear time evolution of perturbations and the time evolution in the TL model. It is shown that J-point Gaussian quadrature can be used to derive the exact adjoint-based observation impact equations and furthermore that it is straightforward to account for the effect of multiple outer loops in these equations if the proposed 4D-Var method is used. The method is illustrated using a three-level quasi-geostrophic model and the Lorenz (1996) model.

Measures of observation impact in non-Gaussian data assimilation

ABSTRACT Non-Gaussian/non-linear data assimilation is becoming an increasingly important area of research in the Geosciences as the resolution and non-linearity of models are increased and more and more non-linear observation operators are being used. In this study, we look at the effect of relaxing the assumption of a Gaussian prior on the impact of observations within the data assimilation system. Three different measures of observation impact are studied: the sensitivity of the posterior mean to the observations, mutual information and relative entropy. The sensitivity of the posterior mean is derived analytically when the prior is modelled by a simplified Gaussian mixture and the observation errors are Gaussian. It is found that the sensitivity is a strong function of the value of the observation and proportional to the posterior variance. Similarly, relative entropy is found to be a strong function of the value of the observation. However, the errors in estimating these two measures using a Gaussian approximation to the prior can differ significantly. This hampers conclusions about the effect of the non-Gaussian prior on observation impact. Mutual information does not depend on the value of the observation and is seen to be close to its Gaussian approximation. These findings are illustrated with the particle filter applied to the Lorenz ’63 system. This article is concluded with a discussion of the appropriateness of these measures of observation impact for different situations.

Optimal linearization trajectories for tangent linear models

We examine differential equations where nonlinearity is a result of the advection part of the total derivative or the use of quadratic algebraic constraints between state variables (such as the ideal gas law). We show that these types of nonlinearity can be accounted for in the tangent linear model by a suitable choice of the linearization trajectory. Using this optimal linearization trajectory, we show that the tangent linear model can be used to reproduce the exact nonlinear error growth of perturbations for more than 200 days in a quasi-geostrophic model and more than (the equivalent of) 150 days in the Lorenz 96 model. We introduce an iterative method, purely based on tangent linear integrations, that converges to this optimal linearization trajectory. The main conclusion from this article is that this iterative method can be used to account for nonlinearity in estimation problems without using the nonlinear model. We demonstrate this by performing forecast sensitivity experiments in the Lorenz 96 model and show that we are able to estimate analysis increments that improve the two-day forecast using only four backward integrations with the tangent linear model. Copyright © 2011 Royal Meteorological Society

On the approximation of local and linear radiative damping in the middle atmosphere

The validity of approximating radiative heating rates in the middle atmosphere by a local linear relaxation to a reference temperature state (i.e., ‘‘Newtonian cooling’’) is investigated. Using radiative heating rate and temperature output from a chemistry–climate model with realistic spatiotemporal variability and realistic chemical and radiative parameterizations, it is found that a linear regressionmodel can capture more than 80% of the variance in longwave heating rates throughout most of the stratosphere and mesosphere, provided that the damping rate is allowed to vary with height, latitude, and season. The linear model describes departures from the climatological mean, not from radiative equilibrium. Photochemical damping rates in the upper stratosphere are similarly diagnosed. Threeimportant exceptions, however, are found.The approximation of linearity breaks down near the edges of the polar vortices in both hemispheres. This nonlinearity can be well captured by including a quadratic term. The use of a scale-independentdamping rate is not well justified in the lower tropical stratosphere because of the presence of a broad spectrum of vertical scales. The local assumption fails entirely during the breakup of the Antarctic vortex, where large fluctuations in temperature near the top of the vortex influence longwave heating rates within the quiescent region below. These results are relevant for mechanistic modeling studies of the middle atmosphere, particularly those investigating the final Antarctic warming.

Microcontroller-based peak current mode control using digital slope compensation

Microcontroller-based peak current mode control of a buck converter is investigated. The new solution uses a discrete time controller with digital slope compensation. This is implemented using only a single-chip microcontroller to achieve desirable cycle-by-cycle peak current limiting. The digital controller is implemented as a two-pole, two-zero linear difference equation designed using a continuous time model of the buck converter and a discrete time transform. Subharmonic oscillations are removed with digital slope compensation using a discrete staircase ramp. A 16 W hardware implementation directly compares analog and digital control. Frequency response measurements are taken and it is shown that the crossover frequency and expected phase margin of the digital control system match that of its analog counterpart.

Classification and fault detection methods for fuel cell monitoring and quality control

In this paper, various types of fault detection methods for fuel cells are compared. For example, those that use a model based approach or a data driven approach or a combination of the two. The potential advantages and drawbacks of each method are discussed and comparisons between methods are made. In particular, classification algorithms are investigated, which separate a data set into classes or clusters based on some prior knowledge or measure of similarity. In particular, the application of classification methods to vectors of reconstructed currents by magnetic tomography or to vectors of magnetic field measurements directly is explored. Bases are simulated using the finite integration technique (FIT) and regularization techniques are employed to overcome ill-posedness. Fisher's linear discriminant is used to illustrate these concepts. Numerical experiments show that the ill-posedness of the magnetic tomography problem is a part of the classification problem on magnetic field measurements as well. This is independent of the particular working mode of the cell but influenced by the type of faulty behavior that is studied. The numerical results demonstrate the ill-posedness by the exponential decay behavior of the singular values for three examples of fault classes.

The effect of adding stinging nettle (Urtica dioica) haylage to a total mixed ration on performance and rumen function of lactating dairy cows

In vitro studies found that inclusion of dried stinging nettle (Urtica dioica) at 100 mg/g dry matter (DM) increased the pH of a rumen fluid inoculated fermentation buffer by 30% and the effect was persistent for 7 days. Our objective was to evaluate the effects of adding stinging nettle haylage to a total mixed ration on feed intake, eating and rumination activity, rumen pH, milk yield, and milk composition of lactating dairy cows. Six lactating Holstein-Friesian cows were used in a replicated 3 × 3 Latin Square design experiment with 3 treatments and 3 week periods. Treatments were a control (C) high-starch (311 g/kg DM) total mixed ration diet and two treatment diets containing 50 (N5) and 100 (N10) g nettle haylage (DM/kg) as a replacement for ryegrass silage (Lolium perenne). There was an increase (linear, P < 0.010) in the proportion of large particles and a reduction in medium (linear, P = 0.045) and fine particles (linear, P = 0.026) in the diet offered with increasing nettle inclusion. A numerical decrease (linear, P = 0.106) in DM intake (DMI) was observed as nettle inclusion in the diet increased. Milk yield averaged 20.3 kg/day and was not affected by diet. There was a decrease (quadratic, P = 0.01) in the time animals spent ruminating as nettle inclusion in the diet increased, in spite of an increase in the number of boli produced daily for the N5 diet (quadratic, P = 0.031). Animals fed the N10 diet spent less time with a rumen pH below 5.5 (P < 0.05) than cows fed the N5 diet. Averaged over an 8.5 h sampling period, there were no changes in the concentration or proportions of acetate or propionate in the rumen, but feeding nettle haylage reduced the concentrations of n-butyrate (quadratic, P < 0.001), i-butyrate (linear, P < 0.009) and n-caproate (linear, P < 0.003). Milk and fat and protein corrected milk yield were not affected when nettles replaced ryegrass silage in the diet of lactating dairy cows, despite a numerical reduction in feed intake. Rumination activity was reduced by the addition of nettle haylage to the diet, which may reflect differences in fibre structure between the nettle haylage and ryegrass silage fed. Changes observed in rumen pH suggest potential benefits of feeding nettle haylage for reducing rumen acidosis. However, the extent to which these effects were due to the fermentability and structure of the nettle haylage compared to the ryegrass silage fed, or a bioactive component of the nettles, is not certain

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.

Thermospheric Control of the Auroral Source of O+Ions for the Magnetosphere

Linear theory, model ion-density profiles and MSIS neutral thermospheric predictions are used to investigate the stability of the auroral, topside ionosphere to oxygen cyclotron waves: variations of the critical height, above which the plasma is unstable, with field-aligned current, thermal ion density and exospheric temperature are considered. In addition, probabilities are assessed that interactions with neutral atomic gases prevent O+ ions from escaping into the magnetosphere after they have been transversely accelerated by these waves. The two studies are combined to give a rough estimate of the total O+ escape flux as a function of the field-aligned current density for an assumed rise in the perpendicular ion temperature. Charge exchange with neutral oxygen, not hydrogen, is shown to be the principle limitation to the escape of O+ ions, which occurs when the waves are driven unstable down to low altitudes. It is found that the largest observed field-aligned current densities can heat a maximum of about 5×1014 O+ ions m−2 to a threshold above which they are subsequently able to escape into the magnetosphere in the following 500s. Averaged over this period, this would constitute a flux of 1012 m−2 s−1 and in steady-state the peak outflow would then be limited to about 1013 m−2 s−1 by frictional drag on thermal O+ at lower altitudes. Maximum escape is at low plasma density unless the O+ scale height is very large. The outflow decreases with decreasing field-aligned current density and, to a lesser extent, with increasing exospheric temperature. Upward flowing ion events are evaluated as a source of O+ ions for the magnetosphere and as an explanation of the observed solar cycle variation of ring current O+ abundance.

Organic soils promote the efficacy of entomopathogenic nematodes, with different foraging strategies, in the control of a major forest pest: a meta-analysis of field trial data

The large pine weevil, Hylobius abietis, is a serious pest of reforestation in northern Europe. However, weevils developing in stumps of felled trees can be killed by entomopathogenic nematodes applied to soil around the stumps and this method of control has been used at an operational level in the UK and Ireland. We investigated the factors affecting the efficacy of entomopathogenic nematodes in the control of the large pine weevil spanning 10 years of field experiments, by means of a meta-analysis of published studies and previously unpublished data. We investigated two species with different foraging strategies, the ‘ambusher’ Steinernema carpocapsae, the species most often used at an operational level, and the ‘cruiser’ Heterorhabditis downesi. Efficacy was measured both by percentage reduction in numbers of adults emerging relative to untreated controls and by percentage parasitism of developing weevils in the stump. Both measures were significantly higher with H. downesi compared to S. carpocapsae. General linear models were constructed for each nematode species separately, using substrate type (peat versus mineral soil) and tree species (pine versus spruce) as fixed factors, weevil abundance (from the mean of untreated stumps) as a covariate and percentage reduction or percentage parasitism as the response variable. For both nematode species, the most significant and parsimonious models showed that substrate type was consistently, but not always, the most significant variable, whether replicates were at a site or stump level, and that peaty soils significantly promote the efficacy of both species. Efficacy, in terms of percentage parasitism, was not density dependent.

Estimation of Gaussian process regression model using probability distance measures

A new class of parameter estimation algorithms is introduced for Gaussian process regression (GPR) models. It is shown that the integration of the GPR model with probability distance measures of (i) the integrated square error and (ii) Kullback–Leibler (K–L) divergence are analytically tractable. An efficient coordinate descent algorithm is proposed to iteratively estimate the kernel width using golden section search which includes a fast gradient descent algorithm as an inner loop to estimate the noise variance. Numerical examples are included to demonstrate the effectiveness of the new identification approaches.

Impact of variation in structure of condensed tannins from sainfoin (Onobrychis viciifolia) on in vitro ruminal methan production and fermentation characteristics

Our study investigated the effects of condensed tannins (CT) on rumen in vitro methane (CH4) production and fermentation characteristics by incubating lucerne in buffered rumen fluid in combination with different CT extracts at 0 (control), 40, 80 and 120 g CT/kg of substrate DM. Condensed tannins were extracted from four sainfoin accessions: Rees ‘A’, CPI63763, Cotswold Common and CPI63767. Gas production (GP) was measured using a fully automated GP apparatus with CH4 measured at distinct time points. Condensed tannins differed substantially in terms of polymer size and varied from 13 (Rees ‘A’) to 73 (CPI63767) mean degree of polymerization, but had relatively similar characteristics in terms of CT content, procyanidin: prodelphinidin (PC: PD) and cis:trans ratios. Compared to control, addition of CT from CPI63767 and CPI63763 at 80 and 120 g CT/kg of substrate DM reduced CH4 by 43% and 65%, and by 23% and 57%, respectively, after 24-h incubation. Similarly, CT from Rees ‘A’ and Cotswold Common reduced CH4 by 26% and 46%, and by 28% and 46% respectively. Addition of increasing level of CT linearly reduced the maximum rates of GP and CH4 production, and the estimated in vitro organic matter digestibility. There was a negative linear and quadratic (p < 0.01) relation between CT concentration and total volatile fatty acid (VFA) production. Inclusion of 80 and 120 g CT/kg of substrate DM reduced (p < 0.001) branched-chain VFA production and acetate: propionate ratio and was lowest for CPI63767. A decrease in proteolytic activity as indirectly shown by a change in VFA composition favouring a shift towards propionate and reduction in branched-chain VFA production varied with type of CT and was highest for CPI63767. In conclusion, these results suggest that tannin polymer size is an important factor affecting in vitro CH4 production which may be linked to the CT interaction with dietary substrate or microbial cells.

Optimising mirid control on cocoa farms through complementary monitoring systems

Mirids (Sahlbergella singularis and Distantiella theobroma) are the most important insect pests affecting cocoa production across West Africa. Understanding the population dynamics of mirids is key to their management, however, the current recommended hand-height assessment method is labour intensive. The objective of the study was to compare recently developed mirid sex pheromone trapping and visual hand-height assessment methods as monitoring tools on cocoa farms and to consider implications for a decision support system. Ten farms from the Eastern and Ashanti regions of Ghana were used for the study. Mirid numbers and damage were assessed fortnightly on twenty trees per farm, using both methods, from January 2012 to April 2013. The mirid population increased rapidly in June, reached a peak in September and began to decline in October. There was a significant linear relationship between numbers of mirids sampled to hand-height and mirid damage. High numbers of male mirids were recorded in pheromone traps between January and April 2012 after which there was a gradual decline. There was a significant inverse relationship between numbers of trapped adult mirids and mirids sampled to hand-height (predominantly nymphs). Higher temperatures and lower relative humidities in the first half of the year were associated with fewer mirids at hand-height but larger numbers of adult males were caught in pheromone traps. The study showed that relying solely on one method is not sufficient to provide accurate information on mirid population dynamics and a combination of the two methods is necessary.

Exact Search Directions for Optimization of Linear and Nonlinear Models Based on Generalized Orthonormal Functions

A novel technique for selecting the poles of orthonormal basis functions (OBF) in Volterra models of any order is presented. It is well-known that the usual large number of parameters required to describe the Volterra kernels can be significantly reduced by representing each kernel using an appropriate basis of orthonormal functions. Such a representation results in the so-called OBF Volterra model, which has a Wiener structure consisting of a linear dynamic generated by the orthonormal basis followed by a nonlinear static mapping given by the Volterra polynomial series. Aiming at optimizing the poles that fully parameterize the orthonormal bases, the exact gradients of the outputs of the orthonormal filters with respect to their poles are computed analytically by using a back-propagation-through-time technique. The expressions relative to the Kautz basis and to generalized orthonormal bases of functions (GOBF) are addressed; the ones related to the Laguerre basis follow straightforwardly as a particular case. The main innovation here is that the dynamic nature of the OBF filters is fully considered in the gradient computations. These gradients provide exact search directions for optimizing the poles of a given orthonormal basis. Such search directions can, in turn, be used as part of an optimization procedure to locate the minimum of a cost-function that takes into account the error of estimation of the system output. The Levenberg-Marquardt algorithm is adopted here as the optimization procedure. Unlike previous related work, the proposed approach relies solely on input-output data measured from the system to be modeled, i.e., no information about the Volterra kernels is required. Examples are presented to illustrate the application of this approach to the modeling of dynamic systems, including a real magnetic levitation system with nonlinear oscillatory behavior.