Model:Lake-Permafrost with Subsidence: Difference between revisions
From CSDMS
No edit summary |
No edit summary |
||
| Line 15: | Line 15: | ||
|Categories=Terrestrial | |Categories=Terrestrial | ||
|Spatial dimensions=1D | |Spatial dimensions=1D | ||
|One-line model description=1-D lake-permafrost thermal model with subsidence. | |One-line model description=1-D lake-permafrost thermal model with subsidence. | ||
|Extended model description=Models temperature of 1-D lake-permafrost system through time, given input surface temperature and solar radiation. Model is fully implicit control volume scheme, and cell size can vary with depth. Thermal conductivity and specific heat capacity are dependent on cell substrate (% soil and % ice) and temperature using the apparent heat capacity scheme where freezing/thawing occurs over a finite temperature range and constants are modified to account for latent heat. Lake freezes and thaws depending on temperature; when no ice is present lake is fully mixed and can absorb solar radiation. Upper 10 m substrate contains excess ice and, if thawed, can subside by this amount (lake then deepens by amount of subsidence). "Cell type" controls whether cell has excess ice, only pore space ice, or is lake water. | |Extended model description=Models temperature of 1-D lake-permafrost system through time, given input surface temperature and solar radiation. Model is fully implicit control volume scheme, and cell size can vary with depth. Thermal conductivity and specific heat capacity are dependent on cell substrate (% soil and % ice) and temperature using the apparent heat capacity scheme where freezing/thawing occurs over a finite temperature range and constants are modified to account for latent heat. Lake freezes and thaws depending on temperature; when no ice is present lake is fully mixed and can absorb solar radiation. Upper 10 m substrate contains excess ice and, if thawed, can subside by this amount (lake then deepens by amount of subsidence). "Cell type" controls whether cell has excess ice, only pore space ice, or is lake water. | ||
}} | }} | ||
{{Start model keyword table}} | |||
{{End a table}} | |||
{{Model technical information | {{Model technical information | ||
|Supported platforms=Windows | |Supported platforms=Windows | ||
| Line 26: | Line 28: | ||
|End year development=2009 | |End year development=2009 | ||
|Model availability=As code | |Model availability=As code | ||
|Program license type=GPL v2 | |Program license type=GPL v2 | ||
|Memory requirements=-- | |Memory requirements=-- | ||
| Line 32: | Line 33: | ||
}} | }} | ||
{{Input - Output description | {{Input - Output description | ||
|Describe input parameters=Daily average solar radiation for location (at surface). | |Describe input parameters=Daily average solar radiation for location (at surface). | ||
|Other input format=Matlab .mat file | |Other input format=Matlab .mat file | ||
|Describe output parameters=Matlab variables, Matlab graphs | |Describe output parameters=Matlab variables, Matlab graphs | ||
| Line 41: | Line 42: | ||
}} | }} | ||
{{Process description model | {{Process description model | ||
|Describe processes represented by the model=Heat conduction in permafrost, lake ice growth-decay, permafrost subsidence due to excess ice | |Describe processes represented by the model=Heat conduction in permafrost, lake ice growth-decay, permafrost subsidence due to excess ice | ||
|Describe key physical parameters and equations=Heat conduction equations, lake ice growth-decay equations | |Describe key physical parameters and equations=Heat conduction equations, lake ice growth-decay equations | ||
|Describe length scale and resolution constraints=Near surface, length scale should be 5-10 cm if subsidence occurring. Deeper in permafrost, 1+ m okay. | |Describe length scale and resolution constraints=Near surface, length scale should be 5-10 cm if subsidence occurring. Deeper in permafrost, 1+ m okay. | ||
|Describe time scale and resolution constraints=Model needs daily (or smaller) timesteps. | |Describe time scale and resolution constraints=Model needs daily (or smaller) timesteps. | ||
|Describe any numerical limitations and issues=Runs slowly - iterates implicit scheme. Some sort of matrix algebra might improve speed. | |Describe any numerical limitations and issues=Runs slowly - iterates implicit scheme. Some sort of matrix algebra might improve speed. | ||
}} | }} | ||
{{Model testing | {{Model testing | ||
|Describe available calibration data sets=-- | |Describe available calibration data sets=-- | ||
|Describe available test data sets=Permafrost model can be tested with ground temperature time series at several depths. Lake model can be tested with known ice depths (though snow, not accounted for in model, can affect depth by several 10s of cm). | |Describe available test data sets=Permafrost model can be tested with ground temperature time series at several depths. Lake model can be tested with known ice depths (though snow, not accounted for in model, can affect depth by several 10s of cm). | ||
|Describe ideal data for testing=Lake temperature time series, known lake depth, underlying permafrost time series. Surface temperature time series as input (not air temperature). | |Describe ideal data for testing=Lake temperature time series, known lake depth, underlying permafrost time series. Surface temperature time series as input (not air temperature). | ||
}} | }} | ||
{{Users groups model | {{Users groups model | ||
| Line 56: | Line 57: | ||
}} | }} | ||
{{Documentation model | {{Documentation model | ||
|Provide key papers on model if any=Model based on previous models by Ling and Zhang (2004), Plug and West (2008), Liston and Hall (1995). Model described in Matell et al. (in preparation) and Matell (2009, masters thesis). | |Provide key papers on model if any=Model based on previous models by Ling and Zhang (2004), Plug and West (2008), Liston and Hall (1995). Model described in Matell et al. (in preparation) and Matell (2009, masters thesis). | ||
|Manual model available=No | |Manual model available=No | ||
}} | }} | ||
| Line 64: | Line 65: | ||
|CCA component=No but possible | |CCA component=No but possible | ||
|IRF interface=No but possible | |IRF interface=No but possible | ||
|CMT component=Not yet | |||
}} | }} | ||
{{Start coupled table}} | |||
{{End a table}} | |||
{{End headertab}} | {{End headertab}} | ||
<!-- Edit the part above to update info on other papers --> | <!-- Edit the part above to update info on other papers --> | ||
Revision as of 17:33, 23 June 2011
