Subsea permafrost stability is the key to whether pre-performed methane sequestered in hydrate deposits escapes to the overlying strata. By making use of the 1D numerical modeling and field data, we analyze the capabilities of the time-domain (transient) electromagnetic method (TDEM) when being applied for subsea permafrost mapping, and study the effect of the background resistivity structure on the inversion models’ accuracy for a series of settings typical for the East Siberian Arctic Shelf—the broadest and shallowest shelf in the world ocean, which represents more than 70% of the subsea permafrost. The synthetic response analysis included the construction of a series of resistivity models corresponding to different settings (presence/absence of ice-bonded permafrost layer, different position of its top and bottom boundaries, different width and thickness of thawed bodies or taliks, variable seawater depth and its resistivity), and calculation of synthetic apparent resistivity responses used to assess their sensitivity to changes in the target parameters of the resistivity structure. This was followed by regularized inversion of synthetic responses and comparing resulting models with original (true) ones, which allowed us to understand the possible uncertainties in the geometry and resistivity of the reconstructed permafrost layer, depending on seawater depth and unfrozen layer thickness, as well as confirm the overall efficacy of TDEM technology for the subsea permafrost imaging. That is crucially important for understanding the current state of the subsea permafrost-hydrate system and possible future dynamics.

Catastrophic megaearthquakes (M > 8) occurring in the subduction zones are among the most devastating hazards on the planet. In this paper we discuss the seismic cycles of the megathrust earthquakes and propose a blockwise geomechanical model explaining certain features of the stress-deformation cycle revealed in recent decades from seismological and satellite geodesy (GNSS) observations. Starting with an overview of the so-called keyboard model of the seismic cycle by L. Lobkovsky, we outline mathematical formalism describing the motion of seismogenic block system assuming viscous rheology beneath and between the neighboring elastic blocks sitting on top of the subducting slab. By summarizing the GNSS-based evidence from our previous studies concerning the transient motions associated with the 2006–2007 Simushir earthquakes, 2010 Maule earthquake, and 2011 Tohoku earthquake, we demonstrate that those data support the keyboard model and reveal specific effect of the postseismic oceanward motion. However, since the seismogenic blocks in subduction systems are mostly located offshore, the direct analysis of GNSS-measured displacements and velocities is hardly possible in terms of the original keyboard model. Hence, the generalized two-segment keyboard model is introduced, containing both frontal offshore blocks and rear onshore blocks, which allows for direct interpretation of the onshore-collected GNSS data. We present a numerical computation scheme and a series of simulated data, which exhibits the consistency with measured motions and enables estimating the seismic cycle characteristics, important for the long-term earthquake forecasting.

A correlation is observed between changes in the level of Earth’s seismic activity and increments of the atmospheric methane concentration over the past 40 years. Trigger mechanisms are proposed for methane emissions and glacier collapse in polar regions. These mechanisms are due to deformation waves caused by large earthquakes in subduction zones located near the polar regions: the Aleutian and Kuril–Kamchatka subduction zones, closest to the Arctic, and the Antarctica–Chilean and Tonga–Kermadec–Macquarie subduction zones. Disturbances of the lithosphere are transmitted over the distances of 3000–4000 km and more at a speed of about 100 km/year. Additional associated stresses come to the Arctic and Antarctica several decades after the occurrence of large earthquakes. In the Arctic zone, additional stresses affect the low-permeability structure of gas bearing sedimentary strata, causing increased methane emission and climate warming. In West Antarctica, deformation waves could trigger the acceleration and intensive collapse of West Antarctic glaciers, which has been observed since the 1970s. These waves are also capable of activating dormant volcanoes located under the sheet glaciers of West Antarctica, leading to an increase in heat flux, to the melting of ice at the glaciers’ base, and to their accelerated sliding towards the ocean, as is happening with the Thwaites Glacier.

Late Mesozoic and Cenozoic geodynamics of the Arctic region is discussed in the context of possible mechanisms which provide multistage cyclic transformations and transport of carbon through crust and mantle. Geodynamic processes control the abiogenic generation of hydrocarbons and the patterns of their localization. Possible mechanisms of abiotic hydrocarbon generation are explained in the context of carbon transport from subduction zones to rifts and serpentinization of ultramafic rocks in the rifts in the case of the Laptev Sea and Gakkel Ridge areas. The carbon of shallow crust origin migrates with encapsulated fragments of marine sediments which are consumed in the Pacific subduction zone where they become destroyed and transformed by different chemical and physical processes. The resulting C-species are involved in mantle convection flows and reach the continental rifts of the Laptev Sea and the Gakkel mid-ocean ridge. Thus, the hydrocarbons formed in the crust and in the mantle acquire signatures of abiotic origin. According to the authors, the scale of manifestation of abiogenic methanogenesis in the lower parts of the lithosphere and in the upper mantle is not so wide. Numerous small (mm and fractions of the mm) particles of exogenous matter and dispersed carbon pulled into the mantle can only form a stable crustal geochemical plume that propagates in the plane of movement of convective flows. Indirectly, the scale of manifestation of this process can be judged by the volumes of degassing of hydrocarbon and carbon dioxide gases, as well as hydrogen and its compounds in the rift systems of the earth’s crust, which are extremely insignificant. However, in the cold seas of the Eastern Arctic, massive emissions of bubble methane of mixed genesis were found. As shown in the literature, the range of variability of stable isotopes of carbon and ¹⁴C of methane in certain areas of discharge associated with rifting demonstrates values (anomalously heavy ¹³C, and young ¹⁴C) that can be considered as examples of presumably abiogenic origin. Our work is mostly theoretical and suggests further discussion and improvement of the mechanism of formation of abiogenic hydrocarbons and the processes of their transformation.

We discussed the peculiarities of the seismic cycle in Aleutian subduction zone, characterized by an oblique subduction setting. It was shown that the orientation of the plate convergence vector relative to the subduction zone axis can have a significant impact on the preparation and occurrence of the largest earthquakes in subduction zones. In particular, from the analysis of the seismic activity occurring in the western part of the Aleutian island arc, it was found that the seismic cycles here are shorter than in the eastern part of the arc. It was revealed that the strongest earthquakes, repeating in the same areas of the western part of the Aleutian subduction zone, differ both in magnitude and length of the fault zone. Taking into account the oblique subduction setting, we proposed the keyboard model of the largest megathrust earthquakes generation as a mechanism potentially capable of explaining the reduction in the seismic cycle duration and noticeable differences in the spatial extent and localization of the fault zones of events with similar magnitudes occurring in the same segment of the western half of the Aleutian subduction zone.