A new approach to predict soil temperature under vegetated surfaces

In this article, the setup and the application of an empirical model, based on Newton’s law of cooling, capable to predict daily mean soil temperature (Tsoil) under vegetated surfaces, is described. The only input variable, necessary to run the model, is a time series of daily mean air temperature. The simulator employs 9 empirical parameters, which were estimated by inverse modeling. The model, which primarily addresses forested sites, incorporates the effect of snow cover and soil freezing on soil temperature. The model was applied to several temperate forest sites, managing the split between Central Europe (Austria) and the United States (Harvard Forest, Massachusetts; Hubbard Brook, New Hampshire), aiming to cover a broad range of site characteristics. Investigated stands differ fundamentally in stand composition, elevation, exposition, annual mean temperature, precipitation regime, as well as in the duration of winter snow cover. At last, to explore the limits of the formulation, the simulator was applied to non-forest sites (Illinois), where soil temperature was recorded under short cut grass. The model was parameterized, specifically to site and measurement depth. After calibration of the model, an evaluation was performed, using ~50 % of the available data. In each case, the simulator was capable to deliver a feasible prediction of soil temperature in the validation time interval. To evaluate the practical suitability of the simulator, the minimum amount of soil temperature point measurements, necessary to yield expedient model performance was determined. In the investigated case 13–20 point observations, uniformly distributed within an 11-year timeframe, have been proven sufficient to yield sound model performance (root mean square error <0.9 °C, Nash–Sutcliffe efficiency >0.97). This makes the model suitable for the application on sites, where the information on soil temperature is discontinuous or scarce.


Introduction
Various biotic, as well as abiotic processes in the soil are temperature dependent (Rankinen et al. 2004). Usually, these dependencies are assumed to have a non-linear nature (Bond-Lamberty et al. 2005;Davidson et al. 2006;Macdonald et al. 1995;Wagle and Kakani 2014), meaning that the response of the process to changes of temperature, strongly depends on the temperature range it is occurring in. Especially for high temperatures, small changes in temperature might yield big changes in the processes response. For the assessment of temperature dependent soil processes, it is therefore crucial to have expedient knowledge about spatial, as well as temporal fluctuations of soil temperature (Bond-Lamberty et al. 2005). The most reliable source of information would be the permanent monitoring of subsurface ground temperature. But in practice it is often hard to measure continuously. Usually, the modeler has to deal with fragmentary timelines of soil temperature, scarce point observations or even no records of T soil at all (Lei et al. 2011). To fill these gaps or to extend the timeline beyond the measurement timeframe, the researcher has to consider the application of a soil temperature model.
The approaches to predict subsurface ground temperature can be coarsely divided in 2 categories; (1) process based models, and (2) empirical models (Kang et al. 2000). Process based approaches to predict soil temperature generally use meteorological input variables (primarily temperature and solar radiation) to calculate energy balance of the soil surface, and heat transport in the soil, by solving the heat equation (Paul et al. 2004). The applicability of these models is often limited by their high complexity, high demand of input data, and specific model parameters, which are often not available for the investigated site (Lei et al. 2011;Svensson et al. 2008). Empirical models, presented in the work of Brown et al. (2000), Kang et al. (2000), or Paul et al. (2004), rely on the statistical relationship between meteorological parameters and soil temperature. More recently, there have been successful attempts to predict T soil using combinations of artificial neural networks and fuzzy logic (Bilgili et al. 2013;Kim and Singh 2014;Kisi et al. 2015;Talaee 2014).
Soil thermal regimes are controlled by various environmental drivers. The most important meteorological factors are air temperature and radiation, laying the base for heat exchange at the soil surface (Hu and Feng 2003). In the latter, forested sites differ substantially from other types of land-cover: The radiation driven heat exchange between soil surface and atmosphere, is limited due to the shielding effect of the canopy (Paul et al. 2004). Therefore, forested sites show strongly dampened T soil fluctuations, compared to sites with sparse vegetation or bare soil (Balisky and Burton 1993). Only a few models exist, which explicitly address the soil thermal conditions of forested ecosystems. Zheng et al. (1993) set up a dynamical T soil model based on Newton's law of cooling, assuming the change of T soil proportional to the temperature difference between air and soil. The fact, that the vegetation cover limits radiation driven heat flux, is taken into account by utilizing a heat transfer coefficient, which depends on the stands leaf area. They assume, that the canopy's damping effect is more pronounced for incoming radiation, than for emission from the ground. This is incorporated, by applying different heat transfer coefficients, whether the soil is warming or cooling. The damping term, dependent on LAI, only comes into effect for soil warming conditions. Based on this work, Kang et al. (2000) set up a spatially resolved T soil model. To describe the soil thermal regimes of South Korean forest sites, they extended the latter approach by introducing a more 'mechanistic' element, based on Fourier's law of heat transport. Besides the spatial and temporal variability of the leaf area, this approach also accounts for the effect of the stands litter layer on soil heat flux. The authors assumed, that T soil does not fall below freezing for most Korean forest sites. As well as in the latter approach, T soil estimates below 0°C were replaced with 0°C. Brown et al. (2000), predicted daily mean T soil of 4 different Northern Hardwood stands, utilizing a statistical relationship between T soil and the average air temperature of the previous day. As a correction term, accounting for the phase shift or 'lagging behind' of the annual course of T soil compared to air temperature, they introduced a cosine function of the Julian day. Despite the simple model structure, the predictions of T soil were quite precise (disregarding the cold season).
To predict daily T soil of various Australian forest sites, Paul et al. (2004) used daily average air temperature and stand parameters like leaf area, understory growth, and litter mass. They assumed T soil oscillating around an annual mean soil temperature, which is calculated from annual mean air temperature, modified with a correction factor, derived from information about the stands' vegetation cover and litter layer. The resulting temperature wave is then offset by a term describing daily fluctuations of T soil , which again, is derived from air temperature. The model specifically addresses the thermal conditions of the topsoil. Therefore, phase shift and attenuation of the temperature oscillation, which become relevant with increasing soil depth, were not considered. Bond-Lamberty et al. (2005) examined the spatiotemporal dynamics of soil thermal regimes during stand development of a disturbed boreal forest. To accompany this investigation and for laying the base to simulate forest dynamics, they implemented an empirical T soil model. Accounting for the influence of recent past air temperature conditions on present T soil , they calculate running averages of the daily mean air temperature. T soil is then calculated as a linear function of multiple running averages, centered to different days in the past. The authors report difficulties to predict T soil close to the freeze/thaw transition.
To evaluate the suitability of T soil as a predictor for the treeline position in the Swiss Alps, Gehrig-Fasel et al. (2008) presented an approach, which strongly differs from others described in this section. To satisfy the statistical requirements for regression modeling, the data was first detrended and then transformed for first differences. After performing the regression analysis, the data was transformed back. Considering that daily mean air temperature was the only input parameter, the model showed high performance in the validation timeframe. Assuming only an insignificant influence of winter soil temperatures on the treeline position (Körner and Paulsen 2004), the validation could be limited to the warm season.
Most approaches presented here disregard T soil dynamics of the cold season. The decoupling of the soil from the atmosphere by a fluctuating snowpack (Betts et al. 2001), the heat transformation processes at the phase change from liquid to frozen (Beltrami 2001;Viterbo et al. 1999), or changes in heat capacity and conductivity seem difficult to be captured in the framework of an empirical approach. In cases where winter T soil is assumed to reach or fall below 0°C, process based approaches, presented by e.g. Rankinen et al. (2004), should be preferred. But, even though this model could be described 'simple' from a mechanistic point of view, solely the empirical snow accumulation/melt module, upstream to the T soil model, requires the assignment of 11 free parameters. An alternative might be the semi-empirical model presented by Katterer and Andren (2009). Making the approach suitable for colder temperature conditions, the formulation presented by Kang et al. (2000) was modified. They interposed a surface temperature term, which acts as link between air and soil temperature. In this term the influence of air temperatures below 0°C is attenuated by a constant factor. This way, they account for the low thermal conductivity of snow.
The objective of this article is the presentation of a model to predict soil temperature of forest stands, which aims to perform like a 'well-tuned' mechanistic simulator, using the straightforwardness of an empiric formulation. The model enables the transformation of fragmentary records of forest soil temperature, into a complete time series of T soil , using average daily air temperature as only input. In this specific case, the created time series is laying the base for the modeling of temperature dependent, biogeochemical soil processes. Due to the fact that many biotic soil processes are sensitive to winter conditions (Campbell et al. 2005), emphasis is laid on an expedient representation of the temperature dynamics of the cold season.
Running the simulation requires the adjustment of nine empirical parameters, which are not defined in a strict physical sense. This is making it hard to deduce parameter values directly from site information. For a proper site specific parameterization, at least some snapshot measurements of T soil are recommended. Therefore, this model primarily aims to sites were T soil data is available, but the time series are inconsistent, or have to be extended beyond the timeframe of measurement.

Model description
The model describes T soil as a function of daily mean air temperature (T air,t ). It employs a daily time step. The formulation is based on Newton's law of cooling (Bergman et al. 2011), which is applied 2 times consecutively.
Utilizing a relatively small heat transfer coefficient (k shift ), the first application of Newton's law provides a phase shifted temperature time series (T shift,t ) which lacks the high frequency fluctuations of T air,t .
A fictive environmental temperature (T env,t ) is postulated as the weighted mean of the elements T air,t , T shift,t , and a constant correction temperature (T corr ). pc air , pc shift , and pc corr are partitioning coefficients, which define the relative weight of the specific element.
The partitioning coefficients sum up to 1, so 2 have to be defined as model parameters, one can be deduced.
DT states the difference of the soil temperature to T env,t .
Taking into account the insulating effect of the snow cover and the heat release/consumption due to the phase change of soil water from liquid to solid and vice versa (Beltrami 2001), a variable heat transfer coefficient (k eff ) is implemented (Fig. 1). k max represents the transfer coefficient above the upper threshold temperature (T 1 ). Below T 1 k eff gets reduced, reaching the minimum (k min ) at the lower threshold (T 0 ), where different k min are applied for soil warming and cooling.
The transition of the transfer coefficient in betweenT 1 andT 0 is described, using a third order polynomial. Fig. 1 Polynomial transition of the heat compensation coefficient (k eff ), between 2 threshold soil temperatures (T 0 , T 1 ), close to soil freezing. The reduction of the coefficient pays respect to the energy release/demand of phase changes, from liquid to solid and vice versa. High model performance was achieved, using different minimal compensation coefficients for soil cooling (k frost ) (solid line) and warming (k thaw ) (dashed line) respectively T soil,t-1 has to be transformed into an auxiliary variable inside the interval 0-1.
At last, Newton's law is applied the 2nd time. The actual daily mean soil temperature calculates as: Study sites/input data

Austria
In the framework of the International Co-operative Programme on Assessment and Monitoring of Air Pollution Effects on Forests (ICP Forests), the Austrian Research Centre for Forests operates several, intensively monitored, forest sites (Level II) (Neumann et al. 2001). In addition to various other environmental parameters, meteorological conditions are monitored continuously. Soil temperature records exist for soil depths, ranging from 5 to 60 cm. The model was originally set up on data from the Level II Plot Klausen-Leopoldsdorf, which is located in the Vienna Woods (48°07 0 16 00 N, 16°02 0 52 00 E), at an elevation of 510 m a. s. l. The research site is a pure beech (Fagus sylvatica L.) stand, which was planted in the late thirties of the last century. The location is facing NE with an inclination of 20 %. The actual forest vegetation coincides with the potential natural one, and can be classified as Hordylemo-Fagetum (Mucina et al. 1993).
Subsequently data from 5 other Level II forest stands were accessed (Fig. 2, Table 1). The selection aims to cover a broad range of site characteristics. Investigated sites show a strong altitudinal and climatic gradient. The elevation of the investigated stands ranges from 290 (Unterpullendorf) to 1540 m a.s.l. (Murau), leading to annual mean temperatures from 9.6 to 5°C, respectively. Austria lies in the transition zone between oceanic and continental climate. Progressing from west to east, investigated locations therefore experience a strong decline in annual precipitation sums, ranging from 1521 mm for mountainous stands in the north-west, affected by orographic precipitation (Mondsee), to 630 mm in the continentally influenced east of the country (Unterpullendorf).
To fill gaps in the record of average daily air temperature, data were accessed, provided by the European Climate Assessment (ECA&D) (Tank et al. 2002). Missing values were replaced, using linear regression with available neighboring stations.

East Coast of the United States
Intending to test the models over regional validity, the continent was switched. Data were accessed from 2 Fig. 2 Location of study sites in the United States and in Austria. The sites used for parameterization of the forest soil temperature simulator cover a broad range of characteristics. For a brief site description see Table 1 intensive long-term ecological research areas in New England; (1) the Hubbard Brook Experimental Forest (HBEF), and (2) Harvard Forest (Fig. 2, Table 1).
The HBEF is located in the White Mountain National Forest in north-central New Hampshire (43°56 0 N, 71°42 0 W). The elevation of the investigated watershed ranges from 250 m to 1000 m. The forest type can be classified as Northern Hardwood, dominated by Sugar maple (Acer saccharum Marsh). The climate is cool, continental, and humid, with mean annual precipitation sums around 1400 mm (Bailey et al. 2003). Approximately onethird of the precipitation is falling as snow, leading to a snowpack, typically lasting from December to April. Soils can be classified as well-drained Spodosols (WRB: Podzol), developed on glacial till (Campbell et al. 2010).
Daily T soil (depth: 5 cm) data, recorded within the framework of the project 'Snow Depth & Soil Freezing as a Regulator of Microbial Processes' (Duran et al. 2014), were obtained. Data of 3 intensive high elevation plots (mean elevation: 560 m, exposition: North) and 3 intensive low elevation plots (mean elevation: 430 m, exposition: South) were used. For each altitude class, one mean timeseries of T soil was calculated.
Because of their proximity to the investigated stands, records of air temperature (Bailey et al. 2003) from meteorological station 23 and 1, for high and low elevation plots respectively were obtained. Missing data were replaced, using offset temperatures of highly correlated neighboring stations. To fill remaining gaps in the air temperature record, the GHCN-Daily dataset was accessed, provided by the NOAA (Menne et al. 2012a, b), utilizing data from the station Wentworth, New Hampshire (43°52 0 22 00 N, 71°54 0 31 00 W).
The Harvard Forest Research Station is located in Central Massachusetts (42°32 0 N, 72°11 0 W). The climate is cool, temperate, and humid. Precipitation is distributed evenly through the year, with annual sums in the range of 1080 mm. The annual mean temperature is 8.5°C (Berbeco et al. 2012). The elevation of the investigated locations is approximately 350 m a. s. l. Soils can be classified as Typic Dystrudepts (WRB: Dystric Cambisol). After a severe disturbance in the beginning of the last century, the forest regrew naturally, resulting in an even aged stand of mixed hardwood species, with Red oak (Quercus rubra L.) dominating (Butler et al. 2012).
Within the forest site, the simulator was applied to 2 sub-sites: (1) Barre Woods (Melillo et al. 2003), and (2) Prospect Hill (Melillo et al. 1999). Both locations were set up to study the effect of soil warming on carbon and nitrogen turnover, by artificially heating the ground (Berbeco et al. 2012;Melillo et al. 2002). The model was adjusted to the topsoil (depth 5 cm) of the undisturbed control plots, whereat on the Prospect Hill site data from 6 control plots were combined, calculating a mean timeseries of T soil . Daily air temperature was obtained from the EMS tower (Munger and Wofsy 1999), where the record 7.6 m above ground was selected. Data gaps were closed, using offset temperature measurements at other heights, or data from the Fisher meteorological station (Boose 2001). If no other source was available, the GHCN-Daily dataset was again accessed, applying offset air temperature data from the Municipal Airport station at Orange, Massachusetts (42°33 0 46 00 N, 72°16 0 59 00 W).

Non-forested sites in Illinois
At last, to explore the limits of the formulation, the model was applied to 6 sites which lack the shielding properties of a dense forest canopy. Therefore, data were obtained from the Illinois Climate Network, (ICN), which operates several open field meteorological stations in Illinois (Fig. 2, Table 1). Air temperature was measured 2 m above ground. Gaps in the air temperature record were closed, using offset temperature measurements of, highly correlated, and neighboring stations. Soil temperature was recorded in 10 cm and 20 cm depth (Hollinger et al. 1994) under sod covered ground. Soil texture was assessed as silt loam, throughout all studied locations. The elevation of the investigated sites ranges from 133 to 265 m a. s. l. Illinois' climate is typically continental with cold winters and warm summers. Moving from north to south, mean annual air temperatures increase from 8.9 to 14.5°C. Also annual precipitation sums reveal a strong north-south gradient, ranging from 810 to 1220 mm. Stations in the north-west of the state are climatically influenced by Lake Michigan, which is attenuating temperature extremes and enhancing winter precipitation (lake effect snow) (Changnon et al. 2008).

Parameterization
The model was applied to each site and depth specifically. Emphasis was laid on its application on longest possible records of T soil , to cover the broadest possible range of different environmental states, which might have a potential influence on soil thermal regimes. On the other hand it seems obvious, that due to changes in leaf area, undergrowth, litter layer, water consumption, etc., forest T soil regimes undergo a certain shift during stand development (compare Kang et al. 2000). In cases where, for reasons unknown, an obvious change in the soil thermal regime was observed, the time frame of the investigation was manually narrowed down. Both T soil,t and T shift,t were initialized at 8°C. The simulator ran a 150 day spin-up prior to the analysis time frame. For model parameterization a simulated annealing algorithm (Kirkpatrick et al. 1983) was applied, selecting an exponential cooling schedule. Optimization/evaluation criterion was in every case the Nash-Sutcliffe model efficiency (NSE) (Nash and Sutcliffe 1970).
NSE ¼ 1 À P n i¼1 T soil;obs;i À T soil;sim;i À Á 2 P n i¼1 T soil;obs;i À T soil;obs Enabling a balanced split, the calibration was conducted on data from odd years, data from even years served in the evaluation. Making the simulation result comparable to other works, other performance indices like Root Mean Squared Error (RMSE), mean absolute error (MAE) and mean bias error (MBE) were calculated.
MAE ¼ n À1 X n i¼1 T soil;obs;i À T soil;sim;i ð11Þ

Parameterization on limited input data
To test the simulators practical suitability to cope with limited input data, the T soil record of Klausen Leopoldsdorf (15 cm depth) was used, ranging from November 2001 to June 2013 (*11 years, 4053 valid observations). The dataset was split into n sectors of approximately equal size. The parameterization (simulated annealing) was performed, drawing only one random observation per sector. The remaining observations served in the evaluation. This step was repeated 12 times per n, each time with different random observations, to generate a distributed result. After 12 iterations, n was incremented, starting with n = 4, gradually progressing to n = 2000. This way, the minimum number of point observations was determined, necessary to yield satisfactory model performance.

Results and discussion
The model was applied to various sites and depths. In this work, a representative selection of 36 simulation runs is displayed ( Table 2). The simulator delivered good estimates of T soil on all investigated forest sites. NSE values above 0.979 and RMSE consistently below 1°C underline the outcome (Table 3), whereat good results were not limited to the topmost soil horizons. Increasing phase shift and the attenuation of the temperature wave with increasing soil depth, were also captured by the simulation (Fig. 3b). Winter T soil dynamics are strongly affected by (1) heat transformations at the freeze/thaw transition and (2) the insulating by the snowpack (Beltrami 2001). The presented model does not specifically address these effects, but it is capable, to account for both effects combined. In most cases, the description of the winter soil thermal regime was successful. Figures 3a, d and 4a clearly show the decoupling of ground temperature from air temperature under snow cover. The simulator was able to track this behavior, where in some cases it failed to predict the exact time when soil temperature rises in spring (Fig. 4a): The melting of the snow cover causes a sharp increase in T soil due to the ceasing insulating effect, hand in hand with an abrupt decrease in surface albedo, making the forest ground susceptible for short wave radiation inputs, which are already considerable in early spring. Rankinen et al. (2004) solved this problem by incorporating a snow dynamics routine into the calculations, but this would require the embedding of more model parameters and meteorological input data. In consideration of the models practical applicability, this was set aside.     One year section of observed and simulated T soil time series plus the corresponding performance scatterplot. Note that the scatterplots cover the whole investigation timeframe! (a, b) winter snow cover decouples the course of air and soil temperature. The melting of the snowpack in the end of March causes T soil to escalate, due to the ceasing insulation plus the abrupt decrease in surface albedo, making the soil susceptible for short wave radiation inputs, which are already considerable in early spring. As the snowpack is not modeled explicitly, the simulator fails to predict the exact time when T soil rises in spring (c, d). Failure to predict a major soil frost event, due to limitations in the model structure: Temperature fluctuations in early winter indicate the absence of a snow pack. When in midwinter all latent heat is released due to the freezing of soil water, T soil suddenly drops. In the formulation the transfer coefficient below the lower threshold temperature (T 0 ) remains constant. As a consequence, our formulation applies best, to sites where severe soil frost plays only a subordinate role (e, f) warm locations, located at low elevations, in the south of Illinois (Belleville, Brownstown, Olney). In contrast to forested sites, open field sites, lack the attenuating properties of a dense canopy, or a thick litter layer. Especially for cold, but snow-free winters, these locations were prone to soil frost (Fig. 4c, St. Charles). Temperature fluctuations in early winter indicate the absence of a thick insolating snow pack. When in midwinter all latent heat is released, due to the freezing of soil water, T soil suddenly drops. Due to the structure of the model, this behavior could not be tracked: In the presented formulation the transfer coefficient below the lower threshold temperature (T 0 ) remains constant at a reduced level, suppressing further soil cooling. This model limitation could be tackled by letting the transfer coefficient rise at temperatures below T 0 . On the other hand, that would require the segregation of the effects of freeze/thaw processes and snow cover insulation, making the model again more complex and input data demanding.
The examination, to determine the minimum amount of point observations of soil temperature, necessary to yield suitable results, was performed on, an 11-years time series, of air and soil temperature at the Level II plot Klausen-Leopoldsdorf (15 cm depth). The time frame was divided in n sectors. Only one observation was selected randomly by sector. All other observations served in the evaluation. Disregarding single outlier runs, good results (NSE [0.97, RMSE \ 0.9°C) were achieved with n B 13. Having available 50 or more daily observations, there was only little difference to the result, compared to utilizing *50 % (n = 2000) of the available data in the calibration process (Fig. 5).
Two considerations led to the implementation of decreasing transfer coefficients with decreasing soil temperature: (1) The heat release/consumption at the freeze/ thaw transition (Beltrami 2001), and (2) the insulating effect of the winter snow cover. So intentionally, values for T 0 and T 1 were searched around 0°C. Surprisingly, in most cases the optimization process led to T 1 values much higher, meaning that the attenuation of the transfer coefficient starts already at higher temperatures. The idea behind utilizing different responses for soil warming and cooling, was the assumption, that soil warming in spring is strongly driven by incoming solar radiation, which is accelerating the temperature rise.
As this model is primarily of an empirical nature, used parameters lack a specific meaning, in a strict physical sense. Nevertheless, it was noted that parameter values were strongly affected by certain site characteristics: k max values clearly decreased with increasing soil depth (Fig. 6a). Meaning, the time demand, to compensate a fraction of the temperature difference between soil layer and air, rose with increasing soil depth. Also the relative partition of the correction temperature (pc corr ), in the calculation of the environmental temperature, increased in deeper soil layers (Fig. 6c). In contrast, the direct influence of air temperature (pc air ) showed a decrease downwards.
Investigated open field sites differed strongly from forest sites, in parameter values of the correction temperature (T corr ). Where on forest locations T corr resided closely to the stands annual mean air temperature, open field sites revealed T corr values, around and above 100°C (Fig. 6d). On the other hand, their relative weight (pc corr ) in the calculation of the environmental temperature, is much lower than on forested sites. It is assumed, that in these cases, they correct for direct radiation energy inputs, which are obviously much higher without the presence of a shielding canopy. The reason that, even under such conditions, the simulator (which does not particularly address radiative heat flux) delivers good estimates of T soil , might be found in the strong correlation between energy balance components, and the air temperature itself (Hock 2003).

Conclusion
The primary intention of this work was the provision of a tool, which enables the transformation of fragmentary records of forest soil temperature, into a complete time series of T soil , using average daily air temperature as only input. In this specific case, the created time series is laying the base for the modeling of temperature dependent, biogeochemical soil processes.
To test the resilience of this model, it was applied to various locations and depths, covering a broad amplitude of site characteristics. The simulator delivered accurate predictions of the temperature of the topsoil, as well as of deeper layers. The high performance was not limited to the warm season. The combination of the insulating effect of the snow cover plus the effect of heat transformations at the freeze/thaw transition, on soil thermal regimes were captured sufficiently. The formulation was applied to forested, as well as open to field locations, where in the open field it failed to reproduce some major soil frost events. Bearing this limitation in mind, this simulator seems to be well applicable to other land use types. b Fig. 6 Four selected parameters and their change with increasing soil depth. a Transfer coefficient values showed a clear decreasing trend with increasing soil temp. b Also the fraction of the air temperature in the calculation of the environmental temperature showed, almost linear, decrement. d Where the correction temperature on forested sites was in a close range to the annual mean air temperature, the open field locations (star symbol) revealed much higher values. On the other hand, the relative weighting (c) of these temperatures was much smaller on non-forested sites. It is assumed, that on these locations, both parameters combined compensate for direct shortwave radiation inputs The model parameters lack a specific meaning in a strict physical sense. Therefore, currently the parameterization requires at least a modest amount of T soil observations, to yield sufficient results. A challenging impulse for future work, would be the attempt to derive model parameters, directly from more easily obtainable site characteristics. This also would enable the capability of the simulator to deal with a changing soil thermal regime, during stand development.