Earth’s Cryosphere, 2024, Vol. XXVIII, No. 5, p. 32-49.
METHODS OF CRYOSPHERIC RESEARCH
NUMERICAL MODELLING OF SUBMARINE PERMAFROST ON THE EURASIAN ARCTIC SHELF CONSIDERING MODERN CLIMATE ZONALITY
Yu.Yu. Smirnov1,2,*, T.V. Matveeva1, N.A. Shchur1,3, A.A. Shchur1, A.V. Bochkarev1
1 All-Russia Research Institute for Geology and Mineral Resources of the Ocean (VNIIOkeangeologia), Angliysky prosp. 1, St. Petersburg, 190121 Russia
2 Russian State Hydrometeorological University, Voronezhskaya St. 79, St. Petersburg, 192007 Russia
3 Peter the Great St. Petersburg Polytechnic University, Politekhnicheskaya St. 29, St. Petersburg, 195251 Russia
*Corresponding author; e-mail: y.y.smirnov@mail.ru
The study and prediction of subsea permafrost conditions hold profound significance in the context of global climate change, the formation of gas hydrates, and the assessment of natural hazards arising from permafrost thawing. This study presents a forecast of relict subsea permafrost distribution across the Eurasian Arctic shelf through the resolution of a one-dimensional, nonstationary Stefan problem using finite difference methods, facilitated by contemporary software tools and libraries. The mathematical description of the model is presented. Particular attention is paid to the boundary conditions of the model and their synthesis. Calculations take into account the climate clustering of the extensive Eurasian shelf, considering zonal variations in both near-surface air temperature and the temperature and salinity of subsea waters. For the first time, model estimates of the position of the upper boundary of subsea permafrost are presented, and its position is compared to drilling data. Analysis of the influence of boundary conditions on the simulation results and comparison of our results with seismoacoustic and drilling data confirm the good quality of the model. The evolution of subsea relict permafrost in time is represented by a set of three maps for the Late Pleistocene–Holocene time spans. As a result, a high-resolution global map of the thickness of underwater permafrost on the Eurasian shelf was created. Detailed modelling and forecast mapping of the distribution of permafrost table was carried out in well-studied areas. According to our forecasts, a wide distribution of subsea relict permafrost is predicted throughout the Eurasian Arctic shelf. The greatest thickness of relict frozen deposits is confined to the littoral zone of the Laptev and East Siberian seas, as well as the offshore areas of the New Siberian Islands. In the Pechora and Chukchi seas, mostly poorly developed, close to complete degradation, relict subsea permafrost is predicted.
Keywords: Eurasian Arctic shelf, submarine permafrost, Python, Julia, numerical permafrost modelling, forecast submarine permafrost mapping.
Recommended citation: Smirnov Yu.Yu., Matveeva T.V., Shchur N.A., Shchur A.A., Bochkarev A.V., 2024. Numerical modelling of submarine permafrost on the Eurasian Arctic shelf considering modern climate zonality. Earth’s Cryosphere XXVIII (5), 32–49.
Received February 5, 2023
Revised May 27, 2024
Accepted July 23, 2024
Translated by A.V. Muravyev