# Modelling the effects of stocking rate, soil type, agroclimate location and nitrogen input on the grass DM yield and forage self-sufficiency of Irish grass-based dairy production systems

Elodie Ruelle, Luc Delaby, Laurence Shalloo, Michael O'Donovan, Deirdre Hennessy, M. Egan, B. Horan, P. Dillon

► To cite this version:

Elodie Ruelle, Luc Delaby, Laurence Shalloo, Michael O'Donovan, Deirdre Hennessy, et al. Modelling the effects of stocking rate, soil type, agroclimate location and nitrogen input on the grass DM yield and forage self-sufficiency of Irish grass-based dairy production systems. Journal of Agricultural Science, 2022, pp.1 - 15. 10.1017/s0021859622000272. hal-03741389

HAL Id: hal-03741389

https://hal.inrae.fr/hal-03741389

Submitted on 1 Aug 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
<figure><img src="image_2.png" /></figure>

Modelling Animal Systems
Research Paper

**Cite this article:** Ruelle E, Delaby L, Shalloo L, O'Donovan M, Hennessy D, Egan M, Horan B, Dillon P (2022). Modelling the effects of stocking rate, soil type, agroclimate location and nitrogen input on the grass DM yield and forage self-sufficiency of Irish grass-based dairy production systems. *The Journal of Agricultural Science* 1–15. https://doi.org/ 10.1017/S0021859622000272

Received: 19 July 2021

Revised: 16 May 2022

Accepted: 31 May 2022

**Key words:**

Farm system modelling; forage self-sufficiency;
grass nitrogen response; Irish dairy system;
stocking rate

**Author for correspondence:**

E. Ruelle,

E-mail: Elodie.ruelle@teagasc.ie

Modelling the effects of stocking rate, soil type,
agroclimate location and nitrogen input on the
grass DM yield and forage self-sufficiency of
Irish grass-based dairy production systems

E. Ruelle¹ ID, L. Delaby², L. Shalloo¹, M. O'Donovan¹, D. Hennessy¹ ID, M. Egan¹, B. Horan¹ and P. Dillon¹

¹Teagasc, Animal & Grassland Research and Innovation Centre, Moorepark, Fermoy, Cork, Ireland and ²INRAE, Physiologie, Environnement et Génétique pour l'Animal et les Systèmes d'Elevage, 35590 Saint-Gilles, France

Abstract

In pasture-based dairy production systems, identifying the appropriate stocking rate (SR; cows/ha) based on the farm grass growth is a key strategic decision for driving the overall farm business. This paper investigates a number of scenarios examining the effects of SR (2–3 cows/ha (0.25 unit changes)), annual nitrogen (N) fertilizer application rates (0–300 kg N/ha (50 kg/ha unit changes)), soil type (heavy and a free-draining soil) and agroclimate location ((south and northeast of Ireland) across 16 years) on pasture growth and forage self-sufficiency using the pasture-based herd dynamic milk model merged with the Moorepark St Gilles grass growth model. The modelled outputs were grass growth, grass dry matter intake, silage harvested and offered, overall farm forage self-sufficiency and N surplus. The model outputs calculated that annual grass yield increased from 9436 kg DM/ha/year when 0 kg N/ha/year was applied to 14 996 kg DM/ha/year when 300 kg N/ha/year were applied, with an average N response of 18.4 kg DM/kg N applied (range of 9.9–27.7 kg DM/kg N applied). Systems stocked at 2.5 cows/ha and applying 250–300 kg N fertilizer/ha/year were self-sufficient for forage. As N input was reduced from 250 kg N/ha/year, farm forage self-sufficiency declined, as did farm N surplus. The results showed that a reduction in N fertilizer application of 50 kg/ha/year will require a reduction in an SR of 0.18 cows/ha to maintain self-sufficiency (R² = 0.90).

Introduction

The sustainability of Irish dairy and drstock farms is dependent on increased productivity and improved efficiency of the conversion of grazed pasture to animal products (Hanrahan *et al.*, 2017). The removal of European Union milk quota regime has allowed the Irish dairy industry to increase cow numbers facilitated by increased land area associated with dairying (Kelly *et al.*, 2020). The Department of Agriculture, Food and the Marine Food Wise 2025 strategy set out as its guiding principal that environmental protection and economic competitiveness would be equally considered and complementary; one was not to be achieved at the expense of the other. While the Irish dairy industry expands, the European Green Deal announced by the European Commission in December 2019 has also identified the requirement to reduce nutrient losses by 50% and fertilizer inputs by 20% by 2030 compared to the base year of 2018.

Nitrogen (N) is essential for plants and it plays a key role in food production being among the most important crop yield-limiting factors in the world together with water (Mueller *et al.*, 2012). At the same time, N losses contribute to climate change and lead to environmental pollution which is harmful for the functioning of ecosystems and human health (Galloway *et al.*, 2008; Fowler *et al.*, 2013). Despite the positive effects of N fertilizer on agricultural production (Whitehead, 1995), the efficiency with which N is used is variable (Watson and Atkinson, 1999). Recent studies have suggested that current N losses from agriculture to the environment are too high (Steffen *et al.*, 2015) and continuing population and consumption growth in the next few decades will increase the demand for N fertilizer and increase N losses unless significant improvements are made in the whole food production–consumption chain (Godfray *et al.*, 2010; Sutton *et al.*, 2013; Mogollón *et al.*, 2018).

Within temperate grass-based dairy systems, there is a reliance on matching feed supply (grass growth) and feed demand (managed through calving date, stocking rate (SR) and strategic supplementation). There is a high reliance on climate conditions associated with the production of animal feed harvested by ruminants. The foundation of this production system is productive grassland and efficient grass utilization over a prolonged grazing season (Shalloo

© The Author(s), 2022. Published by
Cambridge University Press. This is an Open
Access article, distributed under the terms of
the Creative Commons Attribution licence
(http://creativecommons.org/licenses/by/4.0/),
which permits unrestricted re-use, distribution
and reproduction, provided the original article
is properly cited.

<figure><img src="image_3.png" /></figure>

CAMBRIDGE
UNIVERSITY PRESS
and Hanrahan, 2020). Such systems are unique in balancing pas-
ture utilization and maintaining high animal dry matter intake
(DMI) and performance, while reducing and minimizing
imported feed and fertilizer. The beneficial impacts of high SR
on grazing system productivity are widely reported (McCarthy
et al., 2016; Coffey et al., 2018); the impacts of SR on the environ-
ment need to be fully considered. Previously increased SR was
associated with increased chemical fertilizer use and reliance on
imported feed, greater nutrient surpluses and reduced
nutrient-use efficiency, resulting in increased losses to ground-
water and the general environment (Hennessy et al., 2020).
Both McCarthy et al. (2015) and Roche et al. (2016) investigated
the direct effect of SR on nitrate-N leaching and both studies
reported either stable or declining nitrate-N leaching with increas-
ing SR; however, strictly no additional N fertilizer or supplements
were introduced at higher SR.

Currently, Irish dairy farms that have a derogation from the Nitrates Directive (91/676/EEC) can legally apply up to 250 kg N fertilizer/year in conjunction with an organic N loading of a maximum of 250 kg N/ha/year, corresponding to an annual SR of 2.75 cows/ha. However, water quality has been declining in recent years (EPA, 2019) and so there is an overall requirement to reduce chemical N fertilizer use, the impact of variation in the rate of N applied needs to be investigated. The objective of this paper is to investigate the influence of chemical N fertilizer application rate on grass yield, stock carrying capacity (through forage self-sufficiency at different SRs) and N surplus. The pasture-based herd dynamic milk (PBHDM) model (Ruelle et al., 2015) merged with the Moorepark St Gilles grass growth (MoSt GG) model (Ruelle et al., 2018a) was used to simulate a 40 ha farm. The number of cows, soil type, fertilizer application and weather were changed for each simulation, resulting in 140 simulations over 16 climate years. The main model outputs were daily grass growth, cow DMI and silage fed, and farm N surplus.

Materials and methods

Overall summary of the models

The PBHDM model (Ruelle et al., 2015, 2016) simulates a dairy farm at the paddock and cow level. The model simulates each individual animal at grazing and is dependent on animal characteristics but also on grass availability with a reduction in animal DMI during the defoliation process included in the model. Within the model, grassland management-based decision rules are included to ensure realistic simulations. Animals move from one paddock to another based on a fixed residency time in the paddocks or based on an objective post-grazing height, which can be either fixed or dependent on the pre-grazing height. The farm grass cover is evaluated daily and is compared with the animal requirement. In situations where there is a farm grass cover excess on the overall farm, some paddocks can be allocated for silage conservation rather than grazed and vice versa. In the case of a grass deficit (subject to the management rules defined by the user), forage or concentrate supplementation can be added to the diet (Ruelle et al., 2015).

The PBHDM model was evaluated by comparing model out-
puts with two grazing systems in France and in Ireland (Ruelle
et al., 2015). Model outputs that were compared to the actual
experimental data were milk production, pre- and post-grazing
sward height and feed supplementation levels. These outputs
were all compared on a weekly basis while paddock residence

time and total grass harvested and conserved as grass silage was
evaluated over the entire grazing season.

Grass growth in the PBHDM model is predicted by the MoSt
GG model (Ruelle et al., 2018b) which is fully integrated within
the PBHDM (Ruelle et al., 2015). The MoSt GG model (Ruelle
et al., 2018a) is a dynamic model developed in C++ describing
the grass growth of perennial ryegrass swards and the N and
water fluxes of a paddock. The model operates with a daily
time step simulating soil N mineralization, immobilization and
water balance, grass growth, grass N uptake and grass N content.
The model operates through a daily potential growth which is
dependent on solar radiation and the total green biomass. To cal-
culate the actual daily growth, this potential growth is then multi-
plied by parameters depending on environmental conditions
(temperature, water in the soil and radiation) and a parameter
depending on the soil mineral N availability compared to the N
demand associated with potential grass growth. The MoSt
model has been validated against various historical data (Ruelle
et al., 2018a, 2018b; Ruelle and Hennessy, 2019) and is currently
being used to predict grass growth weekly on 78 commercial
farms in Ireland (Ruelle and Hennessy, 2019; Ruelle et al., 2021).

Model simulations

Currently, Ireland has a derogation from the Nitrates Directive (91/676/EEC) allowing a maximal annual organic N loading of 250 kg N/ha (equivalent to an SR of 2.75 cows/ha) in conjunction with a maximal chemical fertilizer application of 250 kg N/ha/year for the farms in derogation. This paper is focused on grazing dairy farm with SRs between 2.0 and 3.0 cows/ha. The descriptions of the simulations are summarized in Table 1. Farm size was fixed at 40 ha and included 40 paddocks. A total of 140 main scenarios were analysed:

(1) Stocking rate – five stocking rates – 2, 2.25, 2.5, 2.75 and 3 cows/ha.

(2) Fertilizer N application rate – seven rates – 0 to 300 kg N/ha/year increasing by a step of 50.

(3) Soil type – two soil types – a free-draining soil (FDS: 6% organic matter (OM), proportion of 0.60 sand and 0.15 clay, soil depth of 100 cm) corresponding to a Moorepark soil type and a heavy soil (HS: 8% OM, proportion of 0.28 sand and 0.36 clay, soil depth of 100 cm) corresponding to a Ballyhaise soil type, those two soil types are the extremes of the soil types found on intensive Irish dairy farms.

(4) Agroclimate location – two weather locations were used – Ballyhaise (Bhy: 2003–18, 54°03′01.3″ N, 7°18′40.3″ W) corresponding to the north of the country and Moorepark (Mpk: 2003–18, 52°09′52.3″ N, 8°15′36.6″ W) corresponding to the south of the country with data being provided by Met Eireann (the Irish Meteorological Service; www.met.ie) synoptic weather stations on the farm sites. The weather data are presented in Table 2.

For each scenario, simulations were completed for 16 years on a
continuous cycle, meaning that some events occurring in year 1
could have consequences in year 2 and so on.

Annual *feeding* management

The feeding system used is the INRA feeding system (INRA, 2010)
where the quality of feed is characterized by its energy value (UFL
<table><caption>**Table 1.** Description of input parameters of the simulations to the PBHDM model</caption><tr><td>Farm size</td><td>40 ha divided into 40 paddocks of 1 ha</td></tr><tr><td>Yearly average grass quality</td><td>FV of 0.93 UEL/kg DM, 0.98 UFL/kg DM and PDI of 100 g/kg DM</td></tr><tr><td>Silage supplementation</td><td>6 kg DM/day/cow maximum (quality FV of 1.1 UEL/kg DM, 0.80 UFL/kg DM, PDI of 75 g/kg DM)</td></tr><tr><td>Concentrate quality</td><td>1.03 UFL/kg DM, PDI of 120 g/kg DM</td></tr><tr><td>Grazing season</td><td>10 February to 20 November</td></tr><tr><td>Grazing severity</td><td>'Very severe' before 15 of April<sup>a</sup><br>'Severe' after the 15 of April<sup>a</sup></td></tr><tr><td>Concentrate fed</td><td>4 kg DM/cow/day of concentrate for the first 3 months of lactation - total of 600 kg fresh/cow/lactation</td></tr><tr><td>Slurry application</td><td>Available slurry is 30 kg N/cow/year, 3/5 is organic 2/5 is mineral<br>Proportion of 0.22 applied in January, 0.22 applied in February, 0.22 applied in March and 0.33 applied in April</td></tr><tr><td>Nitrogen fertilization (kg/ha)</td><td>From January to September<sup>b</sup> from 0 kg N/ha/year to 300 kg N/ha/year in steps of 50 kg N/ha</td></tr><tr><td>SR (cows/ha)</td><td>2 2.25 2.5 2.75 3</td></tr><tr><td>Average number of lactating animals</td><td>80 90 100 110 120</td></tr><tr><td>Agroclimate location</td><td>Ballyhaise (Bhy)<br>GPS: 54°03'01.3" N, 7°18'40.3" W<br>Moorepark (Mpk)<br>GPS: 52°09'52.3" N, 8°15'36.6" W<br>16 years (2003-2018)</td></tr></table>

UFL, 'Unité Fourragère Lait' 1 UFL = 7.1 MJ of net energy for lactation (NEL); PDI, 'Protéine Digestible dans l'intestin' protein unit, UEL, 'Unité d'Encombrement Lait', Fill value unit (Faverdin et al., 2011; INRA, 2010).

<sup>a</sup>Calculation of the objective post grazing sward height (obj_postGH) depending on the pre grazing sward height (preGH) and the severity of the grazing:

Severe: obj_postGH = 4 - 0.1 × preGH + 0.015 × preGH².

Very severe: obj_postGH = 3.5 - 0.1 × preGH + 0.015 × preGH².

<sup>b</sup>Proportion of 0.06 applied in January, 0.05 applied in February, 0.17 applied in March, 0.21 applied in April, 0.12 applied in May, 0.14 applied in June, 0.09 applied in July, 0.09 applied in August and 0.07 applied in September based on Teagasc recommendations.

'Unité Fourragère Lait'; 1 UFL is the equivalent of 7.1 MJ of net energy for lactation, corresponding to the net energy for maintenance and lactation available from 1 kg of air-dried rolled barley), protein (PDI 'protéine digestible dans l'intestin') and fill value (FV expressed in UEL 'Unité d'Encombrement Lait'; 1 UEL is the equivalent to an intake of 140 g DM/kg of metabolic body-weight (BW<sup>0.75</sup>)). Concentrate does not have a fixed FV as it depends on the substitution rate between concentrate and silage. In the simulation, concentrate supplementation (1.03 UFL/kg DM, PDI of 120 g/kg DM) was fed at the rate of 4 kg DM/day/cow for the first 5 months of lactation irrespective of the grass available on the farm resulting in a concentrate feeding level of 600 kg fresh/cow/lactation following Teagasc recommendations. No additional concentrate has been fed during the rest of the lactation. The grazing season was from 10 February to 20 November, irrespective of the amount of grass on the farm, and only lactating cows were allowed to graze. While indoors, lactating cows were

fed grass silage *ad libitum* (quality FV of 1.1 UEL/kg DM, 0.80 UFL/kg DM, PDI of 75 g/kg DM), dry cows were offered 80% of the *ad libitum* intake of the lactating cows (around 10 kg of DM of silage/cow/day). During the grazing season the only reason the herd could return indoors was if the paddocks were deemed ungrazable due to waterlogged soil (i.e. if the soil water content was greater than 90% of saturation level), in that case cows were fed *ad libitum* grass silage indoors.

As the model does not predict the grass nutritive value, the grass quality was constant throughout the simulations at 0.98 UFL/kg DM, PDI of 100 g/kg DM and FV of 0.93 UEL/kg DM. As the model describes the reduced intake during the defoliation process, the value of 0.93 UEL/kg DM corresponds to the FV at the start of grazing in a paddock and reduces as the grass height in the paddock declines (Ruelle *et al.*, 2015).

## Grazing management and supplementation

The lactating herd moves from one paddock to another based on an objective post-grazing sward height which is dependent on the pre-grazing sward height (Table 1). The farm grass cover (i.e. quantity of grass available on the farm; kg DM/ha) is evaluated daily and is compared with the herd requirement. In situations where there is an excess farm grass cover on the overall farm, some paddocks are allocated for silage conservation rather than grazing and vice versa. In the case of a grass deficit, grass silage can be fed up to 6 kg DM/cow. In the simulations, priority is always given to grazing over silage; it is assumed that silage would always be available to purchase onto the farm.

## Fertilizer management

The different annual levels of N fertilizer application simulated were 0–300 kg N/ha with steps of 50 kg N/ha. For each simulation the mineral N fertilizer was applied from 15 January to 15 September with the repartition of application as follows: proportion of 0.06 in January, 0.05 in February, 0.17 in March, 0.21 in April, 0.12 in May, 0.14 in June, 0.09 in July, 0.09 in August and 0.07 in September based on pre-2022 Teagasc recommendations (Wall and Plunkett, 2021). As an example, the N fertilizer application rates for the months of January, February, March, April, May, June, July, August and September used to simulated 100 and 300 kg N/ha were 6 and 18 kg N/ha, 5 and 15 kg N/ha, 17 and 51 kg N/ha, 21 and 63 kg N/ha, 12 and 36 kg N/ha, 14 and 42 kg N/ha, 9 and 27 kg N/ha, 9 and 27 kg N/ha, 7 and 28 kg N/ha, respectively. Slurry was applied to the farm from January to April as follows: proportion of 0.22 in January, 0.22 in February, 0.22 in March and 0.33 in April. The quantity of slurry produced was 30 kg/cow/year (based on 20 kg collected during the housing period and 10 kg during milking throughout the year), assuming 2/5 was mineral N and 3/5 was organic N.

## Model outputs

The main outputs from the simulations are grass growth (kg DM/ha), N response (kg DM/kg N applied), number of grazing days, number of days at grazing without silage supplementation, grass intake (kg/cow and kg/ha), silage fed (at grazing, indoor for lactating and indoor for dry cow in kg DM/cow and kg DM/ha), milk, protein and fat produced (kg/cow and per ha), amount of silage produced (kg/ha) and yearly surplus or deficit of silage (kg/ha). All outputs are extracted per day and are then summarized by week, season and year for the full simulation. Annual farm gate N surplus was calculated for all years for each simulation. The calculation included all purchased N and
<table><caption>**Table 2.** Monthly description of the meteorological data at the two agroclimate locations (sa)</caption><thead><tr><th></th><th>Location</th><th>Jan</th><th>Feb</th><th>March</th><th>April</th><th>May</th><th>June</th><th>Jul</th><th>Aug</th><th>Sept</th><th>Oct</th><th>Nov</th><th>Dec</th><th>Total</th></tr></thead><tbody><tr><td rowspan="2">Mean temperature (°C)</td><td>Bhy</td><td>4.7 (1.7)</td><td>4.7 (1.1)</td><td>5.9 (1.4)</td><td>8.2 (1.2)</td><td>10.7 (0.8)</td><td>13.8 (0.9)</td><td>15.0 (1.2)</td><td>14.5 (0.7)</td><td>12.7 (0.8)</td><td>10.1 (1.1)</td><td>6.6 (1.3)</td><td>5.0 (2.0)</td><td>9.3 (0.4)</td></tr><tr><td>Mpk</td><td>5.6 (1.3)</td><td>5.4 (1.1)</td><td>6.6 (1.2)</td><td>8.8 (1.0)</td><td>11.4 (0.7)</td><td>14.3 (0.9)</td><td>15.6 (1.1)</td><td>15.0 (0.8)</td><td>13.3 (0.8)</td><td>10.6 (1.0)</td><td>7.5 (1.3)</td><td>6.1 (1.8)</td><td>10.0 (0.4)</td></tr><tr><td rowspan="2">Cumulative rainfall (mm)</td><td>Bhy</td><td>100 (43)</td><td>75 (36)</td><td>73 (26)</td><td>60 (24)</td><td>71 (30)</td><td>64 (38)</td><td>86 (32)</td><td>95 (34)</td><td>80 (43)</td><td>92 (36)</td><td>98 (45)</td><td>111 (57)</td><td>1006 (116)</td></tr><tr><td>Mpk</td><td>107 (39)</td><td>78 (56)</td><td>79 (28)</td><td>67 (39)</td><td>71 (23)</td><td>76 (47)</td><td>80 (35)</td><td>73 (51)</td><td>74 (33)</td><td>97 (43)</td><td>107 (58)</td><td>105 (66)</td><td>1014 (132)</td></tr><tr><td rowspan="2">Mean solar radiation (J/cm²)</td><td>Bhy</td><td>202 (24)</td><td>422 (39)</td><td>778 (97)</td><td>1281 (139)</td><td>1588 (126)</td><td>1644 (187)</td><td>1531 (165)</td><td>1236 (113)</td><td>916 (59)</td><td>529 (60)</td><td>275 (33)</td><td>159 (20)</td><td>882 (33)</td></tr><tr><td>Mpk</td><td>250 (30)</td><td>459 (55)</td><td>821 (97)</td><td>1292 (137)</td><td>1555 (129)</td><td>1715 (218)</td><td>1515 (188)</td><td>1315 (110)</td><td>956 (86)</td><td>550 (59)</td><td>314 (34)</td><td>189 (23)</td><td>911 (46)</td></tr></tbody></table>

N sales, but excluded net mineralization and atmospheric deposition. The equation for the annual N surplus (kg N/ha) was *Nin* – *Nout* where:

$$Nin = Nconcentrate + fertilization + Ninsilage bought$$

$$Nconcentrate = concentrate fed per ha (kg DM/ha/year)$$

$$\times \frac{18.1}{6.25 \times 100}$$

With 6.25 as the kjeldahl conversion factor for feed and 18.1 the average protein percentage of concentrate feed (kg DM, DAFM, private report):

$$NinSilage = Silage bought per ha (kg DM/ha/year)$$

$$\times \frac{13.5}{6.25 \times 100}$$

With 6.25 as the kjeldahl conversion factor for feed and with 13.5 the average protein content of silage feed (kg DM, DAFM, private report):

*fertilization* was the chemical N fertilizer applied (kg N/ha/year):

$$Nout = Nmilk + meatSold + Ninsilagesold$$

$$\text{where } Nmilk = milkprotein sold (kg/ha/year)/6.38$$

with 6.38 as the kjeldahl conversion factor for milk.

Cow weight gain was included at 35 kg/year to calculate an annualized net increase in liveweight from when a cow enters the herd as a replacement heifer and leaves as a cull cow. The liveweight of a calf when the calf leaves the farm was assumed to be 40 kg. It was assumed that calves had a N content of 0.029 kg N/kg body weight and the cow liveweight gain had a N content of 0.024 kg N/kg body weight (ARC, 1994) resulting in the following equation:

$$Nmeat = (40 \times 0.029 + 35 \times 0.024) \times SR$$

$$NinSilage = Silage sold per ha (kg DM/ha/year) \times \frac{13.5}{6.25 \times 100}$$

with 6.25 as the kjeldahl conversion factor for feed and with 13.5 the average protein content of silage feed (kg DM, DAFM, private report).

## Results

### Grass dry matter yield

The average annual grass DM yield across scenarios and years was 12 424 kg DM/ha, with a minimum grass yield of 4971 kg DM/ha (SR 2 cows/ha, 0 kg/ha/year fertilizer, FDS, Mpk, 2018) and a maximum grass yield of 17 452 kg DM/ha (SR 3 cows/ha, 300 kg/ha/year fertilizer, HS, Mpk, 2011) (Table 3).

On average, across all years, at 250 kg/ha/year of chemical N and SR 2.5 cows/ha, the average annual grass yield was of 13 816 kg DM/ha for MpkFDS, 13 841 kg DM/ha for BhyFDS, 14 978 kg DM/ha for MpkHS and 14 441 for BhyHS. While the average yield for the two FDS types was close, variation was higher for the MpkFDS: ranging from 9420 kg DM/ha in 2018
<table><caption>**Table 3.** Impact of soil type, agroclimate location and N application rate (kg N/ha) on seasonal grass yield kg DM/ha (SD) and N response (kg DM/kg N)</caption><thead><tr><th>Soil type</th><th>Location</th><th>N level</th><th>Spring<sup>a</sup></th><th>Summer</th><th>Autumn</th><th>Winter</th><th>Annual grass growth</th><th>N response</th></tr></thead><tbody><tr><td rowspan="10">Free-draining</td><td rowspan="5">Bhy</td><td>0</td><td>1396 (383)</td><td>4473 (411)</td><td>2632 (246)</td><td>374 (81)</td><td>8875 (818)</td><td></td></tr><tr><td>50</td><td>1635 (409)</td><td>5141 (400)</td><td>2876 (255)</td><td>406 (86)</td><td>10 059 (854)</td><td>23.7</td></tr><tr><td>100</td><td>1844 (439)</td><td>5691 (403)</td><td>3141 (243)</td><td>431 (90)</td><td>11 108 (867)</td><td>21.0</td></tr><tr><td>150</td><td>2043 (468)</td><td>6194 (426)</td><td>3400 (250)</td><td>459 (95)</td><td>12 095 (879)</td><td>19.7</td></tr><tr><td>200</td><td>2224 (501)</td><td>6664 (447)</td><td>3617 (259)</td><td>483 (98)</td><td>12 988 (899)</td><td>17.9</td></tr><tr><td rowspan="5">Mpk</td><td>250</td><td>2386 (522)</td><td>7101 (464)</td><td>3814 (269)</td><td>510 (102)</td><td>13 810 (921)</td><td>16.5</td></tr><tr><td>300</td><td>2529 (548)</td><td>7513 (484)</td><td>4001 (274)</td><td>531 (105)</td><td>14 575 (956)</td><td>15.3</td></tr><tr><td>0</td><td>1611 (470)</td><td>4452 (517)</td><td>2414 (774)</td><td>442 (94)</td><td>8918 (1384)</td><td></td></tr><tr><td>50</td><td>1865 (503)</td><td>5072 (526)</td><td>2639 (866)</td><td>475 (95)</td><td>10 051 (1510)</td><td>22.7</td></tr><tr><td>100</td><td>2098 (531)</td><td>5612 (537)</td><td>2881 (925)</td><td>511 (100)</td><td>11 102 (1581)</td><td>21.0</td></tr><tr><td rowspan="5">Heavy</td><td rowspan="5">Bhy</td><td>150</td><td>2323 (562)</td><td>6110 (564)</td><td>3102 (965)</td><td>546 (104)</td><td>12 081 (1641)</td><td>19.6</td></tr><tr><td>200</td><td>2524 (583)</td><td>6571 (610)</td><td>3313 (983)</td><td>583 (113)</td><td>12 990 (1687)</td><td>18.2</td></tr><tr><td>250</td><td>2700 (603)</td><td>7002 (637)</td><td>3491 (1026)</td><td>611 (119)</td><td>13 803 (1763)</td><td>16.3</td></tr><tr><td>300</td><td>2853 (627)</td><td>7380 (668)</td><td>3630 (1087)</td><td>630 (121)</td><td>14 494 (1866)</td><td>13.8</td></tr><tr><td>0</td><td>1645 (455)</td><td>5082 (490)</td><td>2782 (369)</td><td>187 (142)</td><td>9695 (894)</td><td></td></tr><tr><td rowspan="10">Mpk</td><td>50</td><td>1876 (480)</td><td>5707 (500)</td><td>3068 (386)</td><td>209 (145)</td><td>10 860 (880)</td><td>23.3</td></tr><tr><td>100</td><td>2085 (525)</td><td>6277 (528)</td><td>3296 (402)</td><td>233 (149)</td><td>11 892 (886)</td><td>20.6</td></tr><tr><td>150</td><td>2277 (563)</td><td>6796 (562)</td><td>3498 (413)</td><td>253 (151)</td><td>12 823 (897)</td><td>18.6</td></tr><tr><td>200</td><td>2431 (598)</td><td>7262 (576)</td><td>3705 (431)</td><td>272 (154)</td><td>13 670 (915)</td><td>16.9</td></tr><tr><td>250</td><td>2550 (618)</td><td>7685 (601)</td><td>3921 (446)</td><td>290 (159)</td><td>14 446 (952)</td><td>15.5</td></tr><tr><td>300</td><td>2651 (634)</td><td>8049 (630)</td><td>4133 (457)</td><td>307 (161)</td><td>15 141 (973)</td><td>13.9</td></tr><tr><td>0</td><td>1851 (607)</td><td>5214 (528)</td><td>2888 (405)</td><td>286 (180)</td><td>10 238 (835)</td><td></td></tr><tr><td>50</td><td>2083 (645)</td><td>5816 (559)</td><td>3156 (407)</td><td>313 (182)</td><td>11 368 (837)</td><td>22.6</td></tr><tr><td>100</td><td>2301 (695)</td><td>6371 (601)</td><td>3388 (416)</td><td>340 (187)</td><td>12 401 (840)</td><td>20.7</td></tr><tr><td>150</td><td>2497 (724)</td><td>6885 (639)</td><td>3602 (426)</td><td>365 (191)</td><td>13 348 (869)</td><td>18.9</td></tr><tr><td>200</td><td>2654 (752)</td><td>7343 (670)</td><td>3811 (440)</td><td>386 (193)</td><td>14 194 (898)</td><td>16.9</td></tr><tr><td>250</td><td>2784 (768)</td><td>7757 (685)</td><td>4030 (454)</td><td>408 (198)</td><td>14 979 (924)</td><td>15.7</td></tr><tr><td>300</td><td>2895 (782)</td><td>8128 (707)</td><td>4234 (464)</td><td>429 (203)</td><td>15 686 (954)</td><td>14.1</td></tr></tbody></table>

<footnote><sup>a</sup>Spring: 1 February to 30 April; summer: 1 May to 31 July; autumn: 1 August to 31 October.</footnote>
<figure><img src="image_4.png" /><figcaption>**Fig. 1.** Representation of the variability of the mean grass growth by week of the year. In black average grass growth across the year, light grey grass growth in the year of 2011 (good growing year), dark grey average grass growth in 2018 (drought year). Mpk, Moorepark; Bhy, Ballyhaise; FDS, free-draining soil; HS, heavy soil.</figcaption></figure>

to 15 573 kg DM/ha in 2011 (±1787) compared to 12 479 kg DM/ha in 2018 and 155 521 kg DM/ha in 2017 (±926) for BhyFDS. This higher variation with the MpkFDS was associated with drought years of 2006, 2013 and 2018, coupled with poor autumn growth (average 1660 kg DM/ha), while for those years, at the Bhy location there was significant autumn growth (average 3772 kg DM/ha).

For the HSs, the annual grass yield variation was similar for the two weather patterns and there was higher yield for MpkHS ranging from 13 645 kg DM/ha in 2012 to 16 500 kg DM/ha (±945) in 2007 and from 12 604 to 16 351 kg DM/ha (±971) in BhyHS.

On average for each 0.25 point increase in SR, there was an increase in annual grass yield of 182 kg DM/ha. This increase in yield reduced as SR increased (193 kg DM/ha when SR increased from 2 to 2.25 cows/ha and 162 kg DM/ha when it increased from 2.75 to 3 cows/ha). There was a greater increase in annual grass yield with increased SR when no fertilizer was applied (average 275 kg DM/ha) compared to when fertilizer N was applied (108 kg DM/ha at 300 kg N/ha/year).

The growth during winter (1 November to 31 January) was lower under Bhy conditions (average 353 kg DM/ha) than Mpk conditions (average 452 kg N/ha) due to colder weather conditions. It was lower for the HS (average 307 kg DM/ha) than for the FDS (average 501 kg DM/ha) due to the HS being easily flooded during this period. In spring (1 February to 30 April) growth was lower under Bhy conditions (average 2112 kg DM/ha) than Mpk conditions (average 2360 kg DM/ha). The HS had a higher grass growth during the spring than the FDS (average of 2331 and 2150 kg DM/ha, respectively). In summer (1 May to 31 July), on average, growth was very similar for the two agroclimate locations, however it was higher for the HS (average 6742 kg DM/ha) than for the FDS (average 6070 kg DM/ha) due to being less prone to soil moisture deficit. In autumn (1 August to 31 October) growth was more variable among soil types and weather conditions. On average, growth was highest for MpkHS (3587 kg DM/ha) followed closely by BhyHS (3486 kg DM/ha). It was the opposite for the FDS with the MpkFDS (3071 kg DM/ha) being lower than the BhyFDS (3354 kg DM/ha). This is due to a warmer weather in the autumn

at Mpk, better for grass growth, but also linked to a lower rainfall, resulting in a drought in August leading to a growth decrease on the MpkFDS. In autumn, the MpkFDS had higher growth variability with a standard deviation (SD) of 1035 kg DM/ha compared to 528 kg DM/ha for the BhyFDS and around 611 kg DM/ha for the HS simulations.

While, on average across all the years, weather and annual grass yield at both locations are similar (Table 2), Fig. 1 highlights the possible differences in growth pattern for the same year at the two locations. Figure 1 presents the weekly average grass growth for all years at the two locations and for the two soil types, which highlights two extreme years of 2011 and 2018. Figure 1 highlights the higher water retention by the HS but also the difference in the weather events between Bhy and Mpk with a greater drought at Mpk than at Bhy.

Annual grass yield, number of days without silage supplementation, grass intake per cow, silage fed per cow, grass intake per ha and grass surplus or deficit per ha for all SRs are presented in Tables 4 to 6 for the annual fertilizer application rates 250, 200 and 150 kg N/ha/year, respectively.

## Cow DMI

Across all simulations, the average DMI per cow annually was 3050 kg DM grazed grass and 1649 kg DM grass silage, with the average amount of concentrate offered fixed at 600 kg. The corresponding production level was of 5456 kg milk, 246 kg fat and 205 kg protein per cow. The maximum grass DMI per cow was 3864 kg DM for MpkFDS at an SR of 2 cows/ha and 300 kg N/ha and the minimum grass DMI per cow was 1774 kg DM for the BhyHS in 2012 and 0 kg N/ha. While 2012 did not correspond to one of the worst years in terms of annual grass yield, the high amount of rainfall led to paddocks being ungrazable for a significant proportion of the grazing season on the HS farm. On average across all of the years, for the SR 2.5 cows/ha and the fertilization rate 250 kg N/ha/year (Table 4), the grass DMI/cow was 3502 kg DM (±130), 2761 kg DM (±382), 3496 kg DM (±308) and 2950 kg DM (±338) for the BhyFDS, BhyHS, MpkFDS and MpkHS
<table><caption>**Table 4.** Impact of soil type, agroclimate location and SR on mean annual grass yield (kg DM/ha), grazing indicators and dairy cow DM intake (kg DM/cow) at an annual fertilizer application rate of 250 kg N/ha</caption><thead><tr><th>Soil type</th><th>Location</th><th>SR (cows/ha)</th><th>Grass yield (kg DM/ha)</th><th>Full grazing days (no silage)</th><th>Grass intake (kg DM/cow)</th><th>Silage fed (kg DM/cow)</th><th>Grass intake (kg DM/ha)</th><th>Surplus/deficit (kg DM/ha)</th><th>N surplus (kg N/ha)</th></tr></thead><tbody><tr><td rowspan="12">Free-draining</td><td rowspan="4">Bhy</td><td>2</td><td>13 499 (914)</td><td>230 (27)</td><td>3607 (163)</td><td>1149 (163)</td><td>7213 (327)</td><td>1724 (738)</td><td>174 (16)</td></tr><tr><td>2.25</td><td>13 648 (893)</td><td>222 (24)</td><td>3570 (134)</td><td>1191 (144)</td><td>8032 (301)</td><td>870 (753)</td><td>188 (16)</td></tr><tr><td>2.5</td><td>13 841 (926)</td><td>208 (23)</td><td>3502 (130)</td><td>1264 (134)</td><td>8755 (325)</td><td>20 (726)</td><td>202 (16)</td></tr><tr><td>2.75</td><td>13 958 (940)</td><td>197 (24)</td><td>3404 (146)</td><td>1327 (147)</td><td>9361 (402)</td><td>-816 (842)</td><td>216 (19)</td></tr><tr><td rowspan="4">Mpk</td><td>3</td><td>14 106 (927)</td><td>189 (32)</td><td>3327 (181)</td><td>1373 (193)</td><td>9981 (543)</td><td>-1609 (883)</td><td>229 (19)</td></tr><tr><td>2</td><td>13 488 (1791)</td><td>229 (43)</td><td>3592 (244)</td><td>1160 (245)</td><td>7184 (487)</td><td>1732 (1300)</td><td>174 (28)</td></tr><tr><td>2.25</td><td>13 651 (1794)</td><td>218 (51)</td><td>3541 (310)</td><td>1214 (306)</td><td>7968 (697)</td><td>886 (1354)</td><td>187 (29)</td></tr><tr><td>2.5</td><td>13 816 (1787)</td><td>210 (50)</td><td>3496 (308)</td><td>1268 (300)</td><td>8740 (770)</td><td>15 (1355)</td><td>202 (29)</td></tr><tr><td rowspan="4">Heavy</td><td>2.75</td><td>13 948 (1813)</td><td>196 (51)</td><td>3367 (307)</td><td>1359 (299)</td><td>9260 (844)</td><td>-820 (1354)</td><td>216 (30)</td></tr><tr><td>3</td><td>14 114 (1793)</td><td>184 (54)</td><td>3282 (338)</td><td>1408 (323)</td><td>9846 (1014)</td><td>-1609 (1436)</td><td>229 (32)</td></tr><tr><td rowspan="4">Bhy</td><td>2</td><td>14 251 (945)</td><td>192 (26)</td><td>2768 (354)</td><td>1900 (344)</td><td>5536 (707)</td><td>1918 (772)</td><td>172 (17)</td></tr><tr><td>2.25</td><td>14 351 (966)</td><td>188 (30)</td><td>2752 (375)</td><td>1919 (357)</td><td>6192 (843)</td><td>1013 (785)</td><td>187 (17)</td></tr><tr><td>2.5</td><td>14 441 (971)</td><td>188 (28)</td><td>2761 (382)</td><td>1917 (354)</td><td>6903 (955)</td><td>104 (835)</td><td>202 (19)</td></tr><tr><td>2.75</td><td>14 535 (970)</td><td>183 (28)</td><td>2718 (363)</td><td>1938 (347)</td><td>7476 (997)</td><td>-769 (920)</td><td>216 (20)</td></tr><tr><td rowspan="4">Mpk</td><td>3</td><td>14 651 (986)</td><td>184 (30)</td><td>2693 (369)</td><td>1940 (348)</td><td>8079 (1106)</td><td>-1597 (904)</td><td>230 (20)</td></tr><tr><td>2</td><td>14 766 (909)</td><td>206 (23)</td><td>2986 (337)</td><td>1709 (314)</td><td>5973 (673)</td><td>2342 (667)</td><td>162 (14)</td></tr><tr><td>2.25</td><td>14 871 (925)</td><td>201 (23)</td><td>2963 (337)</td><td>1730 (312)</td><td>6666 (759)</td><td>1446 (742)</td><td>177 (16)</td></tr><tr><td>2.5</td><td>14 978 (945)</td><td>197 (25)</td><td>2950 (338)</td><td>1751 (308)</td><td>7376 (846)</td><td>543 (785)</td><td>192 (17)</td></tr><tr><td rowspan="2">Heavy</td><td>2.75</td><td>15 081 (947)</td><td>192 (20)</td><td>2896 (305)</td><td>1777 (282)</td><td>7964 (840)</td><td>-300 (808)</td><td>206 (17)</td></tr><tr><td>3</td><td>15 201 (946)</td><td>190 (19)</td><td>2855 (295)</td><td>1791 (276)</td><td>8564 (885)</td><td>-1129 (809)</td><td>220 (17)</td></tr></tbody></table>

Cows were fed with 600 kg of concentrate fresh per year and the only available supplementation available at grazing was grass silage.
<table><caption>**Table 5.** Impact of soil type, agroclimate location and SR on mean annual grass yield (kg DM/ha), grazing indicators and dairy cow DM intake (kg DM/cow) at an annual fertilizer application rate of 200 kg N/ha</caption><thead><tr><th>Soil type</th><th>Location</th><th>SR (cows/ha)</th><th>Grass yield (kg DM/ha)</th><th>Full grazing days (no silage)</th><th>Grass intake (kg DM/cow)</th><th>Silage fed (kg DM/cow)</th><th>Grass intake (kg DM/ha)</th><th>Surplus/deficit (kg DM/ha)</th><th>N surplus (kg N/ha)</th></tr></thead><tbody><tr><td rowspan="12">Free-draining</td><td rowspan="4">Bhy</td><td>2</td><td>12 630 (867)</td><td>222 (26)</td><td>3559 (146)</td><td>1188 (147)</td><td>7119 (293)</td><td>1139 (733)</td><td>137 (16)</td></tr><tr><td>2.25</td><td>12 802 (878)</td><td>208 (30)</td><td>3479 (166)</td><td>1269 (166)</td><td>7827 (374)</td><td>278 (780)</td><td>151 (17)</td></tr><tr><td>2.5</td><td>13 007 (888)</td><td>197 (25)</td><td>3429 (139)</td><td>1330 (146)</td><td>8572 (347)</td><td>-582 (775)</td><td>165 (17)</td></tr><tr><td>2.75</td><td>13 170 (892)</td><td>191 (30)</td><td>3370 (165)</td><td>1367 (171)</td><td>9267 (454)</td><td>-1405 (826)</td><td>179 (18)</td></tr><tr><td rowspan="4">Mpk</td><td>3</td><td>13 330 (913)</td><td>182 (30)</td><td>3260 (187)</td><td>1424 (173)</td><td>9780 (562)</td><td>-2160 (916)</td><td>191 (21)</td></tr><tr><td>2</td><td>12 613 (1687)</td><td>223 (48)</td><td>3557 (267)</td><td>1190 (276)</td><td>7114 (534)</td><td>1137 (1220)</td><td>137 (26)</td></tr><tr><td>2.25</td><td>12 797 (1680)</td><td>207 (51)</td><td>3471 (307)</td><td>1282 (309)</td><td>7811 (691)</td><td>272 (1252)</td><td>151 (27)</td></tr><tr><td>2.5</td><td>12 985 (1705)</td><td>199 (52)</td><td>3427 (320)</td><td>1332 (314)</td><td>8567 (800)</td><td>-576 (1264)</td><td>165 (27)</td></tr><tr><td rowspan="4">Heavy</td><td>2.75</td><td>13 189 (1716)</td><td>187 (53)</td><td>3341 (339)</td><td>1392 (326)</td><td>9187 (933)</td><td>-1369 (1349)</td><td>178 (30)</td></tr><tr><td>3</td><td>13 366 (1754)</td><td>180 (58)</td><td>3236 (387)</td><td>1443 (349)</td><td>9709 (1160)</td><td>-2119 (1413)</td><td>190 (32)</td></tr><tr><td rowspan="4">Bhy</td><td>2</td><td>13 441 (911)</td><td>192 (29)</td><td>2757 (380)</td><td>1905 (356)</td><td>5514 (759)</td><td>1385 (808)</td><td>133 (18)</td></tr><tr><td>2.25</td><td>13 555 (925)</td><td>189 (30)</td><td>2755 (384)</td><td>1917 (360)</td><td>6200 (864)</td><td>484 (913)</td><td>148 (20)</td></tr><tr><td>2.5</td><td>13 678 (942)</td><td>185 (30)</td><td>2736 (385)</td><td>1935 (354)</td><td>6841 (962)</td><td>-410 (899)</td><td>163 (20)</td></tr><tr><td>2.75</td><td>13 821 (933)</td><td>180 (29)</td><td>2707 (373)</td><td>1951 (347)</td><td>7443 (1026)</td><td>-1266 (872)</td><td>177 (19)</td></tr><tr><td rowspan="4">Mpk</td><td>3</td><td>13 853 (925)</td><td>183 (32)</td><td>2681 (388)</td><td>1947 (376)</td><td>8043 (1164)</td><td>-2128 (982)</td><td>192 (22)</td></tr><tr><td>2</td><td>13 946 (882)</td><td>202 (19)</td><td>2959 (322)</td><td>1729 (294)</td><td>5918 (644)</td><td>1785 (673)</td><td>124 (15)</td></tr><tr><td>2.25</td><td>14 073 (893)</td><td>196 (20)</td><td>2919 (299)</td><td>1768 (282)</td><td>6567 (673)</td><td>892 (712)</td><td>139 (15)</td></tr><tr><td>2.5</td><td>14 202 (908)</td><td>193 (23)</td><td>2920 (343)</td><td>1776 (303)</td><td>7300 (857)</td><td>26 (749)</td><td>153 (16)</td></tr><tr><td rowspan="2">Heavy</td><td>2.75</td><td>14 308 (915)</td><td>191 (21)</td><td>2889 (314)</td><td>1783 (289)</td><td>7945 (864)</td><td>-834 (751)</td><td>167 (16)</td></tr><tr><td>3</td><td>14 441 (921)</td><td>185 (18)</td><td>2838 (283)</td><td>1810 (268)</td><td>8515 (848)</td><td>-1659 (846)</td><td>181 (18)</td></tr></tbody></table>

Cows were fed with 600 kg of concentrate fresh per year and the only available supplementation available at grazing was grass silage.
<table><caption>Table 6. Impact of soil type, agroclimate location and SR on mean annual grass yield (kg DM/ha), grazing indicators and dairy cow DM intake (kg DM/cow) at an annual fertilizer application rate of 150 kg N/ha</caption><thead><tr><th>Soil type</th><th>Location</th><th>SR (cows/ha)</th><th>Grass yield (kg DM/ha)</th><th>Full grazing days (no silage)</th><th>Grass intake (kg DM/cow)</th><th>Silage fed (kg DM/cow)</th><th>Grass intake (kg DM/ha)</th><th>Surplus/deficit (kg DM/ha)</th><th>N surplus (kg N/ha)</th></tr></thead><tbody><tr><td rowspan="12">Free-draining</td><td rowspan="4">Bhy</td><td>2</td><td>11 658 (841)</td><td>205 (32)</td><td>3486 (172)</td><td>1264 (183)</td><td>6971 (344)</td><td>434 (707)</td><td>102 (15)</td></tr><tr><td>2.25</td><td>11 878 (834)</td><td>192 (32)</td><td>3402 (185)</td><td>1346 (193)</td><td>7654 (415)</td><td>-382 (777)</td><td>116 (17)</td></tr><tr><td>2.5</td><td>12 128 (849)</td><td>188 (27)</td><td>3374 (141)</td><td>1381 (152)</td><td>8436 (354)</td><td>-1200 (679)</td><td>129 (15)</td></tr><tr><td>2.75</td><td>12 321 (837)</td><td>174 (32)</td><td>3268 (174)</td><td>1465 (188)</td><td>8988 (480)</td><td>-2037 (834)</td><td>142 (18)</td></tr><tr><td rowspan="4">Mpk</td><td>3</td><td>12 491 (881)</td><td>174 (27)</td><td>3230 (174)</td><td>1461 (160)</td><td>9689 (522)</td><td>-2759 (857)</td><td>155 (19)</td></tr><tr><td>2</td><td>11 656 (1608)</td><td>213 (47)</td><td>3496 (268)</td><td>1249 (271)</td><td>6991 (536)</td><td>468 (1107)</td><td>102 (24)</td></tr><tr><td>2.25</td><td>11 886 (1634)</td><td>196 (48)</td><td>3417 (303)</td><td>1341 (301)</td><td>7689 (682)</td><td>-380 (1133)</td><td>115 (25)</td></tr><tr><td>2.5</td><td>12 121 (1634)</td><td>183 (56)</td><td>3315 (366)</td><td>1427 (347)</td><td>8287 (914)</td><td>-1189 (1266)</td><td>128 (28)</td></tr><tr><td rowspan="4">Heavy</td><td>2</td><td>12 552 (888)</td><td>192 (30)</td><td>2759 (375)</td><td>1911 (359)</td><td>5518 (749)</td><td>775 (811)</td><td>96 (18)</td></tr><tr><td>2.25</td><td>12 688 (887)</td><td>186 (30)</td><td>2737 (375)</td><td>1934 (352)</td><td>6158 (843)</td><td>-116 (830)</td><td>111 (18)</td></tr><tr><td>2.5</td><td>12 836 (910)</td><td>185 (33)</td><td>2738 (410)</td><td>1942 (380)</td><td>6845 (1026)</td><td>-989 (882)</td><td>125 (20)</td></tr><tr><td>2.75</td><td>12 965 (916)</td><td>182 (29)</td><td>2711 (363)</td><td>1951 (343)</td><td>7456 (999)</td><td>-1838 (863)</td><td>139 (19)</td></tr><tr><td rowspan="4">Mpk</td><td>3</td><td>13 077 (906)</td><td>179 (27)</td><td>2648 (365)</td><td>1972 (345)</td><td>7944 (1096)</td><td>-2673 (874)</td><td>154 (20)</td></tr><tr><td>2</td><td>13 042 (848)</td><td>201 (21)</td><td>2943 (320)</td><td>1737 (305)</td><td>5886 (641)</td><td>1177 (724)</td><td>87 (16)</td></tr><tr><td>2.25</td><td>13 185 (847)</td><td>196 (24)</td><td>2930 (341)</td><td>1761 (312)</td><td>6593 (768)</td><td>295 (743)</td><td>102 (16)</td></tr><tr><td>2.5</td><td>13 347 (862)</td><td>192 (27)</td><td>2918 (347)</td><td>1775 (317)</td><td>7295 (868)</td><td>-565 (810)</td><td>116 (18)</td></tr><tr><td rowspan="2">Heavy</td><td>2.75</td><td>13 517 (875)</td><td>187 (21)</td><td>2870 (312)</td><td>1808 (290)</td><td>7891 (857)</td><td>-1393 (788)</td><td>129 (17)</td></tr><tr><td>3</td><td>13 650 (882)</td><td>184 (23)</td><td>2827 (304)</td><td>1823 (282)</td><td>8480 (913)</td><td>-2204 (786)</td><td>143 (17)</td></tr></tbody></table>

Cows were fed with 600 kg of concentrate fresh per year and the only available supplementation available at grazing was grass silage.
<figure><img src="image_5.png" /><figcaption>**Fig. 2.** Mean grass intake (kg DM) per ha and per cow depending on soil type, agroclimate location, SR and fertilizer level. Mpk, Moorepark; Bhy, Ballyhaise; FDS, free-draining soil; HS, heavy soil.</figcaption></figure>

systems, respectively (Fig. 2). The corresponding silage fed per cow per year was 1264 kg DM (SD 134 kg DM), 1917 kg DM (SD 354 kg DM), 1268 kg DM (SD 300 kg DM) and 1751 kg DM (SD 308 kg DM) for the BhyFDS, BhyHS, MpkFDS and MpkHS systems, respectively. The lower grass DMI and greater grass DMI variability per cow on the HS was due to paddocks being flooded more often leading to some paddocks not being suitable for grazing during the summer months.

At the SR 2.5 cows/ha, the increase in grass DMI/cow/kg of N applied was of 1.6, 0.5, 1.6 and 0.6 kg DM/kg N for the BhyFDS, BhyHS, MpkFDS and MpkHS, respectively (Fig. 2). These results show that an increase in N fertilizer application rate on an FDS leads to an increase in grass DMI/cow, while on HS the extra fertilizer was used to increase silage produced. Grass DMI/cow declined as SR increased, while grass intake per ha increased (Figs 2 and 3).

## Nitrogen response

The N response was calculated based on the previous level of N fertilizer (increasing in steps of 50 kg N/ha/year) (Table 3). The average N response due to increasing N fertilizer application by 50 kg/ha/year was 18.5 kg DM/ha, with a minimum of 9.4 kg DM/ha (SR 3 cows/ha and increasing fertilizer from 250 to 300 kg N/ha/year, BhyFDS in 2018) and a maximum of 28.0 kg N/ha (SR 2.25 cows/ha and increasing fertilizer from 0 to 50 kg N/ha/year, BhyFDS, 2014). The N response decreased as N fertilizer application rate increased, with an average response across all the simulations of 22.9 kg DM/kg N applied between the level 0 and 50 kg N/ha/year, and a response of 14.5 kg DM/kg N applied as application rate increased from 250 to 300 kg N/ha/year. On

average, across all years, the response to increasing fertilizer application rate from 200 to 250 kg N/ha/year at an SR of 2.5 cows/ha was of 16.6 (±2.0), 16.7 (±1.4), 15.5 (±1.5) and 15.3 (±1.4) for MpkFDS, BhyFDS, MpkHS and BhyHS, respectively.

Figure 4 presents the average N responses by season, fertilizer level, soil type and weather. On average, across all years, the lowest response was in the spring (average 9.2 kg DM/kg N applied from 200 to 250 kg N/ha/year) followed by the autumn (average 12.9 kg DM), while the highest response was always in the summer (average 23.9 kg DM).

## Effect of nitrogen fertilizer rate and SR on feed self-sufficiency

When fertilizer N application rate was reduced, the number of SRs that were forage self-sufficient declined (Fig. 3). When 250 kg N/ha/year were applied SR 2, 2.25 and 2.5 were forage self-sufficient, while the systems stocked at 2.75 and 3 cows/ha had to purchase in 26 920 and 59 280 kg DM grass silage, respectively, to balance energy requirements on the farm (Table 4). When N fertilizer was reduced to 200 kg N/ha/year, only the systems stocked at 2 and 2.25 cows/ha were forage self-sufficient, and the 2.5, 2.75 and 3 cows/ha systems had to purchase in 15 280, 48 600 and 80 480 kg DM of grass silage, respectively, to balance energy requirements on the farm (Table 5). When N fertilizer application rate was reduced to 150 kg N/ha/year, only the system stocked at 2.0 cows/ha was self-sufficient for forage, and the 2.25, 2.5, 2.75 and 3 cows/ha had to purchase in 5680, 39 280, 72 400 and 104 000 kg DM of grass silage, respectively, to balance feed requirements (Table 6). A linear regression analysis was undertaken to examine if there was a relationship between N applied, SR and
<figure><img src="image_6.png" /><figcaption>**Fig. 3.** On farm surplus or deficit depending on soil type, agroclimate location, SR and fertilizer level. Mpk, Moorepark; Bhy, Ballyhaise; FDS, free-draining soil HS, heavy soil.</figcaption></figure>

farm self-sufficiency. Other parameters included in the regression were the soil type, location and year. The results showed that a reduction in N fertilizer application rate of 50 kg/ha/year required a decrease in SR of 0.18 cows/ha to maintain self-sufficiency ($R^2 = 0.90$).

### Effect of nitrogen fertilizer rate and SR on farm gate annual N surplus

Across all simulations, the average annual N surplus was 127 kg of N/ha with a lowest value of -40 kg N/ha (MpkHS in 2011, 0 kg/ ha/year of N applied and SR 2.0 cows/ha) and a highest value of 341 kg N/ha (MpkFDS in 2018, 300 kg/ha/year N applied and SR 3.0 cows/ha). Overall the annual N surplus increased as annual N fertilizer application rate and SR increased. On average, across all simulations, the FDS had a higher annual N surplus (131 kg N/ha) than the HS (123 kg N/ha) and the Bly agroclimate location led to a slightly higher annual N surplus (130 kg N/ha) than the Mpk (125 kg N/ha).

### Discussion

Overall, the simulations conducted showed realistic agronomic output responses. The average response of 727 kg DM/ha per

unit increase in SR in this study aligns well with the previously published research conducted at Mpk between 2013 and 2016 where the average response was of 671 kg of DM/ha (Coffey et al., 2018). McCarthy et al. (2016) found a higher response of 1206 kg DM/ha. McCarthy et al. (2016) also showed that the marginal response decreased as SR increased. The average grass intake for FDS was 3499 kg DM at SR 2.5 cows/ha and N fertilizer application rate of 250 kg N/ha/year. This number compares well with the results of McCarthy et al. (2016) which had an average grass intake per cow of 3718 kg DM at an SR of 2.5 cows/ha.

In the simulations, a reduced grass DM response to increasing fertilizer N application was observed; an average response of 23.0 kg of DM/kg N applied was recorded between 0 and 50 kg N/ha/year, and a response of 15.0 kg of DM/kg N applied was recorded when N fertilizer application increased from 250 to 300 kg N/ha/year. Previous studies (Enriquez-Hidalgo et al., 2016) found similar results with a response ranging from 8.0 to 28.3 kg DM/kg N applied at an application level of 120 kg N/ha/year compared to 0 kg N/ha/year and Hennessy et al. (2008) found responses varying between 15.0 and 33.1 kg DM/kg N applied (when compared to 0 treatment). Previously published studies show a wide range of N responses to N fertilizer application in spring from 2.8 to 22.9 kg DM/kg N applied (O'Donovan et al.,
<figure><img src="image_7.png" /><figcaption>**Fig. 4.** Variation in nitrogen response (kg DM/kg N applied) by season depending on soil type, agroclimate location and nitrogen level. Mpk, Moorepark; Bhy, Ballyhaise; FDS, free-draining soil; HS, heavy soil.</figcaption></figure>

2004). In that study the N response decreased as N applied increased. The average response in early April was of 19.4 kg DM/kg N when 30 kg N/ha was applied which decreased to 11.6 kg DM/kg N when 90 kg N/ha was applied. In the O'Donovan *et al.* (2004) study, cut plots were used rather than grazed swards which could explain the higher response in this study. Some older studies such as MacCarthy (1982) showed lower responses in spring (5.2 kg DM/kg N applied), which was calculated in early March and not in April. The selection of perennial ryegrass cultivars for grazed grassland in Ireland has a strong focus on seasonal growth (spring and autumn), which probably increased the N responses in the past number of years (McEvoy *et al.*, 2011). The grass growth response to fertilizer N application is greatest when grass is actively growing and when conditions of light, temperature and water supply are optimum (Whitehead, 1995; Hennessy, 2005). Differences in these factors are largely responsible for the differences in herbage DM yield and the seasonal response to fertilizer N (Whitehead, 1995; Hennessy, 2005).

Due to the high variability in N response within season, depending on the weather conditions in a given year, developing an appropriate N fertilizer application strategy is challenging. For example, in the spring of 2018 the N response was low compared to the other years, at approximately 7.2 kg DM/kg N applied for both soil types and both sites at annual N fertilizer application rates of 200–250 kg N/ha. In that year, the N application strategy

could have been broadly the same in all regions and soil types due to the low temperatures, high rainfall and poor grass growth, and the best recommendation would have been to delay N application until March or early April. In contrast, in 2007, under the Mpk conditions, the average N response in spring at annual N application rates of 200–250 kg N/ha was approximately 16.1 kg DM/kg N applied for the FDS and 10.9 kg DM/kg N applied for the HS largely due to a milder spring. In that same year under Bhy conditions, the response was only 12.0 kg DM/kg N applied for the FDS and 4.8 kg DM/kg N applied for the HS, which shows the impact of having local and farm specific management for N fertilizer. These examples illustrate the use of precision management based on weather conditions and soil type to guide N fertilizer advice during the low growth periods.

In a farm, soil type and weather are the two parameters a farmer has little control over. Future published modelling work on grass growth should always take into account weather variability (location and year) as well as soil type as it has been shown in this study that a conclusion drawn from 2007 findings would be the opposite to that found in 2018. Precision grassland management techniques that restrict the application of chemical and organic fertilizer N applications when severe weather events are forecasted would also be vital in terms of reducing N loss. This study focused on a small geographical area with differences in soil type and climatic conditions in Ireland. Grazing systems globally will be dealing with much more diverse weather conditions
within a region as well as across soil types. Management needs to
move away from prescriptive structures with singular advice to a
situation where dynamic systems and guidelines interact with the
reality on the ground. Systems that capture data and allow
dynamic management similar to PastureBase Ireland (Hanrahan
et al., 2017) or HappyGrass (https://www.happygrass.fr/) should
be used in other regions and jurisdictions. This would facilitate
more tailored management within and across grass-based systems
of milk production.

This study has shown that there is an increase in N surplus as SR and nitrogen fertilizer application rate increases, which is in agreement with previously published studies (Ledgard et al., 1997; Gourley et al., 2012). Increased N surplus is often linked to an increase in N leaching (Stott and Gourley, 2016). There is an increased risk of N being leached on well-drained and moderately drained soils as surplus N increases (Scholefield et al., 1991, 1993; Fraters et al., 2015). However, this loss is not apparent in HSs (Scholefield et al., 1993), therefore further emphasizing the requirement for tailored management depending on soil type.

The introduction of the farm to fork strategy (2021) requires a reduction in chemical N fertilizer use and N loss. In this study, it was demonstrated that, on average, with each 50 kg/ha/year reduction in N fertilizer application rate, SR had to be reduced by 0.18 cows/ha in order to maintain feed self-sufficiency (Fig. 3). Grass-based systems are focused on maximizing grass production and utilization and minimizing the amount of feed imported to the farm (Shalloo and Hanrahan, 2020). This is more sustainable, both profitably and environmentally, than focusing on maximizing milk yield per cow which tends to move away from grass based systems. With a move to lower grass production there is a risk of greater feed importation to the farm which will lead to reduced profitability and a deterioration in environmental sustainability (O'Brien et al., 2015). In the simulations described, when all systems were compared regardless of soil type or agroclimate location, as SR increased the level of surplus feed produced in the system declined.

Longer-term strategies exist that could be used to maintain grass production with reduced N fertilizer application rate, for example achieving and maintaining optimum soil fertility while incorporating and maintaining clover in the swards (Dillon et al., 2020). The simulations described in the current study assume high soil fertility levels. Increasing soil fertility on farms with sub-optimal soil fertility could result in maintenance of the current level of grass production despite reduced N fertilizer application. Indeed, increasing soil fertility (pH, P and K) has been shown to increase N use efficiency (Wall, 2019). White clover (Trifolium repens L.) is the most important sown legume species in grazed grassland in temperate regions including Ireland (Peyraud et al., 2009). The incorporation of white clover in grazing swards takes advantage of their capacity to fix atmospheric N and make it available for plant growth allowing the reduction in chemical N fertilizer use (Egan et al., 2018; Enriquez-Hidalgo et al., 2018). The use of white clover (T. repens L.) in recent decades in Ireland has been low due to the availability of cheap N fertilizer (Gilliland et al., 2009). Incorporating white clover into perennial ryegrass swards can reduce N fertilizer requirement without affecting herbage production (Egan et al., 2018), while sometimes increasing herbage production at a given N fertilizer application rate (Enriquez-Hidalgo et al., 2018; McClearn et al., 2020a), increasing herbage quality (Ulyatt and Morley, 1981; Ayres et al., 2002) and increasing DM intake (Ribeiro Filho et al., 2003, 2005) and milk production (Ribeiro Filho et al.,

2003, 2005; McClearn et al., 2020a). All of these components
are associated with increases in farm profitability (McClearn
et al., 2020b).

Conclusion

This paper showed that increasing either SR or annual N fertilizer
application rate increased annual grass yield but also lead to an
increase in annual N surplus. While the average annual grass
yield looked similar between the two soil types the FDS had
more grazed grass by the animals while in the HSs, there was
more grass silage conserved and consumed as the ground condi-
tions were not always fit for grazing. Weather variability and its
impact on grass yield highlighted the importance of having flexi-
bility and dynamic grassland management. This paper clearly
illustrates the challenges in balancing feed supply with lower
N input at varying SRs. This study shows that irrespective of
the soil type, a 50 kg reduction in N fertilizer application per ha
would lead to a requirement to reduce SR by 0.18 cows/ha if no
additional measures are taken on farm to maintain herbage
production and feed supplementation is not increased.

**Author contributions.** E. R., M. O. D., P. D., L. S., L. D., B. H., M. E. and D. H. conceptualized the study. E. R., M. O. D., P. D., L. S. and L. D. designed the modelled scenarios. E. R. ran the simulation. D. H., B. H., M. E., L. S., M. O. D. and P. D. analysed the industry consequences of the model outputs. E. R. wrote the article. L. S., D. H. and P. D. critically revised the paper.

**Financial support.** The authors are grateful for the funding provided by the Irish Dairy Levy Trust (Dublin, Ireland).

**Conflict of interest.** The authors declare there are no conflicts of interest.

References

ARC (1994) *The Nutrient Requirements of Nutrient Livestock. Technical Review by an Agricultural Research Council Working Party.* Wallingford, UK: CAB International.

Ayres J, Nandra K and Turner A (2002) A study of the nutritive value of white clover (*Trifolium repens* L.) in relation to different stages of phenological maturity in the primary growth phase in spring. *Grass and Forage Science* **53**(3), 250–259. https://doi.org/10.1046/j.1365-2494.1998.00137.x.

Coffey E, Delaby L, Fleming C, Pierce K and Horan B (2018) Multi-year evaluation of stocking rate and animal genotype on milk production per hectare within intensive pasture-based production systems. *Journal of Dairy Science* **101**, 2448–2462.

Dillon P, Shalloo L, Ruelle E, Delaby L, Fenton O, Crosson P, Creighton P, O'Donovan M, Horan B and Wall D (2020) *Review of the Influence of Chemical Nitrogen Application Rate, Soil Type and Agroclimate Location on Grass Production, Feed Budgets, Nitrogen use Efficiency, Environmental Impact and Farm Profitability*. Oak Park, Carlow, Ireland: Teagasc. ISBN: 978-1-84170-667-2.

Egan M, Galvin N and Hennessy D (2018) Incorporating white clover (*Trifolium repens* L.) into perennial ryegrass (*Lolium perenne* L.) swards receiving varying levels of nitrogen fertilizer: effects on milk and herbage production. *Journal of Dairy Science* **101**, 3412–3427.

Enriquez-Hidalgo D, Gilliland T and Hennessy D (2016) Herbage and nitrogen yields, fixation and transfer by white clover to companion grasses in grazed swards under different rates of nitrogen fertilization. *Grass and Forage Science* **71**, 559–574.

Enriquez-Hidalgo D, Gilliland TJ, Egan M and Hennessy D (2018) Production and quality benefits of white clover inclusion into ryegrass swards at different nitrogen fertilizer rates. *The Journal of Agricultural Science* **156**, 378–386.

EPA (2019) *Water Quality in Ireland 2013–2018*. Wexford, Ireland: EPA.
Faverdin P, Baratte C, Delagarde R and Peyraud JL (2011) GrazeIn: a model of herbage intake and milk production for grazing dairy cows. 1. Prediction of intake capacity, voluntary intake and milk production during lactation. *Grass and Forage Science* **66**, 29–44.

Fowler D, Coyle M, Skiba U, Sutton MA, Cape JN, Reis S, Sheppard LJ, Jenkins A, Grizzetti B and Galloway JN (2013) The global nitrogen cycle in the twenty-first century. *Philosophical Transactions of the Royal Society B: Biological Sciences* **368**, 20130164.

Fraters D, Van Leeuwen T, Boumans L and Reijs J (2015) Use of long-term monitoring data to derive a relationship between nitrogen surplus and nitrate leaching for grassland and arable land on well-drained sandy soils in the Netherlands. *Acta Agriculturae Scandinavica, Section B – Soil & Plant Science* 65(suppl. 2), 144–154.

Galloway JN, Townsend AR, Erisman JW, Bekunda M, Cai Z, Freney JR, Martinelli LA, Seitzinger SP and Sutton MA (2008) Transformation of the nitrogen cycle: recent trends, questions, and potential solutions. *Science* (New York, N.Y.) **320**, 889–892.

Gilliland T, McGilloway D and Conaghan P (2009) Influence of testing procedure on evaluation of white clover (Trifolium repens L.). *Irish Journal of Agricultural and Food Research* **48**, 227–242.

Godfray HCJ, Beddington JR, Crute IR, Haddad L, Lawrence D, Muir JF, Pretty J, Robinson S, Thomas SM and Toulmin C (2010) Food security: the challenge of feeding 9 billion people. *Science* (New York, N.Y.) **327**, 812–818.

Gourley CJP, Dougherty WJ, Weaver DM, Aarons SR, Awty IM, Gibson DM, Hannah MC, Smith AP and Peverill KI (2012) Farm-scale nitrogen, phosphorus, potassium and sulfur balances and use efficiencies on Australian dairy farms. *Animal Production Science* **52**, 929–944.

Hanrahan L, Geoghegan A, O'Donovan M, Griffith V, Ruelle E, Wallace M and Shalloo L (2017) PastureBase Ireland: a grassland decision support system and national database. *Computers and Electronics in Agriculture* **136**, 193–201.

Hennessy D (2005) *Manipulation of Grass Supply to Meet the Feed Demand of Beef Cattle and Dairy Cows*. Belfast, UK: PhD thesis, Queen's University of Belfast.

Hennessy D, O'Donovan M, French P and Laidlaw A (2008) Manipulation of herbage production by altering the pattern of applying nitrogen fertilizer. *Grass and Forage Science* **63**, 152–166.

Hennessy D, Delaby L, Van Den Pol-Van Dasselaar A and Shalloo L (2020) Increasing grazing in dairy cow milk production systems in Europe. *Sustainability* **12**, 2443.

INRA (ed.) (2010) *Alimentation des bovins, ovins et caprins*. France: Editions Quae.

Kelly P, Shalloo L, Wallace M and Dillon P (2020) The Irish dairy industry – recent history and strategy, current state and future challenges. *International Journal of Dairy Technology* **73**, 309–323.

Ledgard S, Penno J and Sprosen M (1997) Nitrogen balances and losses on intensive dairy farms. In *Proceedings of the New Zealand Grassland Association*, Auckland, pp. 49–53. https://www.nzgajournal.org.nz/index.php/ProNZGA/issue/view/vol59.

Maccarthy D (1982) Nitrogen use for dairy farming. In Moorepark Farmers Dairy Conference, May 1982. Ireland: Teagasc

McCarthy J, Delaby L, Hennessy D, McCarthy B, Ryan W, Pierce K, Brennan A and Horan B (2015) The effect of stocking rate on soil solution nitrate concentrations beneath a free-draining dairy production system in Ireland. *Journal of Dairy Science* **98**, 4211–4224.

McCarthy B, Delaby L, Pierce K, McCarthy J, Fleming C, Brennan A and Horan B (2016) The multi-year cumulative effects of alternative stocking rate and grazing management practices on pasture productivity and utilization efficiency. *Journal of Dairy Science* **99**, 3784–3797.

McClearn B, Gilliland T, Guy C, Dineen M, Coughlan F and McCarthy B (2020a) The effect of perennial ryegrass poity and white clover inclusion on milk production of dairy cows. *Animal Production Science* **60**, 143–147.

McClearn B, Shalloo L, Gilliland T, Coughlan F and McCarthy B (2020b) An economic comparison of pasture-based production systems differing in sward type and cow genotype. *Journal of Dairy Science* **103**, 4455–4465.

McEvoy M, O'donovan M and Shalloo L (2011) Development and application of an economic ranking index for perennial ryegrass cultivars. *Journal of Dairy Science* **94**, 1627–1639.

Mogollón J, Lassaletta L, Beusen A, Van Grinsven H, Westhoek H and Bouwman A (2018) Assessing future reactive nitrogen inputs into global croplands based on the shared socioeconomic pathways. *Environmental Research Letters* **13**, 044008.

Mueller ND, Gerber JS, Johnston M, Ray DK, Ramankutty N and Foley JA (2012) Closing yield gaps through nutrient and water management. *Nature* **490**, 254–257.

O'Brien D, Hennessy T, Moran B and Shalloo L (2015) Relating the carbon footprint of milk from Irish dairy farms to economic performance. *Journal of Dairy Science* **98**, 7394–7407.

O'Donovan M, Delaby L, Stakelum G and Dillon P (2004) Effect of autumn/ spring nitrogen application date and level on dry matter production and nitrogen efficiency in perennial ryegrass swards. *Irish Journal of Agricultural and Food Research* **43**, 31–41.

Peyraud JL, Le Gall A and Lüscher A (2009) Potential food production from forage legume-based-systems in Europe: an overview. *Irish Journal of Agricultural and Food Research* **48**, 115–135.

Ribeiro Filho HMN, Delagarde R and Peyraud J-L (2003) Inclusion of white clover in strip-grazed perennial ryegrass swards: herbage intake and milk yield of dairy cows at different ages of sward regrowth. *Animal Science* **77**, 499–510.

Ribeiro Filho HMN, Delagarde R and Peyraud JL (2005) Herbage intake and milk yield of dairy cows grazing perennial ryegrass swards or white clover/perennial ryegrass swards at low- and medium-herbage allowances. *Animal Feed Science and Technology* **119**, 13–27.

Roche J, Ledgard S, Sprosen M, Lindsey S, Penno J, Horan B and Macdonald K (2016) Increased stocking rate and associated strategic dry-off decision rules reduced the amount of nitrate-N leached under grazing. *Journal of Dairy Science* **99**, 5916–5925.

Ruelle E and Hennessy D (2019) Predicting grass growth on farm, utilisation of the MoSt grass growth model. In O'Brien B., Hennessy D and Shallo L. (eds). *ECPLF*. Cork, Ireland: Teagasc, pp. 262–267.

Ruelle E, Shalloo L, Wallace M and Delaby L (2015) Development and evaluation of the pasture-based herd dynamic milk (PBHDM) model for dairy systems. *European Journal of Agronomy* **71**, 106–114.

Ruelle E, Delaby L, Wallace M and Shalloo L (2016) Development and evaluation of the herd dynamic milk model with focus on the individual cow component. *Animal: An International Journal of Animal Bioscience* **10**, 1986–1997.

Ruelle E, Hennessy D and Delaby L (2018a) Development of the Moorepark St Gilles grass growth model (MoSt GG model): a predictive model for grass growth for pasture based systems. *European Journal of Agronomy* **99**, 80–91.

Ruelle E, Hennessy D and Delaby L (2018b) The Moorepark Saint-Gilles grass growth model (MoSt GG) predicts grass growth under grazing conditions. *Fourrages* **236**, 255–262.

Ruelle E, Delaby L and O'Donovan M (2021) La prévision de la croissance d'herbe en Irlande. In *Journées de l'AFPF 2021, France: AFPF*. pp. 33–39.

Scholefield D, Lockyer D, Whitehead D and Tyson K (1991) A model to predict transformations and losses of nitrogen in UK pastures grazed by beef cattle. *Plant and Soil* **132**, 165–177.

Scholefield D, Tyson K, Garwood E, Armstrong A, Hawkins J and Stone A (1993) Nitrate leaching from grazed grassland lysimeters: effects of fertilizer input, field drainage, age of sward and patterns of weather. *Journal of Soil Science* **44**, 601–613.

Shalloo L and Hanrahan L (2020) Setting targets for the Irish dairy industry. *Animal Production Science* **60**, 159–163.

Steffen W, Richardson K, Rockström J, Cornell SE, Fetzer I, Bennett EM, Biggs R, Carpenter SR, De Vries W and De Wit CA (2015) Planetary boundaries: guiding human development on a changing planet. *Science* (New York, N.Y.) **347**(6223), 1259855.

Stott KJ and Gourley CJP (2016) Intensification, nitrogen use and recovery in grazing-based dairy systems. *Agricultural Systems* **144**, 101–112.

Sutton MA, Bleeker A, Howard CM, Erisman J, Abrol Y, Bekunda M, Datta A, Davidson E, De Vries W and Oenema O (2013) Our nutrient world. The challenge to produce more food & energy with less pollution. Centre for Ecology & Hydrology. https://library.wur.nl/WebQuery/wurpubs/full-text/249094

Ulyatt M and Morley F (1981) *The Feeding Value of Temperate Pastures Grazing Animals*. Amsterdam: Elsevier.
**Wall D** (2019) Nitrogen use efficiency – a key production and environmental indicator. In *National Dairy Conference Proceedings 2019*. Ireland Teagasc.

**Wall D** and **Plunkett M** (2021) *Major and Micro Nutrient Advice for Productive Agricultural crops*. Ireland: Teagasc

**Watson C** and **Atkinson D** (1999) Using nitrogen budgets to indicate nitrogen use efficiency and losses from whole farm systems: a comparison of three methodological approaches. *Nutrient Cycling in Agroecosystems* 53, 259–267.

**Whitehead DC** (1995) *Grassland Nitrogen*. UK: CAB international.