Next Article in Journal
Reprogramming Modernist Heritage: Enhancing Social Wellbeing by Value-Based Programming Approach in Architectural Design
Next Article in Special Issue
Comparison of Mechanical, Assisted and Manual Harvest of Origanum vulgare L.
Previous Article in Journal
Effect of Adding Emotion Recognition to Film Teaching—Impact of Emotion Feedback on Learning through Puzzle Films
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dairy Buffalo Life Cycle Assessment (LCA) Affected by a Management Choice: The Production of Wheat Crop

1
Council for Agricultural Research and Economics, Research Centre for Engineering and Agro-Food Processing, CREA-IT, 24047 Treviglio, Italy
2
Department of Medicine Veterinary, University of Bari “Aldo Moro”, 70010 Valenzano, Italy
3
Council for Agricultural Research and Economics, Research Centre Animal Production and Aquaculture, CREA-ZA, 26900 Lodi, Italy
*
Author to whom correspondence should be addressed.
Sustainability 2021, 13(19), 11108; https://0-doi-org.brum.beds.ac.uk/10.3390/su131911108
Submission received: 18 September 2021 / Revised: 1 October 2021 / Accepted: 4 October 2021 / Published: 8 October 2021

Abstract

:
Life cycle assessment (LCA) was performed in dairy buffalo farms representative of Southern Italian farming systems, similar due to several characteristics, with the exception of wheat production. This work evaluated the impacts derived from this management choice, comparing farms with wheat crop (WWC) or not (NWC). In agreement with the literature, economic allocation was chosen as a useful strategy to attribute equivalents to by-products, i.e., culled animals; the same criterion was also adopted to assign pollutants to wheat grain, limited to WWC farms. Environmental impacts in terms of Global Warming Potential (GWP, kg CO2 eq), Acidification Potential (AC, g SO2 eq), Eutrophication Potential (EU, g PO43-eq), Agricultural Land Occupation (ALO, m2y) and Water Depletion (WD, m3) were estimated. The production of wheat crop significantly affected (p < 0.05) the Agricultural Land Occupation (ALO) category as WWC farms need adequate land. WWC farms could allow a significant reduction in eutrophication (EU) compared to NWC farms (p < 0.05).

1. Introduction

Water buffalo (Bubalus bubalis) is a species that provides multiple products worldwide such as draught power, meat, milk, skin and manure. In Italy, buffalo farming has been conducted for centuries in extensive conditions based on marshland environments. Nowadays, most of the Mediterranean Italian Buffalo breed buffaloes are reared under intensive conditions to produce milk, almost entirely processed in mozzarella cheese. The production of buffalo mozzarella cheese is steadily increasing, +26.1% from 2013 to 2017, driven by the constant export increase, especially in Europe, the United States of America and Great Britain, and by the growing international interest in this product, as witnessed by the recent introduction of dairy buffalo in Germany and Great Britain [1,2,3].
Most of these animals are reared in the Protected Designation of Origin (PDO) area comprising the Campania, Lazio, Apulia and Molise regions, where recently buffalo husbandry has moved to more intensive farming conditions with a feeding system based on three different rations corresponding to the three main buffalo productive stages, i.e., lactating cows, dry cows and growing heifers [4]. Berlese et al. [1] also highlighted that in these conditions the animals have no access to pasture and water for wallowing.
Additionally, those areas are historically vocated to produce wheat, mainly to manufacture pasta, which is exported worldwide.
In Northern Italy, mozzarella cheese is sold directly to consumers through a short supply chain. In this context, dairy plants started transforming buffalo milk not exclusively in mozzarella cheese but also in other products, such as ripened cheeses and ricotta, which are appreciated by consumers [1].
Although there is growing attention on the environmental impact of livestock farming, to our knowledge, few studies focused on milk buffalo are available. In particular, few LCA studies have been developed on the allocation of a by-product such as wheat grain (Triticum durum Desf.).
As indicated by several authors [5,6,7], dairy systems often produce crop commodities (cereals), thus in cases of multifunctional processes the environmental impact should be shared among products. The aim of this study was to compare the environmental impact of buffalo milk, provided with two different management systems, integrated with economic allocation.

2. Materials and Methods

2.1. The Farms

All the farms of this study, located in the Apulia and Basilicata regions, are specialized in buffalo dairy farming, with animals kept in confinement. They adopt sexed semen mainly for the heifers to keep high-value animals, and their profiles are described in Table 1. The primary data were obtained from six farms distributed in two groups: No Wheat Crop (NWC) and With Wheat Crop (WWC). Despite the fact that the feeding of both groups is balanced with a 2:1 ratio = corn silage:wheat/straw, only the WWC group does not purchase supplementary forages, such as hay and straw. The WWC farms had different arable lands but a similar incidence of wheat income. The NWC1, WWC2 and WWC3 farms are located in the Protected Designation of Origin (PDO) area of Foggia Province (Apulia Region), whereas Basilicata (Potenza and Matera Provinces) is excluded from the PDO. According to widespread knowledge, Mediterranean Italian Buffalo breed animals show a noticeable heterogeneity in size, and this characteristic would affect management, feeding and productions, also conditioning dry matter intake of feed and milk yields.

2.2. Life Cycle Assessment (LCA) Methodology

The farms were analyzed with an LCA approach, based on ISO 14040 and ISO 14044 methodology [8,9]. The principles and framework for LCA include four distinct phases: (1) definition of the goal and scope (including functional unit and limits of the system); (2) life cycle inventory (LCI) analysis (including input and output data collection for all processes); (3) life cycle impact assessment (LCIA); (4) life cycle interpretation.

2.2.1. Goal and Scope Definition

In this study, the environmental impacts of milk obtained from two different farming managements in terms of Global Warming Potential (GWP, kg CO2 eq), Acidification Potential (AC, g SO2 eq), Eutrophication Potential (EU, g PO43-eq), Agricultural Land Occupation (ALO, m2y) and Water Depletion (WD, m3) were assessed. The functional unit was 1 kg of normalized buffalo milk, with a reference milk fat and protein (fat and protein corrected milk, FPCM) content of 8.3 and 4.73%, respectively. Raw milk was transformed into FPCM with the following equation [10,11]:
FPCM (kg/yr) = ({[(g of fat/L − 83) + (g of protein/L − 47.3)] × 0.00687} + 1) × milk production (kg/yr)

2.2.2. System Boundary Definition

The system boundaries considered in the study were comprised “from cradle to farm gate” (Figure 1 and Figure 2). All the on-farm operations (e.g., animal feeding and care, milking procedures) and sources pertaining to forage production (e.g., arable land, agrochemicals, water) were considered. The consumption of energy and the emissions from activities conducted off-farm were retrieved from databases provided by SimaPro 8.03 [12]. The transport of off-farm feeds, fossil fuels and bedding materials and their emissions were also included in the assessment.

2.2.3. Allocation Criterion

Dairy systems are typically multifunctional processes that coproduce meat from culled animals and surplus calves, for example, in addition to milk. In particular, the farms involved in this study are located in Southern Italian regions (Apulia and Basilicata), where the cultivation of wheat grain (Triticum durum Desf.) is largely widespread. This research focused on the presence of this peculiar coproduct, considered a relevant managing factor, to obtain two advantages for the three wheat-based farms:
  • A significant income in addition to milk.
  • The wheat also provides straw, adopted as fodder (mainly for dry cows) and litter. The With Wheat Crop (WWC) farms benefit from full forage self-sufficiency.
For these reasons, the economic allocation criterion was applied to both the by-product wheat grain and culled cows, in agreement with Pirlo et al. [11].
Young females are kept as replacers in the amount described by the owners, whereas male calves are sold at the 15th day of life; thus, the inputs and outputs provided by these animals were not considered.

2.2.4. Inventory Analysis and Input Data

The inventory data on livestock production, crop cultivations, straw for litter, inputs of purchased feed, electricity, diesel consumption and farm extension (divided in crop Ha year−1) were processed as primary data. The consumption and emission factor of natural gas adopted in the farms (water heating) was excluded due to the low impact and because it was not possible to have precise information about it. The feed consumption was retrieved by interviews with the farmers, the kind of water (i.e., well, tape) used in the farms was also obtained with interviews, whereas the consumption was assessed. In addition, a brief description about milk performances is reported in Table 2.
The inventory analysis involved the buildings, as sheds and concrete paved areas, whereas the structures and the bunkers (for corn silage) to stock the forages were excluded. A productive life of 50 years was assumed for these constructions, as suggested by the Ecoinvent 3 allocation database. The on-farm feed was assumed to be transported for 1 km with a tractor (Transport, tractor and trailer, agricultural (GLO), market for, Alloc, Def, S) whereas the bought feed, carried by truck (Transport, freight lorry > 32 metric ton, EURO 5, RER), was computed based on the distance from the farm, as suggested by Bragaglio et al. [13]. The sunflower meal feed, in most cases, is also imported [14]; thus, a 2700 km travel distance by truck from Ukraine was assessed. According to other authors [15], a 10,000 km journey by ship (Transport, freight, sea, transoceanic ship (GLO), market for, Alloc, Def, S) for soybean and cotton seeds from South America and a 12,000 km journey for palm oil from Malaysia were assumed.
Finally, the journey from the Italian harbor to farm gate was added to this computation. For fossil fuel provision, we also considered transport by truck: transport freight lorry of 3.5–7.5 metric ton, EURO 5, RER.
The water consumption was distinguished in specific categories, provided by the software Simapro 8.03: the water needed for transoceanic crops, the Italian well water mainly adopted for maize and the water adopted for the industrial processes (for example, wheat flour shorts). A description of the formulation and composition of the diets administered by the six farms is shown in Table 3 and Table 4, respectively, for cows (lactating and dry, LC and DC), heifers (HF) and young animals (91–365 days old, YA). The diets of the calves until the 90th day, based on milk replacers, weaning meal and a small amount of hay as sources of fiber, are very similar and were not reported in the table. Their amounts were loaded into the software and processed.

2.3. Emissions

Total emissions were estimated for the farms, according to the main difference that is wheat production, i.e., NWC and WWC. The fuel combustion, electricity consumption, enteric emissions, crop residue emissions, manure management and the emissions due to chemical fertilization were considered.

2.3.1. Enteric Emissions

The Refinement (2019) of the Intergovernmental Panel on Climate Change (IPCC) [16,17] methodology was adopted, considering the relationship between gross energy intake and emissions. According to the updated IPCC guidelines for the specific factors within the Tier 2 method, Equation (2) was adopted.
EF = [GE × (Ym/100) × x/55.56]
where:
EF = emission factor, kg, CH4 head−1 x−1;
GE = gross energy intake, MJ head−1, provided by different ingredients, expressed as dry matter;
Ym = methane conversion factor, percent of GE in feed converted to methane. Ym in dairy buffaloes is assumed to be 6.5, as indicated in Table 10.12 (updated);
x = the mean period, expressed in days, assumed for each livestock category, 90, 270 and 365 days, respectively, for dry cows, lactating cows and heifers. The enteric emissions of calves are not considered;
55.65 = the energy content of methane (MJ/kg CH4).
The data on the GE supplied by different feed sources (e.g., hay, straw and soybean meal) were provided by INRAE [18]. The different percentages of feed were also considered for the assessment of the GE.

2.3.2. Methane Emissions from Manure Management

As suggested by the Refinement (2019) of the IPCC methodology, the Tier 2 method relies on two primary types of inputs that affect the calculation of methane emission factors from manure: the manure characteristics and the animal waste management system characteristics (AWMS). Equation (3) was also adopted to assess these emissions:
EF = (VST*x) × [ B0(T) × 0.67 * ΣS,k*MCFS,k/100 × AWMS(T,S,k)]
where:
EF = CH4 emission factor for livestock category T, kg CH4 animal−1 days−1;
VST = daily volatile solid excreted for livestock category T, kg dry matter animal−1 days−1;
x = basis for calculating annual VS production, days year−1;
B0(T) = maximum methane producing capacity for manure produced by livestock category T, m3 CH4 kg−1 of VS excreted. IPCC guidelines suggest a default value of 0.10 for dairy buffaloes (Table 10.16 (updated));
0.67 = conversion factor of m3 CH4 to kilograms CH4;
MCF(S,k) = methane conversion factors for each manure management system S by climate region k, %. IPCC guidelines suggest a default value of 47% for animals kept in a paddock without distinguishing by climatic areas and relative humidity (Table 10.17 (updated)).
AWMS(T,S,k) = fraction of livestock category T’s manure handled using manure management system S in climate region k, dimensionless. IPCC guidelines suggest default values of 3, 34 and 63%, respectively, for paddock, slurry and solid management in Western Europe (Table 10A.6 (updated)).
The VS content of manure and urine was estimated, in accordance with the Tier 2 method, with Equation (4) (updated):
VS = [GE × (1 − DE%/100) + (UE × GE)] × (1 − ASH/18.45)
where:
VS = volatile solid excretion per day on a dry-organic matter basis, kg day−1;
GE = gross energy intake, MJ day−1;
DE% = digestibility rate of the feed. Different amounts of feed were also considered for the assessment of the DE, also considering the livestock categories. Different feeding periods were considered: from birth to the 90th day (weaning ration) and up to a year for the calves; 365, 95 and 270 days for the heifers, dry and lactating cows. Similarly for the GE, the data provided by INRAE [18] were adopted for the digestibility;
(UE × GE) = urinary energy expressed as fraction of GE. Typically, 0.04 GE can be considered urinary energy excretion by most ruminants, and this value was adopted in the current study;
ASH = the ash content of manure calculated as a fraction of the dry matter feed intake, specifying the different sources, as indicated in Table 3 and Table 4;
18.45 = conversion factor for dietary GE per kg of dry matter (MJ kg−1).

2.3.3. N2O Emissions from Manure Management

Although the IPCC guidelines indicate some criteria useful to assess nitrogen excretion, we followed the approach reported by Romano et al. [7] because some studies showed that buffaloes have a greater efficiency of N utilization compared with cattle [19,20]. Consequently, we adopted the equations suggested by Patra et al. [21] aimed to estimate the nitrogen excretion (urinary and fecal) in buffaloes. The N intake was counted by knowing the crude protein amount distinguished by livestock categories and farms.
The N2O emissions were then assessed adopting the IPCC 2019 guidelines [16,17].

Direct N2O Emissions

As suggested by IPCC 2019 [15], the Tier 1A method was applied, adopting Equation (5) (updated).
N2OD (mm) = [ΣST,P((NT,P × NexT,P) × AWMST,S,P) + Ncdg(s)] × EF3S] × 44/28
where:
N2OD (mm) = direct N2O emissions from manure management in the country, kg N2O year−1, assuming for each livestock category 95, 270, and 365 days, respectively, for dry and lactating cows, heifers and female calves;
N(T,P) = number of heads in each livestock category T on the farm, for production system P;
Nex(T,P) = annual average N excretion per head of livestock category T, on the farm, for production system P, previously assumed with the equation before;
AWMS(T,S,P) = fraction of total annual nitrogen excretion for each livestock category T that is managed in the manure management system S in the country, also considering the production system P (Tier 1A approach), dimensionless. These default values are provided by Table 10A.6 (updated);
Ncdg(s) = annual nitrogen input via co-digestate in the country, kg N yr−1, where the system (s) refers exclusively to anaerobic digestion. Not included in this study;
EF3(S) = emission factor for direct N2O emissions from the manure management system S in the country, kg N2O-N/kg N in the manure management system S; i.e., 0.01 and 0.002 for solid storage and paddock. These values are provided by Table 10.21 (updated) for solid storage [16] and by Table 11.1 (updated) for paddock [17]. In particular, for this manure management a value suggested for dry climates was adopted;
S = manure management system;
T = category of livestock;
44/28 = conversion of (N2O-N) mm emissions to N2O mm emissions.

Indirect N2O Emissions

Specific equations are useful to estimate nitrogen losses; subsequently, these data will be used to assess the indirect emissions due to volatilization and leaching. The Tier 1A method is applied to estimate both losses and the following equations (Equations (6) and (7) (updated)) are adopted to assess N due to leaching and volatilization, respectively:
NVOLATILIZATION-MMS = [ΣST,P((NT,P × NexT,P) × AWMST,S,P) + Ncdg(s) × FracgasMS(T,S)]]
NLEACHING-MMS = [ΣST,P((NT,P × NexT,P) × AWMST,S,P) + Ncdg(s) × FracleachMS(T,S)]]
where:
NVOLATILIZATION-MMS = amount of manure nitrogen that is lost due to the volatilization of NH3 and NOx, assuming for each livestock category 95, 270 and 365 days, respectively, for dry and lactating cows, heifers and female calves; kg N year−1;
NLEACHING-MMS = amount of manure nitrogen that is lost due to leaching, assuming for each livestock category 95, 270 and 365 days, respectively, for dry and lactating cows, heifers and female calves; kg N year−1;
N(T,P) = number of heads in each livestock category T on the farm, for production system P;
Nex(T,P) = annual average N excretion per head of livestock category T, on the farm, for production system P;
AWMS(T,S,P) = fraction of total annual nitrogen excretion for each livestock category T that is managed in the manure management system S in the country, also considering the production system P (Tier 1A approach), dimensionless;
Ncdg(s) = annual nitrogen input via co-digestate in the country, kg N yr−1, where the system (s) refers exclusively to anaerobic digestion. Not included in this study;
P = productivity class, high or low, of the system (Tier 1A approach);
FracgasMS(T,S) = fraction of managed manure nitrogen for livestock category T that volatilizes as NH3 and NOx in the manure management system S. The default value 0.12 is provided by Table 10.22 (updated) and is adopted for solid storage for “other animals” [16];
FracleachMS(T,S) = fraction of managed manure nitrogen for livestock category T that is leached from the manure management system S. The default value 0.02 is provided by Table 10.22 (updated) and is adopted for solid storage for “other animals” [16].
The emissions of N2O due to the volatilization and leaching of manure were assessed with the following equations (Equations (8) and (9)), adopting the previously calculated N losses:
N2OG (mm) = (NVOLATILIZATION-MMS × EF4) × 44/28
N2OL (mm) = (NLEACHING-MMS × EF5) × 44/28
where:
N2OG (mm) = indirect N2O emissions due to volatilization of N from manure management in the farm, assuming for each livestock category 95, 270 and 365 days, respectively, for dry and lactating cows, heifers and female calves, kg N2O year−1;
EF4 = 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. The default value is 0.005 (dry climate) kg N2O-N (kg NH3-N + NOx-N volatilized)−1 and is given in Chapter 11, Table 11.3 (updated) [17].
N2OL (mm) = indirect emissions due to leaching and runoff from manure management in the farm, assuming for each livestock category 95, 270 and 365 days, respectively, for dry and lactating cows, heifers and female calves, kg N2O year−1;
EF5 = emission factor for N2O emissions from nitrogen leaching and runoff, kg N2O-N/kg N leached and runoff. The default value of 0.011 kg N2O-N (kg N leaching/runoff)−1 is given in Chapter 11, Table 11.3 (updated) [16].

2.3.4. Ammonia

The ammonia emissions provided by the livestock were also calculated. In agreement with other studies [7,22], the emission factor 17/14 was adopted to estimate the NH3 amount pertinent to each farm.

2.3.5. CO2 Emissions from Livestock, Emissions from Crop, Soil Residues and Synthetic Fertilizers

The CO2 emissions from livestock were not estimated because the annual net CO2 emissions are assumed to be zero as the CO2 photosynthesized by plants is returned to the atmosphere as respired CO2 [16]. N2O and CO2 emissions from soils and crop residues are included in each crop input selected in the SimaPro database, although several categories were modified and accordingly loaded. The N2O emissions provided by urea and ammonium nitrate were also considered with Tier 2 and Tier 1 methods, respectively, for direct and indirect emissions, distinguished by the following equations:
N2ODIRECT-N = (FSN × EF1i)
where:
N2ODIRECT-N = direct N2O–N emissions from N inputs provided by urea and ammonium nitrate applied to the soils, kg N2O–N year−1;
FSN = annual amount of synthetic fertilizer N applied to soils, kg N year−1;
EF1i = emission factor for N2O emissions from N inputs, kg N2O–N (kg N input)−1. The default value for dry climates is available in Table 11.1 (updated) [17].
N2OATD-N = (FSN × FracGASF) × EF4
where:
N2OATD-N = indirect amount of N2O–N produced from atmospheric deposition of N volatilized from managed soils, kg N2O–N year−1;
FSN = annual amount of synthetic fertilizer N applied to soils, kg N year−1;
FracGASF = fraction of synthetic fertilizer N that volatilizes as NH3 and NOx, kg N volatilized (kg of N applied)−1. In Table 11.3 (updated), the default values for urea and ammonium nitrate-based emissions are available (0.15 and 0.05, respectively) [17];
EF4 = emission factor for N2O emissions from atmospheric deposition of N on soils and water surfaces, [kg N–N2O (kg NH3–N + NOx–N volatilized)−1], available in Table 11.3 (updated) for dry climate, i.e., 0.005 [16].
N2OLEACH-N = (FSN × Frac LEACH) × EF5
N2O LEACH-N = indirect amount of N2O–N produced from leaching and runoff of N additions to managed soils where leaching/runoff occurs, kg N2O–N year−1;
FSN = annual amount of synthetic fertilizer N applied to soils, kg N year−1;
Frac LEACH-N = fraction of all N added to managed soils in regions where leaching/runoff occurs that is lost through leaching and runoff, kg N (kg of N additions)−1. In Table 11.3 (updated), the default value is 0.24 [17];
EF5 = emission factor for N2O emissions from N leaching and runoff, kg N2O–N (kg N leached and runoff)−1. The default value is 0.11 in Table 11.3 (updated) [17].
The conversion of N2O–N emissions to N2O emissions was obtained by the following equation:
N2O = N2O–N × 44/28

2.3.6. Emissions from Electricity and Diesel Fuel

The CO2 equivalents, provided by the combustion of fossil fuels and from electricity use, were estimated considering the amount of diesel fuel and the kWh of electricity consumed for farm operations. The consumption and emission factor of natural gas adopted in the farms was excluded. The amount of purchased diesel was quantified through interviews, such as the electricity consumption. As suggested by ENAMA [23], a standard value of 0.85 kg per liter as diesel density and a 3.13 eq. emission factor to estimate CO2 release from the combustion of 1 kg of diesel were adopted. As for the electricity mix, we used the Italian emission factor (0.47 eq.) adopted in other studies [1,4,7].

2.4. Impact Assessment and Software

The software SimaPro 8.01 PhD, Pré Consultants 2015 was employed to estimate the environmental impacts. Two methods were adopted to assess the impact categories: (i) EPD 2013 for GWP, computed according to the CO2 equivalent factors in a 100-year time horizon, Acidification Potential (AP, g SO2 eq) and Eutrophication Potential (EP, g PO43-eq), and (ii) ReCiPe Midpoint (H) for Agricultural Land Occupation (ALO, m2y) and Water Depletion (WD, m3) [13,24,25]. Table 5 shows the impact categories and the characterization factors used to evaluate the potential environmental burden assessed with the EPD 2013 method.
As indicated by Gerssen-Gondelach et al. [27], Agricultural Land Occupation (ALO) is defined as the area of land needed to produce the FU chosen—in our research, 1 kg FPCM provided by buffalo. When the land occupation is not directly available with the inventory analysis of a study, these authors suggested the following equation:
m2/kg FPCM = total grassland or cropland requirement (Ha) × 10,000 (m2)/kg FPCM
SimaPro Pré Consultants, thanks to the ReCiPe method, allows us to investigate the Water Depletion (WD), processing the raw materials based on studies about water scarcity (WS), water stress index (WSI) and water productivity. The main algorithms adopted in the software, developed by Hoekstra [28,29], considered the environmental relevance of water productivity and the water footprint in water-rich areas.

3. Statistical Analysis

We calculated the GWP, AC, EU, ALO and WD of each farm using the LCA model. Then, these data were analyzed by one-way (general linear model procedure) ANOVA, using the production system as an independent variable. The data were analyzed using the “stats” package of R software [30], and then the Tukey test was adopted, developed with the LTukey function of the “laercio” R package [31]. In agreement with Silva and Azevedo [32], the Tukey test was chosen for carrying out the comparison test as it is characterized by greater rigor than other post hoc statistical comparison tests and has a greater control of type I error.

4. Results and Discussion

Buffalo dairy farms also produce meat and crop commodities, thus the environmental impact should be shared among the co-products [33]. According to other studies focused on bovine milk and meat [7,24,34], the allocation criterion was adopted. In particular, the economic allocation resulted in a more suitable strategy applicable for buffalo farms. In agreement with Pirlo et al. [11], in this study, the economic allocation criterion was applied, which was preferred to mass allocation because it better represents the societal cause of buffalo farm environmental impacts.

4.1. Role of the Allocation

The economic allocation mitigated the impacts, mainly in the WWC system. This criterion was applied to determine the weight of the co-products, avoiding the expansion of system boundaries. The inputs related to the crop (seeds, land, fuels and agrochemicals) were loaded into the software without lightening them when considering straw, which is always reused (in WWC1, WWC2 and WWC3) as feed or litter. In the three WWC farms, all the grain is sold as food for human consumption and all the straw is reused, and the overproduced straw is stored as farm stock. As previously indicated, the meat provided by male calves was excluded, whereas culled cows were involved in both systems. As reported by Mahath et al. [35], the following equation was adopted to subtract the equivalents from the functional unit (FU):
AFeco = (Veco × P)milk/Σ (Veco × P)milk, wheat grain, live weight culled cows
where:
AF is the allocation factor;
Veco is the economic value (EUR/kg);
P is the total production on-farm (kg/year), milk as FPCM, wheat as harvested grain without straw; live-weight cattle as culled cows only.
The economic values of the milk, wheat grain and culled cows were obtained from interviews with the owners and were applied to evaluate the economic allocation. Economic values, mass amounts and percentages are reported in Table 6.

4.2. LCA Categories. Results and Discussion

Table 7 shows the cradle-to-farm-gate life cycle results, related to the two different systems, i.e., NWC and WWC. The WWC group showed a total self-sufficiency in hay and straw supply, and in particular, the absence of purchased straw is related to the wheat crop. This arable area has significantly affected the Agricultural Land Occupation category as the production of wheat crop needs adequate land requirements. Additionally, WWC systems could also allow a significant reduction in eutrophication. These impacts exert a relevant effect on farm management, as shown by the PCA biplot (Figure 3), where the three WWC farms are grouped close together in the second quadrant. Although Berlese et al. [1] attributed to eutrophication a larger contribution of purchased feed than crops, in this study, the higher PO43- equivalents shown by the NWC system seem to find an explanation in the wider use of maize silage. Indeed, the comparison between NWC2/3 and WWC2/3 highlights the administration of corn silage in the NWC farms to all the livestock categories (Table 3 and Table 4).
The multivariate analysis of the main components has shown how the choice of the variables that describe the present study, and which are represented by the LCA descriptors (GWP, AC, EU, ALO, WD), allows us to obtain two main components that explain altogether almost 90% of the variability of the overall dataset. In particular, the first component, set as the horizontal axis of the biplot, explains 65.2% of the overall variability, while the second component, set as the vertical axis of the biplot, explains 24.3% of the overall variability.
The biplot obtained from the PCA analysis (Figure 3) allows us to observe that the two variables that clearly discriminate the observed farms are the ALO and the GWP, placing the WWC farms on the right of the floor, well separated from the NWC farms on the left of the floor. The first farms are in fact characterized by higher ALO and GWP values than the second group.
Regarding the second component, we observed an influential effect of the EU variable between two NWC farms for higher values than WWC farms, creating a vertical separation. Therefore, in the upper part of the biplot there are farms with lower EU values, while in the lower part of the biplot there are farms with higher EU values. In the upper left part, there is a farm belonging to the NWC group, which differs for higher values of WD and AC, but it must be remembered that the statistical analysis of variance did not show such statistically significant differences. The analysis then places the position of the centroids, defined as points whose coordinates are the average of the coordinates of the group elements obtained from the diametrically opposite components, top right for WWC farms and bottom left for NWC farms.
GWP: The impact on climate change was assessed using the IPCC characterization factors for the 100-year time horizon, according to the EDIP (2013) method, adopted in the study. As for climate stability, permanent grasslands store nearly as much carbon as forests (EIP-AGRI 2014) [36], and the carbon sequestration potential of permanent pastures was estimated between 0.01 and 0.3 Gt (gross tons of volume) C year−1 [37,38], but all the farms involved in our research were without pasture. On other hand, the crop, field and fertilizer emissions were assessed with equations provided by IPCC guidelines [17].
Usually, the impacts are influenced by the weight of the functional unit (FU); for example, Berlese et al. [1] found a GWP of 6.4 kg CO2 eq per kg of FPCM, higher than values reported by Pirlo et al. [11] (5.1 kg CO2 eq per kg of FPCM) and also higher than our values (4.91 and 5.19 kg CO2 eq for NWC and WWC, respectively). Indeed, the milk production observed by Berlese et al. [1] was 1409 kg FPCM/buffalo per lactation, in comparison to the mean value showed by Pirlo et al. [11] of 2251 kg FPCM/buffalo per lactation—a comparable amount with ours (2866 and 2000 kg FPCM/buffalo per lactation, respectively, for NWC and WWC systems). As suggested by Gerber et al. [39], the milk production of a farm is a crucial parameter because, in general, high productivity reduces environmental impact. Although NWC farms showed a higher milk yield than WWC farms, the CO2 equivalents were not affected.
AC: In this study, no significant differences were found between the two systems. The values obtained in this study were very similar to those reported by Berlese et al. [1] (37.3 ± 3.97 g SO2 eq), and lower than those reported by Pirlo et al. [11] (65 g SO2), where all the farms were corn silage-based.
A study focused on dairy buffaloes [4] compared a pasture-based system with a confined rearing system, showing that the main contribution for SO2 equivalents was derived from the production of maize silage in the confined system. In agreement with these authors, our research showed that, despite the absence of differences between NWC and WWC farms for acidification potential, the higher trend for mean, minimum and maximum values in NWC farms could be explained by significant inputs related to the maize crop.
In the no-pasture-based system, the ammonium was identified by the authors as the first polluting substance followed by SO2 and NOX. Even though the research investigated bovine milk, Guerci et al. [40] found for grazing cattle (pasture only) very low SO2 equivalents—7.44 g. This result may be attributed to the low use of inorganic fertilizers in pasture-based systems.
In addition, although the emissions of ammonia from manure management are recognized as significant contributors to acidification [4,41], the farms investigated in this study showed comparable manure handling, and the highest values observed in the NWC1 farm suggested a clarification. Research focused on dairy cattle [42] highlighted that the contributions to acidification were spilt among feed productions (both on- and off-farm), with a share of 31–57%. Moreover, Bragaglio et al. [12] found that the largest source of SO2 equivalents is the production of concentrates. The purchasing of off-farm feed (total amount of concentrates, alfalfa hay, straw and more than 50% of hay meadow) concerning the NWC1 farm would explain the high AC value. In Figure 4 and Figure 5, an explanatory comparison of feed sources between NWC1 and WWC1, the two farms without corn silage, is also given.
EU: This LCA descriptor, related to acidification, is affected by the conduction system and by environmental conditions, such as temperature, relative humidity or rainfall. Table 7 suggests that NWC farms showed significantly (p < 0.05) higher values (14.43 g PO43-eq) than WWC farms (11.70 g PO43-eq), which showed a more homogenous profile as a system (Figure 3). Sabia et al. [41] identified a functional unit (FU) of 1 kg of weight gain in the period needed to reach the age of puberty in buffalo heifers and compared animals kept in the pasture or in confinement. They reported higher values of eutrophication potential (g PO43-) in the confined system. Similar findings, although focused on marine eutrophication (g N equivalents), were also reported by other authors [4]. The authors explained the different impacts because the main source of pollution was the production of corn silage, as also observed in dairy cattle by Bartl et al. [43].
As indicated in the third PCA biplot quadrant (Figure 3), the NWC2 and NWC3 farms, producing corn silage, are characterized by high PO43- equivalents. In comparison with NWC corn silage-based farms, WWC farms showed a different administration of this forage. In NWC2 and NWC3 farms, all the livestock categories are feed with corn silage, whereas in the WWC2 farm this is provided to lactating cows, heifers and young animals. Finally, in the WWC3 farm its administration is reserved only for lactating cows.
ALO: Table 7 shows that lower values were found in the NWC system (p < 0.05) compared to the WWC system. Many studies estimated the land use (LU), land use change (LUC) or the agricultural land occupation (ALO) impacts of organic vs. conventional or pasture-based vs. confined farming modes. Often, impact categories show higher values in organic or pasture-based systems, as the low yields (milk and forages) and management efficiency partially explain these results [7,43,44]. In this research, mean yields in NWC and WWC farms of 2866 and 2000 kg FPCM/buffalo per lactation, respectively, were recorded. In addition, WWC farms showed the lowest milk production, as shown in Table 2. The arable land of these farms is characterized by a high extent required for wheat crop (38, 66 and 35% for WWC1, WWC2 and WWC3, respectively); moreover, the fodder self-supply is also ensured by a significant hay area. Finally, the low number of lactating cows in the three WWC farms (Table 1), in addition to explaining the low productions, would suggest a less specialized management compared to the NWC mode. This aspect needs to be confirmed with studies involving a larger number of farms.
The adoption of economic allocation, useful to subtract the pollutants attributable to the co-produced wheat grain, has not allocated enough m2 year−1 equivalents, such as to identify the WWC system neither equally nor less impacting than the NWC one.
WD: This LCA category should be considered from a global perspective, since freshwater is a global resource, with growing global freshwater demand while global freshwater availability is limited [29]. In agreement with this framework, several studies are focused on the water footprint, distinguishing green water, blue water and grey water.
In our research, no significant differences were found comparing the NWC and the WWC system (1.79 ± 0.41 vs. 1.69 ± 0.14 m3). Sabia et al. [4] stated that no data are available on WD in dairy buffaloes, while conflicting results are available for dairy cattle. In Noord-Brabant (Netherlands), De Boer et al. [45] observed that 66 l of consumptive water was needed for 1 kg of FPCM, whereas at world level a water footprint of 1207 m3/ton of milk [46] and a consumption of 544 l per 1 kg di FPCM [45] have been reported. In the other impact categories, dairy buffaloes [11] showed impacts 4-fold higher than dairy cattle [40] due to their lower milk production (roughly 4-fold lower even in terms of FPCM), even though input levels and categories and animal and farm dimensions were roughly similar.
In agreement with these studies [11,41], other authors found [47,48,49,50], in a similar Mediterranean environment, comparable results (0.38–0.55 m3 and 0.52 m3 per kg bovine FPCM, respectively), with a 4-fold ratio.
In the PCA biplot (Figure 3), the NWC1 farm showed, as indicated in Table 7, a different profile from the other NWC farms, with lower m3 equivalents. Despite the absence of irrigated crops (corn), WD could be explained by information acquired from the interviews: drinking water (human use) was indicated as an exclusive source, also adopted to clean the milking parlor. Although the NWC1 and WWC1 farms are characterized by a total absence of crop sources of grain (e.g., oat, barley, maize), additionally the NWC1 farm suffers from a lack of self-supply of forages (all the alfalfa and almost 50% of meadow hay were purchased).

5. Conclusions

This work performed an LCA of buffalo dairy farms, selected for similar general characteristics but different farming managements. A proportion of 1:2 of no maize silage and with maize silage was chosen as a system profile for both NWC and WWC farms. The main goal was the environmental impact assessment of dairy buffalo farming, according to forage self-sufficiency and combination with wheat crop production. WWC farms were characterized by a high availability of wheat straw, adopted as litter and administrated as fodder. In addition, a large arable area ensured a significant amount of hay; therefore, it was necessary to only purchase concentrates. The excess of straw and hay was managed with storage and, considering their low economic value, the allocation of the surplus would probably not have subtracted pollutants from the functional unit, i.e., FPCM. The productions of wheat crop entail adequate land requirements and the low milk yield observed in the WWC farms probably exacerbated the high values of Agricultural Land Occupation in agreement with other studies.
Interestingly, a high land occupation and low production, widespread in Mediterranean pasture-based systems, are often characterized by poor nutritional values of grazing grass.
The ALO results suggest that full self-sufficiency of forages would not be enough to describe a virtuous profile of the farms. Probably higher hay yields would have allowed the farms to obtain lower values of land occupation.
In order to estimate the possible reduction in the emissions due to the modification of some input variables, the performed analysis considers the possibility of reducing the fertilizers used in the agricultural stage. In particular, precision agriculture would allow targeted fertilization, leading to a reduction in the amount of agrochemicals for similar yields; a possible additional step could be to distinguish the different stages (tillage, fertilization, harvesting) with respect to the overall impact (midpoint categories). Afterwards, the base case could be compared with the improved scenario according to the normalized ReCiPe endpoint categories.
The PO43 equivalents (EU) seemed to be mainly conditioned by the feeding strategies. Although in the life cycle assessment the impacts are often mitigated by high FU yields, this trend has not been highlighted in eutrophication.
As reported, dairy systems are often multifunctional processes; however, the allocation of by-products does not always mitigate the impacts. In particular, this study showed that a high availability of forages may have mitigated eutrophication (WWC farms), but it may have more heavily affected another disputed impact category: Agricultural Land Occupation.

Author Contributions

Conceptualization, P.D.P. and A.B.; methodology, A.B. and F.T.; software, A.B., E.R. and A.M.; validation, P.D.P. and E.R.; formal analysis, F.T. and E.R.; investigation, A.B.; resources, P.D.P. and A.M.; data curation, E.R.; writing—original draft preparation, E.R., P.D.P., F.T., A.M. and A.B.; writing—review and editing, F.T., E.R. and A.B.; visualization, P.D.P.; supervision, A.M. and F.T.; project administration, P.D.P. and A.B.; funding acquisition, A.B. All authors have read and agreed to the published version of the manuscript.

Funding

Andrea Bragaglio’s research activity is granted by the European Union and Italian Ministry of Education, University and Research in the program PON 2014–2020 Research and Innovation, framework Attraction and International Mobility-1839894, Activity 1. The project was approved by the Italian Ministry of Education, University and Research.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

a)
CLAL, 2019. Production of Mozzarella di Bufala Campana PDO. Italy. https://www.clal.it/index.php/section=mozzarella_bufala_campana (accessed on 5 August 2021).
b)
FAOSTAT, 2020. www.fao.org/faostat/en/#data (accessed on 5 August 2021).
c)
INRAE, www.inrae.fr. https://mediatheque.inrae.fr/ (accessed on 5 August 2021).
d)
ENAMA, 2005. Handbook of fuel consumption for agricultural subsidized use (Prontuario dei consumi di carburante per l’impiego agevolato in agricoltura). Available from: http://www.enama.it/php/pageflip.php_pdf_enama_int_prontuario.pdf&dir/it/pdf/monografie (accessed on 5 August 2021).
e)
R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Wien, Österreich, 2019; Available online: https://www.R-project.org/ (accessed on 15 August 2021).
f)
Laercio Junio da Silva. laercio: Duncan Test, Tukey Test and Scott-Knott Test. R Package Version 1.0-1. 2010. Available online: https://CRAN.R-project.org/package=laercio (accessed on 15 August 2021).

Acknowledgments

The authors are grateful to Giovanna Calzaretti and Francesco Giannico for their technical support. The authors are grateful to Laura Maresca, veterinary surgeon.

Conflicts of Interest

All authors declare that they do not have any conflicts of interest that could inappropriately influence this manuscript.

References

  1. Berlese, M.; Corazzin, M.; Bovolenta, S. Environmental sustainability assessment of buffalo mozzarella cheese production chain: A scenario analysis. J. Clean. Prod. 2019, 238, 117922. [Google Scholar] [CrossRef]
  2. Borghese, A. Buffalo livestock and products in Europe. Buffalo Bull. 2013, 32, 50–74. [Google Scholar]
  3. CLAL, 2019. Production of Mozzarella di Bufala Campana PDO. Italy. Available online: https://www.clal.it/index.php/section=mozzarella_bufala_campana (accessed on 5 August 2021).
  4. 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]
  5. Pirlo, G. Cradle-to-farm gate analysis of milk carbon footprint: A descriptive review. Ital. J. Anim. Sci. 2012, 11, e20. [Google Scholar] [CrossRef]
  6. Mu, W.; Van Middelaar, C.E.; Bloemhof, J.M.; Engel, B.; De Boer, I.J.M. Benchmarking the environmental performance of specialized milk production systems: Selection of a set of indicators. Ecol. Indic. 2017, 72, 91–98. [Google Scholar] [CrossRef]
  7. Romano, E.; Roma, R.; Tidona, F.; Giraffa, G.; Bragaglio, A. Dairy farms and Life Cycle Assessment (LCA). The allocation criterion useful to estimate undesirable products. Sustainability 2021, 13, 4354. [Google Scholar] [CrossRef]
  8. ISO 14040, Environmental Management—Life Cycle Assessment—Principles and Framework; International Organization for Standardization ISO Central Secretariat Chemin de Blandonnet, 8CP 401—1214 Vernier; ISO: Geneva, Switzerland, 2006.
  9. ISO 14044, Environmental Management—Life Cycle Assessment—Requirements and Guidelines; International Organization for Standardization ISO Central Secretariat Chemin de Blandonnet, 8CP 401—1214 Vernier; ISO: Geneva, Switzerland, 2006.
  10. Di Palo, R. Produzione di Latte Nella Bufala con Diete Tradizionali e con l’impiego di Acidi Grassi. Ph.D. Thesis, University of Naples, Naples, Italy, 1992. [Google Scholar]
  11. Pirlo, G.; Carè, S.; Fantin, V.; Buttol, P.; Terzano, G.M.; Masoni, P.; Pacelli, C. Factors affecting life cycle assessment of milk produced on 6 Mediterranean buffalo farms. J. Dairy Sci. 2014, 97, 6583–6593. [Google Scholar] [CrossRef]
  12. PRe Consultants, 2013. Introduction to LCA with SimaPro. Available online: www.pre-sustainability.com/download/SimaPro8IntroductionToLCA.pdf (accessed on 16 August 2021).
  13. 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]
  14. FAOSTAT. 2020. Available online: www.faofaostat/en/#data (accessed on 5 August 2021).
  15. Nguyen, T.L.T.; Hermansen, J.E.; Mogensen, L. Environmental consequences of different beef production systems in the EU. J. Clean. Prod. 2010, 18, 756–766. [Google Scholar] [CrossRef]
  16. IPCC. Chapter 10: Emissions from livestock and manure management. In Agriculture, Forestry and Other Land Use; Refinement to the 2006 IPCC Guidelines for National Greenhouse Gas Inventories; IPCC: Geneva, Switzerland, 2019; Volume 4, pp. 1–87. Available online: https://www.ipcc-nggip.iges.or.jp/public/2019rf/pdf/4_Volume4/19R_V4_Ch10_Livestock.pdf (accessed on 5 August 2021).
  17. IPCC. Chapter 11: N2O emissions from managed soils, and CO2 emissions from lime and urea application. In Agriculture, Forestry and Other Land Use; Refinement to the 2006 IPCC Guidelines for National Greenhouse Gas Inventories; IPCC: Geneva, Switzerland, 2019; Volume 4, pp. 1–54. Available online: https://www.ipcc_nggip.iges.or.jp/public/2019rf/pdf/4_Volume4/19R_V4_Ch11_Soils_N2O_CO2.pdf (accessed on 5 August 2021).
  18. INRAE. Available online: www.inrae.fr. https://mediatheque.inrae.fr/ (accessed on 5 August 2021).
  19. Singh, K.; Kundu, S.S.; Kushwaha, B.P.; Maity, S.B. Dietary energy levels response on nutrient utilization, nitrogen balance and growth in Bhadawari buffalo calves. Livest. Res. Rural Dev. 2009, 21, 8. [Google Scholar]
  20. Neglia, G.; Balestrieri, A.; Gasparrini, B.; Cutrignelli, M.I.; Bifulco, G.; Salzano, A.; Cimmino, R.; Varricchio, E.; D’Occhio, M.J.; Campanile, G. Nitrogen and Phosphorus Utilisation and Excretion in Dairy Buffalo Intensive Breeding. Ital. J. Anim. Sci. 2014, 13, 3362. [Google Scholar] [CrossRef] [Green Version]
  21. Patra, A.K.; Pal, K.; Lalhriatpuii, M. Prediction of nitrogen excretion in buffalo production systems using dietary and animal variables. Agric. Syst. 2020, 182, 102845. [Google Scholar] [CrossRef]
  22. Addison, P.F.E.; De Bie, K.; Rumpff, L. Setting conservation management thresholds using a novel participatory modeling approach. Conserv. Biol. 2015, 29, 1411–1422. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. ENAMA. Handbook of Fuel Consumption for Agricultural Subsidized Use (Prontuario dei Consumi di Carburante per L’impiego Agevolato in Agricoltura). 2005. Available online: http://www.enama.it/php/pageflip.php_pdf_enama_int_prontuario.pdf& dir/it/pdf/monografie (accessed on 5 August 2021).
  24. Bragaglio, A.; Braghieri, A.; Pacelli, C.; Napolitano, F. Environmental impacts of beef as corrected for the provision of ecosystem services. Sustainability 2020, 12, 3828. [Google Scholar] [CrossRef]
  25. Huijbregts, M. Life-cycle Impact Assessment of Acidifying and Eutrophying Air Pollutants. Calculation of Equivalency Factors with RAINS-LCA; Interfaculty Department of Environmental Science, Faculty of Environmental Science, University of Amsterdam: Amsterdam, The Netherlands, 1999. [Google Scholar]
  26. Heijungs, R.; Guinee, J.B.; Huppes, G.; Lankreijer, R.M.; Udo de Haes, H.A.; Wegener-Sleeswijk, A.; Ansems, A.M.M.; Eggels, P.G.; van Duin, R.; de Goede, H.P. Environmental Life Cycle Assessment of Products: Guide and Backgrounds; Centre of Environmental Science (CML), Leiden University: Leiden, The Netherlands, 1992; Available online: http://hdl.handle.net/1887/8061 (accessed on 16 August 2019).
  27. Gerssen-Gondelach, S.J.; Lauwerijssen, R.B.G.; Havlík, P.; Herrero, M.; Valin, H.; Faaij, A.P.C.; Wicke, B. Intensification pathways for beef and dairy cattle production systems: Impacts on GHG emissions, land occupation and land use change. Agric. Ecosyst. Environ. 2017, 240, 135–147. [Google Scholar] [CrossRef]
  28. Hoekstra, A.Y. A critique on the water-scarcity weighted water footprint in LCA. Ecol. Indic. 2016, 66, 564–573. [Google Scholar] [CrossRef] [Green Version]
  29. Hoekstra, A.Y. Water Footprint Assessment: Evolvement of a New Research Field. Water Resour. Manag. 2017, 31, 3061–3081. [Google Scholar] [CrossRef] [Green Version]
  30. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Wien, Österreich, 2019; Available online: https://www.R-project.org/ (accessed on 15 August 2021).
  31. Da Silva, J.L.; Duncan Test, Tukey Test and Scott-Knott Test. R Package Version 1.0-1. 2010. Available online: https://CRAN.R-project.org/package=laercio (accessed on 15 August 2021).
  32. De Assis Santos e Silva, F.; De Azevedo, C.A.V. Comparison of means of agricultural experimentation data through different tests using the software Assistat. Afr. J. Agric. Res. 2016, 11, 3527–3531. [Google Scholar] [CrossRef]
  33. Cederberg, C.; Stadig, M. System expansion and allocation in life cycle assessment of milk and beef production. Int. J. Life Cycle Assess. 2003, 8, 350–356. [Google Scholar] [CrossRef]
  34. 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]
  35. Mahath, C.S.; Mophin Kani, K.; Dubey, B. Gate-to-gate environmental impacts of dairy processing products in Thiruvananthapuram. India Resour. Conserv. Recycl. 2019, 141, 40–53. [Google Scholar] [CrossRef]
  36. EIP-AGRI Focus Group. Profitability of Permanent Grassland. How to Manage Permanent Grassland in a Way That Combines Profitability, Carbon Sequestration and Biodiversity. 2014. Available online: http://ec. europa.eu/eip/agriculture/sites/agrieip/files/fg9_permanent_grassland_profitability_starting_ paper_2014_en.pdf (accessed on 16 August 2019).
  37. Lal, R. Soil carbon sequestration impacts on global climate change and food security. Science 2004, 304, 1623–1627. [Google Scholar] [CrossRef] [Green Version]
  38. Braghieri, A.; Pacelli, C.; Bragaglio, A.; Sabia, E.; Pacelli, C. 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: Cham, Switzerland, 2015. [Google Scholar] [CrossRef] [Green Version]
  39. Gerber, P.; Vellinga, T.; Opio, C.; Steinfeld, H. Productivity gains and greenhouse gas emissions intensity in dairy systems. Livest. Sci. 2011, 139, 100–108. [Google Scholar] [CrossRef]
  40. Guerci, M.; Knudsen Trydeman, M.; Bava, L.; Zucali, M.; Schönbach, P.; Kristensen, T. Parameters affecting the environmental impact of a range of dairy farming systems in Denmark, Germany and Italy. J. Clean. Prod. 2013, 5, 133–141. [Google Scholar] [CrossRef]
  41. 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]
  42. Baldini, C.; Bava, L.; Zucali, M.; Guarino, M. Milk production Life Cycle Assessment: A comparison between estimated and measured emission inventory for manure handling. Sci. Total Environ. 2018, 625, 209–219. [Google Scholar] [CrossRef] [PubMed]
  43. Bartl, K.; Gomez, C.A.; Nemecek, T. Life cycle assessment of milk produced in two smallholder dairy systems in the highlands and the coast of Peru. J. Clean. Prod. 2011, 13, 1494–1505. [Google Scholar] [CrossRef]
  44. 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]
  45. De Boer, I.J.M.; Hoving, I.E.; Vellinga, T.V.; Van de Wen, G.W.J.; Leffelaar, P.A.; Gerber, P.J. Assessing environmental impacts associated with freshwater consumption along the life cycle of animal products: The case of Dutch milk production in Noord-Brabant. Int. J. Life Cycle Assess. 2013, 18, 193–203. [Google Scholar] [CrossRef] [Green Version]
  46. Mekonnen, M.M.; Hoekstra, A.Y. A Global Assessment of the Water Footprint of Farm Animal Products. Ecosystems 2012, 15, 401–415. [Google Scholar] [CrossRef] [Green Version]
  47. Mekonnen, M.M.; Hoekstra, A.Y. The green, blue and grey water footprint of crops and derived crop products. Hydrol. Earth Syst. Sci. 2011, 15, 1577–1600. [Google Scholar] [CrossRef] [Green Version]
  48. Noya, I.; González-García, S.; Berzosa, J.; Baucells, F.; Feijoo, G.; Moreira, M.T. Environmental and water sustainability of milk production in Northeast Spain. Sci. Total Environ. 2018, 616–617, 1317–1329. [Google Scholar] [CrossRef] [PubMed]
  49. Cosentino, C.; Adduci, F.; Musto, M.; Paolino, R.; Freschi, P.; Pecora, G.; D’Adamo, C.; Valentini, V. Low vs high “water footprint assessment” diet in milk production: A comparison between triticale and corn silage based diets. Emir. J. Food Agric. 2015, 27, 312–317. [Google Scholar] [CrossRef] [Green Version]
  50. Freschi, P.; Musto, M.; Paolino, R.; Cosentino, C. Grazing and biodiversity conservation: Highlights on a natura 2000 network site. In The Sustainability of Agro-Food and Natural Resource Systems in the Mediterranean Basin; Vastola, A., Ed.; Springer International Publishing AG: Basel, Switzerland, 2015; pp. 271–288. [Google Scholar] [CrossRef] [Green Version]
Figure 1. System boundaries of No Wheat Crop (NWC) system.
Figure 1. System boundaries of No Wheat Crop (NWC) system.
Sustainability 13 11108 g001
Figure 2. System boundaries of With Wheat Crop (WWC) system.
Figure 2. System boundaries of With Wheat Crop (WWC) system.
Sustainability 13 11108 g002
Figure 3. Biplot of the PCA analysis with the centroids of the two groups. Centroids are shown by bold boxes as WWC and NCW, respectively.
Figure 3. Biplot of the PCA analysis with the centroids of the two groups. Centroids are shown by bold boxes as WWC and NCW, respectively.
Sustainability 13 11108 g003
Figure 4. The graphs show the different percentage compositions of diets (lactating cows, LC; heifers, HF; dry cows, DC) in the No Wheat Crop 1 (NWC1) farm.
Figure 4. The graphs show the different percentage compositions of diets (lactating cows, LC; heifers, HF; dry cows, DC) in the No Wheat Crop 1 (NWC1) farm.
Sustainability 13 11108 g004
Figure 5. The graphs show the different percentage compositions of diets (lactating cows, LC; heifers, HF; dry cows, DC) in the With Wheat Crop 1 (WWC1) farm.
Figure 5. The graphs show the different percentage compositions of diets (lactating cows, LC; heifers, HF; dry cows, DC) in the With Wheat Crop 1 (WWC1) farm.
Sustainability 13 11108 g005
Table 1. Profile of the farms and of the two groups. Urea 46%N and ammonium nitrate 27%N.
Table 1. Profile of the farms and of the two groups. Urea 46%N and ammonium nitrate 27%N.
No Wheat Crop (NWC)With Wheat Crop (WWC)
NWC1NWC2NWC3WWC1WWC2WWC3
Geographical place (Province)FoggiaPotenzaPotenzaMateraFoggiaFoggia
Total crop area, Ha20658065225270
Hay, Ha2050404050140
Barley, Ha--10---
Maize silage, Ha-1530-1515
Maize grain, Ha----1020
Wheat, Ha---2515095
Herd, heads n.197303446203479613
Lactating cows, n.6012015052160185
Dry cows, n.25508060180185
Heifers, n.908018550100175
Young < 365 days, n. 204525353060
Bulls, n.286698
Urea, t y−1-10.020.037.545.059.5
Ammonium nitrate, t y−1--7.0---
Phosphate, t y−1----0.81.6
Potassium chloride, t y−1--5.2---
Concrete area (shed, services), m21500400035001000600011,000
Milking parlor size, m2200200400150300300
Milk tank, liters140060005000250040006000
Diesel, liters y−118,80023,50076,50021,20064,70070,500
Electricity, kWh y−145,70065,00097,00050,00076,00087,600
Table 2. Dairy performances and milk components of No Wheat Corn (NWC) and With Wheat Corn (WWC) farms.
Table 2. Dairy performances and milk components of No Wheat Corn (NWC) and With Wheat Corn (WWC) farms.
ItemUnitNWC1NWC2NWC3WWC1WWC2WWC3
FPCM yieldkg/year254,405370,400481,800146,000438,000474,850
Average FPCM per lactation *kg/head313730872374207620251897
Fat mean ± SD%7.38 ± 0.317.46 ± 0.457.83 ± 0.407.54 ± 0.418.30 ± 0.418.57 ± 0.53
Protein mean ± SD%4.63 ± 0.124.25 ± 0.164.48 ± 0.074.42 ± 0.214.65 ± 0.094.68 ± 0.14
Bacterial count mean ± SDcfu/mL76,670 ± 50,796316,224 ± 25,215131,667 ± 113,90091,167 ± 37,64973,077 ± 67,66197,947 ± 47,339
* Assuming 270 and 95 days for lactation and dry period, respectively [4].
Table 3. No Wheat Crop (NWC) diets. Animal categories: lactating cows (LC), dry cows (DC), heifers (HF), young animals (YA); # Vicia faba minor; * kg/head/day.
Table 3. No Wheat Crop (NWC) diets. Animal categories: lactating cows (LC), dry cows (DC), heifers (HF), young animals (YA); # Vicia faba minor; * kg/head/day.
NWC1NWC2NWC3
CategoryLCDCHFYALCDCHFYALCDCHFYA
Forage kg/head/day
Meadow hay8.05.53.52.010.07.05.03.03.2-3.21.5
Alfalfa hay2.0-1.8---------
Straw-5.5---4.02.01.0-7.02.0-
Maize silage----8.05.03.02.019.06.013.02.0
Raw concentrate kg/head/day
Maize flour/grain4.0-1.3-----4.0--0.5
Barley1.5-0.41.0----2.02.0-1.0
Soybean meal1.1-----------
Wheat flour shorts-2.0----------
Market concentrate g/head/day
Soybean seeds (roasted)310-----------
Soybean seeds (dehulled/flaked)240-4804002600-8004001650260200150
Sunflower meal480-4403602000-600260450380300220
Cotton seeds540---1600-500220450---
Maize flour450---300-10080-130100-
Fava bean #480-----------
Wheat flour shorts280-480400-----70056080
Beet pulp280----------350
Linseeds240-----------
Wheat flour220-----------
Bran--440360----800680540400
Maize germ meal--220180-----130100-
Maize distillers--6050----80013010080
Palm oil60-----------
Molasses50-6050-----13010080
Chemical composition (%)
Dry matter (DM) *17.811.78.24.418.111.48.35.018.212.110.54.6
Crude protein (%DM)12.459.6515.5015.0014.009.7012.3514.2014.259.0511.8014.50
Ether extract (%DM)5.803.503.453.703.152.752.903.404.403.253.154.15
Crude fiber (%DM)22.7036.6028.0018.7027.6035.9026.6532.2021.4029.8525.2018.95
Ash (%DM)7.607.008.358.357.253.905.454.107.207.057.456.75
Table 4. With Wheat Crop (WWC) diets. Animal categories: lactating cows (LC), dry cows (DC), heifers (HF), young animals (YA); § Avena sativa L. with Vicia sativa L.; # Vicia faba minor; ‡ Pisum sativum L.; * kg/head/day.
Table 4. With Wheat Crop (WWC) diets. Animal categories: lactating cows (LC), dry cows (DC), heifers (HF), young animals (YA); § Avena sativa L. with Vicia sativa L.; # Vicia faba minor; ‡ Pisum sativum L.; * kg/head/day.
WWC1WWC2WWC3
CategoryLCDCHFYALCDCHFYALCDCHFYA
Feed kg/head/day
Oat hay----3.0-1.52.03.58.010.03.5
Mixed meadow hay §10.0-5.02.8--------
Straw-7.5--1.58.52.01.02.03.0--
Maize silage----15.0-3.01.515.0---
Raw concentrate kg/head/day
Maize grain/flour----6.52.04.01.04.0---
Barley------------
Soybean meal----2.90.82.00.51.0---
Fava bean #--------2.0---
Pea ‡--------1.0---
Market concentrate g/head/day
Maize flour32007501100660-------400
Soybean seeds (dehulled/flaked)1800400400240-------100
Soybean seeds (roasted)--1100660-------80
Barley1400-----------
Molasses4005012070-------50
Bran200900360220-------270
Palm oil100-----------
Sunflower meal-300---------250
Beet pulp-100900550--------
Wheat flour shorts-----------380
Chemical composition (%)
Dry matter (DM) *16.09.57.94.617.510.39.64.517.19.88.94.5
Crude protein (%DM)13.508.2012.1512.4014.108.3014.5512.6013.309.9012.0013.80
Ether extract (%DM)3.602.153.203.305.303.405.755.554.302.753.003.45
Crude fiber (%DM)22.2532.6022.2021.8022.3033.5021.6023.0025.3536.7035.0029.65
Ash (%DM)4.855.804.254.255.204.254.905.005.153.303.004.10
Table 5. Characterization factors of the main elementary flows of the impact categories investigated with EPD 2013.
Table 5. Characterization factors of the main elementary flows of the impact categories investigated with EPD 2013.
Category Main Elementary FlowCharacterization FactorSource
Global Warming Potential, kg CO2 eqCO21
CH428IPCC 2019 [15]
N2O265
Acidification, kg SO2 eqNH31.6
NOx0.76Huijbregts 1999 [25]
SO21.2
Eutrophication, kg PO43 eqNO30.1
P2O33.06Heijungs et al. 1992 [26]
Table 6. Yields and economic incomes (%) derived from milk, wheat grain and culled cows (No Wheat Crop, NWC; With Wheat Crop, WWC).
Table 6. Yields and economic incomes (%) derived from milk, wheat grain and culled cows (No Wheat Crop, NWC; With Wheat Crop, WWC).
NWC1NWC2NWC3NWC MeanWWC1WWC2WWC3WWC Mean
Milk income EUR/kg1.601.601.551.581.501.501.601.53
Milk yield/year kg254,405370,400481,800368,868146,000438,000470,850351,617
Total milk income EUR 407,050592,640746790582,160219,000657,000752930542,977
Wheat income EUR/kg----0.330.300.270.30
Wheat yield/year kg----75,000450,000332,000285,670
Total wheat income EUR ----24,750135,00089,77583,175
Culled cows income EUR/head300300300300300300300300
Culled cows head/year82025185161612
Total culled cows income EUR 24006000750053001500480048003700
Milk income %99.4199.0099.0199.1489.3082.4588.8586.86
Wheat income % 0.000.000.000.0010.0916.9410.5912.54
Culled cows income %0.591.000.990.860.610.600.570.59
Table 7. Descriptive characteristics and Tukey’s test.
Table 7. Descriptive characteristics and Tukey’s test.
SYSTEMFARMSGWP
kg CO2 eq
AC
g SO2 eq
EU
g PO43-eq
ALO
m2y−1
WD
m3
No Wheat CropNWC14.6550.2313.278.522.26
NWC25.1439.9614.4210.641.55
NWC34.9537.8115.6211.561.56
mean NWC4.91a42.66a14.43a10.24b1.79a
standard deviation0.246.631.171.560.41
With Wheat CropWWC15.3038.0711.7313.701.84
WWC25.5130.7511.2115.691.57
WWC34.7539.2412.1715.041.65
mean WWC5.19a36.02a11.70b14.81a1.69a
standard deviation0.394.600.481.010.14
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Romano, E.; De Palo, P.; Tidona, F.; Maggiolino, A.; Bragaglio, A. Dairy Buffalo Life Cycle Assessment (LCA) Affected by a Management Choice: The Production of Wheat Crop. Sustainability 2021, 13, 11108. https://0-doi-org.brum.beds.ac.uk/10.3390/su131911108

AMA Style

Romano E, De Palo P, Tidona F, Maggiolino A, Bragaglio A. Dairy Buffalo Life Cycle Assessment (LCA) Affected by a Management Choice: The Production of Wheat Crop. Sustainability. 2021; 13(19):11108. https://0-doi-org.brum.beds.ac.uk/10.3390/su131911108

Chicago/Turabian Style

Romano, Elio, Pasquale De Palo, Flavio Tidona, Aristide Maggiolino, and Andrea Bragaglio. 2021. "Dairy Buffalo Life Cycle Assessment (LCA) Affected by a Management Choice: The Production of Wheat Crop" Sustainability 13, no. 19: 11108. https://0-doi-org.brum.beds.ac.uk/10.3390/su131911108

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