海气耦合模式中挟卷参数化的比较

利用美国国家气象中心的综合性大气-海洋资料序列(COADs),采用三种海气耦合模式对海洋挟卷速度的表达方法,计算了实测资料影响下模式挟卷速度的分布特征。利用求得的挟卷速度和由资料序列求得的热量通量分布,采用Niiler—Kraus一维混合层SST预报模式,估计了各种模式扶卷作用对SST的影响。通过比较各个模式的计算结果,提出: 1.挟卷作用是影响SST变化分布的一个重要因子。挟卷速度的分市取决于所采用参数化形式中主要影响因素的分布情况。不考虑平流作用时,模式SST的分布也取决于挟卷参数化形式中主要影响因素的分布。2.由海气耦合不稳定模态的分析得出三个模式对不稳定海气耦合波影响的特征。     

非粘结柔性立管螺旋键缠绕角度敏感性分析

非粘结柔性立管结构形式复杂,层与层之间涉及摩擦、接触等非线性特性,使其数值模拟方面很难同时保证效率与精度。针对柔性立管特殊的结构形式,本文建立不同角度螺旋键的非粘结柔性立管高效分析简化模型,通过Ansys workbench软件分析非粘结柔性立管在内压和拉伸、扭转、弯曲的组合荷载工况下螺旋键缠绕角度对其刚度的影响。结果表明:本文所建立的立管简化分析模型与立管原型的力学性能吻合较好,并且可以提高计算效率;螺旋键缠绕角度越大,立管的扭转刚度越大,拉伸刚度与弯曲刚度越小;内压的施加对拉伸刚度、扭转刚度和弯曲刚度是有利的,内压对于弯曲刚度有较明显的增加,对拉伸刚度和扭转刚度有较小程度的增加。     

关于海气耦合模式气候漂移及敏感性的一点探讨

本文从概念性的大气或海洋系统及海气耦合系统出发,通过简单的理论分析指出:用观测气候场作为边界条件来调试大气或海洋模式并非最佳选择;海气耦合漂移是由两部分组成的,其中一部分源于未耦合模式模拟得到的气候状态与观测气候状态之间的系统误差,而另一部分则源于海气系统的非线性相互作用;海面“通量修正”可以消除因模式气候与观测气候之间的系统误差而引起的那部分漂移,但仍保留了海气系统的大部分非线性相互作用项.本文最后利用中国科学院大气物理研究所发展的海气耦合模式进行了敏感性试验,指出耦合漂移受到模式海洋的垂真发辨的率、海洋温盐扩散方案、海气耦合强度等诸多因素的影响.     

一个用于全球气候变化研究的海洋大气环流耦合模式Ⅰ.模式的形成及性能

本文初步建立了一个可用于模拟几十年到上百年气候变化的全球大气-海洋-海冰耦合的环流模式(CGCM).该模式采用“预估-校正”的月平均距平耦合方案,把一个两层大气环流模式(AGCM)和一个20层海洋环流模式(OGCM)耦合在一起;由于北冰洋尚未被包括在海洋模式内,因此我们在大气模式中嵌入了一个北冰洋薄层热力学海洋模式.130年的积分试验表明,尽管耦合模式的长时期积分中存在着明显的气候漂移现象,但模式模拟的全球气候的空间分布特征基本上是合理的.这为进一步利用该模式模拟全球增强温室效应提供了数值实验依据.     

有限区域大气-海浪耦合模式的建立及海表粗糙度参数化试验

利用LINUX操作系统下的进程通讯(IPC)技术将中尺度大气模式MM5(V3)与第三代海浪模式WW3进行双向耦合,建立考虑大气-海浪相互作用的风浪耦合模式,在耦合模式中引入3种海表粗糙度参数化方案,通过对一次热带气旋过程的模拟,研究大气-海浪相互作用对热带气旋系统的影响及耦合模式对海表粗糙度参数化方案的敏感性。结果表明:LINUX系统下的进程通讯技术可以方便有效地实现大气和海浪模式的双向耦合,模式运行稳定;耦合模式能够较好的模拟热带气旋的发展和演变过程及其影响下海浪场的分布和演变,模拟结果对海表粗糙度参数化方案较敏感;海浪的反馈作用同时影响了海气间的动力和热力作用过程,不同的海表粗糙度参数化方案下,海浪对两种作用过程不同的影响程度决定了其对气旋系统强度的影响。     

波浪-海流-微地形耦合的沉积动力模式建立及应用

沙纹微地形普遍存在于海底,沙纹的消长能改变底部应力进而影响泥沙的运移。以往研究较多侧重于波致沙纹,并已应用于波浪模式的底摩擦计算,而较少考虑波流联合效应产生的沙纹,也未将其应用于综合的水动力模式和沉积物输运模式。本文在POM水动力模式中嵌入新南威尔士大学泥沙模式,通过耦合波流共同作用的微地形模型与波流相互作用底边界层模型,发展了波浪-海流-微地形（沙纹）耦合的沉积动力模式。本文将该模式应用于澳大利亚Jervis湾,针对波主导和波流联合主导沙纹两种类型,分别进行了沙纹发展状态、几何形态的分布及悬浮泥沙的模拟。结果表明:波致沙纹比波流联合作用的沙波具有更大的波高和波长,因此当波主导时沙纹对悬浮泥沙起着关键作用。通过考虑随沙纹变化的粗糙度,相比于以往模式设置均一的粗糙度,该模型能对悬浮物浓度的骤升过程进行更精细的预测。     

用于厄尔尼诺-南方涛动(ENSO)研究的海气耦合模式综述：中间型和混合型模式

厄尔尼诺-南方涛动(El Niño-Southern Oscillation, ENSO)是地球气候系统中最强的年际变率信号, 起源于热带太平洋海气相互作用过程, 并对全球的天气和气候等产生显著的影响。过去几十年来, 广泛、深入而细致的海气相互作用研究致力于发展和改进海气耦合模式以用于ENSO模拟和预测, 各种类型的海气耦合模式应运而生。经过半个多世纪的努力, ENSO数值模式及其应用已经取得了巨大进展, 包括已发展了一些高度理想化的概念(concept)模型来解释ENSO准周期性循环(包括正负反馈机制等); 同时也已发展了几类复杂程度不同的海气耦合模式并用于对ENSO的真实模拟和实时预测等研究, 尤其是已能提前6个月或更长时间对ENSO事件的发生和发展等进行有效的实时预测。其中最为复杂的模式是基于原始方程组的大气环流模式(Atmospheric General Circulation Models, AGCMs)与海洋环流模式(Oceanic General Circulation Models, OGCMs)等所组成的环流型耦合模式(Coupled General Circulation Models, CGCMs), 这类模式变量取为完全变量的形式(如总的海表温度场, 其可以分解为气候态部分和年际异常部分), 还考虑了尽可能详尽的物理过程及其参数化方案。中间型耦合模式(Intermediate Coupled Models, ICMs)是一类介于高度理想化概念模型与复杂的环流型耦合模式之间的简化模式, 其对应的控制方程组采用距平形式, 直接取大气和海洋年际异常场作为预报变量(如海表温度年际异常), 而相应的气候平均态部分则由对应的观测资料来给定; 大气与海洋模式间的耦合采用异常耦合(anomaly coupling)。混合型耦合模式(Hybrid Coupled Models, HCMs)是另一类简化的海气耦合模式, 其中海洋或大气模式有一个分量模式采用了简化的距平类模式(类似于ICMs),而另一个分量模式则采用环流型模式(General Circulation Models, GCMs); 如可采用统计的大气风应力年际异常模式与OGCM间的耦合而构建一种HCM_OGCM,也可采用简化的海洋距平类模式(如ICM中的海洋分量模式)与AGCM间的耦合而构建另一种HCM^AGCM。历史上, ICMs、HCMs和CGCMs等这几类耦合模式都在ENSO理论体系的发展、数值模拟和实时预测等方面都起到了重要作用。本文主要回顾作者与合作者所研发的ICMs和HCMs的构建、特点和应用例子等。     

一种计算斜坡处海流的诊断模式

本文提出一种考虑底斜效应,适合斜坡处海流诊断计算的模式,并且分别计算了日本以南黑潮区部分测线和对马暖流流源区的流速。其结果与实测符合较好。     

基于局地相似理论的蒸发波导计算方案及敏感性试验

蒸发波导是海洋环境中最常出现的异常大气折射结构,具有很高的利用价值.利用局地相似理论代替 Monin-Obukhov 相似理论来确定近地层温度、湿度垂直分布廓线,并借助热带海洋全球大气-海气耦合响应试验(TOGA COARE) 中提出的总体通量算法,对 Local 模式和 New 的推导过程和算法设计进行了详细的阐述和说明.通过比较不同模式的计算结果和敏感性试验发现,在较强逆温和低风速(强稳定层结) 条件下,除 P-J 模式外,基于 Monin-Obukhov 相似理论的蒸发波导模式计算的高度值均出现异常突变;而基于局地相似理论的 Local 模式能够有效抑制这种异常突变的现象,且与加入人为修正的 P-J 模式结果比较接近.New 模式集 Babin 模式与 Local 模式的优点于一身,应该说是目前理论上最为完善的模式.     

一种计算斜坡处海流的诊断模式

本文提出一种考虑底斜效应，适合斜坡处海流诊断计算的模式，并且分别计算了日本以南黑潮区部分测线和对马暖流流源区的流速。其结果与实测符合较好。     

Temporal electronic structure of non-resonant Raman excited virtual state of P-nitroaniline by 514 nm excitation via bond polarisabilities

We have studied the temporal bond polarisabilities of para-nitroaniline from the Raman intensities by the algorithm proposed by Wu et al. in 1987 (Tian B, Wu G, Liu G 1987 J. Chem. Phys. 87 7300). The bond polarisabilities provide much information concerning the electronic structure of the non-resonant Raman excited virtual state. At the initial moment by the 514.5 nm excitation, the tendency of the excited electrons (mapped out by the bond polarisabilities) is to spread to the molecular periphery, and the electronic structure of the Raman virtual state is close to the pseudo-quinonoidic state. When the final stage of relaxation is approached, the bond polarisabilities of those peripheral bonds relax faster than those closer to the molecular core, the phenyl ring. The molecule is in the benzenoidic form as demonstrated by the bond polarisabilities after relaxation.     

Weak-mode identification and time-series reconstruction from high-level noisy measured data of offshore structures

The identification of true weak modes buried in high-level, noisy, measured data from offshore structures is a practical but challenging problem because weak modes are typically eliminated as noise and rarely, yield a discrete time series. This study proposes a weak-mode identification and time-series reconstruction method for offshore structures when high-level noise is present. A theoretical development proposed in this study extends the traditional modal analysis to reconstructing the discrete time series of weak modes, thereby removing its previous limitations to only frequencies, damping ratios and mode shapes. Additionally, a second development proposed in this study makes the reconstructed time series not simply a combination of harmonic components from a Fourier transform but rather complex exponentials; the damping of the test structure is thus estimated with a better accuracy. A third theoretical development avoids variations in the results from different original signals by handling multiple signals simultaneously. The proposed approach primarily includes three steps: (1) estimate the poles and corresponding residues of high-level, noisy, measured data by converting high-order difference equations to first-order difference equations; (2) isolate the poles of weak modes by assigning multiple rough-pole windows, and subsequently extract the corresponding residues based on the row number of the isolated pole vector; and (3) identify and reconstruct the time series of the weak modes of interest in the form of complex exponentials. The most primary advantage of the proposed process in engineering applications is that the pole windows can be easily obtained and assigned from the relationship between the frequencies and their poles. Three numerical examples are studied: the first presents the detailed numerical operation of the proposed method, the second extends the proposed method from managing one signal to managing multiple signals, and the third demonstrates the advantage of the approach compared with traditional methods. The numerical results indicate that the original signals can be decomposed into multiple complex exponentials with representative poles and corresponding residues, and that the new signals representing weak modes could be reconstructed by assigning a range of frequencies in terms of their relations with the poles. To study the performance of the proposed method when applied to offshore structures such as offshore platforms and marine risers, the experimental data from the high mode VIV experiments sponsored by the Norwegian Deepwater Programme (NDP) are used firstly. The results show that two dominant frequencies corresponding to the in-line and cross-flow directions can be identified simultaneously even one mode is very weak compared with the other, and the time series of the weak mode could be reconstructed with a rough frequency window. Then sea-test data of two offshore platforms are used: one was collected from the JZ20-2MUQ offshore platform when it was excited by ice, and the other was collected from the WZ11-4D platform when it was excited by waves. The results further demonstrated that a large model order is required to estimate all poles and residues of the original noisy signals, and that the row number corresponding to a weak mode of the isolated pole matrix could be easily determined via finite element analysis or engineering experiences. Therefore, the proposed approach provides not only modal parameters, such as frequencies and damping ratios of true weak modes buried in high-level noise, but also the discrete time series of the weak mode.     

一种模态参数识别的虚假模态剔除技术

针对背景噪声下结构的模态参数识别结果存在虚假模态的情况,分析虚假模态产生的原因,通过理论推导提出辨别虚假模态的判断条件,根据判断条件剔除识别模态参数中包含的虚假模态。以比例阻尼海洋平台结构数值模型为例,在模拟噪声水平5%的工况下,获得结构脉冲响应信号,对基于奇异值分解定阶消噪后的信号用复指数法进行模态参数识别,对识别的模态频率和阻尼比进行虚假模态剔除。结果表明:根据判断条件可有效剔除所识别参数中的虚假模态。     

Modal analysis of geoacoustic influences on shallow-waterpropagation

Features observed in horizontal wavenumber spectra from low-frequency acoustic measurements taken off the New Jersey Shelf are analyzed through forward modeling using PE models. Environmental models which produce very good agreement with relative intensity data may only match wavenumbers of lower order modes. A prominent feature in contours of modal amplitude versus horizontal wavenumber and depth is the presence of double images of some higher order modes. Two possible causes of this phenomenon are examined. Range dependence in the environment can produce multiple modal images, but the magnitudes of variations required are much larger than those expected at the experiment site. The more likely cause of the double modes is shown to be a duct deep (two acoustic wavelengths) within the sediment. It produces modes with wavenumbers that are close together and with shapes that are very similar in the water column. Comparisons of model predictions with 50-Hz experimental wavenumbers and mode shapes show good agreement, except for amplitudes of some higher order modes. Predictions at 75 Hz typically require more detailed information about the sediment structure in order to accurately match all the features of higher order modes in the data     

Single-mode excitation in the shallow-water acoustic channel usingfeedback control

The shallow-water acoustic channel supports far-field propagation in a discrete set of modes. Ocean experiments have confirmed the modal nature of acoustic propagation, but no experiment has successfully excited only one of the suite of mid-frequency trapped modes propagating in a coastal environment. The ability to excite a single mode would be a powerful tool for investigating shallow-water ocean processes. A feedback control algorithm incorporating elements of adaptive estimation, underwater acoustics, array processing, and control theory to generate a high-fidelity single mode is presented. This approach also yields a cohesive framework for evaluating the feasibility of generating a single mode with given array geometries, noise characteristics, and source power limitations. Simulations and laboratory wave guide experiments indicate the proposed algorithm holds promise for ocean experiments     

The complex band structure for armchair graphene nanoribbons

Using a tight binding transfer matrix method, we calculate the complex band structure of armchair graphene nanoribbons. The real part of the complex band structure calculated by the transfer matrix method fits well with the bulk band structure calculated by a Hermitian matrix. The complex band structure gives extra information on carrier's decay behaviour. The imaginary loop connects the conduction and valence band, and can profoundly affect the characteristics of nanoscale electronic device made with graphene nanoribbons. In this work, the complex band structure calculation includes not only the first nearest neighbour interaction, but also the effects of edge bond relaxation and the third nearest neighbour interaction. The band gap is classified into three classes. Due to the edge bond relaxation and the third nearest neighbour interaction term, it opens a band gap for N=3M-1. The band gap is almost unchanged for N=3M+1, but decreased for N=3M. The maximum imaginary wave vector length provides additional information about the electrical characteristics of graphene nanoribbons, and is also classified into three classes.     

浅海桶形基础平台水平承载力与抗滑稳定分析

平台在海上受风、浪、流、冰等水平载荷作用产生滑移,而土壤粘结力、摩擦力和土抗力是平台的抗滑力,桶基平台应满足抗滑稳定的要求,以保证平台的整体稳定性.介绍了浅海桶形基础平台受水平力作用时桶形基础的破坏模式及其计算方法,并阐述了抗滑稳定计算方法.     

Coastal Disturbance in Sea Level Propagating along the South Coast of Japan and Its Impact on the Kuroshio

The coastal sea level propagating westward along the south coast of Japan and the impact of the disturbance on the generation of the Kuroshio small meander have been examined. The propagation occurs in sea level variations for periods shorter than 10 days and is remarkable for periods of 4–6 days. Characteristics of the 4–6 day component have been studied using the extended empirical orthogonal function (EEOF). The first and second modes of EEOF are almost in-phase throughout the south coast of Japan. The higher four modes of EEOF are significantly excited when the Kuroshio takes the non-large-meander path, and propagate westward with phase speeds of 2.8 m s⁻¹ (third and fourth modes) and 1.6 m s⁻¹ (fifth and sixth modes) in the Kuroshio region west of Mera in the Boso Peninsula. The analysis shows that more than 70% of the small meanders generate in two months after a significant propagating disturbance reaches south of Kyushu when the velocity of the Kuroshio is high. This effect of coastal disturbance is examined by numerical experiments with a 2.5-layer model in which coastal disturbance is excited by vertical displacement of the upper interface. The result is that offshore displacement of the Kuroshio occurs southeast of Kyushu only in the case of significant upward displacement of the interface under the influence of a high Kuroshio velocity. The significant coastal disturbance, which is associated with upward displacement of the density interface, and a high Kuroshio velocity can therefore be important factors in generating small meanders.     

北半球冬季南太平洋海表温度异常的主要模态及其与ENSO的关系

利用新版本的NOAA扩展重构的海温资料,研究了1950/1951年-2008/2009年北半球冬季(12月-翌年2月平均)南太平洋海温异常的主要模态,并探讨了海温主要模态在年际和年代际时间尺度上与ENSO的关系。研究结果表明:北半球冬季南太平洋海温变化主要存在两个模态,第1模态基本反映了南太平洋海温异常南－北的反相变化特征,第2模态基本反映了南太平洋海温异常东－西的反相变化特征。两个模态的年代际周期均比其年际周期显著,其中第1模态的主要年代际周期为11~16 a,第2模态的主要年代际周期为18~24 a。主要模态与ENSO关系的分析表明,第1模态与同期ENSO事件的关系最显著,主要显示其对ENSO的响应,而第2模态对滞后其1 a的ENSO事件存在一定的影响。为进一步研究它们的年际和年代际关系特征,分别对它们进行滤波后的分析表明,在年际时间尺度上,南太平洋海温两个主模态同ENSO事件的关系与滤波前的基本一致;在年代际时间尺度上,第1模态的年代际变化可能对ENSO的年代际变化存在一定的响应,而第2模态可能不仅对前期的ENSO事件存在一定的响应,还可能对后期的ENSO事件存在一定程度的影响作用。     

太平洋中的主要起伏模态

对海洋中起伏运动(heaving)信号的时空分布研究能够帮助我们更好地了解气候系统中的年际和年代际变率。文章通过再分析资料和模式对太平洋区域的heaving主要模态进行了研究。研究结果表明: 太平洋区域主要存在两种heaving模态: 第一模态主要表现为赤道东西两侧的温跃层异常信号反位相; 第二模态表现为赤道区域和副热带区域的温跃层异常信号呈现反位相变化的规律。本文对这两个主要heaving模态所涉及的物理过程进行详细讨论, 结果表明: 东西反位相模态主要是受赤道波动调节的结果; 而经向结构模态则主要是由赤道地区的波动和副热带区域的风应力旋度异常作用共同导致。此外, 我们还讨论了heaving模态可以通过海洋波动以及Ekman输送等过程对海盆尺度的热输送(振幅约为5×10¹⁴W)以及海洋热含量(振幅约为1.5×10²⁰J)的再分配起到了关键的调制作用, 进一步表明heaving模态对全球气候变化有着重要的作用。