Simulation of Long-Term Carbon and Nitrogen Dynamics in Grassland-Based Dairy Farming Systems to Evaluate Mitigation Strategies for Nutrient Losses

Many measures have been proposed to mitigate gaseous emissions and other nutrient losses from agroecosystems, which can have large detrimental effects for the quality of soils, water and air, and contribute to eutrophication and global warming. Due to complexities in farm management, biological interactions and emission measurements, most experiments focus on analysis of short-term effects of isolated mitigation practices. Here we present a model that allows simulating long-term effects at the whole-farm level of combined measures related to grassland management, animal housing and manure handling after excretion, during storage and after field application. The model describes the dynamics of pools of organic carbon and nitrogen (N), and of inorganic N, as affected by farm management in grassland-based dairy systems. We assessed the long-term effects of delayed grass mowing, housing type (cubicle and sloping floor barns, resulting in production of slurry and solid cattle manure, respectively), manure additives, contrasting manure storage methods and irrigation after application of covered manure. Simulations demonstrated that individually applied practices often result in compensatory loss pathways. For instance, methods to reduce ammonia emissions during storage like roofing or covering of manure led to larger losses through ammonia volatilization, nitrate leaching or denitrification after application, unless extra measures like irrigation were used. A strategy of combined management practices of delayed mowing and fertilization with solid cattle manure that is treated with zeolite, stored under an impermeable sheet and irrigated after application was effective to increase soil carbon stocks, increase feed self-sufficiency and reduce losses by ammonia volatilization and soil N losses. Although long-term datasets (>25 years) of farm nutrient dynamics and loss flows are not available to validate the model, the model is firmly based on knowledge of processes and measured effects of individual practices, and allows the integrated exploration of effective emission mitigation strategies.


Introduction
During the last century, in many agricultural systems the inputs of nitrogen (N) bound by the Haber-Bosch process have largely replaced N from sources like symbiotic fixation and mineralization from manures, crop residues and soil organic matter [1]. The large amount of artificial fertilizers used in agroecosystems has resulted in high concentrations of reactive N in the biosphere, which caused negative effects on soil, water and air quality with detrimental consequences for ecosystems, food supply chains and human health [2,3]. Large improvements have been reached in the environmental performance of dairy farming systems through improved farming practices, underpinning research and supporting policies since the 1980s. Nevertheless, effective integrated approaches to reduce these negative effects of agriculture are still urgently needed.
Flows of N on grassland-based dairy and mixed crop-livestock farming systems can be conceptualized as a cycle from soil N uptake by grassland and crops, which are supplied to animals as feed, the ingested feed is partly incorporated into products but the largest proportion is excreted, and the excreta can be used to fertilize the soil [4][5][6][7]. Gaseous emissions and losses to soil and water can occur at various points in the N cycle, and increase when the total amount of N cycling in the system is enhanced by larger inputs [8]. Therefore, many environmental policies have focused on reduction of inputs to decrease the amount of N cycling in the farming system [9,10]. When artificial N inputs are diminished, or even completely abandoned as in organic farming systems, the dependence on natural sources of N increases. Then management should focus more on incorporation of legumes like clovers to fix atmospheric N 2 , cropping and animal housing systems that optimize crop residue and manure utilization, and on slow processes of build-up of organic matter (OM) and N stocks in soils [11]. The interactions among these biological processes are complex and prone to environmental variability, and as a consequence farmers often struggle to develop a coherent new management strategy at lower input levels [10].
At a given level of N cycling, N use efficiency can be increased through mitigation practices. Many policies and practices have been developed to reduce losses from manure management chains for slurry and solid cattle manure (SCM). In the Netherlands, farms with cubicle housing of livestock that produce slurry are obliged to use covered slurry storage facilities and should apply the slurry into the soil to reduce exposure to air and concomitant losses [12]. In contrast, when animals are housed in deep litter or sloping floor barns, a mixture of faeces, urine and bedding material (mostly wheat straw) is produced. For such straw-based systems, Shah et al. [13,14] showed that application of bedding additives like zeolite, farm topsoil and lava meal inside the barn have potential to reduce N losses and enhance N utilization by crops. In another study, Shah et al. [15] reported that anaerobic storage of solid manure followed by 10 mm of irrigation immediately after its application diminished ammonia (NH 3 ) emission rate by 92% while herbage apparent N recovery increased by 33% as compared to non-irrigated manure.
However, identification of appropriate measures to apply on farms in order to reduce losses is complicated because reduction of losses at one point in the N cycle by a mitigating practice are often compensated by higher losses at other points in the cycle [16,17]. Moreover, the consequences of adjustments to farming practices in the long term should be evaluated. To avoid compensatory loss pathways and negative impacts in the long run, a systems-oriented analysis of the whole farm and the N cycle is needed to construct a coherent long-term strategy of mitigation of losses [18,19]. Simulation models can support the evaluation of measures and the development of effective strategies (see [20] for a review of models of livestock systems). Here, we employ an extended version of the Farm DANCES eco-mathematical model [9,21] to evaluate and compare the long-term productive, environmental and economic performance of dairy systems utilizing different manure types and loss mitigation practices. The model simulates the dynamics of organic carbon (C) and N and of inorganic N available for plant uptake on grassland-based dairy farms. It quantifies mineralization and immobilization, NH 3 volatilization and combined soil N losses (through runoff, leaching and denitrification). Farm management decisions on grassland management, the type of manures produced, adjustment of storage methods, and use of low-emission techniques influence these processes at various points of N cycle.
The objectives of the current study were (i) to explore long-term effects of adapting NH 3 mitigation practices such as use of manure bedding additives, contrasting manure storage methods and irrigation after application of covered manure, on the time course of soil organic C and N contents, soil N mineralization, farm productivity and economics, (ii) to compare effects of strategic adjustments of solid cattle manure and slurry-based systems on farm performance, and (iii) to define a coherent strategy combining effective practices, to mitigate losses and to improve farm productive, environmental and economic performance.

Model Description
The model simulates the dynamics of three state variables that quantify the amount of organic carbon (c) and nitrogen (s), and the integrated amount of inorganic nitrogen that is available for plant uptake throughout the year (n). The state variables and the aggregate flows of N and C on the farm are presented in Fig. 1. Hence, the differential equations of the model are: Where I n = inputs of inorganic nitrogen from fertilizers, deposition and fixation (kg ha -1 year -1 ).
M n = mineralization of organic nitrogen (kg ha -1 year -1 ). E n = losses of inorganic nitrogen through NH 3 volatilization and soil N losses by leaching, runoff and denitrification (kg ha -1 year -1 ). U s = net uptake of inorganic nitrogen into organic material by plants, corrected for mineralization from decay of plant biomass, manure and animal digestion (kg ha -1 year -1 ).
I s = inputs of organic nitrogen in feeds (kg ha -1 year -1 ). P s = export of organic nitrogen in crops, manure and animal products (kg ha -1 year -1 ).
I c = inputs of organic carbon in feeds (kg ha -1 year -1 ). U c = net uptake of organic carbon into organic material by plants, corrected for respiration from decay of manure and animal digestion (kg ha -1 year -1 ). P c = export of organic carbon in crops, manure and animal products (kg ha -1 year -1 ).
R c = respiration of organic carbon through decay by soil biota (kg ha -1 year -1 ).
The model is target-oriented, based on a production level of milk and meat that is defined by the size and productivity of the herd (see Table S1 for parameter values). Energy and protein requirements were calculated on the basis of the Dutch feed evaluation systems [22,23]. The animals are fed with on-farm produced grass and feed crop products (in this case silage maize), and supplementary feed is imported when the amount of feeds produced on the farm is insufficient to cover the energy and protein requirements of the herd. If there is a surplus of on-farm produced feed crop export occurs.
Grassland production is described by the response of N uptake (U) to available inorganic nitrogen (n), and the relation between U and biomass yield (Y) (Fig. 2). These relations are defined by adjusted expo-linear equations [9,24]. N is taken up in harvestable and unharvested biomass, because farm animals can only harvest part of the total amount of plant biomass produced, the remainder staying behind in the field as organic material. Therefore, we distinguished total and harvested amounts of N uptake (U T and U H ) and total and harvested biomass (Y T and Y H ), which resulted in four equations. Equation (4) shows the general form of the expolinear equation.
Where y = the dependent variable, representing U T , Y T , U H or Y H (kg ha -1 year -1 ).
x = the independent variable, representing n, U T , or U H (kg ha -1 year -1 ).
y MAX = maximum value of y, representing U MAX,T , Y MAX,T , U MAX,H or Y MAX,H (kg ha -1 year -1 ). r = initial response of y to x, representing r U,T , r Y,T , r U,H and r Y,H (kg kg -1 ha -1 year -1 ). l = the decline of the response of y to x (kg -1 ).
The initial response of N uptake to available inorganic N and of biomass production to N uptake is linear, with an initial slope r. This initial slope declines with a rate depending on l until the maximum y MAX is reached. The ratio between U MAX,H and U MAX,T is denoted h N : the fraction of harvested N in biomass. Maximum dry matter yield is calculated as Y MAX = Y MAX /a MAX , where a MAX is the maximum N content of grass. The initial response of both total and harvested plant biomass yield to N uptake is calculated from the minimum N content in herbage: r Y = 1/a MIN , with a MIN as the minimum N content of grass. The grassland production curves were calibrated for mixed use by mowing and grazing. For silage maize a yield level Y MAIZE is defined. The parameters for maize production and the grassland production curves can be found in Table S2.
The harvested biomass is fed to animals and can be partly exported in case of a feed surplus. The feed is partly digested by the animals (k D , year -1 ) and the undigested fraction enters the manure, where it is subjected to further degradation (k M , year -1 ) during storage before application to the field. Manure degradability is dependent on the feed quality, therefore the k M is proportional to k D using constant g M , so that k M = k D Ng M . The parameters regarding feed quality are in Table S3.
Soil processes include degradation and additions of OM, and soil N losses. Mineralization results in a decline of the organic carbon pool s, with a fractional rate k S (year -1 ). Moreover, a fraction of the unharvested biomass is degraded in the year of production (k B , year -1 ). The organic carbon in manure and unharvested biomass that remains undegraded after the year or production is added to s. A proportion of available inorganic N is converted to biomass and lost through volatilization, leaching, denitrification and runoff, resulting in a relative rate of withdrawal of inorganic nitrogen (k W , year -1 ). Nitrogen soil losses are calculated as the difference between plant uptake and total withdrawal (Fig. 2). The withdrawal fraction is assumed to be 95% of the available inorganic nitrogen n [9], which is applicable on well-drained soils in temperate regions where winter precipitation exceeds evapotranspiration resulting in complete loss of nitrate N [25]. Only small residues of inorganic N in NH 4 + persist in winter when not subject to nitrification [26]. Parameters of soil processes are presented in Table S4.
Nitrogen in animal excreta is present in organic and inorganic forms. An overview of the processes and conversions of N in manure is presented in Fig. 3. The calculations are largely based on calculation procedures presented by Dä mmgen and Hutchings [27]; we followed the same steps in the calculations but adjusted the calculation of mineralization of manure organic N (see below). Part of the inorganic N can be adsorbed to straw and manure additives. The inorganic N is prone to emission by NH 3 volatilization after excretion in the barn (f E , g g -1 ), during storage (f S , g g -1 ) and after application to the field (f A , g g -1 ). The loss fractions f E , f S and f A are dependent on the barn and storage conditions, and the method of manure application and extra emission mitigating measures after application such as irrigation or application during rainfall. The values of manure parameters as used in the model are presented in Table S5.
To estimate the mineralization or immobilization of N due to degradation of OM by microorganisms we use the following equation: Where M = net mineralization (kg year -1 ). k = fractional degradation rate of the substrate (year -1 ). C S = amount of carbon in the substrate (kg C). e = growth efficiency of the microorganisms (kg kg -1 ). q S = C:N ratio of the substrate (kg C kg -1 N). q M = C:N ratio of the microorganisms (kg C kg -1 N). Microorganisms break down the OM with apparent fractional degradation rate k. However, because microorganisms grow due to this degradation process with growth efficiency e, this results in an addition to the OM, so the observed degradation rate should be corrected for their growth efficiency, and the true degradation rate is C S Nk/(1-e). The degradation of the OM is associated with release of N, determined by the C:N ratio of the substrate (q S ). The micro-organisms will incorporate part or all of this N, dependent on their C:N ratio (q M ). When the C:N ratio of the substrate is high, the release of N from OM is lower than the incorporation into microbial biomass, and as a consequence mineral N from the soil solution may be incorporated. This results in negative values for N release in the equation above, indicating net immobilization. Mineralization occurs when the C:N ratio of the substrate is lower than q S /e. The gross margin (revenues minus variable costs) was used as indicator for the economic performance of the farm. Gross margin is considered the most appropriate economic indicator in a fixed milk quota system [28] and is more sensitive to changes in farm management than total farm results, which also include fixed costs [29]. The revenues included returns from milk, meat and crop sales and other sources, and the variable costs were made for purchases of feeds and bedding material, veterinary care, breeding association and miscellaneous costs (Table S6).

Case Study Farm and Scenarios
We defined a typical grassland-based organic dairy farm on a sandy soil in the province of Gelderland, The Netherlands. The farm area is 66 ha, of which 60 ha is grassland and 6 ha is used for cultivation of silage maize. The herd consists of 85 Holstein Frisian cows with a replacement rate of 25%. The herd is housed in a sloping floor barn, wherein a mixture of faeces, urine and bedding material (mostly wheat straw; SCM) is produced. The SCM is removed regularly from the barn and is stockpiled for storage. The animals graze during a period of 200 days per year (from mid April until the end of October) in a day-and-night grazing system wherein the cows spend 20 hours per day outdoors. The mown grass is conserved as silage and fed to cows indoors during the 165 days winter season.
We developed scenarios to evaluate the effects of changes in parameters due to adjustment in individual farm management practices regarding grassland and manure management, and a more integrated strategy combining various effective measures. With these scenarios we evaluated long-term productive, environmental and economic farm performance as affected by: -Different animal housing systems that produce either SCM (faeces and urine mixed with wheat straw) from the sloping floor barn that is standard for the farm (scenario M), or slurry (mixed faeces and urine) from a cubicle housing system (scenario S). The latter results in lower fractional loss rates during storage and after application (f S and f A ) and higher yield of silage maize (Y MAIZE ). -Delayed mowing of grass (scenarios DS and DM) resulting in harvesting of more mature grass with lower feed quality (k D,GRASS ) and N content (lower a MIN,H , a MAX,H and a MAX,T ) [9,30,31]. Due to the proportionality between feed quality and manure degradability, also k M will decline in these scenarios. -The use of the additives zeolite, lava meal and farm topsoil that are applied on SCM bedding inside the barn (scenarios MZ, ML and MT). The impacts of these additives on N losses after excretion in the barn, during storage and after application to the field and on N uptake and dry matter yield of grass and silage maize (at physiological maturity) have been quantified by [13,14]. Emission factors f E , f S and f A , grassland production ( Fig. 2) and maize yields (Y MAIZE ) were derived from these experimental results. Costs for additives were included in the price of bedding material. -Alternative SCM storage systems of composting, roofing or covering by an impermeable sheet (scenarios MC, MR and MU). Shah et al. [32] quantified the consequences of these measures for emissions during storage (f S ). -Combining covering of the manure by an impermeable sheet (anaerobic storage) with 10 mm irrigation immediately after manure application to the field (scenario MUI). These combined measures affect NH 3 volatilization after application (f A ) and N recovery by the grassland (Fig. 2) as analysed experimentally by [15]. -A combination of measures of SCM handling that appeared most promising for productive, environmental and economic indicators from the previously described scenarios. This integrated strategy contained practices of delayed mowing,  zeolite additive, covering with an impermeable sheet and irrigation after application (scenario DMZUI).
A complete overview of the parameters settings for all scenarios is provided in Table S7.
For the starting conditions for all scenarios we assumed the current situation on the farm, which has been under the management described above that is comparable to scenario M for almost 20 years, hence there is still build-up of soil organic matter and a steady state has not been reached. The long-term dynamics of the state variables s, c and n were evaluated for each of the scenarios for simulation duration of 200 years assuming constant farm management. Moreover, the productivity of the farm under the different scenarios was determined with the feed self-supply rate, which at the target level of animal outputs reflects the production of on-farm feeds (grass and silage maize). The farm gate N balance (inputs minus outputs in products) reflects the total farm N losses and was used as an indicator for environmental performance. The gross margin was the indicator of economic performance for each scenario.

Slurry and SCM-based Systems
The slurry and SCM-based systems without mitigation measures in scenarios M and S contrasted strongly in C and N dynamics (Figs. 4a-4e). Initially the slurry-based system resulted in a larger amount of available inorganic N than the SCM-based system (Fig. 4c), which could support larger grassland productivity. However, the slurry-based system resulted in slightly declining soil organic C and N pools, whereas for the SCM system these pools gradually increased (Figs. 4a and 4b) due to the straw inputs for bedding. As a consequence, after 75 years of simulated management, the inorganic N availability was larger for the SCM system than for the slurry system (Fig. 4c), due to increased mineralization of the large organic N pool. This resulted in increased grassland production and a higher feed self-supply rate for the SCM-based system ( Table 1). The NH 3 emissions were higher from the SCM system (Fig. 4d). N soil losses were strongly linked to available inorganic N, therefore these soil losses from SCM system were initially lower than from the slurry-based system, but were larger than for the slurry system after 75 years (Fig. 4e). Thus, in the equilibrium situation total N losses (NH 3 volatilization and N soil losses) were higher from the SCM-based system than from the slurry-based system (62 vs. 55 kg N ha -1 year -1 ).

Effects of Delayed Mowing
Later mowing of grass results in the on-farm production of more grassland biomass with a lower N content in scenarios DM and DS. This led to reductions of N volatilization of ca. 4.5 kg N ha -1 year -1 (Fig. 4d), while N soil losses also declined slightly (Fig. 4e). In the long term the simulated accumulation of soil organic C and N was larger when more mature grass of lower degradability was fed, both in slurry and SCM-based systems (Figs. 4a and 4b).

Effects of Storage Measures and Irrigation
The storage treatments for SCM-based systems of roofing (scenario MR) or sealing with an impermeable sheet (MU) of the manure helped to reduce NH 3 volatilization losses, although the reductions of losses during storage were largely compensated by extra emissions after application, in particular for the sealing treatment (Fig. 4i). The treatments had no effects on soil organic C and N dynamics (Figs. 4f and 4g), feed self-supply, whole farm Nefficiency and gross margin (Table 1). However, when sealing was combined with irrigation in scenario MUI, the application losses were avoided, and N volatilization could be reduced to only 6 kg N ha -1 year -1 . As a compensation loss for lower volatilization, the combined treatment of sealing and irrigation caused higher N soil losses (Fig. 4j). Nevertheless, the inorganic N availability was larger due to sealing and irrigation (Fig. 4 h), resulting in better grassland production so that the productive, environmental and economic indicators of feed self-supply rate, farm N-efficiency and gross margin of scenario MUI were better than untreated SCM or only storage measures ( Table 1).

Effects of Bedding Additives
The addition of farm topsoil, zeolite or lava meal to the manure bedding (scenarios MT, MZ and ML) resulted in reduced emission through volatilization, which was proportionally but only partly  compensated by more N soil losses, so that inorganic N availability was higher (Figs. 4m-4o). These effects were strongest for the zeolite additive. The three additives resulted in similar improvements of feed self-supply rate, farm N-efficiency and gross margin when compared to the SCM-based system without treatments (Table 1). There were no differences among the additives on the long-term soil C and N dynamics (Figs. 4k, 4l).

Combined Effects in a Coherent Strategy
A strategy of emission mitigation was defined by combining the most successful practices for management of SCM: delayed mowing of grassland, using zeolite as bedding additive, storage under an impermeable sheet, and irrigation after application (scenario DMZUI). This strategy reduced NH 3 volatilization with limited compensatory N soil losses, so that total losses were reduced to 40 kg N ha -1 year -1 and inorganic N availability was enhanced (Figs. 4k-4o; Fig. 5). Also the long-term increments in soil C and N were larger than for untreated SCM due to the lower degradability of mature grass after delayed mowing. Both the feed self-supply rate and the whole farm N efficiency were superior to all alternative systems, and the gross margin was comparable to that of the slurry-based systems (Table 1).

Discussion
SCM-based systems are often associated with larger N losses than slurry-based systems [33]. Many of these losses seem to be avoidable in the short term through appropriate management practices for manure after excretion in the barn and during storage as observed in experiments. These practices could be evaluated in the model simulations for their impacts on N losses and soil C and N pools in the long term: -Delayed mowing results in a higher C:N ratio in the feed, and more organic C and less NH3 in the manure [10,21]. The high C:N ratio of the manure OM may cause immobilization of N upon application to soil, followed by a slow rate of mineralization [34]. Simulation results revealed that this practice contributes to increased soil OM build-up, but has no long-term effects on N losses. -The application of additives like zeolite, farm soil or lava meal to the bedding material results in increased NH 4 + adsorption, which reduces NH 3 volatilization [35][36][37]. The model outputs showed that this will result in some compensatory soil losses, but soil N availability will improve in a long run. Soil organic C and N pools were not affected when compared to SCM system without any treatment.
-Covering manure heaps with an impermeable sheet creates a physical barrier that avoids exposure to air and prevents NH 3 diffusion to the atmosphere [38][39][40]. This only affects the storage phase, but after application to the field the simulated volatilization losses were higher. -Moreover, covering results in anaerobic conditions that slow down OM degradation during storage [34,41], but the labile OM will be rapidly degraded after application to the field. Thus no effects on soil organic C and N pools can be expected, as reflected in the simulation results.
Thus, the simulations demonstrated that individual mitigation measures to reduce losses resulted in compensatory loss pathways. Moreover, in contrast to slurry, SCM cannot be injected in grasslands, while shallow injection of cattle slurry can reduce NH 3 volatilization by up to 74% [42,43]. Conserving inorganic N in the manure during the housing and storage phases leads to higher concentrations in applied manure, and can result in increased emission rates during surface application. Consequently, for SCM additional measures like irrigation or application shortly before rainfall are needed after application to the field to enhance infiltration of total ammoniacal N into the soil [15]. Therefore, for effective mitigation of N losses at the farming systems level and in the long term, a strategy composed of a series of techniques would be needed to address the various potential outflows of N from the system. The model simulations for scenario DMZUI demonstrated that combined management practices of applying zeolite on the SCM bedding inside the barn, anaerobic storage of this manure under impermeable plastic sheet, 10 mm of irrigation immediately after surface application of the manure on grassland and its delayed mowing is the most effective combination to increase soil C and N stocks and to reduce N losses. Such a strategy could result in lower losses, higher productivity and similar economic results as slurry-based systems.
A large advantage of SCM-based systems in the long term is the increased soil organic C and organic N contents (Figs. 4a and 4b) as compared to slurry-based systems due to larger inputs of OM. For SCM scenarios, annual rate of increase in soil OM was greatest in the early phases of the simulation and very low near the end as the soil approached an equilibrium state. This is in agreement with findings in a long-term simulation study [44] and experimental data [45,46]. Besides the contribution to C sequestration, increasing the OM content of soils is important for physical and biological soil properties and processes that support many ecosystem functions. OM contributes to the water holding capacity, cation exchange capacity and infiltration capacity of soils [47,48]. Moreover, there is in general a positive relationship between soil C content and soil microbial biomass [49], and any practice that increases the amount of soil OM improves its biological activity, e.g. [50]. These biota activities can enhance mineralization of soil OM, and hence the supply of inorganic N for plant growth.
In modelling we search for a balance between the level of detail, the precision required, the model's flexibility and the data requirements [51][52][53]. The Farm DANCES model used in this study can be characterized as an eco-mathematical summary model that quantifies the dynamics of organic N and C, and inorganic N, as an instrument to evaluate management strategies. By combining all the relevant processes in the farm N cycle, the model allows to assess interactions among these processes and to identify emergent system properties such as compensatory loss pathways. It offers a quantitative framework for evaluating both short-term and long-term effects of management interventions aimed at improving nutrient use efficiency [54]. This framework supports formulation of scenarios describing future developments, rather than exact prediction (cf. [55]). We are not aware of any empirical data that would enable validation of the whole farm model over a substantial period of time of 25 years or more. Therefore, model validity must be inferred from validity of its components and the plausibility of its results. The model constitutes a complement to studies that emphasize short-term optimization of performance of farm system components, such as emission from barns or N leaching at given soil management, and studies that focus on empirical relations between production factors, such as fertilizer and outputs [9].
From a model user's perspective, the focus of the Farm DANCES model is on quantifying interactions among farm components and biological processes and to provide insight into these interactions to its users, which are currently predominantly researchers and students. The model builds on existing knowledge of biological processes, is data-sparse, can be parameterized with experimental data, and the graphical user interface of the model is intuitive and easy to use. There is scope to improve the model by adding flexibility to the scenarios over time, by including the impact of varying environmental conditions (e.g. temperature and moisture dependence of degradation processes), by further specification of N soil loss pathways and greenhouse gas emissions (e.g. [56]), and by compartmentalizing soil organic matter pools with distinctive degradation dynamics (e.g. [57]). However, these extensions would sacrifice the insightfulness, whereas various modelling studies have demonstrated that relatively simple dynamic models that are based on the correct process representation and data can be extremely accurate and useful (e.g., for soil processes: [58][59][60]; for plant growth: [61,62]; review for livestock systems: [20]).

Conclusions
The simulation results demonstrated that individual emission mitigation measures were often insufficient to reduce N losses at the farming systems level. Practices that reduced NH 3 emissions from animal excreta in the barn or during storage resulted in larger losses after application of manure to the field, through either volatilization of NH 3 or soil losses, i.e. the aggregated flows of runoff, leaching and denitrification. The integrated strategy combining the most effective practices resulted in build-up of soil organic C and N pools, sufficient nutrient availability for plants and low emission rates. This strategy of combined grassland and manure management practices included delayed mowing of grass and fertilization with solid cattle manure that is treated with zeolite, stored under an impermeable sheet and irrigated after application. This strategy can reduce losses to the environment, improve soil properties by larger organic C and N stocks, and increase availability of N for plants uptake, grassland productivity, enhance the feed self-supply. We conclude that SCM-based systems employing a coherent strategy of manure utilization practices can contribute to improved productive, environmental and economic performance of dairy farming systems.

Author Contributions
Conceived and designed the experiments: GAS JCJG EAL. Performed the experiments: GAS JCJG. Analyzed the data: GAS JCJG GMS EAL. Contributed reagents/materials/analysis tools: GAS JCJG GMS EAL. Wrote the paper: GAS JCJG GMS EAL. Developed the simulation model used in analysis: JCJG.