European biorefineries:implications for land, trade and employment

Biorefineries are expected to play a major role in a future low carbon economy and substantial investments are being made to support this vision. However, it is important to consider the wider socio-economic impacts of such a transition. This paper quantifies the potential trade, employment and land impacts of economically viable European biorefinery options based on indigenous straw and wood feedstocks. It illustrates how there could be potential for 70-80 European biorefineries, but not hundreds. A single facility could generate tens of thousands of man-years of employment and employment creation per unit of feedstock is higher than for biomass power plants. However, contribution to national GDP is unlikely to exceed 1% in European member states, although contributions to national agricultural productivity may be more significant, particularly with straw feedstocks. There is also a risk that biorefinery development could result in reduced rates of straw incorporation into soil, raising concerns that economically rational decisions to sell rather than reincorporate straw could result in increased agricultural land-use or greenhouse gas emissions. © 2013.

A tapestry of land uses and land covers typify our study site. Here, a commercial farm surrounds a new residential development.

This material is based upon work supported by the National Science Foundation through the Florida Coastal Everglades Long-Term Ecological Research program under Cooperative Agreements #DBI-0620409 and #DEB-9910514. This image is made available for non-commercial or educational use only.

Land-cover detection and landscape structure analysis in the Pachitea Basin, Peruvian Amazon

Classification procedures, including atmospheric correction satellite images as well as classification performance utilizing calibration and validation at different levels, have been investigated in the context of a coarse land-cover classification scheme for the Pachitea Basin. Two different correction methods were tested against no correction in terms of reflectance correction towards a common response for pseudo-invariant features (PIF). The accuracy of classifications derived from each of the three methods was then assessed in a discriminant analysis using crossvalidation at pixel, polygon, region, and image levels. Results indicate that only regression adjusted images using PIFs show no significant difference between images in any of the bands. A comparison of classifications at different levels suggests though that at pixel, polygon, and region levels the accuracy of the classifications do not significantly differ between corrected and uncorrected images. Spatial patterns of land-cover were analyzed in terms of colonization history, infrastructure, suitability of the land, and landownership. The actual use of the land is driven mainly by the ability to access the land and markets as is obvious in the distribution of land cover as a function of distance to rivers and roads. When considering all rivers and roads a threshold distance at which disproportional agro-pastoral land cover switches from over represented to under represented is at about 1km. Best land use suggestions seem not to affect the choice of land use. Differences in abundance of land cover between watersheds are more prevailing than differences between colonist and indigenous groups.

Half-degree gridded nitrogen and phosphorus fertilizer use for global agriculture production during 1900-2013

In addition to enhance agricultural productivity, synthetic nitrogen (N) and phosphorous (P) fertilizer application in croplands dramatically altered global nutrient budget, water quality, greenhouse gas balance, and their feedbacks to the climate system. However, due to the lack of geospatial fertilizer input data, current Earth system/land surface modeling studies have to ignore or use over-simplified data (e.g., static, spatially uniform fertilizer use) to characterize agricultural N and P input over decadal or century-long period. We therefore develop a global time-series gridded data of annual synthetic N and P fertilizer use rate in croplands, matched with HYDE 3,2 historical land use maps, at a resolution of 0.5º latitude by longitude during 1900-2013. Our data indicate N and P fertilizer use rates increased by approximately 8 times and 3 times, respectively, since the year 1961, when IFA (International Fertilizer Industry Association) and FAO (Food and Agricultural Organization) survey of country-level fertilizer input were available. Considering cropland expansion, increase of total fertilizer consumption amount is even larger. Hotspots of agricultural N fertilizer use shifted from the U.S. and Western Europe in the 1960s to East Asia in the early 21st century. P fertilizer input show the similar pattern with additional hotspot in Brazil. We find a global increase of fertilizer N/P ratio by 0.8 g N/g P per decade (p< 0.05) during 1961-2013, which may have important global implication of human impacts on agroecosystem functions in the long run. Our data can serve as one of critical input drivers for regional and global assessment on agricultural productivity, crop yield, agriculture-derived greenhouse gas balance, global nutrient budget, land-to-aquatic nutrient loss, and ecosystem feedback to the climate system.

ELUM: a spatial modelling tool to predict soil greenhouse gas changes from land conversion to bioenergy in the UK

Acknowledgements This work is based on the Ecosystem Land Use Modelling & Soil Carbon GHG Flux Trial (ELUM) project, which was commissioned and funded by the Energy Technologies Institute (ETI). The authors are grateful to Niall McNamara (Centre for Ecology & Hydrology, Lancaster) for coordinating the project and to Dagmar Henner (University of Aberdeen) for project assistance. We are also grateful to staff at the ETI, particularly to Geraldine Newton-Cross, Geraint Evans and Hannah Evans for constructive advice and feedback, and to Jonathan Oxley for project support. The ELUM Software Package contains Ordnance Survey data © Crown copyright and database right 2012.

Cogestion des ressources naturelles : une approche structurale pour quantifier la contribution des réseaux d'acteurs à la résilience des systèmes socio-écologiques

Alors que les activités anthropiques font basculer de nombreux écosystèmes vers des régimes fonctionnels différents, la résilience des systèmes socio-écologiques devient un problème pressant. Des acteurs locaux, impliqués dans une grande diversité de groupes — allant d’initiatives locales et indépendantes à de grandes institutions formelles — peuvent agir sur ces questions en collaborant au développement, à la promotion ou à l’implantation de pratiques plus en accord avec ce que l’environnement peut fournir. De ces collaborations répétées émergent des réseaux complexes, et il a été montré que la topologie de ces réseaux peut améliorer la résilience des systèmes socio-écologiques (SSÉ) auxquels ils participent. La topologie des réseaux d’acteurs favorisant la résilience de leur SSÉ est caractérisée par une combinaison de plusieurs facteurs : la structure doit être modulaire afin d’aider les différents groupes à développer et proposer des solutions à la fois plus innovantes (en réduisant l’homogénéisation du réseau), et plus proches de leurs intérêts propres ; elle doit être bien connectée et facilement synchronisable afin de faciliter les consensus, d’augmenter le capital social, ainsi que la capacité d’apprentissage ; enfin, elle doit être robuste, afin d’éviter que les deux premières caractéristiques ne souffrent du retrait volontaire ou de la mise à l’écart de certains acteurs. Ces caractéristiques, qui sont relativement intuitives à la fois conceptuellement et dans leur application mathématique, sont souvent employées séparément pour analyser les qualités structurales de réseaux d’acteurs empiriques. Cependant, certaines sont, par nature, incompatibles entre elles. Par exemple, le degré de modularité d’un réseau ne peut pas augmenter au même rythme que sa connectivité, et cette dernière ne peut pas être améliorée tout en améliorant sa robustesse. Cet obstacle rend difficile la création d’une mesure globale, car le niveau auquel le réseau des acteurs contribue à améliorer la résilience du SSÉ ne peut pas être la simple addition des caractéristiques citées, mais plutôt le résultat d’un compromis subtil entre celles-ci. Le travail présenté ici a pour objectifs (1), d’explorer les compromis entre ces caractéristiques ; (2) de proposer une mesure du degré auquel un réseau empirique d’acteurs contribue à la résilience de son SSÉ ; et (3) d’analyser un réseau empirique à la lumière, entre autres, de ces qualités structurales. Cette thèse s’articule autour d’une introduction et de quatre chapitres numérotés de 2 à 5. Le chapitre 2 est une revue de la littérature sur la résilience des SSÉ. Il identifie une série de caractéristiques structurales (ainsi que les mesures de réseaux qui leur correspondent) liées à l’amélioration de la résilience dans les SSÉ. Le chapitre 3 est une étude de cas sur la péninsule d’Eyre, une région rurale d’Australie-Méridionale où l’occupation du sol, ainsi que les changements climatiques, contribuent à l’érosion de la biodiversité. Pour cette étude de cas, des travaux de terrain ont été effectués en 2010 et 2011 durant lesquels une série d’entrevues a permis de créer une liste des acteurs de la cogestion de la biodiversité sur la péninsule. Les données collectées ont été utilisées pour le développement d’un questionnaire en ligne permettant de documenter les interactions entre ces acteurs. Ces deux étapes ont permis la reconstitution d’un réseau pondéré et dirigé de 129 acteurs individuels et 1180 relations. Le chapitre 4 décrit une méthodologie pour mesurer le degré auquel un réseau d’acteurs participe à la résilience du SSÉ dans lequel il est inclus. La méthode s’articule en deux étapes : premièrement, un algorithme d’optimisation (recuit simulé) est utilisé pour fabriquer un archétype semi-aléatoire correspondant à un compromis entre des niveaux élevés de modularité, de connectivité et de robustesse. Deuxièmement, un réseau empirique (comme celui de la péninsule d’Eyre) est comparé au réseau archétypique par le biais d’une mesure de distance structurelle. Plus la distance est courte, et plus le réseau empirique est proche de sa configuration optimale. La cinquième et dernier chapitre est une amélioration de l’algorithme de recuit simulé utilisé dans le chapitre 4. Comme il est d’usage pour ce genre d’algorithmes, le recuit simulé utilisé projetait les dimensions du problème multiobjectif dans une seule dimension (sous la forme d’une moyenne pondérée). Si cette technique donne de très bons résultats ponctuellement, elle n’autorise la production que d’une seule solution parmi la multitude de compromis possibles entre les différents objectifs. Afin de mieux explorer ces compromis, nous proposons un algorithme de recuit simulé multiobjectifs qui, plutôt que d’optimiser une seule solution, optimise une surface multidimensionnelle de solutions. Cette étude, qui se concentre sur la partie sociale des systèmes socio-écologiques, améliore notre compréhension des structures actorielles qui contribuent à la résilience des SSÉ. Elle montre que si certaines caractéristiques profitables à la résilience sont incompatibles (modularité et connectivité, ou — dans une moindre mesure — connectivité et robustesse), d’autres sont plus facilement conciliables (connectivité et synchronisabilité, ou — dans une moindre mesure — modularité et robustesse). Elle fournit également une méthode intuitive pour mesurer quantitativement des réseaux d’acteurs empiriques, et ouvre ainsi la voie vers, par exemple, des comparaisons d’études de cas, ou des suivis — dans le temps — de réseaux d’acteurs. De plus, cette thèse inclut une étude de cas qui fait la lumière sur l’importance de certains groupes institutionnels pour la coordination des collaborations et des échanges de connaissances entre des acteurs aux intérêts potentiellement divergents.

The relationship between vegetation and modern pollen assemblages on Mount Paggeo (NE Greece)

In this study, we investigated the relationship between vegetation and modern-pollen rain along the elevational gradient of Mount Paggeo. We apply multivariate data analysis to assess the relationship between vegetation and modern-pollen rain and quantify the representativeness of forest zones. This study represents the first statistical analysis of pollen-vegetation relationship along an elevational gradient in Greece. Hence, this paper improves confidence in interpretation of palynological records from north-eastern Greece and may refine past climate reconstructions for a more accurate comparison of data and modelling. Numerical classification and ordination were performed on pollen data to assess differences among plant communities that beech (Fagus sylvatica) dominates or co-dominates. The results show a strong relationship between altitude, arboreal cover, human impact and variations in pollen and nonpollen palynomorph taxa percentages.

Setting priorities for private land conservation in fire-prone landscapes: Are fire risk reduction and biodiversity conservation competing or compatible objectives?

Although wildfire plays an important role in maintaining biodiversity in many ecosystems, fire management to protect human assets is often carried out by different agencies than those tasked for conserving biodiversity. In fact, fire risk reduction and biodiversity conservation are often viewed as competing objectives. Here we explored the role of management through private land conservation and asked whether we could identify private land acquisition strategies that fulfill the mutual objectives of biodiversity conservation and fire risk reduction, or whether the maximization of one objective comes at a detriment to the other. Using a fixed budget and number of homes slated for development, we simulated 20 years of housing growth under alternative conservation selection strategies, and then projected the mean risk of fires destroying structures and the area and configuration of important habitat types in San Diego County, California, USA. We found clear differences in both fire risk projections and biodiversity impacts based on the way conservation lands are prioritized for selection, but these differences were split between two distinct groupings. If no conservation lands were purchased, or if purchases were prioritized based on cost or likelihood of development, both the projected fire risk and biodiversity impacts were much higher than if conservation lands were purchased in areas with high fire hazard or high species richness. Thus, conserving land focused on either of the two objectives resulted in nearly equivalent mutual benefits for both. These benefits not only resulted from preventing development in sensitive areas, but they were also due to the different housing patterns and arrangements that occurred as development was displaced from those areas. Although biodiversity conflicts may still arise using other fire management strategies, this study shows that mutual objectives can be attained through land-use planning in this region. These results likely generalize to any place where high species richness overlaps with hazardous wildland vegetation.

Perpetual private land conservation: the case for outdoor recreation and functional leisure

As natural areas, agricultural lands, and open spaces continue to be developed at unprecedented rates, it is important for land conservation professionals to understand the individuals who might play a role in permanently protecting these lands and their ecological services. Many factors have been shown to influence land protection decisions among private owners, including land-use activities, demographic characteristics, and environmental intention and behavior. With the hypothesis that individuals already involved in land conservation programs would be candidates for permanent protection, we set out to model conservation easement decisions within a group of participants in southern Indiana’s Classified Forest and Wildlands Program (ICFWP). We used a mailed questionnaire to survey 500 landowners, garnering 308 responses, about their interest in conservation easements. Our results indicated significant positive relationships between interest in conservation easements with variables representing perception of landscape change, outdoor recreation behavior as an adult, and environmental organization membership. By better understanding the ways these factors promote permanent land-use decisions, land conservation professionals can better allocate limited resources through strategic investments in targeting and outreach.

Sediment Budget for Cane Land on the Lower Herbert River Floodplain, North Queensland, Australia

Land use in the river catchments of tropical North Queensland appears to have increased the export of sediment and nutrients to the coast. Although evidence of harmful effect of sediment on coastal and riverine ecosystems is limited, there is a growing concern about its possible negative impacts. Sugarcane cultivation on the floodplains of the tropical North Queensland river catchments is thought to be an important source of excess sediment in the river drainage systems. Minimum-tillage, trash blanket harvesting has been shown to reduce erosion from sloping sugarcane fields, but in the strongly modified floodplain landscape other elements (e.g. drains, water furrows and headlands) could still be important sediment sources. The main objectives of this thesis are to quantify the amount of sediment coming from low-lying cane land and identify the important sediment sources in the landscape. The results of this thesis enable sugarcane farmers to take targeted measures for further reduction of the export of sediment and nutrients. Sediment budgets provide a useful approach to identify and quantify potential sediment sources. For this study a sediment budget is calculated for a part of the Ripple Creek catchment, which is a sub-catchment of the Lower Herbert River. The input of sediment from all potential sources in cane land and the storage of sediment within the catchment have been quantified and compared with the output of sediment from the catchment. Input from, and storage on headlands, main drains, minor drains and water furrows, was estimated from erosion pin and surface profile measurements. Input from forested upland, input from fields and the output at the outlet of the catchment was estimated with discharge data from gauged streams and flumes. Data for the sediment budget were collected during two ‘wet’-seasons: 1999-2000 and 2000-2001. The results of the sediment budget indicate that this tropical floodplain area is a net source of sediment. Plant cane fields, which do not have a protective trash cover, were the largest net source of sediment during the 1999-2000 season. Sediment input from water furrows was higher, but there was also considerable storage of sediment in this landscape element. Headlands tend to act as sinks. The source or sink function of drains is less clear, but seems to depend on their shape and vegetation cover. An important problem in this study is the high uncertainty in the estimates of the sediment budget components and is, for example, likely to be the cause of the imbalance in the sediment budget. High uncertainties have particularly affected the results from the 20002001 season. The main source of uncertainty is spatial variation in the erosion and deposition processes. Uncertainty has to be taken into consideration when interpreting the budget results. The observation of a floodplain as sediment source contradicts the general understanding that floodplains are areas of sediment storage within river catchments. A second objective of this thesis was therefore to provide an answer to the question: how can floodplains in the tropical North Queensland catchments can be a source of sediment? In geomorphic literature various factors have been pointed out, that could control floodplain erosion processes. However, their importance is not 'uniquely identified'. Among the most apparent factors are the stream power of the floodwater and the resistance of the floodplain surface both through its sedimentary composition and the vegetation cover. If the cultivated floodplains of the North Queensland catchments are considered in the light of these factors, there is a justified reason to expect them to be a sediment source. Cultivation has lowered the resistance of their surface; increased drainage has increased the drainage velocity and flood control structures have altered flooding patterns. For the Ripple Creek floodplain four qualitative scenarios have been developed that describe erosion and deposition under different flow conditions. Two of these scenarios were experienced during the budget study, involving runoff from local hillslopes and heavy rainfall, which caused floodplain erosion. In the longer term larger flood events, involving floodwater from the Herbert River, may lead to different erosion and deposition processes. The present study has shown that the tropical floodplain of the Herbert River catchment can be a source of sediment under particular flow conditions. It has also shown which elements in the sugarcane landscape are the most important sediment sources under these conditions. This understanding will enable sugarcane farmers to further reduce sediment export from cane land and prevent the negative impact this may have on the North Queensland coastal ecosystems.

Forests as landscapes of social inequality: tropical forest cover and land distribution among shifting cultivators

Can social inequality be seen imprinted in a forest landscape? We studied the relationship between land holding, land use, and inequality in a peasant community in the Peruvian Amazon where farmers practice swidden-fallow cultivation. Longitudinal data on land holding, land use, and land cover were gathered through field-level surveys (n = 316) and household interviews (n = 51) in 1994/1995 and 2007. Forest cover change between 1965 and 2007 was documented through interpretation of air photos and satellite imagery. We introduce the concept of “land use inequality” to capture differences across households in the distribution of forest fallowing and orchard raising as key land uses that affect household welfare and the sustainability of swidden-fallow agriculture. We find that land holding, land use, and forest cover distribution are correlated and that the forest today reflects social inequality a decade prior. Although initially land-poor households may catch up in terms of land holdings, their use and land cover remain impoverished. Differential land use investment through time links social inequality and forest cover. Implications are discussed for the study of forests as landscapes of inequality, the relationship between social inequality and forest composition, and the forest-poverty nexus.

Environmental impacts of freshwater use and suspended solids in life cycle assessment

Perturbation of natural ecosystems, namely by increasing freshwater use and its degradative use, as well as topsoil erosion by water of land-use production systems, have been emerging as topics of high environmental concern. Freshwater use has become a focus of attention in the last few years for all stakeholders involved in the production of goods, mainly agro-industrial and forest-based products, which are freshwater-intensive consumers, requiring large inputs of green and blue water. This thesis presents a global review on the available Water Footprint Assessment and Life Cycle Assessment (LCA)-based methods for measuring and assessing the environmental relevance of freshwater resources use, based on a life cycle perspective. Using some of the available midpoint LCA-based methods, the freshwater use-related impacts of a Portuguese wine (white ‘vinho verde’) were assessed. However, the relevance of environmental green water has been neglected because of the absence of a comprehensive impact assessment method associated with green water flows. To overcome this constraint, this thesis helps to improve and enhance the LCA-based methods by providing a midpoint and spatially explicit Life Cycle Impact Assessment (LCIA) method for assessing impacts on terrestrial green water flow and addressing reductions in surface blue water production caused by reductions in surface runoff due to land-use production systems. The applicability of the proposed method is illustrated by a case study on Eucalyptus globulus conducted in Portugal, as the growth of short rotation forestry is largely dependent on local precipitation. Topsoil erosion by water has been characterised as one of the most upsetting problems for rivers. Because of this, this thesis also focuses on the ecosystem impacts caused by suspended solids (SS) from topsoil erosion that reach freshwater systems. A framework to conduct a spatially distributed SS delivery to freshwater streams and a fate and effect LCIA method to derive site-specific characterisation factors (CFs) for endpoint damage on aquatic ecosystem diversity, namely on algae, macrophyte, and macroinvertebrates organisms, were developed. The applicability of this framework, combined with the derived site-specific CFs, is shown by conducting a case study on E. globulus stands located in Portugal as an example of a land use based system. A spatially explicit LCA assessment was shown to be necessary, since the impacts associated with both green water flows and SS vary greatly as a function of spatial location.

Aspects of Land Take in the Metropolitan Area of Naples

Land take is a phenomenon of great concern nowadays because of the large number of its negative impacts regarding biological, economic and social balance. In Italy, the development of urban and other artificial land has been irreversibly transforming a nonrenewable resource such as soil, regardless the almost constant population rate, with different speed depending of the region considered. The aim of this paper is to analyze the phenomenon in the metropolitan area of Naples, which is an area highly affected by territorial aggression of human matrix. The data used are both by the Institute for Environmental Protection and Research (ISPRA) Report 2015 on the usage of the land and by ISTAT relating to the resident population up to the 1st of January 2015 and the extension of land for agricultural use (Census 2010). The mathematical combination of this data creates a new indicator that can be referred to as “residual land”; this residual area is of great extension with many different characteristics and it could represent the area where the phenomenon of land take most occurs. The identification, measurement and analysis of “residual land” provide new insights on the evolution of land take and this new indicator can represent a critical element to work on to prevent future land transformation and protect natural and agricultural areas within the Italian context.

Resource use and greenhouse gas emissions from grain-finishing beef cattle in seven Australian feedlots: a life cycle assessment

Grain finishing of cattle has become increasingly common in Australia over the past 30 years. However, interest in the associated environmental impacts and resource use is increasing and requires detailed analysis. In this study we conducted a life cycle assessment (LCA) to investigate impacts of the grain-finishing stage for cattle in seven feedlots in eastern Australia, with a particular focus on the feedlot stage, including the impacts from producing the ration, feedlot operations, transport, and livestock emissions while cattle are in the feedlot (gate-to-gate). The functional unit was 1 kg of liveweight gain (LWG) for the feedlot stage and results are included for the full supply chain (cradle-to-gate), reported per kilogram of liveweight (LW) at the point of slaughter. Three classes of cattle produced for different markets were studied: short-fed domestic market (55–80 days on feed), mid-fed export (108–164 days on feed) and long-fed export (>300 days on feed). In the feedlot stage, mean fresh water consumption was found to vary from 171.9 to 672.6 L/kg LWG and mean stress-weighted water use ranged from 100.9 to 193.2 water stress index eq. L/kg LWG. Irrigation contributed 57–91% of total fresh water consumption with differences mainly related to the availability of irrigation water near the feedlot and the use of irrigated feed inputs in rations. Mean fossil energy demand ranged from 16.5 to 34.2 MJ lower heating values/kg LWG and arable land occupation from 18.7 to 40.5 m2/kg LWG in the feedlot stage. Mean greenhouse gas (GHG) emissions in the feedlot stage ranged from 4.6 to 9.5 kg CO2-e/kg LWG (excluding land use and direct land-use change emissions). Emissions were dominated by enteric methane and contributions from the production, transport and milling of feed inputs. Linear regression analysis showed that the feed conversion ratio was able to explain >86% of the variation in GHG intensity and energy demand. The feedlot stage contributed between 26% and 44% of total slaughter weight for the classes of cattle fed, whereas the contribution of this phase to resource use varied from 4% to 96% showing impacts from the finishing phase varied considerably, compared with the breeding and backgrounding. GHG emissions and total land occupation per kilogram of LWG during the grain finishing phase were lower than emissions from breeding and backgrounding, resulting in lower life-time emissions for grain-finished cattle compared with grass finishing.

Assessment of grazing selectivity in relation to land cover: a comparison between sheep and goats

Portugal’s Northeast production of sheep and goats are based on the exploitation of landscape by-products such as spontaneous native vegetation and agriculture leftovers. Shepherds tend the flocks throughout grazing itineraries every day, crossing a mosaic of patches of varied land uses. During the journey, the shepherd acts together with the sheep and goats to select each patch in creating an ordered sequence of land uses. The focus of the research is on the land-use composition of the grazing itineraries; determinate how they depend on the patterns of the landscape mosaic. It is utilized a data set of 26 monthly herd’s itineraries, 13 of sheep and 13 of goats, to investigate the relationship of the land uses crossed by the flocks and the land uses of the landscape, evaluating the land-use preferences and selectivity of the sheep and goats. It is utilized the divergences in the time spent and distance travelled by the herds and the area of the land uses in the landscape, the chi-square test to relate the preferred land used and the season, and the discriminate analysis to distinguish the preferences and the selectivity of the herd of sheep and the herd of goats. The herds of the sheep and the goats presented different land-use preferences over the seasons and the discriminant analysis shows that they have different landscape preferences. The herd of sheep has the highest selectivity indexes for the annual irrigated crops, the agricultural complex systems and the agroforestry land uses. The highest selectivity indexes for the herd of goats were found for the deciduous forest, the agriculture with natural and semi-natural spaces and the shrublands land uses. It was concluded that the landscape management for sheep and goats herding has to be different: the agricultural land uses are essential to the flocks of sheep and the forest land uses are decisive to the flocks of goats.