海冰形变对北极冬季CICE模式冰厚模拟结果的影响

海冰数值模式是研究海冰动力热力状态参量及之间联系的有效途径。目前对冰厚数值模拟结果的分析远远少于对海冰范围/面积和密集度的研究,对冰速与海冰形变对冰厚分布影响的研究也尚欠缺。本文利用Los Alamos sea ice model（CICE）海冰模式模拟了1980−2018年的北极海冰变化,并使用遥感、同化冰厚数据进行比对验证,分析了模拟冰速和海冰形变对冰厚的影响,计算了冰速的散度和切变偏差对冰厚偏差的贡献。结果显示,CICE对北极70°N以北区域平均冰厚和冰速的年际变化模拟基本合理,但模拟的平均冰厚和冰速多年变化趋势均小于同化数据的变化率;模拟和观测冰厚的空间分布差异与冰速和形变率的偏差有密切联系,主要表现为波弗特海的正偏差和北极中央区至弗拉姆海峡的负偏差。泛北极区域散度和切变偏差在3月之前对冰厚偏差的贡献在13%～16%之间变化,3−4月则由16%跃变至27%。散度偏差主导了11月、12月波弗特海区域的冰厚正偏差,切变偏差主导了冬季加拿大群岛以北海域和穿极流区域的冰厚负偏差。     

渤海和北极海冰组构及晶体结构特征分析

渤海一年冰的组构和晶体结构表现出的冰的非连续生长速率使得冰层底部具有气泡含量不同的许多分层;静水生长的灰冰具有典型的粒状冰和柱状冰;融化-重冻结冰和重叠冰具有特定的晶体结构.位于北冰洋72°24.037'N,153°33.994'W,长2.2m的冰心样为三年冰,并且发现其内有一种新的晶体类型,定义为碎屑凝聚冰.由位于北冰洋74°58.614'N,160°31.830'W,长4.86m的冰心样的晶体结构分析发现它过去为冰脊,后经融化-冻结改造形成冰丘的证据.     

海冰热力模式的辐射参数化方案在渤海和波罗的海应用的比较

研究了海冰热力模式中的各种辐射参数化方案,对比了模式计算的太阳短波辐射、大气长波辐射以及海冰热力变化,并利用渤海和波罗的海观测资料进行比较和误差分析.冬季大部分时间太阳短波辐射对海冰热力过程的作用有限.简单计算方案一般满足海冰模式要求.误差主要受云和冰雪表面与大气之间的多重反射影响.长波辐射对表面热平衡和海冰质量变化起重要作用.长波辐射参数化方案的计算结果受环境因素影响.云量参数化有待进一步改进.海冰模式计算结果的精度与长波辐射计算精度有一致性.     

渤海海冰特征厚度分析

通过海冰生消机理和数值试验,讨论了渤海海冰特征厚度的存在条件;对不同厚度的海冰表面温度、冰面热量收支、冰面下热传导和太阳辐射透射量进行了对比分析,分析了渤海海冰向特征冰厚的动态演化过程;在不同气温、风速、相对湿度和海洋热通量等气象和水文条件下,对渤海特征冰厚进行了计算;讨论了海冰生消的动态平衡过程,分析了1997/1998年冬季辽东湾JZ20-2海域实测冰厚与特征冰厚的相互关系。对渤海特征冰厚分析将有助于渤海海冰数值模拟工作的开展和对不同重现期设计冰厚的推算。     

气候因子对渤海冰情影响的统计分析

运用相关分析、偏相关分析等统计分析方法研究了北极涛动、北极海冰和西伯利亚高压等气候因子对渤海海冰年际变化的影响。统计表明西伯利亚高压是影响渤海海冰年际变化最重要的因子。秋季北极海冰密集度、冬季北极涛动与渤海冰情显著相关,它们可能通过影响西伯利亚高压的强度,从而影响渤海海冰的严重程度。利用回归分析得到渤海冰情等级与西伯利亚高压、西太平洋副高之间的统计关系,回归结果较好的反映了渤海冰情年际变化。     

北极巴伦支海和拉普捷夫海海冰中细菌和微藻的垂直分布

通过计算机图像测定系统测定、计算获得海冰细菌的大小、面积和碳含量,研究了北极巴伦支海和拉普捷夫海的浮冰中细菌大小、丰度和生物量的垂直分布及其与海冰微藻生物量的关系.细菌在冰心中是非均匀分布的,至少有两种不同分布类型:第一类有单一的丰度峰值,在冰心底部或在冰心中部;第二类有两个丰度峰值.海冰中细菌的丰度为0.4×105～36.7×105个/cm3.不同类型海冰中的细菌大小变化极大,在多年冰中,最大的细菌与叶绿素a最大值出现在同一层,而在一年冰中,细菌大小几乎没有垂直变化.整柱冰心的细菌生物量变化为19.2～79.2mg/cm2,细菌与海冰微藻生物量之比为0.43～10.00.对固定冰比较和研究的结果发现,海冰冰心中细菌大小、丰度和生物量的垂直分布差别极大.据此,分析了对目前海冰研究采样方法的局限性,并提出了规范采样方法的设想.     

Geophysics of sea ice in the Baltic Sea: A review

With improved observation methods, increased winter navigation, and increased awareness of the climate and environmental changes, research on the Baltic Sea ice conditions has become increasingly active. Sea ice has been recognized as a sensitive indicator for changes in climate. Although the inter-annual variability in the ice conditions is large, a change towards milder ice winters has been detected from the time series of the maximum annual extent of sea ice and the length of the ice season. On the basis of the ice extent, the shift towards a warmer climate took place in the latter half of the 19th century. On the other hand, data on the ice thickness, which are mostly limited to the land-fast ice zone, basically do not show clear trends during the 20th century, except that during the last 20 years the thickness of land-fast ice has decreased. Due to difficulties in measuring the pack-ice thickness, the total mass of sea ice in the Baltic Sea is, however, still poorly known. The ice extent and length of the ice season depend on the indices of the Arctic Oscillation and North Atlantic Oscillation. Sea ice dynamics, thermodynamics, structure, and properties strongly interact with each other, as well as with the atmosphere and the sea. The surface conditions over the ice-covered Baltic Sea show high spatial variability, which cannot be described by two surface types (such as ice and open water) only. The variability is strongly reflected to the radiative and turbulent surface fluxes. The Baltic Sea has served as a testbed for several developments in the theory of sea ice dynamics. Experiences with advanced models have increased our understanding on sea ice dynamics, which depends on the ice thickness distribution, and in turn redistributes the ice thickness. During the latest decade, advance has been made in studies on sea ice structure, surface albedo, penetration of solar radiation, sub-surface melting, and formation of superimposed ice and snow ice. A high vertical resolution has been found as a prerequisite to successfully model thermodynamic processes during the spring melt period. A few observations have demonstrated how the river discharge and ice melt affect the stratification of the oceanic boundary layer below the ice and the oceanic heat flux to the ice bottom. In general, process studies on ice–ocean interaction have been rare. In the future, increasingly multidisciplinary studies are needed with close links between sea ice physics, geochemistry and biology.     

基于高光谱遥感的渤海海冰厚度半经验模型

Sea ice thickness is one of the most important input parameters for the prevention and mitigation of sea ice disasters and the prediction of local sea environments and climates. Estimating the sea ice thickness is currently the most important issue in the study of sea ice remote sensing. With the Bohai Sea as the study area, a semiempirical model of the sea ice thickness(SEMSIT) that can be used to estimate the thickness of first-year ice based on existing water depth estimation models and hyperspectral remote sensing data according to an optical radiative transfer process in sea ice is proposed. In the model, the absorption and scattering properties of sea ice in different bands(spectral dimension information) are utilized. An integrated attenuation coefficient at the pixel level is estimated using the height of the reflectance peak at 1 088 nm. In addition, the surface reflectance of sea ice at the pixel level is estimated using the 1 550–1 750 nm band reflectance. The model is used to estimate the sea ice thickness with Hyperion images. The first validation results suggest that the proposed model and parameterization scheme can effectively reduce the estimation error associated with the sea ice thickness that is caused by temporal and spatial heterogeneities in the integrated attenuation coefficient and sea ice surface. A practical semi-empirical model and parameterization scheme that may be feasible for the sea ice thickness estimation using hyperspectral remote sensing data are potentially provided.     

CICE5.0与BCC_CSM2.0模式的耦合及对北极海冰的模拟评估

本文将美国Los Alamos国家实验室发展的最新海冰模式CICE5.0引入国家气候中心气候系统模式BCC_CSM2.0,替代原有的海冰模式SIS,形成一个新的耦合模式。在此基础上,评估新耦合模式对1985-2009年北极海冰的模拟性能,检验引入CICE5.0后对耦合模式中北极海冰、海洋和大气模拟结果的改进。结果表明,引入CICE5.0后,模式能较好地模拟出北极海冰的空间分布、季节以及年际变化特征。相比于旧版本耦合模式,新耦合模式模拟的北极多年冰增多、一年冰减少,同时,海冰增厚、海冰流速减慢,模拟效果得到显著改进,对波弗特涡流模拟的改善尤为明显。进一步分析发现,相比于SIS,CICE5.0对北极海冰特别是海冰厚度模拟性能的提升,在耦合进入BCC_CSM2.0后,会触发冰-温的正反馈机制,改进了模式对海平面气压场、表层气温和海表温度的模拟,由此进一步提高了模式对北极海冰的模拟能力。     

北极海冰密集度预报对大气强迫敏感性的个例研究

A regional Arctic configuration of the Massachusetts Institute of Technology general circulation model （MIT-gcm） is used as the coupled ice-ocean model for forecasting sea ice conditions in the Arctic Ocean at the Na-tional Marine Environmental Forecasting Center of China （NMEFC）, and the numerical weather prediction from the National Center for Environmental Prediction Global Forecast System （NCEP GFS） is used as the atmospheric forcing. To improve the sea ice forecasting, a recently developed Polar Weather Research and Forecasting model （Polar WRF） model prediction is also tested as the atmospheric forcing. Their forecasting performances are evaluated with two different satellite-derived sea ice concentration products as initializa-tions： （1） the Special Sensor Microwave Imager/Sounder （SSMIS） and （2） the Advanced Microwave Scanning Radiometer for EOS （AMSR-E）. Three synoptic cases, which represent the typical atmospheric circulations over the Arctic Ocean in summer 2010, are selected to carry out the Arctic sea ice numerical forecasting experiments. The evaluations suggest that the forecasts of sea ice concentrations using the Polar WRF atmo-spheric forcing show some improvements as compared with that of the NCEP GFS.     

海冰动力过程的改进离散元模型及在渤海的应用

海冰的断裂、重叠和堆积等离散分布特性广泛地存在于极区和副极区的不同海域,并对海冰的生消、运移过程有着重要影响。针对海冰在不同尺度下的离散分布特点,发展海冰动力过程的离散元方法有助于完善海冰数值模式,提高海冰数值模拟的计算精度。为此,本文针对海冰生消运移过程中的非连续分布和形变特性,发展了适用于海冰动力过程的改进离散元模型(MDEM)。不同于传统离散元方法,该模型将海冰离散为具有一定厚度、尺寸和密集度的圆盘单元。海冰单元设为诸多浮冰块的集合体,其在运移和相互接触碰撞过程中,依照质量守恒发生单元尺寸、密集度和厚度的相应变化。基于海冰离散性和流变性的特点,该模型采用黏弹性接触本构模型计算单元间的作用力,并依据Mohr-Coulomb准则计算海冰法向作用下的塑性变形及切向摩擦力。为验证该模型的可靠性,本文对海冰在规则水域内的运移和堆积过程进行了分析,离散元计算结果与解析值相一致;此外,对旋转风场下海冰漂移规律的模拟进一步验证了本文方法的精确性。在此基础上,对渤海辽东湾的海冰动力过程进行了48h数值分析,计算结果与卫星遥感资料和油气作业区的海冰现场监测数据吻合良好。在下一步工作中将考虑海冰离散元模拟中的热力因素影响,发展具有冻结、断裂效应的海冰离散元模型,更精确地模拟海冰动力-热力耦合作用下的生消和运移过程。     

北极楚科奇海海冰面积多年变化的研究

北极气候系统正在发生显著变化,其中,海冰面积和厚度的减小是其最主要的特征.楚科奇海是海冰面积变化最有代表性的区域.文章利用积累了9a的高分辨率海冰分布数据研究海冰面积的多年变化特征.结果表明,各年的冰情有显著的季节内变化,海冰面积距平曲线体现了不同时期海冰面积变化的动态过程.在1997~2005年间,楚科奇海海冰面积经历了轻(1997年)-重(2000~2001年)-轻(2002~2005年)的变化过程.9a的数据总体上体现了海冰面积减小的趋势,2005年的冰情呈现了历史新低.每年融冰期的长短与冰情轻重有密切的关系,冰轻年份融冰开始时间早,冻结结束时间晚.各年海冰面积最小值发生在9月下旬至10月初,各个年份海冰最小面积差别很大.有的年份只有4%,而重冰年可以大于50%.文章采用4个重要参数表达海冰多年变化.其中海冰面积指数反映了当年总体平均的海冰面积距平;海冰最小面积反映了融冰期海冰的极限情况;上一个冬季的气温积温也与翌年海冰面积有良好的关联;分析了风场对海冰的影响,表明风场在融冰期能够在短时间内改变海冰的覆盖面积.     

渤海及黄海北部冰情长期变化趋势分析

统计分析上世纪50年代—2010年渤海及黄海北部海冰资料,对其年代际变化特征进行分析。上世纪50年代—90年代冰情总体呈缓解的趋势,2000年以来冰情略有加重。研究发现太阳活动与渤海及黄海北部冰情变化关系密切,太阳活动可能是渤海及黄海北部冰情长期变化重要影响因素。如果太阳黑子的周期长度比上一个周期长,那么周期内冰情较上一个周期严重。反之亦然。     

Sediment transport by sea ice in the Chukchi and Beaufort Seas: Increasing importance due to changing ice conditions?

Sediment-laden sea ice is widespread over the shallow, wide Siberian Arctic shelves, with off-shelf export from the Laptev and East Siberian Seas contributing substantially to the Arctic Ocean's sediment budget. By contrast, the North American shelves, owing to their narrow width and greater water depths, have not been deemed as important for basin-wide sediment transport by sea ice. Observations over the Chukchi and Beaufort shelves in 2001/02 revealed the widespread occurrence of sediment-laden ice over an area of more than 100,000 km² between 68 and 74°N and 155 and 170°W. Ice stratigraphic studies indicate that sediment inclusions were associated with entrainment of frazil ice into deformed, multiple layers of rafted nilas, indicative of a flaw-lead environment adjacent to the landfast ice of the Chukchi and Beaufort Seas. This is corroborated by buoy trajectories and satellite imagery indicating entrainment in a coastal polynya in the eastern Chukchi Sea in February of 2002 as well as formation of sediment-laden ice along the Beaufort Sea coast as far eastward as the Mackenzie shelf. Moored upward-looking sonar on the Mackenzie shelf provides further insight into the ice growth and deformation regime governing sediment entrainment. Analysis of Radarsat Synthetic Aperture (SAR) imagery in conjunction with bathymetric data help constrain the water depth of sediment resuspension and subsequent ice entrainment (>20 m for the Chukchi Sea). Sediment loads averaged at 128 t km^–2, with sediment occurring in layers of roughly 0.5 m thickness, mostly in the lower ice layers. The total amount of sediment transported by sea ice (mostly out of the narrow zone between the landfast ice edge and waters too deep for resuspension and entrainment) is at minimum 4×10⁶ t in the sampling area and is estimated at 5–8×10⁶ t over the entire Chukchi and Beaufort shelves in 2001/02, representing a significant term in the sediment budget of the western Arctic Ocean. Recent changes in the Chukchi and Beaufort Sea ice regimes (reduced summer minimum ice extent, ice thinning, reduction in multi-year ice extent, altered drift paths and mid-winter landfast ice break-out events) have likely resulted in an increase of sediment-laden ice in the area. Apart from contributing substantially to along- and across-shelf particulate flow, an increase in the amount of dirty ice significantly impacts (sub-)ice algal production and may enhance the dispersal of pollutants.     

Improvement of short-term sea ice forecast in the southern Okhotsk Sea

In this study, a numerical model of 7-day forecast of sea ice produced by the Japan Meteorological Agency was improved by the following approaches. First, a new ice dynamic model was introduced: the distributed mass/discrete floe model. The model takes account of discrete characteristics of ice floes and well simulates the ice edge location at low computational cost. Secondly, the grid size was reduced to 5 × 5 km for the future high resolution forecasts. Next, the sea surface current data was examined because it significantly influences sea ice movement. We applied two new datasets of HINO and Okhotsk Ocean General Circulation Model (Okhotsk OGCM), which are estimated by numerical simulations, for the 7-day forecast of sea ice. Ice southward speed in January and the whorl formations in February and March were well reproduced with Okhotsk OGCM datasets. Finally, the ocean heat flux at the ice-ocean interface was refined. As a result, we achieved an ice edge error reduction from 30.8 km to 23.5 km.     

Properties of sea ice and overlying snow in the Southern Sea of Okhotsk

The general properties of sea ice and overlying snow in the southern Sea of Okhotsk were examined during early February of 2003 to 2005 with the P/V “Soya”. Thin section analysis of crystal structure revealed that frazil ice (48% of total core length) was more prevalent than columnar ice (39%) and that stratigraphic layering was prominent with a mean layer thickness of 12 cm, indicating that dynamic processes are essential to ice growth. The mean thickness of ice blocks and visual observations suggest that ridging dominates the deformation process above thicknesses of 30 to 40 cm. As for snow, it was found that faceted crystals and depth hoar are dominant (78%), as which is also common in the Antarctic sea ice, and is indicative of the strong vertical temperature gradients within the snow. Stable isotope measurements (δ¹⁸O) indicate that snow ice occupies 9% of total core length and that the mass fraction of meteoric ice accounts for 1 to 2% of total ice volume, which is lower than the Antarctic sea ice. Associated with this, the effective fractionation coefficient during the freezing of seawater was also derived. Snow ice was characterized by lower density, higher salinity, and nearly twice the gas content of ice of seawater origin. In addition, it is shown that the surface brine volume fraction and freeboard are well correlated with ice thickness, indicating some promise for remote sensing approaches to the estimation of ice thickness.     

The effects of resolving the Canadian Arctic Archipelago in a finite element sea ice model

Though narrow straits may have a strong influence on the large-scale sea ice mass balance, they are often crudely represented in coarse resolution sea ice models. Unstructured meshes, with their natural ability to fit boundaries and locally increase the mesh resolution, propose an alternative framework to capture the complex oceanic areas formed by coasts and islands. In this paper, we develop a finite element sea ice model to investigate the sensitivity of the Arctic sea ice cover features to the resolution of the narrow straits constituting the Canadian Arctic Archipelago. The model is a two-level dynamic-thermodynamic sea ice model, including a viscous-plastic rheology. It is run over 1979–2005, forced by daily NCEP/NCAR reanalysis data. Confronting qualitatively numerical experiments with observations shows a good agreement with satellite and buoys measurements. Due to its simple representation of the oceanic interactions, the model overestimates the sea ice extent during winter in the southernmost parts of the Arctic, while the Baffin Bay and Kara Sea remain ice-covered during summer. In order to isolate the benefits from resolving the Canadian Arctic Archipelago, a numerical experiment is performed where we artificially close the archipelago. Focusing on the large-scale sea ice thickness pattern, no significant change is found in our model, except in the close surroundings of the archipelago. However, the local and short-term influences of the ice exchanges are nonnegligible. In particular, we show that the ice volume associated to the Canadian Arctic Archipelago represents 10% of the Northern Hemisphere sea ice volume and that the annual mean ice export towards Baffin Bay amounts to 125 km³ yr⁻¹, which may play an important role on the convective overturning in the Labrador Sea.     

基于MODIS热红外数据的渤海海冰厚度反演

Level ice thickness distribution pattern in the Bohai Sea in the winter of 2009–2010 was investigated in this paper using MODIS night-time thermal infrared imagery.The cloud cover in the imagery was masked out manually.Level ice thickness was calculated using MODIS ice surface temperature and an ice surface heat balance equation.Weather forcing data was from the European Centre for Medium-Range Weather Forecasts(ECMWF) analyses.The retrieved ice thickness agreed reasonable well with in situ observations from two off-shore oil platforms.The overall bias and the root mean square error of the MODIS ice thickness are –1.4 cm and 3.9 cm,respectively.The MODIS results under cold conditions(air temperature –10°C) also agree with the estimated ice growth from Lebedev and Zubov models.The MODIS ice thickness is sensitive to the changes of the sea ice and air temperature,in particular when the sea ice is relatively thin.It is less sensitive to the wind speed.Our method is feasible for the Bohai Sea operational ice thickness analyses during cold freezing seasons.     

On the future navigability of Arctic sea routes: High-resolution projections of the Arctic Ocean and sea ice

The rapid Arctic summer sea ice reduction in the last decade has lead to debates in the maritime industries on the possibility of an increase in cargo transportation in the region. Average sailing times on the North Sea Route along the Siberian Coast have fallen from 20 days in the 1990s to 11 days in 2012–2013, attributed to easing sea ice conditions along the Siberian coast. However, the economic risk of exploiting the Arctic shipping routes is substantial. Here a detailed high-resolution projection of ocean and sea ice to the end of the 21st century forced with the RCP8.5 IPCC emission scenario is used to examine navigability of the Arctic sea routes. In summer, opening of large areas of the Arctic Ocean previously covered by pack ice to the wind and surface waves leads to Arctic pack ice cover evolving into the Marginal Ice Zone. The emerging state of the Arctic Ocean features more fragmented thinner sea ice, stronger winds, ocean currents and waves. By the mid 21st century, summer season sailing times along the route via the North Pole are estimated to be 13–17 days, which could make this route as fast as the North Sea Route.     

两种盐度参数化方案在海冰模式中的应用比较

本文详细介绍了SIS海冰模式中引进两种盐度参数化方案即等盐度方案和盐度廓线方案对海冰模拟所存在的差异。利用盐度廓线方案导出的表征盐度与海冰温度间关系的方程比等盐度方案多出一项,将定义为盐度差异项。盐度差异项对海冰厚度的热力作用表现为：在海冰厚度增长季节（11月到次年5月）,盐度差异项通过升高海冰内部温度,抑制海冰增长;在消融的第一阶段（6．8月）,盐度差异项通过升高海冰内部温度加快海冰消融;在消融的第二阶段（9．10月）,盐度差异项通过降低海冰内部的温度抑制海冰消融。但尺度分析表明,盐度差异项要比方程中队海冰温度作用最大项小1．2个量级,如果采用一级近似,可以略去盐度差异项,因此盐度差异项对海冰增长和消融影响很小。同时利用冰洋耦合模式（ModularOceanModel,MOM4）,分别采用两种盐度参数化方案模拟北极海冰厚度和海冰密集度的季节性变化,模拟结果也表明两种方案模拟得到的海冰厚度和海冰密集度的季节性变化相差甚小。