A data assimilation method for using low-resolution Earth observation data in heterogeneous ecosystems

We present an approach for dealing with coarse-resolution Earth observations (EO) in terrestrial ecosystem data assimilation schemes. The use of coarse-scale observations in ecological data assimilation schemes is complicated by spatial heterogeneity and nonlinear processes in natural ecosystems. If these complications are not appropriately dealt with, then the data assimilation will produce biased results. The “disaggregation” approach that we describe in this paper combines frequent coarse-resolution observations with temporally sparse fine-resolution measurements. We demonstrate the approach using a demonstration data set based on measurements of an Arctic ecosystem. In this example, normalized difference vegetation index observations are assimilated into a “zero-order” model of leaf area index and carbon uptake. The disaggregation approach conserves key ecosystem characteristics regardless of the observation resolution and estimates the carbon uptake to within 1% of the demonstration data set “truth.” Assimilating the same data in the normal manner, but without the disaggregation approach, results in carbon uptake being underestimated by 58% at an observation resolution of 250 m. The disaggregation method allows the combination of multiresolution EO and improves in spatial resolution if observations are located on a grid that shifts from one observation time to the next. Additionally, the approach is not tied to a particular data assimilation scheme, model, or EO product and can cope with complex observation distributions, as it makes no implicit assumptions of normality.

The boreal forest and the polar front

The analysis presented in this paper suggests that the larger heating over the boreal forest in the spring and summer, as contrasted with weaker heating over the adjacent tundra, results in a preferred position of the polar front along the northern edge of the boreal forest. This positioning is well documented in the literature (see, for example, Bryson, 1966; Barry and Hare, 1974; Kreps and Barry, 1970). This heating results from the lower albedo of the boreal forest which is not compensated by an increase in transpiration, even with the larger leaf area index of the forest. The warmer temperatures are mixed upward by the deep boundary layer over the forest and mesoscale circulations which result from the patchiness of heating associated with the heterogeneous landscapes of the forest. Thus in contrast to previous assumptions in which the arctic front position in the summer determines the northern limit of the boreal tree line, our study suggests the boreal forest itself significantly influences the preferred position of the front. This conclusion reinforces the findings of Bonan et al. (1992) and Foley et al. (1994) on the important role of boreal forest-tundra interactions with climate.

The simulation of the opposing fluxes of latent heat and CO2 over various land-use types: coupling a gas exchange model to a mesoscale atmospheric model

A mesoscale meteorological model (FOOT3DK) is coupled with a gas exchange model to simulate surface fluxes of CO2 and H2O under field conditions. The gas exchange model consists of a C3 single leaf photosynthesis sub-model and an extended big leaf (sun/shade) sub-model that divides the canopy into sunlit and shaded fractions. Simulated CO2 fluxes of the stand-alone version of the gas exchange model correspond well to eddy-covariance measurements at a test site in a rural area in the west of Germany. The coupled FOOT3DK/gas exchange model is validated for the diurnal cycle at singular grid points, and delivers realistic fluxes with respect to their order of magnitude and to the general daily course. Compared to the Jarvis-based big leaf scheme, simulations of latent heat fluxes with a photosynthesis-based scheme for stomatal conductance are more realistic. As expected, flux averages are strongly influenced by the underlying land cover. While the simulated net ecosystem exchange is highly correlated with leaf area index, this correlation is much weaker for the latent heat flux. Photosynthetic CO2 uptake is associated with transpirational water loss via the stomata, and the resulting opposing surface fluxes of CO2 and H2O are reproduced with the model approach. Over vegetated surfaces it is shown that the coupling of a photosynthesis-based gas exchange model with the land-surface scheme of a mesoscale model results in more realistic simulated latent heat fluxes.

Validation of Canopy Height Profile methodology for small-footprint full-waveform airborne LiDAR data in a discontinuous canopy environment

A Canopy Height Profile (CHP) procedure presented in Harding et al. (2001) for large footprint LiDAR data was tested in a closed canopy environment as a way of extracting vertical foliage profiles from LiDAR raw-waveform. In this study, an adaptation of this method to small-footprint data has been shown, tested and validated in an Australian sparse canopy forest at plot- and site-level. Further, the methodology itself has been enhanced by implementing a dataset-adjusted reflectance ratio calculation according to Armston et al. (2013) in the processing chain, and tested against a fixed ratio of 0.5 estimated for the laser wavelength of 1550nm. As a by-product of the methodology, effective leaf area index (LAIe) estimates were derived and compared to hemispherical photography-derived values. To assess the influence of LiDAR aggregation area size on the estimates in a sparse canopy environment, LiDAR CHPs and LAIes were generated by aggregating waveforms to plot- and site-level footprints (plot/site-aggregated) as well as in 5m grids (grid-processed). LiDAR profiles were then compared to leaf biomass field profiles generated based on field tree measurements. The correlation between field and LiDAR profiles was very high, with a mean R2 of 0.75 at plot-level and 0.86 at site-level for 55 plots and the corresponding 11 sites. Gridding had almost no impact on the correlation between LiDAR and field profiles (only marginally improvement), nor did the dataset-adjusted reflectance ratio. However, gridding and the dataset-adjusted reflectance ratio were found to improve the correlation between raw-waveform LiDAR and hemispherical photography LAIe estimates, yielding the highest correlations of 0.61 at plot-level and of 0.83 at site-level. This proved the validity of the approach and superiority of dataset-adjusted reflectance ratio of Armston et al. (2013) over a fixed ratio of 0.5 for LAIe estimation, as well as showed the adequacy of small-footprint LiDAR data for LAIe estimation in discontinuous canopy forests.

JULES-crop: a parametrisation of crops in the Joint UK Land Environment Simulator

Studies of climate change impacts on the terrestrial biosphere have been completed without recognition of the integrated nature of the biosphere. Improved assessment of the impacts of climate change on food and water security requires the development and use of models not only representing each component but also their interactions. To meet this requirement the Joint UK Land Environment Simulator (JULES) land surface model has been modified to include a generic parametrisation of annual crops. The new model, JULES-crop, is described and evaluation at global and site levels for the four globally important crops; wheat, soybean, maize and rice. JULES-crop demonstrates skill in simulating the inter-annual variations of yield for maize and soybean at the global and country levels, and for wheat for major spring wheat producing countries. The impact of the new parametrisation, compared to the standard configuration, on the simulation of surface heat fluxes is largely an alteration of the partitioning between latent and sensible heat fluxes during the later part of the growing season. Further evaluation at the site level shows the model captures the seasonality of leaf area index, gross primary production and canopy height better than in the standard JULES. However, this does not lead to an improvement in the simulation of sensible and latent heat fluxes. The performance of JULES-crop from both an Earth system and crop yield model perspective is encouraging. However, more effort is needed to develop the parametrisation of the model for specific applications. Key future model developments identified include the introduction of processes such as irrigation and nitrogen limitation which will enable better representation of the spatial variability in yield.

Modelling spectral reflectance of open cork oak woodland: a simulation analysis of the effects of vegetation structure and background

This study analyses the influence of vegetation structure (i.e. leaf area index and canopy cover) and seasonal background changes on moderate-resolution imaging spectrometer (MODIS)-simulated reflectance data in open woodland. Approximately monthly spectral reflectance and transmittance field measurements (May 2011 to October 2013) of cork oak tree leaves (Quercus suber) and of the herbaceous understorey were recorded in the region of Ribatejo, Portugal. The geometric-optical and radiative transfer (GORT) model was used to simulate MODIS response (red, near-infrared) and to calculate vegetation indices, investigating their response to changes in the structure of the overstorey vegetation and to seasonal changes in the understorey using scenarios corresponding to contrasting phenological status (dry season vs. wet season). The performance of normalized difference vegetation index (NDVI), soil-adjusted vegetation index (SAVI), and enhanced vegetation index (EVI) is discussed. Results showed that SAVI and EVI were very sensitive to the emergence of background vegetation in the wet season compared to NDVI and that shading effects lead to an opposing trend in the vegetation indices. The information provided by this research can be useful to improve our understanding of the temporal dynamic of vegetation, monitored by vegetation indices.

A tillering inhibition gene influences root-shoot carbon partitioning and pattern of water use to improve wheat productivity in rainfed environments

Genetic modification of shoot and root morphology has potential to improve water and nutrient 19 uptake of wheat crops in rainfed environments. Near-isogenic lines (NILs) varying for a tillering 20 inhibition (tin) gene and representing multiple genetic backgrounds were investigated in contrasting 21 controlled environments for shoot and root growth. Leaf area, shoot and root biomass were similar 22 until tillering whereupon reduced tillering in tin-containing NILs produced reductions of up to 60% in 23 total leaf area and biomass, and increases in total root length of up to 120% and root biomass to 24 145%. Together, root-to-shoot ratio increased two-fold with the tin gene. The influence of tin on shoot 25 and root growth was greatest in the cv. Banks genetic background, particularly in the biculm-selected 26 NIL, and was typically strongest in cooler environments. A separate de-tillering study confirmed 27 greater root-to-shoot ratios with regular tiller removal in non-tin containing genotypes. In validating 28 these observations in a rainfed field study, the tin allele had a negligible effect on seedling growth but 29 was associated with significantly (P<0.05) reduced tiller number (-37%), leaf area index (-26%) and 30 spike number (-35%) to reduce plant biomass (-19%) at anthesis. Root biomass, root-to-shoot ratio at 31 early stem elongation and root depth at maturity were increased in tin-containing NILs. Soil water use 32 was slowed in tin-containing NILs resulting in greater water availability, greater stomatal 33 conductance, cooler canopy temperatures and maintenance of green leaf area during grain-filling. 34 Together these effects contributed to increases in harvest index and grain yield. In both the controlled 35 and field environments, the tin gene was commonly associated with increased root length and biomass 36 but the significant influence of genetic background and environment suggests careful assessment of 37 tin-containing progeny in selection for genotypic increases in root growth.

CHP toolkit: case study of LAIe sensitivity to discontinuity of canopy cover in fruit plantations

This paper presents an open-source canopy height profile (CHP) toolkit designed for processing small-footprint full-waveform LiDAR data to obtain the estimates of effective leaf area index (LAIe) and CHPs. The use of the toolkit is presented with a case study of LAIe estimation in discontinuous-canopy fruit plantations. The experiments are carried out in two study areas, namely, orange and almond plantations, with different percentages of canopy cover (48% and 40%, respectively). For comparison, two commonly used discrete-point LAIe estimation methods are also tested. The LiDAR LAIe values are first computed for each of the sites and each method as a whole, providing “apparent” site-level LAIe, which disregards the discontinuity of the plantations’ canopies. Since the toolkit allows for the calculation of the study area LAIe at different spatial scales, between-tree-level clumpingcan be easily accounted for and is then used to illustrate the impact of the discontinuity of canopy cover on LAIe retrieval. The LiDAR LAIe estimates are therefore computed at smaller scales as a mean of LAIe in various grid-cell sizes, providing estimates of “actual” site-level LAIe. Subsequently, the LiDAR LAIe results are compared with theoretical models of “apparent” LAIe versus “actual” LAIe, based on known percent canopy cover in each site. The comparison of those models to LiDAR LAIe derived from the smallest grid-cell sizes against the estimates of LAIe for the whole site has shown that the LAIe estimates obtained from the CHP toolkit provided values that are closest to those of theoretical models.

Fluxes of CO(2) above a plantation of Eucalyptus in southeast Brazil

Eddy-covariance measurements of net ecosystem exchange of CO(2) (NEE) and estimates of gross ecosystem productivity (GEP) and ecosystem respiration (R(E)) were obtained in a 2-4 year old Eucalyptus plantation during two years with very different winter rainfall In the first (drier) year the annual NEE GEP and RE were lower than the sums in the second (normal) year and conversely the total respiratory costs of assimilated carbon were higher in the dry year than in the normal year Although the net primary production (NPP) in the first year was 23% lower than that of the second year the decrease in the carbon use efficiency (CUE = NPP/GEP) was 11% and autotrophic respiration utilized more resources in the first dry year than in the second normal year The time variations in NEE were followed by NPP because in these young Eucalyptus plantations NEE is very largely dominated by NPP and heterotrophic respiration plays only a relatively minor role During the dry season a pronounced hysteresis was observed in the relationship between NEE and photosynthetically active radiation and NEE fluxes were inversely proportional to humidity saturation deficit values greater than 0 8 kPa Nighttime fluxes of CO(2) during calm conditions when the friction velocity (u) was below the threshold (0 25 ms(-1)) were estimated based on a Q(10) temperature-dependence relationship adjusted separately for different classes of soil moisture content which regulated the temperature sensitivity of ecosystem respiration (C) 2010 Elsevier B V All rights reserved

Estimativa da área foliar de crotalaria juncea l. a partir de dimensões lineares do limbo foliar

Knowledge of the leaf area plant are needed for agronomic and physiological studies involving plant growth. The aim of this study was to obtain a mathematical model using linear measures of leaf dimensions, which will allow the estimation of leaf area of Crotalaria juncea L. Correlation studies were conducted involving real leaf area (Sf) and leaf length (C), maximum leaf width (L) and the product between C and L. All tested models (linear, exponential or geometric) provided good estimation of leaf area (above 87%). The better fit was attained using linear model, passing or not through the origin. From a practical viewpoint, it is suggested to use the linear model involving the C and L product, using a linear coefficient equal to zero. Estimation of leaf area of Crotalaria juncea L. can be obtained using the model Sf = 0.7160 x (C*L) with a determination coefficient of 0.9712.

Análise de crescimento em amendoinzeiro (Arachis hypogaea L.) submetido à matocompetição

Foram verificados os efeitos de diferentes períodos de controle e convivência de plantas daninhas (0; 14; 28; 42; 56; 70; 84 e 98 dias) na cultura do amendoinzeiro (Arachis hypogaea L.) em Moji-Guaçu (SP). As principais espécies invasoras da área experimental foram: Acanthospermum australe (Loef.) O. Kuntz, Waltheria indica L., Borreria alata DC., Richardia scabra L., Rhynchelytrum repens (Nees) Stapf. et Hubb. e Diodia teres Walt. de maneira geral, durante quase todo o ciclo, a espécie de maior importância relativa foi A. australe, sendo suplantada por outras espécies apenas na colheita, quando já estava encerrando seu ciclo de desenvolvimento. O período inicial de controle das plantas daninhas, além de reduzir a densidade da comunidade infestante, diminuiu sua riqueza em espécies e alterou a importância relativa das populações. A matocompetição resultou na redução do crescimento do amendoinzeiro, principalmente a partir de 42 dias de convivência, afetando a altura das plantas, a área foliar, o índice de área foliar e o acúmulo epígeo de matéria seca. Não foram afetados os números de ramos e de folíolos por planta, nem a qualidade dos grãos e o vigor das sementes do amendoim. A produção de frutos foi reduzida ao redor de 31,5% em decorrência da interferência das plantas daninhas. O período total de prevenção da interferência na produção foi de 14 dias e o período anterior à interferência foi de 70 dias.

Validation of the Deardorff model for estimating energy balance components for a sugarcane crop

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

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

A melancia é uma espécie tradicionalmente conduzida em campo no sistema rasteiro. As cultivares de frutos pequenos (1 a 3 kg), que adquirem melhores preços de mercado, vêm sendo cultivadas também em ambiente protegido, onde são conduzidas no sistema vertical, com poda de ramos e raleio de frutos. Essas práticas possibilitam aumentar o adensamento das plantas, a qualidade e a produtividade de frutos em comparação ao sistema rasteiro. Objetivou-se com este trabalho avaliar a influência de três alturas de condução (1,7; 2,2 e 2,7 m) e duas densidades de plantas (3,17 e 4,76 plantas m-2) sobre as características produtivas e qualitativas da mini melancia Smile cultivada em ambiente protegido. A poda da haste principal foi realizada aos 43, 55 e 66 dias após o transplante (DAT) para as alturas de condução de 1,7; 2,2 e 2,7 m, respectivamente. A massa seca dos ramos, dos pecíolos, das folhas e total foram afetados pela altura de condução, cujos maiores valores foram obtidos para as plantas conduzidas a 2,2 e 2,7 m de altura. A área foliar, a área foliar específica e o índice de área foliar não foram influenciados pela altura de condução das plantas. A altura de condução de 2,7 m elevou a produtividade total. Entretanto, a produtividade comercial, a massa média dos frutos e todas as características qualitativas não foram significativamente diferentes das obtidos pela altura de poda de 2,2 m. em relação à densidade de plantas, a melhor opção foi a de 4,76 plantas m-2, pois elevou a produtividade comercial em 37,4% sem reduzir a massa média dos frutos.

Desenvolvimento da cultura do feijoeiro submetida a dois sistemas de manejo de irrigação e de cultivo

A cultura do feijoeiro apresenta elevada importância no contexto agrícola nacional. Além de ser uma cultura que atende às características da agricultura familiar, também é cultivada em áreas extensivas com alta tecnologia. Pesquisas relacionadas ao desenvolvimento do feijoeiro frente às variantes edafoclimáticas regionais vêm assumindo papel de destaque, principalmente por se tratar de uma planta sensível tanto ao déficit quanto ao excesso hídrico. O objetivo deste trabalho foi avaliar a cultura do feijoeiro de 'inverno' submetida a dois métodos de manejo de irrigação, em pivô central, um via solo por meio de tensiômetros, outro via clima por meio do tanque Classe 'A', em dois sistemas de plantio, convencional e direto, no primeiro ano. Foram avaliados o índice de área foliar (IAF), massa seca total (MS) de plantas e a variação do potencial mátrico de água no solo. Concluiu-se que os manejos de irrigação e os sistemas de plantio não resultaram em diferenças no IAF e nem no acúmulo de matéria seca total ao longo do ciclo da cultura; o manejo da irrigação por tensiometria acarretou maior variação no potencial mátrico de água no solo do que o método via clima por meio do tanque Classe 'A'.

Espaçamentos entrelinhas e entre plantas no crescimento e na produção de repolho roxo

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)