The Cryosphere www.the-cryosphere.net 1994-0416 1994-0424 1 1 2007 10.5194/tc-1-41-2007 http://www.the-cryosphere.net/1/41/2007/ http://www.the-cryosphere.net/1/41/2007/tc-1-41-2007.html http://www.the-cryosphere.net/1/41/2007/tc-1-41-2007.pdf 41 58 2007-11-22 Using in-situ temperature measurements to estimate saturated soil thermal properties by solving a sequence of optimization problems 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 curves in horizons where no direct temperature measurements are available. To determine physical properties of ground material, optimization-based inverse techniques are employed to fit the simulated temperatures to the measured ones. Two major ingredients of these techniques are 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 apply the proposed algorithm to recover the soil properties for the Happy Valley site in Alaska using one-year temperature dynamics. The determined initial approximation is 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, 139 pp., 2004. Alifanov, O.: Inverse Heat Transfer Problems, Springer, Berlin, 348 pp., 1995. 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, 306 pp., 1996. Andersland, O. and Anderson, D.: Geotechnical Engineering for Cold Regions, McGraw-Hill, 566 pp., 1978. Anderson, D. and Morgenstern, N.: Physics, chemistry and mechanics of frozen ground: a review, in: Proceedings of the 2nd International Conference on Permafrost, Yakutsk, USSR, 257–288, 1973. Avriel, M.: Nonlinear Programming: Analysis and Methods, Dover Publications, 554 pp., 2003. Bazaraa, M., Sherali, H., and Shetty, C M.: Nonlinear Programming: Theory and Algorithms, 2nd edn., John Wiley & Sons, 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 Periglac., 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 University Press, London, 520 pp., 1959. Carson, J.: Analysis of soil and air temperatures by Fourier techniques, J. Geophys. Res., 68, 2217–2232, 1963. Ciarlet, P.: The Finite Element Method for Elliptic Problems, North-Holland, Amsterdam, 530 pp., 2002. 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. Meth. Eng., 8, 613–624, 1974. Dalhuijsen, A. and Segal, A.: Comparison of finite element techniques for solidification problems, Int. J. Numer. Meth. Eng., 23, 1807–1829, 1986. de~Vries, D.: Physics of the Plant Environment, chap. Thermal properties of soils, edited by: van Wijk, W. R., Wiley, New York, 210–235, 1963. Dennis, J. and Schnabel, R.: Numerical Methods for Unconstrained Optimization and Nonlinear Equations, Prentice-Hall, 394 pp., 1987. Fletcher, R.: Practical Methods of Optimization, John Wiley & Sons, 450 pp., 2000. Galushkin, Y.: Numerical simulation of permafrost evolution as a part of sedimentary basin modeling: permafrost in the Pliocene-Holocene climate history of the Urengoy field in the West Siberian basin, Can. J. Earth Sci., 34, 935–948, 1997. Goldberg, D.: Genetic Algorithms in Search, Optimization, and Machine Learning, Addison Wesley, 432 pp., 1989. Goodrich, W.: The influence of snow cover on the ground thermal regime, Can. Geotech. J., 19, 421–432, 1982. Gupta, S.: The Classical Stefan Problem, Elsevier, Amsterdam, 404 pp., 2003. Hinkel, K.: 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 Alaskan Arctic, Cold Reg. Sci. Technol., 19, 95–110, 1991. Hobbs, P.: Ice Physics, Claredon Press, Oxford, 856 pp., 1974. Hurley, S. and Wiltshire, R J.: Computing thermal diffusivity from soil temperature measurements, Computers and Geosciences, 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, J. Comput. Appl. Math., 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 Reg. Sci. Technol., 19, 111–122, 1991. Kane, D., Hinkel, K., Goering, D., Hinzman, L., and Outcalt, S.: Non-conductive heat transfer associated with frozen soils, Global Planet. Change, 29, 275–292, 2001. Kolmogorov, A. and Fomin, S.: Introductory Real Analysis, Prentice-Hall, New York, 403 pp., 1975. Lagarias, J., Reeds, J., Wright, M., and Wright, P.: Convergence properties of the Nelder-Mead simplex method in low dimension, SIAM J. Optimiz., 9, 112–147, 1998. Lemmon, E.: Numerical Methods in Thermal Problems, chap. Phase change technique for finite element conduction code, edited by: Lewis, R. W. and Morgan, K., Pineridge Press, Swansea, UK, 149–158, 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 Periglac., 14, 141–150, 2003. Lovell, C.: Temperature effects on phase composition and strength of partially frozen soil, Highway Research Board Bulletin, 168, 74–95, 1957. Lunardini, V.: Freezing of soil with an unfrozen water content and variable thermal properties, CRREL Report 88-2, US Army Cold Regions Research and Engineering Lab, 30 pp., 1987. \mboxMölders, N. and Romanovsky, V.: Long-term evaluation of \mboxHTSVS frozen ground/permafrost component using observations at Barrow, Alaska., J. Geophys. Res., 111, D04105, doi:10.1029/2005JD005957, 2006. \mboxThermal Logic: Thermal conductivity sensor user's manual, PO Box 781, Pullman, WA, 99163, 3 edn., 25 pp., 2001. McGaw, R W., Outcalt, S I., and Ng, E.: Thermal properties of wet tundra soils at Barrow, Alaska, in: Proceedings of 3rd International Conference on Permafrost, vol 1, National Research Council of Canada, Ottawa, Canada, 47–53, 1978. Morgan, K., Lewis, R., and Zienkiewicz, O.: An improved algorithm for heat conduction problems with phase change, Int. J. Numer. Meth. Eng., 12, 1191–1195, 1978. Mundim, M. and Fortes, M.: Numerical Methods in Thermal Problems, vol 6, chap. Evaluation of finite element method utilized in the solution of solid-liquid phase change problems, edited by: Lewis, R. W. and Morgan, K., Pineridge Press, Swansea, UK, 90–100, 1979. Nelson, F. and Outcalt, S.: Anthropogenic geomorphology in northern Alaska, Polar Geography, 3, 17–48, 1987. Oleson, K., 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, 173 pp., 2004. Osterkamp, T. and Romanovsky, V.: Characteristics of changing permafrost temperatures in the Alaskan Arctic, USA, Arctic Alpine Res., 28, 267–273, 1996. Osterkamp, T. and Romanovsky, V.: Freezing of the active layer on the coastal plain of the Alaskan Arctic, Permafrost Periglac., 8, 23–44, 1997. Permyakov, P.: Methods of determining the characteristics of dispersed media at a phase transition, Russian Physics Journal, 47(3), 240–246, 2004. Pham, Q.: Comparison of general-purpose finite element methods for the Stefan problem, Numer. Heat Tr. B-Fund., 27, 417–435, 1995. Pinder, G. and Gray, W.: Finite Element Simulation in Surface and Subsurface Hydrology, Academic Press, New York, 295 pp., 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. Meth. Eng., 19, 1771–1782, 1983. Robert, C. and Casella, G.: Monte Carlo Statistical Methods, Springer-Verlag, 645 pp., 2004. Romanovsky, V. and Osterkamp, T.: Interannual variations of the thermal regime of the active layer and near-surface permafrost in Northern Alaska, Permafrost Periglac., 6, 313–335, 1995. Romanovsky, V. and Osterkamp, T.: Thawing of the active layer on the coastal plain of the Alaskan Arctic, Permafrost Periglac., 8, 1–22, 1997. Romanovsky, V. and Osterkamp, T.: Effects of unfrozen water on heat and mass transport processes in the active layer and permafrost, Permafrost Periglac., 11, 219–239, 2000. Samarskii, A. and Vabishchevich, P.: Computational Heat Transfer, Mathematical Modeling, vol 1, Wiley, 418 pp., 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 \mboxNMR and \mboxTDR methods, CRREL Report 88-18, US Army Cold Regions Research and Engineering Lab, 16 pp., 1988. Stafford, J.: Remote, non-contact and in-situ measurement of soil moisture content: a review, J. Agr. Eng. Res., 14, 151–172, 1988. Thacker, W.: The role of the Hessian matrix in fitting models to measurements, J. Geophys. Res., 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, 12 pp., 1982. Tikhonov, A. and Samarskii, A.: Equations of Mathematical Physics, Pergamon, 776 pp., 1963. Tikhonov, A., Leonov, A., and Yagola, A. G.: Nonlinear Ill-Posed Problems, Chapman and Hall, London, 392 pp., 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, Vol. II, 1064 pp., 1982. Voller, V. and Swaminathan, C.: Fixed grid techniques for phase change problems: a review, Int. J. Numer. Meth. Eng., 30, 875–898, 1990. Watanabe, K. and Mizoguchi, M.: Amount of unfrozen water in frozen porous media saturated with solution, Cold Reg. Sci. Technol., 34, 103–110, 2002. Williams, P.: Properties and behaviour of freezing soils, Tech. Rep 72, Norwegian Geotechnical Institute, 120 pp., 1967. Yershov, E.: General Geocryology, Cambridge University Press, Cambridge, 604 pp., 1998. Yoshikawa, K., Overduin, P., and Harden, J.: Moisture content measurements of moss (Sphagnum spp.) using commercial sensors, Permafrost Periglac., 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, J. Geophys. Res., 106, 33 649–33 670, 2001. Zienkiewicz, O. and Taylor, R.: The Finite Element Method, vol 1, McGraw-Hill, London, 470 pp., 1991.