Advanced
Variations of Soil Temperatures in Winter and Spring at a High Elevation Area (Boulder, Colorado)
Variations of Soil Temperatures in Winter and Spring at a High Elevation Area (Boulder, Colorado)
Journal of Soil and Groundwater Environment. 2015. Oct, 20(5): 16-25
Copyright © 2015, Korean Society of Soil and Groundwater Environment
  • Received : July 08, 2015
  • Accepted : October 05, 2015
  • Published : October 31, 2015
Download
PDF
e-PUB
PubReader
PPT
Export by style
Article
Author
Metrics
Cited by
TagCloud
About the Authors
Jin-Yong, Lee
Department of Geology, Kangwon National University
hydrolee@kangwon.ac.kr
Hyoun Soo, Lim
Department of Geological Sciences, Pusan National University
Ho Il, Yoon
Korea Polar Research Institute
Poongsung, Kim
Department of Geology, Kangwon National University

Abstract
The City of Boulder is located at an average elevation of 1,655 m (5,430 feet), the foothills of the Rocky Mountains in Colorado. Its daily air temperature is much varying and snow is very frequent and heavy even in spring. This paper examines characteristics of shallow (surface and depth = 10 cm) soil temperatures measured from January to May 2015 in the high elevation city Boulder, Colorado. The surface soil temperature quickly responded to the air temperature with the strongest periodicity of 1 day while the subsurface soil temperatures showed a less correlation and delayed response with that. The short-time Fourier of the soil temperatures uncovered their very low frequencies characteristics in heavy snow days while it revealed high frequencies of their variations in warm spring season. The daily minimum air temperature exhibited high cross-correlations with the soil temperatures without lags unlike the maximum air temperature, which is derived from its higher and longer auto-correlation and stronger spectrums of low frequencies than the maximum air temperature. The snow depth showed an inverse relationship with the soil temperature variations due to snow's low thermal conductivity and high albedo. Multiple regression for the soil temperatures using the air temperature and snow depth presented its predicting possibility of them even though the multiple r 2 of the regression is not that much satisfactory (r 2 = 0.35-0.64).
Keywords
1. Introduction
The soil temperatures are affected by a variety of natural and anthropogenic factors including solar radiation, wind, vegetation, rain, snow cover, moisture content, land cover, buildings, etc. ( Davidoff and Selim, 1988 ; Zhang et al., 2008 ; Kim et al., 2015 ; Teubner et al., 2015 ; Yoshioka et al., 2015 ) and the affecting magnitude of each factor can be varying with location (longitude, latitude, and topographic elevation) and interconnected surrounding environments ( Wundram et al., 2010 ). The varying soil temperatures derived from changing climate, influence nutrient supply, soil respiration, and productivity, hence consequently the biological diversity of the soil environment ( Chapin III et al., 2000 ; Rustad et al., 2000 ). The changing soil environment eventually affects human life because the thin soil layer of the Earth provides nearly all of biomass production, and in turn impacts the changing climate ( Nikolaidis and Ragnarsdottir, 2015 ).
The alpine or high elevation areas have some different weather characteristics from those of low lands. Drastic change in daily weather is quite common and daily air temperature variation (maximum-minimum) is quite large (on average over 15℃ in Boulder). It is severely cold in winter and unbearably hot in summer. When it rains, it is subject to be torrential and snow can occur even in summer. In the meanwhile, the impact of changing climate is more profound in alpine and high elevation areas compared with other areas ( Beniston et al., 1997 ; Beniston, 2006 ). Thus, it is essential to understand response and changes of soil temperatures responding to outer natural factors like air temperature and snow in such high elevation areas to adapt and mitigate the changing climate.
The objective of this study is to unravel the temporal and spectral characteristics of soil temperatures measured in a highland city, Boulder, Colorado, USA. For the purpose, the soil temperatures using automatic temperature sensors were monitored and then analyzed using time and frequency domain functions, and multiple regression functions.
2. Methods and Materials
- 2.1. Study Area
The study area is located in the city of Boulder, Colorado, USA ( Fig. 1a ) and it is situated at the foothills of the Rocky Mountains. The mean elevation is 1,655 m and the city is home to the University of Colorado (CU), the state’s largest university. The city covers 66.5 km 2 and its population is 102,420 as of 2014 ( City of Boulder, 2014 ). The highest air temperature was 36.7℃ in 2013 and the lowest was −37.8℃ in 1975 ( NOAA, 2015 ). Generally, July is the hottest (monthly mean = 24℃) and February is the coldest (monthly mean = 0℃). The average annual precipitation for 1893-2013 is 4,823 mm, much of it occurs in April and May ( NOAA, 2015 ), and occasionally causes flash floods in the city ( Morss et al., 2015 ). Mean annual snowfall for 1897-2014 is 199.76 cm (snow depth, not water equivalent) and snow occurs in most of months only except July and August ( NOAA, 2015 ). The freezing index of this city is 368 °F·days ( NOAA, 2015 ).
PPT Slide
Lager Image
(a) Location of the City of Boulder, (b) Soil temperature measuring point near the Boulder Creek, (c) Soil temperature logger, and (d) Scene of a snow day (February 28, 2015) near the measuring point.
There are two great geologic provinces in the Boulder area; the eastern one is the Great Plains and the western one is the Rock Mountains ( BASIN, 2015 ). The main city area is in the flatland, east of the Rockies. The geology of the city is mainly comprised of sedimentary rocks, from Lyons Sandstone of Permian age in the west to Pierre Shale of Cretaceous age in the east, and on the top of the Shale, there is the gravel cap in the Boulder valley (main city area), from which the name of the city is originated ( Fennenman, 1905 ; Baker, 1974 ; Bilodeau et al., 1987 ). Especially, the sandstones of Fountain Formation of Pennsylvanian age (called Flatirons) colored in red mainly due to hematite weathering ( Wells, 1967 ; Baker, 1974 ; Bilodeau et al., 1987 ; BASIN, 2015 ) feature a spectacular landscape and attract many travelers from all over the world.
According to Rigg(1993) the soil properties around the study area in the city are as follows: soil moisture 17.19-20.38%, pore space 29.99-32.04%, pH 6.13-7.42, organic matter 5.91-8%, bulk density 0.81-1 g/cm 3 , NH 4 0.15-0.51 meq/100 g, and SO 4 0.06-0.22 meq/100 g. The soil textures are sand 68.15-83.36%, silt 13.95-20.18%, and clay 2.69-13.45% ( Rigg, 1993 ).
- 2.2. Data Collection
The soil temperatures were monitored every hour in the campus of CU Boulder ( Fig. 1b ) from January 12 to May 31, 2015. The three points were selected in the backyard of Faculty Staff Court and they formed a triangle shape with distances of 1 m each other. One waterproof temperature sensor (HOBO temperature loggers, UA-001-64, ONSET, USA) ( Fig. 1c ) was buried at surface soil (1 cm depth) and two at 10 cm depths (named 10R (right) and 10L (left)). The sensor can measure the temperatures from −20 to 70℃ with accuracy of ±0.53℃ and resolution of 0.14℃. I also measured air temperature on site using a commercial thermometer and the snow depth using a steel ruler (see Figs. 1b and 1d ) for January 12 to March 31, every day around 3 pm. Resolutions of the thermometer and ruler are 1℃ and 1 mm, respectively.
Daily air temperature and snow depth data were collected for the period of January 1 to May 31 from a NOAA observatory ( http://www.esrl.noaa.gov/psd/boulder/ ), 2 km east from the soil temperature measuring points. The NOAA provides daily maximum and minimum air temperatures data (daily mean values are not provided) and maximum snow depth. These data were compared with the manually measured air temperature and snow depth data on site. There is a stream (Boulder Creek), 20 m apart from the measuring points, running west to east (see Fig. 1b ). Occasionally the stream water temperatures were also measured for some reference using the thermometer.
- 2.3. Time and Frequency Domains Analysis
In this study, in order to unravel variation characteristics of air and soil temperatures, some representative time and frequency functions were used, including autocorrelation, spectral density, cross-correlation, and short-time Fourier ( Larocque et al., 1998 ; Lee and Lee, 2000 ; Kim et al., 2005 ). Details on these functions can be found in many literatures including Duffy and Gelhar(1986) , Box et al.(1994) , and Padilla and Pulido-Bosch(1995) . In the short-time Fourier analysis, window and its size were Welch and 32 ( Hammer, 2012 ), respectively. The soil temperature is the main focus of this study, so to model the soil temperature, a multiple regression method was applied ( Zhang et al., 1997 ; Adlam et al., 2010 ) using air temperature and snow depth. All these calculations were done with the help of Hammer (2012) . Before the time and spectral analyses, trends were removed from the original time series data using the liner regression ( Lim et al., 2011 ; Hammer, 2012 ).
3. Results and Discussion
- 3.1. Weather Conditions
Weather data for January 1 to May 31, 2015 are presented in Fig. 2 . During the period, the maximum daily air temperature (from NOAA data) ranged −7.8 to 30℃ (mean = 14.1℃) while the minimum air temperature was between −18.3 and 12.2℃ (mean = −0.3℃). Around end of February, the coldest air temperature was recorded. The daily variation of air temperature (maximum-minimum) was large, ranging from 2.2 to 30℃ (mean = 14.5℃). The manually measured on site air temperature data are generally between the maximum and minimum temperatures ( Fig. 2a ) and their Pearson correlation coefficient are 0.88 and 0.72, respectively. It is considered the NOAA air temperature data (2 km apart from the site) represent well the air temperature condition on site and thus they were used in the further analysis instead of the shorter manually measured air temperature data.
PPT Slide
Lager Image
(a) Daily maximum and minimum of air temperatures (January 1 to May 31, 2015) measured at a NOAA site, 2 km east from the soil temperature measuring point where air temperature and nearby stream water temperatures were also manually measured (January 12 to March 31), and (b) Daily snowfall and snow depth measured at the NOAA site, also showing the snow depth manually measured (using a ruler) at the soil temperature point.
Fig. 2b shows the daily snowfall and maximum snow depth for the period. The daily snowfall ranged from 0 to 241.3 mm/day with a mean of 13.4 mm/day. For 151 days of monitoring, snow occurred in 27 days (17.9%). Coincident with the coldest weather in late February (see Fig. 2a ), a large quantity of snow was also recorded. Among the total of precipitation (2,455.4mm) during the period, snow accounts for 82.7% (2,029.5mm). The snow depth in the same period ranged from 0 to 35.56 cm with a mean of 3.89 cm, which is very similar to that manually measured on site (r = 0.94 at p = 0.00). So it is can be said that the weather conditions at the two sites (the soil temperature monitoring site and the NOAA site, 2 km apart) are very similar.
- 3.2. Time and Spectrum Characteristics of Soil Temperatures
Fig. 3 shows the recorded soil temperatures at surface (a), depth of 10 cm (b and c), and temperature difference between at surface and depth of 10 cm (d and e). The surface soil temperature was between −1.5 and 33.8℃ with a mean of 7.1℃ during the period, with a increasing trend (0.10℃/day, r 2 = 0.47). However, in the late February and early March (heavy snow days), it was nearly constant, not changed. The daily fluctuation of the surface soil temperature was marked and it generally was small in January, was increasing with time. The shallow soil temperatures (depth = 10 cm) showed ranges of 0.1-16.3℃ (mean = 6.3℃; right) and −0.1-16.9℃ (mean = 6.4℃; left) and they also showed the increasing trends (0.10℃/day at r 2 = 0.80 and 0.11℃ at r 2 = 0.79, respectively).
PPT Slide
Lager Image
Soil temperatures (January 12 to May 31, 2015) measured at (a) surface, (b) depth of 10 cm, 1 m right of the surface soil temperature measuring point, (c) depth of 10 cm, 1 m left of the surface soil temperature point, (d) temperature difference between surface and depth 10 cm (right), and (e) that between surface and depth 10 cm (left).
However, the amplitudes of daily fluctuations were much smaller than that of the surface soil temperature, but also in these cases, the very stable soil temperatures in late February and early March were observed when snow was heavy and its depth was large. Figs. 3d and 3e show the difference between surface and shallow subsurface soil temperatures. The differences were between −5.7 and 23.3℃ (mean = 0.8℃; Fig. 3d ), and between −5.9 and 22.4℃ (mean = 0.7℃, Fig. 3e ). The positive and negative differences were alternating with time but the positive ones are more distinctive and predominant, indicating primarily downward heat flux in the shallow soil during the period. Fig. 4 shows the bivariate plots of the surface and shallow soil temperatures. As expected, the surface and shallow soil temperatures are closely related (r 2 = 0.72 and 0.73). Considering the temperature distribution around the 1 : 1 line, it is noted that the surface soil temperature was mostly higher than the shallow subsurface soil temperatures and the temperature turnover was not significant.
PPT Slide
Lager Image
Bivariate plots of soil temperature at surface and those of depth of 10 cm (a: right (10R), b: left (10L)).
Fig. 5 shows auto-correlation and spectral density of the soil temperatures. Auto-correlation indicates inter-dependency of successive values ( Larocque et al., 1998 ; Lee and Lee, 2000 ). The auto-correlation of the surface soil temperature gradually decreases with lag times with periodic fluctuation and reach the null value (0) at a lag time of 82 hours (3.4 days) while the shallow soil temperatures show slower decreasing trends with similar fluctuations and reach the null value at 273 hours (11.4 days) and 250 hours (10.4 days), which indicate the subsurface soil temperatures have a longer memory effect than the surface soil temperature.
PPT Slide
Lager Image
Auto-correlation and spectrum of hourly soil temperatures.
The spectrums of the soil temperatures ( Fig. 5b ) reveal the periodicities of their signals. The highest power (381.2) of the surface temperature at f = 0.0416 means a daily periodicity (cyclic variation) and the second peak (69.3 at f = 0.0834) in high frequencies is indicative of the semidiurnal periodicity. The spectrums of the subsurface temperatures (depth = 10 cm) show somewhat different distributions. The powers of low frequencies (< 0.005) are greatly elevated while those of high frequencies are largely lessened. Even though the daily cyclic behavior is still distinctive (powers 139.3 and 173.7 at f = 0.0416), the semidiurnal variation almost disappeared (powers 8.4 and 10.2 at f = 0.0834). All these spectrums indicate that the surface soil temperatures quickly respond to variations in air temperature or solar radiation while the subsurface soil temperatures slowly respond to that due to heat signal attenuation during transmission through soil.
Fig. 6 shows the short-term Fourier of the soil temperatures during the monitoring period. In the first 800 hours (33.3 days), the spectrum of the surface air temperature shows discernible peaks only in low frequencies, and at around 800-1,500 hours (coincident with heavy rain and thick snow period), the spectrum is greatly strengthened in the lowest frequencies, which is well explained by very stable soil temperature for the time period (see Fig. 3a ). After that, some significant spectrums also in high frequencies (> 0.1) can be found and this is closely related with very frequent fluctuations of the surface soil temperatures for the period (after 1,500 hours) due to outer stresses such as changes in air temperature and solar radiation. The above phenomena are also found in the Fourier of the subsurface soil temperatures ( Figs. 6b and 6c ). However, due to delayed and attenuated heat transmission through soil, spectrums of lower frequencies are highly strengthened but those of high frequencies nearly completely disappeared, which is well consistent with the original soil temperatures behaviors (see Fig. 3b and 3c ). So the short-term Fourier is very useful to reveal variation characteristics of certain signals with time.
PPT Slide
Lager Image
Shot-time Fourier transforms of soil temperatures. The window function is Welch and its size is 32 (Hammer, 2012).
- 3.3. Correlations with Air Temperatures and Snow Depth
Fig. 7 shows correlations of the soil temperatures with outer air temperatures. The maximum air temperature shows very high correlation (r = 0.80) with the surface soil temperature at null lag time, which means the surface soil temperature responds quickly to the air temperature. The correlation reached the null value at a lag time of 15 days. However, the subsurface soil temperatures show somewhat delayed and lowered peak correlations (r = 0.69 and 0.69 at a lag time of 1 day), which is due to heat flux attenuation during its propagation through soil. The correlations reached the null value at lag times of somewhat delayed 17 days. However, compared with the maximum air temperature, the daily minimum air temperature shows different behaviors, showing lowered peak correlations (r = 0.72, 0.64, and 0.64, respectively) without lag times. These differences are derived from differences in time and spectral characteristics of the two air temperatures ( Figs. 7c and 7d ). The minimum air temperature has a longer memory effect like the soil temperatures and it has higher powers of low frequencies (longer periodicity) while the maximum air temperature has a shorter memory and higher powers of high frequencies ( Fig. 7d ).
PPT Slide
Lager Image
Correlations of soil temperatures with (a) maximum air temperature and (b) minimum air temperature, and (c) auto-correlations and (d) spectrums of max. and min. air temperatures.
The correlation of the daily soil temperature variations (maximum-minimum) with the snow depth was also examined. Fig. 8a shows relationships among the snow depth (SD) and daily variations in air temperature (ATV), and soil temperatures (STV). The daily air temperature variation was between 2.2 and 23.9℃ (mean = 14.2℃) but it was not much varying with time (coefficient of variation (CV) = 0.33). The daily variation in the surface soil temperature was between 0 and 30.3℃ (mean = 9.2℃) and it is generally increasing towards the warming months (CV = 0.85) except for the heavy snow days (late February to early March).
PPT Slide
Lager Image
Daily soil temperatures variations (maximum-minimum; STV) with those of air temperatures (ATV) and snow depth (SD), and their correlations. Surf: surface soil, 10R: depth = 10 cm (right), 10L: depth = 10 cm (left).
The subsurface temperature variations were very similar to those of the surface soil temperature except the amplitudes of the variations were much smaller. Also in these cases, the practically no variations of the subsurface soil temperatures are found in the snow days. As one of main influencing factor on the soil temperatures, the air temperature variation is substantially positively linked with the soil temperatures variation (r = 0.48, 0.37, and 0.38, respectively; Figs. 8b - 8d ). The snow depth, however, is inversely related with the soil temperature variations (r = −0.49, −0.46, and −0.47, respectively; Figs. 8e - 8f ). Therefore, it is considered the snow depth strongly affects the soil temperature variations due to its insulation effect by low thermal conductivity and high albedo ( Zhang, 2005 ; Zhang et al., 2008 ).
- 3.4. Estimation of Soil Temperatures using Multiple Regressions
The possibility to modeling the soil temperatures using multiple regression of the air temperature (maximum and minimum applied separately) and the snow depth, most influencing factors that are available, was examined. Fig. 9 shows the multiple regression results for the surface and subsurface soil temperatures. It is seen that the surface soil temperature is comparatively well represented by both the parameters even though they cannot explain well that for the heavy snow days ( Fig. 9a ). Interestingly, in the case of using the minimum air temperature, the surface soil temperature is better explained (r 2 = 0.64) than the maximum air temperature (r 2 = 0.56). In both cases of the subsurface soil temperatures ( Figs. 9b and 9c ), the soil temperatures in the cold days (around before March) were poorly represented by the multiple regression while those in warmer days (after March) are relatively better explained by them. Considering the low multiple r 2 (0.35 and 0.52) of those regression, however, there is a room for large improvement regarding the multiple regression by incorporating other important parameters even though the minimum air temperature shows still a better performance. It is noted that there are many influencing factors on soil temperatures ( Zhang, 2005 ; Kim et al., 2015 ). Therefore, the somewhat unsuccessful analyses in this study using two parameters in estimating the soil temperatures are not discouraging. If more data such as solar radiation and moisture content are available, the prediction can be more accurate and reliable.
PPT Slide
Lager Image
Estimated soil temperatures using multiple linear regression based on air temperature (ATmax and ATmin) and snow depth (SD).
4. Conclusion
The temporal and frequency characteristics of the soil temperatures measured from a high elevation city Boulder, Colorado were examined. The weather condition in the city is much different from low lands in that there is frequent snow and a large daily change in air temperatures. Main findings from this study are as follows: 1) The surface and shallow subsurface soil temperatures well respond to changes in the air temperature but the daily variation amplitude is larger in warmer days, 2) The short-time Fourier is useful to uncover the spectrum characteristics of the soil temperatures with elapsing time, 3) The soil temperatures variations are inversely related with the snow cover depth due to its insulation effect by low thermal conductivity and high albedo, and 4) the shallow soil temperatures are not well explained by the multiple regression using air temperature and snow depth, so further incorporation of other important parameters is needed for a better prediction. Finally additional monitoring at various points considering land use in the city is essentially required for more reliable analysis.
Acknowledgements
This research was supported by the Basic Science Research Program through the National Research Foundation of Korea (NRF), funded by the Ministry of Education (NRF-2011-0007232), also supported by the Korea Polar Research Institute (PP15010), and Pusan National University Research Grant, 2012. Mr. Woo-Hyun Jeon supported the installation of the measuring equipment, Mr. Sang-Woong Yun helped the recovery of the equipment, and the manual measurements were conducted with the first author’s son Jeong Woo. We also appreciate helpful comments by reviewers.
References
Adlam L.S. , Balks M.R. , Seybold C.A. , Campbell D.I. 2010 Temporal and spatial variation in active layer depth in the McMurdo Sound Region, Antarctica Antarc. Sci. 22 (1) 45 - 52    DOI : 10.1017/S0954102009990460
Baker V.R. 1974 Urban geology of Boulder, Colorado: A progress report Environ. Geol. 1 75 - 88
2015 Geology of the Boulder area http://bcn.boulder.co.us/basin/watershed/
Beniston M. 2006 Mountain weather and climate: A general overview and a focus on climatic change in the Alps Hydrobiol. 562 3 - 16    DOI : 10.1007/s10750-005-1802-0
Beniston M. , Diaz H.F. , Bradley R.S. 1997 Climate change at high elevation sites: An overview Clim. Change 36 233 - 251    DOI : 10.1023/A:1005380714349
Bilodeau S.W. , Buskirk D.V. , Bilodeau W. 1987 Geology of Boulder, Colorado, United States of America Bull. Assoc. Engrg. Geol. 24 (3) 289 - 332
Box G.E.P. , Jenkins G.M. , Reinsel G.C. 1994 Time Series Analysis: Forecasting and Control Prentice-Hall Englewood Cliffs, NJ 598 -
Chapin III F.S. , Zavaleta E.S. , Eviner V.T. , Naylor R.L. , Vitousek P.M. , Reynolds H.L. , Hooper D.U. , Lavorel S. , Sala O.E. , Hobbie S.E. , Mack M.C. , Diaz S. 2000 Consequences of changing biodiversity Nature 405 234 - 242    DOI : 10.1038/35012241
2014 2014 Community Profile https://www-static.bouldercolorado.gov/
Davidoff B. , Selim H.M. 1988 Correlation between spatially variable soil moisture content and soil temperature Soil Sci. 145 (1) 1 - 8    DOI : 10.1097/00010694-198801000-00001
Duffy C.J. , Gelhar L.W. 1986 A frequency domain analysis of groundwater quality fluctuations: interpretation of field data Water Resour. Res. 22 (7) 1115 - 1128    DOI : 10.1029/WR022i007p01115
Fennenman N.M. 1905 Geology of the Boulder district, Colorado Government Printing Office Washington DC USGS Bulletin 265 101 -
Hammer Ø. 2012 PAST: PAleontological STatistics Version 2.17 Natural History Museum, University of Oslo 229 -
Kim J.H. , Lee J.H. , Cheong T.J. , Kim R.H. , Koh D.C. , Ryu J.S. , Chang H.W. 2005 Use of time series analysis for the identification of tidal effect on groundwater in the coastal area of Kimje, Korea J. Hydrol. 300 188 - 198    DOI : 10.1016/j.jhydrol.2004.06.004
Kim H. , Lee J.Y. , Jeon W.H. , Lee K.K. , Shrestha S. , Pandey V.S. 2015 Groundwater Environment in Asian Cities: Concepts, Methods and Case Studies Elsevier The Netherlands Evaluation of groundwater environment of Seoul, Korea (in press)
Larocque M. , Mangin A. , Razack M. , Banton O. 1998 Contribution of correlation and spectral analyses to the regional study of a karst aquifer (Charente, France) J. Hydrol. 205 217 - 231    DOI : 10.1016/S0022-1694(97)00155-8
Lee J.Y. , Lee K.K. 2000 Use of hydrologic times series data for identification of recharge mechanism in a fractured bedrock aquifer system J. Hydrol. 229 (3-4) 190 - 201    DOI : 10.1016/S0022-1694(00)00158-X
Lim H.G. , Choi H.M. , Lee J.Y. 2011 Comparison of hydrological time series analysis results before and after detrending J. Soil Ground. Env. (in Korean with English abstract) 16 (2) 30 - 40    DOI : 10.7857/JSGE.2011.16.2.030
Morss R.E. , Demuth J.L. , Bostrom A. , Lazo J.K. , Lazrus H. 2015 Flash flood risks and warning decisions: A mental models study of forecasts, public officials, and media broadcasters in Boulder, Colorado Risk Anal.    DOI : 10.1111/risa.12403
Nikolaidis N. , Ragnarsdottir K.V. , Ragnarsdottir K.V. , Banwart S.A. 2015 Soil: The Life Supporting Skin of Earth University of Sheffield Sheffield, UK Soil and Earth’s critical zone 9 - 15
2015 Boulder, Colorado climate and weather data http://www.esrl.noaa.gov/psd/boulder
Padilla A. , Pulido-Bosch A. 1995 Study of hydrographs of karstic aquifers by means of correlation and cross-spectral analysis J. Hydrol. 168 73 - 89    DOI : 10.1016/0022-1694(94)02648-U
Rigg L. 1993 Study of soil properties along a hillslope in Parson’s Parcel Boulder, CO University of Colorado Geography 4411 40 -
Rustad L.E. , Huntington T.G. , Boone R.D. 2000 Controls on soil respiration: Implications for climate change Biogeochem. 48 (1) 1 - 6    DOI : 10.1023/A:1006255431298
Teubner I.E. , Haimberger L. , Hantel M. 2015 Estimating snow cover duration from ground temperature J. Appl. Meteor. Climatol. 54 959 - 965    DOI : 10.1175/JAMC-D-15-0006.1
Wells J.D. 1967 Geology of the Eldorado Springs Quadrangle Boulder and Jefferson Counties, Colorado Government Printing Office Washington USGS Bulletin 1221-D 91 -
Wundram D. , Pape R. , Loffler J. 2010 Alpine soil temperature variability at multiple scales Arc. Antarc. Alp. Res. 42 (1) 117 - 128    DOI : 10.1657/1938-4246-42.1.117
Yoshioka M. , Takakura S. , Ishizawa T. , Sakai N. 2015 Temporal changes of soil temperature with soil water content in an embankment slope during controlled artificial rainfall experiments J. Appl. Geophy. 114 134 - 145    DOI : 10.1016/j.jappgeo.2015.01.009
Zhang T. 2005 Influence of the seasonal snow cover on the ground thermal regime: An overview Rev. Geophys. 43 RG4002 -
Zhang T. , Osterkamp T.E. , Stamnes K. 1997 Effects of climate on the active layer and permafrost on the North Slope of Alaska, USA Permafrost Periglacial Process. 8 45 - 67    DOI : 10.1002/(SICI)1099-1530(199701)8:1<45::AID-PPP240>3.0.CO;2-K
Zhang Y. , Wang S. , Barr A.G. , Black T.A. 2008 Impact of snow cover on soil temperature and its simulation in a boreal aspen forest Cold Regions Sci. Technol. 52 355 - 370    DOI : 10.1016/j.coldregions.2007.07.001