Multivariate analysis of the dielectric response of materials modeled using networks of resistors and capacitors

We discuss the modeling of dielectric responses of electromagnetically excited networks which are composed of a mixture of capacitors and resistors. Such networks can be employed as lumped-parameter circuits to model the response of composite materials containing conductive and insulating grains. The dynamics of the excited network systems are studied using a state space model derived from a randomized incidence matrix. Time and frequency domain responses from synthetic data sets generated from state space models are analyzed for the purpose of estimating the fraction of capacitors in the network. Good results were obtained by using either the time-domain response to a pulse excitation or impedance data at selected frequencies. A chemometric framework based on a Successive Projections Algorithm (SPA) enables the construction of multiple linear regression (MLR) models which can efficiently determine the ratio of conductive to insulating components in composite material samples. The proposed method avoids restrictions commonly associated with Archie’s law, the application of percolation theory or Kohlrausch-Williams-Watts models and is applicable to experimental results generated by either time domain transient spectrometers or continuous-wave instruments. Furthermore, it is quite generic and applicable to tomography, acoustics as well as other spectroscopies such as nuclear magnetic resonance, electron paramagnetic resonance and, therefore, should be of general interest across the dielectrics community.

A fast two-grid and finite section method for a class of integral equations on the real line with application to an acoustic scattering problem in the half-plane

We consider the numerical treatment of second kind integral equations on the real line of the form ∅(s) = ∫_(-∞)^(+∞)▒〖κ(s-t)z(t)ϕ(t)dt,s=R〗 (abbreviated ϕ= ψ+K_z ϕ) in which K ϵ L_1 (R), z ϵ L_∞ (R) and ψ ϵ BC(R), the space of bounded continuous functions on R, are assumed known and ϕ ϵ BC(R) is to be determined. We first derive sharp error estimates for the finite section approximation (reducing the range of integration to [-A, A]) via bounds on (1-K_z )^(-1)as an operator on spaces of weighted continuous functions. Numerical solution by a simple discrete collocation method on a uniform grid on R is then analysed: in the case when z is compactly supported this leads to a coefficient matrix which allows a rapid matrix-vector multiply via the FFT. To utilise this possibility we propose a modified two-grid iteration, a feature of which is that the coarse grid matrix is approximated by a banded matrix, and analyse convergence and computational cost. In cases where z is not compactly supported a combined finite section and two-grid algorithm can be applied and we extend the analysis to this case. As an application we consider acoustic scattering in the half-plane with a Robin or impedance boundary condition which we formulate as a boundary integral equation of the class studied. Our final result is that if z (related to the boundary impedance in the application) takes values in an appropriate compact subset Q of the complex plane, then the difference between ϕ(s)and its finite section approximation computed numerically using the iterative scheme proposed is ≤C_1 [kh log⁡〖(1⁄kh)+(1-Θ)^((-1)⁄2) (kA)^((-1)⁄2) 〗 ] in the interval [-ΘA,ΘA](Θ<1) for kh sufficiently small, where k is the wavenumber and h the grid spacing. Moreover this numerical approximation can be computed in ≤C_2 N log⁡N operations, where N = 2A/h is the number of degrees of freedom. The values of the constants C1 and C2 depend only on the set Q and not on the wavenumber k or the support of z.

An exploration of the equivalent weights particle filter

Particle filters are fully non-linear data assimilation techniques that aim to represent the probability distribution of the model state given the observations (the posterior) by a number of particles. In high-dimensional geophysical applications the number of particles required by the sequential importance resampling (SIR) particle filter in order to capture the high probability region of the posterior, is too large to make them usable. However particle filters can be formulated using proposal densities, which gives greater freedom in how particles are sampled and allows for a much smaller number of particles. Here a particle filter is presented which uses the proposal density to ensure that all particles end up in the high probability region of the posterior probability density function. This gives rise to the possibility of non-linear data assimilation in large dimensional systems. The particle filter formulation is compared to the optimal proposal density particle filter and the implicit particle filter, both of which also utilise a proposal density. We show that when observations are available every time step, both schemes will be degenerate when the number of independent observations is large, unlike the new scheme. The sensitivity of the new scheme to its parameter values is explored theoretically and demonstrated using the Lorenz (1963) model.

Numerical analyses of Runge–Kutta implicit–explicit schemes for horizontally explicit, vertically implicit solutions of atmospheric models

With the prospect of exascale computing, computational methods requiring only local data become especially attractive. Consequently, the typical domain decomposition of atmospheric models means horizontally-explicit vertically-implicit (HEVI) time-stepping schemes warrant further attention. In this analysis, Runge-Kutta implicit-explicit schemes from the literature are analysed for their stability and accuracy using a von Neumann stability analysis of two linear systems. Attention is paid to the numerical phase to indicate the behaviour of phase and group velocities. Where the analysis is tractable, analytically derived expressions are considered. For more complicated cases, amplification factors have been numerically generated and the associated amplitudes and phase diagnosed. Analysis of a system describing acoustic waves has necessitated attributing the three resultant eigenvalues to the three physical modes of the system. To do so, a series of algorithms has been devised to track the eigenvalues across the frequency space. The result enables analysis of whether the schemes exactly preserve the non-divergent mode; and whether there is evidence of spurious reversal in the direction of group velocities or asymmetry in the damping for the pair of acoustic modes. Frequency ranges that span next-generation high-resolution weather models to coarse-resolution climate models are considered; and a comparison is made of errors accumulated from multiple stability-constrained shorter time-steps from the HEVI scheme with a single integration from a fully implicit scheme over the same time interval. Two schemes, “Trap2(2,3,2)” and “UJ3(1,3,2)”, both already used in atmospheric models, are identified as offering consistently good stability and representation of phase across all the analyses. Furthermore, according to a simple measure of computational cost, “Trap2(2,3,2)” is the least expensive.

Representation errors and retrievals in linear and nonlinear data assimilation

This article shows how one can formulate the representation problem starting from Bayes’ theorem. The purpose of this article is to raise awareness of the formal solutions,so that approximations can be placed in a proper context. The representation errors appear in the likelihood, and the different possibilities for the representation of reality in model and observations are discussed, including nonlinear representation probability density functions. Specifically, the assumptions needed in the usual procedure to add a representation error covariance to the error covariance of the observations are discussed,and it is shown that, when several sub-grid observations are present, their mean still has a representation error ; socalled ‘superobbing’ does not resolve the issue. Connection is made to the off-line or on-line retrieval problem, providing a new simple proof of the equivalence of assimilating linear retrievals and original observations. Furthermore, it is shown how nonlinear retrievals can be assimilated without loss of information. Finally we discuss how errors in the observation operator model can be treated consistently in the Bayesian framework, connecting to previous work in this area.

Uptake of agri-environmental schemes in the Less-Favoured Areas of Greece: The role of corruption and farmers' responses to the financial crisis

Vine-growing in the Less-Favoured Areas of Greece is facing multiple challenges that might lead to its abandonment. In an attempt to maintain rural populations, Rural Development Schemes have been created that offer the opportunity to rural households to maintain or expand their farming businesses including vine-growing. This paper stems from a study that used data from a cross-sectional survey of 204 farmers to investigate how farming systems and farmers’ perception of corruption, amongst other socio-economic factors, affected their decisions to continue vine-growing through participation in Rural Development Schemes, in three remote Less-Favoured Areas of Greece. The Theory of Planned Behaviour was used to frame the research problem with the assumption being that an individual’s intention to participate in a Scheme is based on their prior beliefs about it. Data from the survey were reduced and simplified by the use of non-linear principal component analysis. The ensuing variables were used in selectivity corrected ordered probit models to reveal farmers’ attitudes towards viticulture and rural development. It was found that economic factors, perceived corruption and farmers’ attitudes were significant determinants on whether to participate in the Schemes. The research findings highlight the important role of perceived corruption and the need for policies that facilitate farmers’ access to decision making centres.

Potential of novel dextran oligosaccharides as prebiotics for obesity management through in vitro experimentation.

The energy-salvaging capacity of the gut microbiota from dietary ingredients has been proposed as a contributing factor for the development of obesity. This knowledge generated interest in the use of non-digestible dietary ingredients such as prebiotics to manipulate host energy homeostasis. In the present study, the in vitro response of obese human faecal microbiota to novel oligosaccharides was investigated. Dextrans of various molecular weights and degrees of branching were fermented with the faecal microbiota of healthy obese adults in pH-controlled batch cultures. Changes in bacterial populations were monitored using fluorescent in situ hybridisation and SCFA concentrations were analysed by HPLC. The rate of gas production and total volume of gas produced were also determined. In general, the novel dextrans and inulin increased the counts of bifidobacteria. Some of the dextrans were able to alter the composition of the obese human microbiota by increasing the counts of Bacteroides–Prevotella and decreasing those of Faecalibacterium prausnitzii and Ruminococcus bromii/R. flavefaciens. Considerable increases in SCFA concentrations were observed in response to all substrates. Gas production rates were similar during the fermentation of all dextrans, but significantly lower than those during the fermentation of inulin. Lower total gas production and shorter time to attain maximal gas production were observed during the fermentation of the linear 1 kDa dextran than during the fermentation of the other dextrans. The efficacy of bifidobacteria to ferment dextrans relied on the molecular weight and not on the degree of branching. In conclusion, there are no differences in the profiles between the obese and lean human faecal fermentations of dextrans.

The transaction cost economics theory of the family firm: family-based human asset specificity and the bifurcation bias

We develop a transaction cost economics theory of the family firm, building upon the concepts of family-based asset specificity, bounded rationality, and bounded reliability. We argue that the prosperity and survival of family firms depend on the absence of a dysfunctional bifurcation bias. The bifurcation bias is an expression of bounded reliability, reflected in the de facto asymmetric treatment of family vs. nonfamily assets (especially human assets). We propose that absence of bifurcation bias is critical to fostering reliability in family business functioning. Our study ends the unproductive divide between the agency and stewardship perspectives of the family firm, which offer conflicting accounts of this firm type's functioning. We show that the predictions of the agency and stewardship perspectives can be usefully reconciled when focusing on how family firms address the bifurcation bias or fail to do so.

A new calibrated sunspot group series since 1749: statistics of active day fractions

Although the sunspot-number series have existed since the mid-19th century, they are still the subject of intense debate, with the largest uncertainty being related to the "calibration" of the visual acuity of individual observers in the past. Daisy-chain regression methods are applied to inter-calibrate the observers which may lead to significant bias and error accumulation. Here we present a novel method to calibrate the visual acuity of the key observers to the reference data set of Royal Greenwich Observatory sunspot groups for the period 1900-1976, using the statistics of the active-day fraction. For each observer we independently evaluate their observational thresholds [S_S] defined such that the observer is assumed to miss all of the groups with an area smaller than S_S and report all the groups larger than S_S. Next, using a Monte-Carlo method we construct, from the reference data set, a correction matrix for each observer. The correction matrices are significantly non-linear and cannot be approximated by a linear regression or proportionality. We emphasize that corrections based on a linear proportionality between annually averaged data lead to serious biases and distortions of the data. The correction matrices are applied to the original sunspot group records for each day, and finally the composite corrected series is produced for the period since 1748. The corrected series displays secular minima around 1800 (Dalton minimum) and 1900 (Gleissberg minimum), as well as the Modern grand maximum of activity in the second half of the 20th century. The uniqueness of the grand maximum is confirmed for the last 250 years. It is shown that the adoption of a linear relationship between the data of Wolf and Wolfer results in grossly inflated group numbers in the 18th and 19th centuries in some reconstructions.

Rossby wave propagation on potential vorticity fronts with finite width

The horizontal gradient of potential vorticity (PV) across the tropopause typically declines with lead time in global numerical weather forecasts and tends towards a steady value dependent on model resolution. This paper examines how spreading the tropopause PV contrast over a broader frontal zone affects the propagation of Rossby waves. The approach taken is to analyse Rossby waves on a PV front of finite width in a simple single-layer model. The dispersion relation for linear Rossby waves on a PV front of infinitesimal width is well known; here an approximate correction is derived for the case of a finite width front, valid in the limit that the front is narrow compared to the zonal wavelength. Broadening the front causes a decrease in both the jet speed and the ability of waves to propagate upstream. The contribution of these changes to Rossby wave phase speeds cancel at leading order. At second order the decrease in jet speed dominates, meaning phase speeds are slower on broader PV fronts. This asymptotic phase speed result is shown to hold for a wide class of single-layer dynamics with a varying range of PV inversion operators. The phase speed dependence on frontal width is verified by numerical simulations and also shown to be robust at finite wave amplitude, and estimates are made for the error in Rossby wave propagation speeds due to the PV gradient error present in numerical weather forecast models.

Urinary metabolomic profiling to identify biomarkers of a flavonoid-rich and flavonoid-poor fruits and vegetables diet: the FLAVURS trial in adults: the FLAVURS trial

The present study aims to investigate the dose dependent effects of consuming diets enriched in flavonoid-rich and flavonoid-poor fruits and vegetables on the urine metabolome of adults who had a C1.5 fold increased risk of cardiovascular diseases. A single-blind, dose-dependent, parallel randomized controlled dietary intervention was conducted where volunteers (n = 126) were randomly assigned to one of three diets: high flavonoid diet, low flavonoid diet or habitual diet as a control for 18 weeks. High resolution LC– MS untargeted metabolomics with minimal sample cleanup was performed using an Orbitrap mass spectrometer. Putative biomarkers which characterize diets with high and low flavonoid content were selected by state-of-the-art data analysis strategies and identified by HR-MS and HR-MS/MS assays. Discrimination between diets was observed by application of two linear mixedmodels: one including a diet-time interaction effect and the second containing only a time effect. Valerolactones, phenolic acids and their derivatives were among sixteen biomarkers related to the high flavonoid dietary exposure. Four biomarkers related to the low flavonoid diet belonged to the family of phenolic acids. For the first time abscisic acid glucuronide was reported as a biomarker after a dietary intake, however its origins have to be examined by future hypothesis driven experiments using a more targeted approach. This metabolomic analysis has identified a number of dose dependent urinary biomarkers (i.e. proline betaine or iberin-N-acetyl cysteine), which can be used in future observation and intervention studies to assess flavonoids and nonflavonoid phenolic intakes and compliance to fruit and vegetable intervention.

Effect of inbreeding on growth and reproductive traits of Nellore cattle in Brazil

The objective of this study was to identify and quantify the influence of F (inbreeding coefficient) on weaning weight (WW), weight gain from weaning to 18 months of age (WG345), finishing visual score (precocity) at 18 months of age, muscling visual score at 18 months of age (MUS), hip height (HH), scrotal circumference at 18 months of age (SC), heifer probability of pregnancy at 14 months of age (PP14), and stayability (STAY) in Brazilian Nellore cattle. The complete pedigree included 417,552 animals born between 1984 and 2007 on 12 farms located in the states of Mato Grosso do Sul, Sao Paulo and Bahia. Following the observation of a statistically significant effect (P<0.05) of the covariates individual inbreeding coefficient (F) and maternal inbreeding coefficient, regression analysis of each trait, adjusted for all other effects, was performed as a function of the linear and quadratic effect of F and maternal F (when significant). Inbreeding negatively affected all traits studied (P<0.05), except for muscling. A quadratic effect of individual F on WW, WG345, HH and PP14, and a quadratic effect of maternal F on WG345 and HH were observed. Levels of inbreeding higher than 7-11% affected negatively growth and reproductive performance of Nellore cattle. Therefore, inbreeding should be avoided, except for purposes of genetic breeding whose main objective is the fixation of certain alleles in the population. (C) 2010 Elsevier B.V. All rights reserved.

Exposure to ambient levels of particles emitted by traffic worsens emphysema in mice

Objectives: We investigated effects of chronic exposure (2 months) to ambient levels of particulate matter (PM) on development of protease-induced emphysema and pulmonary remodeling in mice. Methods: Balb/c mice received nasal drop of either papain or normal saline and were kept in two exposure chambers situated in an area with high traffic density. One of them received ambient air and the other had filters for PM. Results: mean concentration of PM10 was 2.68 +/- 0.38 and 33.86 +/- 2.09 mu g/m(3), respectively, in the filtered and ambient air chambers (p<0.001). After 2 months of exposure, lungs from papain-treated mice kept in the chamber with ambient air presented greater values of mean linear intercept, an increase in density of collagen fibers in alveolar septa and in expression of 8-isoprostane (p = 0.002, p < 0.05 and p = 0.002, respectively, compared to papain-treated mice kept in the chamber with filtered air). We did not observe significant differences between these two groups in density of macrophages and in amount of cells expressing matrix metalloproteinase-12. There were no significant differences in saline-treated mice kept in the two chambers. Conclusions: We conclude that exposure to urban levels of PM worsens protease-induced emphysema and increases pulmonary remodeling. We suggest that an increase in oxidative stress induced by PM exposure influences this response. These pulmonary effects of PM were observed only in mice with emphysema. (C) 2009 Elsevier Inc. All rights reserved.

ON INTEGRABLE CODIMENSION ONE ANOSOV ACTIONS OF R(k)

In this paper, we consider codimension one Anosov actions of R(k), k >= 1, on closed connected orientable manifolds of dimension n vertical bar k with n >= 3. We show that the fundamental group of the ambient manifold is solvable if and only if the weak foliation of codimension one is transversely affine. We also study the situation where one 1-parameter subgroup of R(k) admits a cross-section, and compare this to the case where the whole action is transverse to a fibration over a manifold of dimension n. As a byproduct, generalizing a Theorem by Ghys in the case k = 1, we show that, under some assumptions about the smoothness of the sub-bundle E(ss) circle plus E(uu), and in the case where the action preserves the volume, it is topologically equivalent to a suspension of a linear Anosov action of Z(k) on T(n).

An integer linear programming approach for approximate string comparison

We introduce a problem called maximum common characters in blocks (MCCB), which arises in applications of approximate string comparison, particularly in the unification of possibly erroneous textual data coming from different sources. We show that this problem is NP-complete, but can nevertheless be solved satisfactorily using integer linear programming for instances of practical interest. Two integer linear formulations are proposed and compared in terms of their linear relaxations. We also compare the results of the approximate matching with other known measures such as the Levenshtein (edit) distance. (C) 2008 Elsevier B.V. All rights reserved.