Modeled Daily Thaw Depth and Freezing Depth


This data set contains modeled daily thaw depth and freezing depth for the Arctic terrestrial drainage basin. Thaw and freezing depths were calculated over the study area with a resolution of 25 km x 25 km in the NSIDC EASE-Grid projection. Data extend from September 1998 through December 2000. A finite-difference model for one-dimensional heat conduction with phase change calculates daily thaw and freezing depths, based on data of snow water equivalent, soil temperature, soil dry-bulk density, and soil water content. Data are in space-delimited ASCII text format and are available via FTP.

Citing These Data:

Oelke, C., R. Armstrong, M. Serreze, and T. Zhang. 2003. Modeled Daily Thaw Depth and Freezing Depth. Boulder, CO: National Snow and Ice Data Center. Digital Media.

Overview Table

Category Description
Data format Data are in space-delimited ASCII text format, compressed with gzip.
Spatial coverage and resolution This data set covers the Arctic terrestrial drainage basin. Thaw and freezing depths were calculated over the study area with a resolution of 25 km x 25 km.
Temporal coverage and resolution Daily data are available from 01 September 1998 through 31 December 2000, grouped in monthly files.
Grid type and size Data are gridded to the Northern Hemisphere EASE-Grid projection.
File naming convention
File size Data files compressed with gzip: 627 KB - 1.8 MB
Data files uncompressed: 10.9 MB
ease721_drainage_lakglac.txt: 1.3 MB
Parameter(s) Data values include daily thaw depth and freezing depth in cm.
Procedures for obtaining data Data are available via FTP.

Table of Contents

1. Contacts and Acknowledgments
2. Detailed Data Description
3. Data Access and Tools
4. Data Acquisition and Processing
5. References and Related Publications
6. Document Information

1. Contacts and Acknowledgments:

Investigator(s) Name and Title:

Christoph Oelke
Principal investigator
National Snow and Ice Data Center

Richard Armstrong, Mark Serreze, and Tingjun Zhang
National Snow and Ice Data Center

Technical Contact:

NSIDC User Services
National Snow and Ice Data Center
University of Colorado
Boulder, CO 80309-0449  USA
phone: +1 303.492.6199
fax: +1 303.492.2468
form: Contact NSIDC User Services


The National Aeronautics and Space Administration (NASA) Interdisciplinary Study (IDS), grant number NAG5-9568, provided funding for this research.

2. Detailed Data Description:


Data are in space-delimited ASCII text format, compressed with gzip.

File and Directory Structure:

Data are separated into 28 thaw depth files and 28 freezing depth files; each monthly file contains daily data for that month. The file "ease721_drainage_lakglac.txt" provides the latitude and longitude of each EASE-Grid pixel within the Arctic drainage basin (39,926 total pixels) and a permafrost code that indicates if the cell represents continuous, discontinuous, or isolated permafrost; seasonally-frozen ground; or glacier/lake.

File Naming Convention:

The file-naming convention is "" and "," where:

ALD = thaw depth
FGD = freezing depth
yyyy = four-digit year
mm = two-digit month
gz = compressed with gzip

File Size:

Data files compressed with gzip: 627 KB - 1.8 MB
Data files uncompressed: 10.9 MB
ease721_drainage_lakglac.txt: 1.3 MB

Spatial Coverage:

This data set covers the entire Arctic terrestrial drainage basin. This area includes land with ice or water that empty into the Arctic Ocean, Hudson Bay, James Bay, Hudson Strait, and the Bering Sea. The Arctic drainage basin spans overs regions of seasonally-frozen ground in the south, to permafrost regions in the north. Continuous permafrost underlies nearly half of the Arctic drainage basin (Brown et al. 1998, Zhang et al. 2000).

Spatial Coverage Map:

Extent of thaw depth measurements
Extent of thaw depth. Black indicates areas of seasonally-frozen ground.
Extent of freezing depth measurements
Extent of freezing depth. Black indicates areas of continuous permafrost.

Spatial Resolution:

Thaw depth and freezing depth are calculated with a resolution of 25 km x 25 km.

Grid Description:

Input data are interpolated to the NSIDC Northern Hemisphere EASE-Grid. Modeled thaw depth data and freezing depth data are gridded to a spatial subset of the Northern Hemisphere EASE-Grid projection that represents the Arctic drainage basin (39,926 total pixels). The file "ease721_drainage_lakglac.txt" provides the corresponding EASE-Grid latitude and longitude of each thaw depth or freezing depth value within the Arctic drainage basin subset. Data are output to ASCII text format.

Temporal Coverage:

Daily data are available from 01 September 1998 through 31 December 2000.

Parameter or Variable:

Parameter Description:

This data set provides daily thaw depth and freezing depth. The thaw layer is the layer of ground subject to annual thawing and freezing in areas underlain by permafrost. Freezing depth refers to soil or rock in which part or all of the pore water has turned into ice in areas not underlain by permafrost. See the Frozen Ground Glossary for more details.

Unit of Measurement:

Depth values are in cm.

Sample Data Record:

Each data file consists of a one-line header followed by columns for EASE-Grid cell numbers, latitude, longitude, and depth (m) for each day of the month. Missing data, or when days 29, 30 and 31 do not exist, are indicated by values of "-9999.0." Following are the first four lines in a sample thaw-depth file called "ALD.1998.09.txt":

1998 09 Active Layer Depth (ALD) in [m/100], C. Oelke, NSIDC
"CellID" "Lat" "Long" "01" "02" "03" "04" "05" "06" "07" "08" "09" "10" "11" "12" "13" "14" "15" "16" "17" "18" "19" "20" "21" "22" "23" "24" "25" "26" "27" "28" "29" "30" "31"
1 56.5398 201.6560 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0
2 56.6259 201.2912 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0 -9999.0
3 56.7107 200.9245 79.8 80.4 81.4 82.3 83.2 84.2 85.1 85.9 86.8 87.7 88.5 89.3 90.2 91.0 91.8 92.6 93.4 94.2 95.0 95.8 96.5 97.3 98.0 98.7 99.4 100.1 100.8 101.5 102.1 102.8 -9999.0

The file "ease721_drainage_lakglac.txt" provides latitude and longitude of each EASE-Grid cell. Columns represent an index number, EASE-Grid i,j coordinates, latitude, longitude, and a permafrost code, where:

1 = continuous permafrost
2 = discontinuous permafrost
3 = isolated permafrost
4 = seasonally-frozen ground
5 = glacier or lake

Grid cells with a permafrost code of "5" are inside the Arctic drainage basin, but they should be masked because the model does not account for heat transfer through glaciers, ice caps, or lakes. Grid cells for seasonally-frozen ground ("4") are set to the minimum value in thaw depth files.

Following are the first five lines of "ease_721_drainage_lakglac.txt":

1 225 307 56.5398 201.6560 4
2 225 308 56.6259 201.2912 4
3 225 309 56.7107 200.9245 3
4 226 306 56.6703 202.1664 4
5 226 307 56.7584 201.8014 4

Error Sources:

Modeled thaw depths for sporadic permafrost areas, primarily in the southern parts of the Arctic drainage basin, can be spuriously high. In these regions, permafrost is very isolated and occurs at sub-grid scales. The model input data (surface air temperature, snow cover, and soil bulk density) are likely not representative of the true forcing conditions for these small areas, and can produce an unrealistic increase of thaw depth with time. Frozen-ground depths of small non-permafrost areas within discontinuous permafrost are likely too high, as forcing parameters are more representative of colder permafrost climate at these grid cells (Oelke et al. 2003).

Quality Assessment:

Oelke et al. (2003) compared modeled thaw depth measurements for the years 1999 and 2000 with actual measurements from 60 Circumpolar Active Layer Measurement (CALM) sampling locations. The modeled values agree reasonably well with the measured values within their standard deviation. A Root Mean Square (RMS) error of 33.5 cm between measured depth from the 2000 CALM samples and modeled thaw depth results primarily from the difference in grid resolution between 10-m CALM grid values (100 m2 x 100 m2 areas) and the 25-km grid interval used in the model.

3. Data Access and Tools:

Data Access:

Data are available via FTP.

Related Data Collections:

Please visit the Frozen Ground Data Center site at NSIDC for a complete listing of data sets.

4. Data Acquisition and Processing:

Data Source:

Oelke et al. (2003) obtained model input data from the following sources:

National Center for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) reanalysis data (Kalnay et al. 1996) generate air surface temperatures at the lowest signal level (0.995 σ) with modifications. The reanalysis data have a horizontal resolution of 2.5° x 2.5°. Oelke et al. adjusted the topography, based on NCEP/NCAR tropospheric lapse rates, the 0.995 σ level temperature, and topography from a digital elevation model (DEM). Elevation is a first-order determinant of the spatial variation in surface air temperature. The topography adjustments effectively improve the resolution of the NCEP/NCAR data set, making it more compatible with the snow cover data also used for model forcing.

Snow Cover
Oelke et al. (2003) calculated snow water equivalent from the DMSP SSM/I-SSMIS Daily Polar Gridded Brightness Temperatures data set, available from NSIDC. The Chang (1987) algorithm, originally used for calculating snow water equivalent from Scanning Multichannel Microwave Imager (SMMR) passive-microwave data from 1978 to 1987, was modified for SSM/I data.

Oelke et al. derived snow height from snow water equivalent by dividing by a climatological snow density at various locations and seasons. A 45-year time series (1955-1999) of Canadian snow data from the Meteorological Service of Canada (MSC 2000) defined the climatological seasonal cycles of snow density for tundra, taiga, prairie, alpine and maritime regions. Sturm et al. (1995) defined these snow classes, based on climatological values of temperature, precipitation, and wind speed. Over 90% of snow in the Arctic drainage basin is from tundra and taiga regions. More information about the snow classes is available from the Global Seasonal Snow Classification System data set, available from the Arctic System Science Data Coordination Center (ARCSS-ADCC) at NSIDC.

Passive-microwave sensors cannot detect very thin snow cover because this type of snow does not provide a sufficiently strong scattering signal. Oelke et al. used the Northern Hemisphere EASE-Grid Weekly Snow Cover and Sea Ice Extent Version 2 data set, available from NSIDC, to identify snow. This data set represents a digital version of National Oceanic and Atmospheric Administration (NOAA) / National Environmental Satellite, Data and Information Service (NESDIS) snow-cover charts, which were based on information from several visible-band satellites. For grid cells where SSM/I does not detect snow cover but the NOAA charts do, Oelke et al. assume a snow thickness of 3 cm. The NOAA charts are most useful at the beginning of the winter season and for the southern margin of snow cover. Erroneous SSM/I detection of snow, for example, in the middle of summer, is eliminated through comparison with the NOAA snow charts.

Soil Properties
The three major model layers use soil bulk density data derived from the pedon database in the International Geosphere-Biosphere Program's (IGBP) SoilData System (Global Soil Data Task 2000). This database generates maps of soil parameters at any depth and spatial resolution. Since the SoilData System accounts only for mineral soil types, Oelke et al. parameterized a percentage of organic soil (peat) for the top two major soil layers. They also extracted relative compositions of clay, silt, sand and gravel from the SoilData System for each grid cell. These concentrations are used to weight the thermal conductivities of fine-grained and coarse-grained soils in both frozen and thawed states (Kersten 1949). Finally, the University of New Hampshire Permafrost/Water Balance Model (Holden et al. 2001) provided daily soil water content, derived from a 20-year model climatology (1981-2000).

Processing Steps:

Oelke et al. calculated daily thaw depth and freezing depth using a finite-difference model for one-dimensional heat conduction with phase change. Details of this model are described in (Goodrich 1982). This model has shown excellent previous results for active-layer depth and soil temperatures when driven by well-known boundary conditions and forcing parameters at specific locations (Zhang, Osterkamp, and Stamnes 1996).

The researchers applied the frozen-ground model to the entire Arctic drainage area on a daily, 25-km, EASE-grid scale. They divided soil into three major layers: 0-30 cm, 30-80 cm, and 80-1500 cm, and calculated 54 model nodes ranging from a thickness of 10 cm near the surface to 1 m at 15-m depth. They determined thermal properties of mineral soils from soil dry-bulk density and water content according to Kersten (1949) and Lunardini (1988) for peat.

Oelke et al. chose initial temperatures for the model, based on a grid cell's permafrost classification from the International Permafrost Association (IPA) Circum-Arctic Map of Permafrost and Ground-Ice Conditions (Brown et al. 1998). They ran the model for a simulated 365 days to obtain realistic start conditions for temperatures at all model layers. All soil down to the lower model boundary was set to sub-freezing initial temperatures for the thaw depth. The upper model boundary represents the snow surface when snow cover is present. Daily mean air temperature defines the upper boundary condition. When snow is absent, the model's upper boundary is the ground surface, represented by surface temperature (°C). In this model, heat transfer within the snow layer is by conduction only, and snow densities and thermal conductivities are climatological. Oelke et al. assigned an initial geothermal heat flux for every grid cell using a temperature gradient of 3 K per 100 m. They ran the model for 20 years with this condition, and they used the resulting temperature gradient between the lowest two model layers together with the bottom-layer thermal conductivity to calculate the new geothermal heat flux.

Freezing depth was calculated by setting initial soil temperatures to above-freezing values at all depths. Freezing depth increases during winter and decreases with spring thawing.

Oelke et al. calculated thaw depth for areas identified from the IPA map as continuous permafrost, and the frozen parts of discontinuous and sporadic permafrost. Thaw depth was not determined for seasonally frozen-ground areas. In this case, the corresponding model grid cells were coded as "missing" with values of -9999.0. Freezing depth was calculated for areas of seasonally-frozen ground and for the non-frozen parts of sporadic and discontinuous permafrost. Freezing depth was not calculated for areas of continuous permafrost, based on the IPA map. Corresponding grid cell values were coded as missing (-9999.0) in this case. The researchers simulated thaw depth and freezing depth for the regions of discontinuous permafrost and for sporadic permafrost. In areas of continuous permafrost, Oelke et al. showed only simulated thaw depths; only one freezing depth was simulated for seasonally-frozen ground.

5. References and Related Publications:

Armstrong, R.L., and M.J. Brodzik. 2002. Northern Hemisphere EASE-Grid weekly snow cover and sea ice extent, Version 2. Boulder, CO: National Snow and Ice Data Center. CD-ROM.

Brodzik, M.J. 1997. EASE-Grid: A versatile set of equal-area projections and grids. Boulder, CO, USA: National Snow and Ice Data Center.

Brown, J., O.J. Ferrians, Jr., J.A. Heginbottom, and E.S. Melnikov. 1998. Digital Circum-Arctic map of permafrost and ground-ice conditions. In International Permafrost Association, Data and Information Working Group, comp. Circumpolar Active-Layer Permafrost System (CAPS), version 1.0. Boulder, CO: National Snow and Ice Data Center. CD-ROM.

Chang, A.T.C., J.L. Foster, and D.K. Hall. 1987. Nimbus-7 SMMR derived global snow cover parameters. Annals of Glaciology 9:39-44.

Global Soil Data Task. 2000. Global gridded surfaces of selected soil characteristics. International Geosphere-Biosphere Programme - Data and Information Services (IGBP-DIS). Oak Ridge National Laboratory (ORNL) Distributed Active Archive Center. Available online at

Goodrich, L.E. 1982. Efficient numerical technique for one-dimensional thermal problems with phase change. International Journal of Heat and Mass Transfer 21:615-621.

Holden, J.B., C.J. Vörösmarty, S. Frolking, and R. Lammers. 2001. A large-scale water balance model for permafrost terrain. Arctic and Alpine Research. In press.

Kalnay, E., M. Kanamitsu, R. Kistler, W. Collins, D. Deaven, L. Gandin, M. Iredell, S. Saha, G. White, J. Woolen, Y. Zhu, M. Chelliah, W. Ebisuzaki, W. Higgens, J. Janowiak, K.C. Mo, C. Ropelewski, J. Wang, A. Leetma, R. Reynolds, R. Jenne, and D. Joseph. 1996. The NCEP/NCAR 40-year reanalysis project. Bulletin of the American Meteorological Society 77:437-471.

Kersten, M.S. 1949. Laboratory research for the determination of the thermal properties of soils. Final report. University of Minnesota, Engineering Experiment Station.

Lunardini, V.J. 1988. Heat conduction with freezing and thawing. Monograph, 88-1. US Army Corps of Engineers Cold Regions Research and Engineering Laboratory.

Meteorological Service of Canada (MSC). 2000. Canadian Snow Data CD-ROM. Downsview, Ontario, Canada: Meteorological Service of Canada, CRYSYS Project, Climate Processes and Earth Observation Division.

Oelke, C., T. Zhang, M. Serreze, and R. Armstrong. 2003. Regional-scale modeling of soil freeze/thaw over the Arctic Drainage Basin. Journal of Geophysical Research 108(D10) (May).

Sturm, M., J. Holmgren, and G.E. Liston. 1995. A seasonal snow cover classification system for local to global applications. Journal of Climate 8:1261-1283.

Zhang, T., J.A. Heginbottom, R.G. Barry, and J. Brown. 2000. Further statistics on the distribution of permafrost and ground-ice in the Northern Hemisphere. Polar Geography 24(2):126-131.

Zhang, T., R.G. Barry, K. Knowles, J.A. Heginbottom, and J. Brown. 1999. Statistics and characteristics of permafrost and ground ice distribution in the Northern Hemisphere. Polar Geography 23(2):147-169.

Zhang, T., T.E. Osterkamp, and K. Stamnes. 1996. Influence of the depth hoar layer of the seasonal snow cover on the ground thermal regime. Water Resources Research 32(7):2075-2086.

6. Document Information:

List of Acronyms

The following acronyms are used in this document:

ARCSS-ADCC: Arctic System Science Data Coordination Center
ASCII: American Standard Code for Information Interchange
CALM: Circumpolar Active Layer Monitoring project
EASE-Grid: Equal-Area Scalable Earth Grid
FTP: file transfer protocol
IGBP: International Geosphere-Biosphere Programme
IDS: Interdisciplinary Study
IPA: International Permafrost Association
MSC: Meteorological Service of Canada
NASA: National Aeronautics and Space Administration
NCAR: National Center for Atmospheric Research
NCEP: National Center for Environmental Prediction
NESDIS: National Environmental Satellite, Data and Information Service
NOAA: National Oceanic and Atmospheric Administration
RMS: Root Mean Square
SMMR: Scanning Multichannel Microwave Radiometer
SSM/I: Special Sensor Microwave/Imager

Document Creation Date:

February 2003

Document Revision Date:

February 2003

Document Review Date:

February 2003

Document URL: