Agronomy Publications Agronomy
4-2015 Understanding the DayCent model: Calibration, sensitivity, and identifiability through inverse modeling Magdalena Necpálová University of Wisconsin–Madison
Robert P. Anex University of Wisconsin–Madison
Michael N. Fienen United States Geological Survey
Stephen J. Del Grosso Colorado State University
Michael J. Castellano IFoowlalo Swta tthie Usn iaverndsit ay,dd casitteionlmj@ial wasorktates.e duat: http://lib.dr.iastate.edu/agron_pubs Part of the Agricultural Science Commons, Agriculture Commons, and the Agronomy and Crop See next page for additional authors Sciences Commons The ompc lete bibliographic information for this item can be found at http://lib.dr.iastate.edu/ agron_pubs/99. For information on how to cite this item, please visit http://lib.dr.iastate.edu/ howtocite.html.
This Article is brought to you for free and open access by the Agronomy at Iowa State University Digital Repository. It has been accepted for inclusion in Agronomy Publications by an authorized administrator of Iowa State University Digital Repository. For more information, please contact [email protected]. Understanding the DayCent model: Calibration, sensitivity, and identifiability through inverse modeling
Abstract The ba ility of biogeochemical ecosystem models to represent agro-ecosystems depends on their correct integration with field observations. We report simultaneous calibration of 67 DayCent model parameters using multiple observation types through inverse modeling using the PEST parameter estimation software. Parameter estimation reduced the total sum of weighted squared residuals by 56% and improved model fit ot − crop productivity, soil carbon, volumetric soil water content, soil temperature, N2O, and soilNO3 compared to the default simulation. Inverse modeling substantially reduced predictive model error relative to the default − + model for all model predictions, except for soil NO3 and NH4 . Post-processing analyses provided insights into parameter–observation relationships based on parameter correlations, sensitivity and identifiability. Inverse modeling tools are shown to be a powerful way to systematize and accelerate the process of biogeochemical model interrogation, improving our understanding of model function and the underlying ecosystem biogeochemical processes that they represent.
Keywords DayCent model, Inverse modeling, PEST, Sensitivity analysis, Parameter identifiability, Parameter correlations
Disciplines Agricultural Science | Agriculture | Agronomy and Crop Sciences
Rights © 2014 The Authors. Published by Elsevier Ltd. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).
Authors Magdalena Necpálová, Robert P. Anex, Michael N. Fienen, Stephen J. Del Grosso, Michael J. Castellano, John E. Sawyer, Javed Iqbal, Jose L. Pantoja, and Daniel W. Barker
This article is available at Iowa State University Digital Repository: http://lib.dr.iastate.edu/agron_pubs/99 Environmental Modelling & Software 66 (2015) 110e130
Contents lists available at ScienceDirect
Environmental Modelling & Software
journal homepage: www.elsevier.com/locate/envsoft
Understanding the DayCent model: Calibration, sensitivity, and identifiability through inverse modeling
* Magdalena Necpalov a a, , Robert P. Anex a, Michael N. Fienen b, Stephen J. Del Grosso c, Michael J. Castellano d, John E. Sawyer d, Javed Iqbal d, Jose L. Pantoja d, Daniel W. Barker d a Department of Biological System Engineering, University of WisconsineMadison, Madison, WI, USA b United States Geological Survey, Wisconsin Water Science Center, Middleton, WI, USA c Natural Resource Ecology Laboratory, Colorado State University, Fort Collins, CO, USA d Dept. of Agronomy, Iowa State University, Ames, IA, USA article info abstract
Article history: The ability of biogeochemical ecosystem models to represent agro-ecosystems depends on their correct Received 11 August 2014 integration with field observations. We report simultaneous calibration of 67 DayCent model parameters Received in revised form using multiple observation types through inverse modeling using the PEST parameter estimation soft- 12 December 2014 ware. Parameter estimation reduced the total sum of weighted squared residuals by 56% and improved Accepted 12 December 2014 model fit to crop productivity, soil carbon, volumetric soil water content, soil temperature, N2O, and soil Available online 21 January 2015 NO3 compared to the default simulation. Inverse modeling substantially reduced predictive model error þ relative to the default model for all model predictions, except for soil NO3 and NH4 . Post-processing Keywords: e DayCent model analyses provided insights into parameter observation relationships based on parameter correlations, fi Inverse modeling sensitivity and identi ability. Inverse modeling tools are shown to be a powerful way to systematize and PEST accelerate the process of biogeochemical model interrogation, improving our understanding of model Sensitivity analysis function and the underlying ecosystem biogeochemical processes that they represent. Parameter identifiability © 2014 The Authors. Published by Elsevier Ltd. This is an open access article under the CC BY-NC-ND Parameter correlations license (http://creativecommons.org/licenses/by-nc-nd/4.0/).
Software and data availability section
Software DayCent model Developers W. J. Parton, S. J. Del Grosso, S. Ogle, K. Paustian, Contact address Natural Resource Ecology Laboratory, Colorado State University, Fort Collins, CO, USA Telephone (970) 491-2195 Email address [email protected] Abbreviations: ANPP, aboveground net primary productivity; ARS, Agricultural Year first available 1998 Research Service; C, carbon; CEC, cation-exchange capacity; CH4, methane; C/N Hardware required PC with at least 512 K of RAM. A graphics ratio, carbon to nitrogen ratio; d, index of agreement; DEFAC, decomposition factor; DNDC, denitrification decomposition model; EPA, Environmental Protection adapter (CGA, EGA, VGA, or Hercules þ Agency; GHG, greenhouse gas; GML, GausseMarquardteLevenberg; NH4 , monographic) and 2 Mb of disk space are ammonium cation; J, Jacobian matrix; N, nitrogen; N2O, nitrous oxide; NPP, net recommended. primary productivity; NO3 , nitrate anion; PEST, parameter estimation software; Software required Windows MB, mean bias; RMSE, root mean square error; rRMSE, relative root mean square Availability and cost Available on request; Free error; SOC, soil organic carbon; SOM, soil organic matter; SVD, singular value decomposition; SWSR, sum of weighted squared residuals; VSWC, volumetric soil Software PEST version 13.0 water content; UAN, urea ammonium nitrate. Developer John Doherty * Corresponding author. Sustainable Agroecosystems, ETH Zürich, Tannenstrasse Contact address Watermark Numerical Computing, 336 þ 1, 8092 Zürich, Switzerland. Tel./fax: 41 (0)446338607. Cliveden Avenue, Corinda 4075, Australia E-mail addresses: [email protected] (M. Necpalov a), [email protected] (R.P. Anex), mnfi[email protected] (M.N. Fienen), [email protected] Telephone 07 3379 1664 (S.J. Del Grosso), [email protected] (M.J. Castellano), [email protected] Email address [email protected] (J.E. Sawyer), [email protected] (J. Iqbal). http://dx.doi.org/10.1016/j.envsoft.2014.12.011 1364-8152/© 2014 The Authors. Published by Elsevier Ltd. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/). M. Necpalov a et al. / Environmental Modelling & Software 66 (2015) 110e130 111
Year first available 1994 emissions from agricultural soils for the US National Greenhouse Hardware required Desktop or Laptop Gas Inventory compiled by the EPA (Olander and Haugen-Kozyra, Software required Windows or Linux 2011) and reported annually to the UN Framework Convention on Availability and cost down load from: http://www. Climate Change (US EPA, 2014). DayCent consists of sub-models for pesthomepage.org/Downloads.php; soil water content and temperature by layer, plant production and Free. allocation of net primary production (NPP), decomposition of litter and soil organic matter (SOM), mineralization of nutrients, N gas emissions from nitrification and denitrification, and CH4 oxidation in unsaturated soils. The accuracy with which a model represents the natural system 1. Introduction observed in the field depends on how completely the underlying biophysical processes are represented in the model and how well the Greenhouse gases (GHG) released from the soils of terrestrial model parameters are calibrated to field observations. Like other ecosystems are highly variable in space and time due to the inter- biogeochemical process-based models, DayCent is typically cali- action of climatic drivers and ecosystem processes involved in brated manually by adjusting one parameter at a time, thus the carbon (C) and nitrogen (N) transformation associated with pro- calibrated parameters are adjusted in an iterative fashion in multiple duction and consumption of GHGs (Müllera et al., 2002; Rahn et al., stages (Wallach et al., 2014). At each stage, specific processes are 2012; Wrage et al., 2001). Field measurements that capture the high targeted (e.g., plant growth and yield, SOC), and the most influential temporal and spatial variability of N2O fluxes (Bouwman et al., parameters are adjusted to match simulated to observed values (Del 2002; Parkin, 2008; Snyder et al., 2009) or the high spatial vari- Grosso et al., 2011). This approach, however, does not guarantee full ability of soil organic carbon (SOC; Conant and Paustian, 2002; extraction of information from the field observations and it is diffi- Kravchenko and Robertson, 2011) are expensive, time intensive, cult to know when calibration correctly balances the performance of and unable to capture the full range of ecological and environ- all model components (Nolan et al., 2011). It is generally accepted mental conditions. When properly informed by field observations, that manual calibration of complex ecosystem models does not ecosystem process-based models are a powerful way to investigate necessarily yield optimal parameter estimates, is somewhat arbi- the effects of management practices on GHG emissions or SOC from trary, and results in high uncertainty in model parameters and different ecosystems, soils, and climates. simulated variables (Schwarz et al., 2006). Inverse modeling, based A number of biogeochemical models have been developed and on an objective statistical method and mathematical techniques for used to quantify GHG emissions and SOC at both plot and landscape stable parameter estimation, has become a widely accepted way to scales, e.g., Century (Parton et al., 1994; Parton, 1996), DayCent (Del enhance the transfer of information contained in field observations Grosso et al., 2005; Parton et al., 1998), deni- to model parameters (Doherty, 2003; Doherty and Hunt, 2010a; trificationedecomposition (DNDC) (Li et al., 2000), ecosys (Grant Hunt et al., 2007). Despite mathematical objectivity, some subjec- et al., 1993) and EPIC (Wang, 2005). These models are mathemat- tivity is unavoidable: through defining the conceptualization of the ical representations of our understanding of the complicated, inverse problem and making a set of decisions related to regulari- coupled biogeochemical soil processes that allow us to test our zation, parameter bounds, observation weighting strategy, etc. understanding through comparison of model results with obser- (Fienen, 2013). The inverse modeling tool PEST (Doherty, 2010) uses vations, and predict responses to conditions that have not yet been an iterative, nonlinear regression approach that involves simulta- observed, such as ecosystem responses to changing climate. Thus neous adjustment of multiple model parameters and evaluation of these models have become important tools in the study of model fit by the sum of weighted squared residuals between field biogeochemical cycles. Model development is based on a quanti- observations and simulated values. In addition to providing so- tative understanding of the interactions among physical, chemical phisticated estimates of the parameter values that provide the best and biological processes that is critical for predicting the ecosystem possible fit for a given calibration problem, inverse modeling pro- response to land use or climate change. The individual underlying vides a method for comprehensive model analysis through statistical processes are represented by sets of equations in component measures such as the variance/covariance matrix, parameter corre- models that are coupled together to describe a full system (Wallach lations, confidence intervals, sensitivities, identifiability, and pre- et al., 2014). Models usually have a mechanistic structure that re- dictive uncertainty analysis (Moore and Doherty, 2005, 2006). flects our understanding of the processes governing the system These tools can help users recognize model problems that are behavior. Many ecosystem models utilize several hundred param- difficult to identify with manual calibration methods (Hill and eters representing individual physical quantities or combinations of Tiedeman, 2007; Poeter and Hill, 1997). For example, it has been physical quantities that may not be observable through direct repeatedly observed that only a small number of the many pa- measurement. It is thus impossible to measure the sensitivity of rameters used in most environmental models are uniquely esti- system behavior to each of these parameters and information on mable with most datasets (Beck and Halfon, 1991; Beven and Freer, their identifiability through field observations is often not available. 2001; Doherty and Hunt, 2009). The inability to uniquely identify Yet for model users and particularly model developers, an under- certain model parameters can be the result of their high correlation standing of how model parameters influence the simulation of with other parameters, or lack of sensitivity of the model outputs to target ecosystem processes and which field observations are most these parameters. This sort of problem is extremely difficult to useful in defining parameter values is essential. recognize without specialized tools and can lead to misidentifica- The DayCent model is a widely used terrestrial biogeochemical tion of parameter values, model over-fitting, and inaccurate model process-based model of intermediate complexity (Del Grosso et al., projections for conditions outside the range of the calibration 2001, 2002; Parton et al., 1998). It has been used to simulate dataset. Applying inverse modeling tools provides valuable insight ecosystem responses to changes in climate and agricultural man- about parameter dependencies, which parameters are exerting the agement practices in crop, grassland, forest and savanna ecosys- most influence on the simulated values, whether the field obser- tems (Brilli et al., 2013; Cheng et al., 2014; Del Grosso et al., 2008a, vations contain enough information to estimate the model pa- 2009; Hartman et al., 2009; Parton et al., 2007; Parton and rameters, and the uncertainty associated with the predictions Rasmussen, 1994). In the USA, it has been used to quantify N2O based on the estimated parameter values. 112 M. Necpalov a et al. / Environmental Modelling & Software 66 (2015) 110e130
Inverse modeling has been applied to soil physics and ground- rainfall in 2011, 2012, and 2013, respectively. The mean air temperature was 10.0 C water studies where inverse simulations are commonly used to in 2011, 11.6 C in 2012, and 8.5 Cin2013(Fig. 1). parameterize hydraulic functions (e.g., Bitterlich et al., 2005; Kosugi et al., 2001; Lin and Anderson, 2003; Spohrer et al., 2006; Tonkin 2.2. DayCent modeling and Doherty, 2009). It has been also applied to three biogeo- DayCent is a terrestrial ecosystem model designed to simulate fluxes of C and N chemical models, i.e. RZWQM (Malone et al., 2010, 2014; Nolan among the atmosphere, vegetation, and soil (Del Grosso et al., 2001; Parton et al., 1998). DayCent is the daily time-step version of the Century model (Parton et al., et al., 2010) to calibrate soil hydraulic, N leaching, N trans- 1994). The plant growth submodel simulates plant productivity as a function of formation and crop yield parameters; Forest DNDC (Lamers et al., genetic potential, phenology, nutrient availability, water/temperature stress, and 2007) and DayCent (Rafique et al., 2013) to calibrate parameters solar radiation (i.e. energy biomass conversion factor). NPP is allocated to plant associated with C and N trace gas production. These studies per- components (e.g., roots vs. shoots) based on vegetation type, phenology, and water/ fi formed calibration with a reduced number of degrees of freedom to nutrient stress. Nutrient concentrations of plant components vary within speci ed fi limits depending on vegetation type and nutrient availability relative to plant de- increase ef ciency. To our knowledge, no previous study has mand (Del Grosso et al., 2008a). SOM is simulated in the top 20 cm soil layer as a sum demonstrated the full use of inverse modeling by calibrating mul- of dead plant matter and three SOM pools (active, slow, and passive) on the basis of tiple components of an ecosystem biogeochemical model and decomposition rates. The amount of biomass decomposition products entering the explicating model function through the use of associated sensitivity pools and flowing between the pools depend on lignin content and C/N ratio, size of the pools, temperature/water factors, and clay content (Del Grosso et al., 2001, 2011; and identifiability methods. Parton et al., 1994). The decomposition of litter, SOM, and nutrient mineralization The objectives of this study were to demonstrate calibrating the are functions of substrate availability, lignin content, C/N ratio, water/temperature major components of the DayCent model with several types of field stress, and tillage intensity (Del Grosso et al., 2008a). The trace gas model contains both denitrification and nitrification submodels. Daily denitrification rates are observations simultaneously through inverse modeling as imple- mented in the PEST parameter estimation software (Doherty, 2010), calculated for each soil layer based on soil NO3 concentration distributed throughout the soil profile, heterotrophic respiration (i.e., available labile C), soil and to provide insights into the function of the model. Outputs water content, texture, and temperature; while nitrification rates are calculated þ simulated using the estimated parameter values were validated based on soil NH4 concentration, water content, texture, and the temperature in against an independent dataset. Using parameter correlations, the top 15 cm layer (Del Grosso et al., 2001, 2008a; Parton et al., 2001). The soil water sensitivity and identifiability we provide insights into the complex sub-model simulates soil water content and water fluxes (i.e., through the canopy, surface run off, leaching, evaporation and transpiration) for each horizon DayCent model structure and explore the relationships between throughout the defined depth of the soil profile (Parton et al., 1998). When the model parameters and observations. average daily air temperature is freezing, the precipitation is accumulated in the snowpack. Saturated water flow occurs on days that receive rainfall, irrigation, or snow melt every day on a sub-daily time step (Del Grosso et al., 2008a). Each soil 2. Material and methods layer is filled before water flows to the next layer and if the input rate is greater than saturated hydraulic conductivity, the difference goes into surface runoff (Del Grosso 2.1. Study location and data et al., 2001). Unsaturated flow simulated by tipping bucket method is calculated The study was conducted at Iowa State University Agricultural Engineering and using Darcy's law as a function of the hydraulic conductivity of soil layers and the Agronomy Ames Research Farm (42.01N, 93.78W). Long term mean annual tem- difference in hydraulic potential calculated for the centers of adjacent soil layers. The perature and rainfall are 9.4 C and 827 mm yr 1. The soil is predominantly Clarion flux from the top layer depends on potential evapotranspiration rate, water content loam series (fine loamy, mixes, superactive, mesic Typic Hapludolls) with pH 6.62, of the top soil layer, and the minimum water content set for the top layer (Del Grosso 1 cation-exchange capacity (CEC) 20.38, SOC content of 23.3 g kg , total N content et al., 2001). Soil NO3 is distributed throughout the soil profile and available for 1.77 g kg 1 soil at 0e10 cm. The soil texture is loam and clay loam (39.5% sand and leaching into the subsoil. Its movement and leaching are largely controlled by soil þ 38.2% silt at 0e10 cm). water flow and plant N uptake. The nitrification sub-model modifies soil NH4 SOC, crop aboveground net primary productivity (ANPP), grain yields, soil NO3 , concentration which is assumed to be immobile and distributed entirely in the top þ NH4 ,N2O emissions, soil temperature and volumetric soil water content (VSWC) 15 cm layer of soil (Del Grosso et al., 2008a). were measured over a 3 year period (2011e2013) as part of a field experiment DayCent version 4.5 was used to simulate SOC, crop ANPP, N2O fluxes and soil þ studying the effect of winter rye cover crop on soil N2O emissions from a NO3 and NH4 at the plot scale representing the area of the plot (6.1 15.2 m). The cornesoybean cropping system treated with different N fertilizer rates. A detailed initial SOC pools were generated and stabilized using a “spin-up” simulation of description of the sampling strategy and analytical procedures has been reported by native prairie ecosystem (mix of perennial C3 warm and cold grasses species and Mitchell et al. (2013) and Kladivko et al. (2014). This study used one treatment with symbiotic N2 fixing plants), and naturally occurring disturbances (1000e1800) fol- N rate of 135 kg N fertilizer ha 1, no cover crop, managed without tillage. During the lowed by simulation of historical land cover/use data (i.e., grazing until 1880, low corn phase, N fertilizer was applied as urea ammonium nitrate (UAN, 32% N), side- fertilized cornewheat-fallow rotation followed by more intensive cornesoybean dressed in bands at 15 cm depth. All the data, except for SOC and crop ANPP, were rotations) and current management (2010e2020). The dates of the current man- collected from the fertilizer bands which were, on an area basis, calculated to have agement events such as planting/harvest, and fertilizer application were consistent received twice the N fertilizer of surrounding areas. Weather data were collected with the operations in the fields during 2011e2013. The simulation was driven using from a meteorological station. The site received 815.3, 692.4, and 852.1 mm of site-specific measured weather data (i.e., daily rainfall, high and low air
Fig. 1. Daily air temperature and rainfall during calibration (2011e2012) and validation period (2013). M. Necpalov a et al. / Environmental Modelling & Software 66 (2015) 110e130 113 temperatures) from 1950 to 2012, which were used repeatedly over the course of Parameters included in inverse modeling (Table 2) were selected based on our prior simulation. Simulation soil profile comprising of 13 soil layers (0e180 cm) was knowledge of the model as being influential in the DayCent submodels. The default characterized by plot specific soil texture, SOC, bulk density and pH measured at parameter values were used as initial values and upper and lower bounds were 0e10, 10e20, 20e40, and 40e60 cm depth. The field capacity, wilting point, and specified by one of the DayCent developers (personal communication, S. Del Grosso saturated hydraulic conductivity for each depth were calculated using the Soil Water and B. Parton, December 6, 2013, Table 2). All parameters were log-transformed to Characteristics Calculator software (SWCC), version 6.02.74 (USDA Agricultural strengthen the linear relationships between parameters and simulated values Research Service, Washington). “Default simulation” here refers to simulation using (Doherty and Hunt, 2010a). Numerical stability of inverse modeling was ensured site specific soil and weather data with all model parameters set to their default through regularization using truncated singular value decomposition (SVD) of the values as defined by the model developers. The model was calibrated using field weighted Jacobian matrix on an iteration-by-iteration basis. The level of truncation observations, described above, for years 2011 and 2012, and validated against data was automatically calculated based on a stability criterion. This regularization collected in 2013 (Table 1). method transforms the original model parameters into linear combinations (i.e., eigenvectors), determines which are most sensitive (Moore and Doherty, 2005; Tonkin and Doherty, 2005), and truncates the transformed normal equations ma- 2.3. Inverse modelling trix, reducing the number of estimated parameters to maintain numerical stability Inverse modelling of DayCent parameters was accomplished using a model- and maximum reasonableness (Aster et al., 2005). The resulting regularized inver- independent parameter estimation software package PEST (Doherty, 2010). PEST sion process will not include parameters that are unidentifiable with the available allows independent parameter estimation using a nonlinear regression method data. When correlated parameters are included in the inversion, the SVD-based fi grounded in the principles of least-squares minimization (Doherty and Hunt, regression nds the maximum likelihood combination of the parameters that is 2010a); i.e. the model parameters are estimated in an iterative fashion as the code consistent with the observations. systematically varies model inputs, runs the model, reads model output, and eval- Details of the process of coupling the PEST software with DayCent model have uates the model fit using an objective function, which represents weighted least been reported by Rafique et al. (2013). The runtime of the parameter estimation squared difference between observations and simulated values (Doherty, 2010), process was decreased through use of the BeoPEST version of PEST designed for expressed as parallel processing on a distributed grid of processors (Hunt et al., 2010; Schreuder, 2009). h i T Pre-calibration parameter correlations were obtained from the correlation co- fðbÞ¼½y y0ðbÞ Q y y0ðbÞ efficient matrix by using a standard GausseMarquardteLevenberg parameter- where Q is a diagonal matrix with the squared observation weights on the diagonal, estimation method. Relative composite sensitivity of each parameter with respect y is a vector of observations, y′(b) is a vector of model outputs from the DayCent to each observation group and to the entire calibration dataset at the beginning of the model, based on parameter vector b, and collocated with the observations in y, and T parameter estimation process was computed as the magnitude of the column of the indicates matrix transpose. The symbology is adopted from Nolan et al. (2011). Pa- Jacobian matrix corresponding to the ith parameter with each entry in that column rameters that minimize this equation are obtained by solving the normal equations multiplied by the squared weight associated with the corresponding observation, using the GausseMarquardteLevenberg (GML) gradient search algorithm. At the multiplied with the absolute value of the parameter value (Doherty, 2005): start of each iteration, the relationship between the model parameters and the 1=2 ¼ T * v model-generated outputs is linearized through formulation as a Taylor expansion si J QJ i ii based on the current best parameter set. The matrix of all first-order partial de- rivatives of the simulated values that correspond to observations in the calibration where J is the Jacobian matrix and Q is the a diagonal matrix whose elements are dataset to the adjustable model parameters (i.e., the “Jacobian”) is computed using comprised of the squared observation weights and jvij is the absolute value of the fi the nite difference method. The linearized problem is then solved for a better parameter value. Parameters with relative composite sensitivity >10% of the parameter set using the GML algorithm, and the new parameters are tested by maximum relative composite sensitivity obtained for the set of parameters were fi running the model and computing the objective function as de ned above. The considered to be highly sensitive, while parameters with relative composite sensi- parameter changes and objective function improvement are compared with those of tivity <1% were considered to be insensitive to the inversion problem (Hill and fi the previous iteration to determine if another iteration is justi ed. If it is, the entire Tiedeman, 2007). The sensitivity analysis was used as a diagnostic of the inverse process is repeated; if not, the parameter estimation process terminates (Doherty, process, not as a standalone analysis. Therefore, the initial values and the assump- 2010). tions were consistent with those used in the inversion. fi The calibration dataset comprises 111 eld observations (Table 1). A weighted Parameter identifiability, which represents the ability of a calibration dataset to multicomponent objective function was adopted, similar to that used in Nolan et al. constrain model parameters (Doherty and Hunt, 2009), was calculated through SVD (2011) and Lin and Radcliffe (2006). Each observation type was assigned to a of the weighted Jacobian matrix computed on the basis of initial parameter values. different observation group, and each observation group formed a component of the The boundary between the solution and null subspaces was set at a specific singular objective function which was representative of that type of observation (e.g., soil value calculated using the SUPCALC utility (Doherty, 2008) as described by Doherty fi temperature measurements). An inter-group weighting strategy was de ned using and Hunt (2009). The number of singular vectors used to compute identifiability PEST utility PWTADJ1 (Doherty and Welter, 2010) such that each group contributed differed between the observation groups from 5 to 11 by means of the different equally to the objective function at the start of the estimation process. This ensures number of field observations. The identifiability of the ith parameter was calculated that each observation group contributes information to the process, regardless of as a sum of the squared ith components of all eigenvectors spanning the calibration number of observations per group, units of measure, and other confounding factors. solution space (V1)(Doherty, 2010): Individual field observations were weighted equally within each group. Prior to the parameter estimation process, field capacity and wilting point in the ¼ Sð Þ2 Idi V1i DayCent input soil file were adjusted as described by Del Grosso et al. (2011).
Table 1 Field observations included in DayCent calibration using an inverse modeling approach. Sampling strategy and analytical procedures were described by Mitchell et al. (2013) and Kladivko et al. (2014).
Observation group N Units Data source Sampling frequency