Ekman模型时间变化风应力系数的伴随参数估计

本文基于时间分布参数设置,利用伴随同化方法,反演了Ekman模型中随时间变化的风应力拖曳系数,并在孪生实验和实际实验中对该方法进行了验证。在孪生实验中,研究了参数反演结果对不同影响因素的响应,包括:风速分布、风应力系数分布、风应力系数初始猜测值、风应力系数独立变量个数、观测数据误差和观测的深度。孪生实验结果验证了伴随同化方法反演Ekman模型中时变风应力系数的有效性,具体包括如下五个方面结论:1)不同风速分布下均能成功反演出不同风应力拖曳系数分布; 2)反演结果对初始猜测值较为敏感,风应力系数初始猜测值越接近给定值,反演结果越好;3)风应力系数独立点个数的选取会显著影响反演结果,合理的选择有利于提高反演效率及减小观测数据误差;4)观测误差能够影响反演结果,观测数据误差在20%以下时能取得合理的反演结果; 5)反演结果对观测数据的表层和次表层流速更为敏感,这是由Ekman流的物理性质决定的。实际实验,利用百慕大锚系试验平台的风速和流速数据,去除周期性潮流和地转流成分后得到Ekman流成分,并作为观测输入到该同化模型,反演出了适用于该区域和该时段的随时间变化的风应力系数。通过比较模拟流速和观测流速,证明利用伴随同化方法能从实测数据中反演出合理的时变风应力系数,对于海洋模型风应力系数的确定是一项有益的尝试。     

中国物理海洋学研究70年:发展历程、学术成就概览

本文概略评述新中国成立70年来物理海洋学各分支研究领域的发展历程和若干学术成就。中国物理海洋学研究起步于海浪、潮汐、近海环流与水团,以及以风暴潮为主的海洋气象灾害的研究。随着国力的增强,研究领域不断拓展,涌现了大量具有广泛影响力的研究成果,其中包括:提出了被国际广泛采用的“普遍风浪谱”和“涌浪谱”,发展了第三代海浪数值模式;提出了“准调和分析方法”和“潮汐潮流永久预报”等潮汐潮流的分析和预报方法;发现并命名了“棉兰老潜流”,揭示了东海黑潮的多核结构及其多尺度变异机理等,系统描述了太平洋西边界流系;提出了印度尼西亚贯穿流的南海分支(或称南海贯穿流);不断完善了中国近海陆架环流系统,在南海环流、黑潮及其分支、台湾暖流、闽浙沿岸流、黄海冷水团环流、黄海暖流、渤海环流,以及陆架波方面均取得了深刻的认识;从大气桥和海洋桥两个方面对太平洋–印度洋–大西洋洋际相互作用进行了系统的总结;发展了浅海水团的研究方法,基本摸清了中国近海水团的分布和消长特征与机制,在大洋和极地水团分布及运动研究方面也做出了重要贡献;阐明了南海中尺度涡的宏观特征和生成机制,揭示了中尺度涡的三维结构,定量评估了其全球物质与能量输运能力;基本摸清了中国近海海洋锋的空间分布和季节变化特征,提出了地形、正压不稳定和斜压不稳定等锋面动力学机制;构建了“南海内波潜标观测网”,实现了对内波生成–演变–消亡全过程机理的系统认识;发展了湍流的剪切不稳定理论,提出了海流“边缘不稳定”的概念,开发了海洋湍流模式,提出了湍流混合参数化的新方法等;在海洋内部混合机制和能量来源方面取得了新的认识,并阐述了混合对海洋深层环流、营养物质输运等过程的影响;研发了全球浪–潮–流耦合模式,推出一系列海洋与气候模式;发展了可同化主要海洋观测数据的海洋数据同化系统和用于ENSO预报的耦合同化系统;建立了达到国际水准的非地转(水槽/水池)和地转(旋转平台)物理模 型实验平台;发展了ENSO预报的误差分析方法,建立了海洋和气候系统年代际变化的理论体系,揭示了中深层海洋对全球气候变化的响应;初步建成了中国近海海洋观测网;持续开展南北极调查研究;建立了台风、风暴潮、巨浪和海啸的业务化预报系统,为中国气象减灾提供保障;突破了国外的海洋技术封锁,研发了万米水深的深水水听器和海洋光学特性系列测量仪器;建立了溢油、危险化学品漂移扩散等预测模型,为伴随海洋资源开发所带来的风险事故的应急处理和预警预报提供科学支撑。文中引用的大量学术成果文献(每位第一作者优选不超过3篇)显示,经过70年的发展,中国物理海洋学研究培养了一支实力雄厚的科研队伍,这是最宝贵的成果。这支队伍必将成为中国物理海洋学研究攀登新高峰的主力军。     

X波段雷达海面回波谱宽特征研究

海面电磁回波频谱宽度与海浪波高密切相关,可应用频谱宽度进行海浪有效波高反演。本文应用线性滤波法仿真出了海表散射面元在雷达视向上的投影速度,建立了回波谱宽模型,分析了雷达空间分辨率、回波时间序列长度及海洋环境参数等因素对频谱宽度的影响,同时还针对如何在实际观测过程中选择回波时间序列长度、观测方位角等参数进行了讨论。最后还将理论结果与CSIR-X波段雷达实测数据谱宽估计结果进行了比较。结果表明,剔除雷达噪声以及频率泄露的影响后,基于高斯分布标准偏差的谱宽估计方法所得结果与理论结果吻合很好,这从而证明了理论结果的可靠性。本文所得结果对海浪有效波高反演具有一定参考价值。     

Drainage conditions around monopiles in sand

Large diameter monopiles are typical foundation solutions for offshore wind turbines. In design of the monopile foundations in sand, it is necessary to understand the drainage conditions of the foundation soil under the design loading conditions as the soil performance (strength and stiffness) is highly dependent on the drainage conditions. This paper presents a numerical investigation into this issue, with a purpose to develop a simple design criterion for assessing the soil drainage conditions around a monopile in sand. It is found that for typical monopile foundations in sand, the drainage condition during a single load cycle is generally expected to be undrained. However, the current state-of-practice uses p-y springs derived for drained soil responses for monopile design. The impact of this discrepancy on monopile foundation design was evaluated and found to be insignificant due to the relatively low level of loading as compared to the capacity of the soil.     

Obtaining the distribution of quiescent periods directly from the power spectral densities of Sea waves

There is a growing practical interest in the ability to increase the sea states at which marine operations can be safely undertaken by exploiting the quiescent periods that are well known to exist under a wide range of sea conditions. While the actual prediction of quiescent periods at sea for the control of operations is a deterministic process, the long term planning of future maritime tasks that rely on these quiescent periods is a statistical process involving the anticipated quiescence properties of the forecasted sea conditions in the geographical region of interest. It is in principle possible to obtain such data in tabular form either large scale simulation or from field data. However, such simulations are computationally intensive and libraries of appropriate field data are not common. Thus, it is clearly attractive to develop techniques that exploit standard wave spectral models for describing the quiescence statistics directly from such spectra. The present study focuses upon such techniques and is a first step towards the production of a computationally low-cost quiescence prediction tool and compares its efficacy against simulations. Two significant properties emerge for a large class of wave spectral models that encompasses the ubiquitous Neumann and Pierson Moskowitz or Bretschneider forms. Firstly, the auto-correlation function of the wave profile that are required to produce the quiescence property can be obtained analytically in terms of standard special functions. This considerably reduces the computational cost making desktop computer-based planning tools a reality. Secondly, for each class of these parametric spectra, the probability of a given number of consecutive wave heights (normalised to the significant wave heights) less than some critical value is in fact independent of absolute wave height. Thus, for a broad class of practically interesting wave spectra all that is required to obtain the statistical distribution of the quiescent periods is simple rescaling.     

Hydrodynamic forces on a partially submerged cylinder at high Reynolds number in a steady flow

The hydrodynamic forces on the stationary partially submerged cylinder are investigated through towing test with Reynolds number ranging from 5 × 10⁴ to 9 × 10⁵. Three test groups of partially submerged cylinders with submerged depths of 0.25 D, 0.50 D, and 0.75 D and one validation group of fully submerged cylinders are conducted. During the experiments, the hydrodynamic forces on the cylinders are measured using force sensors. The test results show a considerable difference in the hydrodynamic coefficients for the partially submerged cylinders versus the fully submerged cylinders. A significant mean downward lift force is first observed for the partially submerged cylinders in a steady flow. The maximum of the mean lift coefficients can reach 1.5. Two distinct features are observed due to the effects of overtopping: random distributions in the mean drag coefficients and a clear quadratic relationship between the mean lift coefficients and the Froude number appear in the non-overtopping region. However, the novel phenomenon of a good linear relationship with the Froude number for the mean hydrodynamic coefficients is clearly shown in the overtopping region. In addition, fluctuating hydrodynamic coefficients are further proposed and investigated. These results are helpful to have a better understanding of the problem and to improve related structural designs.     

透明胞外聚合颗粒物碳输运新途径

目前大家普遍认为,透明胞外聚合颗粒物(Transparent Exopolymer Particles,TEP)因其独特的凝聚效应导致碳通量向海底输出。但是,近几年的研究表明TEP不仅影响了碳沉降途径,其本身能够悬浮甚至向海水表层迁移,导致其在海洋微表层(Surface Micro-layer, SML)积累,最终显著影响海洋表层碳通量。TEP和其他颗粒物聚集形成凝聚物后,其运动趋势则由凝聚物中TEP的含量占比,即最终颗粒物密度所决定。一个新的值得注意的现象是,密度低的TEP通过与其他微粒聚合形成表面活性物质(Surface-active Substances, SAS),会在海洋–大气界面形成薄膜,显著影响海–气CO₂交换通量,甚至对全球气候变化造成影响。     

Sea Level Prediction in the Yellow Sea From Satellite Altimetry With a Combined Least Squares-Neural Network Approach

Accessible high-quality observation datasets and proper modeling process are critically required to accurately predict sea level rise in coastal areas. This study focuses on developing and validating a combined least squares-neural network approach applicable to the short-term prediction of sea level variations in the Yellow Sea, where the periodic terms and linear trend of sea level change are fitted and extrapolated using the least squares model, while the prediction of the residual terms is performed by several different types of artificial neural networks. The input and output data used are the sea level anomalies (SLA) time series in the Yellow Sea from 1993 to 2016 derived from ERS-1/2, Topex/Poseidon, Jason-1/2, and Envisat satellite altimetry missions. Tests of different neural network architectures and learning algorithms are performed to assess their applicability for predicting the residuals of SLA time series. Different neural networks satisfactorily provide reliable results and the root mean square errors of the predictions from the proposed combined approach are less than 2?cm and correlation coefficients between the observed and predicted SLA are up to 0.87. Results prove the reliability of the combined least squares-neural network approach on the short-term prediction of sea level variability close to the coast.     

一维综合孔径微波辐射计遥感海面温度的敏感性分析

一维综合孔径微波辐射计能够有效提高观测的空间分辨率,其观测入射角通常在0°~55°范围内变化。为了开发适用于一维综合孔径微波辐射计的海面温度反演算法,需要评估其观测亮温对海洋大气环境要素的敏感性。利用海面发射率模型和大气辐射传输模型,构建了适用于一维综合孔径微波辐射计的微波海洋大气辐射传输模式,研究了C波段垂直和水平极化微波辐射亮温在不同入射角下对海洋大气环境要素的敏感性变化情况,并定量计算了相应的敏感系数。结果表明:垂直和水平极化亮温对海洋大气环境要素的敏感性表现出不同的特性。随着入射角的增大,垂直极化亮温对海面温度的敏感性增强,对海面风场的敏感性相对减弱;水平极化亮温则相反。由大气水汽含量和云液态水含量误差引入的垂直和水平极化亮温误差随入射角增大而增大,但是,即使在55°的大入射角下垂直和水平极化亮温误差仍小于0.12 K。对于海面温度反演精度优于1 K的要求,一维综合孔径微波辐射计的测温精度需优于0.6 K。研究结果对于一维综合孔径微波辐射计海面温度反演算法的研究和载荷设计具有一定的理论指导意义。     

A New Wavelet Method for Identification of Eddies and Assessment of Incidents on Islands of the Eastern Caribbean

Using two dimensional continuous wavelet transforms, a novel method for identification of mesoscale eddies is presented to facilitate extraction of characteristics for area, amplitude, type, and location from maps of sea level anomalies. In comparison with the previously established growing method for eddy identification, it is found that the wavelet method identifies more than twice the number of eddies and is particularly better at resolving small eddies down to the 0.25 degree resolution of the data. Such research into eddy identification and tracking is significant to the assessment of eddies with potential to impact on coastlines of small islands. The method is applied to the identification of eddies on tracks towards islands of the Eastern Caribbean over 23?years. Spatial and temporal variation in rate of occurrence and magnitude is established. For Barbados there is an average of 9 anticyclonic incidents a year with maximum amplitude of typically 0.22?m in the dry seasons and 0.16?m in the wet seasons. Seasonal variation is reversed for the other islands with twice the number of anticyclonic incidents having maximum amplitudes of about 0.20?m annually.