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.

The exponentiated generalized gamma distribution with application to lifetime data

A four-parameter extension of the generalized gamma distribution capable of modelling a bathtub-shaped hazard rate function is defined and studied. The beauty and importance of this distribution lies in its ability to model monotone and non-monotone failure rate functions, which are quite common in lifetime data analysis and reliability. The new distribution has a number of well-known lifetime special sub-models, such as the exponentiated Weibull, exponentiated generalized half-normal, exponentiated gamma and generalized Rayleigh, among others. We derive two infinite sum representations for its moments. We calculate the density of the order statistics and two expansions for their moments. The method of maximum likelihood is used for estimating the model parameters and the observed information matrix is obtained. Finally, a real data set from the medical area is analysed.

An estimate of biometric parameters in eucalyptus clone plantations in Southern Bahia: an application of the airborne laser scanning (ALS) technology

The application of airborne laser scanning (ALS) technologies in forest inventories has shown great potential to improve the efficiency of forest planning activities. Precise estimates, fast assessment and relatively low complexity can explain the good results in terms of efficiency. The evolution of GPS and inertial measurement technologies, as well as the observed lower assessment costs when these technologies are applied to large scale studies, can explain the increasing dissemination of ALS technologies. The observed good quality of results can be expressed by estimates of volumes and basal area with estimated error below the level of 8.4%, depending on the size of sampled area, the quantity of laser pulses per square meter and the number of control plots. This paper analyzes the potential of an ALS assessment to produce certain forest inventory statistics in plantations of cloned Eucalyptus spp with precision equal of superior to conventional methods. The statistics of interest in this case were: volume, basal area, mean height and dominant trees mean height. The ALS flight for data assessment covered two strips of approximately 2 by 20 Km, in which clouds of points were sampled in circular plots with a radius of 13 m. Plots were sampled in different parts of the strips to cover different stand ages. The clouds of points generated by the ALS assessment: overall height mean, standard error, five percentiles (height under which we can find 10%, 30%, 50%,70% and 90% of the ALS points above ground level in the cloud), and density of points above ground level in each percentile were calculated. The ALS statistics were used in regression models to estimate mean diameter, mean height, mean height of dominant trees, basal area and volume. Conventional forest inventory sample plots provided real data. For volume, an exploratory assessment involving different combinations of ALS statistics allowed for the definition of the most promising relationships and fitting tests based on well known forest biometric models. The models based on ALS statistics that produced the best results involved: the 30% percentile to estimate mean diameter (R(2)=0,88 and MQE%=0,0004); the 10% and 90% percentiles to estimate mean height (R(2)=0,94 and MQE%=0,0003); the 90% percentile to estimate dominant height (R(2)=0,96 and MQE%=0,0003); the 10% percentile and mean height of ALS points to estimate basal area (R(2)=0,92 and MQE%=0,0016); and, to estimate volume, age and the 30% and 90% percentiles (R(2)=0,95 MQE%=0,002). Among the tested forest biometric models, the best fits were provided by the modified Schumacher using age and the 90% percentile, modified Clutter using age, mean height of ALS points and the 70% percentile, and modified Buckman using age, mean height of ALS points and the 10% percentile.

How rare is too rare to harvest? Management challenges posed by timber species occurring at low densities in the Brazilian Amazon

Tropical forests are characterized by diverse assemblages of plant and animal species compared to temperate forests. Corollary to this general rule is that most tree species, whether valued for timber or not, occur at low densities (<1 adult tree ha(-1)) or may be locally rare. In the Brazilian Amazon, many of the most highly valued timber species occur at extremely low densities yet are intensively harvested with little regard for impacts on population structures and dynamics. These include big-leaf mahogany (Swietenia macrophylla), ipe (Tabebuia serratifolia and Tabebuia impetiginosa), jatoba (Hymenaea courbaril), and freijo cinza (Cordia goeldiana). Brazilian forest regulations prohibit harvests of species that meet the legal definition of rare - fewer than three trees per 100 ha - but treat all species populations exceeding this density threshold equally. In this paper we simulate logging impacts on a group of timber species occurring at low densities that are widely distributed across eastern and southern Amazonia, based on field data collected at four research sites since 1997, asking: under current Brazilian forest legislation, what are the prospects for second harvests on 30-year cutting cycles given observed population structures, growth, and mortality rates? Ecologically `rare` species constitute majorities in commercial species assemblages in all but one of the seven large-scale inventories we analyzed from sites spanning the Amazon (range 49-100% of total commercial species). Although densities of only six of 37 study species populations met the Brazilian legal definition of a rare species, timber stocks of five of the six timber species declined substantially at all sites between first and second harvests in simulations based on legally allowable harvest intensities. Reducing species-level harvest intensity by increasing minimum felling diameters or increasing seed tree retention levels improved prospects for second harvests of those populations with a relatively high proportion of submerchantable stems, but did not dramatically improve projections for populations with relatively flat diameter distributions. We argue that restrictions on logging very low-density timber tree populations, such as the current Brazilian standard, provide inadequate minimum protection for vulnerable species. Population declines, even if reduced-impact logging (RIL) is eventually adopted uniformly, can be anticipated for a large pool of high-value timber species unless harvest intensities are adapted to timber species population ecology, and silvicultural treatments are adopted to remedy poor natural stocking in logged stands. (C) 2008 Elsevier B.V. All rights reserved.

What loggers leave behind: Impacts on big-leaf mahogany (Swietenia macrophylla) commercial populations and potential for post-logging recovery in the Brazilian Amazon

The sustainability of current harvest practices for high-value Meliaceae can be assessed by quantifying logging intensity and projecting growth and survival by post-logging populations over anticipated intervals between harvests. From 100%-area inventories of big-leaf mahogany (Swietenia macrophylla) covering 204 ha or more at eight logged and unlogged forest sites across southern Brazilian Amazonia, we report generally higher landscape-scale densities and smaller population-level mean diameters in eastern forests compared to western forests, where most commercial stocks survive. Density of trees >= 20 cm diameter varied by two orders of magnitude and peaked at 1.17 ha(-1). Size class frequency distributions appeared unimodal at two high-density sites, but were essentially arnodal or flat elsewhere; diameter increment patterns indicate that populations were multi- or all-aged. At two high-density sites, conventional logging removed 93-95% of commercial trees (>= 45 cm diameter at the time of logging), illegally eliminated 31-47% of sub-merchantable trees, and targeted trees as small as 20 cm diameter. Projected recovery by commercial stems during 30 years after conventional logging represented 9.9-37.5% of initial densities and was highly dependent on initial logging intensity and size class frequency distributions of commercial trees. We simulated post-logging recovery over the same period at all sites according to the 2003 regulatory framework for mahogany in Brazil, which raised the minimum diameter cutting limit to 60 cm and requires retention during the first harvest of 20% of commercial-sized trees. Recovery during 30 years ranged from approximately 0 to 31% over 20% retention densities at seven of eight sites. At only one site where sub-merchantable trees dominated the population did the simulated density of harvestable stems after 30 years exceed initial commercial densities. These results indicate that 80% harvest intensity will not be sustainable over multiple cutting cycles for most populations without silvicultural interventions ensuring establishment and long-term growth of artificial regeneration to augment depleted natural stocks, including repeated tending of outplanted seedlings. Without improved harvest protocols for mahogany in Brazil as explored in this paper, future commercial supplies of this species as well as other high-value tropical timbers are endangered. Rapid changes in the timber industry and land-use in the Amazon are also significant challenges to sustainable management of mahogany. (C) 2007 Elsevier B.V. All rights reserved.

Variation of wood density and relationship with the tree-ring width of mahogany trees, Swietenia macrophylla, from Amazonian tropical forest of Peru

Mahogany trees, Swietenia macrophylla, occur in open rainforest, semi deciduous and deciduous and dense rainforest of Peruvian Amazonian tropical forest. They occur, preferentially, in areas with a defined dry season, with typical phenology and seasonal variation activity, forming distinct tree-rings. The present work had as aim to determine the wood density radial variation of 14 mahogany trees, of two populations of the Peruvian Amazonian tropical forest, through the X-ray densitometry and to evaluate their application as methodology, compared to the classic method of measurement table, for the determination of the treering width. The radial wood apparent density of the trees profiles rendered it possible to delimit the areas of juvenile-adult wood and of the heartwood-sapwood, relative to the anatomical structure and chemical composition differences, due to the extractives and the vessels obstruction by tyloses. The mean, minimum and maximum wood apparent density of the mahogany trees for the Populations A and B were of 0.70; 0.29; 1.01 g.cm(-3) and 0.81; 0.29; 1.19 g.cm(-3), respectively. The analysis of the variance and mean test indicate differences of mean wood density among the mahogany trees of each population, probably due to the age of the trees. There was no correlation between mean wood density of mahogany trees among the two populations, as well as, between the tree-ring width and the respective mean density. The X-ray densitometry technique is an important tool in the evaluation of the radial variation of wood apparent density and the delimitation of tree-ring boundaries, with correlations of 0.94 and 0.93 in relation to measurement table, for each sampled population.

Distribution of Aceria guerreronis and Neoseiulus baraki among and within coconut bunches in northeast Brazil

Aceria guerreronis Keifer (Acari: Eriophyidae) is considered a major pest of coconut in many countries in the Americas, Africa and parts of Asia. Neoseiulus baraki Athias-Henriot (Acari: Phytoseiidae) is one of the predatory mites most commonly found in association with A. guerreronis in parts of northeast Brazil. The objective of this work was to study the distribution of A. guerreronis and N. baraki among and within coconut bunches. The hypothesis was tested that A. guerreronis and N. baraki are homogenously distributed over the fruits in a bunch, independent of the fruits` age and position. Five collections of bunches, each corresponding to leaves 12-16 from apex (about 2-6 month-old), were conducted in each of three fields in northeastern Brazil, from February to October, 2007. A total of 1,986 fruits were examined. The number of mites, the percentage of fruits hosting them and the level of damage caused by A. guerreronis were evaluated. The highest density of A. guerreronis was observed on fruits of bunch 4 whereas the highest density of N. baraki was observed on bunch 5. Considering all fruits together, no significant differences were observed between densities of either A. guerreronis or N. baraki among the basal, median and apical thirds of the bunches. In younger bunches, fruits of the apical region tend to have lower densities of both mites than fruits of the basal region. This pattern, in association with a similar pattern for the percentage of fruits hosting N. baraki, suggests that the predator initially reaches the basal bunch region, from where it moves to the apical region. The results of the present study suggest that the pest population reduction in bunches older than bunch 4 could be due to (1) an effect of the predator, (2) reduction of the proportion of undamaged tissues amenable to attack, and/or (3) less favorable characteristics of the fruits to attack by A. guerreronis, as indicated by their increasing lignin content as they get older.

Root Water Extraction under Combined Water and Osmotic Stress

Using a numerical implicit model for root water extraction by a single root in a symmetric radial flow problem, based on the Richards equation and the combined convection-dispersion equation, we investigated some aspects of the response of root water uptake to combined water and osmotic stress. The model implicitly incorporates the effect of simultaneous pressure head and osmotic head on root water uptake, and does not require additional assumptions (additive or multiplicative) to derive the combined effect of water and salt stress. Simulation results showed that relative transpiration equals relative matric flux potential, which is defined as the matric flux potential calculated with an osmotic pressure head-dependent lower bound of integration, divided by the matric flux potential at the onset of limiting hydraulic conditions. In the falling rate phase, the osmotic head near the root surface was shown to increase in time due to decreasing root water extraction rates, causing a more gradual decline of relative transpiration than with water stress alone. Results furthermore show that osmotic stress effects on uptake depend on pressure head or water content, allowing a refinement of the approach in which fixed reduction factors based on the electrical conductivity of the saturated soil solution extract are used. One of the consequences is that osmotic stress is predicted to occur in situations not predicted by the saturation extract analysis approach. It is also shown that this way of combining salinity and water as stressors yields results that are different from a purely multiplicative approach. An analytical steady state solution is presented to calculate the solute content at the root surface, and compared with the outputs of the numerical model. Using the analytical solution, a method has been developed to estimate relative transpiration as a function of system parameters, which are often already used in vadose zone models: potential transpiration rate, root length density, minimum root surface pressure head, and soil theta-h and K-h functions.

Optimal irrigation depth for bean crops considering water constraints and farmer`s risk aversion level

Nowadays, the rising competition for the use of water and environmental resources with consequent restrictions for farmers should change the paradigm in terms of irrigation concepts, or rather, in order to attain economical efficiency other than to supply water requirement for the crop. Therefore, taking into account the social and economical role of bean activity in Brazil, as well as the risk inherent to crop due to its high sensibility to both deficit and excessive water, the optimization methods regarding to irrigation management have become more interesting and essential. This study intends to present a way to determine the optimal water supply, considering different combinations between desired bean yield and level of risk, bringing as a result a graph with the former associated with the latter, depending on different water depths.

Pricing farm-level agricultural insurance: a Bayesian approach

This paper applies Hierarchical Bayesian Models to price farm-level yield insurance contracts. This methodology considers the temporal effect, the spatial dependence and spatio-temporal models. One of the major advantages of this framework is that an estimate of the premium rate is obtained directly from the posterior distribution. These methods were applied to a farm-level data set of soybean in the State of the Parana (Brazil), for the period between 1994 and 2003. The model selection was based on a posterior predictive criterion. This study improves considerably the estimation of the fair premium rates considering the small number of observations.

Crop rotation in center-pivot for phytonematode control: density variation, pathogenicity and crop loss estimation

Crop rotation in center-pivot for phytonematode control: density variation, pathogenicity and crop loss estimation A field study conducted over three consecutive years, on a farm using crop rotation system under center-pivot and infested with the nematodes Pratylenchus brachyurus, P. zeae, Meloidogyne incognita, Paratrichodorus minor, Helicotylenchus dihystera, Mesocriconema ornata and M. onoense, demonstrated that intensive crop systems provide conditions for the maintenance of high densities of polyphagous phytonematodes. Of the crops established on the farm (cotton, maize, soybean and cowpea), cotton and soybean suffered the most severe crop losses, caused respectively by M. incognita and P. brachyurus. Since maize is a good host for both nematodes, but tolerant of M. incognita, its exclusion from cropping system would be favorable to the performance of cotton, soybean and cowpea. Results from experiments carried out in controlled conditions confirmed the pathogenicity of P. brachyurus on cotton. Additional management with genetic resistance was useful in fields infested with M. incognita, although the soybean performance was affected by low resistance of the cultivars used for P. brachyurus. In conclusion, crop rotation must be carefully planned in areas infested with polyphagous nematodes, specifically in the case of occurrence of two or more major pathogenic nematodes.

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.

Soil density evaluated by spectral reflectance as an evidence of compaction effects

Soil compaction, reflected by high bulk density, is an environmental degradation process and new technologies are being developed for its detection. Despite the proven efficiency of remote sensing, it has not been widely used for soil density. Our objective was to evaluate the density of two soils: a Typic Quartzpisament (TQ) and a Rhodic Paleudalf (RP), using spectral reflectance obtained by a laboratory spectroradiometer between 450 and 2500 nm. Undisturbed samples were taken at two depths (0-20 and 60-80 cm), and were artificially compacted. Spectral data, obtained before and after compaction, were compared for both wet and dried compacted samples. Results demonstrated that soil density was greater in RP than in TQ at both depths due to its clayey texture. Spectral data detected high density (compacted) from low density (non-compacted) clayey soils under both wet and dry conditions. The detection of density in sandy soils by spectral reflectance was not possible. The intensity of spectral reflectance of high soil bulk density (compacted) samples was higher than for low density (non-compacted) soils due to changes in soil structure and porosity. Dry samples with high bulk density showed differences in the spectral intensity, but not in the absorption features. Wet samples in equal condition had statistically higher reflectance intensity than that of the low soil bulk density (non-compacted), and absorption differences at 1920 nm, which was due to the altered position of the water molecules. Soil line and spectral reflectance used together could detect soil bulk density variations for the clay soil. This technique could assist in the detection of high soil density in the laboratory by providing new soil information.

Evaluation of Automotive Varnish as a Coating for Density Measurements of Soil Clods

Saran F-310 resin (Dow Chemical Co, Midland, MI) has been widely used to coat soil clods for density and size measurements; however, the manufacturer has recently stopped producing this resin and supplies are difficult to obtain. Hence, we evaluated the feasibility of using Lazzudur 7502 (Sherwin-Williams, Cleveland, OH) automotive varnish to coat soil clods for density measurement. Preliminary evaluations showed that immersion of clods in the varnish did nor affect clod cohesion and that a single immersion in Lazzudur with 30 min of post-immersion drying produced density results nor significantly (P < 0.05) different to those obtained using saran. This technique was tested across seven soils and no significant (P < 0.05) difference was found in the density of the clods measured using the two coating methods. This work suggests that automotive varnish can he used as an alternative to saran resin for clod density measurements.

Parametric correlation functions to model the structure of permanent environmental (co)variances in milk yield random regression models

The objective of the present study was to estimate milk yield genetic parameters applying random regression models and parametric correlation functions combined with a variance function to model animal permanent environmental effects. A total of 152,145 test-day milk yields from 7,317 first lactations of Holstein cows belonging to herds located in the southeastern region of Brazil were analyzed. Test-day milk yields were divided into 44 weekly classes of days in milk. Contemporary groups were defined by herd-test-day comprising a total of 2,539 classes. The model included direct additive genetic, permanent environmental, and residual random effects. The following fixed effects were considered: contemporary group, age of cow at calving (linear and quadratic regressions), and the population average lactation curve modeled by fourth-order orthogonal Legendre polynomial. Additive genetic effects were modeled by random regression on orthogonal Legendre polynomials of days in milk, whereas permanent environmental effects were estimated using a stationary or nonstationary parametric correlation function combined with a variance function of different orders. The structure of residual variances was modeled using a step function containing 6 variance classes. The genetic parameter estimates obtained with the model using a stationary correlation function associated with a variance function to model permanent environmental effects were similar to those obtained with models employing orthogonal Legendre polynomials for the same effect. A model using a sixth-order polynomial for additive effects and a stationary parametric correlation function associated with a seventh-order variance function to model permanent environmental effects would be sufficient for data fitting.