Paula Toro-Mujica- Instituto de Ciencias Agroalimentarias, Animales y Ambientales (ICA3), Universidad de O'Higgins, San Fernando, Chile
One of the main production challenges associated with climate change is the reduction of carbon emissions. Increasing the efficiency of resource utilization is one way to achieve this purpose. The modification of production systems through improved reproductive, genetic, feed, and grazing management practices has been proposed to increase technical–economic efficiency, even though the “environmental viability” of these modifications has not always been evaluated. The objective of this study was to evaluate the use of feeding and management strategies on the carbon footprint (CF) and economic variables in the traditional cow–calf system in southern Chile using a simulation model. The modifications evaluated corresponded to combinations of stocking rate, use of creep feeding practices with different supplementation levels, and the incorporation of feed additives to the supplement, using factorial experiments. Additionally, the scenarios were evaluated with and without carbon sequestration. The CF for the baseline scenarios was 12.5 ± 0.3 kg of CO2−eq/kg of live weight (LW) when carbon sequestration was considered and 13.0 ± 0.4 kg of CO2−eq/kg of LW in the opposite case. Changes in stocking rate, supplementation level, and consideration of carbon sequestration in pasture and soil had a significant effect on the CF in all simulated scenarios. The inclusion of additives in the supplement did not have a significant effect on production costs. With regard to reducing greenhouse gas (GHG) emissions, incorporating canola oil presented the best average results. The model developed made the selection of environmentally viable feed strategies or management adaptations possible.
Introduction
Livestock is one of the main agricultural activities responsible for the anthropogenic emissions of greenhouse gas (GHG). These GHG emissions from global livestock supply chains reach ~7 gigatonnes annually. Of this, 65% are associated with cattle production (1). Modifying livestock systems is a strategy to increase technical and economic efficiency. At the same time, interest exists in decreasing GHG. However, obtaining precise measurements or estimations is the first barrier to such a challenge. The Intergovernmental Panel on Climate Change (IPCC) has generated empirical methodologies with different degrees of precision (Tiers) for estimating GHG for livestock production systems through mathematical equations (2). On the other hand, more recent researchers, such as studies carried out by Cottle and Eckard (3), Niu et al. (4), and Patra and Lalhriatpuii (5), proposed models for beef cattle, dairy cattle, and small ruminants, respectively. Cottle and Eckard (3), through a meta-analysis applied to the MitiGate database of methane-mitigation studies, determined an equation that predicts the methane yield in beef cattle. The equation included the following variables: measurement method, breed, diet, and country. Niu et al. (4) used the fitting linear mixed models to estimate methane emissions, including variables, such as dry matter (DM) intake, dietary composition, milk production and composition, and body weight. Patra and Lalhriatpuii (5) developed linear and non-linear statistical models for small ruminants to predict enteric methane emissions from goats, using the variables dietary nutrient composition, intake of nutrients and energy, digestibility of energy, and organic matter.
Undoubtedly, diet and various associated nutritional management strategies are fundamental variables in any enteric methane emission estimation model. At the same time, it is one of the main routes to increase the efficiency of the livestock production systems (6–10). Then, feeding systems, pastoral resource management (11, 12), and a comprehensive perspective (13) are research areas that have been studied to reduce GHG emissions. Within nutritional management, various research studies have addressed the inclusion of additives to the diet. The proposed additives include nitrates, lipids (e.g., canola and soybean oil), monensin, yeasts, plant extracts (e.g., oregano, green tea, tannins, and saponins tea), seaweed, propolis extract, and mixtures of them. In most of these (14–20), in grazing or confined animals, a decrease in enteric methane emissions has been observed. Evaluating the effects of these changes in the general management and some nutritional management of livestock farms has been estimated through simulation models. Lurette et al. (21) created a model to evaluate the sensitivity of pastoral dairy systems to scenarios of seasonal biomass production variability. White et al. (22) used a whole-farm model to study the relative importance of reproduction, genetics, and nutritional management on minimizing the environmental impact of meat production systems. Furthermore, Rotz et al. (23) developed a simulation model to quantify environmental footprints for beef production systems in the United States. However, the use of simulation models to estimate the effect of the inclusion of feed additives in livestock production cycles, both at the level of total emissions and at the level of economic variables', is scarce (24). Then, as Legesse et al. (25) pointed out, a need exists to include evaluation of the contribution of various feeding systems to GHG emissions that can be covered through simulation models.
In southern Chile, where more than 63% of the cattle production is located, a steady reduction in the number of farms and cattle heads has occurred in recent years (26). In this area, a predominance of cow–calf systems exists, generally within small farms with low technical and economic efficiency (9). The animals are usually kept under extensive grazing regimens. This is supplemented during the winter with hay or silage. The products of these systems correspond to calves of 6–8 months old that are sold for finishing on other farms. In search of increasing the environmental viability of the system, the objective of this study was to evaluate the use of feeding and management strategies on the carbon footprint (CF) and economic variables in the traditional cow–calf system in southern Chile using a simulation model.
Materials and Methods
Study Area and Production System
The southern zone of Chile is constituted by three regions: La Araucanía, Los Ríos, and Los Lagos. The Region of La Araucanía has the second-largest cattle population, reaching 0.7 million heads, representing 17% of the country's cattle (24). This zone of Chile is between 38° and 39°30′S latitude and covers an area of 31,842.3 km2. The cattle farms are small to medium size, with an average herd of 17.7 heads of cattle (standard deviation: 81.2 cattle heads), and where 97% of farmers own < 100 animals (24). Feeding is based on natural and sown grasslands, and some crops cover nearly 10% of the total area (9, 27).
The region's climate is temperate rainy, with an average temperature of 11.2°C, an average rainfall of ~1,200 mm, and 57% of the precipitation between May and August. Forage production varies between 2 and 8 tons/ha, depending on the type of pasture, level of fertilization, and management of grazing. Forage production is seasonal, concentrating between 50 and 60% of production in the spring. This distribution forces farmers to conserve forage as hay or silage to supplement the animal diet during winter and summer (9). In addition, the seasonality of forage production determines the seasonality of the bovine production. The calving season is concentrated between the end of winter and the beginning of spring (28). This ensures the best fit between the supply and demand of forage. Thus, the breeding season is carried out in November and December (29). In the traditional cow–calf production system, the production cycle lasts between 16 and 20 months, depending upon the length of the breeding season and the weaning age. The conventional product of the current breeding systems corresponds to steer calves and heifer calves of 6–8 months with weights of 180–240 kg. Typical stocking rates in the zone studied are around 0.5–1.4 livestock units/ha/year. These stocking rates demonstrate that the systems tend to be extensive. The cow–calf system uses natural or sown grassland, and some farmers have supplementary crops. The most common supplementary crops are corn silage, oats, and alfalfa. The supplementation of cows during the winter months (June to August) is a common practice, as well as the supplementation during the breeding months.
Limits of the Systems
The cradle-to-farm-gate methodology (30) utilized by Toro-Mujica et al. (31), Sykes et al. (32), and Cerri et al. (33) was used to estimate the system CF. This methodology involves assessing the emissions associated with all the inputs and processes within a production cycle. This evaluation is carried out until the animal leaves the farm. Therefore, emissions related to transport, slaughterhouse, livestock fairs, or other farms were not included (34–36). According to Rotz et al. (36), these activities represent <1% of the total energy used. Emissions associated with veterinary and reproductive inputs were not considered (37). Emissions related to inputs required for the production of forage and supplements, such as seed, fertilizers, pesticides, and fuel expenditures on machinery, were assigned to the DM production of the respective crop (Figure 1). Soil carbon sequestration, an ecosystem service, was estimated for soils and pastures (38–41). Quantification of the CF was calculated per kg live weight (LW) of cattle sold.
Figure 1. System boundaries, greenhouse gas (GHG) sources, and storage (sinks) of the cow–calf system.
Model Description
A simulation model was developed in Excel 2016, using spreadsheets and macros of VBA based on the model proposed by Catrileo et al. (9) to evaluate management modifications in cow–calf systems on temperate pastures in Chile. The model contains an initial data sheet where the user defines the characteristics of the farm, including size and general management. With regard to general management, fertility percentage, replacement rate, and cows and calf mortality were also included. For feed management, the model has a configuration sheet where the quantity and type of supplement are defined month by month. The model allows monthly selection and set supplementation for cows and calves. In the case of calves, supplementation is simulated through the creep feeding system. From the model, it is assumed that the farm has three pasture paddocks. In each, the type of pasture present must be selected. Subsequently, the user must select month by month which paddock will be used. The model has a database where information about nutrient composition, the CF of supplements and pasture, climatic data, and price of products, among others, is stored. The model makes daily estimates for forage intake and energy requirements, based on the equations proposed by CSIRO (42, 43) and AFRC (44). A random component was incorporated in the variables voluntary forage intake, milk production, and cow and calf mortality to obtain a stochastic model. The normal reverse distribution was used to incorporate randomness into the variables for voluntary forage intake and milk production, with a standard deviation of 7 and 7.5%, respectively (9).
A variation of Bernoulli's test was used to incorporate randomness into the variable mortality. Then, each animal was considered as an independent test (45), generating for each of them a random number with uniform distribution (0–100). When this number was less than or equal to the mortality variable, it meant that the animal has died.
Pastures and Supplements
The pastures used to simulate the model were tall fescue with subterranean clover (28). Data for monthly growth rates (kg DM/ha), digestibility of DM (%), and crude protein (%) were available. The inputs required for sowing and fertilization were measured in terms of CF based on Edwards-Jones et al. (46) and Saunders and Barber (47). For animals' supplementation with corn silage, energy concentration of 2.5 Mcal/kg DM, 7.5% crude protein, and 72% digestibility of DM (48) were used. The CF of the corn silage was assigned a value of 0.2 kg CO2−eq/kg, as reported by Adom et al. (49). Emissions attributable to additives and their transportation were not included since the low inclusion percentage does not significantly modify the CF of the supplement.
Carbon Sequestration
As pointed out by Toro-Mujica et al. (31), carbon sequestration was assumed to be 10% of the carbon deposited in the soil by fertilization (39) over a time period of 100 years (38). In the case of organic matter deposited in the soil by senescence of the pasture, it was assumed that 15% of the organic matter was converted into the organic matter found in soils (50).
Estimate of the Carbon Footprint
For CF, the original estimation algorithm of the model was complemented with the Tier 2 equations proposed by the IPCC (2, 51) to assess GHG emissions. These equations relate gross energy intake to the production of methane. Carbon dioxide (CO2), methane (CH4), and nitrous oxide (N2O) emissions were transformed into CO2 equivalents (CO2-eq) based on the conversion factors (1 kg of CO2-eq equal to 1 kg of CO2, 25 kg of CO2−eq equal to 1 kg of CH4, and 298 kg of CO2-eq equal to 1 kg of N2O) as proposed by the IPCC (2, 46).
The emissions converted into CO2−eq were added and divided by total LW (kilogram) sold out of the farm gate (culled cows and weaned claves) to calculate CF assuming kg LW at the farm gate as the functional unit (52). The selection of this functional unit was based on the possibility of comparing the results with those of other studies.
Baseline Simulated Scenarios
The baseline scenario corresponded to the traditional cow–calf production system (Group I) identified by Toro-Mujica et al. (53). The definition of a cow–calf production system within the simulation model includes the month of the production cycle initiation, stocking rate, age, weight of cows, pasture areas, type and time of use of the pastures, cow fertility, duration of lactation, and cow and calf mortality, as well as data referring to the use of inputs for grassland fertilization. In the actual system, the herd consisted of cows of different ages and weights; but in the model, average weight and age were used. In the baseline scenario, a stocking rate of 0.7 livestock units (LU/ha), a production cycle of 20 months, the supplementation of cows for 90 days in winter (June to August), and 90 days of postpartum were chosen. The supplementation for cows was carried out with corn silage (2 kg DM per cow per day). The calves were not supplemented.
Proposed Feed Strategies and Experimental Scenarios
The simulation model flexibility allows using combinations of many variables. However, in order to simplify the interpretation of the results, only three variables were evaluated in the experiment. The variables included stocking rate, use of creep feeding with five levels of supplementation (0, 0.25, 0.5, 0.75, and 1 kg DM/calf/day) with corn silage from 2 months of age to age of sale, and use of feed additives in the corn silage for feeding calves (monensin with a dose of 30 mg/kg DM and canola oil with a dose of 46 g/kg DM) (17, 54).
The effect of incorporating feed additives was carried out based on the results reported by Beauchemin and McGinn (17) and Ranga Niroshan Appuhamy et al. (54). The authors suggested that the inclusion of monensin reduced methane emissions from 9 to 15% per unit of gross energy intake. On the other hand, the reduction due to the addition of canola oil was 21% per unit of gross energy intake. With regard to the cost of adding monensin to corn silage, an increase of 3% was estimated over the initial price due to its cost and its mixture with the supplement. In the case of canola oil, the price increase of the supplement was 13%.
The dimensional input variables, management variables, and the levels used in the variables' stocking rate, feed additives, and supplementation to simulate are shown in Table 1. Table 1 also includes the inputs (urea, seed, diesel, and pesticides) associated with the maintenance of the improved pasture.
Table 1. Input data.
Data were analyzed using R software (55). All the experimental scenarios were analyzed with and without consideration of carbon sequestration in pasture and soil in order to evaluate if the variation of CF depends on the incorporation of this factor. The results were evaluated using multifactorial ANOVA (56). Analysis of multiple comparisons of means was conducted through Tukey's honestly significant difference (HSD) test or the Kruskal–Wallis test, after proof equality of variances with Levene's test. The analyses were performed with two data sets:
Baseline stocking rate (0.7 LU/ha) scenario, modified through the calf supplementation and additive uses.
Baseline stocking rate (0.7 LU/ha) scenario, modified through the change of stocking rate, the calf supplementation, and additive uses.
Since the model is stochastic, 20 simulations (reps) were run for each variable's combinations. The number of replicates required for the stochastic model was tested by the methodology described by Toro-Mujica et al. (20).
Results
Baseline Scenario
The CF in the baseline scenario was 12.5 ± 0.3 kg CO2−eq/kg of LW when carbon sequestration was not considered. The incorporation of carbon sequestration increased the CF to 13.0 ± 0.3 kg CO2−eq/kg of LW sold (Table 2). The sale weight and weaning of the calves was 181.5 ± 3.1 kg with 6 months of age. The weight of the cows at weaning was 542.2 ± 10.0 kg. Calves' supplementation had a significant effect on the variables CF, total emissions, weight of calves, and total kg of calves sold (p < 0.01, Table 2). The decrease in CF reached 4% for supplementation scenarios due to the 10.2% increase in the sale weight of calves compared with the non-supplementation scenario. With regard to the total emissions, the use of the supplement increased the total CO2 emissions by 2.5%. However, due to increased total production (kilos sold), the value of the CF decreased. The inclusion of additives to corn silage had a significant effect on the CF, methane emissions, and total CO2 emissions (Table 2). Thus, the scenarios that used monensin decreased CF by 1.6%, methane emissions by 1.7%, and total CO2 emissions by 1.4%, compared with the scenarios without additive. The effect produced by the addition of canola oil was slightly higher, reaching values of 2.4, 3.0, and 2.5% for CF, methane emissions, and total CO2 emissions, respectively.
Table 2. Simulated carbon footprint (CF), methane emissions, total emissions, and productive variables according to calf supplementation (creep feeding) and additives (monensin and canola oil) used in the baseline scenario (stocking rate of 0.7 LU/ha).
Methane emissions represented between 83.7 and 83.3% of total GHG emissions, decreasing as supplementation increased (p < 0.01). The use of additives had an effect on the percentage of methane emissions (p < 0.01). The lowest percentage (83.3%) was obtained with canola oil. GHG emissions associated with pastures, supplementation, and direct N2O represented on average 6.5, 4.2, and 4.1% of total GHG emissions.
These results show that the effect of additives on methane emissions is greater than the increase in total emissions resulting from the inclusion of the supplement (with or without additives).
For the economic variables, supplementation had a significant effect on total costs, raising them by 8.2% at the highest level of supplementation (1 kg DM/calf/day) compared with the baseline scenario. However, no significant differences were observed between supplementation levels for operational and financial income (p < 0.01). The additive/calf supplementation interaction was not significant in any of the variables studied.
Experimental Scenarios
A total of 65 scenarios were evaluated. Of these, 60 (five stocking rates, four levels of supplementation, and two additives + supplement without additive) included the supplementation of the calves. Table 3 displays the results for the three factors evaluated (stocking rate, calf supplementation, and additive). The variables stocking rate and calf supplementation had a significant effect on the weight of the calves and the total kg LW of calves weaned/sold (Figure 2). Regarding the economic variables, a significant effect of calf supplementation was observed on all the variables associated with costs and variable gross income (p < 0.01). The additive/calf supplementation interaction was not significant in any of the variables studied.
Table 3. Cost and income variables concerning calf supplementation (creep feeding) and additives (monensin and canola oil) in the baseline scenario (stocking rate of 0.7 LU/ha).
Figure 2. Calf weight according to stocking rate and level of calf supplementation.
An interaction effect between calf supplementation and stocking rate was observed for CF, calf weight, and total production (p < 0.01) (Table 4). The quantification of carbon sequestration had a significant effect on the CF (p < 0.01) (Figure 3).
Table 4. Simulated CF, total emissions, and productive variables in relation to stocking rate (0.4, 0.55, 0.7, 0.85, and 1 LU/ha), calf supplementation (creep feeding), and additives (monensin and canola oil).
Figure 3. Effect of stocking rate (0.4, 0.55, 0.7, 0.85, and 1 LU/ha), calf supplementation (creep feeding), and carbon sequestration on carbon footprint (CF).
When carbon sequestration was not considered, CF averages decreased by 2.5 and 1.4% for canola oil and monensin, respectively. In the scenario with carbon sequestration, the average decreases in CF were 2.4% for canola oil and 1.3% for monensin (Table 3). Figure 4 shows the effect of the additives on the CF based on the calf supplementation and stocking rate in the estimates with and without consideration of carbon sequestration.
Figure 4. Effect of calf supplementation, additive stocking rate, and carbon sequestration on the carbon footprint (CF) of the simulated scenarios.
Figure 5 illustrates the distribution of total GHG emissions concerning the evaluated variables in the simulated scenarios. In the systems evaluated, the items representing the highest percentages of GHG emissions corresponded to enteric methane emissions, pastures, supplementation, and direct N2O, with average percentages of 83.5, 6.6, 4.2, and 4.1%, respectively. Effects of calf supplementation, additive, and stocking rate were observed on the percentage distribution of total GHG emissions (p < 0.05).
Figure 5. Total greenhouse gas (GHG) emissions according to stocking rate (0.4, 0.55, 0.7, 0.85, and 1 LU/ha), supplementation (creep feeding), and additive (monensin and canola oil).
For the economic variables (Table 5), results showed a significant effect for the supplementation, stocking rate, and the supplementation/stocking rate interaction on total cost, operational cost, total gross income, operational income, average total cost, average operational cost, and financial income (p < 0.01). The increase in stocking rate to 1 LU/ha increased total costs by 40.3% compared with the baseline scenario (0.7 LU/ha). However, due to the higher total production, increases of 39.6, 32.7, and 30.8% by total gross income, financial income, and operational income, respectively, were observed. The use of supplement generated an average increase of 1.9% in the total cost and 2.2% in the operational cost for each increase of 0.25 kg DM/calf/day. Given the increase in kilos sold due to supplementation, an increase occurred in the total income. Despite this, due to the increase in the average cost of production at higher supplementation levels, the operational income on average decreased by 0.3% product of supplementation. A significant effect of the stocking rate/additive interaction was also observed in the financial income and operational income (p < 0.01).
Table 5. Cost and income variables concerning supplementation and additives in the simulated scenarios.
Discussion
The CF values obtained in this research were within the ranges described by Florindo et al. (52), Ogino et al. (57), and Crosson et al. (58), considering that the carcass yield of the common breeds in southern Chile is on average 56% (59). Similarly, the average weights gained by calves at 6 months in the scenarios without supplementation were similar to those reported by Catrileo et al. (9). As observed in Table 2, the CF tended to present higher values when the carbon sequestration was considered. This was confirmed in Figure 6, where part of the CF of the baseline scenario (stocking rate of 0.7 LU/ha) arose from the loss of carbon from the pasture and soil. Furthermore, the scenario with carbon sequestration showed that as the stocking rate increases, the CF increases. These results coincided with those reported by Toro-Mujica et al. (31) in extensive sheep production systems in rain-fed areas. As Sarkar et al. (60) point out, maintaining productivity continues to be a challenge for producers, implying maintaining a slightly higher than optimal system load. Thus, increasing the stocking rate without the corresponding increase in supplementation raises grazing pressure, decreasing the available photosynthetic material, the amount of litter with the potential to incorporate carbon into the soil, and the production of roots (61), resulting in less carbon sequestration to mitigate emissions. As Chen et al. (62) pointed out, improved grazing management regimes that allow the renewal of degraded grassland and promote carbon storage in the soil are required. On the other hand, no clear trend appeared to increase or decrease the CF, as the stocking rate increased in the scenarios without carbon sequestration. This meant that the stocking rate effect on the CF was dependent on the initial value of the stocking rate. However, the expected trend is to increase the CF when stocking rates exceed the level of equilibrium of the production system. This trend was shown slightly when the stocking rate was increased in the scenario without supplementation or with low supplementation levels (Figure 3). With regard to the distribution of total emissions, the percentage of enteric methane was higher than the 44% reported by Samsonstuen et al. (63), consistent with the lower use of inputs in the systems studied.
Figure 6. Carbon footprint (CF) and distribution of greenhouse gas (GHG) emissions in simulated scenarios.
The significant effect of creep feeding on the weight of calves agreed with Carvalho et al. (64). The higher weight of calves with supplementation occurred because, in the supplementation months (summer/autumn), the pasture was in its reproductive stage, decreasing its nutritional value (65). Through a meta-analysis, Carvalho et al. (64) reported that the supplemental weight gain increases quadratically with supplementation (Figure 2). This trend was a consequence of the growth potential of the breed, of the potential feed intake, and the ability to substitute forage for a supplement without altering the functioning of the rumen. The greater gain in weaning weight resulted in reduced CF by kg of LW sold despite the increase in total emissions. In other words, the ratio between the extra production of emissions associated with supplementation and the extra kilos produced was less than the CF of the scenario without supplementation. Furthermore, an interaction between stocking rate and supplementation was observed. Thus, as the supplementation increased, the effect of the increase in the stocking rate on the weight of calves decreased (Figure 2). The interaction stocking rate/supplementation was also observed when the CF was analyzed considering the carbon sequestration. As the level of supplementation increased, the effect of the rise of the stocking rate on the CF decreased. Moreover, the weight of the cows was not affected by the level of supplementation of their calves, as reported by Carvalho et al. (64). This was because the calves decreased their pasture intake and not milk consumption as the level of supplementation rises (66). The pasture intake decreased with increasing supplementation levels (Figure 6). As the supplementation level increased, the proportion that carbon sequestration represented in the CF decreased, according to Crosson et al. (58). The additives incorporated into the supplementation, both monensin and canola oil, had slight but significant effects on decreasing the CF of the simulated systems (Table 3). Thus, the addition of canola oil into the supplement decreased the CF between 2 and 2.8% as the supplementation increased from 0.25 to 1 kg/calf/day. In the case of monensin, the decrease in the CF varied between 1.3 and 1.8%. In addition, the total enteric emissions showed decreases in percentages, ranging from 2.6 to 3.2% for canola oil and from 1.6 to 2.0% for monensin. The additives had a slight effect on the CF due to the fact that the supplement with the additive only represented between 6 and 21% of the total supplementation consumed and between 0.6 and 2.4% of the total feed consumed during the production cycle (for levels of 0.25 and 1 kg DM/calf/day, respectively). Then, a potential exists to reduce GHG emissions by incorporating these additives into the supplementation. For example, in the systems studied, this could be examined by increasing the period or level of supplementation in calves or by incorporating the additives into the supplementation of cows. This is an issue that needs to be analyzed in future research. In these studies, the secondary effects of increasing the level of supplementation need to be incorporated, since supplementation increases the mean retention time of feed in the rumen decreased, increased fractional passage rate, and decreased ruminal fermentation (67). This could alter the relative production of volatile fatty acids in the rumen and feed digestibility (68) and, thus, the efficiency of the feeding strategy. The incorporation of additives into the supplement had no effect on any of the economic variables examined. Three possible explanations emerged: a low level of additive per kg of feed, the periods and levels of supplementation evaluated, and the economic compensation of the cost of additive with the increase of the efficiency in the use of energy of the supplement. The increase in total costs due to supplementation did not give rise to lower financial or operational income. However, marginal profit decreased as the stocking rate increased. The above was consistent with the lower weight of cows and calves, and the subsequent lower milk production by decreasing DM availability in the pastures (69, 70).
Conclusions
The model proposed explained the effect of feed and management strategies on GHG emissions and production costs. The incorporation of additives into the supplementation of calves decreased the CF in the analyzed production system without increasing production costs. The most significant effects on CF were obtained with the highest levels of supplementation. Although the magnitude of the additives' effects was low, it can be increased as the proportion of the feed with additives is increased. In the first instance, including an additive to the cows' supplement emerges as an alternative. The stocking rate used in the baseline scenario was very close to the optimal level of the production system. This shows how extensive systems are managed to reach a point of equilibrium. In this way, the use of new production strategies in the cow–calf production systems requires an adaptation period and a comprehensive adjustment in the factors of production. The adaptation period, as well as the time required to reach an optimal CF, can be reduced using simulation models. For the simulated WCS baseline scenario (0.7 LU/ha), the CF (−5.1%) can be decreased and the financial income (5.8%) increased by supplementing calves with 0.75 kg DM/day without additives. Furthermore, by supplementing calves with 0.25 kg DM/day of monensin, a decrease in CF (−2.9%) and an increase in financial income (8.4%) were obtained. When an increase in stocking rate (0.85 LU/ha) was simulated, a slight increase in the CF (1.9%) allowed a significant increase (15.3%) in the financial income of the farm through the supplementation of calves with 0.75 kg DM/day and canola oil as additives.
Considering the international agreements to reduce GHG emissions, using a simulation model to evaluate feeding or management strategies in cow–calf production systems before its implementation on the farm appears to be an alternative that deserves further exploration. The secondary supplementation effects on rumen function and milk quality need to be examined in future simulation models.
Data Availability Statement
The raw data supporting this study are available on request to the corresponding author.
Author Contributions
PT-M: model design, experimentation, statistical analysis, data interpretation, and wrote the manuscript.
Funding
This work was supported by Comisión Nacional de Ciencia y Tecnología (CONICYT, Chile) through the projects REDI170084 and FONDECYT11190367.
Conflict of Interest
The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
The author also thanks Cecilia Espinoza, PhD, for her technical assistance (ORCID 0000-0002-6975).
References
1. FAO. Tackling Climate Change through Livestock - A global assessment of emissions and mitigation opportunities. (2013). Food and Agriculture Organization of the United Nations Rome, Italy.
2. IPCC. Guidelines for National Greenhouse Gas Inventories. (2006). Greenhouse Gas Inventories Programme, Intergovernmental Panel on Climate Change, IGES, Japan.
3. Cottle D, Eckard R. Global beef cattle methane emissions: yield prediction by cluster and meta-analyses. Anim Prod Sci. (2018) 58:2167–77. doi: 10.1071/AN17832
4. Niu M, Kebreab E, Hristov A, Oh J, Arndt C, Bannink A, et al. Prediction of enteric methane production, yield, and intensity in dairy cattle using an intercontinental database. Glob Change Biol. (2018) 24:3368–89. doi: 10.1111/gcb.14094
5. Patra A, Lalhriatpuii M. Development of statistical models for prediction of enteric methane emission from goats using nutrient composition and intake variables. Agric Ecosyst Environ. (2016) 215:89–99. doi: 10.1016/j.agee.2015.09.018
6. Beck P, Stewart C, Sims M, Gadberry M, Jennings J. Effects of stocking rate, forage management, and grazing management on performance and economics of cow-calf production in Southwest Arkansas. J Anim Sci. (2016) 94:3996–4005. doi: 10.2527/jas.2016-0634
7. Do Carmo M, Claramunt M, Carriquiry M, Soca P. Animal energetics in extensive grazing systems: rationality and results of research models to improve energy efficiency of beef cow-calf grazing Campos systems. J Anim Sci. (2016) 94:84–92. doi: 10.2527/jas.2016-0596
8. Rasby R, Warner J, Jenkins K, Klopfenstein T. Alternative cow-calf production systems: opportunities and challenges. J Anim Sci. (2016) 94:70. doi: 10.2527/msasas2016-150
9. Catrileo A, Toro P, Aguilar C, Vera R. Use of supplements and variation in the stocking rate in cow-calf systems on temperate pastures in Chile: a simulation approach. Anim Prod Sci. (2009) 49:1059–67. doi: 10.1071/EA08113
10. Ramírez-Restrepo C, Vera-Infanzón R, Rao I. Predicting methane emissions, animal-environmental metrics and carbon footprint from Brahman (Bos indicus) breeding herd systems based on long-term research on grazing of neotropical savanna and Brachiaria decumbens pastures. Agric Syst. (2020) 184:102892. doi: 10.1016/j.agsy.2020.102892
11. Basarab J, Baron V, López-Campos Ó, Aalhus J, Haugen-Kozyra K, Okine E. Greenhouse gas emissions from calf- and yearling-fed beef production systems, with and without the use of growth promotants. Animals. (2012) 2:195. doi: 10.3390/ani2020195
12. Chijioke O.gbonna A, Rudinow Saetnan E Meta-analysis of methane mitigation strategies: improved predictions of mitigation potentials and production implications. Iran J Appl Anim Sci. (2018) 8:567–73.
13. Ormond J. Geoengineering super low carbon cows: food and the corporate carbon economy in a low carbon world. Clim Change. (2020) 163:135–53. doi: 10.1007/s10584-020-02766-7
14. Villar M, Hegarty R, Nolan J, Godwin I, McPhee M. The effect of dietary nitrate and canola oil alone or in combination on fermentation, digesta kinetics and methane emissions from cattle. Anim Feed Sci Tech. (2020) 259:114294. doi: 10.1016/j.anifeedsci.2019.114294
15. Beck M, Thompson L, Williams G, Place S, Gunter S, Reuter R. Fat supplements differing in physical form improve performance but divergently influence methane emissions of grazing beef cattle. Anim Feed Sci Tech. (2019) 254:114210. doi: 10.1016/j.anifeedsci.2019.114210
16. Thompson L, Beck M, Gunter S, Williams G, Place S, Reuter R. An energy and monensin supplement reduces methane emission intensity of stocker cattle grazing winter wheat. Appl Anim Behav Sci. (2019) 35:433–40. doi: 10.15232/aas.2018-01841
17. Beauchemin K, McGinn S. Effects of various feed additives on the methane emissions from beef cattle. Int Congr Ser. (2006) 1293:152–5. doi: 10.1016/j.ics.2006.01.042
18. Guyader J, Eugène M, Doreau M, Morgavi D, Gérard C, Martin C. Tea saponin reduced methanogenesis in vitro but increased methane yield in lactating dairy cows. J. Dairy Sci. (2017) 100:1845–55. doi: 10.3168/jds.2016-11644
19. Roque B, Salwen J, Kinley R, Kebreab E. Inclusion of Asparagopsis armata in lactating dairy cows' diet reduces enteric methane emission by over 50 percent. J Clean Prod. (2019) 234:132–8. doi: 10.1016/j.jclepro.2019.06.193
20. Stefenoni H Räisänen S Cueva S Wasson D Lage C Melgar A . Effects of the macroalga Asparagopsis taxiformis and oregano leaves on methane emission, rumen fermentation, and lactational performance of dairy cows. J Dairy Sci. (2021) 104:4157–73. doi: 10.3168/jds.2020-19686
21. Lurette A, Aubron C, Moulin C. A simple model to assess the sensitivity of grassland dairy systems to scenarios of seasonal biomass production variability. Comput Electro Agric. (2013) 93:27–36. doi: 10.1016/j.compag.2013.01.008
22. White R, Brady M, Capper J, McNamara J, Johnson K. Cow–calf reproductive, genetic, and nutritional management to improve the sustainability of whole beef production systems. J Anim Sci. (2015) 93:3197–211. doi: 10.2527/jas.2014-8800
23. Rotz A, Asem-Hiablie S, Place S, Thoma G. Environmental footprints of beef cattle production in the United States. Agric Syst. (2019) 169:1–13. doi: 10.1016/j.agsy.2018.11.005
24. Moraes L, Fadel J, Castillo A, Casper D, Tricarico J, Kebreab E. Modeling the trade-off between diet costs and methane emissions: a goal programming approach. J Dairy Sci. (2015) 98:5557–71. doi: 10.3168/jds.2014-9138
25. Legesse G, Small J, Scott S, Crow G, Block H, Alemu A, et al. Predictions of enteric methane emissions for various summer pasture and winter feeding strategies for cow calf production. Anim Feed Sci Tech. (2011) 166–7:678–87. doi: 10.1016/j.anifeedsci.2011.04.082
26. INE. Instituto Nacional de EstadísticasMicrodatos Encuestas Ganaderas Bovinas 2015. (2015). Available online at: https://www.ine.cl (accessed June 10, 2019).
27. INE. Instituto Nacional de Estadísticas. Data from: Censo Agropecuario 2007. (2007). Base de microdatos. Instituto Nacional de Estadísticas, Santiago, Chile.
28. Rojas C, Romero O. Sistema de crianza de Hereford, utilizando festuca con trébol subterráneo en el valle de la IX Región. Agricultura Técnica (Chile). (1990) 50:379–85.
29. González V, Tapia M. Manual bovino de carne. Instituto de Desarrollo Agropecuario - Instituto de Investigaciones Agropecuarias. Boletín INIA N°04 (2019). Available online at: https://www.inia.cl/wp-content/uploads/ManualesdeProduccion/04%20Manual%20Bovino.pdf (accessed October 5, 2007).
30. FAO. Greenhouse Gas Emissions from the Dairy Sector - A Life Cycle Assessment. (2010). Food and Agriculture Organization of the United Nations, Rome, Italy.
31. Toro-Mujica P, Aguilar C, Vera, Bas RF. Carbon footprint of sheep production systems in semi-arid zone of Chile: a simulation-based approach of productive scenarios and precipitation patterns. Agric Sci. (2017) 157:22–38. doi: 10.1016/j.agsy.2017.06.012
32. Sykes A, Topp C, Wilson R, Reid G, Rees R. A comparison of farm-level greenhouse gas calculators in their application on beef production systems. J Clean Prod. (2017) 164:398–409. doi: 10.1016/j.jclepro.2017.06.197
33. Cerri C, Moreira C, Alves P, Raucci G, de Almeida Castigioni B, Mello F, et al. Assessing the carbon footprint of beef cattle in Brazil: a case study with 22 farms in the State of Mato Grosso. J Clean Prod. (2016) 112:2593–600. doi: 10.1016/j.jclepro.2015.10.072
34. Ripoll-Bosch R, de Boer I, Bernués A, Vellinga T. Accounting for multifunctionality of sheep farming in the carbon footprint of lamb: a comparison of three contrasting Mediterranean systems. Agric Syst. (2013) 116:60–8. doi: 10.1016/j.agsy.2012.11.002
35. British Standards Institute. PAS 2050:2011 Specification for the Assessment of the Life Cycle Greenhouse Gas Emissions of Goods and Services. London: British Standards Institute (2011).
36. Rotz CA, Asem-Hiablie S, Dillon J, Bonifacio H. Cradle-to-farm gate environmental footprints of beef cattle production in Kansas, Oklahoma, and Texas. J Anim Sci. (2015) 93:2509–19. doi: 10.2527/jas.2014-8809
37. Mogensen L, Kristensen T, Nielsen NI, Spleth P, Henriksson M, Swensson C, et al. Greenhouse gas emissions from beef production systems in Denmark and Sweden. Livest Sci. (2015) 174:126–43. doi: 10.1016/j.livsci.2015.01.021
38. Batalla I, Knudsen M, Mogensen L, Hierro O, Pinto M, Hermansen J. Carbon footprint of milk from sheep farming systems in Northern Spain including soil carbon sequestration in grasslands. J Clean Prod. (2015) 104:121–9. doi: 10.1016/j.jclepro.2015.05.043
39. Petersen B, Knudsen M, Hermansen J, Halberg N. An approach to include soil carbon changes in life cycle assessments. J Clean Prod. (2013) 52:217–24. doi: 10.1016/j.jclepro.2013.03.007
40. Soussana J, Tallec T, Blanfort V. Mitigating the greenhouse gas balance of ruminant production systems through carbon sequestration in grasslands. Animal. (2010) 4:334–50. doi: 10.1017/S1751731109990784
41. Lal R. Carbon emission from farm operations. Environ Int. (2004) 30:981–90. doi: 10.1016/j.envint.2004.03.005
44. AFRC. Energy and Protein Requirements of Ruminants: An Advisory Manual. (1993). Agricultural and Food Research Council. Technical Committee on Responses to Nutrients, Great Britain.
45. Grinstead C, Snell L. Introduction to Probability. 2nd revised ed. Providence, RI: American Mathematical Society (1997).
46. Edwards-Jones G, Plassmann K, Harris I. Carbon footprinting of lamb and beef production systems: insights from an empirical analysis of farms in Wales, UK. J Agric Sci. (2009) 147:707–19. doi: 10.1017/S0021859609990165
47. Saunders C, Barber A. Comparative Energy and Greenhouse Gas Emissions of New Zealand's and the UK's Dairy Industry. Philadelphia, PA: Lincoln University (2007).
48. NRC. Nutrient Requirements of Beef Cattle: Eighth Revised Edition. National Research Council. Washington, DC: The National Academies Press (2016).
49. Adom F, Maes A, Workman C, Clayton-Nierderman Z, Thoma G, Shonnard D. Regional carbon footprint analysis of dairy feeds for milk production in the USA. Int J Life Cycle Assess. (2012) 17:520–34. doi: 10.1007/s11367-012-0386-y
50. Lal R. Residue management, conservation tillage and soil restoration for mitigating greenhouse effect by CO2-enrichment. Soil Till Res. (1997) 43:81–107. doi: 10.1016/S0167-1987(97)00036-6
51. IPCC. 2019 Refinement to the 2006 IPCC Guidelines for National Greenhouse Gas Inventories. In: Calvo Buendia E, Tanabe K, Kranjc A, Baasansuren J, Fukuda MS, et al., editors. Geneva: IPCC (2019).
52. Florindo T, de Medeiros Florindo G, Talamini E, da Costa J, Ruviaro C. Carbon footprint and life cycle costing of beef cattle in the Brazilian midwest. J Clean Prod. (2017) 147:119–29. doi: 10.1016/j.jclepro.2017.01.021
53. Toro-Mujica P, Vera R, Vargas-Bello-Pérez E, Pinedo P, Bas F. Trends and drivers of change of pastoral beef production systems in a mediterranean-temperate climate zone of Chile. Animals. (2019) 9:1135. doi: 10.3390/ani9121135
54. Ranga Niroshan Appuhamy J, Strathe A, Jayasundara S, Wagner-Riddle C, Dijkstra J, France J, et al. Anti-methanogenic effects of monensin in dairy and beef cattle: a meta-analysis. J Dairy Sci. (2013) 96:5161–73. doi: 10.3168/jds.2012-5923
55. R Core Team. R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing (2019). Available online at: https://www.R-project.org/ (accessed January 29, 2021).
56. Kleijnen J. Validation of models: statistical techniques and data availability. In: Farrington HBNPA, Sturrock DT, Evans G, editors. Proceedings of the 1999 Winter Simulation Conference (AZ: Phoenix W) (1999). p. 647–54.
57. Ogino A, Sommart K, Subepang S, Mitsumori M, Hayashi K, Yamashita T, et al. Environmental impacts of extensive and intensive beef production systems in Thailand evaluated by life cycle assessment. J Clean Prod. (2016) 112:22e31. doi: 10.1016/j.jclepro.2015.08.110
58. Crosson P, Shalloo L, O'Brien D, Lanigan G, Foley P, Boland T, et al. A review of whole farm systems models of greenhouse gas emissions from beef and dairy cattle production systems. Anim Feed Sci Tech. (2011) 166–7:29–45. doi: 10.1016/j.anifeedsci.2011.04.001
59. Hänsch A, Larraín R. Carbon footprint of a typical grass-based beef production system in Chile. in Journal of Abstracts of the 65th Annual Meeting of the European Federation of Animal Science. Wageningen: Wageningen Academic Publishers (2014). p. 138.
60. Sarkar R, Corriher-Olson V, Long C, Somenahally A. Challenges and potentials for soil organic carbon sequestration in forage and grazing systems. Rangeland Ecol Manag. (2020) 73:786–95. doi: 10.1016/j.rama.2020.04.002
61. Klumpp K, Fontaine S, Attard E, Le Roux X, Gleixner G, Soussana J. Grazing triggers soil carbon loss by altering plant roots and their control on soil microbial community. J Ecol. (2009) 97:876–85. doi: 10.1111/j.1365-2745.2009.01549.x
62. Chen W, Huang D, Liu N, Zhang Y, Badgery W, Wang X, et al. Improved grazing management may increase soil carbon sequestration in temperate steppe. Sci Rep. (2015) 5:10892. doi: 10.1038/srep10892
63. Samsonstuen S, Åby B, Crosson P, Beauchemin K, Wetlesen M, Bonesmo H, et al. Variability in greenhouse gas emission intensity of semi-intensive suckler cow beef production systems. Livest Sci. (2020) 239:104091. doi: 10.1016/j.livsci.2020.104091
64. Carvalho V, Paulino M, Detmann E, Filho S, Lopes S, Rennó L, et al. A meta-analysis of the effects of creep-feeding supplementation on performance and nutritional characteristics by beef calves grazing on tropical pastures. Livest Sci. (2019) 227:175–82. doi: 10.1016/j.livsci.2019.07.009
65. Fonseca M, Valadares Filho S, Henriques L, Paulino P, Detmann E, Fonseca E, et al. Nutritional requirements of nursing Nellore calves. Rev Bras Zoot. (2012) 41:1212–21. doi: 10.1590/S1516-35982012000500019
66. van Wyngaard J, Meeske R, Erasmus L. Effect of concentrate level on enteric methane emissions, production performance, and rumen fermentation of Jersey cows grazing kikuyu-dominant pasture during summer. J Dairy Sci. (2018) 101:9954–66. doi: 10.3168/jds.2017-14327
67. van Gastelen S, Dijkstra J, Bannink A. Are dietary strategies to mitigate enteric methane emission equally effective across dairy cattle, beef cattle, and sheep?. J Dairy Sci. (2019) 102:6109–30. doi: 10.3168/jds.2018-15785
68. Chaokaur A, Nishida T, Phaowphaisal I, Sommart K. Effects of feeding level on methane emissions and energy utilization of Brahman cattle in the tropics. Agric Ecosyst Environ. (2015) 199:225–30. doi: 10.1016/j.agee.2014.09.014
69. Gunter S. Effects of stocking and supplementation rates on the performance of beef steers grazing mixed-grass prairie during the winter. Appl Anim Sci. (2019) 35:641–51. doi: 10.15232/aas.2019-01864
Keywords: simulation model, carbon footprint, feed strategies, stocking rate, additive
Citation: Toro-Mujica P (2021) Evaluation of Feed Strategies and Changes of Stocking Rate to Decrease the Carbon Footprint in a Traditional Cow-Calf System: A Simulation Model. Front. Vet. Sci. 8:587168. doi: 10.3389/fvets.2021.587168
Received: 25 July 2020; Accepted: 30 April 2021;
Published: 09 June 2021.
Edited by:
Paul Crosson, Teagasc, IrelandReviewed by:
Jorge Avila Stagno, University of Concepcion, ChileAklilu W. Alemu, Agriculture and Agri-Food Canada (AAFC), Canada
Copyright © 2021 Toro-Mujica. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Paula Toro-Mujica, paula.toro@uoh.cl