SciELO - Scientific Electronic Library Online

 
vol.28 issueNSPE1Performance evaluation of Uruguayan dairy farming systemsEffect of untreated and stabilized dairy effluent applications on soil fertility and associated health risks author indexsubject indexarticles search
Home Pagealphabetic serial listing  

Services on Demand

Journal

Article

Related links

Share


Agrociencia Uruguay

On-line version ISSN 2730-5066

Agrocienc. Urug. vol.28 no.nspe1 Montevideo  2024  Epub Sep 20, 2024

https://doi.org/10.31285/agro.28.1243 

Articles

Balancing nitrogen at the farm gate: Economic-environmental sustainability trade-off in pastoral dairy systems of Uruguay

Balanceando el nitrógeno en la portera del predio: Sostenibilidad económica y ambiental en los sistemas lecheros pastoriles de Uruguay

Equilibrando o nitrogênio no portão da propriedade: Sustentabilidade econômica e ambiental nos sistemas leiteiros pastoris do Uruguai

1Instituto Nacional de Investigación Agropecuaria (INIA), Colonia, Uruguay, sstirling@inia.org.uy

2Universidad de la República, Facultad de Agronomía, Paysandú, Uruguay

3Instituto Nacional de la Leche, Montevideo, Uruguay

4University of Tennessee, Department of Biosystems Engineering and Soil Science, Knoxville, USA


Abstract:

Uruguay´s dairy can potentially enhance milk productivity competitively, but intensification risks elevating nitrogen (N) surplus, heightening environmental concerns. This study quantified farm-gate N inputs and outputs, calculating N surplus (input-output) and N use efficiency (NUE=output/input) for 17 commercial modal dairy systems identified in the 2014 and 2019 national surveys and six prospective intensified systems based on experimental pastoral farmlets achieving near-maximal rainfed productivity. Current dairy systems maintained N surplus at 71 kg N ha-1 between 2014 and 2019 (range: 44-97 kg N ha-1) while improving NUE from 28.3 to 30.5% (range: 20-35%). Intensification increased N surplus without necessarily reducing NUE. Our analyses highlight three aspects: (i) comparatively low N surplus of current Uruguayan dairy, (ii) nonlinear links between N surplus and stocking rate, feed intake, milk productivity and operating profit, and (iii) inequality between dairy systems in their contribution to national dairy N surplus reflects mainly disparity in farm size. These insights underscore the crucial need for understanding the actual fate of N surpluses: nitrate leaching, ammonia volatilisation, N2 denitrification, or N accumulation in soil organic matter. This is an unavoidable requisite for designing management practices and policies able to effectively optimise the economic and environmental sustainability of Uruguayan dairy.

Keywords: grazing systems; whole-farm balance; N losses; N surplus; N use efficiency

Resumen:

La lechería uruguaya presenta la oportunidad de mejorar la competitividad en la productividad de leche. Sin embargo, la intensificación aumenta el riesgo de elevar el excedente de nitrógeno (N) y la contaminación ambiental. Este estudio evaluó las entradas y salidas de N predial, N excedente (entrada-salida) y eficiencia de uso de N (NUE=salida/entrada) de 17 sistemas modales (encuestas nacionales 2014-2019), y 6 sistemas intensificados basados en estudios experimentales de alta productividad de leche por hectárea en secano. Los sistemas lecheros actuales mantuvieron un N excedente de 71 kg N ha-1 entre 2014 y 2019 (rango: 44-97 kg N ha-1), y mejoraron la NUE de 28,3 a 30,5% (rango: 20-35%). De nuestro análisis se destacan tres puntos: (i) el N excedente comparativamente bajo de la lechería uruguaya actual, (ii) las relaciones no lineales entre N excedente y carga animal, consumo de alimento, productividad de la leche e ingreso de capital, y (iii) la contribución desigual entre los sistemas lecheros al excedente nacional de N de la lechería refleja principalmente la disparidad en el tamaño de los predios. Estos hallazgos evidencian la necesidad crucial de entender el destino real de los excedentes de N: lixiviación de nitratos, volatilización de amoníaco, desnitrificación de N2 o acumulación de N en materia orgánica del suelo. Este es un requisito ineludible para diseñar prácticas y políticas de gestión capaces de optimizar efectivamente la sostenibilidad económica y ambiental de la lechería uruguaya.

Palabras clave: balance predial de N; contaminación; sistemas lecheros

Resumo:

A pecuária leiteira uruguaia apresenta a oportunidade de melhorar a competitividade na produtividade do leite. No entanto, a intensificação aumenta o risco de elevar o excesso de nitrogênio (N) e a poluição ambiental. Este estudo avaliou as entradas e saídas de N na propriedade, o excesso de N (entrada-saída) e a eficiência de uso de N (NUE = saída/entrada) de 17 sistemas modais (pesquisas nacionais 2014-2019) e seis sistemas intensificados com base em estudos experimentais de alta produtividade de leite por hectare em sequeiro. Os sistemas leiteiros atuais mantiveram um excesso de N de 71 kg N ha-1entre 2014 e 2019 (faixa: 44-97 kg N ha-1) e melhoraram a NUE de 28,3% para 30,5% (faixa: 20-35%). De nossa análise, destacam-se três pontos: (i) o excesso de N comparativamente baixo na pecuária leiteira uruguaia atual, (ii) relações não lineares entre o excesso de N e carga animal, consumo de alimentos, produtividade do leite e entrada de capital, e (iii) a contribuição desigual entre os sistemas leiteiros para o excesso nacional de N da pecuária leiteira reflete principalmente a disparidade no tamanho das propriedades. Essas descobertas evidenciam a necessidade crucial de entender o destino real dos excessos de N: lixiviação de nitratos, volatilização de amônia, desnitrificação para N2 ou acumulação de N na matéria orgânica do solo. Este é um requisito inegociável para projetar práticas e políticas de gestão capazes de otimizar efetivamente a sustentabilidade econômica e ambiental da pecuária leiteira uruguaia.

Palavras-chave: balanço de N; poluição; sistemas leiteiros

1. Introduction

Intensive dairy production is debated because of its roles in resource utilisation, nitrogen (N) pollution of the environment, greenhouse gas emissions, and their impacts on landscape and biodiversity1. Concerning N cycling, the intensification of dairy production has been linked to an increase in N surplus, a decrease in nitrogen use efficiency (NUE), and the accumulation of N in the soil. This occurs because NUE in dairy production is limited by the biological capacity of cows to convert N from feed into milk, and by the ability of crops and pastures to transform soil N into grains, forage, and other agricultural products2. The surplus N in the system that is not absorbed by plants or utilised by animals can result in detrimental environmental consequences, including water pollution from nitrate leaching, global warming due to nitrous oxide emissions, and air pollution from ammonia emissions3.

In this regard, the strong association between intensificationviahigher stocking rates and pasture production per hectare and N surplus has been widely reported in intensified pasture-based dairy systems such as those from New Zealand4, Australia5, The Netherlands6 or Ireland7. This situation has led to water pollution due to the leaching of N from urine and fertilizer deposited on pasture in countries with predominantly permeable soils, such as Ireland and New Zealand 8)(9) 10.

The Uruguayan dairy presents the opportunity to increase productivity per hectare but with the challenge of maintaining competitiveness and ensuring the sustainability of natural resources 11)(12) . Uruguay exports around 70% of total milk production to the international dairy market13 and is positioned as the seventh largest milk exporting country worldwide14. Therefore, the future of the dairy sector relies on its ability to remain competitive domestically and internationally. Dairy farm systems in Uruguay present the opportunity to increase milk production through improved utilisation of home-grown forage by increasing the stocking rate11. National research studies show that a two-fold increase in stocking rate from the industry average is biologically feasible and significantly improves dairy systems' biophysical and economic performance 15)(16) 17. These results have encouraged and provided a clear direction for the intensification of commercial dairy systems. These proposed prospective intensification strategies combine stocking rates, cow genotypes, and feeding strategies. Overall, they are characterised by increased stocking rates, home-grown forage harvest per hectare, and reliance on N fertilizers and imported supplements. In this regard, the environmental impact of such strategies in terms of N use remains unclear 11)(18) .

The objective of the present study was to assess the magnitude of farm-gate N surpluses and NUE in dairy farms from Uruguay. The aim was to consolidate information regarding the performance of current and prospective dairy systems in terms of basic N indicators to help optimise the trade-off between economic and environmental sustainability.

2. Materials and methods

The farm-gate balance of 23 representative Uruguayan dairy farm systems was evaluated. The farm systems evaluated were: (i) 17 commercial modal dairy systems classified from two national dairy surveys, and (ii) 6 prospective intensified systems scaled up from experimental farmlets data.

2.1 Data sources

2.1.1 Modal dairy systems

The modal dairy systems are typologies of representative systems of the Uruguayan dairy sector developed by the National Milk Institute (INALE) Technical Team, based on farm scale (annual milk production), land productivity (L ha-1), and pasture harvest (kg DM ha-1) (the latter only applies for the 2014 modal systems) (Table S1, Supplementary material).

Data from national dairy surveys conducted in 2014 and 2019 by INALE and the Ministry of Livestock, Agriculture and Fisheries (MGAP) were used to develop these typologies. The survey included farms located in the main dairy regions of Uruguay (i.e., Canelones, Colonia, Flores, Florida, Paysandú, Río Negro, San José and Soriano). A total of 314 and 294 farms were surveyed in 2014 and 2019, respectively. This accounted for more than 80% of national overall milk production; therefore, we will refer to the results based on this data as the national average 19)(20) . The information surveyed included farm area, animal stocks and reconciliation, land use, milk and live weight productivity, feeding, pasture utilisation, silage and concentrate production and purchases and fertilizer purchases.

Seven modal systems were defined from the 2014 survey (i.e., ML1, ML2, ML4, ML6, ML8, ML10, ML12) and 10 from the 2019 survey (i.e., M1, M2B, M2A, M3B, M3A, M4B, M4A, M5B, M5A, M6). The economic analysis was performed on the modal dairy farms for the 2019-20 financial year. The classification criteria and detailed information about each modal dairy system are presented in Table S1 and Table S2. The national average was calculated as the milk production-weighted of all modal dairy farm systems surveyed in 2014 and 2019.

2.1.2 Experimental farmlets

The intensified systems evaluated in this study arise from the “CRS” and “10MIL” farmlet studies. Both farmlet studies assessed the biophysical feasibility and economic result of intensification options to increase milk production per hectare well above the average national dairy farm. The information from each farmlet (three-year data) was scaled up to the average Uruguayan modal dairy farm (i.e., M4; 145 ha farm) as part of a study carried out by INALE to assess the economics of technological research proposal adoption. It was simulated that the average Uruguayan modal dairy farm (M4) “adopted” the CRS and 10MIL research systems (i.e., farm system design, stocking rate, feeding strategies, labour, infrastructure, fertilization rates, investments, etc.) to achieve similar productive and economic results. This methodology enabled us to compare the modal dairy systems and the experimental farmlets.

A detailed description of the “10MIL” farmlet study is provided by Stirling and others16. Briefly, the farmlets were located at INIA La Estanzuela (34°20′S, 57°41′W, Uruguay) and evaluated for three years (Seasons 2017/18, 2018/19 and 2019/20) the biophysical performance of intensification strategies based on increasing home-grown forage harvest to sustain a two-fold increase in stocking rate relative to the national average. Specifically, this study evaluated two feeding strategies with varying proportions of grazing in the annual feeding budget (grass fixed (GFix) and grass maximum (GMax)) and two Holstein Friesian cow genotypes (New Zealand (NZHF) or North American Holstein Friesian (NAHF)). The farmlets evaluated were GMax-NAHF, GMax-NZHF, GFix-NAHF, and GFix-NZHF.

The “CRS” farmlet study was conducted at the Centro Regional Sur research station, Agronomy Faculty, Canelones, Uruguay21 (34°36.8′ S, 56°13.1′ W). The experiment involved four whole-farm systems evaluated for three years (Seasons 2016/17, 2017/18 and 2018/19). This study evaluated two stocking rates (SR): 1.5 or 2.0 milking cows per ha-1 of milking platform (i.e., 1.4 and 1.7 cows per ha-1 of dairy area) combined by two contrasting residual sward heights: 4 cm residual sward height all year round (Low Residual, LR) or seasonally variable residual sward height (6 cm for autumn and winter, 9 cm for spring and 12 cm for summer) (High Residual, HR). Thus, four treatments resulted from combining the two factors: 1.5LR, 1.5HR, 2.0LR and 2.0HR. For this study, the average SR treatments were analysed as no apparent differences were observed between the LR and HR treatments for the “CRS” farmlets.

2.2 Farm systems characteristics

The modal dairy systems and scaled-up experimental farmlets will be referred to as dairy farm systems. The farm systems descriptions with key productive and economic performance indicators for the modal and intensified dairy systems are given in Table 1. The national average dairy systems (2014 and 2019) present an average of 1.0 cows ha-1 and produce 415 kg milk solids ha-1. The intensified systems (“CRS” and “10MIL”) present, on average, a two fold increase in stocking rate (1.8 cows ha-1) and solids productivity per ha (1,000 kg milk solids ha-1) regarding the national average.

All the systems studied presented common features of typical dairy farming in Uruguay. Around 60% of the total farm area is occupied by long-term mixed pastures that last 3-4 years, and the rest is occupied by annual pastures (i.e., oats, ryegrass) and crops (i.e., sorghum, maize) destined for grazing and silage or grain production. While most of the silage is produced on-farm, these systems rely on purchased concentrate (82-100% bought-in). Generally, supplements are offered on a feed pad or sacrifice paddock where cows are confined during pasture shortage seasons, and when access to pasture is constrained due to weather events. Inorganic fertilizer is imported and applied to pastures and crops in these systems.

Table 1: Description of the national average dairy systems (2014 and 2019) and the “CRS” and “10MIL” intensified dairy systems 

1The national average was calculated as the milk production-weighted of all modal dairy farm systems surveyed in 2014 and 2019.

2.3 Farm-gate balance calculation

The balance presented in this study involved calculating N inputs and N outputs at the farm-gate according to the guidance document developed by the EU N Expert Panel22, as it is widely used in science and policy. Nitrogen surplus was calculated as the difference between total N inputs and total N outputs and expressed on a per-hectare basis. Nitrogen use efficiency was calculated as the ratio of total N outputs divided by total N inputs and expressed as a percentage.

All N input and output data were collected in individual farms and expressed in kg N ha-1 of farm area yr-1. The unit of land considered for this study was the dairy area, defined as the area destined for the milking cows plus the area destined for the dry stock.

2.3.1 Nitrogen imports

Inputs comprised N entering the farm as mineral fertilizer, imported concentrates and silage, biological N fixation and atmospheric deposition. Nitrogen imported in mineral fertilizers was calculated from the amount of chemical fertilizer purchased. Standard N concentrations, as provided by commercial suppliers, were used. Nitrogen entering the farm in concentrate and silage feeds was established by multiplying the total quantity purchased by the crude protein content, using specific values for the modal dairy systems and the farmlet studies, divided by 6.2523.

For biological N2 fixation calculation, a value of 36.5 kg N fixed per Mg of DM legume consumed by animals was used24 (Supplement S3). The amount of consumed legume DM was estimated by calculating total pasture harvested (Mg DM ha-1 yr-1) as the sum of grazed pasture and pasture harvested as forage reserves and then affected by the percentage of the area with grass/legume mixed pastures, which was reported in surveys and measured in farmlets. Afterwards, 40% of legume DM in the mixed pasture was assumed for modal systems, as the survey did not include this data (Supplement S3). For farmlets, the measured proportion of legume DM in mixed pastures was used. Nitrogen atmospheric deposition values were taken for two locations in Uruguay25.

Nitrogen from silage surplus was calculated as the difference between total silage imported and produced and silage consumed and the average N concentration of silage. It was considered a “negative input” rather than an output, according to the guidelines suggested by de Klein and others9.

2.3.2 Nitrogen exports

Nitrogen outputs included milk and live weight production. Nitrogen exported in milk was calculated from the annual amount of milk produced and the crude protein milk concentration divided by 6.3826. Nitrogen export in livestock was calculated from the difference in live weight between the replacement heifers entering the farm and culled cows leaving the farm multiplied by a standard 0.024% N concentration in live weight of adult animals26.

2.3.3 Statistical analysis

Regression analysis was undertaken to describe the relationship between the components of the farm-gate N balance and between N surplus and the productive and economic variables used to characterise the systems. The goodness of fit of fitted models was assessed using the coefficient of determination (R2) adjusted by the number of parameters and the residual mean squared error. Power models fitted the data with similar or better goodness-of-fit than linear models. The power regression model used was y = c + a x b , where y is the dependent variable, x is the independent variable, and the fitted parameters a, b and c are the scaling coefficient, the exponent coefficient, and the intercept, respectively. The exponent coefficient is particularly relevant as it indicates the degree of non-linearity (how far it is from 1) and the shape of the curvature (increasing when >1, decreasing when <1). All analyses were carried out in TBL Curve v5.01 (Systat), using the Levenberg-Marquardt algorithm to iteratively minimise squared residuals.

2.3.4 Contribution of the 2019 modal dairy farms to total national dairy N surplus

The contribution of each modal dairy farm (2019) to the total N surplus in Uruguay was assessed by plotting the cumulative proportion of total milk production and total number of dairy farms in Uruguay vs. the cumulative proportion of total N surplus.

Gini plots were generated by linear interpolation of the cumulative proportions of milk production, number of farms and N surplus. These plots provide a visual and quantitative depiction of the degree of inequality of dairy modal systems to the national dairy N surplus. The Gini coefficient was calculated27, ranging from 0 (perfect equality) to 1 (perfect inequality).

2.3.5 Targets for nitrogen performance indicators

Nitrogen input and output values of the assessed systems were mapped onto the framework proposed by the EU Nitrogen Expert Panel22. Furthermore, we mapped pastoral dairy systems from primary exporting dairy countries in the same plot9. This approach allows the analysis of farm performance simultaneously in terms of its NUE, N output, and N surplus, providing a graphical assessment of the distance between actual and possible target values for each indicator for each system.

For N output, two minimum thresholds were defined based on the milk solids productivity per hectare that a dairy farm system must achieve to be profitable in Uruguay. The long-term sustainability of a dairy farm is achieved when the operating profit can remunerate the land and capital factors. For the calculation of remuneration: (i) land is treated as if it is all leased at an average market value of 160 US$ ha-1; (ii) capital (livestock, machinery, and improvements) is remunerated at an annual interest rate of 4%, which amounts to 108 US$ ha-1. According to this, a dairy farm's minimum acceptable operating profit in the 2021/22 period was calculated as 268 US$ ha-1. For comparison, the operating profit of agricultural systems -which would compete for land with dairy systems- averaged 251 US$ ha-1 over the last ten years (FUCREA, Federación Uruguaya de grupos CREA; pers. comm.). Therefore, dairy systems below the 250-270 US$ ha-1 threshold (i.e., N output below 19 and 21 kg N ha-1) are deemed to have their economic sustainability compromised.

Maximum and minimum thresholds for NUE were set at 20% and 60%, according to targets proposed for mixed crop-livestock systems, where both crop and livestock products are included28. Mixed crop-livestock production systems will have different target values than crop production systems because of the increased risks of N losses from animal manures and the low availability of organically bound N to crops. A NUE higher than 60% risks inducing soil mining (i.e., N removal exceeds N input, declining soil fertility and plant yield), and a NUE below 20% is deemed too low (i.e., N inputs to the systems exceed total N demand), even accounting for the inherently low conversion efficiency of the animal feed protein-N in milk and meat protein-N.

The N surplus target value was not set because there is no information to relate it to the magnitude of actual losses of reactive N species to the environment, particularly water courses. As a tentative reference, we set ranges of N surplus from <80 to 100-160 to >160 kg N surplus ha -1(9) 22. This is based on NUE results (see below): systems with N surplus below 80 kg N ha-1 could increase their NUE, but NUE would decrease for surpluses above 160 kg N ha-1.

3. Results

3.1 Farm-gate balance

The farm-gate balance of the average Uruguayan dairy systems was similar in 2014 and 2019, as no major changes were observed in inputs, outputs, and surplus between the two surveys, although NUE did improve (Table 2, Table S2 and Table S5). In these systems, inputs were similarly distributed between fertilizer, imported feeds and biological N fixation and N atmospheric deposition, which accounted for 32, 31 and 37% of total N inputs, respectively. Outputs were amply dominated by milk production, as live weight production made a minor contribution. The seven modal dairy farms identified in 2019 imported on average 103 kg N ha-1 (ranging from 59 to 137 kg N ha-1) and exported 32 kg N ha-1 (ranging from 13 to 42 kg N ha-1), thus presenting a N surplus of 71 kg N ha-1 (44 and 97 kg N ha-1) (Table 2, Table S5).

There was a disparity between modal dairy farms (2019) in their contribution to the national dairy N surplus, depending on how their contribution is assessed. When considered in terms of produced milk, all modal systems contributed rather proportionally, and thus, the cumulative milk production scaled almost equally with cumulative N surplus (Gini coefficient = 0.21: Figure 1A). However, modal dairy systems differed greatly in size, and therefore just 21% of the number of farms (i.e., those grouped in the M5B, M5A and M6 modal systems) accounted for 74% of the total national dairy N surplus (Gini coefficient = 0.69; Figure 1B).

The farm-gate balance of the “CRS” and “10MIL” intensified systems had a N surplus 79 to 280% higher than the average 2019 modal dairy system, respectively (Table 2). In the “CRS” systems, NUE remained similar to the 2019 average dairy system, but NUE decreased to 22% in the “10MIL” systems. In the intensified systems, fertilizer and imported feeds became the main N inputs, accounting for 37 and 43% of total N inputs, respectively. Although intensification changed little absolute biological N fixation, this source became relatively less important.

Farm-gate N surplus increased with intensification (Table 2, Figure 2 and Figure 3). Close positive relationships were observed between N surplus and N input (R2=0.998) and N output (R2=0.965). However, while N surplus increased proportionally more as N input increased (b=1.26; Figure 2A, Figure 2D), N output increased proportionally less (b=0.33; Figure 2B). Therefore, N surplus and N output had a markedly nonlinear link (b=2.14). As a direct consequence of these allometric relationships, NUE was not associated with N inputs in a simple way (R2=0.00; Figure 2C). Instead, NUE showed a complex relationship with N inputs, increasing up to 100 kg N ha-1 and decreasing above 200 kg N ha-1.

The link between N surplus and all indicators of intensification was also nonlinear. It increased at an accelerating rate as milk solids productivity increased (b=1.77; Figure 3A). This pattern was also clear for stocking rate (b=1.93), total feed intake (b=2.33), as well as for operating profit (b=1.42; Figure 3B-D). This means that increments in all these variables resulted in progressively more significant N surpluses.

Table 2: Farm-gate N balance of the national average dairy systems (2014 and 2019) and the “CRS” and “10MIL” intensified systems 

1The national average was calculated as the milk production-weighted of all modal dairy farm systems surveyed in 2014 and 2019.

Figure 1: Gini plot of cumulative proportional total milk production (A) for total number of dairy farms vs. (B) cumulative proportional N surplus for Uruguayan modal dairy systems in 2019 

Figure 2: Relationships between the components of farm-gate balance: N input, N output, N surplus and N use efficiency (NUE)  

Figure 3: Relationships between N surplus and several parameters associated with intensification: milk solids productivity (A), stocking rate (B), total feed intake (C) and operating profit (D)  

3.2 Mapping results onto the EU Nitrogen Expert Panel input-output framework

When the relationship between N inputs and outputs for each farm system evaluated was mapped onto the EU Nitrogen Expert Panel input-output framework, NUE and N output targets were achieved by some, but not all, assessed systems (Figure 4). On one hand, the most extensive modal systems had adequate NUE but too low N outputs (i.e., N output below the minimum required to be profitable; 19 and 21 kg N ha-1). On the other hand, “10MIL” systems presented high N output and NUE within the desired range, but their N surplus was substantially greater than 160 kg N ha-1. The same pattern was observed in pastoral dairy systems from New Zealand, Australia, Europe, and the US 2)(7) 9.

Figure 4: Annual N input and N output values of dairy systems mapped onto the framework proposed by the EU N Expert Panel22. The two dashed lines encompass the 20-60% range of desired NUE, as lower values risk increasing N pollution and higher values risk mining soil N. The horizontal dotted lines represent the N output associated with the minimum productivity to reach an operating profitable to either (i) remunerate the land and capital factors for current Uruguayan dairy systems or (ii) equal the average operating profit of current agricultural systems that compete for land with dairy systems. Finally, the shaded grey gradient represents the current uncertainty regarding what N surplus would lead to substantial losses of harmful reactive N species. 

4. Discussion

This study amalgamated data from the 2014 and 2019 national surveys of current commercial dairy farms with data from simulated systems adopting intensification strategies based on the outcomes of experimental research farmlets. This unique database of farm-gate N balances of dairy systems constitutes the first assessment of the overall current national dairy situation and the impact of prospective intensified systems.

4.1 Current and future dairy farming systems: the impact of intensification

Current pastoral dairy systems in Uruguay have lower stocking rates and productivity than pasture-based dairy systems from New Zealand, Australia, Northern Europe, or the US (c.f.Table 1vs. McDowell and others29; Luo & Ledgard30; Ros and others31). Likewise, the N surplus of current Uruguayan dairy systems, at 71 kg N ha-1, is at the lower end of the values reported for pastoral dairy systems worldwide9, less than half the national average dairy N surplus of New Zealand30 (186-281 kg N ha-1), Australia32 (156 kg N ha-1), Ireland33 (155 kg N ha-1), or the Netherlands34 (174-208kg N ha-1). However, the NUE of current Uruguayan dairy farms (20-35%) is within the 20-40% range reported for dairy systems worldwide9. Therefore, the performance of the Uruguayan dairy in terms of farm-gate N balance, as measured in 2019, is relatively better than comparable systems. Furthermore, NUE increased by two per cent points between 2014 and 2019 (Table 1).

Intensification of pastoral dairying in Uruguay, via enhanced forage consumption and increased stocking rate, has been shown to improve profitability in commercial11 and experimental dairy farms 16)(17) . For instance, the intensified “CRS” and “10MIL” systems achieved two- to three-fold rises in operating profit compared to the national dairy systems. However, economic performance alone does not guarantee the sustainable intensification of dairy 1)(35) . Increasing milk production per hectare while maintaining low production costs to remain internationally competitive may entail greater environmental risks.

Farm-gate N balance shows that the “CRS” and “10MIL” intensified systems imported on average 187 and 356 kg N ha-1 and exported 59 and 85 kg N ha-1. This resulted in an average N surplus of 128 and 271 kg N ha-1, respectively (Table 2, Figure 4), comparable to those of intensified pasture-based dairy systems in leading exporting countries 29)(30) 31.

These results confirm the general trend reported in the literature that intensification of dairy systems increases stocking rate, milk production and operating profit, but also farm-gate N surplus 2)(7)9)(22) 36. Therefore, it is safe to assume that environmental problems associated with intensified dairy systems in more developed countries may also emerge at some stage due to intensification in Uruguay.

As the average dairy farm intensifies by adopting “CRS” and “10MIL” strategies, N inputs to the farm system increased due to higher use of fertilizer and imported feeds required to support larger feed demand of higher stocking rates. Notably, the NUE of the intensified systems was neither consistently lower nor related to the composition of N inputs. Opposite to negative trends observed in other countries9, in the present study, NUE peaked at intermediate values of N input, a consequence of the shape of the allometric relationship between N outputs and inputs (Figure 2C).

4.2 Mapping Uruguayan dairy farms onto the EU Nitrogen Expert Panel framework

Optimal N performance targets would ideally combine the efficient use of high N inputs to sustain animal productivity while minimising N loss9. Plotting systems in terms of their N input and output allows us to assess their compliance (or not) with three concurrent targets22. First, a N output target is defined by the minimum profitability required to be competitive (e.g., against other land uses). Second, NUE targets set by achievable efficiency under local agroecological conditions. Third, a maximum N surplus target to avoid losses to the environment of harmful reactive N species. This is the most uncertain target and should be defined for the specific conditions according to environmental sensitivity and policy37.

Mapping Uruguayan farms in the N input-output framework shows that current dairy systems present a lower N surplus than intensified systems, and, thus, less potential risk of harmful N losses. However, 40% of these systems (i.e., M1, M2B, M3B, M4B), which account for 20% of national milk production, had N outputs below the minimum 21 kg N ha-1 required to be competitive in Uruguay (Figure 4 ). Notably, de Klein and others9 placed that threshold at 80 kg N ha-1 for other pastoral countries, perhaps due to higher values of land.

At the other end of the distribution, the most intensified systems presented the highest output per hectare and NUE above 20%, but also the highest N surplus. Finally, some current and intensified systems plot in the space comprising the desirable range of outcomes for all targets. This suggests that intensified dairy systems can, in fact, be more profitable than the current average system while complying with N performance targets.

4.3 Strategies for sustainable intensification

Different strategies can improve N performance. These strategies can be grouped into intensification, extensification, increasing efficiency and avoiding soil degradation22 (Figure S2). Our results suggest that, in Uruguay, dairy farms with low productivity could increase N inputs and improve their performance without major increases in risk of environmental impact. These systems -which comprised 50% of dairy farms in 2019- could reach the desirable range of outcomes (unshaded part of the figure) following a “traditional” intensification pathway.

For the most intensified and profitable systems (e.g., “10MIL”), the challenge is to increase their efficiency (i.e., NUE) and reduce N surplus. This requires improving, first, internal N recycling, especially N in animal excreta9, to reduce fertilizer inputs, and second, the amount of harvested grass per unit of N fertilizer applied38. These adjustments allow either greater N output in saleable products for the same N input, or lower N inputs while maintaining N outputs.

Provided appropriate management practices are adopted to account for those two key aspects, intensified pasture-based dairy may be achieved without major increases in N surplus. In the Netherlands, N surplus remained unchanged in commercial dairy farms between 2013 and 2015 due to the implementation of environmental regulations limiting N application and making mandatory the export of manure surpluses34. Similarly, intensified Irish dairy farms reported positive impacts of regulations upon N surplus and NUE due to less inorganic fertilizer N use and improved manure application timing39.

4.4 Need for further research

Farm-gate N balances are the simplest indicator of a dairy farm’s environmental performance. Its simplicity, however, is also its main drawback because N surplus is not linearly or unequivocally related to harmful N losses to the environment, i.e. nitrate leaching or ammonia volatilisation 4)(40)41)(42) 43. Indeed, the relationship between N balances and groundwater concentrations has been found to be poor, unclear, or non-existent 4)(40) 43. Moreover, part of the N surplus might not be lost to the environment but stored in soil organic matter, bringing both environmental and productive benefits44. Therefore, N surplus is best interpreted as an indicator of potential risk 2)(9) .

Pollution from reactive N is arguably the planetary boundary surpassed to a larger extent by food production systems45. Whether the high N surplus associated with intensive dairy systems would end up increasing water pollution -the most problematic local impact of N excess- requires research. Future studies, possibly combining field measurements with simulation models, should quantitatively assess the relationship between N surplus and the biochemical identity, spatial distribution, and temporal dynamics of N losses. Without this, there is no solid basis for designing either (i) effective N management options for the spatial and temporal hotspots that contribute disproportionally to N environmental risk, or (ii) sound targeted policies able to optimise the compromise between economic and environmental sustainability of dairy systems.

5. Conclusions

This study consolidates data on the N balance of current dairy systems in Uruguay and offers insights into the potential effects of intensification. First, current dairy systems in Uruguay present a comparatively low N surplus compared to other pasture-based systems worldwide, but at the same time, many of those current systems have N outputs that compromise their economic sustainability. Second, the intensification of dairy systems achieving higher profits entails increases in farm-gate N surplus. The pronounced nonlinearity observed in the surplus-intensification relationship is particularly concerning.

The insights gleaned from this study may help guide technological advancements and policy formulation, fostering innovative and sustainable farming practices that contribute to achieving a more environmentally responsible and economically robust dairy sector. Specifically, we deem it essential that future research assess the fate of N surplus, i.e., its partitioning into the various N loss pathways, which carry quite different environmental consequences, and accumulation in soil organic matter, which carries environmental and productive benefits. Finally, the large inequality in the contribution of individual modal dairy systems to the national N surplus can be considered if effective policy is sought after.

Acknowledgements:

This work was funded by the National Institute of Agriculture Research (INIA). The authors acknowledge the dairy farm staff and students for their assistance during the farmlet trials at INIA La Estanzuela and CRS Faculty of Agronomy. FL and FAL were financially supported by PROCISUR and FONTAGRO through the project “Uso de leguminosas en sistemas ganaderos sudamericanos”.

References

1. Clay N, Garnett T, Lorimer J. Dairy intensification: Drivers, impacts and alternatives. Ambio. 2020;49(1):35-48. Doi: 10.1007/s13280-019-01177-y. [ Links ]

2. Powell JM, Gourley CJP, Rotz CA, Weaver DM. Nitrogen use efficiency: A potential performance indicator and policy tool for dairy farms. Environ Sci Policy. 2010;13(3):217-28. Doi: 10.1016/j.envsci.2010.03.007. [ Links ]

3. Cameron KC, Di HJ, Moir JL. Nitrogen losses from the soil/plant system: a review. Ann Appl Biol. 2013;162(2):145-73. Doi: 10.1111/aab.12014. [ Links ]

4. Beukes PC, Scarsbrook MR, Gregorini P, Romera AJ, Clark DA, Catto W. The relationship between milk production and farm-gate nitrogen surplus for the Waikato region, New Zealand. J Environ Manage. 2012;93(1):44-51. Doi: 10.1016/j.jenvman.2011.08.013. [ Links ]

5. Gourley CJP, Dougherty WJ, Weaver DM, Aarons SR, Awty IM, Gibson DM, Hannah MC, Smith AP, Peverill KI. Farm-scale nitrogen, phosphorus, potassium and sulfur balances and use efficiencies on Australian dairy farms. Anim Prod Sci. 2012;52(10):929-44. Doi: 10.1071/AN11337. [ Links ]

6. Oenema J, Oenema O. Unraveling feed and nutrient use effciencies in grassland based dairy farms. Front Sustain Food Syst. 2022;6:846561. Doi: 10.3389/fsufs.2022.846561. [ Links ]

7. Treacy M, Humphreys J, Namara KM, Browne R, Teagasc CJW. Farm-gate nitrogen balances on intensive dairy farms in the south west of Ireland. Ir J Agric Food Res. 2008;47:105-17. [ Links ]

8. Clark DA, Caradus JR, Monaghan RM, Sharp P, Thorrold BS. Issues and options for future dairy farming in New Zealand. J Agric Res. 2007;50(2):203-21. Doi: 10.1080/00288230709510291. [ Links ]

9. de Klein CAM, Monaghan RM, Ledgard SF, Shepherd M. A system’s perspective on the effectiveness of measures to mitigate the environmental impacts of nitrogen losses from pastoral dairy farming. In: Proceedings of the 4th Australasian Dairy Science Symposium. Canterbury: Lincoln University; 2010. p. 14-28. [ Links ]

10. Dillon P, Shalloo L, Murphy D, O’Brien D, Richards K, O’Donovan M, Ruelle E. The impact of nitrogen management strategies within grass based dairy systems (Internet). Oak Park: Teagasc; 2021 (cited 2024 Feb 29). 49p. Available from: Available from: https://www.teagasc.ie/publications/2021/the-impact-of-nitrogen-management-strategies-within-grass-based-dairy-systems.php . [ Links ]

11. Fariña SR, Chilibroste P. Opportunities and challenges for the growth of milk production from pasture: the case of farm systems in Uruguay. Agric Syst. 2019;176:102631. Doi: 10.1016/j.agsy.2019.05.001. [ Links ]

12. Chilibroste P. A major challenge for the Uruguayan dairy industry: sustainable growth. Agrocienc Urug. 2021;25(2):e970. Doi: 10.31285/AGRO.25.970. [ Links ]

13. Ministerio de Ganadería, Agricultura y Pesca, DIEA (UY). Anuario estadístico agropecuario 2020. Montevideo: MGAP; 2020 (cited 2024 Feb 29). Available from: Available from: https://bit.ly/3UXE8s6Links ]

14. OECD; FAO. OECD-FAO agricultural outlook 2021-2030. Paris: OECD Publishing; 2021. 337p. Doi: 10.1787/19428846-en. [ Links ]

15. Ortega G, Nuñez T, Custodio D, Mello R, López Y, Chilibroste P. Effect of stocking rate on pasture production and utilisation on a grazing dairy system during winter and spring. In: Abstracts of the 2018 American Dairy Science Association Annual meeting. Knoxville: American Dairy Science Association; 2018. p. 257. [ Links ]

16. Stirling S, Delaby L, Mendoza A, Fariña S. Intensification strategies for temperate hot-summer grazing dairy systems in South America: effects of feeding strategy and cow genotype. J Dairy Sci. 2021;104(12):12647-63. Doi: 10.3168/jds.2021-20507. [ Links ]

17. Stirling S. Estrategias de intensificación para sistemas lecheros pastoriles en Uruguay: análisis biofísico y del potencial impacto ambiental (doctoral’s thesis). Montevideo (UY): Universidad de la República, Facultad de Veterinaria; 2022. 155p. [ Links ]

18. Darré E, Llanos E, Astigarraga L, Cadenazzi M, Picasso V. Do pasture-based mixed dairy systems with higher milk production have lower environmental impacts? A Uruguayan case study. J Agric Res . 2021;64(3):444-62. Doi: 10.1080/00288233.2020.1750433. [ Links ]

19. INALE. Encuesta 2014: resultados definitivos agrupados (Internet). Montevideo: INALE; 2019 (cited 2024 Feb 29). Available from: Available from: https://www.inale.org/estadisticas/encuesta-2014-resultados-definitivos-agrupados/Links ]

20. INALE. Encuesta Lechera 2019: resultados definitivos (Internet). Montevideo: INALE ; 2021 (cited 2024 Feb 29). Available from: Available from: https://www.inale.org/estadisticas/encuesta-lechera-2019-resultados-definitivos/Links ]

21. Ortega G, Chilibroste P, Garrido JM, Waller A, Fariña SR, Lattanzi FA. Monitoring herbage mass and pasture growth rate of large grazing areas: a comparison of the correspondence, cost and reliability of indirect methods. J Agric Sci. 2023;161:502-11. Doi: 10.1017/S0021859623000333. [ Links ]

22. EU Nitrogen Expert Panel. Nitrogen use efficiency (NUE) an indicator for the utilization of nitrogen in agriculture and foodsystems. Wageningen: Wageningen University; 2015. 47p. [ Links ]

23. Mcdonald P, Edwards RA, Greenhalgh JFD, Morgan CA. Animal nutrition. Englewood Cliffs: Prentice Hall; 1995. 696p. [ Links ]

24. Lussich F. Variabilidad de la fijación biológica de nitrógeno de leguminosas forrajeras en Uruguay: posibles causas y consecuencias nutricionales (master’s thesis). Montevideo (UY): Universidad de la República, Facultad de Agronomía; 2020. 102p. [ Links ]

25. Carnelos DA, Portela SI, Jobbágy EG, Jackson RB, Di Bella CM, Panario D, Fagúndez C, Piñeiro-Guerra JM, Grion L, Piñeiro G. A first record of bulk atmospheric deposition patterns of major ions in southern South America. Biogeochemistry. 2019;144(3):261-71. Doi: 10.1007/s10533-019-00584-3. [ Links ]

26. Agricultural Research Council. The nutrient requirements of ruminant livestock. Oxon: CAB International; 1994. 351p. [ Links ]

27. Gini C. Measurement of inequality of incomes. Econ J. 1921;31(121):124-6. Doi: 10.2307/2223319. [ Links ]

28. Brentrup F, Pallière C. Nitrogen use efficiency as an agro-environmental indicator. In: Proceedings of the OECD workshop on OECD agri-environmental indicators: lessons learned and future directions, 23-26 March, 2010. Leysin: OECD; 2010. p. 23-6. [ Links ]

29. McDowell RW, Rotz CA, Oenema J, Macintosh KA. Limiting grazing periods combined with proper housing can reduce nutrient losses from dairy systems. Nat Food. 2022;3(12):1065-74. Doi: 10.1038/s43016-022-00644-2. [ Links ]

30. Luo J, Ledgard S. New Zealand dairy farm systems and key environmental effects. Front Agric Sci Eng. 2021;8(1):148-58. Doi: 10.15302/J-FASE-2020372. [ Links ]

31. Ros MBH, Godber OF, Olivo AJ, Reed KF, Ketterings QM. Key nitrogen and phosphorus performance indicators derived from farm-gate mass balances on dairies. J Dairy Sci . 2023;106(5):3268-86. Doi: 10.3168/jds.2022-22297. [ Links ]

32. Stott KJ, Gourley CJP. Intensification, nitrogen use and recovery in grazing-based dairy systems. Agric Syst . 2016;144:101-12. Doi: 10.1016/j.agsy.2016.01.003. [ Links ]

33. Buckley C, Wall DP, Moran B, O’Neill S, Murphy PNC. Farm gate level nitrogen balance and use efficiency changes post implementation of the EU nitrates Directive. Nutr Cycl Agroecosystems. 2016;104(1):1-13. Doi: 10.1007/s10705-015-9753-y. [ Links ]

34. Oenema J, Oenema O. Intensification of grassland-based dairy production and its impacts on land, nitrogen and phosphorus use efficiencies. Front Agric Sci Eng . 2021;8(1):130-47. Doi: 10.15302/J-FASE-2020376. [ Links ]

35. Garnett T, Appleby MC, Balmford A, Bateman IJ, Benton TG, Bloomer P, Burlingame B, Dawkins M, Dolan L, Fraser D, Herrero M, Hoffmann I, Smith P, Thornton PK, Toulmin C, Vermeulen SJ, Godfray HC. Agriculture: Sustainable intensification in agriculture: premises and policies. Science. 2013;341(6141):33-4. Doi: 10.1126/science.1234485. [ Links ]

36. Haas G, Deittert C, Köpke U. Farm-gate nutrient balance assessment of organic dairy farms at different intensity levels in Germany. Renew Agric Food Syst. 2007;22(3):223-32. Doi: 10.1017/S1742170507001780. [ Links ]

37. Quemada M, Lassaletta L, Jensen LS, Godinot O, Brentrup F, Buckley C, Foray S, Hvid SK, Oenema J, Richards KG, Oenema O. Exploring nitrogen indicators of farm performance among farm types across several European case studies. Agric Syst . 2020;177:102689. Doi: 10.1016/j.agsy.2019.102689. [ Links ]

38. Rawnsley RP, Smith AP, Christie KM, Harrison MT, Eckard RJ. Current and future direction of nitrogen fertiliser use in Australian grazing systems. Crop Pasture Sci. 2019;70(12):1034-43. Doi: 10.1071/CP18566. [ Links ]

39. Mihailescu E, Murphy PNC, Ryan W, Casey IA, Humphreys J. Nitrogen balance and use efficiency on twenty-one intensive grass-based dairy farms in the South of Ireland. J Agric Sci . 2014;152(5):843-59. Doi: 10.1017/S0021859614000045. [ Links ]

40. Humphreys J, O’Connell K, Casey IA. Nitrogen flows and balances in four grassland-based systems of dairy production on a clay-loam soil in a moist temperate climate. Grass Forage Sci. 2008;63(4):467-80. Doi: 10.1111/j.1365-2494.2008.00660.x. [ Links ]

41. Pinxterhuis JB, Beare MH, Edwards GR, Collins RP, Dillon P, Oenema J. Eco-efficient pasture based dairy farm systems: a comparison of New Zealand, The Netherlands and Ireland. In: Grassland and forages in high output dairy farming systems. Wageningen: Grassland Science in Europe; 2015. p. 349-66. [ Links ]

42. Hoekstra NJ, Schulte RPO, Forrestal PJ, Hennessy D, Krol DJ, Lanigan GJ, Müller C, Shalloo L, Wall DP, Richards KG. Scenarios to limit environmental nitrogen losses from dairy expansion. Sci Total Environ. 2020;707:134606. Doi: 10.1016/j.scitotenv.2019.134606. [ Links ]

43. Burchill W, Lanigan GJ, Li D, Williams M, Humphreys J. A system N balance for a pasture-based system of dairy production under moist maritime climatic conditions. Agric Ecosyst Environ. 2016;220:202-10. Doi: 10.1016/j.agee.2015.12.022. [ Links ]

44. Soussana JF, Lemaire G. Coupling carbon and nitrogen cycles for environmentally sustainable intensification of grasslands and crop-livestock systems. Agric Ecosyst Environ . 2014;190:9-17. Doi: 10.1016/j.agee.2013.10.012. [ Links ]

45. Leip A, Bodirsky BL, Kugelberg S. The role of nitrogen in achieving sustainable food systems for healthy diets. Glob Food Sec. 2021;28:100408. Doi: 10.1016/j.gfs.2020.100408. [ Links ]

46. Ledgard SF, Steele KW. Biological nitrogen fixation in mixed legume/grass pastures. Plant and Soil. 1992;141:137-53. Doi: 10.1007/BF00011314. [ Links ]

Transparency of data The entire data set that supports the results of this study was published in the article itself.

Author contribution statement SS: Conceptualization; Data curation; Formal analysis; Investigation; Writing - original draft FL: Data curation; Writing - review & editing GO: Data curation; Writing - review & editing ALM: Data curation AP: Data curation; Writing - review & editing JA: Data curation; Writing - review & editing GG: Data curation SF: Conceptualization; Writing - review & editing PC: Conceptualization; Writing - review & editing FAL: Conceptualization; Supervision; Writing - review & editing

Editor The following editor approved this article: Mariana Carriquiry (https://orcid.org/0000-0002-1600-2591) Universidad de la República, Montevideo, Uruguay.

Supplementary material

Table S1: Farm systems description of the 2014 Modal dairy systems (ML1, ML2, ML4, ML6, ML8, ML10, ML12) and the milk-, area- and farm-weighted average 

Table S2: Farm systems description of the 2019 Modal dairy systems (M1, M2B, M2A, M3B, M3A, M4B, M4A, M5B, M5A, M6) and the milk-, area- and farm-weighted average 

Supplement S3

Estimation of N input via biological N2 fixation by forage legumes To estimate N input via biological N2 fixation, we assumed, first, that legumes accounted for 40% of consumed DM (including pasture grazed and harvested for reserves), and second, that 36.5 kg of N were fixed for each Mg of consumed legume DM. Assumption 1. We simulated the average proportion of legume DM for three frequent sequences in dairy systems of Uruguay (based on information from INALE19). These comprised annual winter forage crops (VI, pure grass, typically ryegrass or oats), mixed perennial pastures (PP, grass/legume) and annual summer forage crops (VV, pure grass, typically sudangrass). The main differences between sequences are the botanical composition and length of the phase with perennial pasture: 3.5 years for alfalfa+dactylis and for tall fescue+white clover/lotus corninulatus mixtures, and 1.5 years for red clover+ryegrass/bromus mixtures. For longer pastures, we simulated evolutions with either high or low legume content. The average proportion of legumes across years and sequences ranged between 31 and 46%. We used 40% as a simplifying factor. Assumption 2. We estimated the amount of N fixed per Mg of consumed DM for the three most frequent legume species in dairy farms in Uruguay: alfalfa, white clover and red clover (based on information from INALE19), assuming a 60-20-20 contribution, respectively. Data for the amount of N fixed per unit shoot DM for each species was taken from Lussich24. This coefficient was then affected by an estimated proportion of consumed shoot DM, an estimated allocation of N to belowground organs46, and an estimated recycling factor of N from senescent leaves so that Kg N fix/Mg DM consumed = kg N fix/Mg produced DM: consumed DM/produced DM × (1+allocation of N belowground) × adjustment N mobilisation.

Table S4: Farm-gate N balance of the 2014 modal dairy systems (ML1, ML2, ML4, ML6, ML8, ML10, ML12) and the milk-, area- and farm-weighted average 

Table S5: Farm-gate N balance of the 2019 Modal dairy systems (M1, M2B, M2A, M3B, M3A, M4B, M4A, M5B, M5A, M6) and the milk-, area- and farm-weighted average 

Figure S1: Farm-gate N balance of the 2014 and 2019 national average dairy systems and the “CRS” and “10MIL” intensified systems 

Figure S2: Ideal pathways that any farm could follow to enter the characteristic operating space (unshaded area): Intensification, extensification, increasing efficiency and avoiding soil degradation within the N input - N output framework proposed by the EU N Expert Panel22  

Received: September 04, 2023; Accepted: February 20, 2024

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License