Review Special Issues

Quantitative characterization of animal social organization: Applications for epidemiological modelling

  • Social organization is a key aspect of animal ecology, closely interlinked with all aspects of animal behaviour. The structure of animal assemblages is highly diverse, both within and between species. The complexity and variety of social systems and the dynamic nature of interactions and dependencies between members of social groups have long been major obstacles for developing operational characterizations of social organization. Here, social network analysis, a set of statistical tools rooted in graph theory, suggests itself as a potential solution for this problem, by offering quantitative measures for various aspects of social relationships. In this review I will first introduce network analysis as a tool to characterize the social organization of animal groups and population and, then, focus on the application of this method for epidemiological modelling, specifically the prediction of spreading patterns of pathogens in livestock and its potential for informing targeted surveillance and planning of intervention measures.

    Citation: Bernhard Voelkl. Quantitative characterization of animal social organization: Applications for epidemiological modelling[J]. Mathematical Biosciences and Engineering, 2020, 17(5): 5005-5026. doi: 10.3934/mbe.2020271

    Related Papers:

    [1] Robinson J. Tanyi, Muyiwa S Adaramola . Bioenergy potential of agricultural crop residues and municipal solid waste in Cameroon. AIMS Energy, 2023, 11(1): 31-46. doi: 10.3934/energy.2023002
    [2] Anil Markandya, Kishore Dhavala, Alessandro Palma . The role of flexible biofuel policies in meeting biofuel mandates. AIMS Energy, 2018, 6(3): 530-550. doi: 10.3934/energy.2018.3.530
    [3] Michael E. Salassi, Lawrence L. Falconer, Tyler B. Mark, Michael A. Deliberto, Brian M. Hilbun, Todd L. Cooper . Economic Potential for Energy Cane Production as a Cellulosic Biofuel Feedstock in the Southeastern United States. AIMS Energy, 2015, 3(1): 25-40. doi: 10.3934/energy.2015.1.25
    [4] Pallav Purohit, Subash Dhar . Lignocellulosic biofuels in India: current perspectives, potential issues and future prospects. AIMS Energy, 2018, 6(3): 453-486. doi: 10.3934/energy.2018.3.453
    [5] Kharisma Bani Adam, Jangkung Raharjo, Desri Kristina Silalahi, Bandiyah Sri Aprilia, IGPO Indra Wijaya . Integrative analysis of diverse hybrid power systems for sustainable energy in underdeveloped regions: A case study in Indonesia. AIMS Energy, 2024, 12(1): 304-320. doi: 10.3934/energy.2024015
    [6] Jonathan B. Norman, Marcelle McManus . The techno-economic viability of bio-synthetic natural gas production utilising willow grown on contaminated land. AIMS Energy, 2019, 7(3): 285-312. doi: 10.3934/energy.2019.3.285
    [7] Romain Akpahou, Marshet M. Admas, Muyiwa S Adaramola . Evaluation of a bioenergy resource of agricultural residues and municipal solid wastes in Benin. AIMS Energy, 2024, 12(1): 167-189. doi: 10.3934/energy.2024008
    [8] Muhammad Rashed Al Mamun, Shuichi Torii . Comparative Studies on Methane Upgradation of Biogas by Removing of Contaminant Gases Using Combined Chemical Methods. AIMS Energy, 2015, 3(3): 255-266. doi: 10.3934/energy.2015.3.255
    [9] Muhammad Rashed Al Mamun, Anika Tasnim, Shahidul Bashar, Md. Jasim Uddin . Potentiality of biomethane production from slaughtered rumen digesta for reduction of environmental pollution. AIMS Energy, 2018, 6(5): 658-672. doi: 10.3934/energy.2018.5.658
    [10] Sebastian Auburger, Richard Wüstholz, Eckart Petig, Enno Bahrs . Biogas digestate and its economic impact on farms and biogas plants according to the upper limit for nitrogen spreading—the case of nutrient-burdened areas in north-west Germany. AIMS Energy, 2015, 3(4): 740-759. doi: 10.3934/energy.2015.4.740
  • Social organization is a key aspect of animal ecology, closely interlinked with all aspects of animal behaviour. The structure of animal assemblages is highly diverse, both within and between species. The complexity and variety of social systems and the dynamic nature of interactions and dependencies between members of social groups have long been major obstacles for developing operational characterizations of social organization. Here, social network analysis, a set of statistical tools rooted in graph theory, suggests itself as a potential solution for this problem, by offering quantitative measures for various aspects of social relationships. In this review I will first introduce network analysis as a tool to characterize the social organization of animal groups and population and, then, focus on the application of this method for epidemiological modelling, specifically the prediction of spreading patterns of pathogens in livestock and its potential for informing targeted surveillance and planning of intervention measures.


    1. Introduction

    According to the EU Renewable Energy Directive (RED), 10% of the fuel used within the transportation sector in the EU should be renewable in 2020, and by 2015, a share of 6.7% was reached [1,2]. The vast majority of these fuels are so called first generations biofuels produced from crops, with biodiesel dominating at 79% [3]. However, the EU's biofuels policy, where the primary goal is reduction of greenhouse gas (GHG) emissions, is much debated. A major criticism is that the use of crop based biofuels causes displacement of food crop production and indirect land use change (iLUC) with impacts on GHG emissions [4]. This discussion has led to an amendment of the EU RED in 2015, with a cap of 7% on the use of crop based biofuels by 2017 and a higher demand on reduction of GHG emissions (60%) for production plants taken in operation after the end of 2015 [5]. In this so called iLUC-directive, iLUC-factors are also introduced, a GHG emission penalty to be included in the reporting for biofuels based on cereal and other starch-rich crops, sugars and oil crops. Together with a higher demand on GHG reduction, some crop based biofuels might not meet the amended sustainability criteria, which might be the reason why a decreased consumption of biodiesel made the EU biofuel consumption as a whole drop with 2% between 2014 and 2015 [3].

    Sweden stands out as having the highest share of renewable energy in domestic transport in the EU, with 24% in 2015 according to the calculation rules in the EU RED [2]. The Swedish biofuels supply is partly based on crops, accounting for 40%, but also on biofuels produced from waste and residues [6]. A unique feature of the Swedish biofuel supply is the use of biogas for transport, representing 9% of the biofuels. Biogas can be produced from various types of feedstock such as sewage sludge, municipal and industrial waste, manure, crop residues and dedicated crops. Biogas can also be utilized for different purposes such as heat and power production, as a vehicle fuel, or as a renewable alternative to natural gas in other applications. Which feedstock and utilization pathways that are applied are highly dependent of which policy instruments that are implemented, as demonstrated by a comparison between Sweden and Germany. In Sweden, 63% of the current biogas production is utilized as vehicle fuel, and 95% of the feedstock is municipal and industrial waste and residues [7]. In Germany, representing 50% of the total biogas production in the EU, other policy instruments are applied which results in a completely different biogas system. In fact, only 5% of the biogas is based on waste while the rest is based on crops (52%) and manure (43%) [8,9]. The biogas produced is utilized for heat and power and less than 1% is used as vehicle fuel [9].

    To reach the national goal of 70% reduction on GHG emissions in the Swedish transportation sector by 2030, a demand of domestic production of biofuels based on agricultural biomass is foreseen [10,11,12,13]. Biofuel production could be based upon a wide range of different crops, which may be more or less suitable depending on which aspects, such as e.g. land use, energy, cost or GHG efficiency, that is prioritized [14,15]. Current EU policy does however prioritize so called advanced biofuels, which are produced from feedstock that does not compete with food/feed crops, such as waste and agricultural residues.

    Although agricultural residues represent a large potential, in Southern Sweden dominated by straw, manure and sugar beet tops, only limited amounts are utilized today. This could indicate that biofuels from such feedstock is non-competitive compared both to fossil fuels and to biofuels produced from crops and waste under existing policy framework. Thus, if an increased production of advanced biofuels from agricultural residues is desired, as described in the iLUC-directive, it is important to evaluate the cost aspects under similar and local conditions. Such evaluation could define whether the use of agricultural residues would increase the production cost as compared to biofuels based on crops, and if so, which systems that should/should not be promoted through policy instruments.

    The study presented in this paper is related to previously published work with the objective to establish and compare land use efficiency, energy efficiency, economic efficiency and GHG efficiency for biogas derived from agricultural feedstock from a coastal agricultural area situated in southern Sweden. In previous publications, energy input in cultivation, crop yields, and overall energy and greenhouse gas efficiency from field to fuel have been presented [14,15].

    In this study, the production costs of biogas of vehicle fuel quality was compared for 8 agricultural crops or residues that are or could potentially be used for biogas production in the investigated region, and which are differently addressed in the iLUC-directive [5];

    (1) Four food crops that are mentioned in the iLUC-directive as to be phased out as biofuel feedstock after 2020, and where iLUC-factors will be added in reporting of GHG-reduction; wheat grain, sugar beet, maize and Triticale. The latter is a high yielding cereal crop that in this study is harvested green as whole crop silage.

    (2) Grass, which has an unclear status at the moment. Grass is not burdened with an iLUC-factor in the iLUC directive, but grass cultivated on arable land for biofuel purposes could potentially fall under the definition as a "crop grown as main crop primarily for energy purposes on agricultural land", and thus according to the iLUC-directive support for this should be phased out. However, "grassy energy crops" is also listed as a feedstock for the production of advanced biofuels, which is encouraged.

    (3) Three agricultural residues: manure, wheat straw and sugar beet tops which are only used in minor amounts for biofuel production in Sweden today. These three are also listed as feedstock for the production of advanced biofuels, and are abundant in the region investigated.

    The overall objective is to evaluate how the biogas production costs depends on feedstock, and which aspects of the production costs that dominate for different feedstocks. The production costs are compared to current vehicle fuel prices in Sweden. The analysis also includes the impact of some current policy instruments and the potential outcome if such instruments were implemented for biogas from residues only and not for biogas from crops.

    It should be noted that prices might change quickly and vary largely depending on several factors. Therefore, the values presented here should not be considered as absolute and care should be taken when comparing the results with data from other sources.


    2. Material and Methods


    2.1. System overview

    Eight different crops and agricultural residues are evaluated as biofuel feedstock from an economic point of view. The evaluation includes the production of biogas, upgrading to vehicle fuel quality, distribution and filling stations. In addition, storage, distribution and application of digestate, the liquid residue after biogas production, on arable land is included.

    Feedstock cost as well as the impact of different features of crops and residues results in different cost to produce and distribute the same amount of biogas. To demonstrate the feedstock-related differences, calculations are based on a process operating on one single feedstock, even if neither biogas production, nor in some cases storage, would be practical for single feedstock.

    All gas volumes are given as dry gas at 0 ℃ and 101 kPa, and the lower heating value for methane is used in the conversion to energy units.


    2.2. Feedstock

    Feedstock properties after storage losses are given in Table 1. Background data, assumptions and references can be found in Appendix A1.

    Table 1. Feedstock properties.
    Wheat grain Beet Maize Triticale Grass Wheat straw Beet tops Manure1
    Dry matter (DM) (%) 86.0 21.5 25.9 35.3 17.1 90.5 12.7 12.0
    Volatile solids (VS) (%) 83.4 20.9 24.8 33.4 15.5 85.2 10.5 9.6
    CH4 (m3/MgDM) 369 372 307 341 243 224 264 170
    Macro and micro nutrients (mg/kg wet weight)
    N 18,000 1588 3499 4437 5654 4860 4213 4536
    P 4452 234 465 581 633 273 396 873
    K 4699 929 1222 2632 3280 9018 2040 4273
    S 1330 83 201 266 358 1140 330 309
    Fe 24.4 14.0 16.2 12.23 21.07 NA 69.9 53.2
    Mo 0.46 0.01 0.06 0.17 0.19 NA 0.06 0.20
    Co NA 0.02 0.01 0.04 0.14 NA 0.07 0.15
    Ni NA 0.11 0.15 0.08 0.33 NA 0.17 0.18
    1 Mix of liquid and solid manure.
     | Show Table
    DownLoad: CSV

    2.3. Biogas process

    The calculations are based on continuous stirred tank reactor (CSTR) processes, operated under mesophilic conditions (37 ℃), which is the dominating solution in German, full-scale, crop-based biogas plants [16]. The reactors are operated as one-stage processes.

    The parameters selected to define the CSTR process are the maximum organic loading rate (OLR), the minimum hydraulic retention time (HRT) and the maximum DM-concentration in the reactor. The output is the active reactor volume. Mineral nutrients of P, S, Fe, Co, Mo, Ni and Se are added if the feedstock provides insufficient amounts to support 5% microbial biomass formation. In additions, a minimum and maximum concentration of total ammoniacal nitrogen (TAN) are set to provide buffering capacity and avoid ammonia inhibition. N is added as carbonic diamide if the concentration is too low. Defining parameters are summarized in Table 2. More details about selected limits and need for mineral nutrients are found in Appendix A2.

    Table 2. Defining parameters for the biogas process.
    Parameter Unit Lower threshold Higher threshold
    Organic loading rate kgDM-rem/(m3*d) - 3
    Hydraulic retention time d 46 -
    DM concentration in reactor % - 10
    TAN in reactor g/L 1.8 4.0
     | Show Table
    DownLoad: CSV

    2.4. Storage and pre-treatment of feedstock

    Maize, Triticale, grass and beet tops are assumed to be harvested using a forage harvester, transported directly to the biogas plant and ensiled in bunker silos at site. Beet and dried wheat grain is stored off site and transported just in time to the biogas plant. Wheat straw is transported and stored as high-density bales. Liquid and solid manure, delivered just in time with trucks, is handled in a closed building with controlled ventilation to avoid odor disturbances in the neighborhood.

    For maize, Triticale, grass and beet tops, no additional pre-treatment is included. Beets and grain requires a minor pre-treatment (slicing and crushing respectively). Regarding straw, the feedstock is pre-treated to make the organic material more available for the microorganisms. Although different kinds of pre-treatment are presented in the literature, economic calculations in this study are based on mechanical methods (extruder and hammer mill) used at Danish biogas plants [17,18].

    Finally, the use of liquid and solid manure must also fulfill the European directive for the treatment of animal byproducts [19]. To meet these requirements, a mixture of liquid and solid manure is pumped with cutting pumps (to ensure that the particle size is less than 12 mm) to hygienization tanks and heated to 70 ℃ for one hour. Thereafter, the manure is pumped to the reactors while crops and crop residues are fed into the biogas reactor with screw conveyors.


    2.5. Digestate

    The reactors are assumed to be ideally stirred, giving effluent concentrations equal to those in the reactor. The digestate is stored in covered storage tanks, with a loss of TAN amounting to 1% of N [20]. Storage capacity for crop digestate sufficient for 12-monthsx production is included. For manure, a digestate volume corresponding to the amount of liquid manure used as feedstock is assumed to be transported back to existing storage tanks at the farms. For the excess digestate, storage capacity for 12-monthsx production is included.


    2.6. Upgrading of biogas to vehicle fuel quality

    Biogas typically contains 50–60% CH4, 40–50% CO2 and some contaminants such as hydrogen sulphide, ammonia and particles. Sweden has implemented a standard for biogas utilized as vehicle fuel that requires the removal of contaminants and a methane concentration of 95–99% [21]. Currently, several technologies for the removal of CO2, also known as upgrading, are commercially available. In Sweden, having 59 biogas upgrading plants in December 2015, 71% of the plants use the water scrubber technology, 8% use pressure swing adsorption (PSA), 19% use chemical scrubbing and 2% use membrane technology [22]. All technologies have different characteristics considering e.g. process energy requirements and methane losses, see Table 3. A more detailed presentation of these and other technologies is e.g. found in Bauer et al. [23] and Hoyer et al. [24].

    Table 3. Process energy requirements per volume of biogas and methane losses for different upgrading technologies [9,24].
    Electricity (kWh/m3) Heat (kWh/m3) Methane losses (%)
    Amin scrubber 0.06–0.15 0.5–0.8 0.1
    PSA 0.20–0.25 - 1–5
    Water scrubber 0.20–0.30 - 0.5–2.0
     | Show Table
    DownLoad: CSV

    Different characteristics regarding for instance process energy requirements, methane losses and methane concentration in the upgraded gas affect the economic outcome. Thus, local conditions such as the price of electricity and heat, possibilities to utilize waste heat and the market price for upgraded biogas determine the most favorable technology in each case. Economy of scale, however, has a more general effect on the production cost. As presented in Hoyer et al. [24], there is a clear efficiency of scale concerning the investment cost for upgrading plants up to an installed capacity of approximately 1000 m3 biogas/h. Increasing the capacity further will have a minor impact on the investment cost per m3 of installed capacity. A similar efficiency of scale is found for the total upgrading cost [25,26,27]. Although larger plants result in lower upgrading cost per energy unit, a larger biogas plant may also results in a higher transportation cost for feedstock and digestate [28]. Thus, the analysis in this study is based on a biogas plant producing 1000 m3/h of biogas. With a methane content of 55%, this corresponds to an annual methane production of 172 TJ. Due to low losses of methane and a low consumption of electricity, calculations are based on a chemical scrubber.


    2.7. Distribution of biogas of vehicle fuel quality

    In Sweden, upgraded biogas is distributed to filling stations by truck (70%) or gas grid (30%) which could be a dedicated biogas grid or the natural gas grid [7]. When distributed by truck, upgraded biogas is compressed to 20–25 MPa and loaded in high-pressure vessels of steel or composite. Biogas can also be liquefied and distributed as LBG (liquefied biogas) although it is not a common solution in Sweden today [29].

    When distributed in gas grids, both local dedicated biogas grids as well as the national natural gas grid are utilized today. Most biogas plants inject gas into the low pressure grid (4 bar) but it is also possible to inject into the high pressure grid (80 bar) [29]. Currently, injection to the national grid normally requires that propane is added to adjust the heating value of the biogas to correspond to the distributed natural gas. In this study, calculations are made for biogas injected into the low-pressure natural gas grid. Based on the discussion presented in Börjesson et al. [29] propane addition is not included, assuming conditions as for a dedicated biogas grid.

    The final step in the chain of distribution is the methane filling station. In Sweden, 211 filling stations was in operation in the end of 2015. The average amount of gas distributed per station annually corresponded to almost 20 GWh for stations dedicated to busses and 4.5 GWh for public stations [29]. In this study, calculations were based on filling stations distributing 10 GWh annually.


    2.8. Energy input

    In the following sections, energy input is presented as data measured at the biogas plant, thus as final energy use. This is also used in the economic calculations.


    2.9. Biogas production

    Process energy input is affected by plant design and operational conditions as well as by characteristics of the feedstock. Process heat requirements depend among other things on process and feedstock temperature, reactor configuration, insulation and heat exchangers as well as on hygienization requirements. In the literature, there are several reports of internal heat demand for biogas plants utilizing different kinds of feedstock. A demand in the range of 70–130 MJ/t has been reported [16,20,30,31,32], but there are also examples where the reported heat demand is as high as 250–330 MJ/Mg [16,33,34]. For comparison, a theoretically calculated heat demand for the feedstocks in the current study (assuming an average feedstock temperature of 8 ℃, a process temperature of 37 ℃ and a specific heat capacity of water and dry matter of 4.18 and 1.0 kJ/kg/K respectively) would be 89–109 MJ/Mg. In addition, heat is required to compensate for losses through the reactor walls and via the biogas produced. These losses are affected by insulation, local weather conditions and reactor configuration etc. Here, the heat demand is increased with 15% to compensate for these losses [31] which gives a calculated heat demand of 102–126 MJ/Mg excluding heat exchangers. For comparison, a recent study of some Danish co-digestion plants found that the net process heat demand was 54 to 72 MJ/Mg of feedstock [17]. Given the numbers calculated here and the ones found in the literature, it is clear that process heat demand varies between different plants depending on feedstock as well as design and operation of the plant. In this study, the calculated heat demand is applied for all biogas plants but the one using manure, where the net heat demand is set to 70 MJ/Mg representing the higher end of the interval presented by Hjort-Gregersen [17].

    Electricity is used for feeding, stirring, pumping etc. and the demand is influenced by feedstock type and reactor design. In the literature, reported electricity input varies from 5–41 kWh/Mg [16,20,31,32,33,34,35]. However, when the biogas is utilized for CHP, it is common that reported figures also include parasitic load for the production of CHP, which could represent 9–40% of the total electricity consumption [35]. In a recent study of some Danish biogas plants, the electricity consumption for large scale co-digestion plants was found to be 6–7 kWh/Mg of feedstock [17]. In this study, the internal electricity demand for the biogas production process only is taken to be 8 kWh/Mg wet weight of added feedstock (including water where applicable).


    2.10. Upgrading

    As presented in Table 3, different upgrading technologies have different energy requirements.

    In this study, calculations are based on a chemical scrubber using 0.4 MJ electricity and 2.0 MJ heat per m3 of biogas.


    2.11. Utilization of waste heat

    Both the biogas process and the upgrading process require heat. However, the upgrading process also produces waste heat that can be utilized to heat the biogas rector. In total, 75–80% of the heat used in the upgrading process is available at a temperature of 55–60 ℃ with a return temperature of 35 ℃. It is also possible to split this heat stream in two, of which 1/3 is available at 85 ℃ and 2/3 is available at 55 ℃ (Karlsson, LE, Läckeby Water AB, personal communication). Thus, the possibility to use this waste heat depends on the temperature required at the biogas plant, where a mesophilic process requires a heat stream at 48 ℃ [36]. When the feedstock requires 70 ℃ hygienization, the temperature should be approximately 85 ℃ [37]. Thus, the waste heat from the upgrading process should be sufficient in both cases. In this study, it is assumed that all available waste heat could be utilized when energy crops are digested, although it is not certain that the biogas plant can reduce the temperature of the heating circuit to the extent required to achieve full utilization. Regarding manure, it is assumed that only the high temperature waste heat could be utilized due to the requirements for hygienization.


    2.12. Distribution and filling stations

    Upgraded biogas is assumed to be injected into the local low-pressure grid, which requires an increase of pressure from 0.1 MPa to 0.4 MPa. At the filling station, the pressure is increased to 20 MPa. Total electricity consumption is set to 0.03 kWh/kWh of upgraded biogas [29].


    2.13. Economic assessment

    Considered parameters are feedstock cost, methane potential as well as feedstock properties and how these affect the biogas process, the need for additives and the design of the biogas plant as such. The assessment also includes storage and handling of digestate as well as its potential value when used as a biofertilizer. Finally, the calculated production cost is compared to current market price of vehicle gas at public filling stations.


    2.14. Feedstock

    The production cost for crops applied in this study, presented in Table 4, includes crop production, harvest, transport to the biogas plant (30 km for grain and 7 km for the other crops) and storage. The production cost for crops does not include the cost for land nor any EU subsidies. Crop yield levels representative for the agricultural area (Gss) in southern Sweden are chosen based on the yields used in the assessments of GHG emissions and energy input that have been performed as part of the same research project [14,15,38]. The cost for agricultural residues includes harvest and handling, storage as well as a cost compensation for the removal of nutrients in crop residues as compared to the current practice of leaving the residues in the field. The cost also includes transportation to the biogas plant corresponding to 30 km (straw) and 7 km (beet tops). Regarding wheat straw, the recoverable amount in the region has been calculated to 2.8 MgDM/ha [39]. According to Valin et al. [40], not more than 33–50% should be removed not to jeopardize soil carbon development. The straw removal used in the present study is 1.6 MgDM/ha, which represents close to 60% of the total amount, following the assumptions in Rosenqvist [38]. For beet tops, calculations are based on a recovery of 3.0 MgDM/ha. For comparison, the mean yield in the county of Skåne has been shown to be 3.4 MgDM/ha, meaning that nearly 90% of the total yield is assumed to be recovered [41].

    Table 4. Production/recovery cost for crops and crop residues (€/MgDM) [38].
    Wheat grain Beet Maize Triticale Grass Wheat Straw Beet tops Manure
    Liquid Solid
    Yield (MgDM/ha) 6.5 11.2 10.0 7.5 10.0 1.6 3.0 - -
    Machinery 59 72 55 65 49 60 33 - -
    Mineral fertilizers 40 26 39 21 47 10 36 - -
    Transport 7 15 20 20 20 14 52 - -
    Other 51 68 46 40 32 6 8 24 15
    Total 156 180 160 146 149 91 129 24 15
     | Show Table
    DownLoad: CSV

    The calculated cost for crops and residues are based on data presented by Rosenqvist [38] representing price levels of 2016. The reason for using one reference instead of a review of multiple studies is to make sure that calculations are made according to the same methodology and using the same assumptions and background data. For comparison, the calculated cost of wheat grain corresponds to 95% of the Swedish market price for wheat grain in 2016 [42]. For straw, the calculated cost of 91 €/Mg DM corresponds to approximately 5 €/GJ (based on the lower heating value) which is 95% of the average market price for wood chips used for district heating in 2016 [43]. The calculated cost is also comparable of the cost for straw used for district heating in Denmark [44]. As mentioned above, the production cost presented in Table 4 is based on yields representative for southern Sweden, which are relatively high in a Swedish context. Assuming more or less suitable agricultural conditions could decrease the production cost up to 16% or increase it up to 54% for some crops or residues (data not shown here) [38]. Given these uncertainties, the impact of different feedstock costs is included in the sensitivity analysis.

    Regarding manure, it is assumed that manure is available at the farm free of charge for the biogas producer, resulting in a transport cost only. The transportation cost is calculated based on assumed loading capacity, time to load and unload, cleaning of the truck, average speed and distance. Based on the assumptions and method presented in Appendix A4, the average one-way distance is calculated to 24.7 and 22.7 km for liquid and solid manure respectively. The resulting transport cost corresponds to 4.3 and 2.9 €/Mg respectively. For liquid manure, the transportation cost also includes return transport of digestate.


    2.15. Pre-treatment and handling of feedstock

    Wheat grain, beet and straw are treated before fed into the reactor. In a recent study, the cost for pre-treatment of straw at Danish biogas plants was calculated to 17–43 €/Mg wet weight depending on technology, scale and level of utilization [17]. The average cost was 31 €/Mg which could be compared to the average cost of 10 €/Mg for pre-treatment of solid manure [17,18]. With a straw DM content of 90.5% (Table 1), the cost for pre-treatment of straw is set to 34 €/MgDM in this study. For further comparison, the cost for extrusion of grass has been calculated to 40 €/MgDM [45]. For wheat grain and beet, the cost for crushing and slicing is set to 0.9 €/Mg [38]. The time required for handling of feedstock at the biogas plant is set to 1.5 minutes/Mg for grain and straw, 3.8 minutes/Mg for maize, Triticale grass and beet tops and 2 minutes/Mg for solid manure [46]. The cost of personnel is set to 30 €/h.


    2.16. Additives

    The costs of additives in the biogas process are 1.1 €/kg (N) and 2.2 €/kg (P) based on mineral fertilizer costs [38]. For S and Fe, the cost are 3.0 €/kg and 4.1 €/kg respectively based on price estimates in April 2017 (Brenntag Nordic AB, Malmö; Kemira Kemi AB, Helsingborg; Yara AB, Köping). Micronutrients (Co, Mo, Ni and Se) are applied as individualised mixtures, for each feedstock (from Schaumann BioEnergy GmbH, Pinneberg) at a fixed cost of 21,600 € per biogas plant and year for 1600 kg of mixture, independent on variations of micronutrients in the mix.


    2.17. Investment costs

    The specific investment costs for biogas plants designed to process energy crops varies in the literature depending on scale, material, feeding technique, etc. Data presented by FNR [16] show that the investment cost for 41 biogas plants varies from 209–895 €/m3 active volume, including the CHP unit, which on average represents 22.5% of the total investments. However, when comparing biogas plants of similar design and features, efficiency of scale is clear, albeit limited for larger biogas plants [47]. For large-scale biogas plants, treating approximately 10,000–100,000 t of maize (90% of the feedstock) and manure (10% of the feedstock), the investment costs, not including production of CHP, vary from 82 to 114 €/t [25]. Assuming a VS content of 33.6% for maize and 6.4% for manure, investment costs were recalculated to €/m3 active reactor volume as shown in Figure 1, by assuming an OLR of 2.5 kgVS added/m3/d [25,47]. The investment cost for a biogas plant with an active reactor volume of 7000 m3 would with this method be calculated to 309 €/m3. For comparison, FNR [48] presents data for biogas plants with an active reactor volume of 7000 and 7200 m3. Not including the cost for biogas utilization and digestate storage, the investment is approximately 320 €/m3. Since the calculated investment cost, as presented in Figure 1 is only a few percentage lower, the equation for the curve shown in Figure 1 is applied for all plants using energy crops and crop residues.

    Figure 1. Calculated investment cost.

    For the biogas plant using manure as feedstock, there will be additional investments for odor control (building, ventilation, filters etc.) and hygienization. The investment cost would therefore be higher. For such biogas plants with a reactor volume of 13–16,000 m3, the investment cost has been calculated to 400–480 €/m3[49]. For comparison, the calculated investment cost for a biogas plant based on manure and with a reactor volume of 32000 m3 is approximately 470 €/m3 (Andersson, S. Biogas Sydöstra Skåne AB, personal communication). Given these numbers, the investment cost for large-scale manure based biogas plants could be 50–100% higher per m3 reactor than biogas plants using crops. In order not to underestimate the investment cost for manure based biogas plants, the investment is calculated with the same method as applied in Figure 1 and then doubled.

    The investment costs for covered digestate storage tanks are set to be 30 €/m3[50].

    For the amine scrubber, with an installed capacity of 1000 m3/h, the investment cost is set to 2.0 million € including a wood chips boiler [29].

    The investment in filling stations is estimated to 830,000 € for a station with an annual capacity of 10 GWh. Thus, the total investment in filling stations is set to 4.2 million € [29].


    2.18. Capital cost

    The weighted average cost of capital (WACC) is set to 6% and the depreciation period is set to 15 years for all installations [29].


    2.19. Process energy

    Process heat is produced on site in a wood chip boiler. The cost for wood chips is set to 21 €/MWh based on the average price paid by district heating plants in Sweden in 2014–2016 [51]. Assuming a boiler efficiency of 85% and a cost for operation and maintenance of 1.7 €/MWh [29] the cost for heat is calculated to 26 €/MWh excluding cost of capital (included in the investment cost). Based on the average electricity price in 2015 and 2016 for industries using 2–20 GWh annually, the cost for electricity used in the biogas system is set to 55 €/MWh [52].


    2.20. Operation and maintenance

    For large-scale biogas plants treating agricultural feedstock, such as energy crops and manure, the costs of operation and maintenance of the biogas plant reported in the literature vary between approximately 5–10% of the investment, corresponding to approximately 5–10 €/t [25,31,53] including personnel, insurance, service and repairs etc. In this study, labor costs for feedstock handling is calculated separately and the costs of operation and maintenance are set to 5% of the investment. For the upgrading, the cost for operation and maintenance, excluding process energy, is set to 3% of the investment [54]. Operation and maintenance of the digestate storages are set to 0.5% of the investment [50].


    2.21. Digestate

    The economic value of the digestate produced depends on the concentration of nutrients and the cost of alternative fertilizers available to the farmer. Different crops and soil type will also require different fertilization, which affects the economic value of the digestate. Based on discussions with industrial actors in Ireland, Smyth et al. [53] estimate the net value to be 4 €/Mg for digestate from the digestion of grass, not including transport and storage. According to Brown et al. [31], a conservative assumption is that the value of the digestate equals the cost of transport and spreading. In this study, the cost of storage and transport of the digestate is included in the biogas production cost.

    The transportation distance is calculated in the same ways as for manure, see Appendix A4, assuming that the biogas plant is in a center of a circle. The size of the circle is calculated based on the average share of arable land in the county of Skåne (41%) and assuming that 30% of the arable land surrounding the biogas plant is fertilized with digestate. The amount digestate spread per hectare will depend on P and TAN content where the maximum amounts added are limited to 22 and 150 kg/ha respectively to comply with Swedish legislation for organic fertilizers [55]. Thus, the transportation distance varies depending of the content of the digestate.

    In order to present the impact of the value of the digestate in the overall economic calculation, two values are applied. As a minimum value, it is assumed that the farmer is given the digestate free of charge. As a maximum value, the nutrients in the digestate (N, P and K) are valued as mineral fertilizers with a market price of 1.1, 2.4 and 0.8 €/kg, respectively [38]. This maximum value is, however, reduced by the cost of spreading the digestate. The cost of spreading liquid manure could vary from 1.7 to 3.8 €/Mg depending on how it is spread [56]. Based on this interval, the cost for spreading of digestate is set to 2 €/Mg. Regarding the digestate from manure, P and K are not affected by the biogas process, and digestion gives no added value. Organically bound nitrogen is, however, mineralized to TAN, so TAN concentration is higher in the digestate than in the untreated manure. The value of this additional TAN is used in the calculation of the maximum economic value. Also, the cost of spreading digestate is reduced by the cost that would have been generated by spreading the untreated manure, which is set to 2 €/Mg for liquid manure and 3.5 €/Mg solid manure [56]. Any additional investments that might be required at individual farms in order to deliver manure and receive digestate are not included.


    2.22. Distribution and filling stations

    The cost for gas grid distribution is set to 8.9 €/MWh and the cost for operation and maintenance of the filling stations is set to 5.6 €/MWh [29].


    2.23. Policy instruments

    The production and utilization of biogas is affected by a wide range of policy instruments [57]. In Sweden, there are energy and carbon dioxide taxes on fossil fuels from which renewable vehicle fuels such as biogas are exempted as long as the sustainability criteria in the EU RED are met and the biofuel is not overcompensated as compared to the market price of fossil fuels. However, these exemptions are only approved by the European Commission until 2018 for liquid biofuels and 2020 for biogas [58,59]. To meet EU regulation, a reduction quota system has been suggested for low-blend biofuels [60]. It is, however, not yet (June 2017) clear which policy instrument that will be applied for high blend or pure biofuels such as biogas.

    In addition to these tax exemptions, there is a production subsidy for biogas from manure which could be as high as 12 €/GJ. However, the subsidy can vary between years, since all producers must share the available funds, which are fixed until the year 2023. It is also possible to apply for different kinds of investment subsidies. However, according to EU regulation, none of these subsidies can be given to biogas plants using energy crops as feedstock [61,62,63]. The potential impact of these subsidies is demonstrated in the sensitivity analysis.

    Finally, the calculated biogas production cost is compared to current average market price of compressed biogas (CBG) and natural gas (CNG) as well as gasoline and diesel including the tax exemption applied until 2020. Since there are restrictions in EU regulations regarding the support of biofuels from food crops, the same comparison is also made with full taxes on biogas from these feedstock.


    2.24. Vehicle fuel market prices

    The market price of vehicle fuels more or less change on a daily basis. The prices presented in Table 5 reflect the situation in 2016, including current policy instruments, at public filling stations. Diesel and CBG/CNG have similar market prices and petrol is approximately 15% more expensive.

    Table 5. Market price of diesel, petrol and CBG excl. VAT (€/GJ) [6,29,64].
    Diesel Petrol CNG CBG
    Fuel price 12.9 15.3 24.2 30.9
    Energy and CO2-tax 17.0 20.5 6.7 -
    Total price 29.9 35.8 30.9 30.9
     | Show Table
    DownLoad: CSV

    3. Results and Discussion


    3.1. Feedstock amounts, reactor dimensions, additives and digestate content

    The feedstock requirements to reach a biogas production corresponding to 1000 m3/h, or 172 TJ/year are shown in Table 6. The result from the dimensioning of the biogas process will apart from the demand for feedstock be a requirement of active reactor volume and in some cases a need for water and additives. The required addition of nutrients reflects the low concentrations of some essential nutrients that are present in some of the crops. The required addition for individual micronutrients can be found in Appendix A3. Apart from biogas, the process will also generate digestate. Calculated amounts, nutrient content and transport distance are also given in Table 6.

    Table 6. Resulting process requirements, feedstock and transport demand.
    Wheat grain Beet Maize Triticale Grass Wheat Straw Beet tops Manure
    Process
    Feed-in, wet [Mg/a] 15,200 60,160 60,720 40,010 119,130 23,730 143,540 236,2001
    Feed-in, DM [Mg/a] 13,070 12,960 15,700 14,130 18,710 21,470 18,270 28,3402
    Water [Mg/a] 41,270 650 90,350
    Active reactor volume [m3] 10,200 10,200 10,200 10,200 15,000 14,400 18,100 29,800
    Additives [kg/a]
    Fe 450 130 70 410 NA
    N 71,340 190,810
    P 5440
    S 330
    Micronutrient mixture 1600 1600 1600 1600 1600 1600 1600 1600
    Digestate
    Amount (t/yr) 45,290 48,980 49,540 29,480 107,950 102,900 132,370 225,000
    N-tot (kg/t)3 6.0/6.0 2.0/3.4 4.3/4.3 6.0/6.0 5.5/5.5 1.1/3.0 4.6/4.5 4.8/4.7
    TAN (kg/t)3 4.0/3.9 0.3/1.8 1.9/1.8 2.8/2.8 3.2/3.1 0.0/1.8 3.0/3.0 3.1/3.1
    P (kg/t) 1.5 0.3 0.6 0.8 0.6 0.1 0.4 0.9
    K (kg/t) 1.6 1.1 1.5 3.6 3.3 2.1 2.2 4.5
    Field application [Mg/ha/a] 14.7 76.5 38.6 27.9 36.8 84.7 50.2 24.0
    One way transport distance [km] 7.7 3.5 5.0 4.5 7.6 4.9 7.2 6.64
    1 Whereof 64,410 t/a is solid manure
    2 Whereof 12,880 t DM/a from solid manure
    3 Nitrogen concentrations are given both as based on concentration in crop after processing and as concentration after addition (if applicable) and storage loss (1% of N-tot at storage under roof cover [65]).
    4 An amount of digestate corresponding to the delivery of liquid manure (171,770 t/a) is assumed to be delivered back to the farm with the return transport. This transport distance is for the additional 53,230 t/a of digestate for which new storage capacity is built and that is applied at fields in the vicinity of the biogas plant.
     | Show Table
    DownLoad: CSV

    3.2. Economy

    The calculated cost for production and distribution of biogas, as presented in Table 7, range from 24 to 33 €/GJ, including the full value of the NPK in the digestate and not including any subsidies. Thus, biogas from all feedstock but grass and straw could meet the average market price of 31 €/GJ CBG. Since this is the average market price and the calculated biogas cost is only slightly higher, it might be possible to find local markets where biogas from all feedstock could be competitive. The economical margins for grass, manure, straw and beet tops are however very small or non-existing. The feedstock cost is based on the amount of feedstock required to produce 172 TJ annually and includes the production cost presented in Table 4 and the methane yield given in Table 1. Due to the difference in methane yields, the feedstock cost per MJ of fuel for straw is i.e. only marginal lower than the cost for wheat grain even if the DM-based feedstock cost is much lover for straw.

    Table 7. Investment and production cost.
    Wheat grain Beet Maize Triticale Grass Straw Beet tops Manure
    Investment cost (k€)
    Biogas plant 2972 2971 2971 2972 4111 3964 4808 14,618
    Digestate storages 1359 1469 1486 884 3239 3087 3971 1597
    Upgrading plant 2000 2000 2000 2000 2000 2000 2000 2000
    Gas grid injection 100 100 100 100 100 100 100 100
    Filling stations 4167 4167 4167 4167 4167 4167 4167 4167
    Annual cost (k€)
    Feedstock 2045 2283 2513 2068 2787 1945 2351 10131
    Feedstock handling 11 120 121 80 238 18 287 64
    Feedstock pre-treatment 14 53 - - - 738 - 659
    Nutrients and water 45 103 23 24 22 301 22 22
    Capital 1091 1102 1104 1042 1402 1371 1549 2315
    Electricity 198 200 200 189 230 228 243 202
    Heat 125 125 125 125 137 130 162 125
    Gas grid distribution 425 425 425 425 425 425 425 425
    Operation and maintenance
    Biogas 149 149 149 149 206 198 240 731
    Digestate storage 7 7 7 4 16 15 20 8
    Upgrading 60 60 60 60 60 60 60 60
    Filling stations 265 265 265 265 265 265 265 265
    Digestate transport 71 60 67 39 168 138 201 -
    Total 4505 4953 5059 4471 5955 5832 5825 5883
    Annual revenue
    Digestate 329 76 129 170 594 192 541 396
    Production cost (€/GJ)2 24.3 28.4 28.7 25.0 31.2 32.8 30.7 31.9
    Production cost (€/GJ)3 26.2 28.8 29.4 26.0 34.6 34.6 33.9 34.2
    1 Transport cost only, including return transport of digestate
    2 Including the value of NPK in the digestate
    3 Excluding the value of NPK in the digestate
     | Show Table
    DownLoad: CSV

    The economic value of the digestate has different impact in different systems, see Figure 2. In general, the value of the digestate depends on the concentration of nutrients and the amount of nutrients in relation to the amount of biogas produced. For beet, maize and straw the value of the digestate reduces the overall cost with a few percent only, mainly due to the relatively low concentration of nutrients. Notice that the value of the digestate from straw is low even if considerable amounts of N and P has been added for process stability reasons. Without these additives, the digestate would not have any economic value. For grain, grass and beet tops, the value of the digestate reduces the cost with up to 10% due to the higher nutrient amounts. Digestate from Triticale has a high concentration of nutrients but because of a high methane potential the total amounts and thus the economic impact is low. Biogas from manure is a special case where not only the increased share of TAN adds to the economic value but also the avoided cost for spreading of solid manure.

    Figure 2. Production and distribution cost for biogas from crops and residues with and without the value of the digestate.

    3.3. Biogas from crops vs residues

    One of the main objectives with this paper is to analyze the difference between biogas systems based on three different kinds of feedstock (crops with an iLUC-factor; grass and residues). The assessment shows that biogas systems based on cereal crops (wheat grain and Triticale) results in the lowest cost while crop residues, grass, and manure gives the highest. Considering the overall cost for the biogas systems analyzed here, feedstock and feedstock pre-treatment represents approximately 40–50% of the total cost for production and distribution of upgraded biogas from crops and crop residues and less than 20% for systems based on manure.

    However, the overall cost presented in Table 7 include upgrading and distribution of biogas, which, with the system design applied here, is not affected by the choice of feedstock. Considering only costs related to biogas production and digestate handling, the feedstock cost represents an even higher share of the production cost, as seen in Figure 3. The choice of feedstock does however not only affect the feedstock cost as such. The features of the feedstock, such as biodegradability, content of micro-and macronutrients and water also affect the need for additives, reactor dimensions as well as storage, transport and economic value of the digestate. Compared to the cost for the feedstock (including pre-treatment), differences in the need for additives and reactor dimensions etcetera has, however, a relatively small impact. Thus, the feedstock with the lowest cost (including pre-treatment) compared to its methane potential is likely to result in the lowest production cost, see Figure 3. Of the analyzed feedstock, manure partly constitutes an exception to the latter two statements. The feedstock cost per produced methane is low, however, due to the low conversion degree to methane a large volume of inert material is treated, making the cost of capital and the cost for operation and maintenance high. This shows the dual importance of high conversion to methane.

    Figure 3. Production cost for biogas from crops and residues excluding upgrading and distribution.

    3.4. Uncertainties in calculated production costs

    The calculated biogas production cost depends on a number of estimations and assumptions. For the biogas systems analyzed in this paper, the production cost is dominated by the cost of feedstock, capital and in some cases pre-treatment, see Table 7 and Figure 3. For energy crops, the feedstock cost depends on local agricultural conditions, which may vary between different regions but also between different years. However, the feedstock cost is not only related to the actual production cost for the crop but also to the market value for the crop and the farmer's alternative use of land which may vary over time. For crop residues, the cost for handling and storage may also vary depending on local agricultural conditions and transport distance to the biogas plant. For manure, the feedstock cost applied in this study consists of transportation costs only. Thus, changing the transportation distance, the average speed or the loading capacity would have a direct impact on the feedstock cost. If there is need for investments in farm installations, this would also increase the cost. Regarding the cost of capital, this is affected by the investment cost as such, which depends on the specific plant design as well as the current market situation for construction projects of that kind. Also, the cost of capital is affected by the interest and depreciation time applied which e.g. depend on the perceived risk seen by the investor. Therefore, the cost of capital could vary over time but also between different projects. Given these uncertainties, the production cost is calculated for scenarios where the cost of capital and feedstock is changed by ± 50%, see Figure 4 and 5. Although these changes have a high impact on the overall production cost, they do not change the conclusion that biogas from energy crops has a lower production cost than biogas from residues and manure. When biogas is produced from straw and manure, pre-treatment and additives also represents a substantial and relatively uncertain part of the total production cost. Reducing these costs would increase the competitiveness of these feedstocks. However, the cost for pre-treatment and additives must be reduced with more than 75% before biogas from straw have a lower production cost than biogas from wheat grain.

    Figure 4. Biogas production cost when capital cost is changed by ± 50%.
    Figure 5. Biogas production cost when feedstock cost is changed by ± 50%.

    In Table 8, the break-even cost of capital and feedstock is also presented for each biogas system. Since the distribution cost applied in this study is based on a relatively well established market, where each filling station distribute 10 GWh annually, calculations are also made for a scenario where the volume is 5 GWh, representing a less developed market. For biogas produced from wheat grain and Triticale, the capital cost could be increased with 41–106% before the production cost equals the current market price. If the feedstock cost is changed instead, the cost could be increased with 29–57%. For biogas from straw, the cost of capital and feedstock must be reduced by 22–41% and 90% respectively to match the current market price.

    Table 8. Calculated break-even cost for capital and feedstock with an annual distributed volume of 10 and 5 GWh per filling station.
    Production cost (€/GJ)1 Capital cost (k€) Feedstock cost (€/t DM)
    Calculated1 Break-even Base case Break-even2
    Wheat grain 24.3–26.8 1091–1520 2247 157 245–212
    Beet 28.4–30.9 1102–1531 1557 180 215–182
    Maize 28.7–31.2 1104–1533 1506 160 186–158
    Triticale 25.0–27.5 1042–1471 2074 146 219–189
    Grass 31.2–33.7 1402–1830 1373 149 147–125
    Straw 32.8–35.3 1371–1800 1063 91 76–56
    Beet tops 30.7–33.2 1549–1978 1597 129 131–108
    Manure 31.9–34.4 2315–2744 2159 36 30–15
    1 Reducing the distributed volume from 10 GWh to 5 GWh at each filling station increase the annual capital cost with 429 k€ or 2.5 €/GJ.
    2 The interval represents the feedstock cost that results in a total cost for production and distribution of 31 €/GJ depending on the distributed volume at each filling station.
     | Show Table
    DownLoad: CSV

    3.5. Impact of policy instruments

    As mentioned earlier, different policy instruments affect biogas production and utilization. In the following, the impact of current investment and production subsidies as well as the exemption from CO2-and energy taxes applied in Sweden is presented.

    The calculated cost presented in Table 7 and the comparison between crops and residues presented in Figure 2 and 3 does not include any subsidies. Including them would reduce the cost and increase the competitiveness for biogas versus other fuels. Since both the investment subsidy as well as the production subsidy are unavailable for biogas plants using food crops they might also affect the comparison between different biogas systems.

    The impact of an investment subsidy on the cost for production and distribution of biogas is presented in Figure 6. The investment subsidy is set to 30% of the total investment for all biogas systems analyzed but the ones using wheat grain, Triticale, maize and beet. For these systems, only filling stations are subsidized.

    Figure 6. Production cost for biogas from crops and residues with and without subsidies and compared to current market prices of vehicle fuels.

    These assumptions reduce the difference between crops and residues, although straw and grass still give the highest cost and wheat grain the lowest. Since the cost of capital represents a high share of the total cost for biogas from manure, this is the feedstock where an investment subsidy has the highest impact. However, the production subsidy for biogas from manure could have an even higher impact. Including a production subsidy of 6 €/GJ (50% of the maximum subsidy) instead of an investment subsidy gives an overall cost of 26 €/GJ. In that case, biogas from manure would be the biofuel with the lowest cost excluding biogas from wheat grain and Triticale.

    In general, biogas from straw, grass and beet tops is not economically competitive to biogas from crops such as wheat grain and Triticale, even if investment subsidies are included. Compared to current market prices of fossil fuels, all biogas systems also seems to be competitive, some with very small margins, if investment subsidies are included, see Figure 7.

    Figure 7. Calculated production cost for biogas including investment subsidies and full CO2-and energy tax for biogas from food crops compared to current market prices of fossil vehicle fuels.

    Implementing these investment subsidies at the same time as biogas from food crops is burdened with the same energy-and CO2-taxes as natural gas used as a vehicle fuel, biogas from beet tops and manure would have the lowest cost. With the feedstock cost applied in this study, biogas from wheat grain would, however, still be competitive as compared to fossil fuels. Since wheat grain is easier and cheaper to transport than the other feedstock analyzed here, these result could also be more generally applicable.

    It should, however, be pointed out that the conditions described in this paper are not representative of Sweden in general. The conditions for producing crops and collect crop residues are likely more favorable than in many other parts of the country. Also, the availability of manure is high in the region analyzed here. The biogas volumes sold on each filling station is also higher than current average representing a more developed biogas market.

    Since the difference between the calculated cost and current market price is relatively small, especially for residues, small changes in the assumptions made here might also change the conclusion that biogas can compete with fossil fuels from an economic point of view. Due to this uncertainty, long term and clearly defined policy instruments will be necessary for actors to invest in the kind of biogas systems analyzed here.


    4. Conclusions

    When producing biogas from crops or residues, the choice of feedstock has a high impact on the overall cost for production and distribution.

    Biogas from Triticale and wheat grain can be competitive to fossil fuels even without subsidies. However, EU directives set out goals to favor the production of biofuels from waste and residues to that from dedicated energy crops. Without investment and production subsidies, biogas from grass, crop residues and manure is more expensive than biogas from several food crops. Even with current subsidies biogas from wheat grain and Triticale are economically more competitive than biogas from grass, wheat straw and beet tops. A shift from food crops to residues as biofuel feedstock, as desired according to EU directives, would thus require additional policy instruments favoring biogas from residues. Current production subsidy for biogas from manure could, however, make manure the most interesting feedstock from an economic point of view.

    An alternative or complementary approach could be to keep the current tax exemption for energy and CO2-taxes for biogas based on desired feedstock and not for biogas from crops. Even with this action, biogas from food crops could be more or less competitive to residues. This policy alone would also result in very small or non-existing economical margins for biogas production from all analyzed feedstocks.

    The reason to produce biofuels is to reduce different kinds of environmental impact. In this case, the main objective is to replace fossil vehicle fuels with biofuels that give large GHG emission reduction. When different policy instruments are discussed, it is important to keep this in mind. Given the policy instruments implemented today, biogas from both food crops and crop residues could, in some cases, compete with fossil fuels from an economic point of view. These instruments are, however, only implemented for a few more years and investments in biogas systems require a longer time perspective. Thus, strong and long-term policy signals are required if a future increase in production of biogas as an advanced biofuel is desired.


    Acknowledgements

    The funding from the Swedish Research Council Formas for the interdisciplinary research project Crops 4 Biogas (Resource efficient production of renewable energy carriers from energy crops) is gratefully acknowledged. We thank our co-workers in that project for interesting and fruitful discussions.


    Conflict of Interest

    All authors declared that, there is no conflict of interest in this paper.




    [1] E. O. Wilson, Sociobiology: The New Synthesis, Belknap Press, 1975.
    [2] R. A. Hinde, Ethology: Its Nature and Relation with Other Sciences, Oxford University Press, 1982.
    [3] H. Whitehead, Analyzing Animal Societies, University of Chicago Press, 2008.
    [4] A. F. Fraser, D. M. Broom, Farm Animal Behaviour and Welfare, CAB International, 1997.
    [5] E. O. Price, Animal Domestication and Behaviour, CABI Publishing, 2002.
    [6] R. A. Hinde, Primate Social Relationships, Blackwell Scientific Publications, 1983.
    [7] T. T. Strusaker, Correlates of ecology and social organization among African cercopithecines, Folia Primatol., 11 (1969), 80-118.
    [8] S. R. Sundaresan, I. R. Fishhoff, J. Dushoff, D. I. Rubenstein, Network metrics reveal differences in social organization between two fission-fusion species, Grevy's zebra and onager, Oecologia, 151 (2007), 140-149.
    [9] D. P. Croft, R. James, J. Krause, Exploring Animal Social Networks, Princeton University Press, 2008.
    [10] T. Wey, D. T. Blumstein, W. Shen, F. Jordan, Social network analysis of animal behaviour: A promising tool for the study of sociality, Anim. Behav., 75 (2008), 333-344.
    [11] C. Kasper, B. Voelkl, A social network analysis of primate groups, Primates, 50 (2009), 343-356.
    [12] J. Moreno, Who Shall Survive?, Beacon, 1934.
    [13] J. Scott, P. J. Carrington, The SAGE Handbook of Social Network Analysis, SAGE Publications, 2011.
    [14] G. A. Lundberg, M. Steel, Social attraction-patterns in a village, Sociometry, 1 (1938), 375-419.
    [15] J. A. Barnes, Class and committee in a Norwegian island parish, Hum. Relat., 7 (1954), 39-58.
    [16] J. H. Levine, The sphere of influence, Am. Sociol. Rev., 37 (1972), 14-27.
    [17] M. Granowetter, The strength of weak ties, Am. J. Sociol., 78 (1973), 1360-1380.
    [18] B. Ryan, N. C. Gross, The diffusion of hybrid seed corn in two Iowa communities, Rural Sociol., 8 (1943), 15-24.
    [19] E. Katz, H. Levine, M. L. Hamilton, Traditions of research on the diffusion of innovation, Am. Sociol. Rev., 28 (1963), 237-253.
    [20] S. Milgram, The small-world problem, Psychol. Today, 2 (1964), 60-67.
    [21] E. M. Rogers, Diffusion of Innovations, Free Press, 2003.
    [22] M. Bond, Social influences on corporate political donations in Britain, Brit. J. Sociol., 55 (2004), 55-77.
    [23] D. Knoke, Policy Networks, in The SAGE Handbook of Social Network Analysis (eds. J. Scott, P. J. Carrington), SAGE Publications, 2011, 210-222.
    [24] M. Diani, Social Movements and collective actions, in The SAGE Handbook of Social Network Analysis (eds. J. Scott, P. J. Carrington), SAGE Publications, 2011, 223-235.
    [25] M. O. Jackson, Social and Economic Networks, Princeton University Press, 2008.
    [26] A. L. Barabâsi, J. Hawoong, N. Zoltan, R. Erzsebet, A. Schubert, T. Vicsek, Evolution of the social network of scientific collaborations, Phys. A, 311 (2002), 590-614.
    [27] M. E. Newman, Coauthorship networks and patterns of scientific collaboration, Proc. Natl. Acad. Sci. USA, 101 (2004), 5200-5205.
    [28] D. J. Watts, S. H. Strogatz, Collective dynamics of "small-world" networks, Nature, 393 (1998), 440-442.
    [29] M. E. Newman, Networks: An Introduction, Oxford University Press, 2010.
    [30] H. Kummer, Soziales Verhalten einer Mantelpavianen-Gruppe, Schweizerische Zeitschr. Psychol., 33 (1957), 1-91.
    [31] D. S. Sade, Some aspects of parent-off spring and sibling relations in a group of rhesus monkeys, with a discussion of grooming, Am. J. Phys. Anthropol., 23 (1965), 1-18.
    [32] R. A. Hinde, Interactions, relationships and social structure, Man, 11 (1976), 1-17.
    [33] R. W. Byrne, A. Whiten, S. P. Henzi, Social relationships of mountain baboons: Leadership and affiliation in a non-female-bonded monkey, Am. J. Primatol., 207 (1989), 191-207.
    [34] B. D. Chepko-sade, K. P. Reitz, D. S. Sade, Sociometrics of Macaca mulatta IV: Network analysis of social structure of a pre-fission group, Soc. Netw., 11 (1989), 293-314.
    [35] C. A. Chapman, Association patterns of spider monkeys: The influence of ecology and sex on social organization, Behav. Ecol. Sociobiol., 26 (1990), 409-414.
    [36] C. P. Yeager, Proboscis monkey (Nasalis larvatus) social organization: Group structure, Am. J. Primatol., 106 (1990), 95-106.
    [37] D. S. Sade, Sociometrics of Macaca mulatta I. Linkage, cliques in grooming matrices, Fol. Primatol., 18 (1972), 196-223.
    [38] D. S. Sade, Sociometrics of Macaca mulatta Ⅲ. n-path centrality in grooming networks, Soc. Netw., 11 (1989), 273-292.
    [39] D. S. Sade, M. Altmann, J. Loy, G. Hausfater, J. A. Breuggeman, Sociometrics of Macaca mulatta: Ⅱ. Decoupling centrality and dominance in rhesus monkey social networks, Am. J. Phys. Anthropol., 77 (1988), 409-425.
    [40] S. Wasserman, K. Faust, Social Network Analysis: Methods and Applications, Cambridge University Press, 1994.
    [41] M. Barthelemy, B. Gondran, E. Guichard, Spatial structure of the internet traffic, Phys. A, 319 (2003), 633-642.
    [42] P. Bonacich, Factoring and weighting approaches to status scores and clique identification, J. Math. Sociol., 2 (1972), 113-120.
    [43] M. E. J. Newman, A measure of betweenness centrality based on random walks, Soc. Netw., 27 (2005), 39-54.
    [44] M. E. J. Newman, M. Girvan, Finding and evaluating community structure in networks, Phys. Rev. E, 69 (2004), 026113.
    [45] A. Clauset, Finding local community structure in networks, Phys. Rev. E, 72 (2005), 026132.
    [46] M. E. J. Newman, Modularity and community structure in networks, Proc. Natl. Acad. Sci. USA, 103 (2006), 8577-8582.
    [47] D. Knoke, S. Yang, Social network analysis, Sage Publications, 2019.
    [48] D. Lusseau, The emergent properties of a dolphin social network, Biol. Lett., 270 (2003), 186-188.
    [49] D. P. Croft, J. Krause, R. James, Social networks in the guppy (Poecilia reticulata), Biol. Lett., 271 (2004), 516-519.
    [50] B. Voelkl, Does group structure influence the social transmission of information?, Fol. Primatol., 75 (2004), 423.
    [51] D. P. Croft, R. James, A. J. W. Ward, M. S. Botham, D. Mawdsley, J. Krause, Assortative interactions and social networks in fish, Oecologia, 143 (2005), 211-219.
    [52] J. C. Flack, M. Girvan, F. B. M. de Waal, D. C. Krakauer, Policing stabilizes construction of social niches in primates, Nature, 439 (2006), 426-429.
    [53] C. Sueur, O. Petit, Organization of group members at departure Is driven by social structure in Macaca, Int. J. Primatol., 29 (2008), 1085-1089.
    [54] S. P. Henzi, D. Lusseau, T. Weingrill, Cyclicity in the structure of female baboon social networks, Behav. Ecol. Sociobiol., 63 (2009), 1015-1021.
    [55] J. Lehmann, C. Boesch, Sociality of the dispersing sex: The nature of social bonds in West African female chimpanzees, Pan troglodytes, Anim. Behav., 77 (2009), 377-387. doi: 10.1016/j.anbehav.2008.09.038
    [56] J. Lehmann, R. I. M. Dunbar, Network cohesion, group size and neocortex size in female-bonded Old World primates, Proc. R. Soc. B, 276 (2009), 4417-4422.
    [57] G. Ramos-Fernández, D. Boyer, F. Aureli, L. G. Vick, Association networks in spider monkeys (Ateles geoffroyi), Behav. Ecol. Sociobiol., 63 (2009), 999-1013.
    [58] N. J. B. Boogert, S. M. Reader, W. Hoppitt, K. N. Laland, The origin and spread of innovations in starlings, Anim. Behav., 75 (2008), 1509-1518.
    [59] B. Voelkl, R. Noë, The influence of social structure on the propagation of social information in artificial primate groups: A graph-based simulation approach, J. Theoret. Biol., 252 (2008), 77-86.
    [60] M. Franz, C. L. Nunn, Network-based diffusion analysis: A new method for detecting social learning, Proc. R. Soc. B, 276 (2009), 1829-1836.
    [61] C. Vital, P. Martins, Using graph theory metrics to infer information flow through animal social groups: A computer simulation analysis, Ethology, 115 (2009), 347-355.
    [62] W. Hoppitt, A. Kandler, J. R. Kendal, K. N. Laland, The effect of task structure on diffusion dynamics: Implications for diffusion curve and network-based analyses, Learn. Behav., 38 (2010), 243-251.
    [63] L. M. Aplin, D. R. Farine, J. Morand-Ferron, A. Cockburn, A. Thornton, B. C. Sheldon, Experimentally induced innovations lead to persistent culture via conformity in wild birds, Nature, 7540 (2015), 538.
    [64] B. Voelkl, C. Kasper, Social structure of primate interaction networks facilitates the emergence of cooperation, Biol. Lett., 5 (2009), 462-464.
    [65] B. Voelkl, The "Hawk-Dove" game and the spread of the evolutionary process in small heterogeneous populations, Games, 1 (2010), 103-116.
    [66] B. Voelkl, The evolution of generalized reciprocity in social interaction networks, Theoret. Popul. Biol., 104 (2015), 17-25.
    [67] B. Mccowan, K. Anderson, A. Heagarty, A. Cameron, Utility of social network analysis for primate behavioral management and well-being, Appl. Anim. Behav. Sci., 109 (2008), 396-405.
    [68] B. A. Beisner, M. E. Jackson, A. Cameron, B. Mccowan, Effects of natal male alliances on aggression and power dynamics in rhesus macaques, Am. J. Primatol., 801 (2011), 790-801.
    [69] V. Dufour, C. Sueur, A. Whiten, The impact of moving to a novel environment on social networks, activity and wellbeing in two new world primates, Am. J. Primatol., 811 (2011), 802-811.
    [70] M. C. Crofoot, D. I. Rubenstein, A. S. Maiya, T. Y. Berger-wolf, Aggression, grooming and group-level cooperation in white-faced capuchins (Cebus capucinus): Insights from social networks, Am. J. Primatol., 833 (2011), 821-833.
    [71] B. Tiddi, F. Aureli, G. Schino, B. Voelkl, Social relationships between adult females and the alpha male in wild tufted capuchin monkeys, Am. J. Primatol., 73 (2011), 812-820.
    [72] A. J. J. MacIntosh, A. Jacobs, C. Garcia, K. Shimizu, K. Mouri, M. A. Huffman, et al., Monkeys in the middle: Parasite transmission through the social network of a wild primate, PLoS One, 7 (2012), e51144.
    [73] L. J. N. Brent, S. Semple, Social capital and physiological stress levels in free-ranging adult female rhesus macaques, Behaviour, 102 (2011), 76-83.
    [74] P. C. Lopes, P. Block, B. König, Infection-induced behavioural changes reduce connectivity and the potential for disease spread in wild mice contact networks, Sci. Rep., 6 (2016), 31790.
    [75] M. Dow, F. B. M. de Waal, Assignment methods for the analysis of network subgroup interactions, Soc. Netw., 11 (1989), 237-255.
    [76] I. Matsuda, P. Zhang, L. Swedell, U. Mori, A. Tuuga, A., H. Bernard, et al., Comparisons of intraunit relationships in nonhuman primates living in multilevel social systems, Int. J. Primatol., 33 (2012), 1038-1053.
    [77] R. Milo, S. Shen-Orr, S. Itzkovitz, N. Kashtan, D. Chklovskii, U. Alon, Network motifs: Simple building blocks of complex networks, Science, 298 (2002), 824-827.
    [78] N. Snyder-Mackler, J. C. Beehner, T. J. Bergman, Defining higher levels in the multilevel societies of geladas (Theropithecus gelada), Int. J. Primatol., 33 (2012), 1054-1068.
    [79] C. Sueur, O. Petit, A. de Marco, A. T. Jacobs, K. Watanabe, B. Thierry, A comparative network analysis of social style in macaques, Anim. Behav., 82 (2011), 845-852.
    [80] B. Voelkl, R. Noë, Simulation of information propagation in real-life primate networks: Longevity, fecundity, fidelity, Behav. Ecol. Sociobiol., 64 (2010), 1449-1459.
    [81] D. I. Rubenstein, Networks of terrestrial ungulates: linking form and function, in Animal Social Networks (eds. J. Krause, R. James, D. W. Franks, D. P. Croft), Oxford University Press, 2015, 184-196.
    [82] E. A. Foster, D. W. Franks, L. J. Morrell, K. C. Balcomb, K. M. Parsons, A. van Ginneken, et al., Social network correlates of food availability in an endangered population of killer whales, Orcinus orca. Anim. Behav., 83 (2012), 731-736.
    [83] R. Albert, A. L. Barabasi, Statistical mechanics of complex networks, Rev. Mod. Phys., 74 (2002), 47-97.
    [84] M. E. J Newman, The structure and function of complex networks, SIAM Rev., 45 (2003), 167-256.
    [85] S. Macdonald, B. Voelkl, Primate social networks, in Animal Social Networks (eds. J. Krause, R. James, D. W. Franks, D. P. Croft), Oxford University Press, 2015,123-136.
    [86] R. R. Kao, L. Danon, D. M. Green, I. Z. Kiss, Demographic structure and pathogen dynamics on the network of livestock movements in Great Britain, Proc. R. Soc. B, 273 (2007), 1999-2007.
    [87] R. R. Kao, D. M. Green, J. Johnson, I. Z. Kiss, Disease dynamics over very different time-scales: Foot-and-mouth disease and scrapie on the network of livestock movements in the UK, J. R. Soc. Interface, 4 (2007), 907-916.
    [88] L. Danon, A. P. Ford, T. House, C. P. Jewell, M. J. Keeling, G. O. Roberts, et al., Networks and the epidemiology of infectious disease, Interdiscipl. Persp. Infect. Dis., 2011 (2011), 284909.
    [89] R. M. Anderson, R. M. May, Population biology of infectious diseases: Part 1, Nature, 280 (1979), 361-367.
    [90] M. J. E. Newman, The spread of epidemic disease on networks, Phys. Rev. E, 66 (2003), 016128.
    [91] R. M. May, Network structure and the biology of populations', Trends Ecol. Evol., 21 (2006), 394-399.
    [92] L. Hufnagel, D. Brockmann, T. Geisel, Forecast and control of epidemics in a globalized world, Proc. Natl. Acad. Sci. USA, 101 (2004), 7794-7799.
    [93] L. A. Meyers, Contact network epidemiology: Bond percolation applied to infectious disease prediction and control, Bull. Am. Math. Soc., 44 (2007), 63-86.
    [94] R. M. May, R. M Anderson, Transmission dynamics of HIV infection, Nature, 326 (1987), 137-142.
    [95] A. S. Klovdahl, J. J. Potterat, D. E. Woodhouse, J. B. Muth, S. Q. Muth, W. W. Darrow, Social networks and infectious disease: The Colorado Springs study, Soc. Sci. Med., 38 (1994), 79-88.
    [96] F. Liljeros, C. R. Edling, L. A. Nunes Amaral, E. Stanley, Y. Åberg, The web of human sexual contacts, Nature, 411 (2001), 907-908.
    [97] M. A. Nowak, Evolutionary Dynamics, Harvard University Press, 2006.
    [98] M. W. Schein, M. H. Fohrman, Social dominance relationships in a herd of dairy cattle, Brit. J. Anim. Behav., 3 (1955), 45-55.
    [99] M. Bigras-Poulin, R. A. Thompson, M. Chriel, S. Mortensen, M. Greiner, Network analysis of Danish cattle industry trade patterns as an evaluation of risk potential for disease spread, Prevent. Vet. Med., 76 (2006), 11-39.
    [100] L. Fiebig, T. Smieszek, J. Saurina, J. Hattendorf, J. Zinsstag, Contacts between poultry farms, their spatial dimension and their relevance for avian influenza preparedness, Geospat. Health, 4 (2009), 79-95.
    [101] B. Martinez-Lopez, A. M. Perez, J. M. Sanchez-Vizcaino, Social network analysis. Review of general concepts and use in preventive veterinary medicine, Transb. Emerg. Dis., 56 (2009), 109-120.
    [102] V. V. Volkova, R. Howey, N. J. Savill, M. E. J. Woolhouse, Sheep movement networks and the transmission of infectious diseases, PLoS One, 5 (2010), e11185.
    [103] R. P. Smith, A. J. C. Cook, R. M. Christley, Descriptive and social network analysis of pig transport data recorded by quality assured pig farms in the UK, Prevent. Vet. Med., 108 (2013), 167-177.
    [104] J. Ribeiro-Lima, E. A. Enns, B. Thompson, M. E. Craft, S. J. Wells, From network analysis to risk analysis-An approach to risk-based surveillance for bovine tuberculosis in Minnesota, US, Prevent. Vet. Med., 118 (2015), 328-340.
    [105] H. H. Lentz, A. Koher, P. Hövel, J. Gethmann, C. Sauter-Louis, T. Selhorst, et al., Disease spread through animal movements: a static and temporal network analysis of pig trade in Germany, PLoS One, 11 (2016), 0155195.
    [106] P. Bajardi, A. Barrat, L. Savini, V. Colizza, Optimizing surveillance for livestock disease spreading through animal movements, J. R. Soc. Interface, 9 (2012), 2814-2825.
    [107] M. M. Mweu, G. Fournié, T. Halasa, N. Toft, S. S. Nielsen, Temporal characterisation of the network of Danish cattle movements and its implication for disease control: 2000-2009, Prevent. Vet. Med., 110 (2013), 379-387.
    [108] S. Nickbakhsh, L. Matthews, J. E. Dent, G. T. Innocent, M. E. Arnold, S. W. Reid, et al., Implications of within-farm transmission for network dynamics: Consequences for the spread of avian influenza, Epidemics, 5 (2013), 67-76.
    [109] B. Vidondo, B. Voelkl, Dynamic network measures reveal the impact of cattle markets and alpine summering on the risk of epidemic outbreaks in the Swiss cattle population, BMC Vet. Res. 14 (2018), 88.
    [110] J. Krause, D. Lusseau, R. James, Animal social networks: An introduction, Behav. Ecol. Sociobiol., 63 (2009), 967-973.
    [111] M. J. Silk, D. P. Croft, R. J. Delahay, D. J. Hodgson, M. Boots, N. Weber, et al., Using social network measures in wildlife disease ecology, epidemiology, and management, BioScience, 67 (2017), 245-257.
    [112] M. E. Craft, Infectious disease transmission and contact networks in wildlife and livestock, Phil. Trans. R. Soc. B, 370 (2015), 20140107.
    [113] R. H. Griffin, C. L. Nunn, Community structure and the spread of infectious disease in primate social networks, Evol. Ecol., 26 (2012), 779-800.
    [114] C. L. Nunn, F. Jordan, C. M. McCabe, J. L. Verdolin, J. H. Fewell, Infectious disease and group size: More than just a numbers game, Phil. Trans. R. Soc. B, 370 (2015), 20140111.
    [115] S. S. Godfrey, C. M. Bull, R. James, K. Murray, Network structure and parasite transmission in a group living lizard, the gidgee skink, Egernia stokesii, Behav. Ecol. Sociobiol., 63 (2009), 1045-1056.
    [116] K. L. VanderWaal, E. R. Atwill, S. Hooper, K. Buckle, B. McCowan, Network structure and prevalence of Cryptosporidium in Belding's ground squirrels, Behav. Ecol. Sociobiol., 67 (2013), 1951-1959.
    [117] T. Porphyre, M. Stevenson, R. Jackson, J. McKenzie, Original article Influence of contact heterogeneity on TB reproduction ratio R0 in a free-living brushtail possum Trichosurus vulpecula population, Vet. Res., 39 (2008), 31.
    [118] J. Rushmore, D. Caillaud, R. J. Hall, R. M. Stumpf, L. A. Meyers, S. Altizer, Network-based vaccination improves prospects for disease control in wild chimpanzees, J. R. Soc. Interface, 11 (2014), 20140349.
    [119] J. A. Drewe, K. T. D. Eames, J. R. Madden, G. P. Pearce, Integrating contact network structure into tuberculosis epidemiology in meerkats in South Africa: Implications for control, Prevent. Vet. Med., 101 (2011), 113-120.
    [120] M. D. J. Blyton, S. C. Banks, R. Peakall, D. B. Lindenmayer, D. M. Gordon, Not all types of host contacts are equal when it comes to E. coli transmission, Ecol. Lett., 17 (2014), 970-978.
    [121] C. R. Webb, Farm animal networks: Unraveling the contact structure of the British sheep population, Prevent. Vet. Med., 68 (2005), 3-17.
    [122] F. Natale, A. Giovannini, L. Savini, D. Palma, L. Possenti, G. Fiore, et al., Network analysis of Italian cattle trade patterns and evaluation of risks for potential disease spread, Prevent. Vet. Med., 92 (2009), 341-350.
    [123] C. Dubé, C. Ribble, D. Kelton, B. Mcnab, A Review of network analysis terminology and its application to foot-and-mouth disease modelling and policy development, Transbound. Emerg. Dis., 56 (2009), 73-85.
    [124] H. Chen, G. Smith, S. Zhang, K. Qin, J. Wang, S. Li, et al., H5N1 virus outbreak in migratory waterfowl, Nature, 436 (2005), 191-192.
    [125] B. J. Hoye, V. J. Munster, H. Nishiura, R. A. M. Fouchier, J. Madsen, M. Klaassen, Reconstructing an annual cycle of interaction: Natural infection and antibody dynamics to avian influenca along a migratory flyway, Oikos, 120 (2011), 748-755.
    [126] K. R. Finn, M. J. Silk, M. A. Porter, N. Pinter-Wollman, The use of multilayer network analysis in animal behaviour, Anim. Behav., 149 (2019), 7-22.
    [127] K. Robert, D. Garant, F. Pelletier, Keep in touch: Does spatial overlap correlate with contact rate frequency?, J. Wildl. Manag., 76 (2012), 1670-1675.
    [128] M. L. Gilbertson, L. A. White, M. E. Craft, Trade‐offs with telemetry‐derived contact networks for infectious disease studies in wildlife, Meth. Ecol. Evol., 2020.
    [129] S. E. Perkins, F. Cagnacci, A. Stradiotto, D. Arnoldi, P. J. Hudson, Comparison of social networks derived from ecological data: implications for inferring infectious disease dynamics, J. Anim. Ecol., 78 (2009), 1015-1022.
    [130] J. Krause, A. D. M. Wilson, D. P. Croft, New technology facilitates the study of social networks, Trends Ecol. Evol., 26 (2011), 5-6.
    [131] C. Rutz, Z. T. Burns, R. James, S. M. H. Ismar, J. Burt, B. Otis, et al., Automated mapping of social networks in wild birds, Curr. Biol., 22 (2012), R669-R671.
    [132] I. Psorakis, B. Voelkl, C. J. Garroway, R. Radersma, L. M. Aplin, R. A. Crates, et al., Inferring social structure from temporal data, Behav. Ecol. Sociobiol., 69 (2015), 857-866.
    [133] J. R. Ginsberg, T. P. Young, Measuring associations between individuals or groups in behavioural studies, Anim. Behav., 44 (1992), 377-379.
    [134] L. Beijder, D. Fletcher, S. Brager, A method for testing association patterns of social animals, Anim. Behav., 56 (1998), 719-725.
    [135] L. A. White, J. D. Forester, M. E. Craft, Using contact networks to explore mechanisms of parasite transmission in wildlife, Biol. Rev., 92 (2017), 389-409.
    [136] R. K. Hamede, J. Bashford, H. McCallum, M. Jones, Contact networks in a wild Tasmanian devil (Sarcophilus harrisii) population: Using social network analysis to reveal seasonal variability in social behaviour and its implications for transmission of devil facial tumour disease, Ecol. Lett., 12 (2009), 1147-1157.
    [137] T. C. Germann, K. Kadau, I. M. Longini, C. A. Macken, Mitigation strategies for pandemic influenza in the United States, Proc. Natl. Acad. Sci. USA, 103 (2006), 5935-5940.
    [138] S. E. Robinson, M. G. Everett, R. M. Christley, Recent network evolution increases the potential for large epidemics in the British cattle population, J. R. Soc. Interface, 4 (2007), 669-674.
    [139] J. C. Gibbens, C. E. Sharpe, J. W. Wilesmith, L. M. Mansley, E. Michalopoulou, J. B., et al., Descriptive epidemiology of the 2001 foot-and-mouth disease epidemic in Great Britain: The first five months, Vet. Rec., 149 (2001), 729-743.
    [140] I. Z. Kiss, D. M. Green, R. R. Kao, The network of sheep movements within Great Britain: Network properties and their implications for infectious disease spread, J. R. Soc. Interface, 3 (2006), 669-677.
    [141] D. M. Green, I. Z. Kiss, R. R. Kao, Modelling the initial spread of foot-and-mouth disease through animal movements, Proc. R. Soc. B, 273 (2006), 2729-2735.
    [142] M. C.Vernon, M. J. Keeling, Representing the UK's cattle herd as static and dynamic networks, Proc. R. Soc. B, 276 (2009), 469-476.
    [143] P. Sah, S. T. Leu, P. C. Cross, P. J. Hudson, S. Bansal, Unravelling the disease consequences and mechanisms of modular structure in animal social networks, Proc. Natl. Acad. Sci. USA, 114 (2017), 4165-4170.
    [144] R. Pastor-Satorras, A. Vespignani, Epidemic spreading in scale-free networks, Phys. Rev. Lett., 86 (2001), 3200-3203.
    [145] A. L. Lloyd, R. M. May, How viruses spread among computers and people, Science, 292 (2001), 1316-1317.
    [146] D. C. Bell, J. S. Atkinson, J. W. Carlson, Centrality measures for disease transmission networks, Soc. Netw., 21 (1999), 1-21.
    [147] M. J. Keeling, The effects of local spatial structure on epidemiological invasions, Proc. R. Soc. B, 266 (1999), 859-867.
    [148] K. T. D. Eames, M. J. Keeling, Modeling dynamic and network heterogeneities in the spread of sexually transmitted diseases, Proc. Natl. Acad. Sci. USA, 99 (2002), 13330-13335.
    [149] C. Buckee, L. Danon, S. Gupta, Host community structure and the maintenance of pathogen diversity, Proc. R. Soc. B, 274 (2007), 1715-1721.
    [150] M. Salathé, J. H. Jones, Dynamics and control of diseases in networks with community structure, PLoS Comp. Biol., 6 (2010), e1000736.
    [151] S. M. Firestone, M. P. Ward, R. M. Christley, N. K. Dhand, The importance of location in contact networks: Describing early epidemic spread using spatial social network analysis, Prevent. Vet. Med., 102 (2011), 185-195.
    [152] J. Frössling, A. Ohlson, C. Björkman, N. Hakansson, M. Nöremark, Application of network analysis parameters in risk-based surveillance-Examples based on cattle trade data and bovine infections in Sweden, Prevent. Vet. Med., 105 (2012), 202-208.
    [153] L. García Álvarez, C. R. Webb, M. A. Holmes, A novel field-based approach to validate the use of network models for disease spread between dairy herds, Epidemiol. Infect., 139 (2011), 1863-1874.
    [154] R. Biek, A. G. Rodrigo, D. Holley, A. Drummond, C. R. Anderson, H. A. Ross, et al., Epidemiology, genetic diversity, and evolution of endemic feline immunodeficiency virus in a population of wild cougars, J. Virol., 77 (2003), 9578-9589.
    [155] B. T. Grenfell, O. G. Pybus, J. R. Gog, J. L. N. Wood, J. M. Daly, J. A. Mumford, et al., Unifying the epidemiological and evolutionary dynamics of pathogens, Science, 303 (2004), 327-333.
    [156] R. Biek, A. Drummond, M. Poss, A virus reveals population structure and recent demographic history of its carnivore host, Science, 311 (2006), 538-542.
    [157] E. A. Archie, G. Luikart, V. O. Ezenwa, Infecting epidemiology with genetics: A new frontier in disease ecology, Trends Ecol. Evol., 24 (2008), 21-30.
    [158] C. M. Bull, S. S. Godfrey, D. M. Gordon, Social networks and the spread of Salmonella in a sleepy lizard population, Mol. Ecol., 21 (2012), 4386-4392.
    [159] K. L. VanderWaal, E. R. Atwill, L. A. Isbell, B. McCowan, Linking social and pathogen transmission networks using microbial genetics in giraffe (Giraffa camelopardalis), J. Anim. Ecol., 83 (2014), 406-414.
    [160] K. L. VanderWaal, E. R. Atwill, L. A. Isbell, B. McCowan, B. Quantifying microbe transmission networks for wild and domestic ungulates in Kenya, Biol. Conserv., 169 (2014), 136-146.
    [161] J. S. Lee, E. W. Ruell, E. E. Boydston, L. M. Lyren, R. S. Alonso, J. L. Troyer, et al., Gene flow and pathogen transmission among bobcats (Lynx rufus) in a fragmented urban landscape, Mol. Ecol., 21 (2012), 1617-1631. doi: 10.1111/j.1365-294X.2012.05493.x
    [162] B. Y. Reis, I. S. Kohane, K. D. Mandl, An epidemiological network model for disease outbreak detection, PLoS Med, 4 (2007), e210.
    [163] F. Schirdewahn, V. Colizza, H. H. Lentz, A. Koher, V. Belik, P. Hövel, Surveillance for outbreak detection in livestock-trade networks, in Temporal Network Epidemiology (eds. M. Naoki, P. Holme), Springer, 2017, 215-240.
    [164] P. Skums, A. Kirpich, P. I. Baykal, A. Zelikovsky, G. Chowell, Global transmission network of SARS-CoV-2: From outbreak to pandemic, medRxiv, 2020.
    [165] D. Lusseau, H. Whitehead, S. Gero, Incorporating uncertainty into the study of animal social networks, Anim. Behav., 75 (2008), 1809-1815.
    [166] R. James, D. P. Croft, J. Krause, Potential banana skins in animal social network analysis, Behav. Ecol. Sociobiol., 63 (2009), 989-997.
    [167] B. Voelkl, C. Kasper, C. Schwab, Network measures for dyadic interactions: Stability and reliability, Am. J. Primatol., 73 (2011), 731-740.
    [168] J. Krause, S. Krause, R. Arlinghaus, I. Psorakis, S. Roberts, C. Rutz, Reality mining of animal social systems, Trends Ecol. Evol., 28 (2013), 541-551.
    [169] M. Berdoy, J. P. Webster, D. W. Macdonald, Fatal attraction in rats infected with Toxoplasma gondii, Proc. R. Soc. B, 267 (2000), 1591-1594.
    [170] A. Vyas, S. Kim, N. Giacomini, J. C. Boothroyd, R. M. Sapolsky, Behavioral changes induced by Toxoplasma infection of rodents are highly specific to aversion of cat odors, Proc. Natl. Acad. Sci. USA, 104 (2007), 6442-6447.
    [171] D. P. Croft, M. Edenbrow, S. K. Darden, I. W. Ramnarine, C. van Oosterhout, J. Cable, Effect of gyrodactylid ectoparasites on host behaviour and social network structure in guppies Poecilia reticulata, Behav. Ecol. Sociobiol., 65 (2011), 2219-2227.
    [172] F. J. Theis, L. V. Ugelvig, C. Marr, S. Cremer, Opposing effects of allogrooming on disease transmission in ant societies, Phil Trans. R. Soc. B, 370 (2015), 20140108.
    [173] V. O. Ezenwa, E. A. Archie, M. E. Craft, D. M. Hawley, L. B. Martin, J. Moore, et al., Host behaviour-parasite feedback: an essential link between animal behaviour and disease ecology, Proc. R. Soc. B, 283 (2016), 20153078.
    [174] L. A. White, J. D. Forester, M. E. Craft, Dynamic, spatial models of parasite transmission in wildlife: Their structure, applications and remaining challenges, J. Anim. Ecol., 87 (2018), 559-580.
    [175] K. Büttner, J. Salau, J., Krieter, Quality assessment of static aggregation compared to the temporal approach based on a pig trade network in Northern Germany, Prevent. Vet. Med., 129 (2016), 1-8.
    [176] M. J. Silk, D. J., Hodgson, C. Rozins, D. P. Croft, R. J. Delahay, M. Boots, et al., Integrating social behaviour, demography and disease dynamics in network models: applications to disease management in declining wildlife populations, Phil. Trans. R. Soc. B, 374 (2019), 20180211.
    [177] M. E. Craft, E. Volz, C. Packer, L. A. Meyers, Disease transmission in territorial populations: The small-world network of Serengeti lions, J. R. Soc. Interface, 8 (2011), 776-786.
    [178] N. Weber, S. P. Carter, S. R. X. Dall, R. J. L. Delahay, J. L. McDonald, S. Bearhop, et al., Badger social networks correlate with tuberculosis infection, Curr. Biol., 23 (2013), R915-R916.
    [179] K. P. Huyvaert, R. E. Russell, K. A. Patyk, M. E. Craft, P. C. Cross, M. G. Garner, et al., Challenges and opportunities developing mathematical models of shared pathogens of domestic and wild animals, Vet. Sci., 5 (2018), 92.
    [180] S. Kraberger, N. M. Fountain-Jones, R. B. Gagne, J. Malmberg, N.G. Dannemiller, K. Logan, et al., Frequent cross-species transmissions of foamy virus between domestic and wild felids, Virus Evol., 6 (2020), vez058.
    [181] R. K. Plowright, C. R. Parrish, H. McCallum, P. J. Hudson, A. I. Ko, A. L. Graham, et al., Pathways to zoonotic spillover, Nat. Rev. Microbiol., 15 (2017), 502.
    [182] B. J. Coburn, B. G. Wagner, S. Blower, Modeling influenza epidemics and pandemics: insights into the future of swine flu (H1N1). BMC Med., 7 (2009), 30.
    [183] C. M. Scoglio, C. Bosca, M. H. Riad, F. D. Sahneh, S. C. Britch, L. W. Cohnstaedt, et al., Biologically informed individual-based network model for Rift Valley fever in the US and evaluation of mitigation strategies, PloS One, 11 (2016), e0162759.
    [184] S. K. Lau, P. C. Woo, K. S. Li, Y. Huang, H. W. Tsoi, B. H. Wong, et al., Severe acute respiratory syndrome coronavirus-like virus in Chinese horseshoe bats, Proc. Natl. Acad. Sci. USA, 102 (2005), 14040-14045.
    [185] C. M. Luo, N. Wang, X. L. Yang, H. Z. Liu, W. Zhang, B. Li, et al., Discovery of novel bat coronaviruses in south China that use the same receptor as Middle East respiratory syndrome coronavirus, J. Virol., 92 (2018), e00116-18.
    [186] N. Wang, S. Y. Li, X. L. Yang, H. M. Huang, Y. J. Zhang, H. Guo, et al., Serological evidence of bat SARS-related coronavirus infection in humans, China, Virol. Sin., 33 (2018),104-107.
    [187] L. E. Escobar, R. Moen, M. E. Craft, K. L. VanderWaal, Mapping parasite transmission risk from white-tailed deer to a declining moose population, Eur. J. Wildl. Res., 65 (2019), 60.
    [188] P. Sah, J. Mann, S. Bansal, Disease implications of animal social network structure: A synthesis across social systems, J. Anim. Ecol., 87 (2018), 546-558.
  • This article has been cited by:

    1. Madhurjya Das, Mainak Chakraborty, Senthil Chinnasamy, Balasubramanian Ramakrishnan, 2020, Chapter 13, 978-3-030-58826-7, 285, 10.1007/978-3-030-58827-4_13
    2. Tong Liu, Anna Schnürer, Johanna Björkmalm, Karin Willquist, Emma Kreuger, Diversity and Abundance of Microbial Communities in UASB Reactors during Methane Production from Hydrolyzed Wheat Straw and Lucerne, 2020, 8, 2076-2607, 1394, 10.3390/microorganisms8091394
    3. Lovisa Björnsson, Thomas Prade, Sustainable Cereal Straw Management: Use as Feedstock for Emerging Biobased Industries or Cropland Soil Incorporation?, 2021, 1877-2641, 10.1007/s12649-021-01419-9
    4. Thomas Prade, Sven-Erik Svensson, Torsten Hörndahl, Emma Kreuger, Impact of Harvest Date and Cutting Length of Grass Ley and Whole-Crop Cereals on Methane Yield and Economic Viability as Feedstock for Biogas Vehicle Fuel Production, 2019, 12, 1939-1234, 137, 10.1007/s12155-018-9952-x
    5. Mikael Lantz, Thomas Prade, Serina Ahlgren, Lovisa Björnsson, Biogas and Ethanol from Wheat Grain or Straw: Is There a Trade-Off between Climate Impact, Avoidance of iLUC and Production Cost?, 2018, 11, 1996-1073, 2633, 10.3390/en11102633
    6. Yanlei Feng, Lorenzo Rosa, Global biomethane and carbon dioxide removal potential through anaerobic digestion of waste biomass, 2024, 19, 1748-9326, 024024, 10.1088/1748-9326/ad1e81
    7. Emma Kreuger, Virginia Tosi, Maja Lindblad, Åsa Davidsson, Co-Digestion and Mono-Digestion of Sewage Sludge and Steam-Pretreated Winter Wheat Straw in Continuous Stirred-Tank Reactors—Nutrient Composition and Process Performance, 2024, 10, 2311-5637, 414, 10.3390/fermentation10080414
  • Reader Comments
  • © 2020 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(5951) PDF downloads(235) Cited by(0)

Article outline

Figures and Tables

Figures(1)  /  Tables(1)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog