Calibration of a Species-Specific Spectral Vegetation Index for Leaf Area Index (LAI) Monitoring: Example with MODIS Reflectance Time-Series on Eucalyptus Plantations

The leaf area index (LAI) is a key characteristic of forest ecosystems. Estimations of LAI from satellite images generally rely on spectral vegetation indices (SVIs) or radiative transfer model (RTM) inversions. We have developed a new and precise method suitable for practical application, consisting of building a species-specific SVI that is best-suited to both sensor and vegetation characteristics. Such an SVI requires calibration on a large number of representative vegetation conditions. We developed a two-step approach: (1) estimation of LAI on a subset of satellite data through RTM inversion; and (2) the calibration of a vegetation index on these estimated LAI. We applied this methodology to Eucalyptus plantations which have highly variable LAI in time and space. Previous results showed that an RTM inversion of Moderate Resolution Imaging Spectroradiometer (MODIS) near-infrared and red reflectance allowed good retrieval performance (R-2 = 0.80, RMSE = 0.41), but was computationally difficult. Here, the RTM results were used to calibrate a dedicated vegetation index (called "EucVI") which gave similar LAI retrieval results but in a simpler way. The R-2 of the regression between measured and EucVI-simulated LAI values on a validation dataset was 0.68, and the RMSE was 0.49. The additional use of stand age and day of year in the SVI equation slightly increased the performance of the index (R-2 = 0.77 and RMSE = 0.41). This simple index opens the way to an easily applicable retrieval of Eucalyptus LAI from MODIS data, which could be used in an operational way.

Effects of irrigation on water use and water use efficiency in two fast growing Eucalyptus plantations

Eucalyptus plantations occupy almost 20 million ha worldwide and exceed 3.7 million ha in Brazil alone. Improved genetics and silviculture have led to as much as a three-fold increase in productivity in Eucalyptus plantations in Brazil and the large land area occupied by these highly productive ecosystems raises concern over their effect on local water supplies. As part of the Brazil Potential Productivity Project, we measured water use of Eucalyptus grandis x urophylla clones in rainfed and irrigated stands in two plantations differing in productivity. The Aracruz (lower productivity) site is located in the state of Espirito Santo and the Veracel (higher productivity) site in Bahia state. At each plantation, we measured stand water use using homemade sap flow sensors and a calibration curve using the clones and probes we utilized in the study. We also quantified changes in growth, leaf area and water use efficiency (the amount of wood produced per unit of water transpired). Measurements were conducted for 1 year during 2005 at Aracruz and from August through December 2005 at Veracel. Transpiration at both sites was high compared to other studies but annual estimates at Aracruz for the rainfed treatment compared well with a process model calibrated for the Aracruz site (within 10%). Annual water use at Aracruz was 1394 mm in rainfed treatments versus 1779 mm in irrigated treatments and accounted for approximately 67% and 58% of annual precipitation and irrigation inputs respectively. Increased water use in the irrigated stands at Aracruz was associated with higher sapwood area, leaf area index and transpiration per unit leaf area but there was no difference in the response of canopy conductance with air saturation deficit between treatments. Water use efficiency at the Aracruz site was also not influenced by irrigation and was similar to the rainfed treatment. During the period of overlapping measurements, the response to irrigation treatments at the more productive Veracel site was similar to Aracruz. Stand water use at the Veracel site totaled 975 mm and 1102 mm in rainfed and irrigated treatments during the 5-month measurement period respectively. Irrigated stands at Veracel also had higher leaf area with no difference in the response of canopy conductance with air saturation deficit between treatments. Water use efficiency was also unaffected by irrigation at Veracel. Results from this and other studies suggest that improved resource availability does not negatively impact water use efficiency but increased productivity of these plantations is associated with higher water use and should be given consideration during plantation management decision making processes aimed at increasing productivity. Published by Elsevier B.V.

Analysis of NDVI time series using cross-correlation and forecasting methods for monitoring sugarcane fields in Brazil

Brazil is the largest sugarcane producer in the world and has a privileged position to attend to national and international market places. To maintain the high production of sugarcane, it is fundamental to improve the forecasting models of crop seasons through the use of alternative technologies, such as remote sensing. Thus, the main purpose of this article is to assess the results of two different statistical forecasting methods applied to an agroclimatic index (the water requirement satisfaction index; WRSI) and the sugarcane spectral response (normalized difference vegetation index; NDVI) registered on National Oceanic and Atmospheric Administration Advanced Very High Resolution Radiometer (NOAA-AVHRR) satellite images. We also evaluated the cross-correlation between these two indexes. According to the results obtained, there are meaningful correlations between NDVI and WRSI with time lags. Additionally, the adjusted model for NDVI presented more accurate results than the forecasting models for WRSI. Finally, the analyses indicate that NDVI is more predictable due to its seasonality and the WRSI values are more variable making it difficult to forecast.

Validation and comparison of two soil-vegetation-atmosphere transfer models for tropical Africa

This study aims to compare and validate two soil-vegetation-atmosphere-transfer (SVAT) schemes: TERRA-ML and the Community Land Model (CLM). Both SVAT schemes are run in standalone mode (decoupled from an atmospheric model) and forced with meteorological in-situ measurements obtained at several tropical African sites. Model performance is quantified by comparing simulated sensible and latent heat fluxes with eddy-covariance measurements. Our analysis indicates that the Community Land Model corresponds more closely to the micrometeorological observations, reflecting the advantages of the higher model complexity and physical realism. Deficiencies in TERRA-ML are addressed and its performance is improved: (1) adjusting input data (root depth) to region-specific values (tropical evergreen forest) resolves dry-season underestimation of evapotranspiration; (2) adjusting the leaf area index and albedo (depending on hard-coded model constants) resolves overestimations of both latent and sensible heat fluxes; and (3) an unrealistic flux partitioning caused by overestimated superficial water contents is reduced by adjusting the hydraulic conductivity parameterization. CLM is by default more versatile in its global application on different vegetation types and climates. On the other hand, with its lower degree of complexity, TERRA-ML is much less computationally demanding, which leads to faster calculation times in a coupled climate simulation.

STRUCTURAL CHARACTERIZATION OF CANOPIES OF Eucalyptus spp. USING RADIOMETRIC DATA FROM TM/Landsat 5

Empirical approaches and, more recently, physical approaches, have grounded the establishment of logical connections between radiometric variables derived from remote data and biophysical variables derived from vegetation cover. This study was aimed at evaluating correlations of dendrometric and density data from canopies of Eucalyptus spp., as collected in Capao Bonito forest unit, with radiometric data from imagery acquired by the TM/Landsat-5 sensor on two orbital passages over the study site (dates close to field data collection). Results indicate that stronger correlations were identified between crown dimensions and canopy height with near-infrared spectral band data (rho(s)4), irrespective of the satellite passage date. Estimates of spatial distribution of dendrometric data and canopy density (D) using spectral characterization were consistent with the spatial distribution of tree ages during the study period. Statistical tests were applied to evaluate performance disparities of empirical models depending on which date data were acquired. Results indicated a significant difference between models based on distinct data acquisition dates.

Parametrização e avaliação do modelo DSSAT/Canegro para variedades brasileiras de cana-de-açúcar

O objetivo deste trabalho foi parametrizar e avaliar o modelo DSSAT/Canegro para cinco variedades brasileiras de cana-de-açúcar. A parametrização foi realizada a partir do uso de dados biométricos e de crescimento das variedades CTC 4, CTC 7, CTC 20, RB 86-7515 e RB 83-5486, obtidos em cinco localidades brasileiras. Foi realizada análise de sensibilidade local para os principais parâmetros. A parametrização do modelo foi feita por meio da técnica de estimativa da incerteza de probabilidade generalizada ("generalized likelihood uncertainty estimation", Glue). Para a avaliação das predições, foram utilizados, como indicadores estatísticos, o coeficiente de determinação (R2), o índice D de Willmott e a raiz quadrada do erro-médio (RMSE). As variedades CTC apresentaram índice D entre 0,870 e 0,944, para índice de área foliar, altura de colmo, perfilhamento e teor de sacarose. A variedade RB 83-5486 apresentou resultados similares para teor de sacarose e massa de matéria fresca do colmo, enquanto a variedade RB 86-7515 apresentou valores entre 0,665 e 0,873, para as variáveis avaliadas.

Vertical growth of mini watermelon according to the training height and plant density

The watermelon is traditionally cultivated horizontally on the ground. The cultivars of small fruits (1 to 3 kg), which reach better market prices, are also being grown in a greenhouse, where the plants are trained upward on vertical supports, with branches pruning and fruits thinning. These practices make possible an increase of the plant density, fruit quality and yield compared to the traditional growth system. The aim of this experiment was to evaluate the influence of three training heights (1.7, 2.2 and 2.7 m) and two planting densities (3.17 and 4.76 plants m-2) over the productive and qualitative characteristics of mini watermelon "Smile" cultivated in greenhouse. The pruning was done at 43, 55 and 66 days after transplanting (DAT), when the plant height reached 1.7, 2.2 and 2.7 m, respectively. The dry mass of branches, petioles, leaves and total were affected by the training height, where the highest values were obtained by the plants pruned at 2.2 and 2.7 m. Leaf area, specific leaf area and leaf area index were not affected by the height of the plants. The training height of 2.7 m raised the total yield, however, marketable yield, average fruit mass and all the quality characteristics did not differ significantly from those obtained by the training height of 2.2 m. Regarding to plant density, the best option was 4.76 plants m-2, due to the increasing of marketable yield in 37.4% without reducing the average weight of fruits.

Assessing biophysical variable parameters of bean crop with hyperspectral measurements

Recently high spectral resolution sensors have been developed, which allow new and more advanced applications in agriculture. Motivated by the increasing importance of hyperspectral remote sensing data, the need for research is important to define optimal wavebands to estimate biophysical parameters of crop. The use of narrow band vegetation indices (VI) derived from hyperspectral measurements acquired by a field spectrometer was evaluated to estimate bean (Phaseolus vulgaris L.) grain yield, plant height and leaf area index (LAI). Field canopy reflectance measurements were acquired at six bean growth stages over 48 plots with four water levels (179.5; 256.5; 357.5 and 406.2 mm) and tree nitrogen rates (0; 80 and 160 kg ha-1) and four replicates. The following VI was analyzed: OSNBR (optimum simple narrow-band reflectivity); NB_NDVI (narrow-band normalized difference vegetation index) and NDVI (normalized difference index). The vegetation indices investigated (OSNBR, NB_NDVI and NDVI) were efficient to estimate LAI, plant height and grain yield. During all crop development, the best correlations between biophysical variables and spectral variables were observed on V4 (the third trifoliolate leaves were unfolded in 50 % of plants) and R6 (plants developed first flowers in 50 % of plants) stages, according to the variable analyzed.

Modified crop model estimation of depleted and potential soybean yield

Despite the great importance of soybeans in Brazil, there have been few applications of soybean crop modeling on Brazilian conditions. Thus, the objective of this study was to use modified crop models to estimate the depleted and potential soybean crop yield in Brazil. The climatic variable data used in the modified simulation of the soybean crop models were temperature, insolation and rainfall. The data set was taken from 33 counties (28 Sao Paulo state counties, and 5 counties from other states that neighbor São Paulo). Among the models, modifications in the estimation of the leaf area of the soybean crop, which includes corrections for the temperature, shading, senescence, CO2, and biomass partition were proposed; also, the methods of input for the model's simulation of the climatic variables were reconsidered. The depleted yields were estimated through a water balance, from which the depletion coefficient was estimated. It can be concluded that the adaptation soybean growth crop model might be used to predict the results of the depleted and potential yield of soybeans, and it can also be used to indicate better locations and periods of tillage.

Estimativa do índice de área foliar e temperatura em uma área de eucalipto

It is known that the presence of large masses of vegetation is a factor that can influence the microclimate of a region. In this paper we analyzed the correlation between leaf area index (LAI) and land surface temperature (LST), both estimated from remote sensing images from Landsat-5 TM in an area of eucalyptus plantation, and these estimates were compared to the observed data. The correlation between LAI and LST was not significant (16%), which indicates that there is no necessarily a direct influence of vegetation in the local temperature. The comparison between estimated and observed data shows that the application of remote sensing techniques in the estimative of interested variables is efficient, because the estimatives followed consistently the observed values.

Predicting Growth of Panicum maximum: An Adaptation of the CROPGRO-Perennial Forage Model

Warm-season grasses are economically important for cattle production in tropical regions and tools to aid in management and research on these forages would be highly beneficial both in research and the industry. This research was conducted to adapt the CROPGRO-Perennial Forage model to simulate growth of the tropical species guineagrass (Panicum maximum Jacq. cv. 'Tanzania') and to describe model adaptation for this species. To develop the CROPGRO parameters for this species, we began with values and relationships reported in the literature. Some parameters and relationships were calibrated by comparison with observed growth, development, dry matter accumulation, and partitioning during a 17-mo experiment with Tanzania guineagrass in Piracicaba, SP, Brazil. Compared with starting parameters for palisadegrass [Brachiaria brizantha (A. Rich.) Stapf. cv. 'Xaraes'], dormancy effects of the perennial forage model had to be minimized, partitioning to storage tissue or root decreased, and partitioning to leaf and stem increased to provide for more leaf and stem growth and less root. Parameters affecting specific leaf area and senescence of plant tissues were improved. After these changes were made to the model, biomass accumulation was better simulated, mean predicted herbage yield was 6576 kg ha(-1), averaged across 11 regrowth cycles of 35 (summer) or 63 d (winter), with a RMSE of 494 kg ha(-1) (Willmott's index of agreement d = 0.985, simulated/observed ratio = 1.014). The model also gave good predictions against an independent data set, with similar RMSE, ratio, and d. The results of the adaptation suggest that the CROPGRO model is an efficient tool to integrate physiological aspects of guineagrass and can be used to simulate growth.

Production and carbon allocation in monocultures and mixed-species plantations of Eucalyptus grandis and Acacia mangium in Brazil

Introducing nitrogen-fixing tree species in fast-growing eucalypt plantations has the potential to improve soil nitrogen availability compared with eucalypt monocultures. Whether or not the changes in soil nutrient status and stand structure will lead to mixtures that out-yield monocultures depends on the balance between positive interactions and the negative effects of interspecific competition, and on their effect on carbon (C) uptake and partitioning. We used a C budget approach to quantify growth, C uptake and C partitioning in monocultures of Eucalyptus grandis (W. Hill ex Maiden) and Acacia mangium (Willd.) (treatments E100 and A100, respectively), and in a mixture at the same stocking density with the two species at a proportion of 1 : 1 (treatment MS). Allometric relationships established over the whole rotation, and measurements of soil CO2 efflux and aboveground litterfall for ages 4-6 years after planting were used to estimate aboveground net primary production (ANPP), total belowground carbon flux (TBCF) and gross primary production (GPP). We tested the hypotheses that (i) species differences for wood production between E. grandis and A. mangium monocultures were partly explained by different C partitioning strategies, and (ii) the observed lower wood production in the mixture compared with eucalypt monoculture was mostly explained by a lower partitioning aboveground. At the end of the rotation, total aboveground biomass was lowest in A100 (10.5 kg DM m(-2)), intermediate in MS (12.2 kg DM m(-2)) and highest in E100 (13.9 kg DM m(-2)). The results did not support our first hypothesis of contrasting C partitioning strategies between E. grandis and A. mangium monocultures: the 21% lower growth (delta B-w) in A100 compared with E100 was almost entirely explained by a 23% lower GPP, with little or no species difference in ratios such as TBCF/GPP, ANPP/TBCF, delta B-w/ANPP and delta B-w/GPP. In contrast, the 28% lower delta B-w in MS than in E100 was explained both by a 15% lower GPP and by a 15% lower fraction of GPP allocated to wood growth, thus partially supporting our second hypothesis: mixing the two species led to shifts in C allocations from above- to belowground, and from growth to litter production, for both species.

Corte foliar e tempo de transplantio para o uso de plântulas do sub-bosque na restauração florestal

O uso de plântulas da regeneração natural tem sido recomendado como estratégia para produção de mudas visando à restauração florestal, contudo muitos aspectos técnicos desse método ainda carecem de investigação científica. O objetivo deste trabalho foi avaliar o efeito da redução da área foliar e do transplantio imediato na sobrevivência e crescimento de mudas de espécies arbóreas produzidas a partir de plântulas obtidas da regeneração natural. Plântulas de Esenbeckia leiocarpa (Rutaceae), Eugenia ligustrina (Myrtaceae) e Maytenus salicifolia (Celastraceae), obtidas em remanescente de vegetação secundária de Floresta Estacional Semidecidual em Bofete, SP, foram extraídas do solo e submetidas aos tratamentos: I) redução de 50% da área de cada folha e transplantio imediato; II) nenhuma redução de área das folhas e transplantio imediato; III) redução de 50% da área de cada folha, manutenção das plântulas em água e transplantio 24 h após a coleta; e IV) nenhuma redução de área das folhas, manutenção das plântulas em água e transplantio 24 h após a coleta. As mudas foram avaliadas com relação à sobrevivência e ao crescimento em altura, ao longo de oito meses. Os resultados evidenciaram que nem o corte das folhas ou a manutenção das plântulas dentro de recipientes com água por 24 h antes do transplantio afetaram os parâmetros avaliados. Assim, para as espécies estudadas a redução da área foliar e o transplantio imediato são desnecessários para a produção de mudas em viveiro a partir de plântulas obtidas da regeneração natural.

Growth of the crabgrass species Digitaria ciliaris and Digitaria nuda

The aim of this research paper was to compare the growth of D. ciliaris and D. nuda crabgrass species under non-competitive conditions. To this end, two experiments were conducted, one from March - July 2010 and the other from February - June 2011. The experimental design of both trials was completely randomized making a factorial (2 seasons x 2 species crabgrass x 12 evaluation periods) with four replications. Assessments began at 15 days after sowing (DAS), and repeated weekly until 92 DAS. The variables evaluated were total dry matter (roots+leaves+stems), leaf area, leaf number and tiller. The results were submitted to analysis of variance and the absolute growth rate, relative growth rate and leaf area ratio were calculated using the means, which were adjusted regression models. The crabgrass species were significantly different in leaf area, leaf number, tiller number and dry matter per plant. D. ciliaris for all variables was statistically higher than D. nuda. Regarding the speed at which the growth of the species occurred, the absolute growth rate and relative growth rate of D. ciliaris was also greater than D. nuda. In addition, D. ciliaris also had a lower leaf area ratio indicating greater efficiency in converting light energy into carbohydrates. It can be concluded that D. ciliaris has a higher growth rate in conditions where there is no limitation of nutrients and water availability in relation to D. nuda, mainly due to D. ciliaris have greater leaf area, number of leaves and dry matter accumulation per plant.

Crescimento inicial de Merremia cissoides, Neonotonia wightii e Stizolobium aterrimum

Este trabalho foi desenvolvido com o objetivo de avaliar o crescimento inicial de três espécies de plantas daninhas que são problemas emergentes para a cultura da cana-de-açúcar: Neonotonia wightii, Stizolobium aterrimum e Merremia cissoides. O delineamento experimental adotado foi inteiramente casualizado, com sete tratamentos e quatro repetições, para cada espécie estudada. Durante todo o experimento foram realizadas sete avaliações de crescimento (tratamentos), espaçadas de 15 dias, totalizando 105 dias de ciclo. As espécies tiveram suas variáveis analisadas, como a área foliar (Af; cm² por planta). Posteriormente, o material amostrado foi seco em estufa, quando se mensurou também a matéria seca (g por planta) das raízes (Mr), da parte aérea (Ma) e total (Mt). O acúmulo de matéria seca observado mostra a elevada capacidade que a planta daninha S. aterrimum possui de produzir fitomassa. Ao final do experimento, observou-se estabilização da matéria seca da parte aérea (Ma) nas três espécies estudadas, acompanhada pela variável matéria seca total (Mt), com estabilização no final do experimento. Para S. aterrimum, obteve-se maior taxa de crescimento inicial e acúmulo de matéria seca total (Mt), além de maior acúmulo de área foliar (Af). N. wightii mostrou ser uma espécie de crescimento inicial lento, quando comparada com as outras estudadas. M. cissoides apresentou acúmulo significativo de matéria seca total (Mt), quando comparado com as outras espécies.