The Scientific World Journal

Volume 2012 (2012), Article ID 590287, 8 pages

doi:10.1100/2012/590287

Research Article

Simulation of Soil Temperature Dynamics with Models Using Different Concepts

Renáta Sándor and Nándor Fodor

Centre for Agricultural Research, Hungarian Academy of Sciences, 2462 Martonvásár, Hungary

Received 20 March 2012; Accepted 14 April 2012

Academic Editors: A. Ferrante, R. M. Mian, T. Rennert, and J. Viers

Copyright © 2012 Renáta Sándor and Nándor Fodor. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

This paper presents two soil temperature models with empirical and mechanistic concepts. At the test site (calcaric arenosol), meteorological parameters as well as soil moisture content and temperature at 5 different depths were measured in an experiment with 8 parcels realizing the combinations of the fertilized, nonfertilized, irrigated, nonirrigated treatments in two replicates. Leaf area dynamics was also monitored. Soil temperature was calculated with the original and a modified version of CERES as well as with the HYDRUS-1D model. The simulated soil temperature values were compared to the observed ones. The vegetation reduced both the average soil temperature and its diurnal amplitude; therefore, considering the leaf area dynamics is important in modeling. The models underestimated the actual soil temperature and overestimated the temperature oscillation within the winter period. All models failed to account for the insulation effect of snow cover. The modified CERES provided explicitly more accurate soil temperature values than the original one. Though HYDRUS-1D provided more accurate soil temperature estimations, its superiority to CERES is not unequivocal as it requires more detailed inputs.

1. Introduction

Soil temperature () is one of the most important variables of the soil. It can significantly influence seed germination [1], plant growth [2], uptake of nutrients [3], soil respiration [4, 5], soil evaporation [6], and the intensity of physical [7], chemical [8, 9], and microbiological processes [10, 11] in the soil.

Solar radiation and air temperature are the main driving forces determining the soil temperature which is influenced by numerous other factors such as precipitation, soil texture, and moisture content as well as the type of surface cover (plant canopy, crop residue, snow, etc.) [12]. Yearly, monthly, or daily means of soil temperature measurements are frequently reported, but the variability of is similarly important [13]. In spite of this, at many meteorological stations only aboveground variables (e.g., air temperature) are observed, or the soil temperature sensors are installed at the station (close to the mast that supports other sensors and the data logger) and not at the plots of the experimental site which could make the measured data unrepresentative.

If soil temperature is not measured, several methods are available to calculate it using meteorological variables and other parameters. As the simple air-temperature-based methods (e.g., [14]) provided inadequate data, an improved formula was introduced that uses precipitation data [15] as well. There are three types of soil temperature models [16]: (1) empirical models that are based on statistical relationships between soil temperature at some depth and climatological and soil variables (e.g., [17]); (2) mechanistic models that focus on physical processes (radiative energy balance as well as sensible, latent, and ground-conductive heat fluxes) to predict the upper boundary temperature and estimate the temperature of deeper layers with Fourier’s equation (e.g., [18]); (3) mixed empirical and mechanistic models that calculate the temperature of different soil layers based on physical principles of heat flow, but the boundary temperature at the soil surface must be provided empirically (e.g., [19]).

Since LAI (leaf area index) and soil water balance strongly influence the, soil temperature dynamics, soil temperature calculating methods function more precisely when those are integrated into hydrological [20, 21] or crop simulation models [22, 23]. The primary purpose of these models is to describe the processes of the very complex atmosphere-soil-plant system, including human activities, using mathematical tools and to simulate them with the help of computers.

The objectives of this paper are as follows: (1) presenting the effect of LAI on soil temperature; (2) comparison of an empirical and a mechanistic soil temperature model using measured data; (3) enhancing the performance of the empirical model.

2. Materials and Methods

Data of the agrometeorological station at Őrbottyán, Hungary were used in the study. The arenosol of the experiment site has the following characteristics: bulk density: 1.67 gcm⁻¹; organic matter content: 0.91%; CaCO₃ content: 5.1%; sand fraction: 86.3%; silt fraction: 8.3%; clay fraction: 5.4% [24]. Saturated hydraulic conductivity (Ks) and characteristic points of the soil water retention curve (SWRC) were measured with Guelph permeameter [25] and Eijkelkamp sand/kaolin box apparatus, respectively. pF-measurements were carried out with 100 cm³ undisturbed samples taken in 5 replicates. The van Genuchten parameters [26] of the SWRC were determined with Soilarium [27]. The above parameters characterize the 0–20 cm layer of the soil. In the 20–60 cm layer, the parameter values are practically the same except for the organic matter content which gradually decreases to zero with depth.

Soil temperature sensors (thermistor type, ±0.5°C accuracy, 0.1°C resolution) were installed at 5 different depths (5, 10, 20, 40, and 60 cm) at the centre of each  m test plot of an experiment with 8 parcels realizing the combinations of the fertilized, nonfertilized, irrigated, and nonirrigated treatments in two replicates. Temperature data were recorded every 15 minutes. A meteorological station was installed next to the experiment where precipitation, relative humidity, wind velocity and direction, global radiation, and air temperature were measured every 5 minutes in 2010 and 2011. In these two years, maize was grown at the site. LAI of every parcel in three-week intervals were determined by direct measurements. Three plants were cut out randomly at every observation time, and the area of the leaves was calculated with Montgomery’s method [28].

Site-specific measured data were used as inputs for the CERES-Maize [29] crop simulation model as well as for the HYDRUS-1D [21] hydrological model. CERES is a daily-step deterministic model that simulates plant (assimilation, biomass accumulation, leaf area, dynamics and root growth) as well as soil (water, temperature and nutrient dynamics) processes using empirical equations. HYDRUS-1D is designed for simulating one-dimensional variably saturated water flow, heat movement, and the transport of solutes in the soil. It numerically solves the Richards equation for saturated-unsaturated water flow (including a sink term to account for water uptake by plant roots) and advection-dispersion type equations for heat and solute transport using Galerkin-type linear finite element schemes. Several studies proved the efficiency of both models [30, 31].

The soil temperature calculation module of CERES belongs to the empirical model group. When calculating the actual temperature () at a given depth (), this model takes into account that the upper soil layers absorb energy, and the heat needs time to reach the lower layers as in (1). The effect of the energy reaching the soil surface appears delayed and decreased in the lower soil layers. The extent of the delay and the decrease is a function of the actual average moisture content () and the average bulk density of the topsoil (BDavg). The model assumes a sinusoidal annual course of the soil surface temperature that is modified by an additive term of a five-day moving average of a factor described by (2) as follows:

where denotes the day of the year; equals 200 on the northern hemisphere, while it is 20 on the southern hemisphere. ALB is the albedo of the surface, and denote the average temperature and the average temperature difference of the site. , , and denote the daily mean and maximum temperature as well as the daily global radiation on the th day of the year, respectively. The term in (1) describes the delay of the effect of energy reaching the surface in deeper layers. The exponential term in (1) is more related to the heat capacity of the topsoil as it governs the decrease of the effect of the incoming energy at the surface in deeper layers.

The soil temperature calculation module of HYDRUS-1D belongs to the mechanistic model group. It numerically solves the convection-dispersion equation describing the one-dimensional heat transfer as follow: (3).

is the volumetric water content; denotes the apparent thermal conductivity of the soil. and are the volumetric heat capacities of the solid and the liquid phases, respectively. is the sink term, and is the Darcian fluid flux density. The apparent thermal conductivity can be expressed with (4) based on the work of de Marsily [32] as well as of Chung and Horton [33]:

where is the thermal dispersivity, while , , and are empirical parameters that can be estimated by using the sand, silt, and clay content of the soil.

Initial conditions of the water flow domain were measured with TDR (IMKO TRIME-FM3) tube access probe in 10 cm increments in three replicates. The initial soil temperature was set to uniform 10°C in the whole profile. HYDRUS-1D requires the setting of boundaries conditions for solving the flow equations. For water flow, atmospheric boundary conditions with surface runoff and free drainage were prescribed at the upper and lower boundary, respectively. For heat flow, the temperature values at both boundaries were provided in the model input file.

The parameters of CERES were calibrated by inverse modeling [34] so that the simulated LAI values would be in good agreement with the observed values (Figure 1). The obtained daily LAI values were used as inputs for HYDRUS-1D, as well. The measured and calculated soil temperature values were compared with simple graphical and statistical tools.

Figure 1: Observed and simulated leaf area index (LAI) values for the 1st parcel (fertilized, nonirrigated) of the experiment.

Originally, the user cannot alter the functions ((1) and (2)) of heat transport in CERES. Though it is an empirical model, it cannot be calibrated. In other words, it is postulated that it works for all soil types. A simple modification of one of its governing equations (1) is proposed to provide greater flexibility and the possibility of site-specific calibration for the following model: (5).

By modifying the value of parameter in (5), the amount of heat reaching the deeper soil layers could be adjusted. Though this parameter has no clear physical meaning, it most likely integrates the effect of soil organic matter, soil structure, and other implicit factors on soil-specific heat.

3. Results

Considerable differences in the leaf area indices were observed in different treatments of the experiment at the end of the canopy development in 2011 (Figure 2). Over 5°C, difference was observed in the daily maximum temperature at 5 cm depth, in the selected parcels. At 20 cm depth the observed difference was still explicitly greater (3.8°C) than the measurement error. The peek temperature at this depth occurred 6 hours later (at sunset) than the peek of the air temperature.

(a)

(b)

Figure 2: Effect of leaf are index (LAI) on the soil temperature dynamics at different depths in a non-fertilized (to the left) and in a fertilized parcel (to the right) on 14/07/2011 at Őrbottyán, Hungary.

The calibration of (5) resulted in for the newly introduced parameter. When the series of measured and calculated temperature data were analyzed, it became obvious that the original CERES considerably underestimated the soil temperature of the deeper soil layers especially in the year 2011 (Figure 3). According to the error indicators (Figure 4), the modified CERES estimated better the than the original CERES especially in the deeper layers. The calculated values of HYDRUS-1D fit the best to the measured data (Figure 4). At the top and bottom layers, the modified CERES presented similar performance indicators than those of HYDRUS-1D.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)

(k)

(l)

(m)

(n)

(o)

Figure 3: Series of measured and calculated soil temperature values at different depths, at Őrbottyán, Hungary. Thick lines—measured, thin lines—calculated.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)

(k)

(l)

(n)

(o)

(p)

Figure 4: Comparison of the measured and calculated soil temperature values at different depths, at Őrbottyán, Hungary.

HYDRUS-1D considerably underestimated the soil temperature of the upper layers in the frosty period when the average air temperature was −7.1°C between 27/12/2010 and 04/01/2011, while the average observed and calculated was −3.4 and −6.5°C, respectively, at 5 cm depth.

During the winter period, all models overestimated the trends of temperature changes and resulted in more pronounced oscillations of soil temperature than that of the observed values.

4. Discussion

The comparison of the course of measured soil temperature at two parcels (different treatments of the experiment) on a summer day highlights the effect of canopy development status on soil temperature dynamics (Figure 2). The vegetation reduces both the average soil temperature and its diurnal amplitude; therefore, considering the LAI is important in modeling.

The measured soil temperature was not below −5°C in the upper layers despite the fact that the average daily air temperature was permanently below −7°C (some days the daily minimum was below −15°C) for longer periods in December 2010 and January 2011. During this period, the site was covered with snow which reduces the effect of freezing since the depth of frost penetration is sensitive to the details of snow cover buildup [35]. As the models underestimated the actual soil temperature and overestimated the temperature oscillation within the winter period, it is obvious that all of them failed to account for the insulation effect of snow cover.

When the series of measured and calculated temperature data were analyzed, it became obvious that the original CERES considerably overestimated the summer soil temperature, while it underestimated during the winter. The minimum of the calculated soil temperature was below −5°C at 60 cm depth, while the corresponding measured temperature was +2°C. The modified CERES let less heat to be transmitted to the deeper layers resulting in lower temperatures. The calculated average error of was reduced with almost 70% compared to the original CERES at 40 and 60 cm depths.

Though HYDRUS-1D gave the most accurate soil temperature estimations, it has to be noted that this model requires more detailed inputs than CERES does. For example, CERES calculates the leaf area development, while HYDRUS-1D requires LAI data to be provided as inputs. Furthermore, HYDRUS-1D requires some thermal properties of the soil (e.g., thermal conductivity) which was estimated from the available textural information in this study. This might explain the relative moderate performance of this model.

5. Conclusions

Two soil temperature models using different concepts were compared in this study. The simpler empirical model was enhanced by introducing an extra parameter in one of its governing equations. The experimental results clearly showed that crop cover significantly influences the soil temperature dynamics of the upper soil layers. Therefore, considering the LAI in model calculation is indispensable. The seasonal snow cover could significantly modify the freezing of soil as it builds up an isolating layer. The simulation of the effect of snow cover should be enhanced in the investigated models. The additional parameter proposed to modify the calculation of the CERES model provided greater flexibility and resulted in better performance, though the comparison was carried out only for a very sandy soil. Further studies should be conducted to investigate the capability of the modified CERES for simulating the heat transport of more structured soils with higher clay and organic matter contents. Though the more sophisticated HYDRUS-1D provided more accurate soil temperature estimations, its superiority to CERES is not unequivocal. The considerable input requirements of HYDRUS-1D may force the users to apply parameter estimation methods which most likely decrease the model accuracy.

Acknowledgment

The authors gratefully acknowledge the financial support of the Hungarian Scientific Research Foundation (K67672).

References

G. Nabi and C. E. Mullins, “Soil temperature dependent growth of cotton seedlings before emergence,” Pedosphere, vol. 18, no. 1, pp. 54–59, 2008.

X. Liu and B. Huang, “Root physiological factors involved in cool-season grass response to high soil temperature,” Environmental and Experimental Botany, vol. 53, no. 3, pp. 233–245, 2005.

S. Dong, C. F. Scagel, L. Cheng, L. H. Fuchigami, and P. T. Rygiewicz, “Soil temperature and plant growth stage influence nitrogen uptake and amino acid concentration of apple during early spring growth,” Tree Physiology, vol. 21, no. 8, pp. 541–547, 2001.

R. D. Boone, K. J. Nadelhoffer, J. D. Canary, and J. P. Kaye, “Roots exert a strong influence on the temperature sensitivity of soil respiration,” Nature, vol. 396, no. 6711, pp. 570–572, 1998.

R. Inclán, D. De La Torre, M. Benito, and A. Rubio, “Soil CO₂ efflux in a mixed pine-oak forest in Valsaín (Central Spain),” TheScientificWorldJournal, vol. 7, no. 1, pp. 166–174, 2007.

G. G. Katul and M. B. Parlange, “Estimation of bare soil evaporation using skin temperature measurements,” Journal of Hydrology, vol. 132, no. 1–4, pp. 91–106, 1992.

G. S. Rahi and R. D. Jensen, “Effect of temperature on soil-water diffusivity,” Geoderma, vol. 14, no. 2, pp. 115–124, 1975.

L. C. Paraíba, A. L. Cerdeira, E. F. Da Silva, J. S. Martins, and H. L. Da Costa Coutinho, “Evaluation of soil temperature effect on herbicide leaching potential into groundwater in the Brazilian Cerrado,” Chemosphere, vol. 53, no. 9, pp. 1087–1095, 2003.

X. Zhou, N. Persaud, D. P. Belesky, and R. B. Clark, “Significance of transients in soil temperature series,” Pedosphere, vol. 17, no. 6, pp. 766–775, 2007.

A. Olness, D. Lopez, D. Archer et al., “Factors affecting microbial formation of nitrate-nitrogen in soil and their effects on fertilizer nitrogen use efficiency,” TheScientificWorldJournal, vol. 1, pp. 122–129, 2001.

J. Pietikäinen, M. Pettersson, and E. Bååth, “Comparison of temperature effects on soil respiration and bacterial and fungal growth rates,” FEMS Microbiology Ecology, vol. 52, no. 1, pp. 49–58, 2005.

K. I. Paul, P. J. Polglase, P. J. Smethurst, A. M. O'Connell, C. J. Carlyle, and P. K. Khanna, “Soil temperature under forests: a simple model for predicting soil temperature under a range of forest types,” Agricultural and Forest Meteorology, vol. 121, no. 3-4, pp. 167–182, 2004.

B. Bond-Lamberty, C. Wang, and S. T. Gower, “Spatiotemporal measurement and modeling of stand-level boreal forest soil temperatures,” Agricultural and Forest Meteorology, vol. 131, no. 1-2, pp. 27–40, 2005.

V. R. Hasfurther and R. D. Burman, “Soil temperature modeling using air temperature as a driving mechanism,” Transactions of the American Society of Agricultural Engineers, vol. 17, no. 1, pp. 78–81, 1974.

D. Zheng, E. R. Hunt Jr., and S. W. Running, “A daily soil temperature model based on air temperature and precipitation for continental applications,” Climate Research, vol. 2, no. 3, pp. 183–191, 1993.

Y. Luo, R. S. Loomis, and T. C. Hsiao, “Simulation of soil temperature in crops,” Agricultural and Forest Meteorology, vol. 61, no. 1-2, pp. 23–38, 1992.

T. J. Toy, A. J. Kuhaida, and B. E. Munson, “The prediction of mean monthly soil temperature from mean monthly air temperature,” Soil Science, vol. 126, pp. 96–104, 1978.

C. H. M. van Bavel and D. I. Hillel, “Calculating potential and actual evaporation from a bare soil surface by simulation of concurrent flow of water and heat,” Agricultural Meteorology, vol. 17, no. 6, pp. 453–476, 1976.

P. J. Wierenga and C. T. De Wit, “Simulation of heat flow in soils,” vol. 34, no. 6, pp. 845–848, 1970.

H.-J. G. Diersch, WASY Software FEFLOW 5.4 Finite Element Subsurface Flow and Transport Simulation System. Reference Manual, WASY Institute for Water Resources Planning and Systems Research Ltd., Berlin, Germany, 2009.

J. Šimůnek, M. T. Van Genuchten, and M. Šejna, “Development and applications of the HYDRUS and STANMOD software packages and related codes,” Vadose Zone Journal, vol. 7, no. 2, pp. 587–600, 2008.

N. Brisson, B. Mary, D. Ripoche et al., “STICS: a generic model for the simulation of crops and their water and nitrogen balances. I. Theory and parameterization applied to wheat and corn,” Agronomie, vol. 18, no. 5-6, pp. 311–346, 1998.

J. W. Jones, G. Hoogenboom, C. H. Porter et al., “The DSSAT cropping system model,” European Journal of Agronomy, vol. 18, no. 3-4, pp. 235–265, 2003.

N. Fodor, R. Sándor, T. Orfanus, L. Lichner, and K. Rajkai, “Evaluation method dependency of measured saturated hydraulic conductivity,” Geoderma, vol. 165, pp. 60–68, 2011.

W. D. Reynolds and D. E. Elrick, “In situ measurement of field-saturated hydraulic conductivity, sorptivity and the alpha-parameter using the Guelph permeameter,” Soil Science, vol. 140, pp. 292–302, 1985.

M. T. van Genuchten, “A closed form equation for predicting the hydraulic conductivity of unsaturated soils,” Soil Science Society of America Journal, vol. 44, no. 5, pp. 892–898, 1980.

N. Fodor and K. Rajkai, “Computer program (SOILarium 1.0) for estimating the physical and hydrophysical properties of soils from other soil characteristics,” Agrokémia és Talajtan, vol. 60, pp. 327–340, 2011.

E. G. Montgomery, “Correlation studies in corn. 24th Annual Report,” Agricultural Experiment Station of Nebraska, pp. 108–159, 1911.

J. T. Ritchie, U. Singh, D. C. Godwin, and W. T. Bowen, “Cereal growth, development and yield,” in Understanding Options for Agricultural Production, G. Y. Tsuji, G. Hoogenboom, and P. K. Thornton, Eds., pp. 79–98, Kluwer Academic Publishers, Dordrecht, The Netherlands, 1998.

P. D. Jamieson, J. R. Porter, J. Goudriaan, J. T. Ritchie, H. Van Keulen, and W. Stol, “A comparison of the models AFRCWHEAT2, CERES-Wheat, Sirius, SUCROS2 and SWHEAT with measurements from wheat grown under drought,” Field Crops Research, vol. 55, no. 1-2, pp. 23–44, 1998.

T. H. Skaggs, T. J. Trout, J. Šimůnek, and P. J. Shouse, “Comparison of HYDRUS-2D simulations of drip irrigation with experimental observations,” Journal of Irrigation and Drainage Engineering, vol. 130, no. 4, pp. 304–310, 2004.

G. de Marsily, Quantitative Hydrogeology, Academic Press, London, UK, 1986.

S. O. Chung and R. Horton, “Soil heat and water flow with a partial surface mulch,” Water Resources Research, vol. 23, no. 12, pp. 2175–2186, 1987.

K. Soetaert and T. Petzoldt, “Inverse modelling, sensitivity and monte carlo analysis in R using package FME,” Journal of Statistical Software, vol. 33, no. 3, pp. 1–28, 2010.

L. E. Goodrich, “The influence of snow cover on the ground thermal regime.,” Canadian Geotechnical Journal, vol. 19, no. 4, pp. 421–432, 1982.