最优插值方法在西北太平洋海温同化中的应用研究

本文开发出一套与三维温盐流数值模式POM相匹配的最优插值数值同化模块,具有将船舶报资料和Argo海温数据加入数值模拟的能力.利用该同化系统,本文对2001～2002年进行了数值模拟试验,通过模拟结果和未加入同化的模拟结果的对比表明,该同化系统能有效地改进三维海温模拟结果,使之更为接近于观测值.     

基于ROMS和4DVAR的沿轨与网格化SSH数据同化效果评价

Remote sensing products are significant in the data assimilation of an ocean model. Considering the resolution and space coverage of different remote sensing data, two types of sea surface height(SSH) product are employed in the assimilation, including the gridded products from AVISO and the original along-track observations used in the generation. To explore their impact on the assimilation results, an experiment focus on the South China Sea(SCS) is conducted based on the Regional Ocean Modeling System(ROMS) and the four-dimensional variational data assimilation(4 DVAR) technology. The comparison with EN4 data set and Argo profile indicates that, the along-track SSH assimilation result presents to be more accurate than the gridded SSH assimilation, because some noises may have been introduced in the merging process. Moreover, the mesoscale eddy detection capability of the assimilation results is analyzed by a vector geometry–based algorithm. It is verified that, the assimilation of the gridded SSH shows superiority in describing the eddy's characteristics, since the complete structure of the ocean surface has been reconstructed by the original data merging.     

基于集合最优插值方法的南海北部海区的数据同化试验

基于ROMS(Regional Ocean Modeling System)模式建立一个南海北部集合最优插值的同化系统,并且利用2008年夏季SCOPE(Northern South China Sea Coastal Oceanographic Process Experiment)航次的温盐数据以及航次期间逐日OS...     

基于ROMS模式的南海SST与SSH四维变分同化研究

卫星遥感观测获得了大量高分辨率的海面实时信息,包括海面温度（SST）和海面高度（SSH）等,同化进入数值模式可有效提升模拟精度。本文基于ROMS模式与四维变分同化方法（4DVAR）,使用AVHRR SST和AVISO SSH数据,开展了南海区域同化实验。为检验同化的效果,分别利用HYCOM再分析资料和Argo温盐实测数据分析了同化结果的海面高度、流场及温盐剖面的精度。对比结果表明,SST和SSH的同化能够改善ROMS的模拟结果:同化后海面高度场能够更为准确地捕捉海洋的中尺度特征,与HYCOM海面高度再分析资料相比,平均绝对偏差和均方根误差分别为0.054 m和0.066 m;与HYCOM 10 m层流场相比,东向与北向流速平均绝对偏差分别为0.12 m/s和0.11 m/s,相比未同化均提升约0.01 m/s;温盐同化结果与Argo温盐实测具有较高的一致性,温度和盐度平均绝对偏差为0.45℃、0.077,均方根误差为0.91℃、0.11,单个的温盐廓线对比说明,同化结果与HYCOM再分析资料精度相当。     

基于地球系统模式的局地化粒子滤波器与集合卡尔曼滤波器同化实验

粒子滤波器(PF)是一种非常具有应用前景的非线性资料同化方法。但由于其算法本身存在的粒子退化问题,目前尚未被广泛地应用于大型地球物理模式。目前主流的集合同化系统仍然倾向于使用集合卡尔曼滤波器(EnKF)及其衍生方法。一种新近被提出的局地化粒子滤波器(LPF)在经典的粒子滤波器算法中引入局地化技术,可以使用较小的计算成本有效地避免退化问题,具有非常大的业务应用潜力。本文在全耦合的通用地球系统模式中开展了LPF和EnKF的同化实验,同化资料为模拟的卫星海表温度资料。着重考察了不同局地化参数对两种方法的不同影响,对比了局地化粒子滤波器与集合卡尔曼滤波器的同化效果差异。比较的结果表明,LPF的同化效果对于局地化参数的选择非常敏感,在使用最优局地化参数的条件下,LPF能达到与EnKF相当甚至优于后者的同化效果,并具有较大的改进空间。     

集合同化方法在黄海海域海表面温度同化中的应用

The effects of sea surface temperature(SST) data assimilation in two regional ocean modeling systems were examined for the Yellow Sea(YS). The SST data from the Operational Sea Surface Temperature and Sea Ice Analysis(OSTIA) were assimilated. The National Marine Environmental Forecasting Center(NMEFC) modeling system uses the ensemble optimal interpolation method for ocean data assimilation and the Kunsan National University(KNU) modeling system uses the ensemble Kalman filter. Without data assimilation, the NMEFC modeling system was better in simulating the subsurface temperature while the KNU modeling system was better in simulating SST. The disparity between both modeling systems might be related to differences in calculating the surface heat flux, horizontal grid spacing, and atmospheric forcing data. The data assimilation reduced the root mean square error(RMSE) of the SST from 1.78°C(1.46°C) to 1.30°C(1.21°C) for the NMEFC(KNU) modeling system when the simulated temperature was compared to Optimum Interpolation Sea Surface Temperature(OISST) SST dataset. A comparison with the buoy SST data indicated a 41%(31%) decrease in the SST error for the NMEFC(KNU) modeling system by the data assimilation. In both data assimilative systems, the RMSE of the temperature was less than 1.5°C in the upper 20 m and approximately 3.1°C in the lower layer in October. In contrast, it was less than 1.0°C throughout the water column in February. This study suggests that assimilations of the observed temperature profiles are necessary in order to correct the lower layer temperature during the stratified season and an ocean modeling system with small grid spacing and optimal data assimilation method is preferable to ensure accurate predictions of the coastal ocean in the YS.     

OSTIA数据在中国近海业务化环流模型中的同化应用

The prediction of sea surface temperature(SST) is an essential task for an operational ocean circulation model. A sea surface heat flux, an initial temperature field, and boundary conditions directly affect the accuracy of a SST simulation. Here two quick and convenient data assimilation methods are employed to improve the SST simulation in the domain of the Bohai Sea, the Yellow Sea and the East China Sea(BYECS). One is based on a surface net heat flux correction, named as Qcorrection(QC), which nudges the flux correction to the model equation; the other is ensemble optimal interpolation(En OI), which optimizes the model initial field. Based on such two methods, the SST data obtained from the operational SST and sea ice analysis(OSTIA) system are assimilated into an operational circulation model for the coastal seas of China. The results of the simulated SST based on four experiments, in 2011, have been analyzed. By comparing with the OSTIA SST, the domain averaged root mean square error(RMSE) of the four experiments is 1.74, 1.16, 1.30 and 0.91°C, respectively; the improvements of assimilation experiments Exps 2, 3 and 4 are about 33.3%, 25.3%, and 47.7%, respectively.Although both two methods are effective in assimilating the SST, the En OI shows more advantages than the QC,and the best result is achieved when the two methods are combined. Comparing with the observational data from coastal buoy stations, show that assimilating the high-resolution satellite SST products can effectively improve the SST prediction skill in coastal regions.     

南海海平面长期变化及其趋势研究

On the basis of the satellite maps of sea level anomaly(MSLA) data and in situ tidal gauge sea level data,correlation analysis and empirical mode decomposition(EMD) are employed to investigate the applicability of MSLA data,sea level correlation,long-term sea level variability(SLV) trend,sea level rise(SLR) rate and its geographic distribution in the South China Sea(SCS).The findings show that for Dongfang Station,Haikou Station,Shanwei Station and Zhapo Station,the minimum correlation coefficient between the closest MSLA grid point and tidal station is 0.61.This suggests that the satellite altimeter MSLA data are effective to observe the coastal SLV in the SCS.On the monthly scale,coastal SLV in the western and northern part of SCS are highly associated with coastal currents.On the seasonal scale,SLV of the coastal area in the western part of the SCS is still strongly influenced by the coastal current system in summer and winter.The Pacific change can affect the SCS mainly in winter rather than summer and the affected area mostly concentrated in the northeastern and eastern parts of the SCS.Overall,the average SLR in the SCS is 90.8 mm with a rising rate of(5.0±0.4) mm/a during1993–2010.The SLR rate from the southern Luzon Strait through the Huangyan Seamount area to the Xisha Islands area is higher than that of other areas of the SCS.     

Ekman抽吸对南海海面升高变化的贡献

利用NECP海表风场资料和T/P卫星高度计资料,研究风场通过Ekman抽吸对南海海面升高变化的贡献。通过分析发现:对于南海200 m以深区域的平均海面升高的变化,Ekman抽吸在1月和2月份起到主要作用,其他月份和海面升高变化反号。在年际尺度上,Ekman抽吸的EOF第二模态对海面高度的第一模态具有紧密的联系。     

中国东海海平面长期变化与趋势研究

From the analyses of the satellite altimeter Maps of Sea Level Anomaly(MSLA) data, tidal gauge sea level data and historical sea level data, this paper investigates the long-term sea level variability in the East China Sea(ECS).Based on the correlation analysis, we calculate the correlation coefficient between tidal gauge and the closest MSLA grid point, then generate the map of correlation coefficient of the entire ECS. The results show that the satellite altimeter MSLA data is effective to observe coastal sea level variability. An important finding is that from map of correlation coefficient we can identify the Kuroshio. The existence of Kuroshio decreases the correlation between coastal and the Pacific sea level. Kurishio likes a barrier or a wall, which blocks the effect of the Pacific and the global change. Moreover, coastal sea level in the ECS is mainly associated with local systems rather than global change. In order to calculate the long-term sea level variability trend, the empirical mode decomposition(EMD) method is applied to derive the trend on each MSLA grid point in the entire ECS. According to the 2-D distribution of the trend and rising rate, the sea level on the right side of the axis of Kuroshio rise faster than in its left side. This result supports the barrier effect of Kuroshio in the ECS. For the entire ECS, the average sea level rose 45.0 mm between 1993 and 2010, with a rising rate of(2.5±0.4) mm/a which is slower than global average.The relatively slower sea level rising rate further proves that sea level rise in the ECS has less response to global change due to its own local system effect.     

The impact of mean dynamic topography on a sea-level anomaly assimilation in the South China Sea based on an eddy-resolving model

The sea-level anomaly (SLA) from a satellite altimeter has a high accuracy and can be used to improve ocean state estimation by assimilation techniques. However, the lack of an accurate mean dynamic topography (MDT) is still a bothersome issue in an ocean data assimilation. The previous studies showed that the errors in MDT have significant impacts on assimilation results, especially on the time-mean components of ocean states and on the time variant parts of states via nonlinear ocean dynamics. The temporal-spatial differences of three MDTs and their impacts on the SLA analysis are focused on in the South China Sea (SCS). The theoretical analysis shows that even for linear models, the errors in MDT have impacts on the SLA analysis using a sequential data assimilation scheme. Assimilation experiments, based on EnOI scheme and HYCOM, with three MDTs from July 2003 to June 2004 also show that the SLA assimilation is very sensitive to the choice of different MDTs in the SCS with obvious differences between the experimental results and observations in the centre of the SCS and in the vicinity of the Philippine Islands. A new MDT for assimilation of SLA data in the SCS was proposed. The results from the assimilation experiment with this new MDT show a marked reduction (increase) in the RMSEs (correlation coefficient) between the experimental and observed SLA. Furthermore, the subsurface temperature field is also improved with this new MDT in the SCS.     

南海北部海面高度的季节内变异及其传播特征 *

利用25年(1993—2017)的卫星高度计资料, 采用复经验正交函数(complex empirical orthogonal function, CEOF)方法, 分析南海北部海区海面高度季节内变异的时空分布及传播特征。标准差分析表明, 南海北部海面高度的季节内变异(intra-seasonal variability of sea level anomalies, SLA-ISV)在沿陆坡外侧区较强, 且SLA-ISV表现出明显的季节性变化, 冬半年强于夏半年。CEOF前两个主要模态能较好地揭示研究海区SLA-ISV的时空分布及其传播特征, 并表明SLA-ISV的强度受到季节性变化和年际变化的调制。全年CEOF的第一模态揭示SLA-ISV从台湾岛西南至西沙群岛以东区域的冬半年西南向传播特征; 而全年CEOF的第二模态则表现了SLA-ISV分别在台湾岛西南和东沙群岛西南的西南向传播特征。南海北部中尺度涡季节变化统计分析表明, CEOF的分解结果与南海北部的涡旋活动一致。     

南海周边相对海平面变化特征及2004年苏门答腊地震影响分析

利用南海周边1989-2014年的潮汐资料和GPS长期观测资料,分析了南海周边相对海平面变化特征,以及2004年苏门答腊地震对该区域相对海平面变化的影响。研究结果表明,南海周边的相对海平面变化以上升为主,平均上升速率（4.53±0.20） mm/a,高于全球平均速率,且2004年后上升趋势加剧;南海周边相对海平面呈现6类较典型的变化特征,并存在与板块构造相对应的分区聚集现象,形成了中国东南和越南沿海、马来半岛、加里曼丹岛北部、菲律宾群岛等4个变化特征区。受2004年苏门答腊大地震的影响,马来半岛、南沙和西沙海域的地壳形变由上升趋势转为下沉,加剧了相对海平面的上升;中国东南沿海和菲律宾群岛受地震影响较小;越南沿海和加里曼丹岛北部区域的地震影响还有待进一步研究。     

全球变暖背景下南海动力海平面变化

由于分辨率不足等原因,当前大部分全球耦合气候模式对南海等海洋区域的模拟能力仍然较低。本文基于超高分辨率（Ultra high-resolution） CESM-UHR耦合模式（大气和海洋水平分辨率分别达到约25 km和约10 km）研究了南海动力海平面对全球变暖的响应。研究发现：（1） CESM-UHR能够较好地模拟出南海冬、夏季节性动力海面高度和表层环流变化;（2）在四倍二氧化碳试验下,冬季南海动力海平面变化呈现出中部低、近岸高的分布特征;夏季则呈现出西北部低、东南部高的分布特征,分别对应冬、夏表层地转流增强趋势;（3）冬、夏动力海平面变化特征与风应力旋度变化具有很好的对应关系;（4）全球变暖下南海海平面变化存在季节循环放大效应,这将增大南海极端水位灾害风险。     

Delving into three-dimensional structure of the West Luzon Eddy in a regional ocean model

The three-dimensional structure and associated dynamics of the prominent cold (cyclonic) West Luzon Eddy (WLE) were investigated by a high-resolution regional ocean model. The WLE was horizontally and vertically heterogeneous, exhibiting asymmetric structures in the circulation, vorticity, vertical motion and energy distributions within the eddy. The asymmetry was mainly attributed to the existence of an eddy dipole formed by a coexisting warm (anti-cyclonic) eddy to the south of the WLE. Analysis of the momentum balance revealed that the coexistence of two eddies intensified barotropic pressure gradients in the southern WLE to locally enhance the eastward jet. The positive (negative) vorticity of the jet strengthened (weakened) the eddy in the southern sector (periphery), which, together with the formation of a subsurface density front, intensified (suppressed) the corresponding upward motion and cooling. The baroclinic pressure gradients opposed the dominant barotropic components and spun down the eddy at greater depths with stronger weakening in the southern sector near the front. Asymmetric energy distributions showed that larger mean kinetic energy (MKE) and eddy available potential energy (EAPE) were stored in the southern sector of the WLE. While the larger MKE was directly linked with the stronger barotropic currents, the larger EAPE in the southern WLE was formed by baroclinic energy conversions due to a strong density gradient at the front.