利用多年卫星资料分析青岛和茂名附近海域海面风场变化

文章围绕青岛海洋经济开发和茂名海洋气象基础科学综合试验对附近海域基础环境数据的需求,利用多年卫星资料开展了青岛和茂名附近海域海面风场统计分析研究。结果表明:青岛附近海域风向频率和风速频率最大的方向均为北,相对应的平均风速值为3.5m/s;茂名附近海域风向频率最大方向为ENE,风速频率最大方向为NNE,茂名附近海域平均风速值全年均大于青岛附近海域平均风速值,研究成果可为开展海洋气象业务观测和科学试验提供基础数据支持,为"一带一路"海上风能资源开发与利用提供决策支撑。     

吕宋海峡黑潮脱落涡旋的特征分析

涡旋脱落在西太平洋和南海的海水属性交换中起到重要作用。为研究吕宋海峡附近海域由黑潮脱落并进入南海的涡旋特征,本文采用1993—2014年法国空间局(AVISO)多卫星融合海面高度距平(SLA)和绝对动力地形(ADT)全球网格化延时数据,美国国家海洋数据中心(NODC)的WOA13年平均温盐剖面气候数据,以及1993—2010年SODA2.2.4月平均海洋同化数据集,并分析了黑潮脱落涡旋与大尺度环流的关系。结果表明:(1)暖涡脱落数量远多于冷涡数量,且脱落的冷涡绝大部分在黑潮西侧边缘生成,而脱落的暖涡则大部分在黑潮控制区生成。(2)冷涡、暖涡脱落时的平均半径、平均振幅相近,但是冷涡的平均生命、平均迁移距离约为暖涡的一半。(3)冷涡不是每年都有脱落,主要在冬季脱落;暖涡则每年均有脱落,主要发生在秋季。(4)脱落涡旋数量与脱落时的黑潮路径类型相关。(5)脱落涡旋的平均西行速度为5.8cm/s,与斜压第一模态长Rossby波波速及大尺度环流的西向平流流速之和相近。     

QuikSCAT卫星遥感风场可靠性分析及其揭示的中国近海风速分布

利用2000—2009年美国国家航空航天局(NASA)在中国近海海域(0°~45°N,105°~135°E)的QuikSCAT卫星遥感风场资料与近海测风塔(位于上海近海)、海上石油平台(位于东海和渤海)、岛屿站(南海珊瑚岛和西沙海边观测塔)的实测风场资料进行对比分析,检验了QuikSCAT卫星遥感风场资料在中国近海海域的可靠性。研究结果如下:各站点实测风速与站点位置以及站点附近的QuikSCAT卫星遥感风场资料相关系数均在0.7以上;QuikSCAT卫星遥感风场资料与海上石油平台的风速均方根误差较小(约1.5 m/s);其年均值均大于实测值,差值范围是0.1~1.3 m/s;其Weibull形状参数K与海上石油平台以及近海测风塔的K值较为接近,表明QuikSCAT卫星遥感风场资料各风速段的频次分布形态与观测站的实测值基本吻合,QuikSCAT卫星遥感风场资料能基本合理地反映出中国近海风速的分布状况。利用QuikSCAT卫星遥感风场资料分析了中国近海及其邻近水域风速的空间分布特征:(1)台湾海峡是中国近海风速最大的区域,从台湾海峡向东北至日本海,往西南至南海北部115°E附近和巴林塘海峡为风速的次大值区;(2)28°N到长江入海口的东海海域年均风速为7.0~7.5 m/s,在黄海和渤海为5.5~7.0 m/s,在南海北部自东向西由8.5 m/s递减为6.0 m/s,北部湾最大风速区位于东方附近海域。     

HY-2B卫星载荷联合观测海面阵风的一种反演方法

国内外对海上阵风的研究并不多,且大多集中在阵风预报和应用研究方面,对于海洋阵风数据的获取技术未见文献系统论述。本文利用HY-2B卫星雷达高度计观测的后向散射系数,结合校正微波辐射计观测的亮度温度信息,提出联合反演阵风风速的方法。两个遥感载荷联合反演得到的阵风风速与2019–2021年美国国家浮标数据中心（NDBC）浮标数据进行真实性检验,结果显示:阵风风速均方根误差（RMSE）为0.98 m/s,相关系数为0.82;基于本方法利用国外同类卫星Jason-3得到的阵风风速与2016–2018年NDBC浮标数据的RMSE为0.96 m/s,相关系数为0.88。本文在HY-2B卫星雷达高度计海面风速观测的基础上,纳入同一卫星平台校正微波辐射计的同步观测信息联合实现了海面阵风的观测,数据的比对结果证明文中方法具有较高的观测精度。同时,该方法对于具有相同观测体制的国内外卫星也适用。     

冬季黑潮延伸体区域涡旋尺度海面温度季节内变化特征

利用ERA_interim再分析资料以及OISST高分辨率海面温度(Sea surface temperature,SST)卫星观测数据,通过小波分析、二维模态相关(Pattern correlation)等方法系统地分析了黑潮延伸体区域涡旋尺度SST信号的季节内变化特征。发现涡旋尺度SST季节内变化信号在冬季最强、夏季最弱。该信号主要分布在黑潮-亲潮海洋锋面区域,在日本沿岸振幅最强、并沿黑潮、亲潮锋面向东延伸,其标准差高达1℃,是该海域冬季SST的重要变化信号。涡旋尺度SST的季节内变化周期以40～100 d周期为主,会引起大气边界层的响应,激发海气界面湍流热通量、海面气温、边界层高度等同位相的季节内变化。在该区域涡旋信号较强个数较多的时间段SST异常的季节内变化信号更明显,边界层大气与涡旋尺度SST季节内变化信号的相关程度更大。     

中国海及邻近海域卫星观测资料同化试验

利用1个基于POMgcs海洋模式和多重网格三维变分同化方法建立的中国海及邻近海域海面高与三维温盐流数值预报模型,通过一系列数值试验,研究了同化卫星测高和卫星遥感海面温度观测资料对该模型预报能力的影响。试验结果表明,同化卫星测高资料可明显改善海面高度与三维温度和盐度的分析预报效果,使1 200 m以上的温度预报误差减小0.16℃,并能有效提高对海洋中尺度现象的预报能力;同化卫星遥感海面温度对100 m以上的温度和盐度的预报效果有所改善,可使海面温度的预报误差减小10%。     

卫星SAR图像的海洋大气边界层特征参数的反演

针对中国南海ERS-2／SAR三维对流涡旋图像,分别利用Young的协方差方法和惯性子域方法,反演了对流海洋大气边界层高度、海洋大气边界层垂直对流尺度速度、表面浮力通量、Obukhov长度、稳定性校正后的海面风速和海面风速稳定性校正因子。结果表明,2种反演方法的结果符合较好,在一定程度上说明利用协方差方法和惯性子域方法从卫星SAR图像反演海洋大气边界层参数的有效性。     

梳状高度计编队探测海洋涡旋仿真分析

单颗星下点高度计卫星可进行一维海面高度测量，但时空分辨率较低，对涡旋的探测能力不足。梳状高度计卫星编队可提升海面高度时空分辨率，增强对涡旋的识别探测能力。为定量分析高度计编队下的涡旋探测能力，本文在黑潮延伸体海域开展观测仿真实验，对海平面高度进行沿轨采样，经优化插值算法处理后得到逐天海面高度异常数据场并进行涡旋识别与统计分析。此外，本文还同步开展了单颗星下点高度计和成像高度计卫星涡旋探测仿真实验及对比分析。结果表明，梳状高度计卫星编队在涡旋数量、极性、半径及振幅等方面识别效果较好，与成像高度计探测结果相近。     

利用卫星数据反演海洋蒸发波导的研究

海洋蒸发波导是海洋大气中的一种异常折射结构,通常是由观测站点的传感器数据计算得出。但这种方法只能得到单点的海洋蒸发波导,无法满足我们对大范围海域的海洋信息需求。文中提出了一种利用卫星遥感数据反演海洋蒸发波导的方法,并利用观测站点的数据对反演结果进行订正,从而得到了完整的大面积海区蒸发波导环境。     

基于HY-2 卫星观测分析南海风浪关系

海面有效波高(H1/3)是表征海浪的重要参数,随着卫星遥感技术的发展,雷达高度计已成为获取海面有效波高的重要手段,但也只能对卫星星下点轨迹处进行有效观测,远无法满足大范围应用的需求.本研究结合2013年10月HY-2雷达高度计观测的海面有效波高和微波散射计观测的海面风场资料,分别对高、低风速下风浪数据进行拟合,建立了适用于0~40 m/s风速范围内的南海海域风浪关系模型,经模型比对和结果验证,结果表明,基于HY-2卫星数据分析建立的南海海域风浪关系模型是可信的,特别是低风速的风浪模型与实测数据建立的风浪模型具有很好的一致性;根据建立的风浪关系模型,从卫星散射计大面观测的海面风场出发,能推算出风浪条件下海面有效波高的大面信息,数据覆盖远高于卫星高度计的星下点观测,能为分析和预报海浪、风暴潮灾害提供大范围的海面有效波高信息.     

MODIS渤海海冰遥感资料反演

鉴于渤海海冰监测和预报对海冰卫星遥感数字化产品的迫切需求,本文利用MODIS的1B级数据进行渤海海冰参数反演,提供海冰遥感图像和海冰密集度、冰厚数值产品,作为渤海海冰监测和海冰数值预报初始场的重要信息来源,以及海冰预报质量检验的参考依据之一。反演结果表明,其各通道对海冰性质有很好的反映,资料信号比较稳定,对不同密集度和厚度的冰有较好的区分,相对NOAA/AVHRR和HY-1A资料有更好的实际应用价值;Terra/MODIS和HY-1A/COCTS海冰遥感反演结果对比也为HY-1A系列卫星海冰遥感的改进和提高提供有益的启示。     

On a translating and transversely oscillating cylinder. Part 1—The effect of the strouhal number on the hydrodynamic forces and the near-wake structure

The problem of unsteady, laminar flow past a circular cylinder which starts translating and oscillating impulsively from rest in a viscous fluid is numerically investigated at a Reynolds number of R = 10³. The flow is incompressible and two-dimensional, and the cylinder oscillations are harmonic. The transverse oscillations are only allowed when the maximum oscillatory-to-translational velocity ratio is 0.5. The investigation is based on an implicit finite difference scheme for integrating the unsteady Navier-Stokes equations together with the mass-conservation equation in their vorticity stream function formulation. A non-inertial coordinate transformation is used so that the grid mesh remains fixed relative to the accelerating cylinder. Present calculations are performed within the range of sufficiently large oscillation amplitude to induce separation. The time variation of the in-line and transverse force coefficients are presented. The study also focuses on the laminar asymmetric flow structure in the near-wake region. In this flow regime, it is found that there is alternate shedding of vortices from either side of the cylinder over an oscillation cycle (as predicted experimentally); this is the classical mode of vortex shedding leading to formation of the Kármán street.     

Island mass effect in the Juan Fernández Archipelago (33°S), Southeastern Pacific

Spatial and temporal variability of the island mass effect (IME; defined as local increases of phytoplankton associated with the presence of islands) at the Juan Fernández Archipelago (JFA) is analyzed using chlorophyll-a (Chl-a) satellite data, altimetry, sea surface temperature, wind, geostrophic currents and net heat flux over a ten year period (2002–2012). The the JFA islands (Robinson Crusoe-Santa Clara (RC-SC) and Alejandro Selkirk (AS)) present wakes with significant Chl-a increases, mainly during spring time. These wakes can reach Chl-a values of one order of magnitude higher (~1 mg m⁻³) than the surrounding oligotrophic waters (<0.1 mg m⁻³). The wakes are similar to von Kármán vortex streets which have been used to explain the impact of IME on Chl-a increases in numerical models. The wakes are formed from a high productivity area in the lee of the island, extending to the oceanic region as high Chl-a patches associated with submesoscale eddies that are detached from the islands and connected by less-productive zones. This pattern coincides with previous models that predict the effects of island-generated flow perturbations on biological production variability. The IME is a recurrent feature of islands that has even been observed in decadal average fields. In such average fields, the Chl-a values in RC-SC and AS islands can exceed values found in a Control Zone (a zone without islands) by ~50% and 30%, respectively. Seasonal and interannual variability reveals that, as a consequence of the IME, the winter Chl-a maximum associated with the development of winter convection and mesoscale eddies that propagate from the continental zone, promote that the Chl-a maximum extends towards spring. The IME has an impact on the island on both a local as well as a more regional scale that affects an area of ~40,000 km² (1°Latitude×4°Longitude) centered on the islands. The transport of high productivity patches associated with submesoscale eddies may be responsible for IME propagation at a regional scale. Around the islands, the presence of a weak oceanic incident flow and strong and recurrent wind-wakes, suggest that the generation of Chl-a wakes result from a combined effect between both forcings.     

On the nature of oceanic eddies shed by the Island of Gran Canaria

We use hydrographic and buoy data to compare the initial temperature fields and Lagrangian evolution of water parcels in two vortices generated by the southward flowing Canary Current passing around the island of Gran Canaria Island. One vortex is anticyclonic, shed in June 1998 as the result of an incident current of about 0.05 m s⁻¹, and the second one is cyclonic, shed in June 2005 with the impinging current estimated as 0.03 m s⁻¹. The two vortices exhibit contrasting characteristics yet display some important similarities. The isopycnals are depressed in the core of the anticyclonic vortex, at least down to a depth of 700 m, whilst they dome up in the core of the cyclonic vortex but only down to 450 m. In the top 300 m the depression/doming of the isotherms is similar for both vortices, with a maximum vertical displacement of the isotherm of about 80 m, which correspond to temperature anomalies of some 2.5 °C at a given depth. A simple method is developed to obtain the initial orbital velocity field from the temperature data, from which we estimate peak values of 0.7 and 0.5 m s⁻¹ for the anticyclonic and cyclonic vortices, respectively. The buoys, three for the anticyclonic vortex and two for the cyclonic one, were drougued at 100 m depth, below the surface mixed layer, and their initial velocities are consistent with the above values. In both vortices, the buoys revolve either within a central core, where the rotation rate remains stable and large for several weeks, or in an outer ring, where the rotation rate is significantly smaller and displays large radial fluctuations. Within the inner core the anticyclonic vortex has significant inward radial velocity, while the cyclonic vortex has near-zero radial mean motions. The cyclonic vortex rotates more slowly than the anticyclonic, their initial periods being 4.5 and 2.5 days, respectively. A simple axisymmetric model with radial diffusion (coefficient K_h≅25 m² s⁻¹) and advection reproduces the observations reasonably well, the diffusive effect being more important than that resulting from the observed radial advection. The model also supports the hypothesis that the rotation rate of cyclonic vortices is less than that of anticyclonic vortices, as otherwise they would become inertially unstable. Both the buoys data and sea surface temperature images confirm that the vortices evolve from youth to maturity, as the cores shrink and the outer rings expands, and then to a decay stage, as the core rotation rates decrease, though frequent interactions with other mesoscale structures result in more accelerated aging. Despite these interaction they last many months as coherent structures south of the Canary Islands.     

HY1C/1D海表温度对Terra/Aqua产品的可替代性分析

通过卫星遥感获取的海表温度（SST）产品已经成为海洋和大气研究中的重要数据源,我国海洋水色遥感卫星（HY1C和HY1D）的海洋水色水温扫描仪（COCTS）具有两个热红外通道,可反演全球SST遥感产品。对比Terra和Aqua卫星的中分辨率成像光谱仪（MODIS）的SST产品,分析COCTS海表温度产品对MODIS相应产品的可替代性。比较了两种卫星的全球SST单日和月平均融合产品的图像空间结构,分析了匹配像元SST值的离散度,统计了HY1C/1D的误差结果,讨论了HY1C与HY1D产品的一致性、不同质量控制方案对SST产品影响以及遥感产品质量对昼夜SST变化研究影响等问题。结果表明,以2020年6月SST(Terra)为真值,HY1C白天SST的单日全球遥感产品的平均偏差、绝对偏差、均方根误差和相关系数分别为0.04℃、0.60℃、0.78℃和0.98,夜晚SST的单日全球遥感产品的平均偏差、绝对偏差、均方根误差和相关系数分别为-0.16℃、0.78℃、0.95℃和0.86。以2020年6月SST(Aqua)为真值,HY1D白天SST的单日全球遥感产品的平均偏差、绝对偏差、均方根误差和相...     

HY-1/COCTS波段信息含量与波段相关性分析--以辽东湾为例

以2003年4月辽东湾的HY-1／COCTS数据为例,采用主成分分析法和相关分析法,从光谱的不同结构、主成分、相关性等方面来揭示COCTS图像的光谱信息。COCTS的10个波段中,9,10两个波段与其他几个波段相比,包含的信息较丰富。在其他的8个波段中,6,7波段的信息较为丰富,6波段对主成分的贡献大。这8个波段与9波段和10波段的相关性为负相关,表明它们间的重叠少。5波段与1波段、7波段与6波段、6波段与2波段的相关性较高,这说明这些数据间的彼此重叠较多,故在进行波段图像的假彩色合成时一般以8,4,3或8,4,1波段的组合为佳。以OIF指数来检查融合方法的可行性。     

中国HY-2卫星雷达高度计有效波高真实性检验

Chinese Haiyang-2(HY-2) satellite is the first Chinese marine dynamic environment satellite. The dual-frequency(Ku and C band) radar altimeter onboard HY-2 has been working effective to provide operational significant wave height(SWH) for more than three years(October 1, 2011 to present).We validated along-track Ku-band SWH data of HY-2 satellite against National Data Buoy Center(NDBC) in-situ measurements over a time period of three years from October 1, 2011 to September 30, 2014, the root mean square error(RMSE) and mean bias of HY-2SWH is 0.38 m and(–0.13±0.35) m, respectively. We also did cross validation against Jason-2 altimeter SWH data,the RMSE and the mean bias is 0.36 m and(–0.22±0.28) m, respectively. In order to compare the statistical results between HY-2 and Jason-2 satellite SWH data, we validated the Jason-2 satellite radar altimeter along-track Ku-band SWH data against NDBC measurements using the same method. The results demonstrate the validation method in this study is scientific and the RMSE and mean bias of Jason-2 SWH data is 0.26 m and(0.00±0.26) m,respectively. We also validated both HY-2 and Jason-2 SWH data every month, the mean bias of Jason-2 SWH data almost equaled to zero all the time, while the mean bias of HY-2 SWH data was no less than –0.31 m before April2013 and dropped to zero after that time. These results indicate that the statistical results for HY-2 altimeter SWH are reliable and HY-2 altimeter along-track SWH data were steady and of high quality in the last three years. The results also indicate that HY-2 SWH data have greatly been improved and have the same accuracy with Jason-2SWH data after April, 2013. SWH data provided by HY-2 satellite radar altimeter are useful and acceptable for ocean operational applications.     

First six months quality assessment of HY-2A SCAT wind products using in situ measurements

The first Chinese microwave ocean environment satellite HY-2A, carrying a Ku-band scatteromenter （SCAT）, was successfully launched in August 2011. The first quality assessment of HY-2A SCAT wind products is presented through the comparison of the first 6 months operationally released SCAT products with in situ data. The in situ winds from the National Data Buoy Center （NDBC） buoys, R/V Polarstern, Aurora Australis, Roger Revelle and PY30-1 oil platform, were converted to the 10 m equivalent neutral winds. The temporal and spatial differences between the HY-2A SCAT and the in situ observations were limited to less than 5 min and 12.5 km. For HY-2A SCAT wind speed products, the comparison and analysis using the NDBC buoys yield a bias of-0.49 m/s, a root mean square error （RMSE） of 1.3 m/s and an increase negative bias with increasing wind speed observation above 3 m/s. Although less accurate of HY-2A SCAT wind direction at low winds, the RMSE of 19.19° with a bias of 0.92° is found for wind speeds higher than 3 m/s. These results are found consistent with those from R/Vs and oil platform comparisons. Moreover, the NDBC buoy comparison results also suggest that the accuracy of HY-2A SCAT winds is consistent over the first half year of 2012. The encouraging assessment results over the first 6 months show that wind products from HY-2A SCAT will be useful for scientific community.     

Wave–current interaction with a vertical square cylinder

A three-dimensional numerical model is developed in this study to investigate the problem of wave–current–body interaction. The model solves the spatially averaged Navier–Stokes equations. Turbulence effects are modeled by a subgrid-scale (SGS) model using the concept of large eddy simulation (LES). The model is employed to study the wave–current interaction with a square cylinder that is mounted on the bottom and vertically pierces the free surface. The force analysis demonstrates that the presence of waves can reduce both the strength and frequency of vortex shedding induced by a uniform current due to the nonlinear wave–current interaction. The free surface elevation, strain rates of the mean flow, and eddy viscosity are found to closely correlate with the mechanism of vortex shedding. It is also shown that when the vortex shedding is neglected in the calculation such as by the potential flow approach, one may significantly underestimate the magnitude of in-line force. The energy spectral analysis reveals that there exist initiating, growing, and decaying regions for shedding vortices around the cylinder. In the vortex initiating region, both coherent and turbulent structures are nearly two-dimensional that become three-dimensional in the vortex growing region. The kinetic energy of both coherent and turbulent motions is dissipated in the vortex decaying region, within which the mean flow gradually returns back to two-dimensional.     

基于NDVI与丰度关系的MODIS影像浒苔混合像元分解方法

MODIS影像是浒苔业务化卫星遥感监测的重要数据源,但其空间分辨率(250m)较低,混合像元效应导致传统的NDVI(Normalized Difference Vegetation Index)阈值法浒苔提取误差较大。为解决此问题,本文以准同步的较高分辨率的HJ-1CCD影像浒苔提取结果为基准,建立基于MODIS NDVI与浒苔像元丰度(浒苔覆盖面积占像元面积百分比)关系的浒苔信息提取模型,利用该模型可提取每个像元的丰度(即浒苔覆盖率),从而提高整景影像浒苔覆盖面积提取精度。该模型可为浒苔灾害等级的确定以及防灾减灾应急决策提供支持。