The Cryosphere Discussions www.the-cryosphere-discuss.net 1994-0432 1994-0440 1 1 2007 10.5194/tcd-1-213-2007 http://www.the-cryosphere-discuss.net/1/213/2007/ http://www.the-cryosphere-discuss.net/1/213/2007/tcd-1-213-2007.html http://www.the-cryosphere-discuss.net/1/213/2007/tcd-1-213-2007.pdf 213 269 2007-08-09 Estimation of thermal properties of saturated soils using in-situ temperature measurements D. J. Nicolsky ftdjn@uaf.edu V. E. Romanovsky G. S. Tipenko Geophysical Institute, University of Alaska Fairbanks, PO Box 757320, Fairbanks, AK 99775, USA Institute of Environmental Geoscience Russian Academy of Sciences, 13-2 Ulansky pereulok, PO Box 145, Moscow, Russia We describe an approach to find an initial approximation to the thermal properties of soil horizons. This technique approximates thermal conductivity, porosity, unfrozen water content curve in horizons where no direct temperature measurements are available. To determine physical properties of ground material, optimization-based inverse modeling techniques fitting the simulated and measured temperatures are commonly employed. Two major ingredients of these techniques is an algorithm to compute the soil temperature dynamics and a procedure to find an initial approximation to the ground properties. In this article we show how to determine the initial approximation to the physical properties and present a new finite element discretization of the heat equation with phase change to calculate the temperature dynamics in soil. We successfully applied the proposed algorithm to recover the soil properties for Happy Valley site in Alaska using one-year temperature dynamics. The determined initial approximation was utilized to simulate the temperature dynamics over several consecutive years; the difference between simulated and measured temperatures lies within uncertainties of measurements. ACIA: Impacts of a Warming Arctic: Arctic Climate Impact Assessment, Cambridge University Press, 139pp, 2004. Alifanov, O.: Inverse Heat Transfer Problems, Springer, Berlin, 1994. Alifanov, O., Artyukhin, E., and Rumyantsev, S.: Extreme Methods for Solving Ill-Posed Problems with Application to Inverse Heat Transfer Problems, Begell House, New York, 1996. Andersland, O. and Anderson, D.: Geotechnical Engineering for Cold Regions, McGraw-Hill, 1978. Anderson, D. and Morgenstern, N.: Physics, chemistry and mechanics of frozen ground: A review, in: Proceedings Second International Conference on Permafrost, 257–288, Yakutsk, USSR, 1973. Avriel, M.: Nonlinear Programming: Analysis and Methods, Dover Publications, 2003. Bazaraa, M., Sherali, H., and Shetty, C M.: Nonlinear Programming: Theory and Algorithms, John Wiley & Sons, 2nd edn., 1993. Boike, J. and Roth, K.: Time Domain Re?ectometry as a Field Method for Measuring Water Content and Soil Water Electrical Conductivity at a Continuous Permafrost Site, Permafrost and Periglacial Processes, 8, 359–370, 1997. Brown, J., Ferrians, O., Heginbottom, J J., and Melnikov, E.: Circum-Arctic map of permafrost and ground-ice conditions, U.S. Geological Survey Circum-Pacific Map CP-45, 1:10 000 000, reston, Virginia, 1997. Carslaw, H. and Jaeger, J.: Conduction of heat in solids, Oxford Univ. Press, London, 1959. Carson, J.: Analysis of soil and air temperatures by Fourier techniques, Journal Geophysical Research, 68, 2217–2232, 1963. Comini, G., Giudice, S D., Lewis, R., and Zienkiewicz, O.: Finite element solution of non-linear heat conduction problems with special reference to phase change, Int. J. Numer. Methods Engng, 8, 613–624, 1974. Dalhuijsen, A. and Segal, A.: Comparison of finite element techniques for solidification problems, Int.j .numer.methods eng., 23, 1807–1829, 1986. Dennis, J. and Schnabel, R.: Numerical Methods for Unconstrained Optimization and Nonlinear Equations, Prentice-Hall, 1987. de~Vries, D.: W. R. van Wijk (editor) Physics of the plant environment, chap. Thermal Properties of Soils, North-Holland, Amsterdam, 1963. Fletcher, R.: Practical methods of optimization, John Wiley & Sons, 2000. Goodrich, W.: The influence of snow cover on the ground thermal regime, Candian Geotechnical Journal, 19, 1982. Gupta, S.: The Classical Stefan Problem, Elsevier, Amsterdam, 2003. Hinkel, K M.: Estimating seasonal values of thermal diffusivity in thawed and frozen soils using temperature time series, Cold Reg. Sci. Technol., 26, 1–15, 1997. Hinzman, L., Kane, D., Gleck, R., and Everett, K.: Hydrological and thermal properties of the active layer in the Alsakan Arctic, Cold Region Science and Technology, 19, 95–110, 1991. Hobbs, P.: Ice physics, Claredon Press, Oxford, 1974. Hurley, S. and Wiltshire, R J.: Computing thermal diffusivity from soil temperature measurements, Comput. Geosci., 19, 475–477, 1993. Jaeger, J C. and Sass, J H.: A line source method for measuring the thermal conductivity and diffusivity of cylindrical specimens of rock and other poor conductors, J. Appl. Phys., 15, 1187–1194, 1964. Javierre, E., Vuik, C., Vermolen, F., and van~der Zwaag, S.: A comparison of numerical models for one-dimensional Stefan problems, Journal of Computational and Applied Mathematics, 192, 445–459, 2006. Kane, D., Hinzman, L., and Zarling, J.: Thermal response of the active layer in a permafrost environment to climatic warming, Cold Region Science Technology, 19, 111––122, 1991. Kane, D., Hinkel, K., Goering, D., Hinzman, L., and Outcalt, S.: Non-conductive Heat Transfer Associated with Frozen Soils, Global and Planetary Change, 29, 275–292, 2001. Lagarias, J., Reeds, J., Wright, M., and Wright, P.: Convergence Properties of the Nelder-Mead Simplex Method in Low Dimension, SIAM Journal of Optimization, 9, 112–147, 1998. Lemmon, E.: R.W. Lewis and K. Morgan (ed.), Numerical Methods in Thermal Problems, chap. Phase change technique for finite element conduction code, pp. 149–158, Pineridge Press, Swansea, U.K., 1979. Ling, F. and Zhang, T.: Impact of the timing and duration of seasonal snow cover on the active layer and permafrost in the Alaskan Arctic, Permafrost Periglacial Processes, 14, 141––150, 2003. Lovell, C.: Temperature effects on phase composition and strength of partially frozen soil, Highway Res. Board Bull., 168, 74–95, 1957. McGaw, R W., Outcalt, S I., and Ng, E.: Thermal properties of wet tundra soils at Barrow, Alaska, Proceedings of 3rd International Conference on Permafrost, 1, 47–53, nat. Res. Counc. of Can., Ottawa,Canada, 1978. Molders, N. and Romanovsky, V.: Long-term evaluation of \mboxHTSVS frozen ground/permafrost component using observations at Barrow, Alaska., J. Geophys. Res., 111, doi:10.1029/2005JD005957, 2006. Morgan, K., Lewis, R., and Zienkiewicz, O.: An improved algorithm for heat conduction problems with phase change, Int. j .numer.methods eng., 12, 1191–1195, 1978. Mundim, M. and Fortes, M.: edited by: Numerical Methods in Thermal Problems, edited by: Lewis, R. W. and Morgan, K., vol 6, chap. Evaluation of finite element method utilised in the solution of solid-liquid phase change problems, pp. 90–100, Pineridge Press, Swansea, U.K., 1979. Muzylev, N V.: The uniqueness of a solution to an inverse nonlinear heat conduction problem, Zhurnal Vychislitel'noi Matematiki i Matematicheskoi Fiziki, 25, 1346–1352, in Russian, 1985. Nelson, F. and Outcalt, S.: Anthropogenic geomorphology in northern Alaska, Physical Geography, 3, 17–48, 1987. Oleson, K W., Dai, Y., Bonan, G., Bosilovich, M., Dickinson, R., Dirmeyer, P., Hoffman, F., Houser, P., Levis, S., Niu, G.-Y., Thornton, P., Vertenstein, M., Yang, Z.-L., and Zeng, X.: Technical description of the Community Land Model (\mboxCLM), \mboxNCAR tech. note \mboxNCAR/TN-461+STR, NCAR, 2004. Osterkamp, T. and Romanovsky, V.: Characteristics of changing permafrost termperatures in the Alaskan Arctic, USA, Arctic and Alpine research, 28, 267–273, 1996. Osterkamp, T. and Romanovsky, V.: Freezing of the active layer on the Coastal Plain of the Alaskan Arctic, Permafrost and Periglacial Processes, 8, 23–44, 1997. Permyakov, P.: Methods of Determining the Characteristics of Dispersed Media at a Phase Transition, Russian Physics Journal, 2004. Pham, Q T.: Comparison of general-purpose finite element methods for the Stefan problem, Numerical Heat Transfer Part B - Fundamentals, 27, 417–435, 1995. Pinder, G. and Gray, W.: Finite Element Simulation in Surface and Subsurface Hydrology, Academic Press, New York, 1977. Rank, E., Katz, C., and Werner, H.: On the importance of the discrete maximum principle in transient analysis using finite element methods, Int. j .numer.methods eng., 19, 1771–1782, 1983. Romanovsky, V. and Osterkamp, T.: Interannual variations of the thermal regime of the active layer and near-surface permafrost in Northern Alaska, Permafrost and Periglacial Processes, 6, 313–335, 1995. Romanovsky, V. and Osterkamp, T.: Thawing of the Active Layer on the Coastal Plain of the Alaskan Arctic, Permafrost and Periglacial processes, 8, 1–22, 1997. Samarskii, A. and Vabishchevich, P.: Computational Heat Transfer, Mathematical Modelling, vol. 1,2, Wiley, first edn., 1996. Sass, J., Lachenbruch, A., and Munroe, R.: Thermal conductivity of rocks from measurements on fragments and its application to heat-flow determinations, J. Geophys. Res., 76, 3391–3401, 1971. Sazonova, T., Romanovsky, V., Walsh, J., and Sergueev, D.: Permafrost dynamics in the 20th and 21st centuries along the East Siberian transect, J. Geophys. Res., 109, D01108, doi:10.1029/2003JD003680, 2004. Schmugge, T., Jackson, T., and McKim, H.: Survey of Methods for Soil Moisture Determination, Water Resour. Res., 16, 961–979, 1980. Smith, M. and Tice, A.: Measurement of the unfrozen water content of soils—comparison of NMR and TDR methods, CRREL Report 88–18, US Army Cold Regions Research and Engineering Lab (CRREL), 1988. Stafford, J.: REMOTE, NON-CONTACT AND IN-SITU MEASUREMENT OF SOIL MOISTURE CONTENT: A REVIEW, The Journal of Agricultural Engineering Research, pp. 151–172, 1988. Thacker, W.: The role of the Hessian matrix in fitting models to measurements, Journal of Geophysical Research, 94, 6177–6196, 1989. Thacker, W. and Long, R.: Fitting dynamics to data, J. Geophys. Res., 93, 1227––1240, 1988. Tice, A., Oliphant, J., Nakano, Y., and Jenkins, T.: Relationship between the ice and unfrozen water phases in frozen soil as determined by pulsed nuclear magnetic resonance and physical desorption data, CRREL Report 82–15, US Army Cold Regions Research and Engineering Lab (CRREL), 1982. Tikhonov, A. and Samarskii, A.: Equations of mathematical physics, Pergamon, 1963. Tikhonov, L. and Leonov, A.: Nonlinear Ill-Posed Problems, vol 14 of \em Applied Mathematics and Mathematical Computation Series\/, Chapman and Hall, London, 1996. Topp, G C., Davis, J L., and Annan, A P.: Electromagnetic determination of soil water content: measurements in coaxial transmission lines, Water Resour. Res., 16, 574–582, 1980. Ulaby, F., Moore, R., and Fung, A.: Microwave Remote Sensing, Addison-Wesley, Reading, MA, 1982. Voller, V. and Swaminathan, C.: Fixed grid techniques for phase change problems: a review, Int.j .numer.methods eng., 30, 875–898, 1990. Watanabe, K. and Mizoguchi, M.: Amount of unfrozen water in frozen porous media saturated with solution, Cold Regions Science and Technology, 34, 103––110, 2002. Williams, P.: Properties and behaviour of freezing soils, Tech. Rep 72, Norwegian Geotechnical Institute, 1967. Yershov, E.: General Geocryology, Cambridge University Press, New York, 1998. Yoshikawa, K., Overduin1, P., and Harden, J.: Moisture Content Measurements of Moss (Sphagnum spp.) using Commercial Sensors, Permafrost and Periglacial Processes, 15, 309-318, 2004. Zhang, T. and Osterkamp, T E.: Considerations in determining the thermal diffusivity from temperature time series using finite difference methods, Cold Reg. Sci. Technol., 23, 333–341, 1995. Zhuang, Q., Romanovsky, V., and McGuire, A.: Incorporation of a permafrost model into a large-scale ecosystem model: Evaluation of temporal and spatial scaling issues in simulating soil thermal dynamics, Journal Geophysical Research, 106, 33 649––33 670, 2001. Zienkiewicz, O. and Taylor, R.: The finite element method, vol. 1,2, McGraw-Hill, London, fourth edn., 1991.