Microclimate buffering of winter temperatures by pine stumps in a temperate forest

In the presented study we evaluated the role of wood in buffering its inhabitants against winter temperatures. We followed the microscale temperature patterns inside and outside decaying pine stumps during two winter periods. We analyzed the data by extracting the minimum, maximum and mean temperatures for the episodes of snow, frost or no frost. We compared the temperature variation by applying an absolute values calculation (=modulus) for the subsequent measurements. Finally, we tested the buffering effect of the horizontal depth inside the stumps by comparing the pattern for 5 and 15 cm depths. The results show that (1) the minimum temperature was significantly higher inside than outside the stumps, while temperature variation tended to be larger outside than inside, (2) the dynamics of temperature varied between weather episodes, i.e. the periods characterized by the presence/absence of snow and of frost, (3) the minimum temperature inside the stumps increased with the horizontal depth, while temperature variation decreased. The results presented are the first truly microclimatic data on wood buffering and they are important for the understanding of the possible effects of climate change on the thermal relationships at the microscale.


Introduction
The role of temperature in different aspects of the lives of organisms is a widely discussed issue (Atkinson and Thorndyke 2001;Johnston and Bennett 2008) and becomes even more important in the face of global warming (Helmuth et al. 2005;Dillon et al. 2010;Salinger 2013). Microclimate is an important driver of performance in organisms; recent studies have focused on habitat response to thermal changes at the microscale level (Helmuth 2009;Helmuth and Hofmann 2001;Helmuth et al. 2010;Suggitt et al. 2011;Potter et al. 2013). The microscale environment was found to affect the physiology of mussels in intertidal habitats (Helmuth 1998(Helmuth , 2009, the population dynamics of herbivores in the canopy (Pincebourde et al. 2007;Pincebourde and Woods 2012), and the critical thermal maxima for frog and lizard species in a primary rainforest in the Philippines (Scheffers et al. 2014). Assessing selective pressures in a set of different microhabitats may provide a wider understanding of the general effect of microscale on performance. The present study contributes to this issue by examining the cold buffering effect in decaying wood habitat.
Concerning the possible response to climatic change, the lower end of the thermal continuum is a target for natural selection modulating physiological adaptation to low temperatures, while tolerance to high temperature is more constrained physiologically (Araújo et al. 2013). Since global temperature increase is driven by the rise of daily minimum rather than daily maximum temperatures (Easterling et al. 1997;Vose et al. 2005;Brown et al. 2008), a stronger organismal reaction is expected at minima rather than at maxima of temperature tolerance.
Tree stumps may serve as temporary or permanent habitat for arthropods of different feeding guilds such as woodfeeders, fungivores and predators. At least in the case of wood-feeders, this type of habitat is evolutionarily demanding because it requires various mechanisms for dealing with nutritionally poor food (Haack and Slansky 1987). On the other hand, it is remarkably safe; the mortality of wood-feeding larvae is lower than the mortality of larvae from other feeding guilds (Walczynska 2010). Interestingly, the life strategy of a species changes with the horizontal depth of its wood boring niche, likely caused by the heterogeneous environment inside the tree (Walczynska 2012). Although it was shown that wood inhabitants are characterized by a relatively low mortality rate (Haack and Slansky 1987;Walczynska 2010), it is unknown whether this trend is driven by low pathogen and predatory pressure or by high winter survival. Predictions from a theoretical model suggest that cold-induced mortality is one of the key factors driving the population dynamics of the mountain pine beetle Dendroctonus ponderosae (Régnière and Bentz 2007), yet our understanding of the factors shaping the life-history strategies of wood inhabitants is poor. In general, overwintering was advocated as a key component of the biology of organisms inhabiting temperate, polar and alpine regions (Williams et al. 2015).
The role of wood in buffering outside winter temperature has previously been investigated mostly in the context of insect pest species, especially the pine beetle. Therefore, most existing studies examined living pine trees, and temperature data were recorded under the bark of trees. The results are contradictory; some studies showed that the buffering effect of wood can be small (<1 °C difference for an annual study; Bolstad et al. 1997), or considerable (up to 12 °C difference for a period of sub-zero temperatures; Beal 1934, after Tran et al. 2007; that it may have a constant effect (Bolstad et al. 1997), or may depend on other environmental factors, as in a two-year study on 60 ash trees by Vermunt et al. (2012). Potter and Andresen (2002) designed a theoretical model of heat flow throughout the tree. However, this model cannot be tested without reliable data on air-and inside-wood temperatures that are generally unavailable. Therefore, more environmental data are needed in this context.
In this study we used dataloggers with two temperature sensors to simultaneously measure the temperature inside and outside partly decayed pine stumps. We monitored temperature changes through two winter seasons. We distinguished 'episodes' within measurement seasons based on the presence/absence of snow and of frost to control for their potentially different effects on wood buffering. For comparison we used different measures of temperature, including extreme and mean values, and a variability measure-a 'temperature modulus' which is a simple statistic based on the absolute differences between consecutive measurements of two thermal datasets. The advantage of this approach is the linearization of the complicated pattern of the rate of temperature change and a straightforward interpretation. Finally, we compared winter temperature patterns at two levels of horizontal depth within the stumps to examine whether wood inhabitants occupying diverse niches (phloem-vs xylem-feeders) may possibly face different selection pressures during winter.

Study site
We measured the temperature inside and outside decaying pine stumps, remaining after felling in 2000, in the 2004/2005 and 2005/2006 winter seasons (different stumps in each season), in the Niepołomice Forest near Krakow, Poland (50°07′N, 20°23′E). The study site was the same as used previously to investigate the life history strategy of the xylem-feeder Aredolpona rubra (Coleoptera: Cerambycidae; Walczynska 2012), when this area was examined for the presence of larvae and adults (Walczynska 2008(Walczynska , 2010. Moreover, this area was inhabited by other beetle species living in decayed wood e.g., Criocephalus rusticus and Spondylis buprestoides.

Study design
Four stumps (two per winter season), about 60 cm in diameter and 20 cm in height, were chosen for the measurements. We recorded the temperature data from 23 November 2004 to 22 February 2005 (for one stump) and to 18 March 2005 (for the other stump), and from 16 December 2005 to 29 April 2006 for the next two stumps. The scheme of the experimental set-up is shown in Fig. 1. We programmed the HOBO H8 dataloggers (Onset, Pocasset, USA) with an accuracy of ±0.7 °C and a resolution of 0.4 °C at 21 °C for hourly measurements. The inside temperature was measured with the external sensors of the dataloggers (ø 4.7 mm), which were placed inside the stump horizontally in holes ø 5 mm in diameter, drilled into the stump to depths of about 5 (in both seasons) and 15 cm (only in 2005/2006 season) at half of the stump height. The outside temperature was measured with the internal sensors of the dataloggers, which were placed on the soil surface close to the stumps (c. 2 cm). Sensors were delicately covered with the plant litter. Because the sensors were under litter or relatively deep in the stump wood, and the stumps were shaded for most of the time, the effect of solar radiation on the measurements can be neglected.

Data organization
We divided the datasets for each winter season into episodes of frost, snow and no frost. Firstly, we categorized days with snow cover as snow episode (S), days without snow cover when the temperature dropped below 0 °C for 24 h for at least one sensor as frost episode (F), and the remaining days as no-frost episode (NF). For each day we isolated the minimum and maximum temperatures separately for inside and for outside. For further analysis we calculated the means for the three lowest and three highest daily values, respectively, as T min and T max , and the mean daily values (T mean ). Finally, we computed the mean daily absolute values between hourly measurements (daily modulus) according to the formula: where T(n)-temperature at time n; T(n − 1)-temperature 1 h before; x-a positive (absolute) value of positive or negative difference between two consecutive hourly measurements. The modulus served as a measure of thermal variability because, contrary to statistical variability measures, it has a straightforward biological interpretation.

Wood buffering at 5 cm horizontal depth
We analysed the variables T min , T max , T mean and modulus for the relationship between inside and outside stumps temperatures in separate analyses. Due to the autocorrelation of all daily measures, the statistical approach was to calculate the means for each episode within each winter season independently; three lowest and three highest values for each episode as T min or T max , respectively, and the means of daily values in the case of T mean and modulus. By such an approach in calculating T min and T max , we obtained the seasonal extreme values. We analyzed these data with Repeated Measures ANOVA (SAS, SAS Institute Inc., v.9.4), with the stump as subject, the season as an independent between-subject factor and the episode (three levels) and the location (inside-outside; two levels) as withinsubject factors.

Horizontal depth effect
To estimate the effect of horizontal depth of the stumps (position of a thermocouple) on the buffering of winter temperatures we compared the T min and modulus inside the stumps for two horizontal depth levels: 5 and 15 cm. We tested this effect only in one season (2005/2006) and only for days without snow cover (32 consecutive days since 27 March 2006). We inspected the data for time autocorrelation and corrected them (Yule-Walker method, PROC AUTOREG, SAS Institute Inc., v.9.4). We analysed the corrected data using a GLM model with the depth and stumps as fixed factors and the respective values outside the stumps (T min or modulus) as a covariate (Statistica 10, StatSoft).

Climatic data
The following climatic data were additionally used in the analyses: daily snow cover and mean temperature, for the time periods respective to our study, and yearly snow cover and minimum temperature, for the period of 1986-2015. All these data derived from the Łazy Meteorological Station (Institute of Geography and Spatial Management, Jagiellonian University) situated about 14.5 km south-east of the study site.

Results
In the 2004/2005 season, the temperature pattern was much more variable (considerable fluctuations) than in the 2005/2006 season (Fig. 2). As the differences might depend on the snow cover, we compared the daily temperature data with the snow cover data.  The comparison of daily means for air temperature (climatic data from the meteorological station) and for our "outside temperature" (sensor on the soil, lightly covered with litter), for the respective periods of time, shows that (1) the general pattern of the frequency and direction of temperature changes is very similar in the two datasets in the absence of snow (in the first half of the 2004/2005 season), and (2) the magnitude of these changes is considerably larger for air temperature (Fig. 3). This effect is discussed below.

Wood buffering at 5 cm horizontal depth
The temperature changes inside the stumps, calculated as absolute values of changes between the hourly measurements, are not only smaller than outside, but also show a more stable and constant pattern ("Appendix"). Because the thermal scale was shifted between winter seasons, the temperature data were standardized prior to the analyses by recalculating to the units of standard deviation, for each season independently. The results of the Repeated Measures ANOVA (Table 1) show a consistent and considerable effect of episode factor and its interaction with the season. The temperature inside and outside the stumps differed statistically in the case of T min . Minimum temperature was lower outside the stumps and this pattern changed with Fig. 2 The mean daily temperatures inside and outside of pine stumps at 5 cm depth in the Niepołomice Forest (left Y axis) and daily snow cover from Łazy meteorological station (right Y axis) episode and season (Table 1, Fig. 4; exact temperature values are given above). A tendency (p < 0.1) towards differences in location was also observed in modulus, which was generally higher outside, but was most probably masked by the considerable between-stump variation in temperature changes outside of the stumps (Fig. 4). Neither T max nor T mean differed in temperature regarding datalogger location (Table 1). However, in the case of T mean , location significantly interacted with season, episode, as well as with both (Table 1; Fig. 4). Overall, the pattern found was driven by differences among seasons, especially by the very short (3 day) frost episode in the 2005/2006 winter.

Horizontal depth effect
The horizontal depth of the internal measurements, analysed for the period without snow cover, with the lowest temperature recorded equal to −0.61 °C, affected both the minimum temperature and the temperature variation. In the former case, the internal temperature was affected by depth (p = 0.018), with the external minimum temperature as a covariate (p < 0.001). Overall, the minimum temperature inside the stumps measured at a horizontal depth of 15 cm and adjusted for the outside temperature was 4.6 °C and was significantly higher than at 5 cm depth (adjusted, 3.2 °C). The temperature variation (modulus) inside the stumps was affected by depth (p < 0.001), stumps (p = 0.001) with the external variation as a covariate (p < 0.001). Temperature changes inside the stumps were significantly higher at 5 cm horizontal depth (adjusted for outside conditions, 8.2) compared to 15 cm (adjusted, 7.1).

Discussion
The results of the temperature measurements inside and outside (in shallow litter) decaying pine stumps in Niepołomice Forest (Poland) during two winter seasons show that (1) a buffering effect in winter operated against minimum values (Table 1; Fig. 4), and (2) that the horizontal depth inside the stumps analysed for the period without snow cover affected both the minimum temperature and variation in temperature. Overall, the minimum outside temperature in shallow litter was about 1 °C lower than inside the stumps at 5 cm horizontal depth. Additionally, we observed a tendency for higher thermal variability (about 1.7 times overall) outside versus inside the stumps. The maximum temperature was not significantly different with respect to sensor location. However, it tended to be consistently higher outside (Fig. 4). This slight difference was most probably not statistically detected because of the small sample size. The complicated temperature dynamics are best illustrated in the case of T mean , which, though not different (as it was standardized) among seasons or sensor locations, significantly Fig. 3 The comparison of daily means for air temperature (meteo data) and for outside stump microclimatic data (our data, mean of two stumps) during the period of temperature measurements  Fig. 4). The representative value of our measurements period within the space of 30 years is shown in Fig. 5. According to this illustration, the winters in our study were relatively cold and relatively snowy. Mainly minimum and maximum temperatures were analysed in similar studies. While the results on the buffering effect of wood against minimum winter temperatures are unequivocal (Bolstad et al. 1997;Tran et al. 2007;Vermunt et al. 2012), maximum temperature values were found to be higher inside (Bolstad et al. 1997;Vermunt et al. 2012), or higher outside the trees (Tran et al. 2007). Apart from examining the extreme and mean values, we also took into account the general thermal variability pattern inside and outside stumps. This important issue has not been examined before, although some theoretical models of insect overwintering strategies showed that including the daily thermal variability statistically improved the model fit (Régnière and Bentz 2007).
The main difference between our approach and previous investigations is that the air temperature in the latter was measured in weather stations close to the study areas. In consequence, comparisons between inside-tree and air temperature might have been inaccurate because of temporal and spatial measurement distance and possibly because of the different resolutions of the devices used. The magnitude of difference in comparing inside-outside data between a microscale and a larger spatial scale is best illustrated in Fig. 3, showing air temperature data derived from the nearest meteorological station contrasted with our outside-stump temperature data. This comparison validates (1) the reliability of our "outside" measurements, because the general pattern of the frequency and direction of temperature variation is very similar for the two datasets in the absence of snow (compare with Fig. 2), (2) the importance of taking the snow cover into account, and (3) the truly microclimatic meaning of our results. The difference in Black dots correspond to the measurement seasons in this study (note that the values presented refer to "years", not to "winters") the magnitude of changes between air and outside stumps temperature may be caused by the contact of our external sensor with the ground or with litter cover, but most probably it is the effect of the buffering by the forest itself. This result strongly supports the recommendation that a reliable study on the buffering effect in a given habitat demands a microscale. This is important if an attempt is made to link temperature changes in time with possible effects on organism growth or survival; in the case of our study, the difference in temperature between air and the wood microsite combines the effects of forest, snow, frost and niche occupied within wood. It is also worth mentioning that wind appeared to be the main factor affecting soil temperature pattern in a study by Ashcroft et al. (2009).
The results of our study, attempting a comprehensive assessment of winter temperature dynamics inside decaying stumps, have important consequences for the overwintering strategies of wood inhabitants. The division into winter episodes supports the already well-known dominating effect of snow (Sinclair 2001;Bale and Hayward 2010). This effect may be critical for woodfeeders exploiting the tree at different heights, such as Marmara sp. species (Lepidoptera: Gracillariidae) or three bark beetles (Coleoptera: Curculionidae): Dendroctonus valens, Ips pini and Pityogenes hopkinsi (Price 1984). Distinctive cold-hardening strategies may be expected in such cases. The interacting effect of temperature and snow in determination of microclimatic conditions and, further, of overwintering strategies of organisms was recently discussed by Williams et al. (2015). Moreover, we show the importance of the no-frost episode, which may interrupt frost and snow episodes. Such interruptions occurred in the 2004/2005 season in our study, and their existence may be an indication of climate warming (e.g. Sinclair 2001). The difference between inside and outside temperatures, for all the variables tested, was the most distinct in this type of episode (Fig. 4). Warm periods during winter may substantially affect the survival of overwintering organisms. A freeze-tolerant organism would be exposed to more freeze-thaw cycles, increasing mortality (Sinclair 2001;Bale and Hayward 2010). On the other hand, a chill-sensitive organism not acclimated to low temperatures, such as an invasive species, would endure cold injury and survive better (Nedvȇd et al. 1998). Considering the buffering effect of decayed wood against minimum temperature and possibly against temperature variation, which we report in this study, we conclude that in the face of climate warming the inside-stump community will be less affected (or at a slower rate) by the negative consequences for native species and by invasions of chill-sensitive foreign species, than the shallow litter community.
In previous studies, the buffering effect of wood was estimated under the bark. Here, we present for the first time a comparison of the temperature inside wood measured at two levels of horizontal depth. The results show that in the absence of snow, deeper parts of the stumps are better protected against both minimum temperature and thermal variability. This conclusion has important consequences for organisms occupying different inside-wood niches; diverse overwintering strategies might be expected in such cases. Sinclair et al. (2003) specified the aspects of the winter environment that affect winter survival. These include extreme temperatures, the rates at which they are approached, and the fluctuation patterns. According to our results, at least the first factor, and possibly the last, are buffered inside decaying wood in comparison to the surrounding shallow litter. This might have consequences for the entire inside-wood community because this type of environment prevents the migration of its inhabitants, especially among strict wood-feeders. The adjustment to external temperatures by behaviour (migration) may be facilitated in shallow litter. For example, beetle species may show different abilities to vertical migration to deeper soil in response to external temperature (Hoshikawa et al. 1988).
Finally, our results contribute to the issue of climate refugia, created by the heterogeneity of temperature pattern at the landscape scale (Ashcroft et al. 2009). Decaying stumps with their relatively mild and stable conditions may act as archipelagos on the cold forest ocean for organisms not obligatory related to wood such as spiders, predaceous beetles, millipedes etc. Incidentally, this constitutes another argument in favour of saving wooden debris in forests of economic importance.
In summary, arthropods overwintering inside tree stumps are well-protected against minimum winter temperatures and possibly against high temperature variability. The horizontal depth in stumps affects the general pattern; deeper parts of the stumps are better protected against lowest temperatures and against temperature variability compared to shallower parts of the stumps. This suggests that phloem-and xylem-feeders face different temperature regimes during winter, at least during snow-free intervals. Finally, the periods with and without frost or snow have different temperature dynamics, and this issue opens a new line of inquiry for testing overwintering strategies. Fig. 6 The absolute values of the differences between the successive temperature values for inside at 5 cm depth (left column) and outside (right column) the pine stumps. For the sake of clarity of the figures, data where the temperature change was 0 for all four datasets within a season were omitted