Fast identification algorithms for Gaussian process model

A class identification algorithms is introduced for Gaussian process(GP)models.The fundamental approach is to propose a new kernel function which leads to a covariance matrix with low rank,a property that is consequently exploited for computational efficiency for both model parameter estimation and model predictions.The objective of either maximizing the marginal likelihood or the Kullback–Leibler (K–L) divergence between the estimated output probability density function(pdf)and the true pdf has been used as respective cost functions.For each cost function,an efficient coordinate descent algorithm is proposed to estimate the kernel parameters using a one dimensional derivative free search, and noise variance using a fast gradient descent algorithm. Numerical examples are included to demonstrate the effectiveness of the new identification approaches.

Development of a food-exchange model to replace saturated fat with MUFAS and n-6 PUFAS in adults at moderate cardiovascular risk

The recommendation to reduce saturated fatty acid (SFA) consumption to ≤10% of total energy (%TE) is a key public health target aimed at lowering cardiovascular disease (CVD) risk. Replacement of SFA with unsaturated fats may provide greater benefit than replacement with carbohydrates, yet the optimal type of fat is unclear. The aim was to develop a flexible food-exchange model to investigate the effects of substituting SFAs with monounsaturated fatty acids (MUFAs) or n-6 (ω-6) polyunsaturated fatty acids (PUFAs) on CVD risk factors. In this parallel study, UK adults aged 21-60 y with moderate CVD risk (50% greater than the population mean) were identified using a risk assessment tool (n = 195; 56% females). Three 16-wk isoenergetic diets of specific fatty acid (FA) composition (%TE SFA:%TE MUFA:%TE n-6 PUFA) were designed using spreads, oils, dairy products, and snacks as follows: 1) SFA-rich diet (17:11:4; n = 65); 2) MUFA-rich diet (9:19:4; n = 64); and 3) n-6 PUFA-rich diet (9:13:10; n = 66). Each diet provided 36%TE total fat. Dietary targets were broadly met for all intervention groups, reaching 17.6 ± 0.4%TE SFA, 18.5 ± 0.3%TE MUFA, and 10.4 ± 0.3%TE n-6 PUFA in the respective diets, with significant overall diet effects for the changes in SFA, MUFA, and n-6 PUFA between groups (P < 0.001). There were no differences in the changes of total fat, protein, carbohydrate, and alcohol intake or anthropometric measures between groups. Plasma phospholipid FA composition showed changes from baseline in the proportions of total SFA, MUFA, and n-6 PUFA for each diet group, with significant overall diet effects for total SFA and MUFA between groups (P < 0.001). In conclusion, successful implementation of the food-exchange model broadly achieved the dietary target intakes for the exchange of SFA with MUFA or n-6 PUFA with minimal disruption to the overall diet in a free-living population. This trial was registered at clinicaltrials.gov as NCT01478958.

Can a climate model reproduce extreme regional precipitation events over England and Wales?

The ability of the HiGEM climate model to represent high-impact, regional, precipitation events is investigated in two ways. The first focusses on a case study of extreme regional accumulation of precipitation during the passage of a summer extra-tropical cyclone across southern England on 20 July 2007 that resulted in a national flooding emergency. The climate model is compared with a global Numerical Weather Prediction (NWP) model and higher resolution, nested limited area models. While the climate model does not simulate the timing and location of the cyclone and associated precipitation as accurately as the NWP simulations, the total accumulated precipitation in all models is similar to the rain gauge estimate across England and Wales. The regional accumulation over the event is insensitive to horizontal resolution for grid spacings ranging from 90km to 4km. Secondly, the free-running climate model reproduces the statistical distribution of daily precipitation accumulations observed in the England-Wales precipitation record. The model distribution diverges increasingly from the record for longer accumulation periods with a consistent under-representation of more intense multi-day accumulations. This may indicate a lack of low-frequency variability associated with weather regime persistence. Despite this, the overall seasonal and annual precipitation totals from the model are still comparable to those from ERA-Interim.

Curl free pressure gradients over orography in a solution of the fully compressible Euler equations with implicit treatment of acoustic and gravity waves

Steep orography can cause noisy solutions and instability in models of the atmosphere. A new technique for modelling flow over orography is introduced which guarantees curl free gradients on arbitrary grids, implying that the pressure gradient term is not a spurious source of vorticity. This mimetic property leads to better hydrostatic balance and better energy conservation on test cases using terrain following grids. Curl-free gradients are achieved by using the co-variant components of velocity over orography rather than the usual horizontal and vertical components. In addition, gravity and acoustic waves are treated implicitly without the need for mean and perturbation variables or a hydrostatic reference profile. This enables a straightforward description of the implicit treatment of gravity waves. Results are presented of a resting atmosphere over orography and the curl-free pressure gradient formulation is advantageous. Results of gravity waves over orography are insensitive to the placement of terrain-following layers. The model with implicit gravity waves is stable in strongly stratified conditions, with N∆t up to at least 10 (where N is the Brunt-V ̈ais ̈al ̈a frequency). A warm bubble rising over orography is simulated and the curl free pressure gradient formulation gives much more accurate results for this test case than a model without this mimetic property.

Plasma free fatty acids do not provide the link between obesity and insulin resistance or β-cell dysfunction: results of the Reading, Imperial, Surrey, Cambridge, Kings (RISCK) study

Aims To investigate the relationship between adiposity and plasma free fatty acid levels and the influence of total plasma free fatty acid level on insulin sensitivity and β-cell function. Methods An insulin sensitivity index, acute insulin response to glucose and a disposition index, derived from i.v. glucose tolerance minimal model analysis and total fasting plasma free fatty acid levels were available for 533 participants in the Reading, Imperial, Surrey, Cambridge, Kings study. Bivariate correlations were made between insulin sensitivity index, acute insulin response to glucose and disposition index and both adiposity measures (BMI, waist circumference and body fat mass) and total plasma free fatty acid levels. Multivariate linear regression analysis was performed, controlling for age, sex, ethnicity and adiposity. Results After adjustment, all adiposity measures were inversely associated with insulin sensitivity index (BMI: β = −0.357; waist circumference: β = −0.380; body fat mass: β = −0.375) and disposition index (BMI: β = −0.215; waist circumference: β = −0.248; body fat mass: β = −0.221) and positively associated with acute insulin response to glucose [BMI: β = 0.200; waist circumference: β = 0.195; body fat mass β = 0.209 (P values <0.001)]. Adiposity explained 13, 4 and 5% of the variation in insulin sensitivity index, acute insulin response to glucose and disposition index, respectively. After adjustment, no adiposity measure was associated with free fatty acid level, but total plasma free fatty acid level was inversely associated with insulin sensitivity index (β = −0.133), acute insulin response to glucose (β = −0.148) and disposition index [β = −0.218 (P values <0.01)]. Plasma free fatty acid concentration accounted for 1.5, 2 and 4% of the variation in insulin sensitivity index, acute insulin response to glucose and disposition index, respectively. Conclusions Plasma free fatty acid levels have a modest negative association with insulin sensitivity, β-cell secretion and disposition index but no association with adiposity measures. It is unlikely that plasma free fatty acids are the primary mediators of obesity-related insulin resistance or β-cell dysfunction.

Assessment of the anthelmintic activity of medicinal plant extracts and purified condensed tannins against free-living and parasitic stages of Oesophagostomum dentatum

Background: Plant-derived condensed tannins (CT) show promise as a complementary option to treat gastrointestinal helminth infections, thus reducing reliance on synthetic anthelmintic drugs. Most studies on the anthelmintic effects of CT have been conducted on parasites of ruminant livestock. Oesophagostomum dentatum is an economically important parasite of pigs, as well as serving as a useful laboratory model of helminth parasites due to the ability to culture it in vitro for long periods through several life-cycle stages. Here, we investigated the anthelmintic effects of CT on multiple life-cycles stages of O. dentatum. Methods: Extracts and purified fractions were prepared from five plants containing CT and analysed by HPLC-MS. Anthelmintic activity was assessed at five different stages of the O. dentatum life cycle; the development of eggs to infective third-stage larvae (L3), the parasitic L3 stage, the moult from L3 to fourth-stage larvae (L4), the L4 stage and the adult stage. Results: Free-living larvae of O. dentatum were highly susceptible to all five plant extracts. In contrast, only two of the five extracts had activity against L3, as evidenced by migration inhibition assays, whilst three of the five extracts inhibited the moulting of L3 to L4. All five extracts reduced the motility of L4, and the motility of adult worms exposed to a CT-rich extract derived from hazelnut skins was strongly inhibited, with electron microscopy demonstrating direct damage to the worm cuticle and hypodermis. Purified CT fractions retained anthelmintic activity, and depletion of CT from extracts by pre-incubation in polyvinylpolypyrrolidone removed anthelmintic effects, strongly suggesting CT as the active molecules. Conclusions: These results suggest that CT may have promise as an alternative parasite control option for O. dentatum in pigs, particularly against adult stages. Moreover, our results demonstrate a varied susceptibility of different life-cycle stages of the same parasite to CT, which may offer an insight into the anthelmintic mechanisms of these commonly found plant compounds.

A numerical model of the ionospheric signatures of time-varying magnetic reconnection: III. Quasi-instantaneous convection responses in the Cowley-Lockwood paradigm

Using a numerical implementation of the Cowley and Lockwood (1992) model of flow excitation in the magnetosphere–ionosphere (MI) system, we show that both an expanding (on a _12-min timescale) and a quasiinstantaneous response in ionospheric convection to the onset of magnetopause reconnection can be accommodated by the Cowley–Lockwood conceptual framework. This model has a key feature of time dependence, necessarily considering the history of the coupled MI system. We show that a residual flow, driven by prior magnetopause reconnection, can produce a quasi-instantaneous global ionospheric convection response; perturbations from an equilibrium state may also be present from tail reconnection, which will superpose constructively to give a similar effect. On the other hand, when the MI system is relatively free of pre-existing flow, we can most clearly see the expanding nature of the response. As the open-closed field line boundary will frequently be in motion from such prior reconnection (both at the dayside magnetopause and in the cross-tail current sheet), it is expected that there will usually be some level of combined response to dayside reconnection.

A numerical model of the ionospheric signatures of time-varying magnetic reconnection: II. Measuring expansions in the ionospheric flow response

A numerical model embodying the concepts of the Cowley-Lockwood (Cowley and Lockwood, 1992, 1997) paradigm has been used to produce a simple Cowley– Lockwood type expanding flow pattern and to calculate the resulting change in ion temperature. Cross-correlation, fixed threshold analysis and threshold relative to peak are used to determine the phase speed of the change in convection pattern, in response to a change in applied reconnection. Each of these methods fails to fully recover the expansion of the onset of the convection response that is inherent in the simulations. The results of this study indicate that any expansion of the convection pattern will be best observed in time-series data using a threshold which is a fixed fraction of the peak response. We show that these methods used to determine the expansion velocity can be used to discriminate between the two main models for the convection response to a change in reconnection.

An isotope dilution model for partitioning phenylalanine and tyrosine uptake by the mammary gland of lactating dairy cows

An isotope dilution model for partitioning phenylalanine and tyrosine uptake by the mammary gland of the lactating dairy cow is constructed and solved in the steady state. The model contains four intracellular and four extracellular pools and conservation of mass principles are applied to generate the fundamental equations describing the behaviour of the system. The experimental measurements required for model solution are milk secretion and plasma flow rate across the gland in combination with phenylalanine and tyrosine concentrations and plateau isotopic enrichments in arterial and venous plasma and free and protein bound milk during a constant infusion of [1-(13)C]phenylalanine and [2,3,5,6-(2)H]tyrosine tracer. If assumptions are made, model solution enables determination of steady state flows for phenylalanine and tyrosine inflow to the gland, outflow from it and bypass, and flows representing the synthesis and degradation of constitutive protein and hydroxylation. The model is effective in providing information about the fates of phenylalanine and tyrosine in the mammary gland and could be used as part of a more complex system describing amino acid metabolism in the whole ruminant.

Modelling signatures of pulsed magnetopause reconnection in cusp ion dispersion signatures seen at middle altitudes

It is shown that the open magnetosphere model can reproduce both the down-going and the up-going magnetosheath ions seen in the cusp and mantle regions by the Polar satellite at middle altitudes. ?he pass studied shows a series of discontinuities in the ion dispersion, most of which are shown to arise from pulses of magnetopause reconnection rate. A total of 9 pulses are detected in an interval estimated to be about 30 min long, giving a mean repetition period of about 3 min: they vary in length between 0.5 min and 3.5 min and are separated by periods of much slower reconnection of duration 1-3 min. One step is not as predicted for reconnection rate pulses but is explained in terms of compressive motions caused by a pulse of solar wind dynamic pressure. The reconnection site is found to be 16 +/- 3 R-E from the ionosphere along the separatrix field line, placing it at low latitudes on the dayside magnetopause.

Energy and pitch-angle dispersions of LLBL/cusp ions seen at middle altitudes: predictions by the open magnetosphere model

Numerical simulations are presented of the ion distribution functions seen by middle-altitude spacecraft in the low-latitude boundary layer (LLBL) and cusp regions when reconnection is, or has recently been, taking place at the equatorial magnetopause. From the evolution of the distribution function with time elapsed since the field line was opened, both the observed energy/observation-time and pitch-angle/energy dispersions are well reproduced. Distribution functions showing a mixture of magnetosheath and magnetospheric ions, often thought to be a signature of the LLBL, are found on newly opened field lines as a natural consequence of the magnetopause effects on the ions and their flight times. In addition, it is shown that the extent of the source region of the magnetosheath ions that are detected by a satellite is a function of the sensitivity of the ion instrument . If the instrument one-count level is high (and/or solar-wind densities are low), the cusp ion precipitation detected comes from a localised region of the mid-latitude magnetopause (around the magnetic cusp), even though the reconnection takes place at the equatorial magnetopause. However, if the instrument sensitivity is high enough, then ions injected from a large segment of the dayside magnetosphere (in the relevant hemisphere) will be detected in the cusp. Ion precipitation classed as LLBL is shown to arise from the low-latitude magnetopause, irrespective of the instrument sensitivity. Adoption of threshold flux definitions has the same effect as instrument sensitivity in artificially restricting the apparent source region.

Ion populations on open field lines within the dayside low-latitude boundary layer: theory and observations during a transient event

We analyze of ion populations observed by the NOAA-12 satellite within dayside auroral transients. The data are matched with an open magnetopause model which allows for the transmission of magnetosheath ions across one or both of the two Alfvén waves which emanate from the magnetopause reconnection site. It also allows for reflection and acceleration of ions of magnetospheric origin by these waves. From the good agreement found between the model and the observations, we propose that the events and the low-latitude boundary precipitation are both on open field lines.

An analysis of the accuracy of magnetopause reconnection rate variations deduced from cusp ion dispersion characteristics

We present an analysis of the accuracy of the method introduced by Lockwood et al. (1994) for the determination of the magnetopause reconnection rate from the dispersion of precipitating ions in the ionospheric cusp region. Tests are made by applying the method to synthesised data. The simulated cusp ion precipitation data are produced by an analytic model of the evolution of newly-opened field lines, along which magnetosheath ions are firstly injected across the magnetopause and then dispersed as they propagate into the ionosphere. The rate at which these newly opened field lines are generated by reconnection can be varied. The derived reconnection rate estimates are then compared with the input variation to the model and the accuracy of the method assessed. Results are presented for steady-state reconnection, for continuous reconnection showing a sine-wave variation in rate and for reconnection which only occurs in square wave pulses. It is found that the method always yields the total flux reconnected (per unit length of the open-closed field-line boundary) to within an accuracy of better than 5%, but that pulses tend to be smoothed so that the peak reconnection rate within the pulse is underestimated and the pulse length is overestimated. This smoothing is reduced if the separation between energy channels of the instrument is reduced; however this also acts to increase the experimental uncertainty in the estimates, an effect which can be countered by improving the time resolution of the observations. The limited time resolution of the data is shown to set a minimum reconnection rate below which the method gives spurious short-period oscillations about the true value. Various examples of reconnection rate variations derived from cusp observations are discussed in the light of this analysis.

Overlapping cusp ion injections: An explanation invoking magnetopause reconnection

An explanation of overlapping cusp ion injections is presented using the pulsating cusp model of the effects of magnetopause reconnection. It is shown that two populations of cusp ions, covering separated energy ranges, can be seen simultaneously by low- or mid-altitude satellites because of the combined effect of the acceleration and the straightening of newly-opened field lines as they evolve away from the reconnection site. Observations of such signatures, recently reported in data from the Viking and Freja satellites, are discussed in terms of pulsed and steady reconnection.

Measuring ion temperatures and studying the ion energy balance in the high-latitude ionosphere

Data are presented for a nighttime ion heating event observed by the EISCAT radar on 16 December 1988. In the experiment, the aspect angle between the radar beam and the geomagnetic field was fixed at 54.7°, which avoids any ambiguity in derived ion temperature caused by anisotropy in the ion velocity distribution function. The data were analyzed with an algorithm which takes account of the non-Maxwellian line-of-sight ion velocity distribution. During the heating event, the derived spectral distortion parameter (D∗) indicated that the distribution function was highly distorted from a Maxwellian form when the ion drift increased to 4 km s−1. The true three-dimensional ion temperature was used in the simplified ion balance equation to compute the ion mass during the heating event. The ion composition was found to change from predominantly O4 to mainly molecular ions. A theoretical analysis of the ion composition, using the MSIS86 model and published values of the chemical rate coefficients, accounts for the order-of-magnitude increase in the atomic/molecular ion ratio during the event, but does not successfully explain the very high proportion of molecular ions that was observed.