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.

Radial variation of the wood anatomical structure of Gmelina arborea trees from different climatic and management conditions in Costa Rica

The evaluations of the effect of the climatic conditions and of the intensity of forest management in the trunk of the Gmelina arborea Linn. Roxb. trees are restricted to its physical-mechanical properties and use. The present work has as objective to study the radial variations of the wood anatomy of the gmelina trees sampled in plantations of 30 sites in Costa Rica, characterized by two climatic conditions (tropical dry and humid) and three intensities of forest management (intensive, moderate and without management). The results of the analyses demonstrated the existence of radial variation of the different anatomical parameters, except for the fiber lumen diameter and multiple vessels in the wood of the gmelina trees. For the wood anatomical elements, fibers (width, lumen diameter, and length), vessels (multiple vessels, diameter and frequency) and radial parenchyma (height) relationships were observed with the climate (tropical humid and dry). The radial variations of the wood anatomical elements were, also, influenced by the management regimes of the gmelina trees.

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.

Distribution of Tetranychus evansi and its predator Phytoseiulus longipes (Acari : Tetranychidae, Phytoseiidae) in southern Brazil

Tetranychus evansi is an important pest of tomato in several countries. The predatory mite Phytoseiulus longipes has been found in association with it in Uruguaiana, State of Rio Grande do Sul, Brazil. The objective of this study was to evaluate the distribution of those two species in Uruguaiana region. Bi-weekly samples of plant parts were taken between January 28 and April 30, 2007 in Uruguaiana and five neighboring counties. Tetranychus evansi was found in all counties, but P. longipes was only found in Uruguaiana, mainly in the urban area. It is conceivable that such restricted distribution of P. longipes is due to its recent introduction to the region. Alternatively, unfavorable winter conditions could prevent its persistence, except in protected places. Buildings and vegetation characteristics of the urban area of Uruguaiana could moderate strong winds, intense rainfall and low temperatures that commonly occur in the area.

Evaluation of FAO Penman-Monteith and alternative methods for estimating reference evapotranspiration with missing data in Southern Ontario, Canada

Grass reference evapotranspiration (ETo) is an important agrometeorological parameter for climatological and hydrological studies, as well as for irrigation planning and management. There are several methods to estimate ETo, but their performance in different environments is diverse, since all of them have some empirical background. The FAO Penman-Monteith (FAD PM) method has been considered as a universal standard to estimate ETo for more than a decade. This method considers many parameters related to the evapotranspiration process: net radiation (Rn), air temperature (7), vapor pressure deficit (Delta e), and wind speed (U); and has presented very good results when compared to data from lysimeters Populated with short grass or alfalfa. In some conditions, the use of the FAO PM method is restricted by the lack of input variables. In these cases, when data are missing, the option is to calculate ETo by the FAD PM method using estimated input variables, as recommended by FAD Irrigation and Drainage Paper 56. Based on that, the objective of this study was to evaluate the performance of the FAO PM method to estimate ETo when Rn, Delta e, and U data are missing, in Southern Ontario, Canada. Other alternative methods were also tested for the region: Priestley-Taylor, Hargreaves, and Thornthwaite. Data from 12 locations across Southern Ontario, Canada, were used to compare ETo estimated by the FAD PM method with a complete data set and with missing data. The alternative ETo equations were also tested and calibrated for each location. When relative humidity (RH) and U data were missing, the FAD PM method was still a very good option for estimating ETo for Southern Ontario, with RMSE smaller than 0.53 mm day(-1). For these cases, U data were replaced by the normal values for the region and Delta e was estimated from temperature data. The Priestley-Taylor method was also a good option for estimating ETo when U and Delta e data were missing, mainly when calibrated locally (RMSE = 0.40 mm day(-1)). When Rn was missing, the FAD PM method was not good enough for estimating ETo, with RMSE increasing to 0.79 mm day(-1). When only T data were available, adjusted Hargreaves and modified Thornthwaite methods were better options to estimate ETo than the FAO) PM method, since RMSEs from these methods, respectively 0.79 and 0.83 mm day(-1), were significantly smaller than that obtained by FAO PM (RMSE = 1.12 mm day(-1). (C) 2009 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.

Yield loss caused by huanglongbing in different sweet orange cultivars in So Paulo, Brazil

Citrus huanglongbing (HLB) reduces an affected orchard`s economic life. This work aimed to characterize yield loss due to HLB for different sweet orange cultivars and determine the relationship between disease severity and yield. Disease severity and yield were assessed on 949 individual trees distributed in 11 different blocks from sweet orange cultivars Hamlin, Westin, Pera and Valencia. In each block, plants showing a range of HLB severity levels and asymptomatic plants were selected. Total yield (weight of harvested fruit), mean weight of asymptomatic and symptomatic fruit, relative yield (symptomatic tree yield/mean yield of asymptomatic trees from the same block) and relative number of fruits (fruit number from symptomatic tree/mean number of fruits from asymptomatic trees from the same block) were determined. The weight of symptomatic fruit was lower than the weight of asymptomatic fruit, but the weights of asymptomatic and symptomatic fruit were not correlated with disease severity, indicating that the effects of HLB were restricted to symptomatic branches. The relationship of the relative yield with HLB severity can be satisfactorily described by a negative exponential model. The rates of yield decrease as a function of disease severity were similar for all assessed cultivars. A relative yield (up to 19%) was observed even for trees where disease severity was 100%. The strong linear relationship between relative number of fruits per tree and the relative yield per tree suggested that the yield reduction was due primarily to early fruit drop or lack of fruit set on affected branches.

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.

An integrated molecular map of yellow passion fruit based on simultaneous maximum-likelihood estimation of linkage and linkage phases

The development of genetic maps for auto-incompatible species, such as the yellow passion fruit (Passiflora edulis Sims f.flavicarpa Deg.) is restricted due to the unfeasibility of obtaining traditional mapping populations based on inbred lines. For this reason, yellow passion fruit linkage maps were generally constructed using a strategy known as two-way pseudo-testeross, based on monoparental dominant markers segregating in a 1:1 fashion. Due to the lack of information from these markers in one of the parents, two individual (parental) maps were obtained. However, integration of these maps is essential, and biparental markers can be used for such an operation. The objective of our study was to construct an integrated molecular map for a full-sib population of yellow passion fruit combining different loci configuration generated from amplified fragment length polymorphisms (AFLPs) and microsatellite markers and using a novel approach based on simultaneous maximum-likelihood estimation of linkage and linkage phases, specially designed for outcrossing species. Of the total number of loci, approximate to 76%, 21%, 0.7%, and 2.3% did segregate in 1:1, 3:1, 1:2:1, and 1:1:1:1 ratios, respectively. Ten linkage groups (LGs) were established with a logarithm of the odds (LOD) score >= 5.0 assuming a recombination fraction : <= 0.35. On average, 24 markers were assigned per LG, representing a total map length of 1687 cM, with a marker density of 6.9 cM. No markers were placed as accessories on the map as was done with previously constructed individual maps.

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.

Random regression models to estimate test-day milk yield genetic parameters Holstein cows in Southeastern Brazil

A total of 152,145 weekly test-day milk yield records from 7317 first lactations of Holstein cows distributed in 93 herds in southeastern Brazil were analyzed. Test-day milk yields were classified into 44 weekly classes of DIM. The contemporary groups were defined as herd-year-week of test-day. The model included direct additive genetic, permanent environmental and residual effects as random and fixed effects of contemporary group and age of cow at calving as covariable, linear and quadratic effects. Mean trends were modeled by a cubic regression on orthogonal polynomials of DIM. Additive genetic and permanent environmental random effects were estimated by random regression on orthogonal Legendre polynomials. Residual variances were modeled using third to seventh-order variance functions or a step function with 1, 6,13,17 and 44 variance classes. Results from Akaike`s and Schwarz`s Bayesian information criterion suggested that a model considering a 7th-order Legendre polynomial for additive effect, a 12th-order polynomial for permanent environment effect and a step function with 6 classes for residual variances, fitted best. However, a parsimonious model, with a 6th-order Legendre polynomial for additive effects and a 7th-order polynomial for permanent environmental effects, yielded very similar genetic parameter estimates. (C) 2008 Elsevier B.V. All rights reserved.

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.

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.