the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Modelling active layer thickness in mountain permafrost based on an analytical solution of the heat transport equation, Kitzsteinhorn, Hohe Tauern Range, Austria
Abstract. The active layer thickness (ALT) refers to the seasonal thaw depth of a permafrost body and is an essential parameter for natural hazard analysis, construction, land-use planning and the estimation of greenhouse gas emissions in periglacial regions. The aim of this study is to model the annual maximum thaw depth for determining ALT based on temperature data measured in four shallow boreholes (SBs, 0.1 m deep) in the summit region of the Kitzsteinhorn (Hohe Tauern Range, Austria, Europe). We set up our heat flow model with temperature data (2016–21) from a 30 m deep borehole (DB) drilled into bedrock at the Kitzsteinhorn north-face. For modeling purposes, we assume 1D conductive heat flow and present an analytical solution of the heat transport equation through sinusoidal temperature waves resulting from seasonal temperature oscillations (damping depth method). The model approach is considered successful: In the validation period (2019–21), modeled and measured ALT differed by only 0.1±0.1 m. We then applied the DB-calibrated model to four SBs and found that the modeled seasonal ALT maximum ranged between 2.5 m (SB 2) and 10.6 m (SB 1) in the observation period (2013–2021). Due to small differences in altitude (~ 200 m) within the study area, slope aspect had the strongest impact on ALT. To project future ALT deepening due to global warming, we integrated IPCC climate scenarios SSP1-2.6 and SSP5-8.5 into our model. By mid-century (~ 2050), ALT is expected to increase by 48 % at SB 2 and by 62 % at DB under scenario SSP1-2.6 (56 % and 128 % under scenario SSP5-8.5), while permafrost will no longer be present at SB 1, SB 3 and SB 4. By the end of the century (~ 2100), permafrost will only remain under scenario SSP1-2.6 with an ALT increase of 51 % at SB 2 and of 69 % at DB.
- Preprint
(62929 KB) - Metadata XML
- BibTeX
- EndNote
Status: open (until 31 Jul 2024)
-
RC1: 'Comment on egusphere-2023-3006', Anonymous Referee #1, 18 Jun 2024
reply
General comments.
This study is on active layer thickness in an alpine rock permafrost environment with borehole temperature measurements and analytical modeling of the heat diffusion equation. The measurements and data set used is very interesting and the underlying problem of using shallow ground temperature measurements (here at 0.1 m) for making predictions of ALT and permafrost warming is very important especially in the context of climate change and its severe impacts for alpine regions.
However, the thermal modeling approach used is not new or novel, and has very limited applicability for active layer modeling. The equations described are well-known analytic solutions to the heat conduction (diffusion) equation and can be found in textbooks on the subject, see e.g. Carslaw and Jaeger (Conduction of heat in solids, 1959), Williams (The Frozen Earth, 1989), Woo (Permafrost Hydrology, 2012), and others. It is severely limited for active layer and permafrost modeling because it does not represent many of the critical freezing and thawing processes which control heat propagation in the subsurface.
Many of the assumptions imposed are mentioned in various parts of the manuscript and include neglecting latent heat exchange, advective heat carried by water migration (both lateral and vertical), snow insulation and snowmelt, and heterogeneity in subsurface thermal and hydraulic properties. Critically, even for homogenous ground, latent heat, thermal conductivity and heat capacity change due to variations in liquid-ice-air phase saturation in the active layer, and these changes need to be accounted for in order to be able to calculate heat propagation through the active layer and to the permafrost table and below. A possible exception could be if the ground is extremely dry. Or if the interest is only long-term and deep temperature responses, say below the depth of zero annual amplitude. Even if the intact/unfractured bedrock contains little pore space and hence may be suitably represented by heat conduction, there are evidently fractures feeding the subsurface with unfrozen water. This is noted in the discussion in Section 5.1 and is clearly seen by the prolongation of temperature near the freezing point (the zero-degree curtain) for the thermal sensors located at 3 m depth in Fig 3a.
Several analytical and semi-analytical approaches have been developed to address the limitations of the analytical solution to the heat conduction equation, other than full numerical solutions of the coupled system of differential equations for heat transport with water flow. For some background and an overview, please see the excellent papers by Riseborough et al. (2008) and Kurylyk et al. (2014). Perhaps some of those more robust methods could be employed to make predictions on ALT based on the shallow boreholes, or even better, numerical models which incorporate even fewer limitations.
Riseborough, D., Shiklomanov, N., Etzelmüller, B., Gruber, S., Marchenko, S., 2008. Recent advances in permafrost modeling. Permafrost and Periglacial Processes 19, 137–156. https://doi.org/10.1002/ppp.615
Kurylyk, B.L., McKenzie, J.M., MacQuarrie, K.T.B., Voss, C.I., 2014. Analytical solutions for benchmarking cold regions subsurface water flow and energy transport models: One-dimensional soil thaw with conduction and advection. Advances in Water Resources 70, 172–184. https://doi.org/10.1016/j.advwatres.2014.05.005
Another problem is the biased implementation of the NSE as an error metric, which seems to be based on modeled values after discarding values which perform poorly (Fig 5a). This is a biased filtering and is probably the reason high values of the NSE metric are obtained. Clearly, if all the modeled data are used, including the freeze-up period, and compared against the measurements, the NSE metric will be much lower.
Finally, the interpretation of the ALT projections for the shallow boreholes SB 2, 3 and 4 in Section 4.3 and 5.3 is questionable. These boreholes seem to be grouped near a peak and the spatial distance between them only seems to be a few 10s of meters. Despite their proximity, the prediction on the ALT/PF table are very different (Table 5), from no PF (“infinite” ALT) in SB3 to relatively shallow ALT of 2-3 meters in SB2. It is unrealistic to infer that the permafrost table varies as extremely as this. The result is probably an artifact of the method which assumes 1D conduction only, and that the ALT prediction is based solely on a single near-surface temperature sensor. In reality of course, lateral subsurface heat transport distributes the thermal response and the permafrost extent will be smoother even if differential warming occurs on different faces of the mountain surface.
Detailed comments.
L125-130. Has the thermal effect of the concrete annulus been evaluated? It could be that the thermal properties of concrete have an effect on heat conduction between the rock and the brass segments and potentially bias sensor measurements. Also, there could be effects of vertical heat conduction in the concrete annulus which may differ from that of the bedrock.
L130-135. Is the surface smooth bedrock or is there overburden and unconsolidated material? It is difficult to judge from the photograph in Fig 2b, and not entirely clear how the depth 0.1 m is obtained, because it seems surface roughness could easily vary within this range of a few dm.
L250-255. Water saturation generally increases with depth and, more critically, changes over time with variable infiltration from snowmelt, rainfall, and seasonal freezing/thawing.
L195-200. The problem with the damping depth calibration approach used is that it applies to the thaw depth of the time period considered. With the climate warming scenarios imposed later, the thaw depth changes, yielding the calibration obsolete. The calibration-validation exercise simply indicates no significant changes occur within the time period.
Fig 3a shows the zero-degree isotherm, indicating presence of water and showing that latent heat exchange is playing a significant role. Thus, the assumptions used for the model are questionable.
Eqn 7. It is not clear how temperature comes in to play for the equation for d_phase?
L195-200. The NSE reference is not included in the bibliography. Also, values can be negative, indicating predictions are worse than the mean of the measured data, i.e., not capturing the variability of measurements.
Fig 4b. Unclear why only two years are shown. Please make the markers more visible, currently they are difficult to see. The colors for 2017 and 2018 are inconsistent with Fig 4c. Is it correct there are only 3 measurements in the interval 0-2 m; if so, then a linear regression is of course expected to fit well. Same goes for Fig 4c.
It would be useful to show a trumpet diagram of the temperature measurements with depth (similar to Fig 4a) for the deep borehole with systematically selected times, e.g., once a month. This serves to present data without interpolations as in Fig 3b and A1, and will help visualize the active layer depth, the depth of zero annual amplitude, if there is an inflection point, and if the temperature trend with depth indicates a transition to the permafrost base. This can readily be combined with annual averages for each depth.
Fig 5a. The modeled values can of course be calculated at all depths but the measurements are only made at specific depths. Therefore, it would be reasonable to plot lines for the simulations and markers for the measurements at their corresponding depths. This would help clarify the comparison between modeled results and measurements. Also, it would be helpful to clarify that the values in brackets indicate the ALT.
Fig 5a and Section 4.2. It seems by ignoring the data which does not fit the model is the reason high NSE and RMSE values are obtained – but this is by biasing the measurements and as such highlights the limitations of the model approach.
Fig 5a. The main reason the freeze-up is not represented well by the model is the absence of latent heat.
Citation: https://doi.org/10.5194/egusphere-2023-3006-RC1
Viewed
HTML | XML | Total | BibTeX | EndNote | |
---|---|---|---|---|---|
190 | 53 | 19 | 262 | 13 | 12 |
- HTML: 190
- PDF: 53
- XML: 19
- Total: 262
- BibTeX: 13
- EndNote: 12
Viewed (geographical distribution)
Country | # | Views | % |
---|
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1