Next Article in Journal
Studying Industrial Decarbonisation: Developing an Interdisciplinary Understanding of the Conditions for Transformation in Energy-Intensive Natural Resource-Based Industry
Next Article in Special Issue
Systematic Mapping of Research on Farm-Level Sustainability in Egg and Chicken Meat Production
Previous Article in Journal
Solving the Multi-Depot Green Vehicle Routing Problem by a Hybrid Evolutionary Algorithm
Previous Article in Special Issue
Silage Fermentation and In Vitro Degradation Characteristics of Orchardgrass and Alfalfa Intercrop Mixtures as Influenced by Forage Ratios and Nitrogen Fertilizing Levels
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effect of Feed Concentrate Intake on the Environmental Impact of Dairy Cows in an Alpine Mountain Region Including Soil Carbon Sequestration and Effect on Biodiversity

1
Faculty of Science and Technology, Free University of Bozen-Bolzano, Piazza Università 5, 39100 Bolzano, Italy
2
Research Institute of Organic Farming, Kasseler Strasse 1a, 60486 Frankfurt am Main, Germany
*
Authors to whom correspondence should be addressed.
Sustainability 2020, 12(5), 2128; https://0-doi-org.brum.beds.ac.uk/10.3390/su12052128
Submission received: 20 February 2020 / Revised: 3 March 2020 / Accepted: 4 March 2020 / Published: 9 March 2020
(This article belongs to the Special Issue Sustainable Livestock Production)

Abstract

:
Several studies on the environmental impacts of livestock enterprises are based on the application of life cycle assessments (LCA). In Alpine regions, soil carbon sequestration can play an important role in reducing environmental impacts. However, there is no official methodology to calculate this possible reduction. Biodiversity plays an important role in the Alpine environment and is affected by human activities, such as cattle farming. Our aim was to estimate the carbon footprint (CF) of four different dairy production systems (different in breeds and feeding intensity) by using the LCA approach. The present study included 44 dairy Alpine farms located in the autonomous province of Bolzano in northern Italy. Half of the farms (n = 22) kept Alpine Grey and the other half (n = 22) Brown Swiss cattle. Within breeds, the farms were divided by the amount of concentrated feed per cow and day into high concentrate (HC) and low concentrate (LC). This resulted in 11 Alpine Grey low concentrate (AGLC) farms feeding an average amount of 3.0 kg concentrated feed/cow/day and 11 Alpine Grey high concentrate (AGHC) farms with an average amount of 6.3 kg concentrated feed/cow/day. Eleven farms kept Brown Swiss cows with an average amount of 3.7 kg concentrated feed/cow/day (BSLC) and another 11 farms feeding on average 7.6 kg concentrated feed/cow/day (BSHC). CF for the four systems was estimated using the LCA approach. The functional unit was 1 kg of fat and protein corrected milk (FPCM). Furthermore, two methodologies have been applied to estimate soil carbon sequestration and effect on biodiversity. The system with the lowest environmental impact in terms of CF was BSHC (1.14 kg CO2-eq/kg of FPCM), while the most impactful system was the AGLC group (1.55 kg CO2-eq/kg of FPCM). Including the CF reduction due to soil carbon sequestered from grassland, it decreased differently for the two applied methods. For all four systems, the main factor for CF was enteric emission, while the main pollutant was biogenic CH4. Conversely, AGLC had the lowest impact when the damage to biodiversity was considered (damage score = 0.41/kg of FPCM, damage to ecosystem diversity = 1.78 E-07 species*yr/kg FPCM). In comparison, BSHC had the greatest impact in terms of damage to biodiversity (damage score = 0.56/kg of FPCM, damage to ecosystem diversity = 2.49 E-07 species*yr/kg FPCM). This study indicates the importance of including soil carbon sequestration from grasslands and effects on biodiversity when calculating the environmental performance of dairy farms.

1. Introduction

With a share of about 18%, livestock farming is one of the strongest contributors to greenhouse gas emissions [1]. Methanogenesis in the rumen is an essential metabolic process to maintain steady state fermentation as this scavenges the molecular hydrogen generated during fermentation. The methane production from enteric fermentation in animal agriculture contributes around 20% of the total global methane [2]. Due to an overall increase in milk production by 30% and a global increase in the number of dairy cows by 11%, the dairy sector’s greenhouse gas (GHG) emissions have increased by 18 percent between 2005 and 2015 [3]. According to the Intergovernmental Panel on Climate Change (IPCC), the main GHG emissions attributed to the agricultural sector are methane (CH4) and nitrous oxide (N2O). Life cycle assessment (LCA) is one of the most common approaches to estimate the environmental impact of a product for its entire life cycle. Several studies used LCA to estimate the impact of dairy cattle production [4,5], beef production [6,7], dairy buffalo farming [8,9,10], and dairy goat farming [11,12]. Carbon footprint (CF) is entirely referred to global warming potential (GWP) and it is commonly used to communicate the contribution of dairy production to climate change to stakeholders [13]. Dairy farming in the Italian Alps has transformed greatly in the last decades; between 1990 and 2010, farms declined by 60%, while the number of cows decreased by 29%, with a growth of the average herd size by 76% [14,15,16]. Many of these farms have abandoned traditional Alpine pasture practices by increasing the amount of concentrated feed purchased off-farm. In the Alpine region of the autonomous province of Bolzano, the indigenous breeds Alpine Grey and Brown Swiss are commonly used. Both abandonment of traditional farms or their transformation into intensive enterprises leads to a loss of open areas and forest re-growth [17], a damage of biodiversity [18,19,20], and a rural area undergoing radical socio-economic changes [21,22,23]. Only a few studies took soil carbon sequestration from grassland into account [24,25,26]. These grasslands may act as a carbon sink, with soil carbon sequestration of about 1 t CO2 ha−1 year−1 having been described [27]. However, there is currently no distinct approach to which methodology should be used to determine soil carbon in livestock farming systems. Another important point to consider for a proper assessment of the environmental impact generated by livestock production systems is the effect on biodiversity [28]. The damage score approach is based on the use of the characteristic factors for ecosystem damage of different land uses and agricultural practices [29]. The ReCiPe endpoint methodology is a life cycle impact assessment developed by [30] that provides harmonized characterization factors at midpoint and endpoint levels. At the endpoint level, 17 endpoint impact classes are identified and then pooled to 3 endpoint damage categories labelled as ‘human health’ damage, ‘ecosystems’, and ‘resources’ [30,31]. This method has been used only in a few studies on dairy farms [32,33,34], and it has never been used for small-scale Alpine farms. In the present study, both of the mentioned methodologies were applied to assess the effect of small-scale dairy farms on biodiversity in an Alpine environment. The aim of the present study was to investigate the effect of the used breed and the amount of concentrated feed on the carbon footprint, with and without soil carbon sequestration, and the effect on biodiversity.

2. Materials and Methods

2.1. Experimental Design and Data Collection

The present study involved 44 dairy Alpine farms located in the autonomous province of Bolzano in northern Italy. Farms were conducted from October 2017 until September 2018. Half of the farms (n = 22) kept Alpine Grey and the other half (n = 22) Brown Swiss cattle. Within breeds, the farms were divided by the amount of concentrated feed per cow and day into high concentrate (HC) and low concentrate (LC). This resulted in 11 Alpine Grey low concentrate (AGLC) farms, feeding an average amount of 3.0 kg concentrated feed/cow/day and 11 Alpine Grey high concentrate (AGHC) farms with an average amount of 6.3 kg concentrated feed/cow/day. Eleven farms kept Brown Swiss cows with an average amount of 3.7 kg concentrated feed/cow/day (BSLC) and another 11 farms feeding on average 7.6 kg concentrated feed/cow/day (BSHC). All farms were run as family businesses. The average herd size was 12 lactating cows/farm and farms were located on average at an altitude of 1000 meters above sea level. Primary data were obtained directly from the farmers by providing a questionnaire. Table 1 shows the main characteristics of the farm groups. Using commercial software [35], the amounts of raw materials, as shown in Table 2, contained in the concentrated feed used by the farms were estimated. Based on the labeling of the individual concentrated feed and their basic analytical composition (e.g., crude protein, crude fiber) and energy composition, the percentage of raw material contained in the concentrated feed was estimated with an error of approximately ± 10%. Crop yield was estimated by using the production crop database 2017 from Italy [36]. Table 2 shows in detail the main ingredients and their chemical composition in the concentrated feeds used on the farms.

2.2. Method for Carbon Footprint Calculation

By using the LCA approach, CF values were calculated. The CF was assessed using the software package SimaPro 8.01 PhD. In particular, the ReCiPe method midpoint/endpoint (H) V1.10/Europe Recipe H/A module of this package was used. CF is the net GWP emission per production unit and needed to express the functional units (FU). To calculate the GWP, the LCA approach suggested by [37] was used, and in addition the indications from [1].

2.2.1. Functional Unit and System Boundaries

Following the indications given by the Food and Agriculture Organization and the International Dairy Federation [38] for the evaluation of the CF of dairy production, 1 kg of fat and protein corrected milk (FPCM) was used as the FU:
1 kg FPCM = raw milk (kg) * (0.337 + 0.116 * fat content (%) + 0.06 * protein content (%)).
System boundaries from “cradle to farm-gate” included all farm activities (on-farm and off-farm), except machinery construction, building construction, and medicines.

2.2.2. Emission Calculations for Carbon Footprint Assessment

The total emissions were estimated for each farm system including enteric fermentation, manure management, fuel and electric consumption, permanent pasture management, meadow hay production, and also the off-farm production of concentrated feed. For estimations of CH4 enteric fermentation and manure emission and management, we applied the equation (10.21) Tier 2 method suggested by [1]:
EF = [GE * (Ym/100) * 365 / 55.65];
where:
EF = emission factor, kg CH4 head−1 y−1;
GE = gross energy intake, MJ head−1 day−1;
Ym = methane conversion factor, % of gross energy in feed converted to methane (4% for young cattle and 6% for lactating cows, [39]).
For estimation of GE, the equation (10.16) suggested by [1] was used.
For estimations of the dry matter intake (DMI) for mature dairy cows, the equation (10.18b) suggested by [1] and [40] was used:
DMI = [((5.4 * BW)/500)/((100-DE%)100))];
where:
BW = live body weight, kg (Alpine Grey was 550 kg; Brows Swiss was 650 kg);
DE% = digestible energy expressed as a percentage of gross energy was 60%.
According to the method suggested by [41] for the estimation of CH4 emission from manure management, we used the equation (10.22). For calculation of direct N2O (kg y−1) emission from manure management, we used the equation (10.25) suggested by [41]. For this equation, the emission factor was 0.005 for solid storage and 0.01 for pasture kg N2O-N/kg N N2O. Indirect emissions of N2O were estimated according to the method suggested by [41], as shown in equation (10.27), based on NO3 leaching-runoff and the emission factor for N2O emissions from atmospheric deposition of nitrogen on soils and water surfaces, kg N2O-N (kg NH3-N + NOx-N volatilized) −1 default value was 0.01 kg N2O-N (kg NH3-N + NOx-N volatilized) −1. The emission generated from diesel and electric consumption for agricultural practice were estimated from farmers’ information, as shown in Table 1. In particular, we used 0.85 kg per liter as a standard value for diesel density. The emission factor for Italy was 3.13 kg CO2-eq for burning 1 kg of diesel [9,25]. Regarding electricity, the emission factor was set at 0.47 kg CO2-eq per kWh using national emission factors [42]. The database from [43] was used to estimate the impact of concentrate production off-farm. The fuel consumption for the delivery of the concentrated feed to the farms was also considered.

2.3. Soil Carbon Sequestration

Soils are the largest carbon reservoir of the terrestrial carbon cycle. However, there is no official method to estimate soil carbon for livestock farms and only a few studies have investigated this factor [24,25]. Different agricultural practices and land use can have significant effects on carbon stored in the soil. In the present study, two different methodologies were used to estimate the reduction of carbon footprint by taking carbon sequestration into account.

2.3.1. Soil Carbon Sequestration

Soussana et al. [44] approached this method based on the net carbon fluxes. These authors studied the role played by the soil and pasture in sequestrating carbon as a sink and mitigating its effects. However, the emissions of CH4 and N20 should be reduced while preserving soil carbon storage. In particular, agricultural practices that prevent carbon sinks should be avoided. Grassland management implies a high level of uncertainty. In our study, we used the mean value reported by [44], corresponding to 22 g C/m2 for our estimates.

2.3.2. Soil Carbon Sequestration

The methodology suggested by [45] is based on a 100-year perspective and is when 10% of the total carbon added to the soil will be sequestered. These authors also give relevance to the potential played by soils for carbon sequestration as a suitable strategy to reduce GHG emissions. As suggested by [44,45], carbon change estimates are based on net carbon fluxes. Herbage residues and manure were included in the computation of the annual inputs of carbon into grassland. The former was calculated according to [24]. These authors considered above- and below-ground residues at 40% and 16% of total yield, respectively, and assumed a 45% content of C on a dry matter basis.

2.4. Damage on Biodiversity

Biodiversity can be defined as the genetic variability among any organisms living in the marine environment, in other aquatic systems, on land, or in the ground. Biodiversity, therefore, concerns the diversity in terms of ecosystems, between species and within species. Agricultural practices and the consequent land use may markedly damage biodiversity. In this study, we used two different methodological approaches to evaluate biodiversity damage.

2.4.1. Biodiversity Damage Score

The loss of biodiversity was estimated as suggested by [29,46]. This method is based on the use of characterization factors (CFs). The damage score (DS) was used to express the impact on biodiversity. This score measures the relative change in the number of species detected in the occupied area as compared with the baseline. Table 3 shows the CFs used to evaluate biodiversity damage in this study. The local DS was estimated for 1 m2 of farmland and then related to 1 kg of FPCM.
DS = CFs * A * t
where:
DS = damage score;
CFs = characterization factors, as shown in Table 3;
A = area occupied (m2);
t = time (years).

2.4.2. Damage to Ecosystem Diversity

The damage to ecosystem diversity is expressed in number of species per year and represents the loss of species over a specified period of time and in a given geographical local area. Carried out on commercial software SimaPro 8.01 PhD, the ReCiPe endpoint method approached the damage on ecosystems and biodiversity per 1 kg FPCM. This is a life cycle impact assessment developed by [30] that provides harmonized characterization factors at midpoint and endpoint levels. At endpoint level, 17 endpoint impact classes were identified, and damage categories labelled as ‘ecosystems’ were studied.

2.5. Statistical Analysis

The data at farm level, carbon footprint, main processes, main pollutants, and indicator of damage to biodiversity were analyzed by using one-way (general linear model procedure) ANOVA, using the farming systems as a factor. A correlation analysis was used to examine the relationship between the two methods of damage score estimation. The farm was used as an experimental unit. Therefore, the correlation coefficients were calculated using either Pearson tests according to the data distribution of the variables. Data were analyzed using SAS software (SAS Institute Inc., Cary, NC). Significance was accepted at P < 0.05.

3. Results and Discussion

3.1. Farm Levels

Table 1 shows the main characteristics of the four farming systems. The number of lactating cows averaged 12 without a significant difference between the four systems. It could be observed that small scale dairy farms in the Alpine area were mainly family businesses. Salvador et al. [25] observed an average number of lactating cows between 8.7 and 40.6 in the Alpine systems, whereas [47] found an average of 14.0 lactating dairy cows per farm, and [15] observed on average 12.8 dairy cows per dairy farm in the Italian Alps. In terms of milk production, a significant difference was found between the systems, as shown in Table 2. The AGLC system shows significantly lower kg of FPCM per cow per year compared to all other three systems (4206 ± 231 kg of FPCM/cow/year), whereas the BSHC system shows higher milk production of FPCM than any other system (7356 ± 231 kg of FPCM/cow/year). The differences observed were most probably due to a breed effect (Alpine Grey vs. Brown Swiss). Furthermore, several studies observed an increase in milk production due to an increased amount of concentrated feed [48,49,50]. For both breeds considered in this study, an average increase of about 50% in the quantity of concentrated feed was accompanied by an increase of on average 30% in the milk quantity, as shown in Table 2. In terms of land use for hay production, no significant difference was found between the groups, while in the AGLC system there was a significantly larger use of permanent pasture grazing resources compared to the BSHC system, and also the average number of days on pasture was higher (P < 0.01). The breeding of the Alpine Grey dairy cow was at one of the oldest and most traditional farms in the province of Bolzano, and most probably the change to other breeds of dairy cows and the increased use of commercial feed have led to a gradual reduction of use of the pasture resource. No significant difference was found for the general consumption of electricity and diesel.

3.2. Carbon Footprint Without Soil Carbon Sequestration

The carbon footprint generated from the four different systems per kg of FPCM is shown in Table 4. The system AGLC has the highest CF, differing significantly from all other systems. The AGLC system had an impact of about 26% greater than the BSHC system (P < 0.001), and about 21% greater than the AGHC system (P < 0.01), while no significant difference was found compared to the BSLC system. The AGHC system showed significantly less impact than the BSLC system (1.22 vs. 1.43 kg CO2-eq per kg FPCM respectively, P < 0.01), whereas no significant difference was observed between the AGHC and the BSHC system. In terms of the global warming potential, recent studies conducted in the Alpine region showed an environmental impact between 1.94 and 1.59 kg CO2-eq per kg FPCM [25], while [32] observed 1.60 kg CO2-eq per kg FPCM. Similar results were observed by [51] in southern Germany (1.53 kg CO2-eq/kg FPCM). Other authors showed lower values in the Alpine area similar to ours. Penati et al. (2013) [52] observed values of 1.14 kg CO2-eq per kg FPCM, whereas [53] in Switzerland observed results of 1.08 kg CO2-eq per kg FPCM. The results obtained from several studies were all dependent on the production of milk. From an environmental point of view, a system that produces more milk is more efficient. The BSHC system showed an average milk production of 7300 kg FPCM per cow, while the least productive system (AGLC) achieves an average milk yield of 4200 kg FPCM per cow and lactation. Most probably, the differences observed compared to other studies conducted in the same geographical area were due to the different milk production levels from different breeds studied and differences in management practices. Nevertheless, our results are in the ranges indicated by [3]. The main processes contributing to CF are shown in Table 5. The enteric emission from different systems was the process with the highest impact, followed by the production of concentrated feed and meadow hay. The two systems with lower amounts of concentrated feed (AGLC and BSLC), showed a greater impact generated from enteric emissions than the two systems with a high use of concentrated feed. This is most probably due to the use of cows fed with a greater amount of raw fiber, which can generate a greater impact due to the process from enteric fermentation. In a recent report by [2], enteric fermentations were responsible for about 58.5% of total emissions, while [32] observed a range between 44.1% and 65.9%, and [54] observed a value of 48.4% in organic dairy cow farms. In cows, as in other ruminants, CH4 is emitted as a direct result of fiber hydrolysis. In particular, in the rumen a number of fermentative processes occur, resulting in the production of important nutrients, such as volatile fatty acids and proteins of bacterial origin deriving from the fermentation of vegetal fibers. Table 6 shows the main pollutants in the different systems. The AGLC system shows a high emission of CH4 biogenic compared to other systems. Several studies tried to understand how to reduce the impact generated by enteric fermentations from ruminants; [55] showed that concentrate supplementation on great quality pasture-only diets have the possibility to decrease CH4 emissions. Jiao et al. (2014) [48] observed a reduction in CH4 enteric emissions in dairy cows fed with an increased amount of concentrate feed in the diet, whereas [56] suggested several feeding strategies to minimize the emissions of enteric methane, such as an increase of the voluntary feed intake on permanent pasture. As reported by [57], estimation of in vitro methane production in buffalo and dairy cow feedstuffs are characterized by high structural glucides of low quality. They consequently result in a high CH4 production in both species. Studies indicate that the use of saponins in ruminant feed rations may reduce the production of biogenic methane [58,59,60], while [61,62,63] indicate that the inclusion of more digestible forages in ruminant diets may reduce emissions from the ruminant’s fermentation. This is in line with the equation (10.18b) suggested by [1]; the digestibility of feed ingested by ruminants plays a key role in the emission factor. The second main process responsible for CF was the off-farm production for concentrated feed, as shown in Table 5. This process was significantly higher in the BSHC system than in other systems; conversely, the AGLC and BSLC systems showed a lower amount of CF in this process (12.2% and 11.9%, respectively). This corresponds to the larger and lower use of purchased feed in systems. The two main pollutant compounds present for this process were CO2 and N2O, as shown in Table 6. Salvador et al. (2016) [64] observed that 14.4% of total CF was due to purchased feed in organic dairy cow rearing systems in the Alpine, while [65] in buffalo farms observed a proportion of 30%. Table 5 shows that AGLC systems produce more hay than the two systems with higher amounts of concentrates; this is most probably due to more permanent grassland present in the AGLC system, as shown in Table 1. No significant differences were observed in terms of general consumption between the different systems.

3.3. Carbon Footprint Including Soil Carbon Sequestration

Table 4 shows the results of the CF including the reduction due to carbon sequestered from the soil following the methodologies by [44,45]. The amount of CO2-eq reduction depends on the method used. The reduction observed varies from 12.2% to 14.2% with the methodology suggested by [44] and between 14.0% and 21.2% using the method suggested by [45]. In a study on dairy sheep farming, [24] observed a reduction between 2% and 43% when using the Soussana et al. [44] model and a reduction of 3% to 41% when using the Petersen et al. [45] model. On dairy farms, [25] observed an average reduction by 29.7%, while [33] stated that the reduction by soil carbon sequestration to be negligible. For both methodologies, the differences observed between the four systems were maintained. In particular, the two systems with the use of low concentrated feeds (AGLC and BSLC) show a greater reduction in terms of CF at the on-farm level than systems with high concentrated feed use, and this is most probably due to a greater use of local resources, such as extensive permanent pasture and meadow hay production in the region. This may indicate that the less resources used off-farm (such as concentrated feed), the greater the benefit from soil carbon sequestration. For both methodologies, BSHC was the most efficient system when it came to the reduction of CF by soil carbon sequestration. With the methodology suggested by Petersen et al. [45], the CF generated by the BSHC system was 0.98 kg CO2-eq per kg of FPCM, while the most impactful system remains the AGLC, with 1.33 kg CO2-eq (with the method in [44]) and 1.22 kg CO2-eq (with the method suggested by [45]). Without these methodological differences, the results show that there is a possibility for permanent soil carbon sequestration in these systems if the same amount of C was added to the grasslands every year. This should be included in the CF methodology and it is possible to add the consequence of soil carbon sequestration in the CF [45]. Recent studies have observed that high levels of organic matter in the soil may contribute to lower GHG emission in the atmosphere and mitigate climate change [66]. This highlights the importance of including soil carbon sequestration in LCA due to the climate mitigation of Alpine grazing systems, and the probability of having a lower CF in extensive systems when soil carbon sequestration is included, as [67] pointed out in their study for dairy systems in Ireland.

3.4. Effect on Biodiversity with Two Different Methods

The results from the indicators of damage on biodiversity in different farm systems per 1 kg of FPCM are reported in Table 7. For both methods considered, the BSHC system showed a significantly higher damage to biodiversity compared to the other systems. In BSHC, the value observed was 0.52 DS per kg of FPCM, while the AGLC system showed the lowest value with 0.41 DS per kg of FPCM (P < 0.001). Guerci et al. [32], in a study conducted in EU farms, observed a DS between 0.25 and 1.25, with lower values observed for organic farms. Battini et al. [29] showed a DS no-allocation between 0.75 and 2.40. The AGLC system has a lower use of concentrated feed. It is well known that the production of concentrated feed needs less land but is accompanied by a greater damage to biodiversity, so that the characterization factor used for intensive arable land was 0.79, while organic fertile grassland has been used with a value equal to −0.01, as shown in Table 3. Therefore, the AGLC system goes along with greater benefits to biodiversity. Tumisto et al. [46] observed that some other types of nature reserves or natural meadows may have a greater level of biodiversity than natural forests, whereas meadows that change into other crop production are a driver of habitat loss and thus have a documented negative effect on general biodiversity [68]. A recent study conducted in the Alpine environment [20] observed that the abandonment of pastures, mountain pastures, and advancing forest reduces the biodiversity of vertebrates living in mountain streams. The results show the same trend when analyzed with the ReCiPe endpoint methodology [30]. BSHC showed a greater impact in terms of damage to biodiversity with a damage score of 2.49 E-07, higher than other systems, as shown in Table 7. Kaenchan et al. [69] claimed that land use can have negative effects on the ecosystem including humans, whereas [70] stated that the preservation of extensive agricultural systems, including species-rich meadows, is crucial for the maintenance of the biodiversity in semi-natural management conditions. Du et al. [71] observed that the production of manually harvested sugarcane generated damage to the species of 3.11 E-07, whereas the estimated damage to ecosystem diversity caused by mechanically harvested sugarcane was 2.38 E-07. In a recent study, Sabia et al. [34] observed that for 1 kg of sheep milk using the ReCiPe endpoint approach, the impacts in terms of damage on biodiversity was equal to 3.29 E-07 species*year. Steinmann et al. [72] observed that this damage can be the result of the combination of land and energy footprints. In addition, [73] observed that agronomic practices may undermine both soil and water quality while also negatively affecting the natural biodiversity through, for example, N deposition in ecosystems. In order to promote a balance between damage and restoration (i.e., compensation), the monetization of the environmental impact is necessary to combine all estimations of the ecological damage into one single value [74]. In a recent study conducted at the European level on the impact of total food consumption on biodiversity using the ReCiPe endpoint approach, Crenna et al. [75] observed that the damage produced in terms of species*year was equal to 3.34 E-05 per capita consumption. Figure 1 shows that a relative positive strong correlation was found between the two methods (r = 0.73, P ≤ 0.001). The two methodologies are very similar and have different points of affinity. The methodology of DS approaches through characterization factors, while the ReCiPe endpoint method is the category ‘ecosystems’ expressed in number of species per year and represents the loss of species over a specified period of time and in a given geographical area. Both methodologies have a direct approach to land use, land use change, and local geographical effects. However, there is no official methodology on how to approach the damage to biodiversity generated by livestock farming systems and, in particular, by dairy cows.

4. Conclusions

This study has taken into account small-scale Alpine dairy farms with different amounts of concentrated feed and either the breed Alpine Grey or Brown Swiss. The BSHC system was the most environmentally efficient system in terms of carbon footprint per kg of FPCM. However, the two systems using the Alpine Grey breed showed satisfactory environmental performance even though the AGLC system was less efficient. For all four systems, the main factor responsible for the carbon footprint was the enteric emission. Thus far, there is no generally accepted methodology to estimate the contribution from soil carbon sequestration in LCA. We took into account the soil carbon sequestered from grassland, which reduced the carbon footprint and depended on the method used. For biodiversity damage, the least impacting system was AGLC, while the system with the highest damage score was BSHC for both applied methods. Strong positive correlations were found between the results of damage to biodiversity for the two methods. More studies are needed to define the correct approach and methodology to determine the reduction due to carbon sequestration in soils and effect on biodiversity.

Author Contributions

Conceptualization, E.S. and M.G.; methodology, E.S. and M.G.; software, E.S.; validation, E.S. and M.G.; formal analysis, E.S.; investigation, E.S. S.K. L.F. C.L. and M.G.; resources, M.G.; data curation, E.S. S.K. L.F. C.L. and M.G; writing—original draft preparation, E.S and M.G..; writing—review and editing, E.S. S.K. L.F. C.L. and M.G; visualization, M.G.; supervision, M.G; project administration, M.G.; funding acquisition, M.G. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the “Provincia Autonoma di Bolzano”. This work was supported by the Open Access Publishing Fund of the Free University of Bozen-Bolzano.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. IPCC. Chapter 10 Emissions from Livestock and Manure Management. 2006. Available online: http://www.ipccggip.iges.or.jp/public/2006gl/pdf/4_Volume4/V4_10_Ch10_Livestock.pdf (accessed on 15 May 2019).
  2. Boussaada, A.; Arhab, R.; Calabro, S.; Grazioli, R.; Ferrara, M.G.; Musco, N.; Cutrignelli, M.I.; Tlidjane, M. Comparison of the Effect of Olives Leaves Extracts (Olea europaea) on In vitro Methane Production, Fermentation Efficiency and Protozoa Activity. Glob. Vet. 2017, 18, 445–453. [Google Scholar]
  3. FAO; GDP. Climate Change and the Global Dairy Cattle Sector—The Role of the Dairy Sector in a Low-Carbon Future; The Food and Agriculture Organization: Rome, Italy, 2018. [Google Scholar]
  4. Thomassen, M.A.; Van Calker, K.J.; Smits, M.C.J.; Iepema, G.L.; De Boer, I.J.M. Life cycle assessment of conventional and organic milk production in The Netherlands. Agric. Syst. 2008, 96, 95–107. [Google Scholar] [CrossRef]
  5. Braghieri, A.; Pacelli, C.; Bragaglio, A.; Sabia, E.; Napolitano, F. The Hidden Costs of Livestock Environmental Sustainability: The Case of Podolian Cattle. In The Sustainability of Agro-Food and Natural Resource Systems in the Mediterranean Basin; Vastola, A., Ed.; Springer: Berlin/Heidelberg, Germany, 2015; pp. 47–56. [Google Scholar] [CrossRef] [Green Version]
  6. Bragaglio, A.; Napolitano, F.; Pacelli, C.; Pirlo, G.; Sabia, E.; Serrapica, F.; Serrapica, M.; Braghieri, A. Environmental impacts of Italian beef production: A comparison between different systems. J. Clean. Prod. 2018, 172, 4033–4043. [Google Scholar] [CrossRef]
  7. Alan Rotz, C.; Asem-Hiablie, S.; Place, S.; Thoma, G. Environmental footprints of beef cattle production in the United States. Agric. Syst. 2019, 169, 1–13. [Google Scholar] [CrossRef]
  8. Sabia, E.; Napolitano, F.; Claps, S.; Braghieri, A.; Piazzolla, N.; Pacelli, C. Feeding. Nutrition and Sustainability in Dairy Enterprises: The Case of Mediterranean Buffaloes. In The Sustainability of Agro-Food and Natural Resource Systems in the Mediterranean Basin; Vastola, A., Ed.; Springer: Berlin/Heidelberg, Germany, 2015; pp. 57–64. [Google Scholar] [CrossRef] [Green Version]
  9. Sabia, E.; Napolitano, F.; Claps, S.; De Rosa, G.; Barile, V.L.; Braghieri, A.; Pacelli, C. Environmental impact of dairy buffalo heifers kept on pasture or in confinement. Agric. Syst. 2018, 159, 42–49. [Google Scholar] [CrossRef]
  10. Sabia, E.; Napolitano, F.; Claps, S.; De Rosa, G.; Braghieri, A.; Pacelli, C. Dairy buffalo life cycle assessment as affected by heifer rearing system. J. Clean. Prod. 2018, 192, 647–655. [Google Scholar] [CrossRef]
  11. Pardo, G.; Martin-Garcia, I.; Arco, A.; Yañez-Ruiz, D.R.; Moral, R.; del Prado, A. Greenhouse-gas mitigation potential of agro-industrial by-products in the diet of dairy goats in Spain: A life-cycle perspective. Anim. Prod. Sci. 2016, 56, 646–654. [Google Scholar] [CrossRef]
  12. Gutiérrez-Peña, R.; Mena, Y.; Batalla, I.; Mancilla-Leytón, J.M. Carbon footprint of dairy goat production systems: A comparison of three contrasting grazing levels in the Sierra de Grazalema Natural Park (Southern Spain). J. Environ. Manag. 2019, 232, 993–998. [Google Scholar] [CrossRef]
  13. Opio, C.; Gerber, P.; Mottet, A.; Falcucci, A.; Tempio, G.; MacLeod, M.; Vellinga, T.; Henderson, B.; Steinfeld, H. Greenhouse Gas Emissions from Rumin. Supply Chains—A Global Life Cycle Assessment; Food and Agriculture Organization of the United Nations (FAO): Rome, Italy, 2013. [Google Scholar]
  14. ISTAT, Agricultural Census at a Glance. Available online: http://censimentoagricoltura.istat.it/inbreve/?QueryId=&lang=en&graph=&subtheme=&cube (accessed on 25 June 2019).
  15. Battaglini, L.; Bovolenta, S.; Gusmeroli, F.; Salvador, S.; Sturaro, E. Environmental Sustainability of Alpine Livestock Farms. Ital. J. Anim. Sci. 2014, 13, 431–443. [Google Scholar] [CrossRef]
  16. Van den Pol-van Dasselaar, A.; Hennessy, D.; Isselstein, J. Grazing of Dairy Cows in Europe—An In-Depth Analysis Based on the Perception of Grassland Experts. Sustainability 2020, 12, 1098. [Google Scholar] [CrossRef] [Green Version]
  17. Cocca, G.; Sturaro, E.; Gallo, L.; Ramanzin, M. Is the abandonment of traditional livestock farming systems the main driver of mountain landscape change in Alpine areas? Land Use Policy 2012, 29, 878–886. [Google Scholar] [CrossRef]
  18. Giupponi, C.; Ramanzin, M.; Sturaro, E.; Fuser, S. Climate and land use changes, biodiversity and agri-environmental measures in the Belluno province, Italy. Environ. Sci. Policy 2006, 9, 163–173. [Google Scholar] [CrossRef]
  19. Marini, L.; Klimek, S.; Battisti, A. Mitigating the impacts of the decline of traditional farming on mountain landscapes and biodiversity: A case study in the European Alps. Environ. Sci. Policy 2011, 14, 258–267. [Google Scholar] [CrossRef]
  20. Scotti, A.; Füreder, L.; Marsoner, T.; Tappeiner, U.; Stawinoga, A.E.; Bottarin, R. Effects of land cover type on community structure and functional traits of alpine stream benthic macroinvertebrates. Freshw. Biol. 2019, 65, 524–539. [Google Scholar] [CrossRef] [Green Version]
  21. Bernués, A.; Ruiz, R.; Olaizola, A.; Villalba, D.; Casasùs, I. Sustainability of pasture-based livestock farming systems in the European Mediterranean context: Synergies and trade-offs. Livest. Sci. 2011, 139, 44–57. [Google Scholar] [CrossRef]
  22. Cillis, G.; Statuto, D.; Picuno, P. Vernacular Farm Buildings and Rural Landscape: A Geospatial Approach for Their Integrated Management. Sustainability 2020, 12, 4. [Google Scholar] [CrossRef] [Green Version]
  23. Kühl, S.; Flach, L.; Gauly, M. Economic assessment of small-scale mountain dairy farms in South Tyrol depending on feed intake and breed. Ital. J. Anim. Sci. 2020, 19, 41–50. [Google Scholar] [CrossRef]
  24. Batalla, I.; Knudsen, M.T.; Mogensen, L.; del Hierro, Ó.; Pinto, M.; Hermansen, J.E. Carbon footprint of milk from sheep farming systems in Northern Spain including soil carbon sequestration in grasslands. J. Clean. Prod. 2015, 104, 121–129. [Google Scholar] [CrossRef]
  25. Salvador, S.; Corazzin, M.; Romanzin, A.; Bovolenta, S. Greenhouse gas balance of mountain dairy farms as affected by grassland carbon sequestration. J. Env. Manag. 2017, 196, 644–650. [Google Scholar] [CrossRef]
  26. Stanley, P.L.; Rowntree, J.E.; Beede, D.K.; DeLonge, M.S.; Hamm, M.W. Impacts of soil carbon sequestration on life cycle greenhouse gas emissions in Midwestern USA beef finishing systems. Agric. Syst. 2018, 162, 249–258. [Google Scholar] [CrossRef]
  27. Janssens, I.A.; Freibauer, A.; Schlamadinger, B.; Ceulemans, R.; Ciais, P.; Dolman, A.J.; Heimann, M.; Nabuurs, G.-J.; Smith, P.; Valentini, R.; et al. The carbon budget of terrestrial ecosystems at country-scale a European case study. Biogeosciences 2005, 2, 15–26. [Google Scholar] [CrossRef] [Green Version]
  28. FAO. Principles for the Assessment of Livestock Impacts on Biodiversity. In Livestock Environmental Assessment and Performance (LEAP) Partnership; FAO: Rome, Italy, 2016. [Google Scholar]
  29. De Schryver, A.M.; Goedkoop, M.J.; Leuven, R.S.E.W.; Huijbregts, M.A.J. Uncertainties in the application of the species area relationship for characterization factors of land occupation in life cycle assessment. Int. J. LCA 2010, 15, 682–691. [Google Scholar] [CrossRef] [Green Version]
  30. Goedkoop, M.; Heijungs, R.; Huijbregts, M.; De Schryver, A.; Struijs, J.; van Zelm, R. ReCiPe A LifeCycle Impact Assessment Method which Comprises Harmonised Category indicators at the Midpoint and the Endpoint Level; Leiden University: Leiden, The Netherlands, 2008. [Google Scholar]
  31. FAO. Environmental Performance of Large Ruminant Supply Chains: Guidelines for Assessment. In Livestock Environmental Assessment and Performance Partnership; FAO: Rome, Italy, 2016. [Google Scholar]
  32. Guerci, M.; Knudsen, M.T.; Bava, L.; Zucali, M.; Scheonbach, P.; Kristensen, T. Parameters affecting the environmental impact of a range of dairy farming systems in Denmark, Germany and Italy. J. Clean. Prod. 2013, 54, 133–141. [Google Scholar] [CrossRef]
  33. Battini, F.; Agostini, A.; Tabaglio, V.; Amaducci, S. Environmental impacts of different dairy farming systems in the Po Valley. J. Clean. Prod. 2016, 112, 91–102. [Google Scholar] [CrossRef]
  34. Sabia, E.; Gauly, M.; Napolitano, F.; Serrapica, F.; Cifuni, G.F.; Claps, S. Dairy sheep carbon footprint and ReCiPe end-point study. Small Rumin. Res. 2020, 185, 106085. [Google Scholar] [CrossRef]
  35. SuperMix Farm Computer Systems, Cremona—Italy. Available online: http://www.fcs.it/supermix_inst.php (accessed on 10 May 2019).
  36. FAO, STAT. 2017. Available online: http://www.fao.org/faostat/en/#home (accessed on 25 November 2019).
  37. Guinée, J.B.; Gorree, M.; Heijungs, R.; Huppes, G.; Kleijn, R.; de Koning, A.; Van Oers, L.; Wegener Sleeswijk, A.; Suh, S.; Udo de Haes, H.A.; et al. Handbook on Life Cycle Assessment. In An Operational Guide to the ISO Standards; Kluwar Academic Publishers: Dordrecht, The Netherlands, 2002. [Google Scholar]
  38. IDF. A common carbon footprint approach for the dairy sector: The IDF guide to standard life cycle assessment methodology; Bulletin International Dairy Federation: Brussels, Belgium, 2010. [Google Scholar]
  39. ISPRA Agricoltura. Inventario nazionale delle emissioni e disaggregazione provinciale. In Rapporto Istituto Superiore per la Protezione e la Ricerca Ambientale (ISPRA); ISPRA: Rome, Italy, 2008. [Google Scholar]
  40. NRC. Nutrient Requirements of Dairy Cattle, 6th ed.; National Academy Press: Washington, DC, USA, 1989. [Google Scholar]
  41. IPCC. Chapter 11 N2O Emissions from Managed Soils, and CO2 Emissions from Lime and Urea Application. 2006. Available online: https://www.ipcc-nggip.iges.or.jp/public/2006gl/pdf/4_Volume4/V4_11_Ch11_N2O&CO2.pdf (accessed on 15 May 2019).
  42. Còndor, R.D. Agricoltura: Emissioni nazionali in atmosfera dal 1990 al Rapporto ISPRA 140/Istituto Superiore per la Protezione e la Ricerca Ambientale (ISPRA); ISPRA: Rome, Italy, 2011. [Google Scholar]
  43. Ecoinvent Centre. Ecoinvent Data v. 2.2 e Final Reports Ecoinvent 2000 No. 1-Swiss Centre for Life Cycle Inventories. 2007. Available online: https://www.ecoinvent.org/database/older-versions/ecoinvent-version-2/reports-on-ecoinvent-2/reports-on-ecoinvent-2.html (accessed on 15 May 2019).
  44. Soussana, J.F.; Tallet, T.; Blanfort, V. Mitigating the greenhouse gas balance of ruminant production systems though carbon sequestration in grasslands. Animal 2010, 4, 334–350. [Google Scholar] [CrossRef] [Green Version]
  45. Petersen, B.M.; Knudsen, M.T.; Hermansen, J.E.; Halberg, N. An approach to include soil carbon changes in life cycle assessments. J. Clean. Prod. 2013, 52, 217–224. [Google Scholar] [CrossRef]
  46. Tuomisto, H.L.; Hodge, I.D.; Riordan, P.; Macdonald, D.W. Comparing energy balances, greenhouse gas balances and biodiversity impacts of contrasting farming systems with alternative land uses. Agric. Syst. 2012, 108, 42–49. [Google Scholar] [CrossRef]
  47. Sturaro, E.; Marchiori, E.; Cocca, G.; Penasa, M.; Ramanzin, M.; Bittante, G. Dairy systems in mountainous areas: Farm animal biodiversity, milk production and destination, and land use. Livest. Sci. 2013, 158, 157–168. [Google Scholar] [CrossRef]
  48. Jiao, H.P.; Dale, A.J.; Carson, A.F.; Murray, S.; Gordon, A.W.; Ferris, C.P. Effect of concentrate feed level on methane emissions from grazing dairy cows. J. Dairy Sci. 2014, 97, 7043–7053. [Google Scholar] [CrossRef]
  49. Lawrence, D.C.; O’Donovan, M.; Boland, T.M.; Lewis, E.; Kennedy, E. The effect of concentrate feeding amount and feeding strategy on milk production, dry matter intake, and energy partitioning of autumn-calving Holstein-Friesian cows. J. Dairy Sci. 2015, 98, 338–348. [Google Scholar] [CrossRef] [PubMed]
  50. Schmitz, R.; Schnabel, K.; von Soosten, D.; Meyer, U.; Spiekers, H.; Rehage, J.; Danicke, S. The effects of energy concentration in roughage and allowance of concentrates on performance, health and energy efficiency of pluriparous dairy cows during early lactation. Arch. Anim. Nutr. 2018, 72, 100–120. [Google Scholar] [CrossRef] [PubMed]
  51. Kiefer, L.R.; Menzel, F.; Bahrs, E. Integration of ecosystem services into the carbon footprint of milk of South German dairy farms. J. Environ. Manag. 2015, 152, 11–18. [Google Scholar] [CrossRef] [PubMed]
  52. Penati, C.A.; Tamburini, A.; Bava, L.; Zucali, M.; Sandrucci, A. Environmental impact of cow milk production in the central Italian Alps using Life Cycle Assessment. Ital. J. Anim. Sci. 2013, 12, 584–592. [Google Scholar]
  53. Schader, C.; Jud, K.; Meir, M.S.; Kuhn, T.; Oehen, B.; Gattinger, A. Quantification of the effectiveness of greenhouse gas mitigation measures in Swiss organic milk production using a life cycle assessment approach. J. Clean. Prod. 2014, 73, 227–235. [Google Scholar] [CrossRef]
  54. Pirlo, G.; Lolli, S. Environmental impact of milk production from samples of organic and conventional farms in Lombardy (Italy). J. Clean. Prod. 2019, 211, 962–971. [Google Scholar] [CrossRef]
  55. Van Wyngaard, J.D.W.; Meeske, R.; Erasmus, L.J. Effect of concentrate feeding level on methane emissions, production performance and rumen fermentation of Jersey cows grazing ryegrass pasture during spring. Anim. Feed Sci. Technol. 2018, 241, 121–132. [Google Scholar] [CrossRef] [Green Version]
  56. Marino, R.; Atzori, A.S.; D’Andrea, M.; Iovane, G.; Trabalza-Marinucci, M.; Rinaldi, L. Climate change: Production performance, health issues, greenhouse gas emissions and mitigation strategies in sheep and goat farming. Small Rumin. Res. 2016, 135, 50–59. [Google Scholar] [CrossRef]
  57. Calabrò, S.; Infascielli, F.; Tudisco, R.; Musco, N.; Grossi, M.; Monastra, G.; Cutrignelli, M.I. Estimation of In vitro Methane Production in Buffalo and Cow. Buff. Bull. 2013, 32, 924–927. [Google Scholar]
  58. Ramírez-Restrepo, C.A.; Tan, C.; O’Neill, C.J.; Lòpez-Villalobos, N.; Padmanabha, J.; Wang, J.; McSweeney, C.S. Methane production, fermentation characteristics, andmicrobial profiles in the rumen of tropical cattle fed tea seedsaponin supplementation. Anim. Feed Sci. Technol. 2016, 216, 58–67. [Google Scholar] [CrossRef] [Green Version]
  59. Junior, P.F.; Cassiano, E.C.O.; Martins, M.F.; Romero, L.A.; Zapata, D.C.V.; Pinedo, L.A.; Marino, C.T.; Rodrigues, P.H.M. Effect of tannins-rich extract from Acacia mearnsii or monensin as feed additives on ruminal fermentation efficiency in cattle. Livest. Sci. 2017, 203, 21–29. [Google Scholar] [CrossRef]
  60. Ugbogu, E.A.; Elghandour, M.M.M.Y.; Ikpeazu, V.O.; Buendìa, G.R.; Molina, O.M.; Arunsi, U.O.; Emmanuel, O.; Salem, A.Z.M. The potential impacts of dietary plant natural products on the sustainable mitigation of methane emission from livestock farming. J. Clean. Prod. 2019, 213, 915–925. [Google Scholar] [CrossRef]
  61. Sabia, E.; Claps, S.; Napolitano, F.; Annicchiarico, G.; Bruno, A.; Francaviglia, R.; Sepe, L.; Aleandri, R. In vivo digestibility of two different forage species inoculated with arbuscular mycorrhiza in Mediterranean red goats. Small Rumin. Res. 2015, 123, 83–87. [Google Scholar] [CrossRef]
  62. Sabia, E.; Claps, S.; Morone, G.; Bruno, A.; Sepe, L.; Aleandri, R. Field inoculation of arbuscular mycorrhiza on maize (Zea mays L.) under low inputs: Preliminary study on quantitative and qualitative aspects. Ital. J. Agron. 2015, 10, 30–33. [Google Scholar] [CrossRef] [Green Version]
  63. Caputo, A.; Morone, G.; Di Napoli, M.A.; Rufrano, D.; Sabia, E.; Paladino, F.; Sepe, L.; Claps, S. Effect of destoned olive cake on the aromatic profile of cows’ milk and dairy products: Comparison of two techniques for the headspace aroma profile analysis. Ital. J. Agron. 2015, 10, 15–20. [Google Scholar] [CrossRef]
  64. Salvador, S.; Corazzin, M.; Piasentier, E.; Bovolenta, S. Environmental assessment of small-scale dairy farms with multifunctionality in mountain areas. J. Clean. Prod. 2016, 124, 94–102. [Google Scholar] [CrossRef]
  65. Pirlo, G.; Terzano, G.M.; Pacelli, C.; Abeni, F.; Carè, S. Carbon footprint of milk produced at Italian buffalo farms. Livest. Sci. 2014, 161, 176–184. [Google Scholar] [CrossRef]
  66. Hoogsteen, M.J.J.; Bakker, E.J.; van Eekeren, N.; Tittonell, P.A.; Groot, J.C.J.; van Ittersum, M.K.; Lantinga, E.A. Do Grazing Systems and Species Composition Affect Root Biomass and Soil Organic Matter Dynamics in Temperate Grassland Swards? Sustainability 2020, 12, 1260. [Google Scholar] [CrossRef] [Green Version]
  67. O’Brien, D.; Capper, J.L.; Garnsworthy, P.C.; Grainger, C.; Shalloo, L. A case study of the carbon footprint of milk from high-performing confinement and grass-based dairy farms. J. Dairy Sci. 2014, 97, 1835–1851. [Google Scholar] [CrossRef] [Green Version]
  68. Assandri, G.; Bogliani, G.; Pedrini, P.; Brambilla, M. Species-specific responses to habitat and livestock management call for carefully targeted conservation strategies for declining meadow birds. J. Nat. Conserv. 2019, 52, 125757. [Google Scholar] [CrossRef]
  69. Kaenchan, P.; Guinée, J.; Gheewala, S.H. Assessment of ecosystem productivity damage due to land use. Sci. Total Environ. 2018, 621, 1320–1329. [Google Scholar] [CrossRef] [PubMed]
  70. Riedener, E.; Rusterholz, H.P.; Baur, B. Effects of different irrigation systems on the biodiversity of species-rich hay meadows. Agric. Ecosyst. Environ. 2013, 164, 62–69. [Google Scholar] [CrossRef]
  71. Du, C.; Dias, L.C.; Freire, F. Robust multi-criteria weighting in comparative LCA and S-LCA: A case study of sugarcane production in Brazil. J. Clean. Prod. 2019, 218, 708–717. [Google Scholar] [CrossRef]
  72. Steinmann, Z.J.N.; Schipper, A.M.; Hauck, M.; Giljum, S.; Wernet, G.; Huijbregts, M.A.J. Resource Footprints are Good Proxies of Environmental Damage. Environ. Sci. Technol. 2017, 51, 6360–6366. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Aguilera, E.; Lassaletta, L.; Sanz-Cobena, A.; Garnier, J.; Vallejo, A. The potential of organic fertilizers and water management to reduce N2O emissions in Mediterranean climate cropping systems. A review. Agric. Ecosyst. Environ. 2013, 164, 32–52. [Google Scholar] [CrossRef] [Green Version]
  74. Chen, S.; Wu, D. Adapting ecological risk valuation for natural resource damage assessment in water pollution. Environ. Res. 2018, 164, 85–92. [Google Scholar] [CrossRef]
  75. Crenna, E.; Sinkko, T.; Sala, S. Biodiversity impacts due to food consumption in Europe. J. Clean. Prod. 2019, 227, 378–391. [Google Scholar] [CrossRef]
Figure 1. Pearson correlation coefficient (r) between two methods used for estimation of damage to biodiversity per kg of FPCM (FPCM: fat and protein corrected milk).
Figure 1. Pearson correlation coefficient (r) between two methods used for estimation of damage to biodiversity per kg of FPCM (FPCM: fat and protein corrected milk).
Sustainability 12 02128 g001
Table 1. Main characteristics of the dairy farms located in the South Tyrolean Alpine region.
Table 1. Main characteristics of the dairy farms located in the South Tyrolean Alpine region.
AGLC (n 11)AGHC (n 11)BSLC (n 11)BSHC (n 11)S.E.M.
Lactating cows, n°12.312.811.511.91.2
Cattle, n°5.54.25.64.51.1
Calves, n°5.2 a4.5 a2.6 b2.0 b0.9
Milk yield, kg FPCM/cow/y4206 Aa6031 Bc5029 b7356 Bd231
Concentrate, kg/cow/day3.0 A6.3 B3.7 A7.6 a0.3
Meadow hay, ha8.98.59.57.41.0
Permanent pasture, ha3.4 a2.52.20.7 b0.8
Pasture duration, day/y93 a517734 b16
Altitude of farms, m.s.l.114712861262120688
Electricity, kWh/y82607272562283171873
Fuel, liter/y97794811381335224
AGLC: Alpine Grey low concentrate; AGHC: Alpine Grey height concentrate; BSLC: Brown Swiss low concentrate; BSHC: Brown Swiss height concentrate; FPCM: fat and protein corrected milk; S.E.M.: standard error of mean. AB P ≤ 0.001; ab P ≤ 0.01; Ab P ≤ 0.001; cd P ≤ 0.001.
Table 2. Main ingredient and chemical composition (mean ± S.D.) of feed concentrates used in dairy farms in the South Tyrolean Alpine region.
Table 2. Main ingredient and chemical composition (mean ± S.D.) of feed concentrates used in dairy farms in the South Tyrolean Alpine region.
Ingredient (% of DM)AGLC (n = 11)AGHC (n = 11)BSLC (n = 11)BSHC (n = 11)
Maize flour 29 ± 329 ± 3 28 ± 129 ± 3
Wheat bran23 ± 926 ± 1320 ± 1126 ± 11
Soy bean meal 8 ± 311 ± 411 ± 513 ± 4
Rapeseed meal 9 ± 48 ± 48 ± 58 ± 5
Alfalfa, artificially dried7 ± 157 ± 155 ± 92 ± 6
Sugar beet pulp 7 ± 64 ± 77 ± 62 ± 6
Sunflower meal 7 ± 33 ± 45 ± 43 ± 4
By-products of distillation7 ± 33 ± 45 ± 43 ± 4
Sugar cane 2 ± 12 ± 1 2 ± 12 ± 1
Barley flour 1 ± 24 ± 71 ± 55 ± 8
Chemical composition (% of DM)
Crude protein 16.8 ± 1.916.5 ± 2.118.0 ± 1.217.9 ± 1.9
Crude fiber 8.2 ± 3.79.1 ± 3.57.6 ± 2.87.0 ± 1.8
Ether extract3.4 ± 0.43.4 ± 0.63.4 ± 0.33.8 ± 0.5
Crude ash5.8 ± 0.55.8 ± 0.56.0 ± 0.46.2 ± 0.6
Net Energy Lactation (MJ)7.0 ± 0.47.1 ± 0.47.0 ± 0.37.1 ± 0.2
AGLC: Alpine Grey low concentrate; AGHC: Alpine Grey high concentrate; BSLC: Brown Swiss low concentrate; BSHC: Brown Swiss high concentrate; DM: dry matter.
Table 3. Characterization factors of land use type [29].
Table 3. Characterization factors of land use type [29].
Land Use TypeMedian95% Confidence Level
Organic fertile grassland −0.01−0.18–0.15
Organic tall grassland0.04−0.12–0.18
Less intensive arable land0.440.31–0.54
Intensive arable land0.790.73–0.83
Baseline: semi-natural woodland0.00n.a.
Table 4. Carbon footprint (kg CO2-eq) per 1 kg of FPCM generated by different concentrated feed systems before and after including methods of carbon soil sequestration.
Table 4. Carbon footprint (kg CO2-eq) per 1 kg of FPCM generated by different concentrated feed systems before and after including methods of carbon soil sequestration.
AGLC (n = 11)AGHC (n = 11)BSLC (n = 11)BSHC (n = 11)S.E.M.
CF (kg CO2-eq/kg FPCM), without soil carbon sequestration1.55 Aa1.22 b1.43 a1.14 Bb0.06
CF (kg CO2-eq/kg FPCM), including soil carbon sequestration [44]1.33 Aa1.06 b1.23 a1.00 Bb0.06
Reduction of environmental impact (kg CO2-eq/kg FPCM) [44]0.23 Aa0.16 b0.21 a0.13 Bb0.01
% ON-Farms [44]79.769.381.757.69.1
CF (kg CO2-eq/kg FPCM), including soil carbon sequestration [45]1.22 Aa1.00 b1.17 a0.98 Bb0.05
Reduction of environmental impact (kg CO2-eq/kg FPCM) [45]0.34 Aa0.23 b0.26 a0.17 Bb0.02
% ON-Farms [45]84.086.183.960.98.2
AGLC: Alpine Grey low concentrate; AGHC: Alpine Grey high concentrate; BSLC: Brown Swiss low concentrate; BSHC: Brown Swiss high concentrate; FPCM: fat and protein corrected milk; S.E.M.: standard error of mean; CF: carbon footprint. AB P ≤ 0.001; ab P ≤ 0.01.
Table 5. Main processes responsible for the carbon footprint generated in the different farms (% on total greenhouse gas emissions).
Table 5. Main processes responsible for the carbon footprint generated in the different farms (% on total greenhouse gas emissions).
Process AGLC (n = 11)AGHC (n = 11)BSLC (n = 11)BSHC (n = 11)S.E.M.
Enteric fermentation 65.2 a56.6 b68.2 a56.2 b1.8
Concentrated feed 12.2 A22.1 Ba11.9 A27.5 Bb1.8
Meadow hay 8.8 Aa6.2 b8.0 4.9 B1.0
General consumptions 8.010.08.18.01.8
AGLC: Alpine Grey low concentrate; AGHC: Alpine Grey high concentrate; BSLC: Brown Swiss low concentrate; BSHC: Brown Swiss high concentrate; S.E.M.: standard error of mean. AB: P ≤ 0.001; ab: P ≤ 0.01.
Table 6. Main pollutants responsible for the carbon footprint generated in the different farms (% on total greenhouse gas emissions).
Table 6. Main pollutants responsible for the carbon footprint generated in the different farms (% on total greenhouse gas emissions).
PollutantAGLC (n = 11)AGHC (n = 11)BSLC (n = 11)BSHC (n = 11)S.E.M.
CH4 biogenic 75.0 Aa65.9 b70.2 ab57.8 B1.8
CO2 fossil11.5 a14.3 a17.2 b19.9 b1.8
N2O 11.6 Aa14.7 b10.5 Aa16.8 B1.0
AGLC: Alpine Grey low concentrate; AGHC: Alpine Grey high concentrate; BSLC: Brown Swiss low concentrate; BSHC: Brown Swiss high concentrate; S.E.M.: standard error of mean. AB: P ≤ 0.001; ab: P ≤ 0.01.
Table 7. Indicators of damage to biodiversity in different farm systems.
Table 7. Indicators of damage to biodiversity in different farm systems.
AGLC (n = 11)AGHC (n = 11)BSLC (n = 11)BSHC (n = 11)S.E.M.
Biodiversity, DS/kg FPCM [30]0.41 A0.46 A0.44 A0.56 B0.03
Damage to Ecosystem Diversity (Species*y)1/kg FPCM [31]1.78 E-07 a2.09 E-07 ABab1.60 E-07 A2.49 E-07 Bb1.75 E-08
AGLC: Alpine Grey low concentrate; AGHC: Alpine Grey high concentrate; BSLC: Brown Swiss low concentrate; BSHC: Brown Swiss high concentrate; S.E.M.: standard error of mean; DS: damage score; 1unit of species*y is the weighted damage on the basis of the total number of species on land and in water bodies [31]; AB: P ≤ 0.001; ab; P ≤ 0.01.

Share and Cite

MDPI and ACS Style

Sabia, E.; Kühl, S.; Flach, L.; Lambertz, C.; Gauly, M. Effect of Feed Concentrate Intake on the Environmental Impact of Dairy Cows in an Alpine Mountain Region Including Soil Carbon Sequestration and Effect on Biodiversity. Sustainability 2020, 12, 2128. https://0-doi-org.brum.beds.ac.uk/10.3390/su12052128

AMA Style

Sabia E, Kühl S, Flach L, Lambertz C, Gauly M. Effect of Feed Concentrate Intake on the Environmental Impact of Dairy Cows in an Alpine Mountain Region Including Soil Carbon Sequestration and Effect on Biodiversity. Sustainability. 2020; 12(5):2128. https://0-doi-org.brum.beds.ac.uk/10.3390/su12052128

Chicago/Turabian Style

Sabia, Emilio, Sarah Kühl, Laura Flach, Christian Lambertz, and Matthias Gauly. 2020. "Effect of Feed Concentrate Intake on the Environmental Impact of Dairy Cows in an Alpine Mountain Region Including Soil Carbon Sequestration and Effect on Biodiversity" Sustainability 12, no. 5: 2128. https://0-doi-org.brum.beds.ac.uk/10.3390/su12052128

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop