Subsea Permafrost Research Articles

Vulnerability of Arctic-Boreal methane emissions to climate change

The rapid warming of the Arctic-Boreal region has led to the concern that large amounts of methane may be released to the atmosphere from its carbon-rich soils, as well as subsea permafrost, amplifying climate change. In this review, we assess the various sources and sinks of methane from northern high latitudes, in particular those that may be enhanced by permafrost thaw. The largest terrestrial sources of the Arctic-Boreal region are its numerous wetlands, lakes, rivers and streams. However, fires, geological seeps and glacial margins can be locally strong emitters. In addition, dry upland soils are an important sink of atmospheric methane. We estimate that the net emission of all these landforms and point sources may be as much as 48.7 [13.3–86.9] Tg CH4 yr−1. The Arctic Ocean is also a net source of methane to the atmosphere, in particular its shallow shelves, but we assess that the marine environment emits a fraction of what is released from the terrestrial domain: 4.9 [0.4–19.4] Tg CH4 yr−1. While it appears unlikely that emissions from the ocean surface to the atmosphere are increasing, now or in the foreseeable future, evidence points towards a modest increase from terrestrial sources over the past decades, in particular wetlands and possibly lakes. The influence of permafrost thaw on future methane emissions may be strongest through associated changes in the hydrology of the landscape rather than the availability of previously frozen carbon. Although high latitude methane sources are not yet acting as a strong climate feedback, they might play an increasingly important role in the net greenhouse gas balance of the Arctic-Boreal region with continued climate change.

Seismic Imaging of the Arctic Subsea Permafrost Using a Least-Squares Reverse Time Migration Method

High-resolution seismic imaging allows for the better interpretation of subsurface geological structures. In this study, we employ least-squares reverse time migration (LSRTM) as a seismic imaging method to delineate the subsurface geological structures from the field dataset for understanding the status of Arctic subsea permafrost structures, which is pertinent to global warming issues. The subsea permafrost structures in the Arctic continental shelf, located just below the seafloor at a shallow water depth, have an abnormally high P-wave velocity. These structural conditions create internal multiples and noise in seismic data, making it challenging to perform seismic imaging and construct a seismic P-wave velocity model using conventional methods. LSRTM offers a promising approach by addressing these challenges through linearized inverse problems, aiming to achieve high-resolution, subsurface imaging by optimizing the misfit between the predicted and the observed seismic data. Synthetic experiments, encompassing various subsea permafrost structures and seismic survey configurations, were conducted to investigate the feasibility of LSRTM for imaging the Arctic subsea permafrost from the acquired seismic field dataset, and the possibility of the seismic imaging of the subsea permafrost was confirmed through these synthetic numerical experiments. Furthermore, we applied the LSRTM method to the seismic data acquired in the Canadian Beaufort Sea (CBS) and generated a seismic image depicting the subsea permafrost structures in the Arctic region.

Passive Seismology: Lightweight and Rapid Detection of Arctic Subsea and Sub‐Aquatic Permafrost

AbstractLow sea levels during the last Ice Age exposed millions of square kilometers of Arctic shelves which have been subsequently submerged, creating subsea permafrost. In onshore settings, permafrost can also exist beneath water bodies such as coastal lagoons, rivers, and thermokarst lakes. We explored passive seismology as a method for mapping unfrozen sediment thickness above subsea and sub‐aquatic permafrost. We present passive seismic data collected with the Mobile Ocean Bottom Seismic Instrument (MOBSI) from the Beaufort Sea near Tuktoyaktuk in Canada, Ivashkina Lagoon on the Bykovsky Peninsula, as well as a lake and river in the Lena Delta, Siberia, Russia. We use borehole data and frost probe measurements to identify permafrost‐related H/V measurement peaks and calibrate shear wave velocities for frequency‐to‐depth conversion. We employ the shortest path and maximum signal amplitude to connect peaks and generate geological profiles. The MOBSI detected the ice‐bonded permafrost table beneath the Beaufort Sea, as well as beneath a Siberian lake and lagoon. At Tuktoyaktuk, an ocean bottom seismometer revealed a 5% scatter about the peak frequency for three‐minute time windows and over 8 hr of recording time. With peak frequencies ranging from 4.9 ± 0.2 Hz to 27.6 ± 1.4 Hz, the depth to subsea permafrost ranged from 1.4 ± 0.1 m bsl at the shoreline to 14.0 ± 0.4 m bsl 240 m offshore. Given an accurate shear wave velocity, our findings highlight that MOBSI deployment times as short as 3 min are adequate for detecting Arctic subsea and sub‐aquatic permafrost.

Evidence for subsea permafrost in subarctic Canada linked to submarine groundwater discharge

Evidence for subsea permafrost in subarctic Canada linked to submarine groundwater discharge

Deep-learning viscoelastic seismic inversion for mapping subsea permafrost

Marine seismic surveys can be used to map ice-bearing subsea permafrost on a large scale. However, current seismic processing technologies have limited capacity to image permafrost distribution at depth, mainly due to the low sensitivity of primary reflections and refractions to the velocity inversion found at the base of permafrost. Guided waves and multiples are more sensitive to the velocity variations below the top of permafrost, but they remain challenging to use in physics-based inversion approaches. A deep-learning-based seismic inversion has the potential to improve seismic imaging below the top of permafrost by automatically extracting information from all wave modes. We develop a multi-input neural network (NN) to estimate seismic velocities from marine seismic data. The network is trained on synthetic data generated from the representative distributions of the seismic properties of subsea permafrost. We find that our network can image large velocity contrasts and reversals in depth, typical of subsea permafrost. We use our network to estimate the P- and S-wave velocity and Q-factor models from a seismic line in the Beaufort Sea. The NN indicates highly discontinuous subsea permafrost with variable thickness in the area. Our work shows that deep-learning-based seismic inversion could become a cost-effective technology to map the distribution of subsea permafrost on a large scale and, more generally, high-velocity geologic layers located in shallow waters.

Glacial isostatic adjustment reduces past and future Arctic subsea permafrost

Sea-level rise submerges terrestrial permafrost in the Arctic, turning it into subsea permafrost. Subsea permafrost underlies ~ 1.8 million km2 of Arctic continental shelf, with thicknesses in places exceeding 700 m. Sea-level variations over glacial-interglacial cycles control subsea permafrost distribution and thickness, yet no permafrost model has accounted for glacial isostatic adjustment (GIA), which deviates local sea level from the global mean due to changes in ice and ocean loading. Here we incorporate GIA into a pan-Arctic model of subsea permafrost over the last 400,000 years. Including GIA significantly reduces present-day subsea permafrost thickness, chiefly because of hydro-isostatic effects as well as deformation related to Northern Hemisphere ice sheets. Additionally, we extend the simulation 1000 years into the future for emissions scenarios outlined in the Intergovernmental Panel on Climate Change’s sixth assessment report. We find that subsea permafrost is preserved under a low emissions scenario but mostly disappears under a high emissions scenario.

Biodegradation of Ancient Organic Carbon Fuels Seabed Methane Emission at the Arctic Continental Shelves

AbstractThis study explores the carbon stability in the Arctic permafrost following the sea‐level transgression since the Last Glacial Maximum (LGM). The Arctic permafrost stores a significant amount of organic carbon sequestered as frozen particulate organic carbon, solid methane hydrate and free methane gas. Post‐LGM sea‐level transgression resulted in ocean water, which is up to 20°C warmer compared to the average annual air mass, inundating, and thawing the permafrost. This study develops a one‐dimensional multiphase flow, multicomponent transport numerical model and apply it to investigate the coupled thermal, hydraulic, microbial, and chemical processes occurring in the thawing subsea permafrost. Results show that microbial methane is produced and vented to the seawater immediately upon the flooding of the Arctic continental shelves. This microbial methane is generated by the biodegradation of the previously frozen organic carbon. The maximum seabed methane flux is predicted in the shallow water where the sediment has been warmed up, but the remaining amount of organic carbon is still high. It is less likely to cause seabed methane emission by methane hydrate dissociation. Such a situation only happens when there is a very shallow (∼200 m depth) intra‐permafrost methane hydrate, the occurrence of which is limited. This study provides insights into the limits of methane release from the ongoing flooding of the Arctic permafrost, which is critical to understand the role of the Arctic permafrost in the carbon cycle, ocean chemistry and climate change.

Subsea permafrost and associated methane hydrate stability zone: how long can they survive in the future?

Subsea permafrost and associated methane hydrate stability zone: how long can they survive in the future?

Распространение субаквальной мерзлоты в море Лаптевых по данным сейсморазведки методом преломленных волн

For the first time, for an area of 454 thousand km2 of the Laptev Sea, processing and comprehensive analysis of a large volume of records of the first arrivals of refracted waves of common shot gathers was carried out along 113 CDP seismic lines of JSC MAGE with a total length of about 20.7 thousand km. Fundamentally new information was obtained on the state of the subsea cryolithozone, and the boundary between the predominant distribution of frozen (Southern zone) and thawed (Northern zone) ground was identified. It was substantiated that a number of identified through taliks in the Southern zone have an endogenous genesis and are associated with large disjunctive displacements. High seismic activity in the central part of the Laptev Sea improves faults permeability, which contributes to the activation of subvertical migration of deep fluid flows and enhances the role of the endogenous factor in the degradation of permafrost (frozen ground) and gas hydrates. A large thawed zone was discovered in the North Western area of JSC MAGE activities, approaching close to the coast of Taimyr. A comparison of the results obtained with the data of the stratigraphic well DL-1, drilled in 2022 by FSBI VSEGEI and JSC Rosgeologia near the eastern border of the Laptev Sea, showed their complete correspondence, clearly indicating completed degradation of subsea permafrost in a significant part of the seas of Eastern Siberia.

The hypothesis of the shape of the permafrost in Hornsund, Spitsbergen and the potential impact of its degradation on the Arctic

The areas covered by permafrost in the polar regions are vulnerable to rapid changes in the current climate. The well-studied near-surface active layer and permafrost zone are in contrast to the unknown exact shape of the bottom permafrost boundary. Therefore, the entire shape of permafrost between the upper and lower boundaries is not identified with sufficient accuracy. Since most of the factors affecting deep cryotic structures are subsurface in nature, their evolution in deeper layers is also relatively unclear. Here, we propose a hypothesis based on the results of geophysical studies regarding the shape of the permafrost in the coastal area of Svalbard, Southern Spitsbergen. In the article, we emphasize the importance of recognizing not only the uppermost active layer but also the bottom boundary of permafrost along with its transition zone, due to the underestimated potential role of its continuity in observing climate change. The lower permafrost boundary is estimated to range from 70 m below the surface in areas close to the shore to 180 m inland, while a continuous layer of an entirely frozen matrix can be identified with a thickness between 40 m and 100 m. We also hypothesized the presence of the possible subsea permafrost in the Hornsund. The influence of seawater intrusions, isostatic uplift of deglaciated areas, and surface-related processes that affect permafrost evolution may lead to extensive changes in the hydrology and geology of the polar regions in the future. For all these reasons, monitoring, geophysical imaging and understanding the characteristics and evolution of deep permafrost structures requires global attention and scientific efforts.

Satellite-Based Evaluation of Submarine Permafrost Erosion at Shallow Offshore Areas in the Laptev Sea

Large areas of the seafloor in the Laptev Sea consist of submarine permafrost, which has experienced intense degradation over the last decades and centuries. Thermal abrasion of the submarine permafrost results in upward advection of suspended matter, which could reach the surface layer in shallow areas. This process is visually manifested through increased turbidity of the sea surface layer, which is regularly detected in optical satellite imagery of the study areas. In this study, satellite data, wind and wave reanalysis, as well as in situ measurements are analyzed in order to reveal the main mechanisms of seafloor erosion in shallow areas of the Laptev Sea. We describe the synoptic variability in erosion at the Vasilyevskaya and Semenovskaya shoals in response to wind and wave conditions. Finally, using reanalysis data, daily suspended matter flux from this area was evaluated during ice-free periods in 1979–2021, and its seasonal and inter-annual variabilities were described. The obtained results contribute to our understanding of subsea permafrost degradation, the sediment budget, and carbon and nutrient cycles in the Laptev Sea.

In situ bottom sediment temperatures in the Siberian Arctic seas: Current state of subsea permafrost in the Kara sea vs laptev and East Siberian seas

Features of sediment temperature on the shelf and continental slope areas of the Russian Arctic seas and its physical properties are important for understanding the current state of subsea permafrost and the gas hydrates stability zone. New data are reported for the Kara Sea region where the bottom sediment temperatures are influenced by warming effects from great Siberian rivers and the Atlantic currents. The data collected during marine expeditions in 2019–2022 are combined with results of earlier marine studies, drilling operations, and geophysical surveys in the Laptev and East Siberian seas, in order to identify major trends of in situ temperature and properties distribution of bottom sediments in the Russian Arctic region.Most (85%) of bottom sediments in the Kara Sea shelf, as well as in the Laptev and East Siberian shelves, consist of water-saturated silty clay and silt with rather uniform particle size distribution. The obtained thermal conductivity and heat capacity values for the Kara Sea sediments agree with the values of 1.0 W/(m·K) and 2900 kJ/m3, respectively, obtained previously from other Arctic seas. Thermal conductivity becomes up to 40% higher depthward from 0 to 2 m subbottom depth, possibly, because of lower moisture content and porosity in more lithified sediments.The bottom sediment temperatures in the Arctic seas are distributed unevenly, especially in the Kara Sea shelf (from +5.0 °C in the west to −1.4 °C in the east), where the high sediment temperatures in the western and central parts of the Kara Sea being due to the effect of warm water inputs. The distribution of bottom sediment temperatures correlates well with distribution of relic subsea permafrost. Ice-bearing permafrost in the Siberian Arctic shelf extends from the shoreline till sea depths of 80–100 m, within the respective offshore distances of ∼800–1000 km in some areas, but permafrost remnants may exist locally at sea depths within 120 m. Buried 100–600 m thick continuous subsea permafrost may occur in the Kara, Laptev, and East Siberian shelves under unfrozen (cryotic) saline shallow sediments. However, subsea permafrost is discontinuous and sporadic at sea depths ∼70 m and more. Thus, the bottom sediment temperature features in the Arctic seas can be used as a proxy of subsea permafrost extent contenting intrapermafrost and subpermafrost gas and gas hydrate accumulations.

Understanding seismic velocity variations of subsea permafrost: A sensitivity study

Monitoring the change in permafrost conditions and distribution is crucial for forecasting global warming. As seismic velocities increase with ice content, marine seismic surveys can map the top of ice-bearing subsea permafrost on a large scale. However, conventional seismic methods cannot fully map internal velocity variations linked to changes in ice content because the removal of guided waves and multiples is challenging for reflection processing. Nevertheless, these arrivals carry information about velocity variations with depth. We investigate if a joint analysis of various wave arrivals could provide information about velocity variations within permafrost and accurately estimate the permafrost thickness. Through a sensitivity analysis, we estimate the feasibility of using different wave arrivals, such as reflections, refractions, multiples, and guided waves, to better characterize permafrost conditions. We find that guided waves are sensitive to velocity variations within permafrost and can detect the depth of the permafrost’s base under certain geologic conditions. Then, we combine the analysis of all seismic arrivals to derive a subsurface permafrost model for a seismic line collected in the Beaufort Sea. The joint analysis reveals the transitions in depth between ice-bonded and partially ice-bonded permafrost and offers a crude estimate of permafrost thickness.

Subsea permafrost organic carbon stocks are large and of dominantly low reactivity

Subsea permafrost carbon pools below the Arctic shelf seas are a major unknown in the global carbon cycle. We combine a numerical model of sedimentation and permafrost evolution with simplified carbon turnover to estimate accumulation and microbial decomposition of organic matter on the pan-Arctic shelf over the past four glacial cycles. We find that Arctic shelf permafrost is a globally important long-term carbon sink storing 2822 (1518–4982) Pg OC, double the amount stored in lowland permafrost. Although currently thawing, prior microbial decomposition and organic matter aging limit decomposition rates to less than 48 Tg OC/yr (25–85) constraining emissions due to thaw and suggesting that the large permafrost shelf carbon pool is largely insensitive to thaw. We identify an urgent need to reduce uncertainty in rates of microbial decomposition of organic matter in cold and saline subaquatic environments. Large emissions of methane more likely derive from older and deeper sources than from organic matter in thawing permafrost.

Time-Domain Electromagnetics for Subsea Permafrost Mapping in the Arctic: The Synthetic Response Analyses and Uncertainty Estimates from Numerical Modelling Data

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.

Subsea Methane Hydrates: Origin and Monitoring the Impacts of Global Warming

The East Siberian Arctic shelf is the area where the largest natural gas reserves are concentrated. The formation of permafrost of the Arctic shelf during the Ice Age contributed to the emergence of a zone of stable existence of gas hydrates in the sedimentary layer, and subsequent flooding of the shelf led to its gradual degradation, the thawing of gas hydrates and the subsequent emissions of methane into the atmosphere. In the first part of the paper, we use mathematical modeling to study the processes of the formation of subsea permafrost on the Arctic shelf considering changes in the sea levels over the past 200 thousand years. Numerical simulations show the influence of climate warming up to 2200 on the degradation of subsea permafrost and the violation of the conditions for the stable existence of methane hydrates in bottom sediments using the example of the East Siberian shelf. The second part of the paper proposes a method for seismic monitoring of the state of gas hydrates based on a solution of multi-parameter inverse seismic problems. In particular, the degree of attenuation of seismic energy is one of the objective parameters for assessing the consolidation of gas hydrates: the closer they are to the beginning of decomposition, the higher the attenuation and, hence, the lower the quality factor. In this publication, we do not solve a multi-parameter inverse seismic problem for a real geological object. This would be impossible due to the lack of necessary data. Instead, we focus on substantiating the possibility of correct solutions for the problem of the reconstruction of the absorption and velocities for a viscoelastic medium in relation to the problem of monitoring the state of gas hydrate deposits. As noted in a range of publications, the thawing of gas hydrates leads to an increase in the fluid saturation of the geological medium followed by an increase in the absorption of seismic energy—that is, a decrease in the quality factor. Thus, the methods of seismic monitoring of the state of gas hydrates to predict the possibility of developing dangerous scenarios should be based on solving a multi-parameter inverse seismic problem. This publication is devoted to the presentation of this approach.

Exploring the capabilities of electrical resistivity tomography to study subsea permafrost

Abstract. Sea level rise and coastal erosion have inundated large areas of Arctic permafrost. Submergence by warm and saline waters increases the rate of inundated permafrost thaw compared to sub-aerial thawing on land. Studying the contact between the unfrozen and frozen sediments below the seabed, also known as the ice-bearing permafrost table (IBPT), provides valuable information to understand the evolution of sub-aquatic permafrost, which is key to improving and understanding coastal erosion prediction models and potential greenhouse gas emissions. In this study, we use data from 2D electrical resistivity tomography (ERT) collected in the nearshore coastal zone of two Arctic regions that differ in their environmental conditions (e.g., seawater depth and resistivity) to image and study the subsea permafrost. The inversion of 2D ERT data sets is commonly performed using deterministic approaches that favor smoothed solutions, which are typically interpreted using a user-specified resistivity threshold to identify the IBPT position. In contrast, to target the IBPT position directly during inversion, we use a layer-based model parameterization and a global optimization approach to invert our ERT data. This approach results in ensembles of layered 2D model solutions, which we use to identify the IBPT and estimate the resistivity of the unfrozen and frozen sediments, including estimates of uncertainties. Additionally, we globally invert 1D synthetic resistivity data and perform sensitivity analyses to study, in a simpler way, the correlations and influences of our model parameters. The set of methods provided in this study may help to further exploit ERT data collected in such permafrost environments as well as for the design of future field experiments.

Detecting subsea permafrost layers on marine seismic data: An appraisal from forward modelling

ABSTRACTDetecting the top and base subsea permafrost from 2D seismic reflection data in shallow marine settings is a non‐trivial task due to the occurrence of strong free surface multiples. The potential to accurately detect permafrost layers on conventional 2D seismic reflection data is assessed through viscoelastic modelling. Reflection imaging of permafrost layers is examined through the evaluation of specific characteristics of the subsurface, acquisition parameters and their impact. Results show that limitations are related to the principles of the method, the intrinsic nature of the permafrost layers, and the acquisition geometry. The biggest challenge is the occurrence of free surface multiples that overprint the base permafrost reflection, with the worst‐case scenario the case of a thin layer of ice‐bonded sand. Wedge models suggest that if the base permafrost is dipping, it would intersect internal and free surface multiples of the seafloor and the top permafrost and be detected. Also, the amplitude ratio of the base permafrost reflection and the multiples decreases with the increasing thickness of permafrost. Therefore, the crosscutting relationship between the reflection at base permafrost reflection and the multiples might not be enough to detect the base permafrost for thicker permafrost layers. Finally, the experiment results show that, for partially ice‐bonded layers, the attenuation combined with the low reflectivity of the basal interface limits the likelihood to resolve the base permafrost, especially for thick permafrost layers.

Climate warming-driven changes in the flux of dissolved organic matter and its effects on bacterial communities in the Arctic Ocean: A review

The warming of the Arctic Ocean impacts the dissolved organic matter (DOM) imports into the Arctic region, which affects the local bacterial communities. This review addressed the current status of DOM inputs and their potential influences on bacteria data (e.g., population, production, and metabolic activity of bacteria), as well as the projected changes of DOM inputs and bacterial communities as a result of climate warming. Microbial communities are likely affected by the warming climate and the transport of DOM to the Arctic Ocean. Imported DOM can alter Arctic bacterial abundance, cell size, metabolism, and composition. DOM fluxes from Arctic River runoff and adjacent oceans have been enhanced, with warming increasing the contribution of many emerging DOM sources, such as phytoplankton production, melted sea ice, thawed permafrost soil, thawed subsea permafrost, melted glaciers/ice sheets, atmospheric deposition, groundwater discharge, and sediment efflux. Imported DOM contains both allochthonous and autochthonous components; a large quantity of labile DOM comes from emerging sources. As a result, the Arctic sea water DOM composition is transformed to include a wider range of various organic constituents such as carbohydrates (i.e., glucose), proteinaceous compounds (i.e., amino acid and protein-like components) and those with terrigenous origins (i.e., humic-like components). Changes to DOM imports can alter Arctic bacterial abundance, cell size, metabolism, and composition. Under current global warming projections, increased inflow of DOM and more diverse DOM composition would eventually lead to enhanced CO2 emissions and frequent emergence of replacement bacterial communities in the Arctic Ocean. Understanding the changes in DOM fluxes and responses of bacteria in the Arctic broadens our current knowledge of the Arctic Ocean’s responses to global warming.

Atmospheric CO2 and CH4 Fluctuations over the Continent-Sea Interface in the Yenisei River Sector of the Kara Sea

Observations of the atmospheric sources and sinks of carbon dioxide (CO2) and methane (CH4) in the pan-Arctic domain are extremely scarce, limiting our knowledge of carbon turnover in this climatically sensitive environment and the fate of the enormous carbon reservoirs conserved in the permafrost. Especially critical are the gaps in the high latitudes of Siberia, covered by the vast permafrost underlain tundra, where only several atmospheric monitoring sites are operational. This paper presents the first two years (September 2018–January 2021) of accurate continuous observations of atmospheric CO2 and CH4 dry mole fractions at the recently deployed tower-based measurement station “DIAMIS” (73.5068° N, 80.5198° E) located on the southwestern coast of the Taimyr Peninsula, Siberia, at the Gulf of the Yenisei River that opens to the Kara Sea (Arctic Ocean). In this paper, we summarized the scientific rationale of the site, examined the seasonal footprint of the station with an analysis of terrestrial vegetation and maritime sector contributing to the captured atmospheric signal, and illustrated temporal patterns of CO2 and CH4 for the daytime mixed atmospheric layer over the continent–sea interface. Along with the temporal variations reflecting a signal caused pan-Arctic and not very much influenced by the local processes, we analyzed the spatiotemporal distribution of the synoptic anomalies representing the atmospheric signatures of regional sources and sinks of CO2 and CH4 for the studied high-arctic Siberian domain of ~625 thousand km2, with nearly equal capturing the land surface (54%) and the ocean (46%) throughout the year. Both for CO2 and CH4, we have observed a sea–continent declining trend, presuming a larger depletion of trace gases in the maritime air masses compared to the continental domain. So far, over the Kara Sea, we have not detected any prominent signals of CH4 that might have indicated processes of subsea permafrost degradation and occurrence of cold seeps–still mainly observed in the eastern Arctic Seas—The Laptev Sea and the East-Siberian Sea.