General results for the beta-modified Weibull distribution

We study in detail the so-called beta-modified Weibull distribution, motivated by the wide use of the Weibull distribution in practice, and also for the fact that the generalization provides a continuous crossover towards cases with different shapes. The new distribution is important since it contains as special sub-models some widely-known distributions, such as the generalized modified Weibull, beta Weibull, exponentiated Weibull, beta exponential, modified Weibull and Weibull distributions, among several others. It also provides more flexibility to analyse complex real data. Various mathematical properties of this distribution are derived, including its moments and moment generating function. We examine the asymptotic distributions of the extreme values. Explicit expressions are also derived for the chf, mean deviations, Bonferroni and Lorenz curves, reliability and entropies. The estimation of parameters is approached by two methods: moments and maximum likelihood. We compare by simulation the performances of the estimates from these methods. We obtain the expected information matrix. Two applications are presented to illustrate the proposed distribution.

Private valuation of carbon sequestration in forest plantations

Approval of the Clean Development Mechanism, provided for in the Kyoto Protocol, enables countries with afforested land to trade in carbon emissions reduction certificates related to carbon dioxide equivalent quantities (CO(2-e)) stored within a certain forest area. Potential CO(2-e) above base line sequestration was determined for two forest sites on commercial eucalyptus plantations in northern Brazil (Bahia). Compensation values for silvicultural regimes involving rotation lengths greater than economically optimal were computed using the Faustmann formula. Mean values obtained were US$8.16 (MgCO(2-e))(-1) and US $7.19 (MgCO(2-e))(-1) for average and high site indexes, respectively. Results show that carbon supply is more cost-efficient in highly productive sites. Annuities of US$18.8 Mg C(-1) and US$35.1 Mg C(-1) and yearly payments of US$4.4 m(-3) and US$8.2 m(-3) due for each marginal cubic meter produced were computed for high and average sites, respectively. The estimated value of the tonne of carbon defines minimum values to be paid to forest owners, in order to induce a change in silvicultural management regimes. A reduction of carbon supply could be expected as a result of an increase in wood prices, although it would not respond in a regular manner. For both sites, price elasticity of supply was found to be inelastic and increased as rotation length moved further away from economically optimal: 0.24 and 0.27 for age 11 years in average- and high-productivity sites, respectively. This would be due to biomass production potential as a limiting factor; beyond a certain threshold value. an increase in price does not sustain a proportional change in carbon storage supply. The environmental service valuation model proposed might be adequate for assessing potential supply in plantation forestry, from a private landowner perspective, with an economic opportunity cost. The model is not applicable to low commercial value forest plantations. (C) 2009 Elsevier B.V. All rights reserved.

Why don`t our stands grow even faster? Control of production and carbon cycling in eucalypt plantations

The growth of Eucalyptus stands varies several fold across sites, under the influence of resource availability, stand age and stand structure. We describe a series of related studies that aim to understand the mechanisms that drive this great range in stand growth rates. In a seven-year study in Hawaii of Eucalyptus saligna at a site that was not water limited, we showed that nutrient availability differences led to a two-fold difference in stand wood production. Increasing nutrient supply in mid-rotation raised productivity to the level attained in continuously fertilised plots. Fertility affected the age-related decline in wood and foliage production; production in the intensive fertility treatments declined more slowly than in the minimal fertility treatments. The decline in stem production was driven largely by a decline in canopy photosynthesis. Over time, the fraction of canopy photosynthesis partitioned to below-ground allocation increased, as did foliar respiration, further reducing wood production. The reason for the decline in photosynthesis was uncertain, but it was not caused by nutrient limitation, a decline in leaf area or in photosynthetic capacity, or by hydraulic limitation. Most of the increase in carbon stored from conversion of the sugarcane plantation to Eucalyptus plantation was in the above-ground woody biomass. Soil carbon showed no net change. This study and other studies on carbon allocation showed that resource availability changes the fraction of annual photosynthesis used below-ground and for wood production. High resources (nutrition or water) decrease the partitioning below-ground and increase partitioning to wood production. Annual foliage and wood respiration and foliage production as a fraction of annual photosynthesis was remarkably constant across a wide range of fertility treatments and forest age. In the Brazil Eucalyptus Productivity Project, stand structure was manipulated by planting clonal Eucalyptus all at once or in three groups at three-monthly intervals, producing a stand where trees did not segregate into dominants and one that had strong dominance. The uneven stand structure reduced production 10-15% throughout the rotation.

Production and carbon allocation in a clonal Eucalyptus plantation with water and nutrient manipulations

We examined resource limitations on growth and carbon allocation in a fast-growing, clonal plantation of Eucalyptus grandis x urophylla in Brazil by characterizing responses to annual rainfall, and response to irrigation and fertililization for 2 years. Productivity measures included gross primary production (GPP), total belowground carbon allocation (TBCA), bole growth, and net ecosystem production (NEP). Replicate plots within a single plantation were established at the midpoint of the rotation (end of year 3), with treatments of no additional fertilization or irrigation, heavy fertilization (to remove any nutrient limitation), irrigation (to remove any water limitation), and irrigation plus fertilization. Rainfall was unusually high in the first year (1769mm) of the experiment, and control plots had high rates of GPP (6.64 kg C m(-2) year(-1)), TBCA (2.14 kg C m(-2) year(-1)), and bole growth (1.81 kg C m(-2) year). Irrigation increased each of these rates by 15-17%. The second year of the experiment had average rainfall (1210 mm), and lower rainfall decreased production in control plots by 46% (GPP), 52% (TBCA), and 40% (bole growth). Fertilization treatments had neglible effects. The response to irrigation was much greater in the drier year, with irrigated plots exceeding the production in control plots by 83% (GPP), 239% (TBCA), and 24% (bole growth). Even though the rate of irrigation ensured no water limitation to tree growth, the high rainfall year showed higher production in irrigated plots for both GPP (38% greater than in drier year) and bole growth (23% greater). Varying humidity and supplies of water led to a range in NEP of 0.8-2.7 kg C m(-2) year. This difference between control and irrigated treatments, combined with differences between drier and wetter years, indicated a strong response of these Eucalyptus trees to both water supply and atmospheric humidity during the dry season. The efficiency of converting light energy into fixed carbon ranged from a low of 0.027 mol C to a high of 0.060 mol C per mol of absorbed photosynthetically active radiation (APAR), and the efficiency of bolewood production ranged from 0.78 to 1.98 g wood per MJ of APAR. Irrigation increased the efficiency of wood production per unit of water used from 2.55 kg wood m(-3) in the rainfed plot to 3.51 kg m(-3) in irrigated plots. Detailed information on the response of C budgets to environmental conditions and resource supplies will be necessary for accurate predictions of plantation yields across years and landscapes. (V) 2007 Elsevier B.V. All rights reserved.

Macroscopic root water uptake distribution using a matric flux potential approach

Hydrological models featuring root water uptake usually do not include compensation mechanisms such that reductions in uptake from dry layers are compensated by an increase in uptake from wetter layers. We developed a physically based root water uptake model with an implicit compensation mechanism. Based on an expression for the matric flux potential (M) as a function of the distance to the root, and assuming a depth-independent value of M at the root surface, uptake per layer is shown to be a function of layer bulk M, root surface M, and a weighting factor that depends on root length density and root radius. Actual transpiration can be calculated from the sum of layer uptake rates. The proposed reduction function (PRF) was built into the SWAP model, and predictions were compared to those made with the Feddes reduction function (FRF). Simulation results were tested against data from Canada (continuous spring wheat [(Triticum aestivum L.]) and Germany (spring wheat, winter barley [Hordeum vulgare L.], sugarbeet [Beta vulgaris L.], winter wheat rotation). For the Canadian data, the root mean square error of prediction (RMSEP) for water content in the upper soil layers was very similar for FRF and PRF; for the deeper layers, RMSEP was smaller for PRF. For the German data, RMSEP was lower for PRF in the upper layers and was similar for both models in the deeper layers. In conclusion, but dependent on the properties of the data sets available for testing,the incorporation of the new reduction function into SWAP was successful, providing new capabilities for simulating compensated root water uptake without increasing the number of input parameters or degrading model performance.

Greenhouse evaluation of the potential of sorghum, pearl millet and crotalaria in the management of Meloidogyne javanica

Crop rotation can play a valuable role in managing plant parasitic nematodes, depending on the availability of profitable non-host or poor host crops. Alternatively, non-host cover crops or green manures can be used in succession to summer cash Crops for this Purpose. The aim of the current study was to evaluate, under greenhouse conditions, the host status of commercial hybrids and cultivars of grain and silage sorghum (Sorghum bicolor) for Meloidogyne javanica, and to assess the effect of sorghum on nematode population in comparison with pearl millet (poor host for M. javanica), showy crotalaria and sunn hemp (both non-hosts). Based on two experiments, it was stated that, as a rule, grain sorghum is a poor host for M. javanica, but silage sorghum is a good host. Silage sorghum `BRS 601` was an exception. In other experiments, grain sorghum, pearl millet (Pennisetum glaucum `BN 2`), showy crotalaria (Crotolaria spectabilis `Comum`) and sunn hemp (C. juncea `IAC-KR-1`) reduced M. javanica population level, while silage sorghum increased the nematode density.

Surface Liming and Zinc Availability in a Long-Term Experiment under No-Till System

No-till (NT) system with crop rotation is one of the most effective strategies to improve agricultural sustainability in tropical and subtropical regions. To control soil acidity in NT, lime is broadcast on the surface without incorporation. The increase in soil pH due to surface liming may decrease zinc (Zn) availability and its uptake by crops. A field experiment was performed in Parana State, Brazil, on a loamy, kaolinitic, thermic Typic Hapludox to evaluate Zn bioavailability in a NT system after surface liming and re-liming. Dolomitic lime was surface applied on the main plots in July 1993 at the rates of 0, 2, 4, and 6 Mg ha-1. In June 2000, the main plots were divided in two subplots to study of the effect of surface re-liming at the rates of 0 and 3 Mg ha-1. The cropping sequence was soybean [Glycine max (L.) Merrill] (2001-2 and 2002-3), wheat (Triticum aestivum L.) (2003), soybean (2003-4), corn (Zea mays L.) (2004-5), and soybean (2005-6). Soil samples were collected at the following depths: 0-0.05, 0.05-0.10, and 0.10-0.20m, 10 years after surface liming and 3 years after surface re-liming. Soil Zn levels were extracted by four extractants: (i) 0.005molL-1 diethylenetriaminepentaacetic acid (DTPA) + 0.1molL-1 triethanolamine (TEA) + 0.01molL-1 calcium chloride (CaCl2) solution at pH7.3 (DTPA-TEA), (ii) 0.1molL-1 hydrochloric acid (HCl) solution, (iii) Mehlich 1 solution, and (iv) Mehlich 3 solution. Zinc concentrations in leaves and grains of soybean, wheat, and corn were also determined. Soil pH (0.01molL-1 CaCl2 suspension) varied from 4.4 to 6.1, at the 0- to 0.05-m depth, from 4.2 to 5.3 at the 0.05- to 0.10-m depth, and from 4.2 to 4.8 at the 0.10- to 0.20-m depth, after liming and re-liming. Zinc concentrations evaluated by DTPA-TEA, 0.1molL-1 HCl, Mehlich 1, and Mehlich 3 solutions were not changed as a result of lime rate application. Re-liming increased Zn concentrations extracted by 0.1molL-1 HCl at 0-0.05m deep and by DTPA-TEA at 0.05-0.10m deep. Surface-applied lime promoted a decrease in Zn concentrations of the crops, mainly in grains, because of increased soil pH at the surface layers. Regardless of the liming treatments, levels of Zn were sufficient to soybean, wheat, and corn nutrition under NT.

Leaf area index estimation with MODIS reflectance time series and model inversion during full rotations of Eucalyptus plantations

The leaf area index (LAI) of fast-growing Eucalyptus plantations is highly dynamic both seasonally and interannually, and is spatially variable depending on pedo-climatic conditions. LAI is very important in determining the carbon and water balance of a stand, but is difficult to measure during a complete stand rotation and at large scales. Remote-sensing methods allowing the retrieval of LAI time series with accuracy and precision are therefore necessary. Here, we tested two methods for LAI estimation from MODIS 250m resolution red and near-infrared (NIR) reflectance time series. The first method involved the inversion of a coupled model of leaf reflectance and transmittance (PROSPECT4), soil reflectance (SOILSPECT) and canopy radiative transfer (4SAIL2). Model parameters other than the LAI were either fixed to measured constant values, or allowed to vary seasonally and/or with stand age according to trends observed in field measurements. The LAI was assumed to vary throughout the rotation following a series of alternately increasing and decreasing sigmoid curves. The parameters of each sigmoid curve that allowed the best fit of simulated canopy reflectance to MODIS red and NIR reflectance data were obtained by minimization techniques. The second method was based on a linear relationship between the LAI and values of the GEneralized Soil Adjusted Vegetation Index (GESAVI), which was calibrated using destructive LAI measurements made at two seasons, on Eucalyptus stands of different ages and productivity levels. The ability of each approach to reproduce field-measured LAI values was assessed, and uncertainty on results and parameter sensitivities were examined. Both methods offered a good fit between measured and estimated LAI (R(2) = 0.80 and R(2) = 0.62 for model inversion and GESAVI-based methods, respectively), but the GESAVI-based method overestimated the LAI at young ages. (C) 2010 Elsevier Inc. All rights reserved.

A Log-Linear Regression Model for the Beta-Weibull Distribution

We introduce the log-beta Weibull regression model based on the beta Weibull distribution (Famoye et al., 2005; Lee et al., 2007). We derive expansions for the moment generating function which do not depend on complicated functions. The new regression model represents a parametric family of models that includes as sub-models several widely known regression models that can be applied to censored survival data. We employ a frequentist analysis, a jackknife estimator, and a parametric bootstrap for the parameters of the proposed model. We derive the appropriate matrices for assessing local influences on the parameter estimates under different perturbation schemes and present some ways to assess global influences. Further, for different parameter settings, sample sizes, and censoring percentages, several simulations are performed. In addition, the empirical distribution of some modified residuals are displayed and compared with the standard normal distribution. These studies suggest that the residual analysis usually performed in normal linear regression models can be extended to a modified deviance residual in the proposed regression model applied to censored data. We define martingale and deviance residuals to evaluate the model assumptions. The extended regression model is very useful for the analysis of real data and could give more realistic fits than other special regression models.

MODIS NDVI time-series allow the monitoring of Eucalyptus plantation biomass

The use of remote sensing is necessary for monitoring forest carbon stocks at large scales. Optical remote sensing, although not the most suitable technique for the direct estimation of stand biomass, offers the advantage of providing large temporal and spatial datasets. In particular, information on canopy structure is encompassed in stand reflectance time series. This study focused on the example of Eucalyptus forest plantations, which have recently attracted much attention as a result of their high expansion rate in many tropical countries. Stand scale time-series of Normalized Difference Vegetation Index (NDVI) were obtained from MODIS satellite data after a procedure involving un-mixing and interpolation, on about 15,000 ha of plantations in southern Brazil. The comparison of the planting date of the current rotation (and therefore the age of the stands) estimated from these time series with real values provided by the company showed that the root mean square error was 35.5 days. Age alone explained more than 82% of stand wood volume variability and 87% of stand dominant height variability. Age variables were combined with other variables derived from the NDVI time series and simple bioclimatic data by means of linear (Stepwise) or nonlinear (Random Forest) regressions. The nonlinear regressions gave r-square values of 0.90 for volume and 0.92 for dominant height, and an accuracy of about 25 m(3)/ha for volume (15% of the volume average value) and about 1.6 m for dominant height (8% of the height average value). The improvement including NDVI and bioclimatic data comes from the fact that the cumulative NDVI since planting date integrates the interannual variability of leaf area index (LAI), light interception by the foliage and growth due for example to variations of seasonal water stress. The accuracy of biomass and height predictions was strongly improved by using the NDVI integrated over the two first years after planting, which are critical for stand establishment. These results open perspectives for cost-effective monitoring of biomass at large scales in intensively-managed plantation forests. (C) 2011 Elsevier Inc. All rights reserved.

Relating MODIS vegetation index time-series with structure, light absorption and stem production of fast-growing Eucalyptus plantations

By allowing the estimation of forest structural and biophysical characteristics at different temporal and spatial scales, remote sensing may contribute to our understanding and monitoring of planted forests. Here, we studied 9-year time-series of the Normalized Difference Vegetation Index (NDVI) from the Moderate Resolution Imaging Spectroradiometer (MODIS) on a network of 16 stands in fast-growing Eucalyptus plantations in Sao Paulo State, Brazil. We aimed to examine the relationships between NDVI time-series spanning entire rotations and stand structural characteristics (volume, dominant height, mean annual increment) in these simple forest ecosystems. Our second objective was to examine spatial and temporal variations of light use efficiency for wood production, by comparing time-series of Absorbed Photosynthetically Active Radiation (APAR) with inventory data. Relationships were calibrated between the NDVI and the fractions of intercepted diffuse and direct radiation, using hemispherical photographs taken on the studied stands at two seasons. APAR was calculated from the NDVI time-series using these relationships. Stem volume and dominant height were strongly correlated with summed NDVI values between planting date and inventory date. Stand productivity was correlated with mean NDVI values. APAR during the first 2 years of growth was variable between stands and was well correlated with stem wood production (r(2) = 0.78). In contrast, APAR during the following years was less variable and not significantly correlated with stem biomass increments. Production of wood per unit of absorbed light varied with stand age and with site index. In our study, a better site index was accompanied both by increased APAR during the first 2 years of growth and by higher light use efficiency for stem wood production during the whole rotation. Implications for simple process-based modelling are discussed. (C) 2009 Elsevier B.V. All rights reserved.

Natural regeneration of Atlantic Forest species in the understory of Eucalyptus saligna Smith. in a former forest production unit at the Parque das Neblinas, Bertioga, SP

This study aimed at characterizing the potential for natural regeneration of native vegetation in the under-story of an earlier Eucalyptus saligna Smith production stand. The study was carried out at the Parque das Neblinas, Bertioga municipality, SP, in a 45 ha third rotation stand; which had been abandoned 15 years ago for natural regeneration to occur. The sampling was done in 24 plots of 20 x 40 m. The sampled area was of 19,200 m(2), with inventory made of 100% of the eucalyptus trees. All regeneration trees with a height >= 1.30 m and DBH >= 5.0 cm were measured, as well as adult individuals with DBH >= 5.0 cm; surveyed in two size classes. 1,417 individuals of E. saligna were measured, with a density of 738,02 individuals/ha and a basal area of 22.69 m(2)/ha. Among 2,763 natural regeneration individuals, 111 species belonged to 66 genera and 34 botanical families. The species represented 43.7% of the tree richness of neighboring native forest fragments. The total estimated density and the basal area were respectively 1,052.6 individuals/ha and 6.4 m(2)/ha of autochthonous trees with DBH >= 5.0 cm (Class 1); while for regeneration there were 3,864.58 individuals/ha, and 2.76 m(2)/ha of individuals with a height >= 1.30 m and DBH <5.0 cm (Class 2). Shannon diversity (H`) was 2.83 and 3.68, respectively, for Classes 1 and 2, and the corrected species richness for a 1000-individual sample (R(1000)) were 75.6 and 87.29 (Fisher`s a index) for the same classes. The majority of the species (34.84%) was typical from the understory of wet tropical forest and had zoochoric fruit dispersal (67.57%). The results indicate that, under these conditions, a eucalyptus forest is able to provide adequate regeneration niches for native vegetation, and may represent a sink habitat for local populations.

Production of finger millet (Eleusine coracana (L.) Gaertn.) dry matter under anticipated fertilization of soybeans

Due to the lack of information about the effects of soybean fertilization anticipation and the use of Eleusine coracana (L.) Gaerm. (ANSB Pe-de-galinha 5352) as cover culture, and the need for new techniques for the management of the agro-ecosystem in a more conservationist, sustainable and functional manner the goal of this research was to study the effects of anticipated fertilization on the production of dry matter in E. coracana. The experiments were carried out in an Oxisol, during the growing seasons of 200112002, 200212003 and 200312004, in Piracicaba, SP Brazil (22`50`25""S and 48`01`65""W). The experimental design was of totally randomized blocks and twelve treatments (lev- els of anticipated fertilization) were used with three repetitions. The base fertilization of soybean culture was partially anticipated, in the sowing for the finger-millet crop. Approximately 70 days after sowing, samples of the finger millet plants were collected in order to evaluate the dry matter production and to desiccate them thereafter It is concluded that the anticipation of phosphorus and potassium fertilization of soybean increases the dry matter production of E. coracana; in addition, E. coracana holds a great potential for the production of plant residues and it can be used in culture rotation or in no-tillage production systems.

QSAR Modeling of a Set of Pyrazinoate Esters as Antituberculosis Prodrugs

Tuberculosis is an infection caused mainly by Mycobacterium tuberculosis. A first-line antimycobacterial drug is pyrazinamide (PZA), which acts partially as a prodrug activated by a pyrazinamidase releasing the active agent, pyrazinoic acid (POA). As pyrazinoic acid presents some difficulty to cross the mycobacterial cell wall, and also the pyrazinamide-resistant strains do not express the pyrazinamidase, a set of pyrazinoic acid esters have been evaluated as antimycobacterial agents. In this work, a QSAR approach was applied to a set of forty-three pyrazinoates against M. tuberculosis ATCC 27294, using genetic algorithm function and partial least squares regression (WOLF 5.5 program). The independent variables selected were the Balaban index (I), calculated n-octanol/water partition coefficient (ClogP), van-der-Waals surface area, dipole moment, and stretching-energy contribution. The final QSAR model (N = 32, r(2) = 0.68, q(2) = 0.59, LOF = 0.25, and LSE = 0.19) was fully validated employing leave-N-out cross-validation and y-scrambling techniques. The test set (N = 11) presented an external prediction power of 73%. In conclusion, the QSAR model generated can be used as a valuable tool to optimize the activity of future pyrazinoic acid esters in the designing of new antituberculosis agents.

Topliss method in the optimization of salicylic acid derivatives as potential antimycobacterial agents

The Topliss method was used to guide a synthetic path in support of drug discovery efforts toward the identification of potent antimycobacterial agents. Salicylic acid and its derivatives, p-chloro, p-methoxy, and m-chlorosalicylic acid, exemplify a series of synthetic compounds whose minimum inhibitory concentrations for a strain of Mycobacterium were determined and compared to those of the reference drug, p-aminosalicylic acid. Several physicochemical descriptors (including Hammett`s sigma constant, ionization constant, dipole moment, Hansch constant, calculated partition coefficient, Sterimol-L and -B-4 and molecular volume) were considered to elucidate structure-activity relationships. Molecular electrostatic potential and molecular dipole moment maps were also calculated using the AM1 semi-empirical method. Among the new derivatives, m-chlorosalicylic acid showed the lowest minimum inhibitory concentration. The overall results suggest that both physicochemical properties and electronic features may influence the biological activity of this series of antimycobacterial agents and thus should be considered in designing new p-aminosalicylic acid analogs.