The watermass variability and southward shift of the Southern Hemispheremid-depth supergyre

The Southern Hemisphere subtropical supergyre at intermediate depths connects all three ocean basins and plays a significant role in responding and conveying the climate-change-related variations in the glob- al ocean. On the basis of the Simple Ocean Data Assimilation/SODA） ocean reanalysis, the thermohaline variability and southward shift of the mid-depth supergyre are demonstrated. The steric height of the sub- surface relative to 1 500 m （400-1 500 m） from the SODA depicts exactly the flow patterns and variability of the oceanic supergyre. During 1958-2007 the water masses in the gyre interiors become cooler/fresher, with the significant exceptions of the Agulhas Current system and Agulhas leakage. The results also exhibit a pronounced strengthening of the inter-basin connection of the supergyre, and the strongest southward shift, by about 2.5° over the whole period, occurs in the central-south Pacific, which is associated with the changes in the basin-scale wind forcing.     

中国近海以及邻近海域海洋再分析产品(CORA)中季节内变化的评估

A regional reanalysis product—China Ocean Reanalysis(CORA)—has been developed for the China's seas and the adjacent areas. In this study, the intraseasonal variabilities(ISVs) in CORA are assessed by comparing with observations and two other reanalysis products(ECCO2 and SODA). CORA shows a better performance in capturing the intraseasonal sea surface temperatures(SSTs) and the intraseasonal sea surface heights(SSHs) than ECCO2 and SODA do, probably due to its high resolution, stronger response to the intraseasonal forcing in the atmosphere(especially the Madden-Julian Oscillation), and more available regional data for assimilation. But at the subsurface, the ISVs in CORA are likely to be weaker than reality, which is probably attributed to rare observational data for assimilation and weak diapycnal eddy diffusivity in the CORA model. According to the comparison results, CORA is a good choice for the study related to variabilities at the surface, but cares have to be taken for the study focusing on the subsurface processes.     

基于ROMS模式利用集合最优插值同化沿轨海表面高度异常数据的研究

The ensemble optimal interpolation （EnOI） is applied to the regional ocean modeling system （ROMS） with the ability to assimilate the along-track sea level anomaly （TSLA）. This system is tested with an eddy-resolving system of the South China Sea （SCS）. Background errors are derived from a running seasonal ensemble to account for the seasonal variability within the SCS. A fifth-order localization function with a 250 km localization radius is chosen to reduce the negative effects of sampling errors. The data assimilation system is tested from January 2004 to December 2006. The results show that the root mean square deviation （RMSD） of the sea level anomaly decreased from 10.57 to 6.70 cm, which represents a 36.6% reduction of error. The data assimilation reduces error for temperature within the upper 800 m and for salinity within the upper 200 m, although error degrades slightly at deeper depths. Surface currents are in better agreement with trajectories of surface drifters after data assimilation. The variance of sea level improves significantly in terms of both the amplitude and position of the strong and weak variance regions after assimilating TSLA. Results with AGE error （AGE） perform better than no AGE error （NoAGE） when considering the improvements of the temperature and the salinity. Furthermore, reasons for the extremely strong variability in the northern SCS in high resolution models are investigated. The results demonstrate that the strong variability of sea level in the high resolution model is caused by an extremely strong Kuroshio intrusion. Therefore, it is demonstrated that it is necessary to assimilate the TSLA in order to better simulate the SCS with high resolution models.     

Wave-Induced Loads on Very Large FPSOs at Restricted Water Depth

The effects of water depth on the wave-induced vertical bending moment and shearing force on a very large FPSO are studied by experiments and computations for regular and irregular waves. The restricted water depth composite Green function is employed to develop a program for the computation of the hydrodynamic coefficients of the very large FPSO at shallow water. A thrce-segment model with 1 : 100 scale is tested in the State Key Laboratory of Ocean Engineering at Shanghai Jiao Tong University for the verification of the nmnerical method. The experimental and computational results show that the water depth has a substantial effect on wave-induced loads. The wave-induced vertical loads increase with the decrease of water depth for shallow water. Especially, for ultra-shallow water these loads increase very evidently with the decrease of water depth. The long-term prediction values of wave-induced vertical loads increase with the decrease of the ratio of water depth to draught. The long-term prediction values of wave-induced vertical loads are about 8% larger than those for deep water when the ratio of water depth to draught is 3.0. However, water depth hardly affects the longterm prediction values of wave-induced loads when the ratio of water depth to draught is larger than 5.0.     

利用热带太平洋TAO浮标阵列评估全球高分辨率浪-潮-流耦合模式同化结果

In order to evaluate the assimilation results from a global high resolution ocean model, the buoy observations from tropical atmosphere ocean(TAO) during August 2014 to July 2015 are employed. The horizontal resolution of wave-tide-circulation coupled ocean model developed by The First Institute of Oceanography(FIOCOM model) is 0.1°×0.1°, and ensemble adjustment Kalman filter is used to assimilate the sea surface temperature(SST), sea level anomaly(SLA) and Argo temperature/salinity profiles. The simulation results with and without data assimilation are examined. First, the overall statistic errors of model results are analyzed. The scatter diagrams of model simulations versus observations and corresponding error probability density distribution show that the errors of all the observed variables, including the temperature, isotherm depth of 20°C(D20), salinity and two horizontal component of velocity are reduced to some extent with a maximum improvement of 54% after assimilation. Second, time-averaged variables are used to investigate the horizontal and vertical structures of the model results. Owing to the data assimilation, the biases of the time-averaged distribution are reduced more than70% for the temperature and D20 especially in the eastern Pacific. The obvious improvement of D20 which represents the upper mixed layer depth indicates that the structure of the temperature after the data assimilation becomes more close to the reality and the vertical structure of the upper ocean becomes more reasonable. At last,the physical processes of time series are compared with observations. The time evolution processes of all variables after the data assimilation are more consistent with the observations. The temperature bias and RMSE of D20 are reduced by 76% and 56% respectively with the data assimilation. More events during this period are also reproduced after the data assimilation. Under the condition of strong 2014/2016 El Ni?o, the Equatorial Undercurrent(EUC) from the TAO is gradually increased during August to November in 2014, and followed by a decreasing process. Since the improvement of the structure in the upper ocean, these events of the EUC can be clearly found in the assimilation results. In conclusion, the data assimilation in this global high resolution model has successfully reduced the model biases and improved the structures of the upper ocean, and the physical processes in reality can be well produced.     

A New Wetting and Drying Method for Moving Boundary in Shallow Water Flow Models

To deal with the moving boundary hydrodynamic problems of the tidal flats in shallow water flow models,a new wetting and drying (WD) method is proposed.In the new method,a "predicted water depth" is evaluated explicitly based on the simplified shallow water equations and used to determine the status (wet or dry) together with the direction of flow.Compared with previous WD method,besides the water elevation,more factors,such as the flow velocity and the surface shear stress,are taken into account in the new method to determine the moving boundary.In addition,a formula is deduced to determine the threshold,as critical water depth,which needs to be preset before simulations.The new WD method is tested with five cases including three 1D ones and two 2D ones.The results show that the new WD method can simulate the wetting and drying process,in both typical and practical cases,with smooth manner and achieves effective estimation of the retention volume at shallow water body.     

A New Wetting and Drying Method for Moving Boundary in Shallow Water Flow Models

To deal with the moving boundary hydrodynamic problems of the tidal flats in shallow water flow models, a new wetting and drying (WD) method is proposed. In the new method, a “predicted water depth” is evaluated explicitly based on the simplified shallow water equations and used to determine the status (wet or dry) together with the direction of flow. Compared with previous WD method, besides the water elevation, more factors, such as the flow velocity and the surface shear stress, are taken into account in the new method to determine the moving boundary. In addition, a formula is deduced to determine the threshold, as critical water depth, which needs to be preset before simulations. The new WD method is tested with five cases including three 1D ones and two 2D ones. The results show that the new WD method can simulate the wetting and drying process, in both typical and practical cases, with smooth manner and achieves effective estimation of the retention volume at shallow water body.     

A comparison of two global ocean-ice coupled models with different horizontal resolutions

A global eddy-permitting ocean-ice coupled model with a horizontal resolution of 0.25 by 0.25 is established on the basis of Modular Ocean Model version 4 (MOM4) and Sea Ice Simulator (SIS). Simulation results are compared with those of an intermediate resolution ocean-ice coupled model with a horizontal resolution of about 1 by 1 . The results show that the simulated ocean temperature, ocean current and sea ice concentration from the eddy-permitting model are better than those from the intermediate resolution model. However, both the two models have the common problem of ocean general circulation models (OGCMs) that the majority of the simulated summer sea surface temperature (SST) is too warm while the majority of the simulated subsurface summer temperature is too cold. Further numerical experiments show that this problem can be alleviated by incorporating the non-breaking surface wave-induced vertical mixing into the vertical mixing scheme for both eddy-permitting and intermediate resolution models.     

Direct measurements of surface and mid-depth circulationin the Shikoku Basin by Argo profiling floats

The circulation in the Shikoku Basin plays a very important role in the pathway of the Kuroshio and the water exchange in the subtropical gyre in the North Pacific Ocean. The Argo profiling floats deployed in the Shikoku Basin are used to study the circulations and water masses in the basin. The trajectories and parking depth velocity fields derived from all Argo floats show an anticyclonic circulation at 2 000 m in the Shikoku Basin. There are inhanced eddy activities in the Shikoku Basin, which have large influence on the Shikoku Basin circulation patterns. The characteristics of temperature-salinity curves indicate that there are North Pacific Ocean tropical water (NPTW), North Pacific Ocean subtropical mode water (NPSTMW) and North Pacific Ocean intermediate water (NPIW) in the Shikoku Basin. The NPTW is only exists south of 32°N. In the middle part of the basin, which is 28°~31°N,133°~135°E, there is a confluence region. Water masses coming from the Kuroshio mix with the water in the Shikoku Basin.     

Foraminifera in surface sediments of the Bering and Chukchi Seas and their sedimentary environment

Based on a quantitative analysis of foraminifera in 39 surface samples of the Bering andChukchi Seas, the nearly absence of planktonic foraminifera in the surface sediments can be related to the low surface primary productivity and strong carbonate dissolution in the study area. It has been revealed that the surface primary productivity, and carbonate dissolution and properties of water masses related to the water depth mainly control the distribution of benthic foraminifera. The shelf of the Chukchi Sea is dominated by the Elphidium spp. assemblage and Nonionella robusta assemblage with low foraminiferal abundance and diversity, which is controlled by the coastal water mass of the Arctic Ocean. The slope of the Bering Sea is dominated by the Uvigerina peregrina - Globobulimina affinis assemblage with abundant N. robusta, and relatively high foraminiferal abundance and diversity, which is controlled by the intermediate and deep water masses of the Pacific Ocean. However, the Bering Sea has relatively sha     

CMIP6模式对北冰洋海洋热含量的模拟能力评估

本文利用PHC、ECCO2、SODA、GECCO3和CMIP6资料,分析了北冰洋热含量的水平分布特征、季节变化和长期变化趋势等,评估了CMIP6模式对北冰洋海洋热含量的模拟能力。研究发现,北冰洋海洋热含量表现出明显的季节变化:热含量在4月份最低,9月份最高;在历史情形下(1850?2014年),相较观测和再分析资料,CMIP6多模式集合平均(MME)的上层500 m热含量在格陵兰海偏暖,在挪威海、巴伦支海和欧亚海盆偏冷,MME的全水深热含量在北冰洋几乎所有区域均偏暖,在格陵兰海偏差最大;CMIP6模式对北冰洋温度剖面模拟偏差较大,MME平均温度在1 000 m以深均高于观测和再分析资料。在未来情形下(2015?2100年),MME表现出明显的北冰洋增暖情形,但绝大多数中国模式没有表现出明显的增暖情形。中国模式中,BCC-CSM2-MR和BCC-ESM1对北冰洋年平均热含量的模拟较差,CIESM对热含量季节和年代际变化模拟较差,FIO-ESM-2-0对北冰洋上层500 m年平均热含量及热含量季节和年代际变化的模拟都比较好。     

基于PHC3.0数据集的北冰洋水团分析

基于PHC3.0极地科学中心水文气候数据集(简称PHC3.0数据集)的温度和盐度资料,使用聚类分析和Bayes判别分析的方法,对北纬70°以北海域的水团结构进行了分析,在北冰洋区域划分出4个水团:北冰洋表层水(ASW)、大西洋中层水(AIW)、太平洋水(PW)和北冰洋深层水(ADW)。北冰洋表层水(ASW)遍布于欧亚海盆和加拿大海盆,以低温低盐为特征。大西洋中层水(AIW)位于约200～900m深度,在北冰洋环极边界流的作用下,其影响可达到加拿大海盆。太平洋水(PW)受经白令海峡进入北冰洋的海水影响,相对高温低盐,夏季时影响显著。北冰洋深层水(ADW)在海盆中相当均匀,几乎没有季节变化,盐度约在34.95psu,温度在加拿大海盆约为-0.3℃,欧亚海盆约为-0.7℃。     

北极冰海耦合模式对两种不同大气再分析资料响应的分析

本文中我们比较了Climate Forecast System Reanalysis(CFSR)高分辨率的再分析数据集和低分辨率的Japanese 25-year Reanalysis Project(JRA25)再分析数据集在向下短波辐射、向下长波辐射、10m风场、近地面气温、降水、湿度上的不同,发现二者差异最大的为降水数据,其次为向下短波辐射数据、向下长波辐射数据。用这两个数据集驱动同一冰海耦合模式,CFSR强迫的海冰、北冰洋中层水和加拿大海盆温盐结构与实测相比有很大差距,等密度面上的地转流速在加拿大海盆和欧亚海盆比JRA25强迫的结果高20%,同时等密度面的深度偏深、位温偏高,在弗拉姆海峡的流通量也比海洋再分析数据Simple Ocean Data Assimilation(SODA)偏多。CFSR的向下辐射数据更加接近实测,采用此数据的敏感性实验模拟结果与实测符合的更好。对于海冰的模拟,云量起着至关重要的作用,降水带来的淡水通量通过影响大西洋入流水携带的热量进而影响到冰区。此外,CFSR过量的降水也是二者对于北冰洋温盐结构、弗拉姆海峡流通量以及地转流强度模拟产生偏差的主要原因。尽管风场的分辨率不同,在海盆尺度上对于海冰和海水温盐结构的影响并不大。     

An assessment of Arctic cloud water paths in atmospheric reanalyses

The role of Arctic clouds in the recent rapid Arctic warming has attracted much attention. However, Arctic cloud water paths(CWPs) from reanalysis datasets have not been well evaluated. This study evaluated the CWPs as well as LWPs(cloud liquid water paths) and IWPs(cloud ice water paths) from five reanalysis datasets(MERRA-2,MERRA, ERA-Interim, JRA-55, and ERA5) against the COSP(Cloud Feedback Model Intercomparison Project Observations Simulator Package) output for MODIS from the MERRA-2 CSP(COSP satellite simulator) collection(defined as M2 Modis in short). Averaged over 1980–2015 and over the Arctic region(north of 60°N), the mean CWPs of these five datasets range from 49.5 g/m~2(MERRA) to 82.7 g/m~2(ERA-Interim), much smaller than that from M2 Modis(140.0 g/m~2). However, the spatial distributions of CWPs, show similar patterns among these reanalyses, with relatively small values over Greenland and large values over the North Atlantic. Consistent with M2 Modis, these reanalyses show larger LWPs than IWPs, except for ERA-Interim. However, MERRA-2 and MERRA underestimate the ratio of IWPs to CWPs over the entire Arctic, while ERA-Interim and JRA-55 overestimate this ratio. ERA5 shows the best performance in terms of the ratio of IWPs to CWPs. All datasets exhibit larger CWPs and LWPs in summer than in winter. For M2 Modis, IWPs hold seasonal variation similar with LWPs over the land but opposite over the ocean. Following the Arctic warming, the trends in LWPs and IWPs during 1980～2015 show that LWPs increase and IWPs decrease across all datasets, although not statistically significant. Correlation analysis suggests that all datasets have similar interannual variability. The study further found that the inclusion of re-evaporation processes increases the humidity in the atmosphere over the land and that a more realistic liquid/ice phase can be obtained by independently treating the liquid and ice water contents.     

2016年夏季北冰洋浮冰站近地层辐射和湍流通量观测

2016年8月7－14日中国第七次北极科学考察期间,在83°N附近设立的长期浮冰站开展了辐射和湍流通量观测研究。结果表明,观测期间反照率变化范围为0.64~0.92,平均反照率为0.78;基于现场观测数据评估了PW79、HIRHAM、ARCSYM和CCSM3 4种不同复杂度的反照率参数化方案在天气尺度的表现,最为复杂的CCSM3结果优于其他参数化方案,但不能体现降雪条件下的反照率快速增长。浮冰区冰雪面平均净辐射为18.10 W/m2,平均感热通量为1.73 W/m2,平均潜热通量为5.55 W/m2,海冰表面消融率为（0.30±0.22） cm/d,表明此时北冰洋浮冰正处于快速消融期。冰面的平均动量通量为0.098（kg·m/s）/（m2·s）,动量通量与风速有很好的对应关系,相关系数达0.80。     

Labrador Sea convection and subpolar North Atlantic Deep Water export in the SODA assimilation model

Labrador Sea convection was most intense and reached the greatest depths in the early 1990s, followed by weaker, shallower, and more variable convection after 1995. The Simple Ocean Data Assimilation (SODA) version 2.0.2/2.0.4 assimilation model is used to explore convective activity in the North Atlantic Ocean for the period from 1992 to 2007. Hydrographic conditions, which are relatively well observed during this period, are used to compare modeled and observed winter mixed-layer depths and water mass anomalies in relation to Deep Western Boundary Current transports and meridional overturning circulation (MOC) changes at the exit of the subpolar basin. The assimilation differs markedly from local observations in the March mixed-layer depth, which represents deep convection and water mass transformation. However, mean MOC rates at the exit of the subpolar gyre, forced by stratification in the mid-latitudes, are similar to estimates based on observations and show no significant decrease during the 1992–2007 period. SODA reproduces the deep Labrador Sea Water formation in the western North Atlantic without any clear indication of significant formation in the Irminger Sea while the lighter upper Labrador Sea Water density range is reached in the Irminger Sea in the 1990s, in agreement with existing assumptions of deep convection in the Irminger Sea and also supported by computed lag correlations with the Labrador Sea. Deep Water transformation mainly takes place in the eastern North Atlantic. The introduction of CFC-11 into the SODA model as a tracer reproduces the mean and multiyear variations of observed distributions.     

基于神经网络的全球三维温盐场重构技术研究

文章利用果蝇优化广义回归神经网络算法FOAGRNN (fruit fly optimization algorithm, FOA; generalized regression neural network, GRNN)对SODA (simple ocean data assimilation)再分析数据进行训练, 构建海表温度、盐度、海面高度与次表层温盐场之间的投影关系模型, 并在全球范围使用SODA和卫星遥感数据评估了模型的应用性能。首先, 利用独立的2016年SODA海表数据作为模型输入进行理想重构试验, 结果显示全球重构温、盐平均均方根误差(MRMSE)分别为0.36℃和0.08‰, 与世界海洋图集WOA13资料相比减小约50%和60%。然后, 利用卫星观测的海表信息作为模型输入进行实际应用试验, 并与Argo观测剖面进行比较评估。试验结果表明, 重构模型能有效表征海水温、盐特征, 其中重构温、盐MRMSE分别为0.79℃和0.16‰, 相比WOA气候态减小27%和11%。误差的垂向分布显示, 重构温度RMSE从海表向下迅速增大, 至100m达到峰值1.35℃, 而后又迅速回落,至250m处为0.81℃, 跃层往下不断减小; 重构盐度RMSE基本随深度增大而减小, 误差峰值位于25m附近, 约为0.25‰。此外, Argo浮标跟踪分析和区域水团统计结果也表明模型能够较好地刻画海洋三维温盐场的内部结构特征。     

西太平洋暖池热状况变化特征及其东传过程

应用JEDAC海温和SODA再分析资料,分析研究了西太平洋暖池(WPWP)的热状况变化特征,探讨了WPWP的热含量变化与不同深度海温变化,特别是与表层和次表层海温场之间的内在联系,分析研究了WPWP异常海温东传过程中与温跃层变化之间的关系,最后讨论了赤道潜流在异常海温东传过程的作用。结果表明,WPWP区域的热含量变化与次表层的海温变化具有相同变化趋势。WPWP区域的热含量变化与次表层海温达到(超过)99.9%信度检验的相关出现在次表层海温变化超前热含量变化的12个月至次表层海温变化滞后热含量变化12个月,最佳相关出现在同期,相关系数r=0.92(样本n=552)。WPWP区域的热含量变化与SST达到99.9%信度检验的相关只出现在热含量滞后SST5个月至同期。这一结果表明,WPWP区域的热含量变化主要体现了次表层海温场的变化特征。由此认为,WPWP区域热含量变化(次表层海温场)对气候变化的贡献主要体现在东向传播对ENSO的贡献上。由WPWP异常海温沿赤道东传过程和温跃层变化之间的关系可知,温跃层变化滞后于WPWP区域东传的异常海温,在WPWP区域异常海温东传过程中,赤道潜流起了重要作用。     

巴芬湾西部和格陵兰岛东部大气经向热量输送协同变化对夏季北极海冰的影响

利用欧洲中心气候再分析资料和美国国家冰雪数据中心北极海冰面积资料,分析了夏季北极海冰面积与前期大气经向热量输送年际变化的联系。结果表明：6月北半球中高纬大气的经向热量输送以瞬变热量形式为主,其中巴芬湾西部（B区）和格陵兰岛东部（G区）是瞬变热量向极区传输的两个通道,二者之间存在反位相的协同变化,且这种协同变化与夏季北极海冰面积变化密切相关。可能的机制为：6月,AD、AO和NAO三种北极大气环流型能够引起巴芬湾西部和格陵兰岛东部瞬变热量输送的协同变化,这种协同变化通过涡旋动力作用激发夏季极区大气表现为AD异常,同时影响途经区域的气温,从而通过热动力作用影响夏季北极海冰。将向极区输送的热量称为暖输送,从极区输出的热量为冷输送,则上述两个区域的瞬变热量协同输送可分为三种情况：B暖G冷、B冷G暖、B和G均冷,而B和G均暖的情况十分罕见。当B区向极区输入、G区输出热量时,有利于太平洋扇区和喀拉海的海冰偏少;当G区输入、B区输出热量时,利于喀拉海和拉普捷夫海海冰偏少;当B区和G区均输出热量时,利于波佛特海南部、喀拉海和拉普捷夫海海冰偏多,反之则相反。