Carbon and nitrogen cycling in a tropical Brazilian soil cropped with sugarcane and irrigated with wastewater

Carbon (C) and nitrogen (N) dynamics in agro-systems can be altered as a consequence of treated sewage effluent (TSE) irrigation. The present study evaluated the effects of TSE irrigation over 16 months on N concentrations in sugarcane (leaves, stalks and juice), total soil carbon (TC), total soil nitrogen (TN), NO(3)(-)-N in soil and nitrate (NO(3)(-)) and dissolved organic carbon (DOC) in soil solution. The soil was classified as an Oxisol and samplings were carried out during the first productive crop cycle, from February 2005 (before planting) to September 2006 (after sugarcane harvest and 16 months of TSE irrigation). The experiment was arranged in a complete block design with five treatments and four replicates. Irrigated plots received 50% of the recommended mineral N fertilization and 100% (T100), 125% (T125), 150% (T150) and 200% (T200) of crop water demand. No mineral N and irrigation were applied to the control plots. TSE irrigation enhanced sugarcane yield but resulted in total-N inputs(804-1622 kg N ha(-1)) greater than exported N (463-597 kg N ha(-1)). Hence, throughout the irrigation period, high NO(3)(-) concentrations (up to 388 mg L(-1) at T200) and DOC (up to 142 mg L(-1) at T100) were measured in soil solution below the root zone, indicating the potential of groundwater contamination. TSE irrigation did not change soil TC and TN. (C) 2009 Elsevier B.V. All rights reserved.

Time stability of soil water storage measured by neutron probe and the effects of calibration procedures in a small watershed

The knowledge of soil water storage (SWS) of soil profiles is crucial for the adoption of vegetation restoration practices. With the aim of identifying representative sites to obtain the mean SWS of a watershed, a time stability analysis of neutron probe evaluations of SWS was performed by the means of relative differences and Spearman rank correlation coefficients. At the same time, the effects of different neutron probe calibration procedures were explored on time stability analysis. mean SWS estimation. and preservation of the spatial variability of SWS. The selected watershed, with deep gullies and undulating slopes which cover an area of 20 ha, is characterized by an Ust-Sandiic Entisol and an Aeolian sandy soil. The dominant vegetation species are bunge needlegrass (Stipa bungeana Trim) and korshinsk peashrub (Carugano Korshinskii kom.). From June 11, 2007 to July 23,2008, SWS of the top1 m soil layer was evaluated for 20 dates, based on neutron probe data of 12 sampling sites. Three calibration procedures were employed: type 1, most complete, with each site having its own linear calibration equation (TrE); type II. with TrE equations extended over the whole field: and type III, with one single linear calibration curve for the whole field (UnE) and also correcting its intercept based on site specific relative difference analysis (RdE) and on linear fitting of data (RcE), both maintaining the same slope. A strong time stability of SWS estimated by TrE equations was identified. Soil particle size and soil organic matter content were recognized as the influencing factors for spatial variability of SWS. Land use influenced neither the spatial variability nor the time stability of SWS. Time stability analysis identified one site to represent the mean SWS of the whole watershed with mean absolute percentage errors of less than 10%, therefore. this site can be used as a predictor for the mean SWS of the watershed. Some equations of type II were found to be unsatisfactory to yield reliable mean SWS values or in preserving the associated soil spatial variability. Hence, it is recommended to be cautious in extending calibration equations to other sites since they might not consider the field variability. For the equations with corrected intercept (type III), which consider the spatial variability of calibration in a different way in relation to TrE, it was found that they can yield satisfactory means and standard deviation of SWS, except for the RdE equations, which largely leveled off the SWS values in the watershed. Correlation analysis showed that the neutron probe calibration was linked to soil bulk density and to organic matter content. Therefore, spatial variability of soil properties should be taken into account during the process of neutron probe calibration. This study provides useful information on the mean SWS observation with a time stable site and on distinct neutron probe calibration procedures, and it should be extended to soil water management studies with neutron probes, e.g., the process of vegetation restoration in wider area and soil types of the Loess Plateau in China. (C) 2009 Elsevier B.V. All rights reserved.

A log-extended Weibull regression model

A bathtub-shaped failure rate function is very useful in survival analysis and reliability studies. The well-known lifetime distributions do not have this property. For the first time, we propose a location-scale regression model based on the logarithm of an extended Weibull distribution which has the ability to deal with bathtub-shaped failure rate functions. We use the method of maximum likelihood to estimate the model parameters and some inferential procedures are presented. We reanalyze a real data set under the new model and the log-modified Weibull regression model. We perform a model check based on martingale-type residuals and generated envelopes and the statistics AIC and BIC to select appropriate models. (C) 2009 Elsevier B.V. All rights reserved.

Bias-corrected Pearson estimating functions for Taylor`s power law applied to benthic macrofauna data

Estimation of Taylor`s power law for species abundance data may be performed by linear regression of the log empirical variances on the log means, but this method suffers from a problem of bias for sparse data. We show that the bias may be reduced by using a bias-corrected Pearson estimating function. Furthermore, we investigate a more general regression model allowing for site-specific covariates. This method may be efficiently implemented using a Newton scoring algorithm, with standard errors calculated from the inverse Godambe information matrix. The method is applied to a set of biomass data for benthic macrofauna from two Danish estuaries. (C) 2011 Elsevier B.V. All rights reserved.

A generalized modified Weibull distribution for lifetime modeling

A four parameter generalization of the Weibull distribution capable of modeling a bathtub-shaped hazard rate function is defined and studied. The beauty and importance of this distribution lies in its ability to model monotone as well as non-monotone failure rates, which are quite common in lifetime problems and reliability. The new distribution has a number of well-known lifetime special sub-models, such as the Weibull, extreme value, exponentiated Weibull, generalized Rayleigh and modified Weibull distributions, among others. We derive two infinite sum representations for its moments. The density of the order statistics is obtained. The method of maximum likelihood is used for estimating the model parameters. Also, the observed information matrix is obtained. Two applications are presented to illustrate the proposed distribution. (C) 2008 Elsevier B.V. All rights reserved.

On the use of correlated beta random variables with animal population modelling

We give reasons why demographic parameters such as survival and reproduction rates are often modelled well in stochastic population simulation using beta distributions. In practice, it is frequently expected that these parameters will be correlated, for example with survival rates for all age classes tending to be high or low in the same year. We therefore discuss a method for producing correlated beta random variables by transforming correlated normal random variables, and show how it can be applied in practice by means of a simple example. We also note how the same approach can be used to produce correlated uniform triangular, and exponential random variables. (C) 2008 Elsevier B.V. All rights reserved.

Log-modified Weibull regression models with censored data: Sensitivity and residual analysis

This paper proposes a regression model considering the modified Weibull distribution. This distribution can be used to model bathtub-shaped failure rate functions. Assuming censored data, we consider maximum likelihood and Jackknife estimators for the parameters of the model. We derive the appropriate matrices for assessing local influence on the parameter estimates under different perturbation schemes and we also present some ways to perform global influence. Besides, for different parameter settings, sample sizes and censoring percentages, various simulations are performed and the empirical distribution of the modified deviance residual is displayed and compared with the standard normal distribution. These studies suggest that the residual analysis usually performed in normal linear regression models can be straightforwardly extended for a martingale-type residual in log-modified Weibull regression models with censored data. Finally, we analyze a real data set under log-modified Weibull regression models. A diagnostic analysis and a model checking based on the modified deviance residual are performed to select appropriate models. (c) 2008 Elsevier B.V. All rights reserved.

Spatio-temporal variability behavior of land surface soil water content in shrub- and grass-land

The knowledge of the relationship between spatial variability of the surface soil water content (theta) and its mean across a spatial domain (theta(m)) is crucial for hydrological modeling and understanding soil water dynamics at different scales. With the aim to compare the soil moisture dynamics and variability between the two land uses and to explore the relationship between the spatial variability of theta and theta(m), this study analyzed sets of surface theta measurements performed with an impedance soil moisture probe, collected 136 times during a period of one year in two transects covering different land uses, i.e., korshinsk peashrub transect (KPT) and bunge needlegrass transect (BNT), in a watershed of the Loess Plateau, China. Results showed that the temporal pattern of theta behaved similarly for the two land uses, with both relative wetter soils during wet period and relative drier soils during dry period recognized in BNT. Soil moisture tended to be temporally stable among different dates, and more stable patterns could be observed for dates with more similar soil water conditions. The magnitude of the spatial variation of theta in KPT was greater than that in ENT. For both land uses, the standard deviation (SD) of theta in general increased as theta(m) increased, a behavior that could be well described with a natural logarithmic function. Convex relationship of CV and theta(m) and the maximum CV for both land uses (43.5% in KPT and 41.0% in BNT) can, therefore, be ascertained. Geostatistical analysis showed that the range in KPT (9.1 m) was shorter than that in BNT (15.1 m). The nugget effects, the structured variability, hence the total variability increased as theta(m) increased. For both land uses, the spatial dependency in general increased with increasing theta(m). 2011 Elsevier B.V. All rights reserved.

On estimation and influence diagnostics for zero-inflated negative binomial regression models

The zero-inflated negative binomial model is used to account for overdispersion detected in data that are initially analyzed under the zero-Inflated Poisson model A frequentist analysis a jackknife estimator and a non-parametric bootstrap for parameter estimation of zero-inflated negative binomial regression models are considered In addition an EM-type algorithm is developed for performing maximum likelihood estimation Then the appropriate matrices for assessing local influence on the parameter estimates under different perturbation schemes and some ways to perform global influence analysis are derived In order to study departures from the error assumption as well as the presence of outliers residual analysis based on the standardized Pearson residuals is discussed The relevance of the approach is illustrated with a real data set where It is shown that zero-inflated negative binomial regression models seems to fit the data better than the Poisson counterpart (C) 2010 Elsevier B V All rights reserved

Regression models for grouped survival data: Estimation and sensitivity analysis

In this study, regression models are evaluated for grouped survival data when the effect of censoring time is considered in the model and the regression structure is modeled through four link functions. The methodology for grouped survival data is based on life tables, and the times are grouped in k intervals so that ties are eliminated. Thus, the data modeling is performed by considering the discrete models of lifetime regression. The model parameters are estimated by using the maximum likelihood and jackknife methods. To detect influential observations in the proposed models, diagnostic measures based on case deletion, which are denominated global influence, and influence measures based on small perturbations in the data or in the model, referred to as local influence, are used. In addition to those measures, the local influence and the total influential estimate are also employed. Various simulation studies are performed and compared to the performance of the four link functions of the regression models for grouped survival data for different parameter settings, sample sizes and numbers of intervals. Finally, a data set is analyzed by using the proposed regression models. (C) 2010 Elsevier B.V. All rights reserved.

Can Eucalyptus invade native forest fragments close to commercial stands?

Some Eucalyptus species are widely used as a plantation crop in tropical and subtropical regions. One reason for this is the diversity of end uses, but the main reason is the high level of wood production obtained from commercial plantings. With the advancement of biotechnology it will be possible to expand the geographical area in which eucalypts can be used as commercial plantation crops, especially in regions with current climatic restrictions. Despite the popularity of eucalypts and their increasing range, questions still exist, in both traditional planting areas and in the new regions: Can eucalypts invade areas of native vegetation, causing damage to natural ecosystems biodiversity? The objective of this study it was to assess whether eucalypts can invade native vegetation fragments in proximity to commercial stands, and what factors promote this invasive growth. Thus, three experiments were established in forest fragments located in three different regions of Brazil. Each experiment was composed of 40 plots (1 m(2) each one), 20 plots located at the border between the forest fragment and eucalypts plantation, and 20 plots in the interior of the forest fragments. In each experimental site, the plots were paired by two soil exposure conditions, 10 plots in natural conditions and 10 plots with soil exposure (no plant and no litter). During the rainy season, 2 g of eucalypts seeds were sown in each plot, including Eucalyptus grandis or a hybrid of E. urophylla x E. grandis, the most common commercial eucalypt species planted in the three region. At 15, 30, 45, 90, 180, 270 and 360 days after sowing, we assessed the number of seedlings of eucalypts and the number of seedlings of native species resulting from natural regeneration. Fifteen days after sowing, the greatest number of eucalypts seedlings (37 m(-2)) was observed in the plots with lower luminosity and exposed soil. Also, for native species, it was observed that exposed soil improved natural germination reaching the highest number of 163 seedlings per square meter. Site and soil exposure were the factors that have the greatest influence on seed germination of both eucalypt and native species. However, 270 days after sowing, eucalypt seedlings were not observed at any of the three experimental sites. The result shows the inability of eucalypts to adapt to condition outside of their natural range. However, native species demonstrated their strong capacity for natural regeneration in forest fragments under the same conditions where eucalypts were seeded. (C) 2011 Elsevier B.V. All rights reserved.

Biogeochemical cycles of nutrients in tropical Eucalyptus plantations Main features shown by intensive monitoring in Congo and Brazil

The sustainability of fast-growing tropical Eucalyptus plantations is of concern in a context of rising fertilizer costs, since large amounts of nutrients are removed with biomass every 6-7 years from highly weathered soils. A better understanding of the dynamics of tree requirements is required to match fertilization regimes to the availability of each nutrient in the soil. The nutrition of Eucalyptus plantations has been intensively investigated and many studies have focused on specific fluxes in the biogeochemical cycles of nutrients. However, studies dealing with complete cycles are scarce for the Tropics. The objective of this paper was to compare these cycles for Eucalyptus plantations in Congo and Brazil, with contrasting climates, soil properties, and management practices. The main features were similar in the two situations. Most nutrient fluxes were driven by crown establishment the two first years after planting and total biomass production thereafter. These forests were characterized by huge nutrient requirements: 155, 10, 52, 55 and 23 kg ha(-1) of N, P, K, Ca and Mg the first year after planting at the Brazilian study site, respectively. High growth rates the first months after planting were essential to take advantage of the large amounts of nutrients released into the soil solutions by organic matter mineralization after harvesting. This study highlighted the predominant role of biological and biochemical cycles over the geochemical cycle of nutrients in tropical Eucalyptus plantations and indicated the prime importance of carefully managing organic matter in these soils. Limited nutrient losses through deep drainage after clear-cutting in the sandy soils of the two study sites showed the remarkable efficiency of Eucalyptus trees in keeping limited nutrient pools within the ecosystem, even after major disturbances. Nutrient input-output budgets suggested that Eucalyptus plantations take advantage of soil fertility inherited from previous land uses and that long-term sustainability will require an increase in the inputs of certain nutrients. (C) 2009 Elsevier B.V. All rights reserved.

A positive growth response to NaCl applications in Eucalyptus plantations established on K-deficient soils

Contrasting responses of Eucalyptus trees to K fertilizer applications have been reported on soils with low K contents. A complete randomized block experiment was set up in Brazil to test the hypothesis that large atmospheric deposits of NaCl in coastal regions might lead to a partial substitution of K by Na in Eucalyptus physiology and enhance tree growth. Treatments with application of 1.5, 3.0, 4.5 kmol K ha(-1) (K(1.5), K(3.0), 1(4.5, respectively) as KCl, 3.0 kmol K ha(-1) applied as K(2)SO(4), 3.0 kmol Na ha(-1) (Na(3.0)) as NaCl commercialized for cattle feeding, and a mixture of 1.5 kmol K + 1.5 kmol Na ha(-1) (K(1.5) + Na(1.5)) were compared to a control treatment (C) with no K and Na applications. All the plots were fertilized with large amounts of the other nutrients. A positive effect of NaCl applications on the growth of E. grandis trees was observed. NaCl and KCl additions in treatments Na(3.0) and K(3.0) increased above-ground biomass by 56% and 130% three years after planting, respectively, in comparison with the C treatment. By contrast, accumulated litterfall up to age 3 years was not significantly modified. NaCl applications in the Na(3.0) treatment significantly increased Na accumulation in above-ground tree components but did not modify K accumulation, whatever the sampling age. A partial substitution of K by Na in tree physiology, as observed for various agricultural crops, might explain this behaviour. Our results suggest the possibility of applying inexpensive K fertilizers, which are less purified in Na, and explain why high yields are achieved without K fertilizer applications in areas with large dry depositions of marine aerosols. Further investigations are necessary to identify the processes involving Na in Eucalyptus tree physiology. (C) 2009 Elsevier B.V. All rights reserved.

Can intensive management accelerate the restoration of Brazil`s Atlantic forests?

Only 7% of the once extensive forest along the eastern coast of Brazil remains, and much of that is degraded and threatened by agricultural expansion and urbanization. We wondered if methods similar to those developed to establish fast-growing Eucalyptus plantations might also work to enhance survival and growth of rainforest species on degraded pastures composed of highly competitive C(4) grasses. An 8-factor experiment was laid out to contrast the value of different intensities of cultivation, application of fertilizer and weed control on the growth and survival of a mixture of 20 rainforest species planted at two densities: 3 m x 1 m, and 3 m x 2 m. Intensive management increased seedling survival from 90% to 98%, stemwood production and leaf area index (LAI) by similar to 4-fold, and stemwood production per unit of light absorbed by 30%. Annual growth in stem biomass was closely related to LAI alone (r(2) = 0.93, p < 0.0001), and the regression improved further in combination with canopy nitrogen content (r(2) =0.99, p < 0.0001). Intensive management resulted in a nearly closed forest canopy in less than 4 years, and offers a practical means to establish functional forests on abandoned agricultural land. (C) 2009 Elsevier B.V. All rights reserved.

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.