Received 8 June 2016; accepted 13 August 2016; published 16 August 2016
Agricultural activity is estimated to be responsible for approximately 14% of global anthropogenic greenhouse gas (GHG) emissions  . In the European Union (EU), agriculture comprises 10% of emissions, with CH4 and N2O contributing 49% and 63%, respectively to sectoral emissions  . In the Republic of Ireland (ROI), agricultural emissions comprise one-third of national emissions and remain a key component of national emissions despite recent decreases, due to the economic downturn  . Agricultural GHGs, particularly N2O, are produced mainly through biological processes and the degree of variation (spatial and temporal) in emissions depends on soil type, land use and climatic factors (e.g.  -  . Agricultural soils may either be net sinks or sources, depending largely on the balance between N2O release and carbon (C) sink-source strength and functional relationships exist between organic C and N, derived from either inorganic or organic sources, to produce GHGs. These relationships are regulated by agricultural activities and associated disturbances (e.g.   . Thus, an understanding of the associated controlling factors and their interactions, including impact of site-specific soil conditions are key requirements  for understanding the GHG balance and for the development and selection of an appropriate model for accounting and reporting.
Most of the Annex-I countries are using IPCC Tier 1 methodologies    for the estimation of agricultural GHGs due to a lack of detailed, spatially-explicit activity data and the absence of disaggregated emission factors (EFs). Some countries (e.g. New Zealand, USA) have moved to Tier 2, with country-specific emission factors and are developing Tier 3 (modelling) methodologies. The Tier 1 approach has several limitations for studies of the GHG balance relevant to Agriculture, Forestry and Other Land Uses (AFOLU)/LULUCF  . Development of higher tiers requires good country/regional-specific activity data allied to extensive GHG emissions datasets. Compared to Tier 2, more additional resources are required for the development of Tier 3, including an appropriate biogeochemical model. A process-based model could take into account functional relationships and provide a flexible and structured way to assess how different scenarios including land-use management and land-use change can affect GHG emissions and soil C and N dynamics. A modelling approach can provide improved estimates of GHG budgets and reflect more robust emissions assessments (sink or source) by reducing the uncertainties associated with the impacts of soil, climate and management activities. The advantages in using a model include an ability to 1) scale GHG emissions from the site-specific to the national/regional level, 2) identify potential mitigation options and the interactions between different gaseous and/or other loss pathways, and 3) provide a better understanding of how agricultural soils can act as C sinks or sources.
In line with commitments under the UNFCCC, the ROI is committed to improving the estimation of GHG budgets by developing Tier 3 approaches. There has been much progress in recent years in developing models to simulate GHG emissions. Modelling is considered a low-cost method of estimating GHG emissions from agricultural soils at different scales and for exploring potential mitigation strategies  . Any discrepancies in predicting field measured GHG fluxes can be used to identify limitations to model use, suggesting further potential developments that could lead to a better representation of the impacts of differences in land use, land use management or environmental factors. However, spatial variability in GHG fluxes due to soil heterogeneity, in particular, has to be considered when comparing model results with field measurements   . There is also a particular need for improved methodologies for up-scaling of GHG emissions from site to regional/national scales  . The use of models for estimating soil GHG emissions is expected to increase, and further improvements in model accuracy and precision will be essential.
Several process-based models are currently used to predict a variety of variables related to different ecosystems. DNDC (DeNitrification-DeComposition) is a process-based model that simulates carbon (C) and nitrogen (N) biogeochemistry in agro-ecosystems and has been used for predicting GHG emissions, soil C dynamics, crop growth and other relevant data  . DailyDayCent is the daily time-step version of the CENTURY biogeochemical model  , which simulates daily N-gas fluxes, C fluxes and other ecosystem parameters  . More recently, the Rothamsted Carbon Model (RothC) and SUNDIAL (SimUlation of Nitrogen Dynamics in Arable Land) have been used in the development of a multi-pool, process-based model, called “ECOSSE” (Estimating Carbon in Organic Soils-Sequestration and Emissions)  . The ECOSSE model simulates soil C and N turnover, including GHG/trace gas emissions, in both mineral and organic soils using only limited meteorological, land-use and soil data, compared to other models. In the ROI, some models (DNDC, DayCent, RothC, PASIM, etc.) have been tested/validated using limited datasets measured from grassland and arable systems  -  . However the results have not been sufficiently robust to allow them to be used in the inventory process without further investigations.
Based on the different characteristics and performances, three process-based models (DNDC v9.4, DailyDayCent and ECOSSE v5+) were chosen to evaluate GHG emissions associated with the major Irish cropland type. The goal was to establish the basis of an emission inventory system using process-based models with the minimum number of commonly available input parameters that reflect the site-specific diversity of management practices that influence GHG emissions. Barley, with dominancy of spring barley, is the major cereal crop in Ireland, comprising 71% of the total cereals in 2014  . Multi-year GHG fluxes data for spring barley, measured at the plot scale, were available to initiate model comparison exercises. The main objectives were: 1) to simulate daily N2O, CH4 and CO2 emissions from conventionally-tilled spring barley fields located in Carlow, Ireland over 8 years using the DNDC v9.4, DailyDayCent and ECOSSE v5+ models; 2) to assess the extent of statistical agreements (R2, RMSE, RE and MD) between model outputs and measured datasets; and 3) to evaluate the differences between the measured and modelled seasonal/annual GHG emissions and their estimated EFs.
2. Materials and Methods
2.1. Experimental Sites and Datasets
Data on inputs and management practices were collected from plot-scale field experiments conducted at the Teagasc Oak Park Research Centre, Carlow (52°86' N and 6°54' W). The soil (0 - 10 cm depth) at Oak Park site is classified as a sandy loam (overlying loam) in texture, free draining, Euteric Cambisol (Grey Brown Podzolics). Detailed site characteristics, which may differ from other published information    due to averaging of samples taken from both small and large plots, are given in Table 1. Thirty years (1982-2011) climate data measured from the nearby weather stations (Oak Park, Carlow and Kilkenny, 30 km away) by Met Eireann were used as inputs to run the models. The meteorological inputs required for DNDC and DailyDayCent are daily minimum and maximum air temperature, rainfall, wind speed, radiation and relative humidity; and for ECOSSE only average air temperature, rainfall and potential evapotranspiration is required.
2.2. Description of Field Experiments
The larger plots used (2.5 ha) for field-scale studies were part of an experiment comparing the effects of conventional and minimum tillage practices    -  . An experiment was also carried out using small plots (6 m × 25 m = 150 m2, each plot containing a 0.27 m2 ground area chamber), which were on the border of the large plots. Only the results from the conventional (CT) tillage treatment were used in the present paper. The CT plots were prepared using a mouldboard plough to a depth of 22 - 25 cm. Subsequently, a light tilling was performed in the CT treatment and seeds (spring barley, cv. Tavern/Quench) were sown using a cultivator drill, followed by rolling. The small plots comprised three randomized. Each main plot divided into two subplots containing different fertilizer treatments: fertilized (N1) and non-fertilized (N0) plots. Treatments were randomly distributed and each treatment was replicated four times.
Following the harvesting of the crop (July or August), crop residues were chopped and left on the field over the autumn and winter period (Table 1). The experiments were conducted under rainfed conditions. N fertilizer was applied in the form of calcium ammonium nitrate (CAN). The amount of N applied varied slightly from year to year. From 2004-2006 fertilizer was applied at 0 and 135 - 159 kg∙N∙ha−1, whilst from 2007-2011 it was applied at 0 and 135 kg∙N∙ha−1. It was split into two applications from 2005 onwards. The unfertilized control started in 2003 and prior to that the whole field had received annually 140 - 160 kg∙N∙ha−1 and had been in spring barley production since the mid-1990s. Herbicide (glyphosate 3 l∙ha−1 of a 360 g product) was applied in January or early February to control over-wintering weeds and volunteer barley seedlings.
Table 1. Site characteristics of experimental field as well as inputs and management practices (EC = Eddy Covariance/large plot received highest N rate).
2.3. Measurements of N2O, CO2 and CH4
Measurements of GHGs from the experimental plots were either made seasonally or annually and for three years commencing from 2009 to 2011, at daily or fortnightly intervals. N2O emissions were measured using the static closed chamber method. Gas was sampled at 0, 30 and 60 min intervals between 9 and 11 am every week and more intensively (twice weekly) following fertilizer application. The gas samples were stored in exetainers (Labco, High Wycombe, UK) prior to the analyses. The gas analyses were performed using a gas chromatography (Varian CP 3800 GC, Varian, USA) fitted with a 63Ni electron capture detector (ECD) for N2O analysis and a Flame Ionisation Detector (FID) for CH4 analysis with high purity helium as a carrier gas. Samples were returned to ambient pressure prior to analysis and fed into the system by a Combi-Pal automatic sampler (CTC Analysis, Switzerland). Following a two-year gap, gas samples for the measurement of both N2O and CH4 were collected from September 2008 to September 2010 and from April 2009 to September 2010, respectively. Gas sampling was carried out weekly during the crop growth period and less frequently (2 - 3 weeks) during the fallow period using static chambers, with 18 replicates.
Eddy Covariance (EC) systems installed in the large plots, consisted of Gill R3 sonic anemometer (Gill Instruments, USA) and Li-7000 infra-red gas analyser (Licor Inc., USA), for net ecosystem exchange (NEE) and ecosystem respiration (Reco) measurements. Estimates of Reco (2003-2007) were based on half-hourly measurements and expressed on a daily basis.
2.4. Determination of Soil Nitrate Concentrations
Soils were sampled during the gas sampling periods and soil nitrate concentrations were determined on 2M KCl extracts using an auto-analyzer (Bran and Luebbe, Norderstedt, Germany)  .
2.5. Description of Models
Three dynamic models (ECOSSE v5 updated in 2012, DNDC v9.4 and DailyDayCent) were selected for this comparative study. Input requirements for each model differ as indicated previously. However, the site characteristics and crop management practices used were the same for all the models which were run for 8 years. Other inputs were either defaults or module-based. A brief description of the models is given below. ECOSSE was mainly calibrated under UK conditions  with subsequent improvements of the N2O and CH4 modules using Irish data  . Both DNDC and the DailyDayCent were calibrated/validated under Irish conditions    .
2.5.1. ECOSSE Model
The ECOSSE model was developed to simulate SOC in highly organic soils from algorithms originally derived for mineral soils in the RothC and SUNDIAL models   . The model uses a pool type approach, and all of the major processes of C and N turnover in the soil are included and are driven by readily available input variables (e.g., SOC, soil water, plant inputs, nutrient applications and timing of management operations). It is a tool for site-specific simulations that apparently does not result in any major loss in accuracy at this scale and makes full use of the limited information that is available to run models whilst still providing accurate simulations of GHGs. The N2O fluxes derive from both nitrification and denitrification, CO2 corresponds to RH and CH4 through a balance between methanogenesis and methanotrophy and changes in SOC stocks. The model can be used with both organic and mineral soils, to provide estimates of the net change in soil C and N in response to changes in land use and climate. This model considers variations for outputs by calculating them on each soil layer for each time step. This model doesn’t use crop growth parameters as inputs but uses a built-in default functional relation.
2.5.2. DNDC Model
The DNDC is a widely used process-based model   , but several modifications/versions have been developed for different production systems. This model couples denitrification and decomposition processes to predict emissions of C, with CH4 oxidation, and N from agricultural soils that are governed by various soil and environmental factors. It contains six sub models: soil climate, crop growth, decomposition, denitrification, nitrification and fermentation, and includes subroutines for cropping practices (fertilization, irrigation, tillage, crop rotation and manure addition) to simulate SOM turnover. The model considers decomposition process as first order kinetics, and the soil is considered as a series of discrete horizontal layers with uniform soil properties within each layer, except for some soil physical properties that are anticipated as being constant across all layers. However, time-dependent variations in soil moisture, temperature, pH, C and N pools are considered for a reliable estimate of C and N fluxes by calculating them for each soil layer for each time step.
2.5.3. DailyDayCent Model
DailyDayCent is a biogeochemical model based on the Century soil C model and, for the most part, the parameter files used are identical to the ones used by Century 4.5 and DayCent 4.5   . This model simulates C and N fluxes between the atmosphere, vegetation, and soil. Major factors (e.g. nutrient availability, water, temperature) controlling plant growth are included in order to simulate GHGs and SOC changes over time. This model considers nutrient supply as a function of SOM decomposition and external nutrient additions. Other model inputs are the timing and description of management events (e.g. fertilization, tillage, harvest), and soil texture data. There are submodels that consider plant production, SOM decomposition, soil water and temperature for each layer, as well as nitrification and denitrification, and CH4 oxidation. Improvements in this model are on-going, and comparison of model results and plot data have shown that DayCent reliably simulates crop yield, SOM levels, and trace gas fluxes for various native and managed systems  .
2.6. Model Run, Statistical Evaluation and Calculation
The datasets were collated and compiled to prepare a list of common input parameters with respect to site characteristics and managements to initialize the models (Table 1). The models were run using the common input parameters and the weather data (data not shown). For DNDC and DailyDaycent, the corresponding simulation spin-up for 30 and 700 years were used to allow the model to reach equilibrium state. For ECOSSE, the soil C pool at steady state equilibrium, with crop residue inputs and N as the concentration measured immediately before the start of experiments was used for the initialization.
The ECOSSE model can predict soil heterotrophic respiration (RH) only whereas the EC provides Reco (soil autotrophic and heterotrophic respiration + crop respiration). For comparison and validation of ECOSSE-simu- lated RH with measured ones, daily Reco measured by EC from the large fertilized plot was transformed to daily RH using DailyDayCent fractions (RH/Reco) obtained from this study. Calculation of the total/cumulative N2O, RH and Reco through integration of the measurement values and the sum of simulated values were performed. Seasonal and annual emission factors (EFs) for N2O over the 8 years were calculated by subtracting cumulative measured and model outputs of the unfertilized control from that of the fertilized treatments and dividing by the respective N inputs.
The outputs were collated and converted into standard units for comparison with measured datasets. The simulated values of GHGs were compared and validated quantitatively with measured values using MS Excel speadsheet (MODEVAL v 2.0)  . Based on the available measured flux data, an evaluation of the consistency of seasonal (4 - 5 months)/annual N and C emissions with simulated values was carried out. The approach was to take a simple mean and standard error (SE) of the values for each dataset and to calculate statistics (e.g., R2 = Coefficient of determination; RMSE = Root Mean Square Error; RE = Relative Error; MD = Mean Difference) that describe the model fits for all data points simulated by placing equal weight on all values. An analysis of variance for significance at the 0.05 level of probability was performed and the 95% confidence intervals calculated using SAS v. 9.3 (SAS Inc.), MODEVAL and MS Excel (v. 2010).
3.1. Simulated and Measured Nitrate-N Concentrations
The measured surface soil concentration was found to reach a maximum of 71.3 kg∙N∙ha−1 following fertilization, decreasing to 0.82 kg∙N∙ha−1 during later periods (Figure 1). In the unfertilized field, the minimum and maximum surface levels (seasonal/annual) measured was 0.20 and 25.2 kg∙N∙ha−1, respectively. The DailyDayCent and DNDC-predicted soil levels were highly variable. For the fertilized field, the average contents predicted over the 8 years by DailyDayCent (71 ± 1.6 kg∙N∙ha−1) and DNDC (193 ± 4.6 kg∙N∙ha−1) were markedly higher than the ECOSSE simulated (20 ± 0.4 kg∙N∙ha−1) values or the measured ones. For the unfertilized field, the corresponding average surface concentrations predicted were 12 ± 0.2, 17 ± 0.4 and 8 ± 0.1 kg∙N∙ha−1. Only the ECOSSE model predicted values that were closer to measured values, which were con- sistent over the 8 years, and to the amount of applied. For the fertilized field, the ECOSSE (R2 = 0.55)
Figure 1. Comparison of field measured (seasonal/annual; open circle/square with vertical bars as standard errors) nitrate-N concentrations (kg∙N∙ha−1) in the 0 - 10 cm soil depth with values simulated (solid line) using the three process-based models over 8 years, commencing from 17 August 2003 (day of harvest), in conventionally-tilled fertilized arable land cropped to spring barley. Solid arrows indicate the days of ploughing and dotted arrows indicate the days of N fertilizer application.
and the DNDC (R2 = 0.31) model estimates correlated significantly (p < 0.05) with the measured values (Table 2). None of the models simulated the surface soil concentrations for the unfertilized field accurately, showing poor coefficients of determination (R2 = −0.07 to 0.13). The total error and bias differences between the simulated and measured soil concentrations were large and significantly (p < 0.05) greater than their 95% confidence intervals for both fields.
3.2. Performance of Models to Simulate GHG Emissions
3.2.1. N2O Emissions
The maximum N2O flux measured across all years was observed in 2004 (56.0 g∙N∙ha−1∙d−1) (Figure 2). For the other years the maximum value was 17.6 g∙N∙ha−1∙d−1 and the minimum −8.0 g∙N∙ha−1∙d−1, demonstrating small differences with the unfertilized plot (16.6 versus −10.4 g∙N∙ha−1∙d−1). Regardless of the models, the simulated N2O fluxes were consistent over the years but differed from the measured values, and none of the models predicted fluxes less than zero. The N2O fluxes varied largely between the fertilized (80.0 - 100.9 g∙N∙ha−1∙d−1) and unfertilized (24.5 - 56.5 g∙N∙ha−1∙d−1) plots, with the highest values predicted by DailyDayCent, including an unusual peak for the unfertilized field (110.1 g∙N∙ha−1∙d−1). The ECOSSE-simulated values correlated well with the measured values (R2 = 0.33, p < 0.05) under fertilized conditions only (Table 2). The total bias and error differences between simulated and measured values did not vary significantly with their 95% confidence levels. Similar results were observed for the unfertilized field but the coefficient of determination was poor (R2 = −0.02 to −0.04, negative).
For the fertilized fields, both DNDC (87%) and DailyDayCent (81%) underestimated the total N2O fluxes
Table 2. Statistical comparison between simulated and the measured daily soil NO3− concentration (kg∙N∙ha−1) and N2O fluxes (g∙N∙ha−1) from a conventionally-tilled plot cropped to spring barley.
*Significant at 5% level of probability. R2 = Coefficient of Determination; RMSE = Root Mean Square Error; RE = Relative Error (Mean); MD = Mean Difference; n = 130.
Figure 2. Comparison of field measured (seasonal/annual; open circle/square with vertical bars as standard errors) N2O emissions (g∙N∙ha−1) with values simulated (line) using the three process-based models over 8 years, commencing from 17 August 2003 (day of harvest), in conventionally-tilled fertilized arable land cropped to spring barley. Solid arrows indicate the days of ploughing and dotted arrows indicate the days of N fertilizer application.
(seasonal/annual), whilst the ECOSSE model overestimated these by 59% (Table 3). Based on an 8-year average, the DNDC simulated total N2O fluxes for the fertilized (207 kg∙N∙ha−1) and unfertilized (81 kg∙N∙ha−1) fields were 2 - 15 times lower than the estimates provided by the other two models. The DNDC-simulated values resulted in significant underestimation of N2O EFs, whilst those derived from DailyDayCent, and the ECOSSE were closer to those calculated from measured values. Compared to the annual EFs derived from measured values, DNDC was 94% lower, DailyDayCent 44% lower, and ECOSSE 35% lower. An estimation discrepancy for total fluxes between integrated values and the corresponding sum of daily fluxes and thereby EFs was observed. On an 8-year average, the simulated N2O EF was 0.09% with DNDC, 0.31% with DailyDayCent and 0.52% with the ECOSSE model.
3.2.2. Ecosystem and Heterotrophic Respiration
The Reco measured using EC from the large fertilized plot reached a maximum flux of 75.6 kg∙C∙ha−1∙d−1 during crop growth that decreased to 0.59 kg∙C∙ha−1∙d−1 during the non-crop period, corresponding to RH (Figure 3). The DNDC simulated values for Reco showed trends similar to the measured values, with an R2 of 0.34 (p < 0.05), and the total bias and error differences between simulated and measured values were ≤34% and ≤91%, respectively (Table 4). The estimated Reco for the DailyDayCent model also showed trends similar to the measured values, with higher fluxes from 2007 onwards, with an R2 of 0.41 (p < 0.05) and relatively small total bias (≤50%) and error (≤85%) difference between simulated and measured values. All models simulated RH satisfactorily, with R2 ranging from 0.44 - 0.62 (p < 0.05), with a small bias (≤50%) and error (≤87%) difference between simulated and measured values.
The annual total Reco measured using the EC was on average 6771 kg∙C∙ha−1, which is closer to the DailyDayCent value (6736) but higher than the DNDC estimate (4455; Table 4). Based on a 4-year average, the estimated RH based on measured values of Reco was 3624 kg∙C∙ha−1, which is closer to the ECOSSE and the DailyDayCent simulated values, but higher than the DNDC estimate (1794 kg∙C∙ha−1). Based on an 8-year average, the RH differed somewhat from the 4-year average though the simulated amount was similar to the measured value.
3.2.3. CH4 Fluxes
The measured CH4 fluxes (emission and oxidation) were small and differed significantly between the fertilized (−040 to 0.36 g∙C∙ha−1∙d−1) and unfertilized (−0.09 to 0.12) plots (Figure 4). The highest simulated oxidation and emission, respectively, were 2.92 and 0 g∙C∙ha−1∙d−1 with DNDC, 4.02 and 0 with DailyDayCent, and 0.24 and 0.31 with ECOSSE, providing values closer to the measured ones for the fertilized plot only. The DNDC and the DailyDayCent simulated values correlated poorly with the measured values (R2 = 0.02), demonstrating
Table 3. Statistical comparison between simulated and the measured seasonal and annual N2O fluxes (g∙N∙ha−1) and emission factors (EFs) derived from a conventionally-tilled field cropped to spring barley.
*Integrated (harvest to harvest); **Sum of daily simulated values (harvest to harvest); EF = Emission factor.
Figure 3. Comparison of Eddy Covariance measured (3 years; open circle) CO2 emissions (Reco, soil respiration only; kg∙C∙ha−1) with values simulated (and/or estimated using DNDC-derived fractions for DailyDayCent and ECOSSE; lines, solid/dotted including RH, heterotrophic respiration) using the three process-based models over 8 years, commencing from 17 August 2003 (day of harvest), in conventionally-tilled fertilized arable land cropped to spring barley. Solid arrows indicate the days of ploughing and dotted arrows indicate the days of N fertilizer application.
Table 4. Validation of daily soil (Reco) and heterotrophic respiration (RH) simulated by three process-based models with values measured from spring barley fields.
*Significant at 5% level of probability. ! = estimated using DailyDayCent derived ratio; f = Reco estimated using a conversion ratio derived from DNDC outputs for ECOSSE and DailyDayCent. R2 = Coefficient of determination; RMSE = Root Mean Square Error; RE = Relative Error (Mean); MD = Mean Difference.
Figure 4. Comparison of field measured (seasonal/annual; open circle/square with vertical bars as standard errors) CH4 emissions/oxidation (g∙N∙ha−1) with values simulated (line) using the three process-based models over 8 years, commencing from 17 August 2003 (day of harvest), in conventionally-tilled fertilized arable land cropped to spring barley. Solid arrows indicate the days of ploughing and dotted arrows indicate the days of N fertilizer application.
large bias and error differences between them (Table 5). The ECOSSE simulated values correlated well (R2 = 0.34, p < 0.05), with the total bias and error differences between simulated and measured values were at < 95% confidence intervals. For the unfertilized plots, either model estimates showed low R2 values, large total biases and errors.
Annual budgets based on the measured data showed the arable land to be a small CH4 source, with an emission of 2.35 g∙C∙ha−1 from the unfertilized plot, increasing to 3.50 g∙C∙ha−1 for the fertilized plot (Table 5). For an 8-year average, the model estimates indicated that cropland was a sink for CH4, with the annual oxidation of 666 g∙C∙ha−1 predicted by DNDC, 704 g∙C∙ha−1 from the DailyDayCent and 28 g∙C∙ha−1 from the ECOSSE models. The integrated and sum of the daily flux approaches taken to calculate total CH4 fluxes showed a small difference.
4.1. Performance of the Models in Simulating Nitrate-N Concentrations
Compared to the unfertilized field, the measured soil concentration in the fertilized plot following CAN fertilizer application was higher. The decrease in soil concentration over time is pre-assumed to have been due to plant uptake and other N loss processes, such as leaching. The DailyDayCent and DNDC-predicted
Table 5. Validation of daily CH4 effluxes (g∙C∙ha−1∙d−1) simulated by three process-based models with values measured from spring barley fields and their total fluxes.
*Significant at 5% level of probability. R2 = Coefficient of determination; RMSE = Root Mean Square Error; RE = Relative Error (Mean); MD = Mean Difference.
soil levels are noisy, and attributed to a mismatch between plant N uptake and other N loss processes, or due to errors associated with the crop growth module, particularly for the latter model. Moreover, the simulated upper or lower limit of levels are generally above the measurement values. Similar large overestimations were reported by others     using the earlier DNDC versions. This may be ascribed to limitations in the ability of the models to accurately account for variable soil depth increments or movement through soils and/or high mineralization of N and rapid nitrification.
The ECOSSE simulated values are closer to the amount of NO3-N applied, in line with the results of other researchers  , and consistent over the 8 years, unlike the simulations predicted by the other two. With the exception of fertilizer-induced peaks in, there were small differences between fertilized and unfertilized plots. Statistical evaluations confirm that the ECOSSE model can simulate soil well (R2 = 0.50) for the fertilized field, and that it performs better than DNDC (R2 = 0.31) and DailyDayCent (R2 = 0.14). However, all the models have difficulties in predicting background levels of soil, though the total bias and error differences are within their 95% confidence intervals, in line with the DNDC v9.2 estimates  .
4.2. Simulation Capacity of the Models for GHG Emissions
4.2.1. N2O Emissions
Simulated N2O emissions using the three models are reasonably consistent over the different years. However, all models are unable to predict N2O fluxes less than zero. This contrasts with the measured values where a sink of N2O under conditions of low oxygen and/or mineral N was observed   . DailyDayCent simulated N2O fluxes reasonably well. The total bias and error differences are somewhat large but within their 95% confidence levels. This indicates a high predictive potentials for the models although only the ECOSSE-simulated values show a significant relationship with the measured ones under fertilized conditions, in agreement with other observations   . Higher R2 values for cumulative N2O fluxes derived from DNDC simulated values were reported by other workers  but these did not correspond with the daily fluxes. Similar strong relationships for total N2O fluxes might be achieved from the other two models; however, the overall performance mainly depends on daily fluxes. Given that the simulated values were within the 95% confidence intervals, none of the models simulated daily N2O fluxes particularly well for the unfertilized field. Thus, an appropriate methodological compromise for the calculation of EFs may not be achievable without considering other factors controlling N2O emissions.
For the fertilized fields, both DNDC and DailyDayCent underestimated, and the ECOSSE model overestimated the total N2O fluxes (seasonal/annual). Based on an 8-year average, DNDC simulated total fluxes are 2 - 15 times lower than the DailyDayCent and the ECOSSE estimates. The variations among the model estimates and their relationship with key driving forces such as soil water and levels are assumed to be functionally related to the production and release of N2O  . These are consistent with the DNDC and the DailyDayCent-simulated peaks. This indicates that both models consider denitrification as the major contributor to N2O production. In contrast, nitrification might be the major pathway in the ECOSSE model and the N2O fluxes are reasonably consistent with the simulated levels; however, further enhancement of the emissions is possible with increasing soil water contents provided that substrate supply is not limiting  . The results from ECOSSE are in general agreement with the literature values for total N2O emissions from crop fields, which range from 0.7 to 3.5 kg∙N∙ha−1∙yr−1  -  . The inconsistencies and uncertainties in the modelled predictions are thought to be partially associated with differences in model version, methods of data analyses, and crop managements between years as well as the use of default values particularly for DNDC and DailyDayCent.
Similarly, a large underestimations of N2O EFs by DNDC as well as by DailyDayCent, compared to the measured data, are evident. Estimation of EFs using simulated values is constrained by total flux differences between the fertilized and unfertilized plots. Replacement of unfertilized values by using background annual N2O emissions of 1 kg∙N∙ha−1  could also be erroneous. Similar overall underestimations, particularly using the earlier versions of DNDC, have been reported    . This may be attributed to the limited number of field measurements, as this could result in large uncertainties in the measured values  . The fact that the DailyDayCent model also underestimated the N2O EF (44%), with similar findings using DayCent (~25%), when compared with the default annual value, was reported by Del Grosso et al. (2005). In contrast, the ECOSSE model, on average, increased the EF by 35%, but was within a closer range of the measured estimate (0.52%). This is in line with the previous version of the model used  , although it is still lower than the IPCC default value (1%).
There is a discrepancy between the integration approach and the corresponding sum of daily fluxes in calculating the total/cumulative N2O fluxes, which may under or overestimate the values, depending on the corresponding peak sizes, and thereby influence the EFs. Nitrous oxide emissions show large temporal and/or spatial variability  , resulting in an EF uncertainty of >50%   . This uncertainty could be more significant over several years of measurements than management-induced variations   . However, the lower measured total N2O fluxes and the corresponding EFs may be explained by the application of CAN during a relatively dry period, leading to less denitrification, and a higher SOC density, which is favourable for complete denitrification to occur under anoxic conditions. The above statement remains equivocal due to intermittent gas sampling, suggesting the need for intensive/continuous sampling to cover the impact of tillage, fertilization, rainfall events and other environmental factors that regulate the degree of N2O emissions. Moreover, further improvements of the models by identifying errors associated with the processes that interactively produce N2O are imperative. The use of more robust measurement protocols are also required for accurate validation and calculation of N2O EFs across disaggregated land use types and in response to different management practices.
4.2.2. Ecosystem and Heterotrophic Respiration
The simulated values for Reco from both DailyDayCent and DNDC demonstrated good correlation with the measured values (R2 = 0.41 versus 0.34, p < 0.05), with relatively small total bias and error differences. Likewise, DNDC simulates well the cumulative CO2 fluxes of cropland sites in Europe, except for some overestimation of net CO2 uptake  . It was found that DayCent provided robust simulated CO2 emissions for various land use systems  . Our study indicates that a further improvement of both models is required to remove the discrepancy with regard to the mis-match between the simulated CO2 peaks and the measured values. Irrespective of the models, this shift cannot be seen for RH in the modelled values, which correlated well (R2 = 0.44 - 0.62; p < 0.05), with small biases and errors. The DailyDayCent and DNDC simulated RH better than ECOSSE. The relatively poor performance of ECOSSE is probably due to the estimation errors associated with the requirement to convert the measured Reco data to RH.
Accordingly, the annual estimates for the total Reco measured using EC (6771 kg∙C∙ha−1) is closer to the global cropland average (5440 ± 800 kg∙C∙ha−1)  . The measured value is similar to the DailyDayCent predictions whilst DNDC underestimated it by 34%, which may be attributed to the poor prediction of soil water content between 50% and 80% WFPS   or around 55% - 60%  . Soil water as a confounding factor together with soil temperature also regulates the processes involved in decomposition  . Decomposition decreases as the soil dries the extent of the decrease is determined by diffusion limitations and the availability of oxygen  . Higher temperatures are often accompanied by low water contents and vice versa  and the strong interdependencies between these two factors make it difficult to separate their effects on soil respiration. Moreover, it was reported that DNDC may not predict Reco perfectly, due to some limitations in the crop growth module  . The crop growth module has an application of a sigmoid curve based upon degree days which require additional parameters, e.g. base temperature, degree days of phenology stages and radiation use efficiency to correctly define the growth curves for crops in terms of temporal carbon take up. Despite having a lower R2, the ECOSSE model simulated total RH values that were closer to the measured values, whilst the DailyDayCent and DNDC values underestimated it. Similar values based on the 8-year average were also observed, implying that the ECOSSE predicts RH better than the other two models. This can be attributed to the absence of crop growth and biomass-related inputs to run the models.
4.2.3. CH4 Emission/Oxidation
The measured data demonstrated both CH4 emissions and oxidation though the magnitude of the fluxes was relatively small. This might be linked to the contribution of RH with simultaneous influence of mainly soil water contents/precipitation events creating aerobic and anaerobic conditions  . The measured peaks for CH4 show a stimulating effect of N fertilization on both emissions and oxidation. The three models are unable to predict similar trends, and are assumed to be constrained by several regulating factors, so that consideration of the flux variations between fertilized and unfertilized fields may not be appropriate for judging the model’s performance. Both DNDC and DailyDayCent predict CH4 oxidation, in line with former versions   , and the ECOSSE also provided values closer to the measured ones, showing small sources, particularly for the fertilized field. The simulated and measured values are poorly correlated, including remarkably high total bias and error differences, particularly with the DNDC and DailyDayCent models. The modified version of ECOSSE used in this study performs better (R2 = 0.34) than the previous one  . Although the CH4 fluxes from arable lands might have less impact in terms of overall GHG accounting the functional relationship between C and N emissions need to be improved for accurate estimations.
The measured data show that cropland is a CH4 source that increases with the application of N fertilizer, indicating fertilizer-induced limitation for CH4 oxidation to occur. Arable soils are mainly considered as a sink rather than a source of CH4   . Indeed, tillage and N fertilization have a tendency to reduce oxidation potentials  . Increased CH4 oxidation in arable soils  may be linked to well aerated conditions with a positive redox potential that limits methanogenic activities through draining, coupled with ploughing  . The three models also demonstrated an annual reduction in CH4 oxidation from the N fertilizer-treated plots compared to the unfertilized ones. However, the higher predicted oxidation with the DNDC and DailyDayCent models compared to ECOSSE indicates that annual estimates are not in agreement with the amounts of mineral N either mineralized and/or applied as fertilizer. These results suggest that functional constraints on the DailyDayCent and DNDC models for predicting CH4 emissions were greater, in comparison to the ECOSSE model and assumed to be mainly due to input parametric limitations.
Compared to the measured values, ECOSSE could simulate nitrate concentration more robustly than DNDC and DailyDayCent. Both DNDC and DailyDayCent underestimated daily and total N2O fluxes compared to ECOSSE, providing an improved prediction of fertilizer-induced N2O fluxes and EFs. All models could simulate soil and/or heterotrophic respiration adequately, except for an underestimation with DNDC that may be related to the greater impact of variations in soil properties compared to other model predictions. Only the ECOSSE model was able to predict field CH4 emissions/oxidation that were closer to the measured ones, and demonstrate the overall dominance of oxidation processes. There are constraints in terms of processes and driving forces in all the models for predicting coupled C and N emissions, leading to the underestimation of GHGs. Thus, refinement and further validation of the models using country-specific activity data are required to better predict GHG emissions. In addition, to avoid a dependency on default inputs that may lead to significant errors in the model outputs more measurements are required that account for temporal and spatial variability. Furthermore, validations and sensitivity tests need to focus more on site-related characteristics, land use differences, management interventions, and climatic factors for providing national GHG estimates and for identification of mitigation options.
The senior author gratefully acknowledges the funding by the Science, Technology, Research and Innovation for the Environment (STRIVE) Programme of the Irish Government under the National Development Plan 2007-2013 and the Department of the Environment, Heritage and Local Government. The authors would like to thanks Phillip O’Brien (EPA) for extending technical and relevant support; Mike Williams, Mike Jones and Matt Saunders (TCD), Komsan Rueangritsarakul and Mohamed Helmy (UCD) for supplying experimental data for modelling work; as well as Tom Bolger and Tommy Gallagher (UCD) for providing administrative support.