A note on the likelihood and moments of the skew-normal distribution

In this paper an alternative approach to the one in Henze (1986) is proposed for deriving the odd moments of the skew-normal distribution considered in Azzalini (1985). The approach is based on a Pascal type triangle, which seems to greatly simplify moments computation. Moreover, it is shown that the likelihood equation for estimating the asymmetry parameter in such model is generated as orthogonal functions to the sample vector. As a consequence, conditions for a unique solution of the likelihood equation are established, which seem to hold in more general setting.

Application of electron paramagnetic resonance Spectroscopy for validation of the novel (AN plus DN) solvent polarity scale

Based on solvation studies of polymers, the sum (1: 1) of the electron acceptor (AN) and electron donor (DN) values of solvents has been proposed as an alternative polarity scale. To test this, the electron paramagnetic resonance isotropic hyperfine splitting constant, a parameter known to be dependent on the polarity/proticity of the medium, was correlated with the (AN+DN) term using three paramagnetic probes. The linear regression coefficient calculated for 15 different solvents was approximately 0.9, quite similar to those of other well-known polarity parameters, attesting to the validity of the (AN+DN) term as a novel ""two-parameter"" solvent polarity scale.

Structure, stability, depolarized light scattering, and vibrational spectra of fullerenols from all-electron density-functional-theory calculations

We have investigated the stability, electronic properties, Rayleigh (elastic), and Raman (inelastic) depolarization ratios, infrared and Raman absorption vibrational spectra of fullerenols [C(60)(OH)(n)] with different degrees of hydroxylation by using all-electron density-functional-theory (DFT) methods. Stable arrangements of these molecules were found by means of full geometry optimizations using Becke's three-parameter exchange functional with the Lee, Yang, and Parr correlation functional. This DFT level has been combined with the 6-31G(d,p) Gaussian-type basis set, as a compromise between accuracy and capability to treat highly hydroxylated fullerenes, e.g., C(60)(OH)(36). Thus, the molecular properties of fullerenols were systematically analyzed for structures with n=1, 2, 3, 4, 8, 10, 16, 18, 24, 32, and 36. From the electronic structure analysis of these molecules, we have evidenced an important effect related to the weak chemical reactivity of a possible C(60)(OH)(24) isomer. To investigate Raman scattering and the vibrational spectra of the different fullerenols, frequency calculations are carried out within the harmonic approximation. In this case a systematic study is only performed for n=1-4, 8, 10, 16, 18, and 24. Our results give good agreements with the expected changes in the spectral absorptions due to the hydroxylation of fullerenes.

Self-organized distribution of periodicity and chaos in an electrochemical oscillator

We report a detailed numerical investigation of a prototype electrochemical oscillator, in terms of high-resolution phase diagrams for an experimentally relevant section of the control (parameter) space. The prototype model consists of a set of three autonomous ordinary differential equations which captures the general features of electrochemical oscillators characterized by a partially hidden negative differential resistance in an N-shaped current-voltage stationary curve. By computing Lyapunov exponents, we provide a detailed discrimination between chaotic and periodic phases of the electrochemical oscillator. Such phases reveal the existence of an intricate structure of domains of periodicity self-organized into a chaotic background. Shrimp-like periodic regions previously observed in other discrete and continuous systems were also observed here, which corroborate the universal nature of the occurrence of such structures. In addition, we have also found a structured period distribution within the order region. Finally we discuss the possible experimental realization of comparable phase diagrams.

Complex kinetics, high frequency oscillations and temperature compensation in the electro-oxidation of ethylene glycol on platinum

Despite the fact that the majority of the catalytic electro-oxidation of small organic molecules presents oscillatory kinetics under certain conditions, there are few systematic studies concerning the influence of experimental parameters on the oscillatory dynamics. Of the studies available, most are devoted to C1 molecules and just some scattered data are available for C2 molecules. We present in this work a comprehensive study of the electro-oxidation of ethylene glycol on polycrystalline platinum surfaces and in alkaline media. The system was studied by means of electrochemical impedance spectroscopy, cyclic voltammetry, and chronoamperometry, and the impact of parameters such as applied current, ethylene glycol concentration, and temperature were investigated. As in the case of other parent systems, the instabilities in this system were associated with a hidden negative differential resistance, as identified by impedance data. Very rich and robust dynamics were observed, including the presence of harmonic and mixed mode oscillations and chaotic states, in some parameter region. Oscillation frequencies of about 16 Hz characterized the fastest oscillations ever reported for the electro-oxidation of small organic molecules. Those high frequencies were strongly influenced by the electrolyte pH and far less affected by the EG concentration. The system was regularly dependent on temperature under voltammetric conditions but rather independent within the oscillatory regime.

Conversion to soy on the Amazonian agricultural frontier increases streamflow without affecting stormflow dynamics

Large-scale soy agriculture in the southern Brazilian Amazon now rivals deforestation for pasture as the region`s predominant form of land use change. Such landscape-level change can have substantial consequences for local and regional hydrology, but these effects remain relatively unstudied in this ecologically and economically important region. We examined how the conversion to soy agriculture influences water balances and stormflows using stream discharge (water yields) and the timing of discharge (stream hydrographs) in small (2.5-13.5 km2) forested and soy headwater watersheds in the Upper Xingu Watershed in the state of Mato Grosso, Brazil. We monitored water yield for 1 year in three forested and four soy watersheds. Mean daily water yields were approximately four times higher in soy than forested watersheds, and soy watersheds showed greater seasonal variability in discharge. The contribution of stormflows to annual streamflow in all streams was low (< 13% of annual streamflow), and the contribution of stormflow to streamflow did not differ between land uses. If the increases in water yield observed in this study are typical, landscape-scale conversion to soy substantially alters water-balance, potentially altering the regional hydrology over large areas of the southern Amazon.

Influence of land-use change on near-surface hydrological processes: Undisturbed forest to pasture

Soil compaction that follows the clearing of tropical forest for cattle pasture is associated with lower soil hydraulic conductivity and increased frequency and volume of overland flow. We investigated the frequency of perched water tables, overland flow and stormflow in an Amazon forest and in an adjacent 25-year-old pasture cleared from the same forest. We compared the results with the frequencies of these phenomena estimated from comparisons of rainfall intensity and soil hydraulic conductivity. The frequency of perched water tables based on rainfall intensity and soil hydraulic conductivity was expected to double in pasture compared with forest. This corresponded closely with an approximate doubling of the frequency of stormflow and overland flow in pasture. In contrast, the stormflow volume in pasture increased 17-fold. This disproportional increase of stormflow resulted from overland flow generation over large areas of pasture, while overland flow generation in the forest was spatially limited and was observed only very near the stream channel. In both catchments, stormflow was generated by saturation excess because of perched water tables and near-surface groundwater levels. Stormflow was occasionally generated in the forest by rapid return flow from macropores, while slow return flow from a continuous perched water table was more common in the pasture. These results suggest that deforestation for pasture alters fundamental mechanisms of stormflow generation and may increase runoff volumes over wide regions of Amazonia. (C) 2009 Elsevier B.V. All rights reserved.

Simulation of sugarcane residue decomposition and aboveground growth

Due to the worldwide increase in demand for biofuels, the area cultivated with sugarcane is expected to increase. For environmental and economic reasons, an increasing proportion of the areas are being harvested without burning, leaving the residues on the soil surface. This periodical input of residues affects soil physical, chemical and biological properties, as well as plant growth and nutrition. Modeling can be a useful tool in the study of the complex interactions between the climate, residue quality, and the biological factors controlling plant growth and residue decomposition. The approach taken in this work was to parameterize the CENTURY model for the sugarcane crop, to simulate the temporal dynamics of aboveground phytomass and litter decomposition, and to validate the model through field experiment data. When studying aboveground growth, burned and unburned harvest systems were compared, as well as the effect of mineral fertilizer and organic residue applications. The simulations were performed with data from experiments with different durations, from 12 months to 60 years, in Goiana, TimbaA(0)ba and Pradpolis, Brazil; Harwood, Mackay and Tully, Australia; and Mount Edgecombe, South Africa. The differentiation of two pools in the litter, with different decomposition rates, was found to be a relevant factor in the simulations made. Originally, the model had a basically unlimited layer of mulch directly available for decomposition, 5,000 g m(-2). Through a parameter optimization process, the thickness of the mulch layer closer to the soil, more vulnerable to decomposition, was set as 110 g m(-2). By changing the layer of mulch at any given time available for decomposition, the sugarcane residues decomposition simulations where close to measured values (R (2) = 0.93), contributing to making the CENTURY model a tool for the study of sugarcane litter decomposition patterns. The CENTURY model accurately simulated aboveground carbon stalk values (R (2) = 0.76), considering burned and unburned harvest systems, plots with and without nitrogen fertilizer and organic amendment applications, in different climates and soil conditions.

Biomass burning in Brazil`s Amazonian ""arc of deforestation"": Burning efficiency and charcoal formation in a fire after mechanized clearing at Feliz Natal, Mato Grosso

Estimates of greenhouse-gas emissions from deforestation are highly uncertain because of high variability in key parameters and because of the limited number of studies providing field measurements of these parameters. One such parameter is burning efficiency, which determines how much of the original forest`s aboveground carbon stock will be released in the burn, as well as how much will later be released by decay and how much will remain as charcoal. In this paper we examined the fate of biomass from a semideciduous tropical forest in the ""arc of deforestation,"" where clearing activity is concentrated along the southern edge of the Amazon forest. We estimated carbon content, charcoal formation and burning efficiency by direct measurements (cutting and weighing) and by line-intersect sampling (LIS) done along the axis of each plot before and after burning of felled vegetation. The total aboveground dry biomass found here (219.3 Mg ha(-1)) is lower than the values found in studies that have been done in other parts of the Amazon region. Values for burning efficiency (65%) and charcoal formation (6.0%, or 5.98 Mg C ha(-1)) were much higher than those found in past studies in tropical areas. The percentage of trunk biomass lost in burning (49%) was substantially higher than has been found in previous studies. This difference may be explained by the concentration of more stems in the smaller diameter classes and the low humidity of the fuel (the dry season was unusually long in 2007, the year of the burn). This study provides the first measurements of forest burning parameters for a group of forest types that is now undergoing rapid deforestation. The burning parameters estimated here indicate substantially higher burning efficiency than has been found in other Amazonian forest types. Quantification of burning efficiency is critical to estimates of trace-gas emissions from deforestation. (C) 2009 Elsevier B.V. All rights reserved.

Gaseous and fluvial carbon export from an Amazon forest watershed

The transfer of carbon (C) from Amazon forests to aquatic ecosystems as CO(2) supersaturated in groundwater that outgases to the atmosphere after it reaches small streams has been postulated to be an important component of terrestrial ecosystem C budgets. We measured C losses as soil respiration and methane (CH(4)) flux, direct CO(2) and CH(4) fluxes from the stream surface and fluvial export of dissolved inorganic C (DIC), dissolved organic C (DOC), and particulate C over an annual hydrologic cycle from a 1,319-ha forested Amazon perennial first-order headwater watershed at Tanguro Ranch in the southern Amazon state of Mato Grosso. Stream pCO(2) concentrations ranged from 6,491 to 14,976 mu atm and directly-measured stream CO(2) outgassing flux was 5,994 +/- A 677 g C m(-2) y(-1) of stream surface. Stream pCH(4) concentrations ranged from 291 to 438 mu atm and measured stream CH(4) outgassing flux was 987 +/- A 221 g C m(-2) y(-1). Despite high flux rates from the stream surface, the small area of stream itself (970 m(2), or 0.007% of watershed area) led to small directly-measured annual fluxes of CO(2) (0.44 +/- A 0.05 g C m(2) y(-1)) and CH(4) (0.07 +/- A 0.02 g C m(2) y(-1)) per unit watershed land area. Measured fluvial export of DIC (0.78 +/- A 0.04 g C m(-2) y(-1)), DOC (0.16 +/- A 0.03 g C m(-2) y(-1)) and coarse plus fine particulate C (0.001 +/- A 0.001 g C m(-2) y(-1)) per unit watershed land area were also small. However, stream discharge accounted for only 12% of the modeled annual watershed water output because deep groundwater flows dominated total runoff from the watershed. When C in this bypassing groundwater was included, total watershed export was 10.83 g C m(-2) y(-1) as CO(2) outgassing, 11.29 g C m(-2) y(-1) as fluvial DIC and 0.64 g C m(-2) y(-1) as fluvial DOC. Outgassing fluxes were somewhat lower than the 40-50 g C m(-2) y(-1) reported from other Amazon watersheds and may result in part from lower annual rainfall at Tanguro. Total stream-associated gaseous C losses were two orders of magnitude less than soil respiration (696 +/- A 147 g C m(-2) y(-1)), but total losses of C transported by water comprised up to about 20% of the +/- A 150 g C m(-2) (+/- 1.5 Mg C ha(-1)) that is exchanged annually across Amazon tropical forest canopies.

Discharge-calcium concentration relationships in streams of the Amazon and Cerrado of Brazil: soil or land use controlled

Stream discharge-concentration relationships are indicators of terrestrial ecosystem function. Throughout the Amazon and Cerrado regions of Brazil rapid changes in land use and land cover may be altering these hydrochemical relationships. The current analysis focuses on factors controlling the discharge-calcium (Ca) concentration relationship since previous research in these regions has demonstrated both positive and negative slopes in linear log(10)discharge-log(10)Ca concentration regressions. The objective of the current study was to evaluate factors controlling stream discharge-Ca concentration relationships including year, season, stream order, vegetation cover, land use, and soil classification. It was hypothesized that land use and soil class are the most critical attributes controlling discharge-Ca concentration relationships. A multilevel, linear regression approach was utilized with data from 28 streams throughout Brazil. These streams come from three distinct regions and varied broadly in watershed size (< 1 to > 10(6) ha) and discharge (10(-5.7)-10(3.2) m(3) s(-1)). Linear regressions of log(10)Ca versus log(10)discharge in 13 streams have a preponderance of negative slopes with only two streams having significant positive slopes. An ANOVA decomposition suggests the effect of discharge on Ca concentration is large but variable. Vegetation cover, which incorporates aspects of land use, explains the largest proportion of the variance in the effect of discharge on Ca followed by season and year. In contrast, stream order, land use, and soil class explain most of the variation in stream Ca concentration. In the current data set, soil class, which is related to lithology, has an important effect on Ca concentration but land use, likely through its effect on runoff concentration and hydrology, has a greater effect on discharge-concentration relationships.

A new monitor for routine thermal and epithermal neutron fluence rate monitoring in k(0) INAA

The Zr-Au set for monitoring the thermal and epithermal neutron fluence rate and the epithermal spectrum parameter a is not always practicable for routine application of INAA in well-thermalized facilities. An alternative set consisting of Cr, Au and Mo provides values for the thermal neutron fluence rate, f and alpha that are not significantly different from those found via the Zr-Au method and the Cd-covered Zr-method. The IRMM standard SMELS-II was analyzed using the (Au-Cr-Mo) monitor and a good agreement was obtained. (C) 2008 Elsevier Ltd. All rights reserved.

Kernel machines for epilepsy diagnosis via EEG signal classification: A comparative study

Objective: We carry out a systematic assessment on a suite of kernel-based learning machines while coping with the task of epilepsy diagnosis through automatic electroencephalogram (EEG) signal classification. Methods and materials: The kernel machines investigated include the standard support vector machine (SVM), the least squares SVM, the Lagrangian SVM, the smooth SVM, the proximal SVM, and the relevance vector machine. An extensive series of experiments was conducted on publicly available data, whose clinical EEG recordings were obtained from five normal subjects and five epileptic patients. The performance levels delivered by the different kernel machines are contrasted in terms of the criteria of predictive accuracy, sensitivity to the kernel function/parameter value, and sensitivity to the type of features extracted from the signal. For this purpose, 26 values for the kernel parameter (radius) of two well-known kernel functions (namely. Gaussian and exponential radial basis functions) were considered as well as 21 types of features extracted from the EEG signal, including statistical values derived from the discrete wavelet transform, Lyapunov exponents, and combinations thereof. Results: We first quantitatively assess the impact of the choice of the wavelet basis on the quality of the features extracted. Four wavelet basis functions were considered in this study. Then, we provide the average accuracy (i.e., cross-validation error) values delivered by 252 kernel machine configurations; in particular, 40%/35% of the best-calibrated models of the standard and least squares SVMs reached 100% accuracy rate for the two kernel functions considered. Moreover, we show the sensitivity profiles exhibited by a large sample of the configurations whereby one can visually inspect their levels of sensitiveness to the type of feature and to the kernel function/parameter value. Conclusions: Overall, the results evidence that all kernel machines are competitive in terms of accuracy, with the standard and least squares SVMs prevailing more consistently. Moreover, the choice of the kernel function and parameter value as well as the choice of the feature extractor are critical decisions to be taken, albeit the choice of the wavelet family seems not to be so relevant. Also, the statistical values calculated over the Lyapunov exponents were good sources of signal representation, but not as informative as their wavelet counterparts. Finally, a typical sensitivity profile has emerged among all types of machines, involving some regions of stability separated by zones of sharp variation, with some kernel parameter values clearly associated with better accuracy rates (zones of optimality). (C) 2011 Elsevier B.V. All rights reserved.

Homogenization in a thin domain with an oscillatory boundary

In this paper we analyze the behavior of the Laplace operator with Neumann boundary conditions in a thin domain of the type R(epsilon) = {(x(1), x(2)) is an element of R(2) vertical bar x(1) is an element of (0, 1), 0 < x(2) < epsilon G(x(1), x(1)/epsilon)} where the function G(x, y) is periodic in y of period L. Observe that the upper boundary of the thin domain presents a highly oscillatory behavior and, moreover, the height of the thin domain, the amplitude and period of the oscillations are all of the same order, given by the small parameter epsilon. (C) 2011 Elsevier Masson SAS. All rights reserved.

Analytical calculation of the transition to complete phase synchronization in coupled oscillators

Here we present a system of coupled phase oscillators with nearest neighbors coupling, which we study for different boundary conditions. We concentrate at the transition to the total synchronization. We are able to develop exact solutions for the value of the coupling parameter when the system becomes completely synchronized, for the case of periodic boundary conditions as well as for a chain with fixed ends. We compare the results with those calculated numerically.