Do woody and herbaceous species compete for soil water across topographic gradients? Evidence for niche partitioning in a Neotropical savanna

Savannas are characterized by sparsely distributed woody species within a continuous herbaceous cover, composed mainly by grasses and small eudicot herbs. This vegetation structure is variable across the landscape, with shifts from open grassland to savanna woodland determined by factors that control tree density. These shifts often appear coupled with environmental variations, such as topographic gradients. Here we investigated whether herbaceous and woody savanna species differ in their use of soil water along a topographic gradient of about 110 m, spanning several vegetation physiognomies generally associated with Neotropical savannas. We measured the delta H-2 and delta O-18 signatures of plants, soils, groundwater and rainfall, determining the depth of plant water uptake and examining variations in water uptake patterns along the gradient. We found that woody species use water from deeper soil layers compared to herbaceous species, regardless of their position in the topographic gradient. However, the presence of a shallow water table restricted plant water uptake to the superficial soil layers at lower portions of the gradient. We confirmed that woody and herbaceous species are plastic with respect to their water use strategy, which determines niche partitioning across topographic gradients. Abiotic factors such as groundwater level, affect water uptake patterns independently of plant growth form, reinforcing vegetation gradients by exerting divergent selective pressures across topographic gradients. (C) 2013 SAAB. Published by Elsevier B.V. All rights reserved.

Root behavior of savanna species in Brazil's Pantanal wetland.

The objective of this study was to determine the maximum depth, structure, diameter and biomass of the roots of common woody species in two savanna physiognomies (savanna woodland and open woody savanna) in Brazil's Pantanal wetland. The root systems of 37 trees and 34 shrubs of 15 savanna species were excavated to measure their length and depth and estimate the total root biomass through allometric relationships with stem diameter at ground level. In general, statistical regression models between root weight and stem diameter at ground level showed a significance of P < 0.05 and R2 values close to or above 0.8. The average depths of the root system in wetland savanna woodland and open woody savanna are 0.8 ± 0.3 m and 0.7 ± 0.2 m, respectively, and differ from the root systems of savanna woody species in non-flooding areas, whose depth usually ranges from 3 to 19 m.Weattribute this difference to the adaptation of woody plant to the shallow water table, particularly during the wet season. This singularity of woody species in wetland savannas is important when considering biomass and carbon stocks for national and global carbon inventories.

Scenario analysis of alternative vegetation management options on the greenhouse gas budget of two grazing businesses in north-eastern Australia

The emerging carbon economy will have a major impact on grazing businesses because of significant livestock methane and land-use change emissions. Livestock methane emissions alone account for similar to 11% of Australia's reported greenhouse gas emissions. Grazing businesses need to develop an understanding of their greenhouse gas impact and be able to assess the impact of alternative management options. This paper attempts to generate a greenhouse gas budget for two scenarios using a spread sheet model. The first scenario was based on one land-type '20-year-old brigalow regrowth' in the brigalow bioregion of southern-central Queensland. The 50 year analysis demonstrated the substantially different greenhouse gas outcomes and livestock carrying capacity for three alternative regrowth management options: retain regrowth (sequester 71.5 t carbon dioxide equivalents per hectare, CO2-e/ha), clear all regrowth (emit 42.8 t CO2-e/ha) and clear regrowth strips (emit 5.8 t CO2-e/ha). The second scenario was based on a 'remnant eucalypt savanna-woodland' land type in the Einasleigh Uplands bioregion of north Queensland. The four alternative vegetation management options were: retain current woodland structure (emit 7.4 t CO2-e/ha), allow woodland to thicken increasing tree basal area (sequester 20.7 t CO2-e/ha), thin trees less than 10 cm diameter (emit 8.9 t CO2-e/ha), and thin trees <20 cm diameter (emit 12.4 t CO2-e/ha). Significant assumptions were required to complete the budgets due to gaps in current knowledge on the response of woody vegetation, soil carbon and non-CO2 soil emissions to management options and land-type at the property scale. The analyses indicate that there is scope for grazing businesses to choose alternative management options to influence their greenhouse gas budget. However, a key assumption is that accumulation of carbon or avoidance of emissions somewhere on a grazing business (e.g. in woody vegetation or soil) will be recognised as an offset for emissions elsewhere in the business (e.g. livestock methane). This issue will be a challenge for livestock industries and policy makers to work through in the coming years.

Efeitos do tipo de floresta e da estrutura de habitat em assembléias de primatas no sudoeste da Amazônia

A heterogeneidade ambiental expressa diferenças naturais entre áreas e é um fator determinante para a riqueza e abundância local de primatas. Neste estudo nós investigamos a composição e estrutura de assembléias de primatas em quatro tipos de floresta: floresta de terra firme, florestas de igapó sazonalmente inundáveis por rios de águas claras (aberta e densa) e cerradão na Reserva Biológica do Guaporé, sudoeste da Amazônia Brasileira. Além disso, avaliamos a associação entre a ocorrência e abundância dos primatas com diferenças estruturais das florestas. Realizamos 617,8 km de censos pelo método de transecção linear (~154 km por tipo de floresta) e avaliamos a estrutura da vegetação em 108 parcelas de 200 m² (0,54 ha por tipo de floresta). Dez espécies de primatas foram registradas durante os 11 meses deste estudo. A floresta de terra firme apresentou o maior número de espécies e a maior densidade de primatas, principalmente devido à presença exclusiva de Callicebus moloch e a maior abundância de Sapajus apella. A elevada densidade de Ateles chamek na floresta aberta inundável foi preponderante para a maior biomassa de primatas neste tipo de floresta. Nas florestas inundáveis e na terra firme, Ateles chamek e Sapajus apella responderam juntas por mais de 70% da biomassa de primatas, e no cerradão apenas Sapajus apella foi responsável por 68% da biomassa. Diferenças entre tipos de floresta na composição específica e abundância relativa de primatas foram associadas com o regime de inundação e com algumas variáveis de estrutura de habitat (densidade de árvores no sub-bosque e no dossel, abertura do dossel, altura total do dossel e densidade de palmeiras e lianas). Nossos resultados reforçam a importância de paisagens heterogêneas na Amazônia, pois estas áreas tendem a contribuir para uma maior diversidade de espécies em uma escala de paisagem.

Bauru: terra de cerrado ou floresta?

Bauru: region of cerrado ou forest? The city of Bauru is located in São Paulo State, where the maps showing native vegetation distribution does not clearly defi ne its nature. In order to clarify which types of native vegetation was found in this region, a bibliographic survey was performed, ranging from ancient documents prepared by naturalists describing this region and recent research results of floristic and phytosociological character. It is concluded that in Bauru, semideciduous seasonal forests overlays the northwestern and that cerrado comes over the southeast region. In riparian areas where the cerrado prevailed, there are still traces of swamp forests and swamp grasslands. In the border areas between forests and cerrado, the occurrence of vegetation transition is common, with two distinct types of it and other typical of these ecotone zones.

Fatores abióticos condicionantes da distribuição de espécies arbóreas em quatro formações florestais do Estado de São Paulo

No estudo das comunidades florestais, estabelecer a importância relativa dos fatores que definem a composição e a distribuição das espécies é um desafio. Em termos de gradientes ambientais o estudo das respostas das espécies arbóreas são essenciais para a compreensão dos processos ecológicos e decisões de conservação. Neste sentido, para contribuir com a elucidação dos processos ecológicos nas principais formações florestais do Estado de São Paulo (Floresta Ombrófila Densa de Terras Baixas, Floresta Ombrófila Densa Submontana, Floresta Estacional Semidecidual e Savana Florestada) este trabalho objetivou responder as seguintes questões: (I) a composição florística e a abundância das espécies arbóreas, em cada unidade fitogeográfica, variam conforme o gradiente edáfico e topográfico?; (II) características do solo e topografia podem influenciar na previsibilidade de ocorrência de espécies arbóreas de ampla distribuição em diferentes tipos vegetacionais? (III) existe relação entre o padrão de distribuição espacial de espécies arbóreas e os parâmetros do solo e topografia? O trabalho foi realizado em parcelas alocadas em unidades de conservação (UC) que apresentaram trechos representativos, em termos de conservação e tamanho, das quatro principais formações florestais presentes no Estado de São Paulo. Em cada UC foram contabilizados os indivíduos arbóreos (CAP ≥ 15 cm), topografia, dados de textura e atributos químicos dos solos em uma parcela de 10,24 ha, subdividida em 256 subparcelas. Análises de correspodência canônica foram aplicadas para estabelecer a correspondência entre a abundância das espécies e o gradiente ambiental (solo e topografia). O método TWINSPAN modificado foi aplicado ao diagrama de ordenação da CCA para avaliar a influência das variáveis ambientais (solo e topografia) na composição de espécies. Árvores de regressão \"ampliadas\" (BRT) foram ajustadas para a predição da ocorrência das espécies segundo as variáveis de solo e topografia. O índice de Getis-Ord (G) foi utilizado para determinar a autocorrelação espacial das variáveis ambientais utilizadas nos modelos de predição da ocorrência das espécies. Nas unidades fitogeográficas analisadas, a correspondência entre o gradiente ambiental (solo e topografia) e a abundância das espécies foi significativa, especialmente na Savana Florestada onde observou-se a maior relação. O solo e a topografia também se relacionaram com a semelhança na composição florística das subparcelas, com exceção da Floresta Estacional Semicidual (EEC). As principais variáveis de solo e topografia relacionadas a flora em cada UC foram: (1) Na Floresta Ombrófila Densa de Terras Baixas (PEIC) - teor de alumínio na camada profunda (Al (80-100 cm)) que pode refletir os teor de Al na superfície, acidez do solo (pH(H2O) (5-25 cm)) e altitude, que delimitou as áreas alagadas; (2) Na Floresta Ombrófila Densa Submontana (PECB) - altitude, fator que, devido ao relevo acidentado, influencia a temperatura e incidência de sol no sub-bosque; (3) Na Savana Florestada (EEA) - fertilidade, tolerância ao alumínio e acidez do solo. Nos modelos de predição BRT, as variáveis químicas dos solos foram mais importantes do que a textura, devido à pequena variação deste atributo no solo nas áreas amostradas. Dentre as variáveis químicas dos solos, a capacidade de troca catiônica foi utilizada para prever a ocorrência das espécies nas quatro formações florestais, sendo particularmente importante na camada mais profunda do solo da Floresta Ombrófila Densa de Terras Baixas (PEIC). Quanto à topografia, a altitude foi inserida na maioria dos modelos e apresentou diferentes influências sobre as áreas de estudo. De modo geral, para presença das espécies de ampla distribuição observou-se uma mesma tendência quando à associação com os atributos dos solos, porém com amplitudes dos descritores edáficos que variaram de acordo com a área de estudo. A ocorrência de Guapira opposita e Syagrus romanzoffiana, cujo padrão variou conforme a escala, foi explicada por variáveis com padrões espaciais agregados que somaram entre 30% e 50% de importância relativa no modelo BRT. A presença de A. anthelmia, cujo padrão também apresentou certo nível de agregação, foi associada apenas a uma variável com padrão agregado, a altitude (21%), que pode ter exercido grande influência na distribuição da espécie ao delimitar áreas alagadas. T. guianensis se associou a variáveis ambientais preditoras com padrão espacial agregado que somaram cerca de 70% de importância relativa, o que deve ter sido suficiente para estabelecer o padrão agregado em todas as escalas. No entanto, a influência dos fatores ambientais no padrão de distribuição da espécie não depende apenas do ótimo ambiental da espécie, mas um resultado da interação espécie-ambiente. Concluiu-se que: (I) características edáficas e topográficas explicaram uma pequena parcela da composição florística, em cada unidade fitogeográfica, embora a ocorrência de algumas espécies tenha se associado ao gradiente edáfico e topográfico; (II) a partir de características dos solos e da topografia foi possível prever a presença de espécies arbóreas, que apresentaram particularidades em relação a sua associação com o solo de cada fitofisionomia; (III) a partir de associações descritivas o solo e a topografia influenciam o padrão de distribuição espacial das espécies, na proporção em que contribuem para a presença das mesmas.

Long-term effects of Macroternies termites, herbivores and annual early fire on woody undergrowth community in Sudanian woodland, Burkina Faso

Termites, herbivores and fire are recognized as major guilds that structure woody plant communities in African savanna and woodland ecosystems. An understanding of their interaction is crucial to design appropriate management regimes. The aim of this study was to evaluate the long-term impacts of herbivore, fire and termite activities on regeneration of trees. Permanent experimental quadrats were established in 1992 in the Sudanian woodland of Burkina Faso subjected to grazing by livestock and annual early fire and the control. Within the treatment quadrats, an inventory of the woody undergrowth community was conducted on termitaria occupied by Macrotermes subhyalinus, extended termitosphere (within 5 m radius from the mound base) and adjacent area (beyond 5 m from the mound base). Hierarchical analysis was performed to determine significant differences in species richness, abundance and diversity indices among vegetation patches within fire and herbivory treatments. Grazed quadrats had significantly (P < 0.001) more species and stem density of woody undergrowth than non-grazed quadrats but maintained similar level of species richness and stem density of woody undergrowth on termitaria. There were not significant differences (P>0.05) in species richness and stem density between burnt and unburnt quadrats. Termitaria supported a highly diverse woody undergrowth with higher stem density than either the extended termitosphere or rest of quadrats. The density of woody undergrowth was significantly related with mature trees of selected species on termitaria (R-2 = 0.593; P<0.001) than that on the extended termitosphere (R-2 = 0.333; P<0.001) and adjacent area (R-2 = 0.197; P<0.001). It can be concluded that termites facilitate the regeneration of woody species while grazing and annual early fire play a minor role in the regeneration of woody species. The current policy that prohibits grazing should be revised to accommodate the interests of livestock herders. (C) 2014 Elsevier GmbH. All rights reserved.

Global fire emissions and the contribution of deforestation, savanna, forest, agricultural, and peat fires (1997-2009)

New burned area datasets and top-down constraints from atmospheric concentration measurements of pyrogenic gases have decreased the large uncertainty in fire emissions estimates. However, significant gaps remain in our understanding of the contribution of deforestation, savanna, forest, agricultural waste, and peat fires to total global fire emissions. Here we used a revised version of the Carnegie-Ames-Stanford-Approach (CASA) biogeochemical model and improved satellite-derived estimates of area burned, fire activity, and plant productivity to calculate fire emissions for the 1997-2009 period on a 0.5° spatial resolution with a monthly time step. For November 2000 onwards, estimates were based on burned area, active fire detections, and plant productivity from the MODerate resolution Imaging Spectroradiometer (MODIS) sensor. For the partitioning we focused on the MODIS era. We used maps of burned area derived from the Tropical Rainfall Measuring Mission (TRMM) Visible and Infrared Scanner (VIRS) and Along-Track Scanning Radiometer (ATSR) active fire data prior to MODIS (1997-2000) and estimates of plant productivity derived from Advanced Very High Resolution Radiometer (AVHRR) observations during the same period. Average global fire carbon emissions according to this version 3 of the Global Fire Emissions Database (GFED3) were 2.0 PgC year-1 with significant interannual variability during 1997-2001 (2.8 Pg Cyear-1 in 1998 and 1.6 PgC year-1 in 2001). Globally, emissions during 2002-2007 were rela-tively constant (around 2.1 Pg C year-1) before declining in 2008 (1.7 Pg Cyear-1) and 2009 (1.5 PgC year-1) partly due to lower deforestation fire emissions in South America and tropical Asia. On a regional basis, emissions were highly variable during 2002-2007 (e.g., boreal Asia, South America, and Indonesia), but these regional differences canceled out at a global level. During the MODIS era (2001-2009), most carbon emissions were from fires in grasslands and savannas (44%) with smaller contributions from tropical deforestation and degradation fires (20%), woodland fires (mostly confined to the tropics, 16%), forest fires (mostly in the extratropics, 15%), agricultural waste burning (3%), and tropical peat fires (3%). The contribution from agricultural waste fires was likely a lower bound because our approach for measuring burned area could not detect all of these relatively small fires. Total carbon emissions were on average 13% lower than in our previous (GFED2) work. For reduced trace gases such as CO and CH4, deforestation, degradation, and peat fires were more important contributors because of higher emissions of reduced trace gases per unit carbon combusted compared to savanna fires. Carbon emissions from tropical deforestation, degradation, and peatland fires were on average 0.5 PgC year-1. The carbon emissions from these fires may not be balanced by regrowth following fire. Our results provide the first global assessment of the contribution of different sources to total global fire emissions for the past decade, and supply the community with an improved 13-year fire emissions time series. © 2010 Author(s).

Birds of a habitat spectrum in the Itirapina Savanna, São Paulo, Brazil (1982-2003)

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Birds from cerradão woodland, an overlooked forest of the Cerrado region, Brazil

O Cerrado ainda recebe pouca atenção no que diz respeito à ornitologia embora seja a única savana tropical do mundo considerada um hotspot de biodiversidade. O cerradão é uma das fisionomias menos conhecidas e mais desmatadas do bioma e poucos levantamentos avifaunísticos foram realizados nessas florestas. Para revisar os estudos sobre aves de cerradão e complementar os poucos inventários já existentes realizados nesse tipo florestal no estado de São Paulo, foi realizado um levantamento bibliográfico dos estudos publicados sobre aves de cerradão. Adicionalmente foi conduzido um levantamento das aves de um fragmento de cerradão de 314 ha localizado na região central do estado de São Paulo, Brasil, entre setembro de 2005 e dezembro de 2006 com a utilização de transecções lineares com raio ilimitado de detecção. de 95 estudos envolvendo aves de cerradão, apenas 17 (18%) discriminaram espécies registradas dentro desta fisionomia daquelas que obtiveram registros em outros ambientes de Cerrado. Exceto por um estudo, nenhuma outra investigação encontrou mais de 64 espécies de aves neste ambiente, resultado compartilhado com diversas regiões do Brasil e também da Bolívia. Diferenças no número de espécies entre cerradões não puderam ser atribuídas à degradação dos ambientes estudados ou tamanho de fragmento. Considerando os registros de cerradões no Brasil e na Bolívia, a compilação de dados acumulou 250 espécies distribuídas em 36 famílias e 15 ordens. Durante nossos trabalhos de campo em localidade do interior paulista foram registradas 48 espécies distribuídas em 20 famílias, incluindo o fruxu-do-cerradão (Neopelma pallescens), ameaçada em São Paulo, e o soldadinho (Antilophia galeata), quase ameaçada no estado e endêmica do Cerrado. Dentre as espécies mais abundantes no fragmento, nenhuma delas é ameaçada ou endêmica do bioma.

Estimating the aboveground biomass and wood volume of savanna woodlands in Brazil's Pantanal wetlands based on allometric correlations

The quantity and distribution of vegetal biomass are important aspects to consider in ecosystem studies. However, little information is available about Brazil's Pantanal woodland savannas. This work involved the development of regression equations of the aerial biomass and wood volume of native tree species in a region of woodland savanna on Rio Negro farm in the Pantanal of Nhecolandia, Brazil. Samples were taken from 10 trees of each of five species: Protium heptaphyllum (Aub1.) Marchand, Magonia pubescens A. St.-Hil., Diptychandra aurantiaca Tul., Terminalia argentea Mart. and Zucc. and Licania minutiflora (Sagot) Fritsch and from a miscellaneous group of I I different species. Linear and nonlinear regression analyses were developed relating the diameter at breast height to the dry weight of wood, branches and leaves, wood volume and total aerial biomass. All the regressions showed a significance of P < 0.05 and an R-2 close to or above 0.8. The biomass curve predicted by linear regression analysis of the studied species was similar to the nonlinear regression, with the exception of L. minutiflora and the miscellaneous group. The breast height diameter proved a good choice for estimating biomass and wood volume. The estimated wood volume and biomass of the Pantanal woodland savanna is crucial information for understanding the carbon cycle and for ensuring the region's conservation and sustainable use. (c) 2006 Elsevier B.V. All rights reserved.

Land use monitoring in the Nigerian savanna using aerial photographs

Aerial photography was used to determine the land use in a test area of the Nigerian savanna in 1950 and 1972. Changes in land use were determined and correlated with accessibility, appropriate low technology methods being used to make it easy to extend the investigation to other areas without incurring great expense. A test area of 750 sq km was chosen located in Kaduna State of Nigeria. The geography of the area is summarised together with the local knowledge which is essential for accurate photo interpretation. A land use classification was devised and tested for use with medium scale aerial photography of the savanna. The two sets of aerial photography at 1:25 000 scale were sampled using systematic dot grids. A dot density of 8 1/2 dots per sq km was calculated to give an acceptable estimate of land use. Problems of interpretation included gradation between categories, sample position uncertainty and personal bias. The results showed that in 22 years the amount of cultivated land in the test area had doubled while there had been a corresponding decrease in the amount of uncultivated land particularly woodland. The intensity of land use had generally increased. The distribution of land use changes was analysed and correlated with accessibility. Highly significant correlations were found for 1972 which had not existed in 1950. Changes in land use could also be correlated with accessibility. It was concluded that in the 22 year test period there had been intensification of land use, movement of human activity towards the main road, and a decrease in natural vegetation particularly close to the road. The classification of land use and the dot grid method of survey were shown to be applicable to a savanna test area.

Viewpoint : social and economic dimensions of involving savanna communities in carbon management systems

In this paper we discuss the social, economic and institutional aspects of the development of carbon management systems within Australia's tropical savannas. Land-use values in savanna landscapes are changing as a result of changing economic markets, greater recognition of native title, and growing social demands and expectations for tourism, recreation and conservation. In addition, there is increasing interest in developing markets and policy arrangements for greenhouse gas abatement, carbon sequestration and carbon trade in savannas. We argue that for carbon management to lead to national greenhouse outcomes, attention must be paid to social, economic and institutional issues in environmental planning and policy arrangements. From an economic perspective, the financial impact of carbon management on savanna enterprises will depend on appropriate and available policy mechanisms, unit price for carbon, landscape condition, existing management strategies and abatement measurements used. Local social and cultural features of communities and regions may enhance or constrain the implementation of carbon abatement strategies, depending on how they are perceived. In terms of institutional arrangements, policies and plans must support and enable carbon management. We identify three areas that require priority investigation and adjustment: regional planning arrangements, property rights, and rules for accounting at enterprise and regional scales. We conclude that the best potential for managing for carbon will be achieved while managing for range of other natural resource management outcomes, especially where managing for carbon delivers collateral benefits to enterprises.

The genome of woodland strawberry (Fragaria vesca)

The woodland strawberry, Fragaria vesca (2n = 2x = 14), is a versatile experimental plant system. This diminutive herbaceous perennial has a small genome (240 Mb), is amenable to genetic transformation and shares substantial sequence identity with the cultivated strawberry (Fragaria Ã- ananassa) and other economically important rosaceous plants. Here we report the draft F. vesca genome, which was sequenced to ×-39 coverage using second-generation technology, assembled de novo and then anchored to the genetic linkage map into seven pseudochromosomes. This diploid strawberry sequence lacks the large genome duplications seen in other rosids. Gene prediction modeling identified 34,809 genes, with most being supported by transcriptome mapping. Genes critical to valuable horticultural traits including flavor, nutritional value and flowering time were identified. Macrosyntenic relationships between Fragaria and Prunus predict a hypothetical ancestral Rosaceae genome that had nine chromosomes. New phylogenetic analysis of 154 protein-coding genes suggests that assignment of Populus to Malvidae, rather than Fabidae, is warranted.

A regional interpretation of rules and good practice for greenhouse accounting : Northern Australian savanna systems

Land-use change, particularly clearing of forests for agriculture, has contributed significantly to the observed rise in atmospheric carbon dioxide concentration. Concern about the impacts on climate has led to efforts to monitor and curtail the rapid increase in concentrations of carbon dioxide and other greenhouse gases in the atmosphere. Internationally, much of the current focus is on the Kyoto Protocol to the United Nations Framework Convention on Climate Change (UNFCCC). Although electing to not ratify the Protocol, Australia, as a party to the UNFCCC, reports on national greenhouse gas emissions, trends in emissions and abatement measures. In this paper we review the complex accounting rules for human activities affecting greenhouse gas fluxes in the terrestrial biosphere and explore implications and potential opportunities for managing carbon in the savanna ecosystems of northern Australia. Savannas in Australia are managed for grazing as well as for cultural and environmental values against a background of extreme climate variability and disturbance, notably fire. Methane from livestock and non-CO2 emissions from burning are important components of the total greenhouse gas emissions associated with management of savannas. International developments in carbon accounting for the terrestrial biosphere bring a requirement for better attribution of change in carbon stocks and more detailed and spatially explicit data on such characteristics of savanna ecosystems as fire regimes, production and type of fuel for burning, drivers of woody encroachment, rates of woody regrowth, stocking rates and grazing impacts. The benefits of improved biophysical information and of understanding the impacts on ecosystem function of natural factors and management options will extend beyond greenhouse accounting to better land management for multiple objectives.