Complex-valued B-spline neural networks for modeling and inverting Hammerstein systems

Many communication signal processing applications involve modelling and inverting complex-valued (CV) Hammerstein systems. We develops a new CV B-spline neural network approach for efficient identification of the CV Hammerstein system and effective inversion of the estimated CV Hammerstein model. Specifically, the CV nonlinear static function in the Hammerstein system is represented using the tensor product from two univariate B-spline neural networks. An efficient alternating least squares estimation method is adopted for identifying the CV linear dynamic model’s coefficients and the CV B-spline neural network’s weights, which yields the closed-form solutions for both the linear dynamic model’s coefficients and the B-spline neural network’s weights, and this estimation process is guaranteed to converge very fast to a unique minimum solution. Furthermore, an accurate inversion of the CV Hammerstein system can readily be obtained using the estimated model. In particular, the inversion of the CV nonlinear static function in the Hammerstein system can be calculated effectively using a Gaussian-Newton algorithm, which naturally incorporates the efficient De Boor algorithm with both the B-spline curve and first order derivative recursions. The effectiveness of our approach is demonstrated using the application to equalisation of Hammerstein channels.

Prediction of nutrient digestibility and energy concentrations in fresh grass using nutrient composition

Improved nutrient utilization efficiency is strongly related to enhanced economic performance and reduced environmental footprint of dairy farms. Pasture-based systems are widely used for dairy production in certain areas of the world, but prediction equations of fresh grass nutritive value (nutrient digestibility and energy concentrations) are limited. Equations to predict digestible energy (DE) and metabolizable energy (ME) used for grazing cattle have been either developed with cattle fed conserved forage and concentrate diets or sheep fed previously frozen grass, and the majority of them require measurements less commonly available to producers, such as nutrient digestibility. The aim of the present study was therefore to develop prediction equations more suitable to grazing cattle for nutrient digestibility and energy concentrations, which are routinely available at farm level by using grass nutrient contents as predictors. A study with 33 nonpregnant, nonlactating cows fed solely fresh-cut grass at maintenance energy level for 50 wk was carried out over 3 consecutive grazing seasons. Freshly harvested grass of 3 cuts (primary growth and first and second regrowth), 9 fertilizer input levels, and contrasting stage of maturity (3 to 9 wk after harvest) was used, thus ensuring a wide representation of nutritional quality. As a result, a large variation existed in digestibility of dry matter (0.642-0.900) and digestible organic matter in dry matter (0.636-0.851) and in concentrations of DE (11.8-16.7 MJ/kg of dry matter) and ME (9.0-14.1 MJ/kg of dry matter). Nutrient digestibilities and DE and ME concentrations were negatively related to grass neutral detergent fiber (NDF) and acid detergent fiber (ADF) contents but positively related to nitrogen (N), gross energy, and ether extract (EE) contents. For each predicted variable (nutrient digestibilities or energy concentrations), different combinations of predictors (grass chemical composition) were found to be significant and increase the explained variation. For example, relatively higher R(2) values were found for prediction of N digestibility using N and EE as predictors; gross-energy digestibility using EE, NDF, ADF, and ash; NDF, ADF, and organic matter digestibilities using N, water-soluble carbohydrates, EE, and NDF; digestible organic matter in dry matter using water-soluble carbohydrates, EE, NDF, and ADF; DE concentration using gross energy, EE, NDF, ADF, and ash; and ME concentration using N, EE, ADF, and ash. Equations presented may allow a relatively quick and easy prediction of grass quality and, hence, better grazing utilization on commercial and research farms, where nutrient composition falls within the range assessed in the current study.

Quartic anharmonic resonances in acetylene molecules

Formulas are derived for the quartic anharmonic resonance coefficients observed to be important between C–H stretching and the combination of one quantum of C≡C stretching and two quanta of H–C≡C bending in a number of acetylene molecules. Examples of this resonance are ν3 with ν2+ν4+ν5 in 12C2H2, ν1 with ν2+2ν5 in 13C2H2, and ν1 with ν2+2ν4 in monofluoroacetylene and monochloroacetylene. The coefficients characterizing the resonances in these examples, which we denote K3,245, K1,255, and K1,244, arise from cubic and quartic terms in the anharmonic force field, in the normal coordinate representation, through second order and first order perturbation treatments respectively, where the second order resonances are calculated by a Van Vleck resonance formalism. The experimentally determined values of these coefficients are compared with values calculated from model anharmonic force fields.

Breeding success of Spotted Flycatchers Muscicapa striata in southern England - is woodland a good habitat for this species?

The UK population of the Spotted Flycatcher Muscicapa striata has declined markedly in the last 30 years but there have been few recent studies of the species. This study examined the relationship between nest success and the predominant habitat type around Spotted Flycatcher nests in two contrasting areas of England. A breeding population in eastern England, a region where numbers of Spotted Flycatchers are known to have decreased dramatically in recent decades, was compared with another in southwest England, where numbers have remained stable or even increased. Whilst there was no difference in breeding success between the two study areas, there were significant differences between habitats, with garden nests more successful than those in farmland or woodland, at both egg and chick stages. Estimates of productivity per nesting attempt were also lower in farmland and woodland, with nests in gardens fledging twice as many chicks as those in either woodland or farmland. The proximate cause of lower success in farmland and woodland was higher nest predation rates during both egg and chick stages. In terms of nesting success, farmland and woodland appear to be similar in quality for this species, but both appear to be suboptimal habitats when compared with gardens, providing evidence of a problem on the breeding grounds for this species, in at least these two habitats.

Breeding success of Spotted Flycatchers Muscicapa striata in southern England - is woodland a good habitat for this species?

The UK population of the Spotted Flycatcher Muscicapa striata has declined markedly in the last 30 years but there have been few recent studies of the species. This study examined the relationship between nest success and the predominant habitat type around Spotted Flycatcher nests in two contrasting areas of England. A breeding population in eastern England, a region where numbers of Spotted Flycatchers are known to have decreased dramatically in recent decades, was compared with another in southwest England, where numbers have remained stable or even increased. Whilst there was no difference in breeding success between the two study areas, there were significant differences between habitats, with garden nests more successful than those in farmland or woodland, at both egg and chick stages. Estimates of productivity per nesting attempt were also lower in farmland and woodland, with nests in gardens fledging twice as many chicks as those in either woodland or farmland. The proximate cause of lower success in farmland and woodland was higher nest predation rates during both egg and chick stages. In terms of nesting success, farmland and woodland appear to be similar in quality for this species, but both appear to be suboptimal habitats when compared with gardens, providing evidence of a problem on the breeding grounds for this species, in at least these two habitats.

Genotypic sex determination enabled adaptive radiations of extinct marine reptiles

Adaptive radiations often follow the evolution of key traits, such as the origin of the amniotic egg and the subsequent radiation of terrestrial vertebrates. The mechanism by which a species determines the sex of its offspring has been linked to critical ecological and life-history traits(1-3) but not to major adaptive radiations, in part because sex-determining mechanisms do not fossilize. Here we establish a previously unknown coevolutionary relationship in 94 amniote species between sex-determining mechanism and whether a species bears live young or lays eggs. We use that relationship to predict the sex-determining mechanism in three independent lineages of extinct Mesozoic marine reptiles (mosasaurs, sauropterygians and ichthyosaurs), each of which is known from fossils to have evolved live birth(4-7). Our results indicate that each lineage evolved genotypic sex determination before acquiring live birth. This enabled their pelagic radiations, where the relatively stable temperatures of the open ocean constrain temperature-dependent sex determination in amniote species. Freed from the need to move and nest on land(4,5,8), extreme physical adaptations to a pelagic lifestyle evolved in each group, such as the fluked tails, dorsal fins and wing-shaped limbs of ichthyosaurs. With the inclusion of ichthyosaurs, mosasaurs and sauropterygians, genotypic sex determination is present in all known fully pelagic amniote groups (sea snakes, sirenians and cetaceans), suggesting that this mode of sex determination and the subsequent evolution of live birth are key traits required for marine adaptive radiations in amniote lineages.

Assessment of the RE(OH)3 Ising magnetic materials as possible candidates for the study of transverse-field-induced quantum phase transitions

The LiHoxY1−xF4 Ising magnetic material subject to a magnetic field perpendicular to the Ho3+ Ising direction has shown over the past 20 years to be a host of very interesting thermodynamic and magnetic phenomena. Unfortunately, the availability of other magnetic materials other than LiHoxY1−xF4 that may be described by a transverse-field Ising model remains very much limited. It is in this context that we use here a mean-field theory to investigate the suitability of the Ho(OH)3, Dy(OH)3, and Tb(OH)3 insulating hexagonal dipolar Ising-type ferromagnets for the study of the quantum phase transition induced by a magnetic field, Bx, applied perpendicular to the Ising spin direction. Experimentally, the zero-field critical (Curie) temperatures are known to be Tc≈2.54, 3.48, and 3.72 K, for Ho(OH)3, Dy(OH)3, and Tb(OH)3, respectively. From our calculations we estimate the critical transverse field, Bxc, to destroy ferromagnetic order at zero temperature to be Bxc=4.35, 5.03, and 54.81 T for Ho(OH)3, Dy(OH)3, and Tb(OH)3, respectively. We find that Ho(OH)3, similarly to LiHoF4, can be quantitatively described by an effective S=1/2 transverse-field Ising model. This is not the case for Dy(OH)3 due to the strong admixing between the ground doublet and first excited doublet induced by the dipolar interactions. Furthermore, we find that the paramagnetic (PM) to ferromagnetic (FM) transition in Dy(OH)3 becomes first order for strong Bx and low temperatures. Hence, the PM to FM zero-temperature transition in Dy(OH)3 may be first order and not quantum critical. We investigate the effect of competing antiferromagnetic nearest-neighbor exchange and applied magnetic field, Bz, along the Ising spin direction ẑ on the first-order transition in Dy(OH)3. We conclude from these preliminary calculations that Ho(OH)3 and Dy(OH)3 and their Y3+ diamagnetically diluted variants, HoxY1−x(OH)3 and DyxY1−x(OH)3, are potentially interesting systems to study transverse-field-induced quantum fluctuations effects in hard axis (Ising-type) magnetic materials.

The economic consequences of mortgage debt

Unsurprisingly, a great deal of attention has been paid to the economic consequences of the credit crunch. However, this paper shows that the credit crunch was preceded by a strong build-up of mortgage debt internationally, which, in the long run, could turn out to be more significant than the credit crunch itself. Indeed, the debt build-up suggests that the credit crunch is more likely to reoccur, because highly-indebted households have weaker buffers to withstand unexpected shocks to their incomes or to interest rates. The paper presents a model that can explain the debt build-up and changes to the distribution of debt between existing owners and first-time buyers, which hinders access to home-ownership for the latter, even amongst those households who would be considered as credit-worthy.

Periodic time series modeling of environmental proxy records with guaranteed positive growth rate estimation

Identifying a periodic time-series model from environmental records, without imposing the positivity of the growth rate, does not necessarily respect the time order of the data observations. Consequently, subsequent observations, sampled in the environmental archive, can be inversed on the time axis, resulting in a non-physical signal model. In this paper an optimization technique with linear constraints on the signal model parameters is proposed that prevents time inversions. The activation conditions for this constrained optimization are based upon the physical constraint of the growth rate, namely, that it cannot take values smaller than zero. The actual constraints are defined for polynomials and first-order splines as basis functions for the nonlinear contribution in the distance-time relationship. The method is compared with an existing method that eliminates the time inversions, and its noise sensitivity is tested by means of Monte Carlo simulations. Finally, the usefulness of the method is demonstrated on the measurements of the vessel density, in a mangrove tree, Rhizophora mucronata, and the measurement of Mg/Ca ratios, in a bivalve, Mytilus trossulus.

Literature for the very young

This chapter looks at books which will support the very young through their earliest childhood and first few years of schooling. Learning to read can be hard but it is a lot harder if you never encounter the sort of literature that will engage you and motivate you to see the whole learning process as worthwhile. The chapter considers ways to share great texts with young children and how to select good books to share.

An Earth Observation Land Data Assimilation System (EO-LDAS)

Current methods for estimating vegetation parameters are generally sub-optimal in the way they exploit information and do not generally consider uncertainties. We look forward to a future where operational dataassimilation schemes improve estimates by tracking land surface processes and exploiting multiple types of observations. Dataassimilation schemes seek to combine observations and models in a statistically optimal way taking into account uncertainty in both, but have not yet been much exploited in this area. The EO-LDAS scheme and prototype, developed under ESA funding, is designed to exploit the anticipated wealth of data that will be available under GMES missions, such as the Sentinel family of satellites, to provide improved mapping of land surface biophysical parameters. This paper describes the EO-LDAS implementation, and explores some of its core functionality. EO-LDAS is a weak constraint variational dataassimilationsystem. The prototype provides a mechanism for constraint based on a prior estimate of the state vector, a linear dynamic model, and EarthObservationdata (top-of-canopy reflectance here). The observation operator is a non-linear optical radiative transfer model for a vegetation canopy with a soil lower boundary, operating over the range 400 to 2500 nm. Adjoint codes for all model and operator components are provided in the prototype by automatic differentiation of the computer codes. In this paper, EO-LDAS is applied to the problem of daily estimation of six of the parameters controlling the radiative transfer operator over the course of a year (> 2000 state vector elements). Zero and first order process model constraints are implemented and explored as the dynamic model. The assimilation estimates all state vector elements simultaneously. This is performed in the context of a typical Sentinel-2 MSI operating scenario, using synthetic MSI observations simulated with the observation operator, with uncertainties typical of those achieved by optical sensors supposed for the data. The experiments consider a baseline state vector estimation case where dynamic constraints are applied, and assess the impact of dynamic constraints on the a posteriori uncertainties. The results demonstrate that reductions in uncertainty by a factor of up to two might be obtained by applying the sorts of dynamic constraints used here. The hyperparameter (dynamic model uncertainty) required to control the assimilation are estimated by a cross-validation exercise. The result of the assimilation is seen to be robust to missing observations with quite large data gaps.

Extended factorizations of exponential functionals of Lévy processes

Pardo, Patie, and Savov derived, under mild conditions, a Wiener-Hopf type factorization for the exponential functional of proper Lévy processes. In this paper, we extend this factorization by relaxing a finite moment assumption as well as by considering the exponential functional for killed Lévy processes. As a by-product, we derive some interesting fine distributional properties enjoyed by a large class of this random variable, such as the absolute continuity of its distribution and the smoothness, boundedness or complete monotonicity of its density. This type of results is then used to derive similar properties for the law of maxima and first passage time of some stable Lévy processes. Thus, for example, we show that for any stable process with $\rho\in(0,\frac{1}{\alpha}-1]$, where $\rho\in[0,1]$ is the positivity parameter and $\alpha$ is the stable index, then the first passage time has a bounded and non-increasing density on $\mathbb{R}_+$. We also generate many instances of integral or power series representations for the law of the exponential functional of Lévy processes with one or two-sided jumps. The proof of our main results requires different devices from the one developed by Pardo, Patie, Savov. It relies in particular on a generalization of a transform recently introduced by Chazal et al together with some extensions to killed Lévy process of Wiener-Hopf techniques. The factorizations developed here also allow for further applications which we only indicate here also allow for further applications which we only indicate here.

Functional model for catecholase-like activity: A mechanistic approach with manganese(III) complexes of salen type Schiff base ligands

Three new Mn(III) complexes [MnL1(OOCH)(OH2)] (1), [MnL2(OH2)(2)][Mn2L22(NO2)(3)] (2) and [Mn2L21(NO2)(2)] (3) (where H2L1 = H(2)Me(2)Salen = 2,7-bis(2-hydroxyphenyl)-2,6-diazaocta-2,6-diene and H2L2 = H(2)Salpn = 1,7-bis(2-hydroxyphenyl)-2,6-diazahepta-1,6-diene) have been synthesized. X-ray crystal structure analysis reveals that 1 is a mononuclear species whereas 2 contains a mononuclear cationic and a dinuclear nitrite bridged (mu-1 kappa O:2 kappa O') anionic unit. Complex 3 is a phenoxido bridged dimer containing terminally coordinated nitrite. Complexes 1-3 show excellent catecholase-like activity with 3,5-di-tert-butylcatechol (3,5-DTBC) as the substrate. Kinetic measurements suggest that the rate of catechol oxidation follows saturation kinetics with respect to the substrate and first order kinetics with respect to the catalyst. Formation of bis(mu-oxo)dimanganese(III,III) as an intermediate during the course of reaction is identified from ESI-MS spectra. The characteristic six line EPR spectra of complex 2 in the presence of 3,5-DTBC supports the formation of manganese(II)-semiquinonate as an intermediate species during the catalytic oxidation of 3,5-DTBC.

Estimates of aerosol radiative forcing from the MACC re-analysis

The European Centre for Medium-range Weather Forecast (ECMWF) provides an aerosol re-analysis starting from year 2003 for the Monitoring Atmospheric Composition and Climate (MACC) project. The re-analysis assimilates total aerosol optical depth retrieved by the Moderate Resolution Imaging Spectroradiometer (MODIS) to correct for model departures from observed aerosols. The reanalysis therefore combines satellite retrievals with the full spatial coverage of a numerical model. Re-analysed products are used here to estimate the shortwave direct and first indirect radiative forcing of anthropogenic aerosols over the period 2003–2010, using methods previously applied to satellite retrievals of aerosols and clouds. The best estimate of globally-averaged, all-sky direct radiative forcing is −0.7±0.3Wm−2. The standard deviation is obtained by a Monte-Carlo analysis of uncertainties, which accounts for uncertainties in the aerosol anthropogenic fraction, aerosol absorption, and cloudy-sky effects. Further accounting for differences between the present-day natural and pre-industrial aerosols provides a direct radiative forcing estimate of −0.4±0.3Wm−2. The best estimate of globally-averaged, all-sky first indirect radiative forcing is −0.6±0.4Wm−2. Its standard deviation accounts for uncertainties in the aerosol anthropogenic fraction, and in cloud albedo and cloud droplet number concentration susceptibilities to aerosol changes. The distribution of first indirect radiative forcing is asymmetric and is bounded by −0.1 and −2.0Wm−2. In order to decrease uncertainty ranges, better observational constraints on aerosol absorption and sensitivity of cloud droplet number concentrations to aerosol changes are required.

Impact of the modal aerosol scheme GLOMAP-mode on aerosol forcing in the Hadley Centre Global Environmental Model

The Hadley Centre Global Environmental Model (HadGEM) includes two aerosol schemes: the Coupled Large-scale Aerosol Simulator for Studies in Climate (CLASSIC), and the new Global Model of Aerosol Processes (GLOMAP-mode). GLOMAP-mode is a modal aerosol microphysics scheme that simulates not only aerosol mass but also aerosol number, represents internally-mixed particles, and includes aerosol microphysical processes such as nucleation. In this study, both schemes provide hindcast simulations of natural and anthropogenic aerosol species for the period 2000–2006. HadGEM simulations of the aerosol optical depth using GLOMAP-mode compare better than CLASSIC against a data-assimilated aerosol re-analysis and aerosol ground-based observations. Because of differences in wet deposition rates, GLOMAP-mode sulphate aerosol residence time is two days longer than CLASSIC sulphate aerosols, whereas black carbon residence time is much shorter. As a result, CLASSIC underestimates aerosol optical depths in continental regions of the Northern Hemisphere and likely overestimates absorption in remote regions. Aerosol direct and first indirect radiative forcings are computed from simulations of aerosols with emissions for the year 1850 and 2000. In 1850, GLOMAP-mode predicts lower aerosol optical depths and higher cloud droplet number concentrations than CLASSIC. Consequently, simulated clouds are much less susceptible to natural and anthropogenic aerosol changes when the microphysical scheme is used. In particular, the response of cloud condensation nuclei to an increase in dimethyl sulphide emissions becomes a factor of four smaller. The combined effect of different 1850 baselines, residence times, and abilities to affect cloud droplet number, leads to substantial differences in the aerosol forcings simulated by the two schemes. GLOMAP-mode finds a presentday direct aerosol forcing of −0.49Wm−2 on a global average, 72% stronger than the corresponding forcing from CLASSIC. This difference is compensated by changes in first indirect aerosol forcing: the forcing of −1.17Wm−2 obtained with GLOMAP-mode is 20% weaker than with CLASSIC. Results suggest that mass-based schemes such as CLASSIC lack the necessary sophistication to provide realistic input to aerosol-cloud interaction schemes. Furthermore, the importance of the 1850 baseline highlights how model skill in predicting present-day aerosol does not guarantee reliable forcing estimates. Those findings suggest that the more complex representation of aerosol processes in microphysical schemes improves the fidelity of simulated aerosol forcings.