The objective of the present study is to analyze the interactions between a methane hydrates stability zone and the ocean temperature variations and to define the hydrate sensitivity to the contemporary warming in the Arctic Ocean. To obtain the spatial–temporary variability of the ocean bottom temperature we employ the ICMMG regional Arctic-North Atlantic ocean model that has been developed in the Institute of Computational Mathematics and Mathematical Geophysics. With the ice-ocean model the Arctic bottom water temperatures were analyzed. The resulting warming ocean bottom water is spatially inhomogeneous, with a strong impact by the Atlantic inflow on shallow regions of 200-500 m depth. Results of the mathematical modeling of the dynamics of methane hydrate stability zone in the Arctic Ocean sediment are reported. We find that the reduction of the methane hydrate stability zone occurs in the Arctic Ocean between 250 and 400 m water depths within the upper 100 m of sediment in the area influenced by the Atlantic inflow. We have identified the areas of the Arctic Ocean where an increase in methane release is probable to occur at the present time.
The vertical distribution of the hydrological characteristics in the upper ocean layer is mostly formed under the influence of turbulent and convective mixing, which are not resolved in the system of equations for large-scale ocean. Therefore it is necessary to include additional parameterizations of these processes into the numerical models. In this paper we carry out a comparative analysis of the different vertical mixing parameterizations in simulations of climatic variability of the Arctic water and sea ice circulation. The 3D regional numerical model for the Arctic and North Atlantic developed in the ICMMG SB RAS (Institute of Computational Mathematics and Mathematical Geophysics of the Siberian Branch of the Russian Academy of Science) and package GOTM (General Ocean Turbulence Model<sup>1,2</sup>, http://www.gotm.net/) were used as the numerical instruments . NCEP/NCAR reanalysis data were used for determination of the surface fluxes related to ice and ocean. The next turbulence closure schemes were used for the vertical mixing parameterizations: 1) Integration scheme based on the Richardson criteria (RI); 2) Second–order scheme TKE with coefficients Canuto-A<sup>3</sup> (CANUTO); 3) First–order scheme TKE with coefficients Schumann and Gerz<sup>4</sup> (TKE-1); 4) Scheme KPP5 (KPP). In addition we investigated some important characteristics of the Arctic Ocean state including the intensity of Atlantic water inflow, ice cover state and fresh water content in Beaufort Sea.