Articles | Volume 16, issue 9
Model description paper
09 May 2023
Model description paper |  | 09 May 2023

Modelling the role of livestock grazing in C and N cycling in grasslands with LPJmL5.0-grazing

Jens Heinke, Susanne Rolinski, and Christoph Müller

To represent the impact of grazing livestock on carbon (C) and nitrogen (N) dynamics in grasslands, we implement a livestock module into LPJmL5.0-tillage, a global vegetation and crop model with explicit representation of managed grasslands and pastures, forming LPJmL5.0-grazing. The livestock module uses lactating dairy cows as a generic representation of grazing livestock. The new module explicitly accounts for forage quality in terms of dry-matter intake and digestibility using relationships derived from compositional analyses for different forages. Partitioning of N into milk, feces, and urine as simulated by the new livestock module shows very good agreement with observation-based relationships reported in the literature. Modelled C and N dynamics depend on forage quality (C:N ratios in grazed biomass), forage quantity, livestock densities, manure or fertilizer inputs, soil, atmospheric CO2 concentrations, and climate conditions. Due to the many interacting relationships, C sequestration, GHG emissions, N losses, and livestock productivity show substantial variation in space and across livestock densities. The improved LPJmL5.0-grazing model can now assess the effects of livestock grazing on C and N stocks and fluxes in grasslands. It can also provide insights about the spatio-temporal variability of grassland productivity and about the trade-offs between livestock production and environmental impacts.

1 Introduction

Grazing lands occupy about 25 % of the global land area (excluding Antarctica; Klein Goldewijk et al.2017) and provide nearly half of the total biomass used in global livestock production (Herrero et al.2013). They also play an important role in Earth's carbon cycle by storing large amounts of soil organic carbon (Conant et al.2017) and contributing to the terrestrial carbon sink (Chang et al.2021). Grazing alters carbon and nitrogen cycling in grasslands (McSherry and Ritchie2013; Conant et al.2017; He et al.2020; Zhou et al.2017), with potential effects on carbon stocks, carbon uptake, N2O emissions, NO3- leaching, and NH3 volatilization. In addition, grazing ruminants produce large amounts of CH4, a potent greenhouse gas (GHG).

Dynamic global vegetation models (DGVMs) are powerful tools for quantifying global biogeochemical cycles, including important aspects such as carbon uptake and loss from the terrestrial biosphere, and their response to changing climate and management (Fernández-Martínez et al.2019; Kondo et al.2020). Yet, most state-of-the-art DGVMs do not account for grazing or only do so in a simplified form by removing a fixed fraction of the above-ground living biomass and, in some cases, returning part of it to grassland litter pools (Friedlingstein et al.2022). A notable exception is ORCHIDEE-GM (Chang et al.2013), which employs the animal module of the process-based grassland model PaSIM (Riedo et al.1998; Vuichard et al.2007) to estimate livestock feed intake from stocking density, animal weight, and biomass availability and to determine the partitioning of C and N into maintenance respiration, products, feces, and urine. However, the effect of feed composition on feed intake and the digestibility of C, N, and energy in feed are not accounted for.

In this paper, we describe the implementation of grazing livestock – represented by lactating dairy cows – into the dynamic global vegetation and agricultural model LPJmL5-tillage (Lutz et al.2019), forming LPJmL5.0-grazing. LPJmL5-tillage incorporates an explicit representation of the N cycle (von Bloh et al.2018) and includes all the improvements of grassland vegetation and carbon dynamics described in Rolinski et al. (2018). The representation of lactating dairy cows is primarily based on established relationships also used in livestock management applications (National Research Council2001) supplemented by relationships from the scientific literature. To account for the effect of forage quality on the digestibility of C and N, as well as digestible energy content, we use compositional data for a wide range of forage plants from Feedipedia (2020) to link these properties to forage N content.

2 Model description

The principal concept of the new livestock module is to determine the mass balances of C and N for grazing dairy cows. Both are arranged in a way that one term can be estimated as the remainder of the balance. In the N balance, that remainder is urinary N excretion (mN,urine), which is calculated as follows:

(1) m N , urine = m N , intake - m N , feces - m N , milk

where mN,intake (Eq. 29) is N intake with forage, mN,feces is N excreted with feces (Eq. 48), and mN,milk is N in milk produced (Eq. 42). In the C balance, the respective remainder is C converted to CO2 through respiration (mC,respiration), which is calculated as follows:

(2) m C , respiration = m C , intake - m C , feces - m C , urine - m C , methane - m C , milk

where mC,intake is C intake with forage (Eq. 30), mC,feces is C excreted with feces (Eq. 49), mC,urine is C excreted with urine (Eq. 50), mC,methane is C in methane from enteric fermentation (Eq. 47), and mC,milk is C in milk produced (Eq. 43).

All terms on the right-hand side of Eqs. (1) and (2) depend directly or indirectly on forage composition, i.e, the mass fractions of crude protein, fatty acids, non-fibre carbohydrates (starch and sugars), fibre carbohydrates (cellulose and hemicellulose), and lignin. However, the only constituents of biomass in LPJmL5.0-grazing are C and N so that all relevant forage properties that depend on forage composition need to be linked to a metric based on C and N content. We use compositional data for a wide range of forage plants from Feedipedia (2020) to calculate the weight fraction of C in forage dry matter and the weight fraction of N in the total mass of C and N (wC,DM and wN,CN, respectively; Sect. 2.1), digestible fractions of C and N in forage (fC and fN; Sect. 2.2), and the digestible energy (DE) content of forage dry matter (Sect. 2.3). We then determine statistical relationships for fC, fN, DE, and wC,DM in relation to wN,CN (Sect. 2.4) to calculate these parameters in LPJmL5.0-grazing.

The remaining part of the method section is structured as follows. Section 2.5 describes the calculation of dry-matter intake, Sect. 2.6 describes the calculation of the energy and protein requirements of livestock, Sect. 2.7 describes the calculation of milk production, Sect. 2.9 describes the calculation of C and N in feces and urine, and Sect. 2.8 describes the calculation of methane emissions from enteric fermentation. Figure 1 provides a schematic representation of the LPJmL5.0-grazing with the new livestock module and its linkages to LPJmL5.0 through fluxes of C and N. In addition, we provide an R script with a fully functional implementation of the livestock model in the Supplement.

Figure 1Schematic representation of LPJmL5.0-grazing. The new livestock module calculates intake of C and N from grazing and determines their partitioning into CO2, methane (CH4), and C and N in milk, urine, and feces. The equations described in the method section can be broadly divided into the following groups: Equations (3)–(19) are used to process compositional data from Feedipedia (2020) to derive Eqs. (20)–(23), which link important variables required in the main model equations (Eqs. 3150) to N in forage (wN,CN). Equations (26)–(30) are used to calculate C and N intake; Eqs. (1)–(2) are the mass balance equations of N and C, respectively; and Eqs. (7)–(10) are weight fractions of C and N in protein, fat, and carbohydrates used to calculate the C and N content of milk.


2.1 Conversion of forage dry matter to C and N

2.1.1 Composition of forage dry matter

Compositional analyses from Feedipedia (2020) provide mass fractions of forage dry matter for crude protein (wCP,DM), neutral-detergent fibre (wNDF,DM), lignin (wL,DM), ether extract (wEE,DM), and ash (wA,DM). wNDF,DM contains a small amount of neutral-detergent-insoluble crude protein (NDICP), which is also included in wCP,DM and needs to be subtracted to obtain the mass fraction of nitrogen-adjusted neutral-detergent fibre (NDF) in forage dry matter (National Research Council2001):

(3) w nNDF , DM = w NDF , DM - w NDICP , DM

where wNDICP,DM is the weight fraction of NDICP in forage dry matter. From data for a wide range of grasses and legumes, Weiss et al. (1992) determined the following relationship for the estimation of the weight fraction of NDICP in forage dry matter:

(4) w NDICP , DM = - 0.0877 + 0.33 w CP , DM + 0.143 w NDF , DM .

From the mass fractions of crude protein, nitrogen-adjusted neutral-detergent fibre, ether extract, and ash, the mass fraction of non-fibre carbohydrates is determined as a residuum (National Research Council2001):

(5) w NFC , DM = 1 - w CP , DM - w nNDF , DM - w EE , DM - w A , DM .

The ether extract obtained from nutritional analysis consists of all lipids in the sample. Some of these lipids, such as pigments and waxes, have almost no nutritional value (Weiss et al.1992). To obtain the weight fraction of highly digestible fatty acids, wEE,DM needs to be adjusted to account for about 1 % of non-fatty-acid components in forage dry matter (National Research Council2001):

(6) w FA , DM = w EE , DM - 0.01 if  w EE , DM > 0.01 0 otherwise .

2.1.2 C and N content of crude protein

Proteins are large macromolecules, which consist of amino acids linked by peptide bonds. There are 20 commonly occurring proteinogenic amino acids with molar weights between 75.1 and 204.2 g mol−1, mass fractions of C between 0.30 and 0.65, and mass fractions of N between 0.08 and 0.32. When amino acids are linked through peptide bonds, one molecule of water is released for each bond. Therefore, anhydrous amino acids in peptide chains have a molar weight that is about 18 g mol−1 lower and a higher mass fraction of C and N of 0.35 to 0.73 and 0.09 to 0.36, respectively.

The C and N content of proteins varies depending on their amino acid composition. Tomé et al. (2019) give a range of 13 % to 19 % for the N content of proteins, but data on C content are rare. Kozlowski (2017) analyzed nearly 14 billion amino acids in nearly 30 million proteins to determine the composition of proteins in 1612 eucaryote species. Using the average amino acid frequencies across all eucaryote species, we estimate a C fraction of 0.53 and an N fraction of 0.17 for average eucaryote protein. Using the amino acid frequencies for the two grass species in the database (Setaria viridis and Eragrostis curvula), we estimate a C fraction of 0.53 and an N fraction of 0.18 for protein in both grass species. Given the close agreement of the estimated C fractions, we assume the following for the weight fraction of C in crude protein:

(7) w C , CP = 0.53 .

The estimates of N content, however, are different and are both larger than the N fraction of 0.16 assumed in compositional analyses (Santos and Huber2002; Feedipedia2020). Because all estimates of crude protein have been calculated from measured N content using that lower value, we also use it to convert crude protein back to N. Thus, the weight fraction of N in crude protein is defined as follows:

(8) w N , CP = 0.16 .

2.1.3 C content of carbohydrates

The basic building blocks of carbohydrates are hexoses (C6H12O6) and pentoses (C5H10O5) with a C fraction of 0.40 each. In polymers, their monomers (hexosans and pentosans) are linked by glycosidic bonds, which are formed by releasing one water molecule for each bond. Therefore, polysaccharides from hexoses and pentoses have a higher C fraction of 0.44 and 0.45, respectively. Disaccharides and oligosaccharides lie in between.

Carbohydrates in plant fibre comprise mainly cellulose and hemicellulose. Cellulose is a polymer of glucose (a hexose), and hemicellulose consists of a mix of hexosans and pentosans (Abu Ghalia and Dahman2017). Information on the composition of non-fibre carbohydrates is scarce, but it can be assumed to be a mix of starch (a polymer of glucose) and various disaccharides and oligosaccharides, which allows for a relatively broad range of possible C fractions. However, the mass fraction of non-fibre carbohydrates in total carbohydrates in the data from Feedipedia (2020) is small (on average 17 %), which limits its relevance for the estimation of C content in forage dry matter. Therefore, we assume the C fraction of starch and celluloses to be a representative value for the weight fraction of all carbohydrates.

(9) w C , CHO = 0.44

2.1.4 C content of fatty acids

The mass fraction of C of fatty acids primarily depends on the fatty-acid chain length. Caprylic acid, a saturated fatty acid with eight C atoms (C8H16O2), has a C fraction of 0.67, while stearic acid, a saturated fatty acid with 18 C atoms (C18H36O2), has a C fraction of 0.76. The degree of saturation also has an effect on C content. For example, linolenic acid, a tri-unsaturated fatty acid with 18 C atoms (C18H30O2), has a C fraction of 0.78. Triglycerides (esters from glycerol and three fatty acids) have a higher C content than single fatty acids of the same type, but the effect becomes negligible with increasing chain length. The triglyceride of linolenic acids has the same C content of 0.78 as linolenic acid.

The composition of ether extract is not reported in the data from Feedipedia (2020). According to a meta-analysis by Glasser et al. (2013), linolenic acid is by far the most abundant fatty acid in grasses, making up more than half of the total lipids. Other important fatty acids are linoleic acid and palmitic acid, each contributing about 10 % to 20 % to the total lipids. The C content of these two fatty acids is slightly lower with 0.77 and 0.75, respectively. Because of the dominance of linolenic acid and the similar C content of other fatty acids, we assume the C fraction of linolenic acid for the weight fraction of all fats and fatty acids as follows:

(10) w C , FA = 0.78 .

2.1.5 C content of lignin

Lignin is a group of large macromolecules derived mainly from three precursors, namely p-coumaryl alcohol, coniferyl alcohol, and sinapyl alcohol (Amthor2003). During polymerization, the monomer units hydroxyphenyl (H), guaiacyl (G), and syringyl (S) are formed from these alcohols, which cross-link through a variety of bonds to form complex three-dimensional macromolecules. The H, G, and S units in lignin have molecular weights of 149.2, 179.2, and 209.2 g mol−1 and C fractions of 0.73, 0.67, and 0.63 (Amthor2003).

The abundance of H, G, and S units in lignin varies depending on the species and tissue type. We obtain the monomer composition for lignin from eight different herbaceous plants from Baucher et al. (1998) and determine their C fraction in lignin using the molecular weight and C fractions of the three monomers. Despite considerable differences in composition, we find all C fractions to be between 0.65 and 0.67, with an overall average of 0.66. Based on that, we define the weight fraction of C in lignin as

(11) w C , L = 0.66 .

2.1.6 C and N content of dry matter

The C content of forage dry matter is calculated from the different forage components multiplied by their respective C fraction:

(12) w C , DM = w C , CP w CP , DM + w C , CHO w NFC , DM + w nNDF , DM - w L , DM + w C , L w L , DM + w C , FA w EE , DM .

Note that, for the calculation of total C content, the weight fraction of total lipids (ether extract) is used, which also comprises pigments and waxes.

The N content in forage dry matter is calculated from crude protein only:

(13) w N , DM = w N , CP w CP , DM .

We also calculate the fraction of N in the sum of C and N:

(14) w N , CN = w N , DM w C , DM + w N , DM .

This variable can also be calculated in LPJmL5.0-grazing, where biomass is represented in terms of C and N only.

2.2 Digestible nutrients

The nutritious value of forage components does not only depend on their energy content but also on how well they can be digested. Weight fractions of digestible nutrients (d) in total dry matter from NFC, CP, FA, and nNDF are calculated using relationships from National Research Council (2001), outlined as follows:


Ash and lignin do not contribute digestible nutrients. wADICP,DM in Eq. (15b) is the mass fraction of acid-detergent-insoluble crude protein (ADICP) in forage dry matter, which is estimated from wNDICP,DM using a relationship from Clipes et al. (2006):

(16) d ADICP = 0.008145 + 0.1131 w NDICP , DM .

For the calculations within LPJmL5.0-grazing, the fraction of digestible C from the total C in dry matter is required, which is calculated as follows:

(17) f C = w C , CP d CP + w C , CHO d NFC + d nNDF + w C , FA w FA , DM w C , DM .

Similarly, the fraction of digestible N from the total N in dry matter is calculated as follows:

(18) f N = w N , CP d CP w N , DM .

2.3 Energy value of forages

The digestible energy (DE) in Mcal kg−1 of forage dry matter is estimated by multiplying the mass fractions of digestible forage components with their energy content (National Research Council2001):

(19) DE = 4.2 d NFC + 5.6 d CP + 9.4 d FA + 4.2 d nNDF - 0.3 .

The heat of combustion is 4.2 Mcal kg−1 for carbohydrates, 5.6 Mcal kg−1 for protein, and 9.4 Mcal kg−1 for long-chain fatty acids (National Research Council2001). Because Eqs. (15a), (15b), (15c), and (15d) give true digestibilities, a correction for metabolic fecal energy is needed, which is assumed to be 0.3 Mcal kg−1 of forage dry matter (National Research Council2001). In a strict sense, this calculation of DE according to Eq. (19) is valid for energy intake at maintenance only. At higher levels of forage intake, the digestibility of diets containing high shares of digestible nutrients (above 60 %) is reduced (National Research Council2001). While this is relevant in highly productive dairy systems, where intake can exceed the maintenance level by a factor of 4, the intake above maintenance in grazing systems is only moderate (see Fig. 3). In addition, the share of digestible nutrients in forages from Feedipedia (2020) is 54 % on average and rarely exceeds 60 %. We therefore assume that this effect is negligible for grazing cattle and do not apply a correction for intake above maintenance.

Figure 2Relationships between forage properties and wN,CN derived from data from Feedipedia (2020). Shaded areas denote 95 % prediction intervals.


2.4 Forage properties in relation to wN,CN

For the application within LPJmL5.0-grazing, we calculate fC, fN, DE, and wC,DM for a wider range of forage plants from Feedipedia (2020) and determine their relationship with wN,CN (Fig. 2). We test three different types of functional relationships: (i) no relationship with wN,CN (i.e. the mean of the dependent variable), (ii) a linear relationship, and (iii) a non-linear exponential relationship with three parameters. All three relationships are determined for each forage property, and the relationship with the lowest value for the Akaike's information criterion (AIC) is selected.

Figure 3Mass balances of nitrogen and carbon as a function of wN,CN for a single livestock unit of 500 kg. Dry matter intake is assumed not to be limited by forage availability, and unfulfilled requirements of NE and MP from previous days are ignored.


For the digestibility of C (fC), we find the lowest AIC for a linear relationship with wN,CN.

(20) f C = 0.561 + 2.190 w N , CN

The estimated parameters are statistically significant, with p values <0.001. The residual standard error (RSE) of the fitted model is about 0.036, which translates into a 95 % prediction interval of ±10.8 % at wN,CN=0.05.

For the digestibility of N (fN), we find the lowest AIC for a concave exponential relationship with wN,CN.

(21) f N = 0.914 - 0.494 exp - 59.559 w N , CN

All estimated parameters are statistically significant, with p values <0.001. The RSE is small (0.0078), which translates into a narrow 95 % prediction interval of ±1.8 % at wN,CN=0.05.

Similarly to the digestibility of C, digestible energy (DE) in forage increases linearly with wN,CN.

(22) DE = 1.952 + 11.438 w N , CN

Also here, the estimated parameters are statistically significant, with p values <0.001. Similarly to the model of C digestibility, the RSE is relatively large (0.15), which results in a similar 95 % prediction interval of ±12.0 % at wN,CN=0.05.

C content of forage dry matter (wC,DM) is found to be independent of wN,CN.

(23) w C , DM = 0.424

However, the sample standard deviation of wC,DM is only 0.01, which translates into a 95 % prediction interval of ±4.7 %.

2.5 Dry-matter intake

Relationships for predicting dry-matter intake of lactating dairy cows usually include milk production as an independent variable because their purpose is to determine the amount and composition of feed required to achieve a desired milk yield (National Research Council2001). To meet the nutritional requirements of high-yielding dairy cows, their diets must contain high proportions of concentrates (e.g. maize and barley) to achieve a high concentration of readily available nutrients. Such relationships are obviously inappropriate for the estimation of voluntary forage intake of grazing cows, which is limited by the capacity of the rumen to digest fibrous materials rather than the metabolic capacity of the animal to utilize the available energy (Tedeschi et al.2019). Several relationships for predicting ad libitum forage intake of grazing cattle have been proposed, but most of them have not been explicitly developed for lactating cows. Forage intake of lactating cows is higher than for dry cows of equal size because lactation causes an increase in the size of the gastrointestinal tract (Coleman et al.2014). Tulloh (1966) has measured the size of the gastrointestinal tract of twin pairs of lactating and dry cows and found that the weight of the whole tract per body weight was on average 33.1 % higher in lactating cows. The weight and the water-filled volume of the reticulo-rumen per body weight were on average about 21.3 % and 44.6 % higher in lactating cows. We choose a relationship from Coleman et al. (2014), obtained from the data of Moore et al. (1999), that predicts voluntary forage intake per kg of body weight as a function of crude protein content in forages. Because the relationship was derived from data of dry dairy cows, we multiply body weight by 1.33 for the calculation of total daily dry-matter intake (DMImax) to account for the larger gastrointestinal tract of lactating cows:

(24) DMI max = 1.33 BW 0.0235 - 0.0385 exp - 32 w CP , DM ,

where DMImax is in kg d−1, BW is the body weight in kg, and wCP,DM is the mass fraction of crude protein in forage dry matter. The equation predicts a maximum daily forage intake of 2.82 % of body weight or 14.1 kg d−1 dry matter for a lactating cow weighing 500 kg. At low levels of crude protein in the feed (wCP,DM<0.10), the ability of ruminal microbes to break down fibrous material is affected by low N availability, which leads to lower passage rates and lower feed intake (Coleman2005).

The protein content in forage dry matter (wCP,DM) required in Eq. (24) for the estimation of DMImax is calculated by dividing the N content of forage dry matter (wN,DM) by the N content of crude protein (CP; see Sect. 2.1.2):

(25) w CP , DM = w N , DM w N , CP .

wN,DM can be calculated from wN,CN by rearranging Eq. (14):

(26) w N , DM = w C , DM w N , CN 1 - w N , CN .

We note that an alternative model for estimating the voluntary feed intake of dairy cows has been proposed by Faverdin et al. (2011). Unlike our approach, which is based on the NRC energy and protein system (National Research Council2001), their model is based on the INRA fill unit system (Institut national de la recherche agronomique2007) and requires iteration between multiple equations of that system. Therefore, it is not applicable here.

The dry-matter intake calculated in Eq. (24) only considers limitations by the animal and does not account for the effect of biomass availability. When biomass availability declines, the grazing efficiency of the animals is reduced because less biomass is acquired with each bite (Hodgson1985). To account for this effect, we adopt the sigmoid grazing function used in the Hurley model (Johnson and Parsons1985), which describes the decline in the proportion of DMImax as a function of leaf area index (LAI):

(27) DMI = DMI max LAI / K q 1 + LAI / K q

where q is a parameter that determines the slope of the curve, and K is a position parameter that can be interpreted as the LAI value at which half of the DMImax is achieved. For q, a value of 3 is suggested to give realistic results (Johnson and Parsons1985; Herrero et al.2000). For the calculation of K, we adapt a relationship from Herrero et al. (2000), which accounts for the effect of animal size (body weight):

(28) K = 0.229 BW 0.36 .

For a cow of 500 kg, a value of 2.15 for K is obtained.

The daily intake of N per head (mN,intake) is calculated by multiplying dry-matter intake (DMI) with the N content of forage dry matter as follows:

(29) m N , intake = DMI w N , DM .

The daily intake of C per head (mC,intake) is calculated by multiplying DMI with the C content of forage dry matter:

(30) m C , intake = DMI w C , DM .

2.6 Energy and protein requirements

Calculations of energy requirements in National Research Council (2001) are based on net energy units, which represent energy available to the animal after all losses have been subtracted. Part of the digestible energy is lost as urine and methane, which reduces the amount of energy that is actually metabolized by the body (Weiss2011). Metabolizable energy (ME) in Mcal kg−1 of forage dry matter is calculated from DE as follows (National Research Council2001):

(31) ME = 1.01 DE - 0.45 .

Part of the metabolizable energy is lost as heat. Net energy (NE) in Mcal kg−1 of forage dry matter that is available for maintenance, activity, and lactation is calculated as follows (National Research Council2001):

(32) NE = 0.703 ME - 0.19 .

Over a range for wN,CN from 0.02 to 0.1, the ratio of NE to DE increases from 50 % to 57 %.

Daily net energy requirements for maintenance are proportional to metabolic body weight BW0.75 (National Research Council2001).

(33) NE M = 0.08 BW 0.75

The value of 0.08 Mcal d−1 kg−0.75 metabolic body weight accounts for the increased maintenance requirements of lactating cows and includes a 10 % activity allowance (National Research Council2001). Although additional energy requirements for grazing and walking are not included in this allowance (National Research Council2001), they are ignored here and will be addressed in a subsequent version of our model.

The requirements of metabolizable protein (MP) consist of urinary protein requirements (MPUP) and metabolic fecal protein (MPMFP). MPUP is proportional to BW0.5 (National Research Council2001):

(34) MP UP = 0.0041 BW 0.5 ,

whereas MPMFP is proportional to DMI (National Research Council2001):

(35) MP MFP = 0.03 DMI .

2.7 Milk production

Net energy requirements for milk production (lactation) equal the energy content of milk (National Research Council2001):

(36) NE milk = 0.36 + 9.69 w fat , milk ,

where wfat,milk is the weight fraction of fat in milk. We assume a constant fat content of 4 % in milk (wfat,milk=0.04), which corresponds to the net energy requirements for milk production of 0.748 Mcal kg−1.

The amount of milk that can be produced from the available net energy above maintenance requirements is calculated as follows:

(37) m milk , NE = DMI NE - NE M NE milk .

The requirements of metabolizable protein for lactation are calculated from milk protein content assuming a conversion efficiency of metabolizable protein to milk protein of 0.67 (National Research Council2001):

(38) MP milk = w protein , milk 0.67 ,

where wprotein,milk is the weight fraction of protein in milk. Assuming a constant protein content of 3.2 % in milk (wprotein,milk=0.032) gives metabolizable protein requirements of 0.048 kg kg−1 for milk production.

The amount of milk that can be produced from the available MP above maintenance is calculated as follows:

(39) m milk , MP = MP avl - MP UP - MP MFP MP milk ,

where MPavl is the total available metabolizable protein calculated from digested N:

(40) MP avl = DMI w N , DM f N w N , CP .

The actual amount of milk that can be produced is the minimum of mmilk,NE and mmilk,MP:

(41) m milk = min ( m milk , NE , m milk , MP ) .

The amount of N contained in milk (mN,milk) is calculated by multiplying with the mass fraction of N in milk (wN,milk):

(42) m N , milk = w N , milk m milk .

In the same way, the amount of C is determined as follows:

(43) m C , milk = w C , milk m milk .

Weight fractions of N (wN,milk) and C wC,milk in milk can be determined from milk composition:

(44) w N , milk = w N , CP w protein , milk ,


(45) w C , milk = w C , CP w protein , milk + w C , CHO w CHO , milk + w C , FA w fat , milk .

Assuming a constant composition of 4 % fat, 3.2 % protein, and 4.85 % sugar in milk, we obtain estimates for wN,milk and wC,milk of 0.00512 and 0.0695, respectively.

2.8 Methane from enteric fermentation

Enteric fermentation in the rumen produces methane as a byproduct. According to IPCC guidelines for greenhouse gas inventories, 6.5±1.0 % of gross energy intake is converted to methane (IPCC2006). The lower and upper bounds of this range are described to be appropriate for good and poorer feed, respectively, but a quantitative relationship is not given by IPCC (2006). Hence, we calculate methane production (mmethane in kg d−1) assuming a constant methane conversion factor of 6.5 % of gross energy intake:

(46) m methane = DMI 18.4 0.065 55.6 ,

where 18.4 and 55.6 are the gross energy content of feed and methane, respectively, in MJ kg−1. Note that, in light of the large uncertainty entailed by assuming a constant emission factor, we do not account for possible variations in gross energy content for forage here.

Since C makes up 75 % of the molar weight of methane (12 g mol−1 out of 16 g mol−1), the amount of C converted to methane is calculated as follows:

(47) m C , methane = 0.75 m methane .

2.9 Feces and urine

N excreted with feces comprises indigestible N in feed and the N contained in metabolic fecal protein:

(48) m N , feces = m N , intake 1 - f N + w N , CP MP MFP .

Analogously, C excreted with feces is calculated as follows:

(49) m C , feces = m C , intake 1 - f C + w C , CP MP MFP .

N excreted with urine is calculated as the residual of the animal's N balance (Eq. 1). The amount of C excreted with nitrogenous components in urine is calculated by multiplying mN,urine with an appropriate C:N ratio. According to Dijkstra et al. (2013), the majority of N in urine (50 %–90 %) is present as urea (Dijkstra et al.2013), which has a C:N ratio of 0.5. Around 5 % of N is present as hippuric acid, with a C:N ratio of 9. The remainder are purine derivatives, creatine, and creatinine, which all have C:N ratios between 1 and 1.33. This implies a plausible range of average C:N ratios in urine of about 0.95 to 1.3. For simplicity and because the amount of C excreted with urine generally makes up a small part of the C balance, we assume a C:N ratio of 1 in urine.

(50) m C , urine = 1 m N , urine

C and N excreted with feces are added to the respective aboveground litter pools in LPJmL. Nitrogenous compounds in urine are assumed to be quickly degraded to ammonium. Thus, mN,urine is added to the ammonium pool of the top soil layer in LPJmL, while mC,urine is added to the aboveground litter pool.

2.10 Integration in LPJmL5.0-grazing

The new module to calculate grazing, digestion, and returning C and N to the soil in form of feces and urine is implemented as a new harvest function, harvest_grass_grazing_ext_livestock, in the source code file harvest_stand.c. The function is called daily in the daily_grassland function when the grassland management option GS_GRAZING_EXT is set in the configuration file lpjml.js. The dairy cow representation is only compatible with the grassland implementation, as described by Rolinski et al. (2018). The new functions of LPJmL5.0-grazing have been implemented in the LPJmL model version LPJmL5.0-tillage, as described by Lutz et al. (2019).

In situations where grass biomass availability exhibits a strong seasonality, it can happen that forage intake is insufficient to fulfil daily maintenance requirements of net energy and protein at all times of the year. In order to prevent overly optimistic milk yield estimates under such circumstances, the unfulfilled daily requirements are tracked by adding the deficit to a buffer (one for net energy and one for protein). When energy and protein intake are above maintenance requirements, the buffers are balanced first; milk production can only occur when both deficit buffers are zero and when energy and protein intake are above maintenance requirements. To prevent the accumulation of a large deficit during long deficit periods (e.g, during spin-up) and the impediment of milk production during succeeding, more productive periods, we constrain the size of the deficit buffers of energy and protein to 365⋅NEM and 365⋅MPUP, respectively.

2.11 Modelling protocol

For the evaluation of the intake and production model alone (without interaction with LPJmL), we calculate C and N uptake and their division into C and N contained in milk, urine, feces, methane, and CO2 for 1000 values of wN,CN between 0.015 and 0.09 (the allowed range for grass biomass in LPJmL). We assume a body weight of 500 kg per livestock unit (LSU), unlimited grass availability, and no unfulfilled requirements of NE and MP from previous days.

For the simulations of the livestock model within LPJmL5.0-grazing, we also assume a body weight of 500 kg per LSU. Simulations are performed at 0.5 arcdeg resolution (about 55 km at the Equator) for all global land cells except Antarctica. The model is driven by climate forcing from the GSWP3-W5E5 dataset (Kim2017; Cucchi et al.2020; Lange et al.2022), historical atmospheric deposition of NO3- and NH4+ (Yang and Tian2020), and historic atmospheric CO2 concentrations (Büchner and Reyer2022). We perform simulations with 45 different LSU densities from 0 to 4 LSU ha−1 of grazing area; the increment between scenarios increases with increasing LSU densities from 0.01 to 0.2 LSU ha−1. For each LSU density setting, the model is run for 7000 years using a random permutation of years from the 1901–1931 period to bring C and N pools into equilibrium. After that, the model is run from 1901–2016 with transient climate, atmospheric deposition, and atmospheric CO2 concentration.

3 Evaluation

Intake of C and N per cow, as well as their division into C and N contained in milk, urine, feces, methane, and CO2, changes as a function of wN,CN (Fig. 3). For a LSU of 500 kg, the intake of C is 3.7 kg d−1 for wN,CN=0.015 and increases steeply with increasing wN,CN until it reaches 6 kg d−1 at wN,CN=0.033, which is about 90 % of the maximum daily C intake. Further increases in wN,CN have a comparatively small impact on C intake. In contrast, N intake increases almost linearly with wN,CN, dominated by the increasing N content of grazed biomass. Because the digestibility of C and N increases with wN,CN (Fig. 2), the proportion of intake excreted with feces decreases, which results in a decline in the absolute amount of C in feces for wN,CN>0.038. Net energy intake is sufficient to meet energy requirements for maintenance over the considered range of wN,CN, but metabolizable protein corresponding to digested N is insufficient to meet urinary and metabolic fecal protein requirements below wN,CN=0.019. Thus, milk production only occurs for wN,CN>0.019, where it is limited by available N up to wN,CN=0.32 and by available energy for higher wN,CN.

To evaluate the validity of our model, we compare the simulated partitioning of dietary N as a function of wN,CN to corresponding relationships from Huhtanen et al. (2008) (Fig. 4). Huhtanen et al. (2008) analyse relationships between feed properties and nitrogen utilization and partitioning in dairy cows obtained from 998 data points from 207 lactation trials. For nitrogen use efficiency (i.e. the amount of N in milk divided by total N intake), Huhtanen et al. (2008) provide nine different relationships, of which we selected the one with the lowest residual mean square error (RMSE) and lowest AIC.

(51) m N , milk m N , intake = 0.627 - 33.9 w CP , DM ME + 650 w CP , DM ME 2

For N in feces, urine, and manure (feces + urine) and the fraction of urinary N in manure N, up to five different relationships are provided in Huhtanen et al. (2008), of which some contained independent variables for which we had no corresponding estimate. From the remaining relationships, we chose the ones with the lowest RMSE and AIC. N in feces is related to DMI and mN,intake.

(52) m N , feces = - 0.021 + 0.00673 DMI + 0.101 m N , intake

N in urine is related to mN,intake, DMI, and mmilk.

(53) m N , urine = 0.04 + 0.879 m N , intake - 0.009 DMI - 0.0039 m milk

N in manure (mN,manure=mN,feces+mN,urine) is related to mN,intake, DMI, and ME.

(54) m N , manure = 0.081 + 0.947 m N , intake - 0.0059 DMI - 0.0059 ME

The fraction of urinary N in manure N is quadratically related to wCP,DM.

(55) m N , urine m N , manure = - 0.241 + 7.11 w CP , DM - 13 w CP , DM 2

To obtain the relationships with wN,CN shown in Fig. 4, wCP,DM, ME, DMI, mN,intake, and mmilk are calculated from wN,CN using the relationships described in Sect. 2.

There is close agreement between the relationships simulated by our model and those determined by Huhtanen et al. (2008). The partitioning of N intake into N in milk, N in feces, and N in urine as a function of wN,CN in our model is strongly supported by the relationships derived from trials by Huhtanen et al. (2008). Unfortunately, we were not able to obtain similar reference data for the evaluation of C partitioning. However, C:N ratios for all elements of the N balance are well determined or closely related to wN,CN in forage. Therefore, the comparison in Fig. 4 also provides an indirect validation of the C balance.

Figure 4Comparison of simulated partition of dietary nitrogen as a function of wN,CN to corresponding relationships from Huhtanen et al. (2008) obtained from experimental data. Simulation are performed assuming forage intake not limited by forage availability and ignoring unfulfilled requirements of NE and MP from previous days.


Figure 5(a) and C (b) budgets for the whole grassland system as a function of livestock density in a grid cell near Potsdam (52.25 N, 13.25 E). The lower boundary represent changes in C and N soil storage. The upper boundaries represent C and N fluxes into the system: gross primary productivity for C and atmospheric deposition and biological N fixation for N. Atmospheric deposition of N at this site is 2.15 g m−2 yr−1. All values represent averages for 1971–2016. The lower end of the C budget (upper right) is shown in greater detail in a separate graph (lower right).


Figure 6(a) and C (b) budgets for the dairy herd as a function of livestock density in a grid cell near Potsdam (52.25 N, 13.25 E). All values represent averages for 1971–2016.


Figure 7GHG emissions (a) and emission intensity (b) as a function of livestock density in a grid cell near Potsdam (52.25 N, 13.25 E). Emission intensity is defined as the net emissions of livestock grazing (i.e. actual emissions minus emissions for 0 LSU ha−1) per kg of protein produced. All values represent averages for 1971–2016.


Figure 8Global maps of milk production (top), emission intensity per kg of milk protein (middle), and stocking density (bottom) for the stocking densities yielding the highest milk production (left) and the lowest emission intensity per kg of milk protein (left). Only grid cells with pastures or rangeland in the year 2000 according to the HYDE3.2 dataset (Klein Goldewijk et al.2017) are shown. All values represent averages for 1971–2016.

4 Results

With the implementation of grazing dairy systems in LPJmL, including the effect of forage quality and quantity on the uptake and partitioning of C and N, the model can now explicitly represent the effects of grazing management on land productivity and the C and N budgets.

In LPJmL5, gross primary productivity (GPP) is computed on a daily basis depending on climate, atmospheric CO2 concentration, leaf area index (LAI), and plant-available water and mineral N in the soil. One direct effect of grazing in LPJmL is a reduction in leaf area, which tends to reduce GPP but also autotrophic respiration. Depending on the relative strength of both effects, this can lead to an increase or a decrease in net primary productivity (i.e. the amount of carbon available for allocation to leaves and roots). Another important direct effect of grazing is the partial withdrawal of grazed N from the ecosystem (as milk) along with a transfer of the remaining N in other compartments (organic litter N and soil ammonium pools). This can result in an increase or a decrease of plant-available mineral N in soils, depending on the relative strength of the two effects. Other processes in LPJmL respond in various ways to these primary impacts of grazing and lead to changes in all stocks and flows of carbon, nitrogen, and water in the model. It should be noted, however, that the current implementation considers neither the mechanical impacts of grazing (e.g. trampling or soil compaction) nor the preferential grazing of certain grass types or compartments.

As examples for site- and management-specific dynamics, we show C and N balances of grassland (Fig. 5), C and N partitioning by the dairy herd (Fig. 6), and GHG emissions (Fig. 7) for Potsdam, Germany in the main text and for three additional sites in the Sahel of Burkina Faso (14.25 N, 0.25 E); in Cordoba, Argentina (31.75 S, 62.75 W); and in Riau, Indonesia (1.25 N, 101.75 E) in the Appendix (Figs. A1A9). The different sites exhibit a wide range of responses to changes in stocking density, which are caused by the interplay between direct grazing impacts and their subsequent effect on other processes in LPJmL5.0-grazing. A robust response across sites is a reduction of autotrophic respiration with increasing stocking density, which is the result of the reduction of standing biomass (Figs. 5, A1A3). The reduction in biomass (leaf area), however, does not always lead to a reduction in GPP, which can be attributed to an increase in plant-available mineral N in the soils as a result of grazing. An increase in N and C intake with stocking density (Figs. 6, A4A6) is to be expected, but the shape of the relationship varies among sites and does not always monotonically increase, reflecting variations in grass availability and quality (N content). Methane emissions increase proportionally with C intake. Milk production first increases with stocking density but then decreases again when increases in intake become smaller than increases in maintenance requirements. At the site in Indonesia (Fig. A6), this point apparently lies beyond 4 LSU ha−1, which is the largest stocking density used in these simulations. The N withdrawal from the ecosystem with milk production is clearly reflected in a reduction in N losses. In fact, NO3 leaching and N2O emissions are usually smallest at or near to the stocking density, for which the highest milk yield is achieved (Figs. 5, A1A3). This is very much in contrast to carbon sequestration in soils in response to increasing atmospheric CO2 concentrations, which is typically largest for much lower stocking densities.

At medium to high stocking densities (above 1–2 LSU ha−1), total GHG emissions (expressed in CO2 equivalents for a 100-year global warming potential) are usually dominated by methane emissions (Figs. 7, A7A9). However, at low stocking densities, enhanced carbon sequestration in soils and reduced N2O emissions through grazing can offset methane emissions and, in some cases, even lead to a decrease in total GHG emissions (e.g. Fig. A8). Thus, there is usually a trade-off between maximizing productivity and minimizing GHG intensity (net emissions per kg of protein) at lower productivity levels. To highlight this tradeoff and to demonstrate that LPJmL5.0-grazing can be used to assess such trade-offs on a global scale, we use the set of simulations described in Sect. 2.11 to determine for each grid cell the stocking densities that correspond to the highest milk production and lowest GHG intensity (Fig. 8). The results show that GHG intensity of livestock production can indeed be negative in many locations at low or very low stocking densities but at the cost of much lower livestock productivity than is possible. It is important to note that these results should be interpreted with care because dairy grazing may not be the most appropriate production system in all locations. However, the change in GHG emissions in response to grazing can be expected to be similar for other livestock grazing systems, and the level of milk production provides an indication of the level of grassland productivity in general.

5 Discussion and conclusions

The new implementation of grazing dairy cows in LPJmL5.0-grazing greatly improves the representation of C and N cycling in grasslands in relation to livestock grazing. The importance of forage quality, in addition to forage quantity, for the uptake and partitioning of C and N by ruminants is well known (e.g. National Research Council2001; Institut national de la recherche agronomique2007), but this has, to our knowledge, never been represented in a dynamic global vegetation model, such as LPJmL. Limitations of the model arise from the generic representation of grazing systems by grazing dairy systems, which are not the most important systems globally (Herrero et al.2013; Heinke et al.2020). However, this system, with its continuous productivity, is much easier to describe and parametrize than grazing beef cattle, which require modelling the weight gain of individual animals and also herd dynamics. In addition, the effect of supplementation with feed crops is currently not included in the model.

Despite these limitations, LPJmL5.0-grazing can be applied to assess the impacts of grazing on C and N cycles – including carbon sequestration – for given levels of grazing intensity. For such analyses, the representation of grazing systems by grazing dairy systems is a reasonable simplification. Other possible applications include the determination of stocking densities to fulfil predefined targets (e.g. maximum productivity or lowest GHG intensity; see Fig. 8). But the results of these kinds of assessments need to be carefully interpreted in light of the generic representation by grazing dairy systems. However, such analyses can provide valuable insights into the spatial variations of grassland productivity and GHG intensity, as well as the synergies and trade-offs entailed by pasture management, which can be generalized to other forms of livestock production.

Appendix A: Additional site dynamics

What follows are plots of C and N dynamics at additional sites for a marginal site in the Sahel of Burkina Faso (14.25 N, 0.25 E); a site in Cordoba, Argentina (31.75 S, 62.75 W); and a productive site in Riau, Indonesia (1.25 N, 101.75E) in comparison to the plots for Potsdam, Germany (52.25 N, 13.25 E) in the main text.

Figure A1Same as Fig. 5 in the main text but for a site in the Sahel of Burkina Faso (14.25 N, 0.25 E). Atmospheric deposition of N at this site is 0.37 g m−2 yr−1.


Figure A2Same as Fig. 5 in the main text but for a site in Argentina (31.75 S, 62.75 W). Atmospheric deposition of N at this site is 0.40 g m−2 yr−1.


Figure A3Same as Fig. 5 in the main text but for a site in Indonesia (1.25 N, 101.75 E). Atmospheric deposition of N at this site is 4.99 g m−2 yr−1.


Figure A4Same as Fig. 6 in the main text but for a site in the Sahel of Burkina Faso (14.25 N, 0.25 E).


Figure A5Same as Fig. 6 in the main text but for a site in Argentina (31.75 S, 62.75 W).


Figure A6Same as Fig. 6 in the main text but for a site in Indonesia (1.25 N, 101.75 E).


Figure A7Same as Fig. 7 in the main text but for a site in the Sahel of Burkina Faso (14.25 N, 0.25 E).


Figure A8Same as Fig. 7 in the main text but for a site in Argentina (31.75 S, 62.75 W).


Figure A9Same as Fig. 7 in the main text but for a site in Indonesia (1.25 N, 101.75 E).


Appendix B: List of symbols
Symbol Unit Description
mN,urine kg d−1 Mass of nitrogen in urine
mN,intake kg d−1 Mass of nitrogen in forage intake
mN,feces kg d−1 Mass of nitrogen in feces
mN,milk kg d−1 Mass of nitrogen in milk
mC,respiration kg d−1 Mass of carbon used for maintenance respiration
mC,intake kg d−1 Mass of carbon in forage intake
mC,feces kg d−1 Mass of carbon in feces
mC,urine kg d−1 Mass of carbon in urine
mC,methane kg d−1 Mass of carbon in methane
mC,milk kg d−1 Mass of carbon in milk
wnNDF,DM kg kg−1 Weight fraction of nitrogen-adjusted neutral-detergent fibre in forage dry matter
wNDF,DM kg kg−1 Weight fraction of neutral-detergent fibre in forage dry matter
wNDICP,DM kg kg−1 Weight fraction of neutral-detergent-insoluble crude protein in forage dry matter
wCP,DM kg kg−1 Weight fraction of crude protein in forage dry matter
wNFC,DM kg kg−1 Weight fraction of non-fibre carbohydrates in forage dry matter
wEE,DM kg kg−1 Weight fraction of ether extract in forage dry matter
wA,DM kg kg−1 Weight fraction of ash in forage dry matter
wFA,DM kg kg−1 Weight fraction of fatty acids in forage dry matter
wL,DM kg kg−1 Weight fraction of lignin in forage dry matter
wADICP,DM kg kg−1 Weight fraction of acid-detergent-insoluble crude protein in forage dry matter
wC,CP kg kg−1 Weight fraction of carbon in crude protein
wN,CP kg kg−1 Weight fraction of nitrogen in crude protein
wC,CHO kg kg−1 Weight fraction of carbon in carbohydrates
wC,L kg kg−1 Weight fraction of carbon in lignin
wC,FA kg kg−1 Weight fraction of carbon in fatty acids
wC,DM kg kg−1 Weight fraction of carbon in forage dry matter
Symbol Unit Description
wN,DM kg kg−1 Weight fraction of nitrogen in forage dry matter
wN,CN kg kg−1 Weight fraction of nitrogen in sum of carbon and nitrogen
dNFC kg kg−1 Weight fractions of digestible non-fibre carbohydrates in total dry matter
dCP kg kg−1 Weight fractions of digestible crude protein in total dry matter
dFA kg kg−1 Weight fractions of digestible fatty acids in total dry matter
dnNDF kg kg−1 Weight fractions of digestible nitrogen-adjusted neutral-detergent fibre in total dry matter
dADICP kg kg−1 Weight fractions of digestible acid-detergent-insoluble crude protein in total dry matter
fC kg kg−1 Digestible fraction of carbon in total dry matter
fN kg kg−1 Digestible fraction of nitrogen in total dry matter
BW kg Body weight
DMImax kg d−1 Dry-matter intake capacity
DMI kg d−1 Dry-matter intake
q Slope parameter in Eq. (27)
K Position parameter in Eq. (27)
DE Mcal kg−1 Digestible energy in forage dry matter
ME Mcal kg−1 Metabolizable energy in forage dry matter
NE Mcal kg−1 Net energy in forage dry matter
NEM Mcal d−1 Net energy requirements for maintenance
NEmilk Mcal kg−1 Net energy requirements for milk production
MPavl kg d−1 Available metabolizable protein
MPUP kg d−1 Metabolizable protein requirements for urinary protein
MPMFP kg d−1 Metabolizable protein requirements for metabolic fecal protein
Symbol Unit Description
MPmilk kg kg−1 Metabolizable protein requirements for milk production
mmilk,NE kg d−1 Potential milk production from net energy
mmilk,MP kg d−1 Potential milk production from metabolizable protein
mmilk kg d−1 Milk production
wfat,milk kg kg−1 Weight fractions of fat in milk
wprotein,milk kg kg−1 Weight fractions of protein in milk
wCHO,milk kg kg−1 Weight fractions of carbohydrates in milk
wC,milk kg kg−1 Weight fractions of carbon in milk
wN,milk kg kg−1 Weight fractions of nitrogen in milk
mmethane kg d−1 Methane production
Code availability

The source code of LPJmL5.0-grazing is archived at Zenodo under (Heinke et al.2022).

Data availability

The historical climate data from the GSWP-W5E5 dataset are available from (Lange et al.2022). The historical data of atmospheric N deposition and atmospheric CO2 concentrations can be obtained from (Yang and Tian2020) and (Büchner and Reyer2022), respectively. All input data, model code, model outputs, and scripts that have been used to produce the results presented in this paper are archived at the Potsdam Institute for Climate Impact Research and are available upon request.


The supplement related to this article is available online at:

Author contributions

JH developed the model, performed simulations, and analysed results. JH, SR, and CM wrote the paper.

Competing interests

At least one of the (co-)authors is a member of the editorial board of Geoscientific Model Development. The peer-review process was guided by an independent editor, and the authors also have no other competing interests to declare.


Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


Jens Heinke acknowledges funding from the EXIMO project (grant no. 01LP1903D) funded through the German Ministry for Education and Research (BMBF). Susanne Rolinski acknowledges funding from the projects AGrEc (grant no. 01DG21039) and CLIMASTEPPE (grant no. 01DJ18012), both funded through the BMBF.

Financial support

This research has been supported by the Bundesministerium für Bildung und Forschung (grant nos. 01LP1903D, 01DJ18012, and 01DG21039).

The publication of this article was funded by the Open Access Fund of the Leibniz Association.

Review statement

This paper was edited by Marko Scholze and reviewed by three anonymous referees.


Abu Ghalia, M. and Dahman, Y.: Synthesis and utilization of natural fiber-reinforced poly (lactic acid) bionanocomposites, in: Lignocellulosic Fibre and Biomass-Based Composite Materials, Elsevier, 313–345,, 2017. a

Amthor, J. S.: Efficiency of Lignin Biosynthesis: a Quantitative Analysis, Ann. Botany, 91, 673–695,, 2003. a, b

Baucher, M., Monties, B., Montagu, M. V., and Boerjan, W.: Biosynthesis and Genetic Engineering of Lignin, Crit. Rev. Plant Sci., 17, 125–197,, 1998. a

Büchner, M. and Reyer, C.: ISIMIP3a atmospheric composition input data (v1.2), ISIMIP Repository [data set],, 2022. a, b

Chang, J., Ciais, P., Gasser, T., Smith, P., Herrero, M., Havlík, P., Obersteiner, M., Guenet, B., Goll, D. S., Li, W., Naipal, V., Peng, S., Qiu, C., Tian, H., Viovy, N., Yue, C., and Zhu, D.: Climate warming from managed grasslands cancels the cooling effect of carbon sinks in sparsely grazed and natural grasslands, Nat. Commun., 12, 118,, 2021. a

Chang, J. F., Viovy, N., Vuichard, N., Ciais, P., Wang, T., Cozic, A., Lardy, R., Graux, A.-I., Klumpp, K., Martin, R., and Soussana, J.-F.: Incorporating grassland management in ORCHIDEE: model description and evaluation at 11 eddy-covariance sites in Europe, Geosci. Model Dev., 6, 2165–2181,, 2013. a

Clipes, R., Detmann, E., Silva, J., Vieira, R., Nunes, L., Lista, F., and Ponciano, N.: Evaluation of acid detergent insoluble protein as an estimator of rumen non-degradable protein in tropical grass forages, Arquivo Brasileiro de Medicina Veterinária e Zootecnia, 58, 694–697,, 2006. a

Coleman, S. W.: Predicting Forage Intake by Grazing Ruminants, Proceedings of 2005 Florida Ruminant Nutrition Symposium, 1–2 February 2005, Gainesville, Florida, 72–90, 2005. a

Coleman, S. W., Gunter, S. A., Sprinkle, J. E., and Neel, J. P.: Difficulties associated with predicting forage intake by grazing beef cows, J. Anim. Sci., 92, 2775–2784,, 2014. a, b

Conant, R. T., Cerri, C. E. P., Osborne, B. B., and Paustian, K.: Grassland management impacts on soil carbon stocks: a new synthesis, Ecol. Appl., 27, 662–668,, 2017. a, b

Cucchi, M., Weedon, G. P., Amici, A., Bellouin, N., Lange, S., Müller Schmied, H., Hersbach, H., and Buontempo, C.: WFDE5: bias-adjusted ERA5 reanalysis data for impact studies, Earth Syst. Sci. Data, 12, 2097–2120,, 2020. a

Dijkstra, J., Oenema, O., van Groenigen, J., Spek, J., van Vuuren, A., and Bannink, A.: Diet effects on urine composition of cattle and N2O emissions, Animal, 7, 292–302,, 2013. a, b

Faverdin, P., Baratte, C., Delagarde, R., and Peyraud, J. L.: 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 Forage Sci., 66, 29–44,, 2011. a

Feedipedia: Animal Feed Resources Information System, (last access: 14 October 2021), 2020. a, b, c, d, e, f, g, h, i, j

Fernández-Martínez, M., Sardans, J., Chevallier, F., Ciais, P., Obersteiner, M., Vicca, S., Canadell, J. G., Bastos, A., Friedlingstein, P., Sitch, S., Piao, S. L., Janssens, I. A., and Peñuelas, J.: Global trends in carbon sinks and their relationships with CO2 and temperature, Nat. Clim. Change, 9, 73–79,, 2019. a

Friedlingstein, P., Jones, M. W., O'Sullivan, M., Andrew, R. M., Bakker, D. C. E., Hauck, J., Le Quéré, C., Peters, G. P., Peters, W., Pongratz, J., Sitch, S., Canadell, J. G., Ciais, P., Jackson, R. B., Alin, S. R., Anthoni, P., Bates, N. R., Becker, M., Bellouin, N., Bopp, L., Chau, T. T. T., Chevallier, F., Chini, L. P., Cronin, M., Currie, K. I., Decharme, B., Djeutchouang, L. M., Dou, X., Evans, W., Feely, R. A., Feng, L., Gasser, T., Gilfillan, D., Gkritzalis, T., Grassi, G., Gregor, L., Gruber, N., Gürses, Ö., Harris, I., Houghton, R. A., Hurtt, G. C., Iida, Y., Ilyina, T., Luijkx, I. T., Jain, A., Jones, S. D., Kato, E., Kennedy, D., Klein Goldewijk, K., Knauer, J., Korsbakken, J. I., Körtzinger, A., Landschützer, P., Lauvset, S. K., Lefèvre, N., Lienert, S., Liu, J., Marland, G., McGuire, P. C., Melton, J. R., Munro, D. R., Nabel, J. E. M. S., Nakaoka, S.-I., Niwa, Y., Ono, T., Pierrot, D., Poulter, B., Rehder, G., Resplandy, L., Robertson, E., Rödenbeck, C., Rosan, T. M., Schwinger, J., Schwingshackl, C., Séférian, R., Sutton, A. J., Sweeney, C., Tanhua, T., Tans, P. P., Tian, H., Tilbrook, B., Tubiello, F., van der Werf, G. R., Vuichard, N., Wada, C., Wanninkhof, R., Watson, A. J., Willis, D., Wiltshire, A. J., Yuan, W., Yue, C., Yue, X., Zaehle, S., and Zeng, J.: Global Carbon Budget 2021, Earth Syst. Sci. Data, 14, 1917–2005,, 2022. a

Glasser, F., Doreau, M., Maxin, G., and Baumont, R.: Fat and fatty acid content and composition of forages: A meta-analysis, Animal Feed Sci. Tech., 185, 19–34,, 2013. a

He, M., Zhou, G., Yuan, T., Groenigen, K. J., Shao, J., and Zhou, X.: Grazing intensity significantly changes the C : N : P stoichiometry in grassland ecosystems, Global Ecol. Biogeogr., 29, 355–369,, 2020. a

Heinke, J., Lannerstad, M., Gerten, D., Havlík, P., Herrero, M., Notenbaert, A. M. O., Hoff, H., and Müller, C.: Water Use in Global Livestock Production – Opportunities and Constraints for Increasing Water Productivity, Water Resour. Res., 56, e2019WR026995,, 2020. a

Heinke, J., Rolinski, S., and Müller, C.: LPJmL5.0-grazing Model Code, Zenodo [code],, 2022. a

Herrero, M., Fawcett, R., Silveira, V., Busqué, J., Bernués, A., and Dent, J.: Modelling the growth and utilisation of kikuyu grass (Pennisetum clandestinum) under grazing. 1. Model definition and parameterisation, Agr. Syst., 65, 73–97,, 2000. a, b

Herrero, M., Havlik, P., Valin, H., Notenbaert, A., Rufino, M. C., Thornton, P. K., Blummel, M., Weiss, F., Grace, D., and Obersteiner, M.: Biomass use, production, feed efficiencies, and greenhouse gas emissions from global livestock systems, P. Natl. Acad. Sci. USA, 110, 20888–20893,, 2013. a, b

Hodgson, J.: The control of herbage intake in the grazing ruminant, P. Nutr. Soc., 44, 339–346,, 1985. a

Huhtanen, P., Nousiainen, J., Rinne, M., Kytölä, K., and Khalili, H.: Utilization and Partition of Dietary Nitrogen in Dairy Cows Fed Grass Silage-Based Diets, J. Dairy Sci., 91, 3589–3599,, 2008. a, b, c, d, e, f, g

Institut national de la recherche agronomique: Alimentation des bovins, ovins et caprins, Éditions Quæ, Paris, 2007. a, b

IPCC: Emissions from livestock and manure management, in: 2006 IPCC Guidelines for National Greenhouse Gas inventories, pp. 10.1–10.87, Prepared by the National Greenhouse Gas Inventories Programme, Hayama, Japan, (last access: 16 September 2021), 2006. a, b

Johnson, I. and Parsons, A.: A theoretical analysis of grass growth under grazing, J. Theor. Biol., 112, 345–367,, 1985. a, b

Kim, H.: Global Soil Wetness Project Phase 3 Atmospheric Boundary Conditions (Experiment 1), Data Integration and Analysis System (DIAS) [data set],, 2017. a

Klein Goldewijk, K., Beusen, A., Doelman, J., and Stehfest, E.: Anthropogenic land use estimates for the Holocene – HYDE 3.2, Earth Syst. Sci. Data, 9, 927–953,, 2017. a, b

Kondo, M., Patra, P. K., Sitch, S., Friedlingstein, P., Poulter, B., Chevallier, F., Ciais, P., Canadell, J. G., Bastos, A., Lauerwald, R., Calle, L., Ichii, K., Anthoni, P., Arneth, A., Haverd, V., Jain, A. K., Kato, E., Kautz, M., Law, R. M., Lienert, S., Lombardozzi, D., Maki, T., Nakamura, T., Peylin, P., Rödenbeck, C., Zhuravlev, R., Saeki, T., Tian, H., Zhu, D., and Ziehn, T.: State of the science in reconciling top‐down and bottom‐up approaches for terrestrial CO2 budget, Global Change Biol., 26, 1068–1084,, 2020. a

Kozlowski, L. P.: Proteome-pI: Proteome isoelectric point database, Nucleic Acids Res., 45, D1112–D1116,, 2017. a

Lange, S., Mengel, M., Treu, S., and Büchner, M.: ISIMIP3a atmospheric climate input data (v1.0), ISIMIP Repository [data set],, 2022. a, b

Lutz, F., Herzfeld, T., Heinke, J., Rolinski, S., Schaphoff, S., von Bloh, W., Stoorvogel, J. J., and Müller, C.: Simulating the effect of tillage practices with the global ecosystem model LPJmL (version 5.0-tillage), Geosci. Model Dev., 12, 2419–2440,, 2019. a, b

McSherry, M. E. and Ritchie, M. E.: Effects of grazing on grassland soil carbon: a global review, Global Change Biol., 19, 1347–1357,, 2013. a

Moore, J. E., Brant, M. H., Kunkle, W. E., and Hopkins, D. I.: Effects of supplementation on voluntary forage intake, diet digestibility, and animal performance, J. Animal Sci., 77, 122,, 1999. a

National Research Council: Nutrient Requirements of Dairy Cattle: Seventh Revised Edition, 2001, National Academies Press, Washington, DC,, 2001. a, b, c, d, e, f, g, h, i, j, k, l, m, n, o, p, q, r, s, t, u, v

Riedo, M., Grub, A., Rosset, M., and Fuhrer, J.: A pasture simulation model for dry matter production, and fluxes of carbon, nitrogen, water and energy, Ecol. Model., 105, 141–183,, 1998. a

Rolinski, S., Müller, C., Heinke, J., Weindl, I., Biewald, A., Bodirsky, B. L., Bondeau, A., Boons-Prins, E. R., Bouwman, A. F., Leffelaar, P. A., te Roller, J. A., Schaphoff, S., and Thonicke, K.: Modeling vegetation and carbon dynamics of managed grasslands at the global scale with LPJmL 3.6, Geosci. Model Dev., 11, 429–451,, 2018. a, b

Santos, J. and Huber, J.: Feeds, Prediction of Energy and Proteins | Feed Proteins, in: Encyclopedia of Dairy Sciences, 2nd edn., edited by: Fuquay, J. W., Academic Press, San Diego, 409–417,, 2002. a

Tedeschi, L. O., Molle, G., Menendez, H. M., Cannas, A., and Fonseca, M. A.: The assessment of supplementation requirements of grazing ruminants using nutrition models, Translational Animal Science, 3, 811–828,, 2019.  a

Tomé, D., Cordella, C., Dib, O., and Péron, C.: Nitrogen and protein content measurement and nitrogen to protein conversion factors for dairy and soy protein-based foods: a systematic review and modelling analysis, World Health Organization and Food and Agriculture Organization of the United Nations, Geneva, (last access: 25 November 2021), 2019. a

Tulloh, N. M.: Physical studies of the alimentary tract of grazing cattle, New Zeal. J. Agr. Res., 9, 999–1008,, 1966. a

von Bloh, W., Schaphoff, S., Müller, C., Rolinski, S., Waha, K., and Zaehle, S.: Implementing the nitrogen cycle into the dynamic global vegetation, hydrology, and crop growth model LPJmL (version 5.0), Geosci. Model Dev., 11, 2789–2812,, 2018. a

Vuichard, N., Soussana, J.-F., Ciais, P., Viovy, N., Ammann, C., Calanca, P., Clifton-Brown, J., Fuhrer, J., Jones, M., and Martin, C.: Estimating the greenhouse gas fluxes of European grasslands with a process-based model: 1. Model evaluation from in situ measurements, Global Biogeochem. Cycles, 21, 1–14,, 2007. a

Weiss, W.: Feeds, Prediction of Energy and Proteins | Feed Energy, in: Encyclopedia of Dairy Sciences, 2nd edn., edited by: Fuquay, J. W., Academic Press, San Diego, 403–408,, 2011. a

Weiss, W., Conrad, H., and St. Pierre, N.: A theoretically-based model for predicting total digestible nutrient values of forages and concentrates, Animal Feed Sci. Tech., 39, 95–110,, 1992. a, b

Yang, J. and Tian, H.: ISIMIP3b N-deposition input data (v1.0), ISIMIP Repository [data set],, 2020. a, b

Zhou, G., Zhou, X., He, Y., Shao, J., Hu, Z., Liu, R., Zhou, H., and Hosseinibai, S.: Grazing intensity significantly affects belowground carbon and nitrogen cycling in grassland ecosystems: a meta-analysis, Global Change Biol., 23, 1167–1179,, 2017. a

Short summary
We develop a livestock module for the global vegetation model LPJmL5.0 to simulate the impact of grazing dairy cattle on carbon and nitrogen cycles in grasslands. A novelty of the approach is that it accounts for the effect of feed quality on feed uptake and feed utilization by animals. The portioning of dietary nitrogen into milk, feces, and urine shows very good agreement with estimates obtained from animal trials.