北冰洋浮冰站近地层参数的观测估算

利用2008年8月20～27日我国第3次北极考察队在85°N附近设立的冰站上进行的湍流通量、辐射观测所获取的相关资料,采用涡动相关法对夏季北冰洋浮冰下垫面的近地层参数进行了估算.结果显示,观测期间浮冰区冰雪面的平均感热、潜热和净辐射通量分别是0.2 W/m2,1.2 W/m2和9.9 W/m2,表明下垫面获得的大部分热...     

基于走航观测数据的ERA5南海北部夏季短波辐射收支评估

海表短波辐射收支是海–气界面能量交换的重要物理过程。本研究利用2019年南海北部夏季科考航次的走航观测数据,评估了ERA5再分析数据的海表短波辐射通量收支。结果表明,ERA5的向下短波辐射相比观测偏小,11时和15时(北京时间)的偏差最大,可达-100 W/m²。与此同时,ERA5的海表反照率整体偏低,其中高太阳高度角时段偏差较小,约为-0.03,低太阳高度角时段偏差较大,约为-0.15。向下短波辐射和反照率的偏差共同造成ERA5白天平均海表净短波辐射通量比观测偏小约25.4 W/m²;其中,反照率低估抵消了约50%向下短波辐射偏差的贡献。研究表明,在不同大气透射率情况下,ERA5的海表辐射收支偏差存在不同表现。ERA5海表反照率的低估可能与其采用的参数化方案在南海北部的适用性不足有关。基于观测本研究也给出了一个简单的参数优化方案。     

厦门近岸210Po和210Pb的大气沉降通量

为揭示中国东南沿海地区210Po和210Pb的大气沉降时空变化特征，探讨该地区气溶胶的停留时间，于2013年1月至2014年12月对厦门地区210Po和210Pb的大气沉降通量进行了时间序列研究。结果发现，210Po和210Pb的平均日沉降通量分别为（65.38±4.79） mBq/（m2·d）（n=54）和（0.78±0.09） Bq/（m2·d）（n=54），表现出明显的周年变化。东北季风期间，210Po和210Pb的沉降通量较高，而西南季风期间其通量较低。2013年和2014年，210Po的年沉降通量分别为19.29 Bq/（m2·a）和9.25 Bq/（m2·a），210Pb的年沉降通量分别为159.2 Bq/（m2·a）和189.6 Bq/（m2·a）。两核素的年沉降通量表现出不同程度的年际差异。210Po与210Pb沉降通量之间存在显著的线性正相关关系，揭示了大气中210Po和210Pb具有相同的迁出机制，降雨和大气中核素含量是影响210Po和210Pb沉降通量的主要因素。该研究结果可以为探求台湾海峡海水中210Po与210Pb的收支平衡提供大气来源项。     

北极夏季海冰反照率的观测和数值模拟试验

在中国第3次北极科学考察浮冰站开展了积雪/海冰反照率观测.本文对观测结果进行了分析,并结合一维高分辨雪/冰模式(HIGHTSI)对3个常用的反照率参数化方案在天气尺度的表现进行了评估.观测期间测站反照率变化范围0.75～0.85,其天气尺度变化同天气和表面冰、雪状况紧密相关,降雪和吹雪过程可改变表面积雪厚度及水平分布,...     

雪热传导系数与穿过海冰的热通量研究

雪热传导系数是海冰质量平衡过程中的重要物理参数，决定了穿透海冰的热传导通量。北冰洋海冰质量平衡浮标观测获得多年冰上冬季温度链剖面可以明显地区分冰雪界面。本文考虑到冰雪界面处温度随时间变化，再根据冰雪界面热传导通量连续假定，提出了新的雪热传导系数计算方法。受不同环境因素影响，多年冰上各个浮标的雪热传导系数在0.23~0.41 W/（m·K）之间，均值为（0.32±0.08） W/（m·K）。北冰洋多年冰上冬季穿过海冰的热传导通量最大发生在11月至翌年3月，约14~16 W/m2。结冰季节，来自海冰自身降温的热量对穿过海冰向大气传输的热量贡献逐月减少，从9月100%减小到12月的35%，翌年的1月至3月稳定在10%左右。夏季，短波辐射通能量通过热传导自上而下加热海冰，海冰上层温度高于下层，热量传播方向与冬季反向，往海冰内部传递。直到9月短波辐射完全消失，气温下降，热量再次转变为自下往上传递。从冰底热传导来看，夏季出现海冰向冰水界面传递热量现象。由于雪较好的绝热性，冰上覆雪极大地削弱了海冰上层热传导通量，从而减缓了秋冬季节的结冰速度。尽管受雪厚影响，多年冰上层热传导通量与气温依旧具有很好的线性关系，气温每降低1℃，热传导通量增加约0.59 W/m2。     

塔基海-气界面动量通量观测：向岸风和离岸风情形比较

本文使用塔基直接观测法研究海洋大气边界层中的海-气界面动量通量。首先，我们收集数据并和前人观测结果比对，其比对结果符合一致。其次，在低风速至中等风速条件下，我们发现海-气界面动量通量的交换系数（又称拖曳系数）对于向岸风和离岸风两种情形有所差异。为此，我们使用一个考虑表面波的参数化方案解释海-气界面动量通量对于表面波的依赖关系。这些结果一方面证实表面波对于海-气界面动量通量的影响，另一方面验证一个考虑表面波参数化方案的有效性。     

琼州海峡冬末春初潮余流场特征

1995年2月26日至3月6日,在琼州海峡的新海一四塘断面上,大小潮期间进行15条船同步测流.该断面上涨潮流速普遍大于落潮流速,实测最大涨潮流速为172cm/s,最大落潮流速为142cm/s;大潮期间北部6个站全层平均余流速度为18.4cm/s.海峡中间3个站0～20m层平均流速为18.3cm/s,南部6个站平均余流速度为10.2m/s.中间和北部诸站余流方向指向W一S范围,南部6个站因受地形影响,流向指向NE.用数值计算方法,再现了大小潮期间琼州海峡整个潮流场.大潮期间,通过西断面(灯楼角──玉苞角)落潮流总通量为4.73×1010m3,涨潮流总通量为4.29×1010m3;通过东断面(东营──龙塘镇)落潮流总通量为5.22×1010m3,涨潮流总通量为4.90×1010m3,其净通量,西、东断面分别为0.43×1010和0.32×1010m3,其方向指向西.     

极区通量观测系统及其在国际极地年(IPY)全球协同观测中的应用

本文对极区通量观测系统作了介绍,在国际极地年(IPY)全球协同观测中,极区通量观测系统在南极中山站进行了连续14个月的观测。结果表明,中山站年净辐射通量为12.9 W/m2。感热通量夏半年(10～2月)为正值,冬半年(3～9月)为负值,年平均1.9 W/m2。潜热通量全年都为正值,年平均11.2 W/m2。总体而言,地表通过净辐射获得热能,又通过感热和潜热方式向大气输送。观测得到的CO2通量全为负值,年平均为-0.031 mg/m2,表明南极中山站是CO2汇。     

南黄海海气热通量观测及其与OAflux数据集比较研究

2007年在南黄海进行了3个航次的热通量观测,包括长、短波辐射,近海表空气温度、湿度,风速,海表皮温等观测数据。依据计算的冬季、春季、秋季三个航次的海气热通量分析了热通量不同季节特征,南黄海海域冬季、春季和秋季平均潜热通量分别为80.7W/m2,5.6W/m2和142.1W/m2,感热通量分别为32.0W/m2,-12.5W/m2和18.9W/m2(海洋向大气传递为正)。将国际较为通用的OAflux数据集与3个季节观测数据做了逐点的比对,作为对OAflux数据集在南黄海海域的评估,结果显示:OAflux数据集热通量结果与观测数据在2006—2007年冬季最为接近,感热和潜热通量均方差是15.3W/m2和21.4W/m2。春季的潜热通量存在明显偏差,均方差为28.4W/m2。秋季的感热和潜热通量均存在显著偏差,均方差分别为20.5W/m2和57.5W/m2。导致春季偏差的主要原因是OAflux数据集和现场观测的近海表空气湿度差异,而秋季偏差则应主要归因于海表温度的偏差。     

基于Icepack海冰柱模式的融池反照率模拟研究

基于Icepack一维海冰柱模式,以2014年中国第6次北极科学考察长期冰站ICE06的3个融池的辐射参量和气象参量的连续观测作为大气强迫数据,对融池反照率及相关参量进行了模拟。本文引入观测的融池深度及海冰厚度作为初始条件,通过考虑融池覆盖率的作用,改进了平整冰融池参数化方案中海冰干舷的计算,修正了冰上可允许的最大融池深度,成功实现了对融池参数变化的模拟;同时,还修正了入射辐射分量比例系数与对应反照率分量权重系数不一致的问题。标准试验中,模拟的3个融池的反照率与观测结果之间的平均误差分别为0.01、0.05和0.13;入射辐射比例的敏感性试验结果表明,当可见光辐射比例增大8%时,融池反照率的模拟结果增大了6%～8%;融池表面再冻结试验的结果显示,当再冻结冰层厚度小于2 cm时,模拟冰面反照率的增加不足0.006,由此引起的表面能量收支减少了约1.1 W/m2。本文研究指出,准确的入射辐射比例对于改善北极海冰反照率模拟是必要的;并指出目前模式仍存在融池表面再冻结参数化、热收支计算、表面吹雪效应等有待解决的问题。     

北冰洋浮冰区湍流通量观测试验及参数化研究

利用2008年8月21～29日我国第3次北极考察期间在北冰洋海区(84°27′N,143°37′W～85°13′N,147°20′W)冰站观测的湍流资料及相关资料,对海冰近地层湍流通量及其特征参数进行了研究.结果表明:观测期间浮冰近地层始终存在逆温和逆湿层.这与我们以前(1999年在75°N和2003年在78°N)的观...     

2010年马卡罗夫海盆冰站观测期间垂向热通量时间变化研究

Based on hydrographic data obtained at an ice camp deployed in the Makarov Basin by the 4th Chinese Arctic Research Expedition in August of 2010, temporal variability of vertical heat flux in the upper ocean of the Makarov Basin is investigated together with its impacts on sea ice melt and evolution of heat content in the remnant of winter mixed layer(r WML). The upper ocean of the Makarov Basin under sea ice is vertically stratified. Oceanic heat flux from mixed layer(ML) to ice evolves in three stages as a response to air temperature changes, fluctuating from 12.4 W/m2 to the maximum 43.6 W/m2. The heat transferred upward from ML can support(0.7±0.3) cm/d ice melt rate on average, and daily variability of melt rate agrees well with the observed results. Downward heat flux from ML across the base of ML is much less, only 0.87 W/m2, due to enhanced stratification in the seasonal halocline under ML caused by sea ice melt, indicating that increasing solar heat entering summer ML is mainly used to melt sea ice, with a small proportion transferred downward and stored in the r WML. Heat flux from ML into r WML changes in two phases caused by abrupt air cooling with a day lag. Meanwhile, upward heat flux from Atlantic water(AW) across the base of r WML, even though obstructed by the cold halocline layer(CHL), reaches0.18 W/m2 on average with no obvious changing pattern and is also trapped by the r WML. Upward heat flux from deep AW is higher than generally supposed value near 0, as the existence of r WML enlarges the temperature gradient between surface water and CHL. Acting as a reservoir of heat transferred from both ML and AW, the increasing heat content of r WML can delay the onset of sea ice freezing.     

东南极中山站陆缘固定冰2011/2012年度的时空变化

Annual observations of first-year ice(FYI) and second-year ice(SYI) near Zhongshan Station, East Antarctica,were conducted for the first time from December 2011 to December 2012. Melt ponds appeared from early December 2011. Landfast ice partly broke in late January, 2012 after a strong cyclone. Open water was refrozen to form new ice cover in mid-February, and then FYI and SYI co-existed in March with a growth rate of 0.8 cm/d for FYI and a melting rate of 2.7 cm/d for SYI. This difference was due to the oceanic heat flux and the thickness of ice,with weaker heat flux through thicker ice. From May onward, FYI and SYI showed a similar growth by 0.5 cm/d.Their maximum thickness reached 160.5 cm and 167.0 cm, respectively, in late October. Drillings showed variations of FYI thickness to be generally less than 1.0 cm, but variations were up to 33.0 cm for SYI in March,suggesting that the SYI bottom was particularly uneven. Snow distribution was strongly affected by wind and surface roughness, leading to large thickness differences in the different sites. Snow and ice thickness in Nella Fjord had a similar "east thicker, west thinner" spatial distribution. Easterly prevailing wind and local topography led to this snow pattern. Superimposed ice induced by snow cover melting in summer thickened multi-year ice,causing it to be thicker than the snow-free SYI. The estimated monthly oceanic heat flux was ～30.0 W/m2 in March–May, reducing to ～10.0 W/m2 during July–October, and increasing to ～15.0 W/m2 in November. The seasonal change and mean value of 15.6 W/m2 was similar to the findings of previous research. The results can be used to further our understanding of landfast ice for climate change study and Chinese Antarctic Expedition services.     

基于中国北极考察布放浮标观测和数值模拟的海冰和积雪厚度变化分析

Sea ice and the snow pack on top of it were investigated using Chinese National Arctic Research Expedition(CHINARE) buoy data.Two polar hydrometeorological drifters,known as Zeno? ice stations,were deployed during CHINARE 2003.A new type of high-resolution Snow and Ice Mass Balance Arrays,known as SIMBA buoys,were deployed during CHINARE 2014.Data from those buoys were applied to investigate the thickness of sea ice and snow in the CHINARE domain.A simple approach was applied to estimate the average snow thickness on the basis of Zeno~ temperature data.Snow and ice thicknesses were also derived from vertical temperature profile data based on the SIMBA buoys.A one-dimensional snow and ice thermodynamic model(HIGHTSI) was applied to calculate the snow and ice thickness along the buoy drift trajectories.The model forcing was based on forecasts and analyses of the European Centre for Medium-Range Weather Forecasts(ECMWF).The Zeno~ buoys drifted in a confined area during 2003–2004.The snow thickness modelled applying HIGHTSI was consistent with results based on Zeno~ buoy data.The SIMBA buoys drifted from 81.1°N,157.4°W to 73.5°N,134.9°W in 15 months during2014–2015.The total ice thickness increased from an initial August 2014 value of 1.97 m to a maximum value of2.45 m before the onset of snow melt in May 2015;the last observation was approximately 1 m in late November2015.The ice thickness based on HIGHTSI agreed with SIMBA measurements,in particular when the seasonal variation of oceanic heat flux was taken into account,but the modelled snow thickness differed from the observed one.Sea ice thickness derived from SIMBA data was reasonably good in cold conditions,but challenges remain in both snow and ice thickness in summer.     

Sea-ice retreat in the Sea of Okhotsk and the ice–ocean albedo feedback effect on it

Sea-ice retreat processes are examined in the Sea of Okhotsk. A heat budget analysis in the sea-ice zone shows that net heat flux from the atmosphere at the water surface is about 77 W m⁻² on average in the active ice melt season (April) due to large solar heating, while that at the ice surface is about 12 W m⁻² because of the difference in surface albedo. The temporal variation of the heat input into the upper ocean through the open water fraction corresponds well to that of the latent heat required for ice retreat. These results suggest that heat input into the ice–upper ocean system from the atmosphere mainly occurs at the open water fraction, and this heat input into the upper ocean is an important heat source for ice melting. The decrease in ice area in the active melt season (April) and the geostrophic wind just before the melt season (March) show a correlation: the decrease is large when the offshoreward wind is strong. This relationship can be explained by the following process. Once ice concentration is decreased (increased) by the offshoreward (onshoreward) wind just before the melt season, solar heating of the upper ocean through the increased (decreased) open water fraction is enhanced (reduced), leading to (suppressing) a further decrease in ice concentration. This positive feedback is regarded as the ice–ocean albedo feedback, and explains in part the large interannual variability of the ice cover in the ice melt season.     

Annual and Seasonal Variations of the Sea Surface Heat Fluxes in the East Asian Marginal Seas

Based on the twice-daily marine atmospheric variables which were derived mostly from the weather maps for 18 years period from 1978 to 1995, the surface heat flux over the East Asian marginal seas was calculated at 0.5°×0.5° grid points twice a day. The annual mean distribution of the net heat flux shows that the maximum heat loss occurs in the central part of the Yellow Sea, along the Kuroshio axis and along the west coast of the northern Japanese islands. The area off Vladivostok turned out to be a heat-losing region, however, on the average, the amount of heat loss is minimum over the study area and the estuary of the Yangtze River also appears as a region of the minimum heat loss. The seasonal variations of heat flux show that the period of heat gain is longest in the Yellow Sea, and the maximum heat gain occurs in June. The maximum heat loss occurs in January over the study area, except the Yellow Sea where the heat loss is maximum in December. The annual mean value of the net heat flux in the East/Japan Sea is −108 W/m² which is about twice the value of Hirose et al. (1996) or about 30% higher than Kato and Asai (1983). For the Yellow Sea, it is about −89 W/m² and it becomes −75 W/m² in the East China Sea. This increase in values of the net heat flux comes mostly from the turbulent fluxes which are strongly dependent on the wind speed, which fluctuates largely during the winter season. This revised version was published online in August 2006 with corrections to the Cover Date.     

2018年北极太平洋区域夏季海冰物理及光学性质的研究

The reduction in Arctic sea ice in summer has been reported to have a significant impact on the global climate. In this study, Arctic sea ice/snow at the end of the melting season in 2018 was investigated during CHINARE-2018, in terms of its temperature, salinity, density and textural structure, the snow density, water content and albedo, as well as morphology and albedo of the refreezing melt pond. The interior melting of sea ice caused a strong stratification of temperature, salinity and density. The temperature of sea ice ranged from –0.8℃ to 0℃, and exhibited linear cooling with depth. The average salinity and density of sea ice were approximately 1.3 psu and 825 kg/m~3, respectively, and increased slightly with depth. The first-year sea ice was dominated by columnar grained ice. Snow cover over all the investigated floes was in the melt phase, and the average water content and density were 0.74% and 241 kg/m~3, respectively. The thickness of the thin ice lid ranged from 2.2 cm to 7.0 cm, and the depth of the pond ranged from 1.8 cm to 26.8 cm. The integrated albedo of the refreezing melt pond was in the range of 0.28–0.57. Because of the thin ice lid, the albedo of the melt pond improved to twice as high as that of the mature melt pond. These results provide a reference for the current state of Arctic sea ice and the mechanism of its reduction.     

Meridional energy flux in the Arctic from data of the radiosonde archive IGRA

The meridional energy transport into high latitudes of the Northern Hemisphere is an important climate-forming factor in the Arctic. This work presents the results of calculating the meridional energy flux across 70° N based on the Integrated Global Radiosonde Archive (IGRA) data from the radio sounding of the atmosphere. The long-term mean energy flux over the period 1992–2007 in the layer from the Earth’s surface to 30 hPa is 70.6 W m^?2. The fraction of the sensible heat flux is 23.2 W m^?2, i.e., 33% of the total energy flux; the fraction of the latent heat flux is 28.0 W m^?2 (40% of the total energy flux); the fraction of the potential energy is 20.0 W m^?2 (27%); and the fraction of the kinetic energy is 0.53 W m^?2, i.e., less than 1% of the total energy flux. The vertical structure of the flux shows that the main energy transport into the Arctic takes place in the middle troposphere-lower stratosphere layer, whereas the energy is transported mainly out of the Arctic in the lower troposphere, which agrees well with the schematic notion about the polar circulation cell. The spatial structure of the flux shows that the key regions with a positive (directed into the Arctic) energy flux are located in the vicinity of 160° E (the northwestern part of Eurasia, Pacific sector) and 50° W (Greenland sector). The regions with a negative (directed out of the Arctic) energy flux are located near 120° W (Canadian Arctic Archipelago) and from 20° E to 90° E (Atlantic sector). In the period from 1992 to 2007, the meridional energy transport into the Arctic weakened by ?0.26 W m^?2 yr^?1. The changes were mutually correlated; namely, positive and negative energy fluxes weakened in amplitude, almost without changing their locations.