Life Cycle Assessment modeling of milk production in Iran

Accepted Manuscript Life Cycle Assessment Modeling of Milk Production in Iran Hamzeh Soltanali, Bagher Emadi, Abbas Rohani, Mehdi Khojastehpour, Amin Nikkhah PII: DOI: Reference:

S2214-3173(15)00030-X http://dx.doi.org/10.1016/j.inpa.2015.06.003 INPA 27

To appear in:

Information Processing in Agriculture

Received Date: Revised Date: Accepted Date:

17 March 2015 12 June 2015 14 June 2015

Please cite this article as: H. Soltanali, B. Emadi, A. Rohani, M. Khojastehpour, A. Nikkhah, Life Cycle Assessment Modeling of Milk Production in Iran, Information Processing in Agriculture (2015), doi: http://dx.doi.org/10.1016/ j.inpa.2015.06.003

This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.

Life Cycle Assessment Modeling of Milk Production in Iran Hamzeh Soltanali a, Bagher Emadi a,*, Abbas Rohani a, Mehdi Khojastehpour a and Amin Nikkhah b a

Department of Biosystems Engineering, Ferdowsi University of Mashhad, Azadi Square, Mashhad, Iran * Corresponding author. Fax: +98 51 38805838, [email protected] b Young Researchers and Elite Club, Rasht Branch, Islamic Azad University, Rasht, Iran [email protected], http://orcid.org/0000-0002-3008-2401

Abstract Livestock units are known as one of the most influential sectors in the environment pollution. Therefore, the aim of this study was to investigate the environmental impacts of milk production in Guilan province of Iran through Life Cycle Assessment (LCA) methodology. The primary data were collected from 45 units of milk production through a field survey with the help of a structured questionnaire. The reliability was assessed using Cronbach's alpha coefficient and was estimated an acceptable value of 0.91. The consumption of resources and emissions were allocated to a Functional Unit (FU) of one ton of milk. Impacts of emissions in five impact categories of global warming, acidification, eutrophication, photochemical oxidation and depletion of resources were investigated. The results showed that the characterization index for these impact categories were 1831 kg CO2 eq, 7.97 kg SO2 eq, 3.42 kg PO4-3 eq, 0.21 kg C2H4 eq and 838.39 MJ, respectively. Final indices for these impact categories were calculated as 0.24, 0.28, 0.076, 0.017 and 0.046, respectively. Environmental index (EcoX) and resources depletion index (RDI) were obtained 0.61 and 0.04, respectively. In this study, the highest potential for environmental impacts of production revealed for acidification and followed by global warming impact category. Keywords: Acidification, Dairy farm, Environmental impact, Global warming

1. Introduction Besides securing food for the rising population, according to predictions will be more than 9 billion in 2075, the environment preservation is one of the most important challenges facing humanities [1]. Agriculture is one of the main sources of greenhouse gas (GHG) emissions such as CO2, CH4 and N2O [2,3]. It is necessary to evaluate environmental impacts of different sectors of agriculture particularly dairy cow breeding units that are one of the main resources of GHG emissions [4]. The FAO report "Livestock long shadow: environmental issues and options’’ claims that livestock units constitute about 18% of the entire GHG emissions of which 3% to 5.1% are from dairy cow breeding units [5,6]. Although there are different methods for assessing the environmental impacts associated with food systems, but the appropriate one is the Life Cycle Assessment (LCA), because this method is suitable for determination of inputs, outputs and environmental consequence in a production process [7,8]. A Literature review disclosed that many researchers have reported the valuable application of LCA model in environmental management of agricultural production [9-14]. Some studies have been conducted in terms of environmental impacts assessment during dairy farm activities, for example Van der Werf et al [15] evaluated the environmental impact of dairy cow breeding in terms of traditional and organic systems. They investigated the impact categories of global warming, acidification, eutrophication, land occupation and fossil resources depletion. They claimed that the impact of land occupation in the organic system was higher in comparison to a traditional system, while the other impact categories had not many differences in the two systems. Castanheira et al. [16] studied environmental impacts of dairy cow breeding units in Portugal. The biggest sources of N2O, NH3 and CH4 were reported to be diesel fuel, manure management and enteric fermentation.

McGeough et al. [17] performed a study on LCA in dairy units of Canada and reported that the amount of GHG emission per one liter of milk was 0.92 kg CO2 eq. O’Brien et al. [18] studied environmental impacts of milk production in Ireland. The results indicated that the effective environmental impacts were global warming, acidification, eutrophication, land occupation and fossil resources depletion for each functional unit of milk production. LCA model was also applied to evaluate environmental impacts resulted from global warming potential in milk production units in the United States. The total GHG per liter of milk was equal to 0.20 kg CO2 eq [19]. Guerci et al. [20] analyzed the environmental impacts of dairy farms activities in Denmark, Germany and Italy through LCA and found that the average annual production of milk per cow was between 6275 to 10964 liters. Effective environmental impact categories were acidification, eutrophication, depletion of fossil sources and global warming. The impacts of all of the effective environmental groups in dairy units of Denmark were less than Germany and Italy. In another study, carried out by Zhang et al. [21] they examined the environmental impacts for dairy cow breeding units in Canada using LCA methodology based on an integrated system. The integrated system in the study was for applying an aerobically digested for production of biogas energy and digested slurry. The results indicated that usage of an integrated system in livestock units can decrease fossil recourses and global warming up to 80% and eutrophication and respiratory effects up to 50%. Some other studies which were conducted on environmental impacts through LCA in dairy cow breeding units are also available [14,22,23,24,25,26,27]. The sustainable production of milk in Guilan province of Iran requires the consideration of environmental management in the production systems. However, to the best of knowledge of authors, no previous analytical work has been reported on the environmental impacts of milk

production in Iran. Therefore, the aim of present study was the environmental impacts assessment of milk production in dairy cow breeding industrial units in Guilan province, Iran based on LCA methodology.

2. Materials and methods 2.1. Site of study and sample selection The study was carried out in industrial units of dairy cow breeding of Guilan province in northern Iran during agricultural year of 2012-2013. Guilan province is located in the north of Iran on the south of Caspian Sea, within 36◦ 34‫ ׳‬and 38◦ 27‫ ׳‬north latitude and 48◦ 53‫ ׳‬and 50◦ 34‫׳‬ east longitude. Guilan has a population of approximately 2.5 million people [28]. From 180 industrial units of cow growing in Guilan at the time of the study, 129 of which were dairy cow breeding units [29]. The average numbers of cow in the studied area were ranged 20 to 200 cows that average was equal to 54.5 head. The most using of machinery were for the process of milking equipment, processing of animal feed, and followed by allocated milk cooling. In addition, many activities such as animal feeding operation were done by human labor. Based on the Cochran formula, 45 active units were selected for this study [30].

(1)

(2) where, n = sample size, N = number of holdings in the target population, t = the reliability coefficient (1.96), s = the variance, and d = precision. Data were obtained from 45 units of milk production using a face to face questionnaire method during 2012-2013. The instrument used in this research was a questionnaire whose validity was

confirmed by university faculty members and agricultural experts. The estimated reliability, using Cronbach’s Alpha, was 0.91, which was an acceptable reliability. Each producer was asked to detail activity as inputs to milk production recorded as diesel fuel (lit), electricity (kWh), natural gas (m3), animal feed (kg) and manure (kg), and as the output yield (kg).

2.2. Life Cycle Assessment This method comprises four sections as: goal and scope definition, life cycle inventory (LCI), life cycle impact assessment (LCIA) and finally interpretation [8,31,32].

2.2.1. Goal and scope definition The purpose of this study was to assess the environmental performance of milk production in dairy farming industries in Guilan province, Iran. The Functional Unit (FU) is connected to the inputs and outputs and provides a condition for comparison, which is usually equivalent to one ton of milk [16,33]. In the current study, the functional unit was chosen as one ton produced milk.

2.2.2. System Boundary In first, inputs including fuel, electricity, animal feed and other activities such as animal manure management (transport and accumulation) and enteric fermentation that potentially have adverse effects on the environment were identified. Then, consumption and emissions for one ton milk production in dairy units were calculated. The system outputs contained all of the pollutants emitted to the environment and were calculated out of using these inputs based on a

functional unit. Fig. 1 shows all of the inputs and outputs for milk production in Guilan dairy system.

Fig. 1. Inputs and outputs of milk production in Guilan province, Iran (system boundary)

2.2.2.1.

The inputs of the system

Fuels, electricity, animal feed, enteric fermentation and manure for production of one ton of milk are illustrated in the Table 1.

Table 1. The inputs and output to produce one ton of milk in Guilan province, Iran

2.2.2.2.

The outputs of the system

In this section, the level of the emitted pollutants from potentially pollutant inputs was assessed. Table 2 displays the emission coefficient from inputs consumption in milk production units.

Table 2. Emission coefficients of contaminants resulting from the use of inputs and activities

2.2.3. Impact Assessment The purpose of this section was to better analyze and interpret outputs and inputs of the system for production of one ton of milk production which was conducted in three phases of characterization, normalization and weighting [37]. In characterization phase, after dividing the impacts into different effective factors, the amount of pollutants resulted from milk production

will be calculated regarding its efficiency in different impact categories. This study investigated five impact categories including: depletion of fossil sources, global warming, eutrophication, acidification and photochemical oxidation. The index for characterization of each impact category is calculated through Eq. (2) [37-39]: (2) where ICIi is indicator value per functional unit for impact category i; Ej or Rj, emission of j mixture or consumption of j resource on each functional unit; CFi is the characterization factor for j mixture in impact category of i. The characterization factor in each impact factor shows the mixture potential for creating the impact. The characterization index for each of the impact categories and the efficiency of each mixture are shown in Table 3.

Table 3. Characterization factor of environmental impacts

2.2.3.1.

Normalization and Weighting

In order to have comprehensive understanding of the calculated amounts, characterization index of impact categories get into non-dimensions. Normalization factors are shown in Table 4. To better understand the level of damage of each impact categories on the environment, weighting stage was used and higher resulted value in this stage shows more potential risk for damaging the environment. Weighting factors for five impact categories such as global warming, acidification, eutrophication, depletion of fossil sources and photochemical oxidation are mentioned in Table 4.

Table 4. Normalization and weighting factors

Groups of global warming, acidification, eutrophication and photochemical oxidation affect environmental impacts (EcoX) and depletion of fossil sources affect the depletion of resources (RDI). The difference between the groups of EcoX and RDI index at resource depletion is that the latter affects future generations of human [37-39].

3. Results and discussion Table 5 illustrates the amount of emitted pollutants from the consumption of inputs and activities for production of one ton of milk in Guilan province, Iran. Inputs of animal feed, electricity and diesel fuel were the key factors in environmental hazards in most impact categories. In a similar study in Portugal, Castanheira et al. [16] showed that animal feed in that concentrate was the main sources of pollutants emission in milk production. Given study in the USA, animal feed and management of animal manure were introduced as the main factor of weather pollution [43]. As shown in Table 5, enteric fermentation and management of manure were the next pollutants of dairy units in Guilan province.

Table 5. The emissions of contaminants resulting from the consumption of inputs and activities for each ton of milk production

Fig. 2 shows the total emissions of pollutants in industrial dairies in Guilan province. The total emissions of pollutants CH4, N2O, NOx, CO2, NH3, SO2, CO, NO3 and PO4-3 for one ton of milk production were obtained 28.06, 1.60, 4.82, 744.71, 2.80, 0.90, 0.77, 16.80 and 0.14 kg eq, respectively. Among the pollutants CO2, CH4 and NO3 had the most roles in environment pollution, respectively. The prime factor for CO2 emission in dairy farms were animal feed,

electricity and diesel fuel which were 676.3, 38.93 and 24.98 kg for one ton of milk, respectively. Among animal feed, concentrate was identified as the main factor for the emission of CO2. Enteric fermentation was the main factor of the emission of CH4 that was equal to 26.57 kg eq for one ton of milk. The reason for high level of methane emission was related to the animal feed consumed. Grainger et al. [44] showed that using appropriate strategies such as: the increase of oil substances in livestock diet in the management of animal feed can minimize the emission of methane gas due to enteric fermentation, 1% increase in oil substances in livestock would lead to 6% decrease in methane emission. In a similar study carried out in dairy farms of Portugal N2O, NH3 and CH4 were introduced as the biggest pollutants. Enteric fermentation was the main source of CH4 emission [16] in another study carried out in Canada 56% of CH4 emission was reported as enteric fermentation [45] Conclusion of study in Guilan farms indicated that concentrate was the main reason for emission of NO3, which was equal to 16.8 kg for one ton of milk production.

Fig. 2. The total emissions of pollutants (kg) in industrial dairies in Guilan province, Iran

Table 6 displays the indices of characterization, normalization and the final index of environmental impact categories of milk production in Guilan province, Iran. The characterization indices of global warming, acidification, eutrophication, photochemical oxidation and depletion of resources for each ton of milk production were computed 1831 kg CO2 eq, 7.97 kg SO2 eq, 3.42 kg PO4-3 eq, 0.21 kg C2H4 eq and 838.39 MJ, respectively. The characterization index of global warming in the studies in Portugal, Ireland, France, Japan for one ton of milk were reported to be equal to 1021, 950.85, 1060, 980 and 1230 kg CO2 eq,

respectively. This impact category contained three important GHG such as CH4, CO2 and N2O. In most studies, the main reasons for its emission were enteric fermentation, diesel fuel and the management of manure [15,16,18,20,35]. The characterization index of acidification for one ton of milk in countries including Portugal, Ireland, France, Japan, and Germany, Italy and Denmark were reported 20, 9.40, 7.20, 7.13, 16.93 kg SO2 eq, respectively. Similar results were found in the current study which showed that SO2 and N2 emitted from concentrate were the main reasons for the emission of this impact category [15,16,18,20,35]. The results revealed that the amount of acidification (5.95 kg SO2 eq) in Guilan province of Iran, was relatively lower than the European countries. One of the main reasons for low amount of this impact category was the limited lands in Guilan. In the other hand, it is not possible to production of plant grass and grain products for animal feed in large scale. In the present study, the highest factor for emission of PO4-3 was concentrate which was equal to 0.14 kg for one ton of milk. Other inputs and activities lacked PO4-3. In some studies, the main reason for emission of this impact category was the emitted N2 from concentrate, corn silage and grass silage [16,20]. Nowadays, there are many solutions to decrease environmental impacts in milk production farms namely: animal feed management. To decrease and control GHG on dairy farms in the Netherlands, it was suggested to replace grass silage with corn silage for feed management. Applying this strategy can decrease negative environmental impacts resulted from enteric fermentation and also can decrease greenhouse gases as 12.8 kg CO2 eq yearly for one ton of milk production [46]. There are some other studies in which adopting the best management strategy in animal feed has been introduced as an appropriate way to decrease negative environmental impacts in dairy farms [25,27]. The index of characterization for photochemical oxidation in this study was 0.21 kg C2H4 eq which was less than the index of dairy farms of

Ireland and Japan [35,18]. The main elements of this group were CH4, SO2 and CO and the main reasons of their emission were enteric fermentation, concentrate, wheat and rice straw. The amount of fossil fuel depletion which contained the energy provided from the consumption of diesel fuel and electricity in Guilan's milk production units was 838.39 MJ which was less than the other countries [15,16,18,20,35]. The reason for the low amount of this impact category is inadequate number of cows in livestock units due to high price of animal feed. This has led to reduction of fuel consumption in livestock units. In most of the developed countries biogas technology is used to manage the pollutions of emission from fossil fuels in dairy farms more efficiently. This is one of the basic methods to control and decrease the GHG emission and it has an economic benefit [21,47].

Table 6. Life cycle assessment of Guilan dairy farms for the production of one ton of milk

Fig. 3 shows the final indices for each impact category; global warming, acidification, eutrophication, fossil fuel depletion and photochemical oxidation which were 0.24, 0.28, 0.076, 0.046 and 0.017 respectively. Acidification and global warming in environmental impact categories (ExoX), had the most disastrous impacts on the environment. The impact of acidification in this study was high due to high consumption of concentrate and electricity which were mostly consumed for heating water and pumps in milk equipment. In addition, animal feed, enteric fermentation, electricity and diesel fuel for tractors were introduced as the main reasons for high global warming.

Fig. 3. The final indices of the environmental impacts of milk production in Guilan province, Iran

4

Conclusion

Results revealed that CO2, CH4 and NO3 caused the most environmental pollutions in most impact categories. The main potential inputs for high CO2 emission in dairy farms was animal feed, electricity and diesel fuel which were respectively 676.3, 38.93 and 24.98 kg for one ton of milk production. Among animal feed, concentrate was selected as the main factor of GHG emissions. Acidification and global warming caused the most amount of damage to the environment. The main point of these environmental impacts were animal feed particularly concentrate, electricity, enteric fermentation and diesel fuel. Other impact categories were not significantly different from each other in case of pollution. In this study, to better environmental management of milk production in Guilan province of Iran, implementation of animal feed management programs and use of biogas systems have been proposed. These systems have not been used widely in Iran due to high costs and lack of adequate knowledge. In addition, it is recommended to take major steps through facilitate livestock units establishment and also encouraging the managers to be familiarize with new technologies like biogas.

Acknowledgment This research was supported by Ferdowsi University of Mashhad, Iran. Further, we are indebted to Surur Khorramdel (Department of Agronomy, Ferdowsi University of Mashhad, Iran) for her critical comments on the manuscript.

References [1] Food and Agriculture Organization of the United Nations (FAO).World Agriculture: towards 2030/2050. Interim Report, Rome, Italy;2006a [2] Van der Maas CWM, Coenen PWHG, Zijlema PJ, Brandes LJ, Baas K, Van den Berghe G, Van den Born GJ, Guis B, Geilenkirchen G, Te Molder R, Nijdam DS, Olivier, JGJ, Peek CJ, Van Schijndel MW, vander Sluis SM. Greenhouse Gas Emissions in the Netherlands 1990–2007 National Inventory Report 2009. Netherlands Environmental Assessment Agency (PBL), Bilthoven; 2009. p. 210. [3] Smith P, MartinoD, CaiZ, Gwary D, Janzen H, KumarP, McCarl B, Ogle S, O’Mara F, Rice C, Scholes B, Sirotenko O. Agriculture. In Climate Change 2007: Mitigation. Contribution of Working Group III to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change [Metz B, Davidson OR, Bosch PR, Dave R, Meyer LA, (eds)], Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA; 2007. [4] Stanford University. Environmental and Social Impact of the Livestock revolution Retrieved from: Available from: http://www.science daily.com. 2010. [5] Food and Agriculture Organization of the United Nations (FAO). Livestock’s long shadowenvironmental issues and options, Food and Agriculture Organization of the United Nations; 2006b. [6] Food and Agriculture Organization of the United Nations (FAO). Greenhouse Gas Emissions from the Dairy Sector. A Life Cycle Assessment, Food and Agriculture Organization of the United Nations; 2010.

[7] Christie KM, Rawnsley RP, Eckard RJ. A whole farm systems analysis of greenhouse gas emissions of 60 Tasmanian dairy farms. Anim Feed Sci and Tech 2011;166-167: 653-662. [8] Guinée JB, Gorrée M, Heijungs R, Huppes G, Kleijn R, De Koning A, Van Oers L, Wegener Sleewijk A, Suh S, Udo De Haes HA, De Brunij H, Huijbregts MAJ, Lindeijer E, Roorda AAH, Van der Ven BL, Weidema BP. Handbook on Life Cycle Assessment; Operational Guide to the ISO Standards. Kluwer Academic Publishers, Dordrecht; 2001. [9] Bojacá CR, Wyckhuys KAG, Schrevens E. Life cycle assessment of Colombian greenhouse tomato production based on farmer-level survey data. J Cleaner Produ 2014; 69: 26-33. [10] Khoshnevisan B, Rafiee SH, Mousazadeh H. Environmental impact assessment of open field and greenhouse strawberry production. European J Agro 2013; 50: 29-37. [11] Khoshnevisan B, Rajaeifar MA, Clark S, Shamahirband S, Anuar NB, Mohd NL, Shuib A, Gani A. Evaluation of traditional and consolidated rice farms in Guilan province, Iran, using life cycle assessment and fuzzy modeling. Sci Total Envir 2014; 481: 242-51. [12] Sabzevari A, Kouchaki-Penchah H, Nabavi-Pelesaraei A. Investigation of life cycle assessment of hazelnut production in Guilan province of I.R. Iran based on orchards size levels. Biological Forum – An International Journal. 2015; 7 (1): 807-813. [13] Meier MS, Stoessel F, Jungbluth N, Juraske R, Schader C, Stolze M. Environmental impacts of organic and conventional agricultural products – Are the differences captured by life cycle assessment? J Envir Manage 2015; 149: 193-208.

[14] Thomassena MA, Van-Calkerb KJ, Smitsc MCJ, Iepemad GL, de IJM, Boera. Life cycle assessment of conventional and organic milk production in the Netherlands. Agricultural Systems 2008; 96 (1–3): 95–107 [15] Van der Werf HMG, Kanyarushoki C, CorsonMS. An operational method for the evaluation of resource use and environmental impacts of dairy farms by life cycle assessment. J Envir Manage 2009; 90: 3643–3652. [16] Castanheira ÉG, Dias AC, Arroja L, Amaroc R. The environmental performance of milk production on a typical Portuguese dairy farm. Agri Sys 2010; 103: 498–507. [17] McGeough EJ, Little, SM, Janzen HH, McAllister TA, McGinn SM, Beauchemin KA. Life-cycle assessment of greenhouse gas emissions from dairy production in Eastern Canada: A case study. J Dairy Sci (JDS) 2012; 95: 5164–5175. [18] O’Brien D, Shalloo L, Patton J, Buckley F, Grainger C, Wallace M. A life cycle assessment of seasonal grass-based and confinement dairy farms. Agri Sys 2012; 107: 33–46. [19] Nutter DW, Kim DS, Ulrich R, Thoma G. Greenhouse gas emission analysis for USA fluid milk processing plants: Processing, packaging, and distribution. Int Dairy J 2012; 31: 57-64. [20] Guerci M, Knudsenb MT, Bava L, Zucali M, Schönbach PH, Kristensen T. Parameters affecting the environmental impact of a range of dairy farming systems in Denmark, Germany and Italy. J Cleaner Prod 2013; 54: 133-141. [21] Zhang S, Bi-Tony X, Clift R. A Life cycle assessment of integrated dairy farm-greenhouse systems in British Columbia. Bio Tech 2013; 150: 496–50.

[22] Beukes PC, Gregorini P, Romera AJ, Levy G, Waghorn GC. Improving production efficiency as a strategy to mitigate greenhouse gas emissions on pastoral dairy farms in New Zealand. Agri, Ecosys & Envir 2010; 136 (3–4): 358-365. [23] Christie KM, Rawnsley RP, Eckard RJ. A whole farm systems analysis of greenhouse gas emissions of 60 Tasmanian dairy farms. Anim Feed Sci and Tech 2011; 166-167: 653-662. [24] Flysjö A, Henriksson M, Cederberg C, Ledgard S, Englund JE. The impact of various parameters on the carbon footprint of milk production in New Zealand and Sweden. Agri Sys 2011; 104: 459–469. [25] Gerber P, Velling T, Opio C, Steinfeld H. Productivity gains and greenhouse gas emissions intensity in dairy systems. Livestock Sci 2011; 139 (1-2): 100-108. [26] Hensen A, Groot TT, vanden Bulk WCM, Vermeulen AT, Olesen JE, Schelde K. Dairy farm CH4 and N2O emissions, from one square meter to the full farm scale. Agri, Ecosys & Envir 2006; 112 (2–3): 146-152. [27] Prado AD, Mas K, Pardo G, Gallejones P. Modelling the interactions between C and N farm balances and GHG emissions from confinement dairy farms in northern Spain. Sci Total Envir 2013; 465: 156-165. [28] The Iranian center of Statistics Centre of Iran. Available at amar.org.ir. 2014. [29] The Iranian center of Statistics. Annual Statistics of Guilan, Iran. Statistics related to agricultural units and natural resources; 2013. [30] Snedecor GW, Cochran WG. Statistical methods. Iowa State University Press; 1980.

[31] International Organisation for Standardisation (ISO). Environmental Management – Life Cycle Assessment: Requirements and Guidelines (ISO 14044:2006). European Committee for Standardization, Brussels, Belgium; 2006a. [32] International Organisation for Standardisation, (ISO). Environmental Management – Life Cycle Assessment: Principles and Framework (ISO, 14040:2006). European Committee for Standardization, Brussels, Belgium; 2006b. [33] Basset-Mens C, Ledgard S, Boyes B. Eco-efficiency of intensification scenario for milk production in New Zealand. Ecological Econ 2009; 68: 127-144. [34] Ecoinvent Centre.Ecoinvent 2.0 database. Swiss centre for life cycle inventories, Dübendorf. http://www.ecoinvent.ch (accessed 30.03.11). 2010. [35] Ogino A, Ishida M, Ishikawa T, Ikeguchi A, Waki M, Yokoyama H, Tanaka Y, Hirooka H. Environmental impacts of a Japanese dairy farming system using whole-crop rice silage as evaluated by life cycle assessment. Anim Sci J 2008; 79: 727-736. [36] Department of Environment (DOE). Iran’s initial national communication to UNFCCC, Climate Change Office Publication. Islamic Republic of Iran; 2003. p.147 [37] Brentrup F, Küsters J, Kuhlmann H, Lammel J. Environmental impact assessment of agricultural production systems using the life cycle assessment methodology: I. Theoretical concept of a LCA method tailored to crop production. European J Agro 2004a; 20: 247-264. [38] Nikkhah A, Khojastehpour M, Emadi B, Taheri-Rad A, Khorramdel S. Environmental impacts of peanut production system using life cycle assessment methodology. J Cleaner Prod 2015; 92: 84-90.

[39] Fallahpour F, Aminghafouri A, Ghalegolab-Behbahani A, Bannayan M. The environmental impact assessment of wheat and barley production by using life cycle assessment (LCA) methodology. Envir, Develop and Sust 2012; 14: 979-992. [40] Intergovernmental Panel on Climate Change (IPCC).Climate Change 1995. The Science of Climate Change. Contribution of Working Group I to the Second Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge University Press, Cambridge, UK; 1996. [41] Derwent RG, Jenkin ME, Saunders SM, Pilling MJ. Photochemical ozone creation potentials for organic compounds in Northwest Europe calculated with a master chemical mechanism. Atmospheric Env1998; 32: 2429-2441. [42] Mirhaji H, Khojastehpour M, Abaspour-Fard M. H. Environmental Effects of wheat production in the Marvdasht region. Journal of Natural Environment. J.Natural Envir 2013; 66 (2): 223-232 (In Persian). [43] Thoma G, Popp J, Shonnard D, Nutter D, Matlock M, Ulrich R, Kellogg W, Kim DS, Neiderman Z, Kemper N, Adomd F, East C. Regional analysis of greenhouse gas emissions from USA dairy farms: A cradle to farm-gate assessment of the American dairy industry circa 2008. Int Dairy J 2013; 31: 29-40. [44] Grainger C, Clarke T, Beauchemin KA, McGinn SM, Eckard RJ. Effect of whole cottonseed supplementation on energy and nitrogen partitioning and rumen function in dairy cattle on a forage and cereal grain diet. Australian J Exp Agri 2008; 48: 860-865.

[45] McGeough EJ, Little, SM, Janzen HH, McAllister TA, McGinn SM, Beauchemin KA. Lifecycle assessment of greenhouse gas emissions from dairy production in Eastern Canada: A case study. J Dairy Sci (JDS). 2012; 95: 5164–5175. [46] Van Middelaar CE, Cederberg C, Vellinga TV, Van Der Werf HMG, De Boer IJM. Exploring variability in methods and data sensitivity in carbon footprints of feed ingredients. Int J Life Cycle Asses 2013; 18: 768–782. [47] Whiting A, Azapagic A. Life cycle environmental impacts of generating electricity and heat from biogas produced by anaerobic digestion. Energy 2014; 70: 181–193.

Figure captions Fig. 1. Inputs and outputs of milk production in Guilan province, Iran (system boundary) Fig. 2. The total emissions of pollutants (kg) in industrial dairies in Guilan province, Iran Fig. 3. The final indices of the environmental impacts of milk production in Guilan province, Iran

Table captions Table 1. The inputs and output to produce one ton of milk in Guilan province, Iran Table 2. Emission coefficients of contaminants resulting from the use of inputs and activities Table 3. Characterization factor of environmental impacts Table 4. Normalization and weighting factors Table 5. The emissions of contaminants resulting from the consumption of inputs and activities for each ton of milk production Table 6. Life cycle assessment of Guilan dairy farms for the production of one ton of milk

Table 1 The inputs and outputs to produce one ton of milk in Guilan province of Iran Inputs & Outputs Unit Mean Standard Deviation (std) Diesel Fuel

lit t-1

9.40

6.60

Electricity

kWh t-1

67.13

12.58

Natural Gas

m3t-1

7.86

6.56

Animal feed

kg t -1

1311.42

437.88

Manure

kg t-1

458.50

42.21

Milk

t head-1 year-1

7.34

1.40

Table 2 Emission coefficients of contaminants resulting from the use of inputs and activities Inputs Diesel Fuel (kgl-1) Natural Gas (kg m3 -1) Electricity (kg kWh-1) Feed (kg kg-1) -Concentrate -Silage -Straw -Alfalfa Enteric Fermentation (GgYear-1 ) Manure Management ( GgYear-1)

PO4

NO3

CO

Emission of pollutants (kg per unit) SO2 NH3 CO2 NOx

-

10-5

-

10-5×34

-

2.63

105×132

-

-

-

-

-

133×10-3

-

-

-

-

-

10-5×27

0.03112

-

-

-

10-3×9

-

-

-

-

-

-

10 5 ×276 105 ×112 10-4× 8 10-5×9

N2O 105 ×181

CH4 105 ×143

7×10-6

26×10-4

[34]

-

[34]

Reference

-

10-5

0.580

10-3×5

0.70

10 5 ×121

-

10 ×21 10-4

10 ×187 2.8 10-4×229

105 ×217 10-6 0.04 10-5×56

-

-

-

-

-

436.53

[36]

-

-

-

-

21.33

23.73

[36]

-5

-4

10-5×73 -5

10 ×2 10-4×3 10-5×7

10×147 10-4 -

5

[34] [35] [35] [35]

Table 3 Characterization factor of environmental impacts Impact Category Global Warming Potential (kg CO2 eq) Acidification (kg SO2 eq) Eutrophication (kg PO43- eq) Depletion of fossil resources (MJ) Photochemical oxidation (kg C2H4 eq)

Potential of compounds CO2=1, CH4=21, N2O=310 SO2=1.2, NOX=0.5, NH3=1.6 NOX=0.13, NH3=0.35, NO3= 0.1, PO4=1 42.86 CH4 = 0.006, SO2= 0.048, CO = 0.027

Reference [38,40] [37] [37] [37] [41]

Table 4 Normalization and weighting factors Impact Category Global Warming Potential (kg CO2 eq) Acidification (kg SO2 eq) Eutrophication (kg PO43- eq) Depletion of fossil resources (MJ) Photochemical oxidation (kg C2H4 eq)

Weighting factor

Reference

8143

1.05

[38,42]

52

1.8

[38,42]

63

1.4

[38,42]

39167

1.14

[37]

9.69

0.8

Normalization factor (unit)

[8]

Table 5 The emissions of contaminants resulting from the consumption of inputs and activities for each ton of milk production Inputs Diesel Fuel Electricity Natural Gas Animal feed: -Concentrate -Silage -Straw -Alfalfa Enteric Fermentation Manure Management

CH4 0.013

pollutions (kg per one ton of milk) N 2O CO SO2 NH3 0.017

0.032 0.18

0.020

0.000023

0.79

0.39 0.009 0.034 0.014

0.008 26.57 0.66

1.14

-6

10 ×6.7

CO2

NOX

NO3

24.98

0.012

0.0066

38.93

0.08

PO4

1.04 0.60 0.77

0.068 0.020

2.70 0.1

378.18 8.87 241.10 48.15

1.17 0.005 3.44 0.11

16.80

0.14

Table 6 Life cycle assessment of Guilan dairy farms for the production of one ton of milk Characterization Impact Category Normalization Index Final Index index t-1 t-1 t-1 Global Warming 1831(kg CO2 eq) 0.22 0.24 Potential Acidification 7.97 (kg SO2 eq) 0.15 0.28 -3 Eutrophication 3.42 (kg PO4 eq) 0.054 0.076 Depletion of fossil 838.39(in MJ) 0.040 0.046 resources Photochemical 0.21(kg C2H4 eq) 0.022 0.017 oxidation

Highlights • • • • •

Iranian milk production was environmentally assessed with a life cycle approach. Five impact categories are assessed. The highest potential for environmental impacts of production revealed for acidification. Environmental index (EcoX) and resources depletion index (RDI) were obtained 0.61 and 0.04, respectively. Implementation of animal feed management programs and use of biogas systems have been proposed.