Kalman滤波风暴潮数值预报四维同化模式研究进展

本文首先介绍了Kalman滤波在风暴潮数值预报中的应用，特别介绍了近年来国际上发展的一些在实际中可行的次优化Kalman滤波算法。并通过一个稳态Kalman滤波风暴潮数值预报模式的实例表明，使用资料同化可以明显改进风暴潮后报结果；资料同化能够提供更为合理的预报初始场，对风暴潮的短期预报有较明显的改进。一旦没有资料同化到模式中去，预报结果很快接近确定性模式。     

变分伴随数据同化方法在断面海温数值计算中的应用研究

以二维断面海温分布模型为例,利用海温实际观测数据,将变分伴随方法应用于断面海温初始场的优化。讨论了变分伴随方法的基本思想,分别从模型方程的连续和离散形式出发推导伴随模型系统,并对这两种途径建立的伴随系统之间的相互关系进行了分析。数值试验的结果表明了变分伴随数据同化方法在海温数值计算和数值预报业务中的良好的应用前景。     

一个稳态Kalman滤波风暴潮数值预报模式

利用Kalman滤波资料同化技术将海洋站水位观测资料融入二维线性风暴潮模式中,研制具有资料同化能力的风暴潮预报模式,改进风暴潮模式计算结果.通过在风暴潮模式的动量方程中加入模式噪声项来修正模式本身和气象强迫力的不确定性.确定性模式的输出通过带有观测噪声的观测方程与可利用的海洋站的潮位观测资料联系起来.假定初始的模式噪声和观测噪声满足均值为0的高斯分布,用迭代法得到计算区域的状态向量的稳态Kalman滤波,进而得到风暴潮模式输出的最优线性校正结果.利用这种资料同化技术,对1956年发生在东海的一次强风暴潮过程进行了后报试验,结果表明,该同化方法对短期风暴潮水位后(预)报有一定的改进.     

MASNUM-WAM海浪模式集合Kalman滤波同化研究——Ⅱ.集合样本对同化效果的影响

背景误差相关结构的确定是影响海浪同化效果的关键因素之一。集合Kalman滤波是一种较为成熟的同化方法,其可以对背景误差进行实时更新和动态估计,现已广泛应用于海洋和大气领域的研究。本文基于MASNUM-WAM海浪模式,分别采用静态样本集合Kalman滤波和EAKF方法,针对2014年全球海域开展海浪数据同化实验,同化资料为Jason-2卫星高度计数据,利用Saral卫星高度计资料对同化实验结果进行检验。结果表明,两组同化方案均有效提高了海浪模式的模拟水平,EAKF方案在风场变化较大的西风带区域表现显著优于静态样本集合Kalman滤波方案,但总体上两者相差不大。综合考虑计算成本和同化效果,静态样本集合Kalman滤波方案更适用于海浪业务化预报。     

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

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

基于BP神经网络方法的近岸数值海温预报释用技术

为了提高近岸精细化海温预报精度,利用神经网络方法,分析了海温数值预报及观测数据在释用中的作用,研究了定点近岸海温影响因子的最优配置方案,建立了定点海温精细化数值预报释用模型,评估了释用模型性能。误差分析结果显示,数值海温产品及其观测在建模中起到了降低和稳定模型误差的作用;释用模型将定点数值预报的误差从2.2°C减少至0.7°C;预报误差较调训误差略高,但考虑到预报误差的稳定性,数值释用与人工经验预报水平持平,因此,该方法具有十分广阔的拓展空间和应用前景。     

集合滤波和三维变分混合数据同化方法研究

发展了一种新的混合数据同化方法——基于集合滤波和三维变分的混合数据同化方法。该方法将集合调整卡尔曼滤波(ensemble adjustment Kalman filter,EAKF)得到的集合样本扰动通过一个转换矩阵的形式直接作用到背景场上,利用顺序滤波的思想得到分析场的一个扰动;然后在三维变分(three dimensional variational analysis,3D-Var)的框架下与观测数据进行拟合,从而给出分析场的最优估计。文中以Lorenz63模型为例,开展了理想数据同化试验,结果表明,相比于集合调整卡尔曼滤波,这种新的混合同化方法可以给出更好的同化结果。     

现代海洋/大气资料同化方法的统一性及其应用进展

海洋／大气资料同化的理论基础是用数值模式作为动力学强迫对观测信息进行提炼，或者说，从包含观测误差(噪声)的空间分布不均匀的实测资料中依据动力系统自身的演化规律(动力学方程或模式)来确定海洋／大气系统状态的最优估计。本文对主要的现代海洋／大气资料同化方法，包括最优插值(()ptimal Interpolation，简称()Ⅰ)、变分方法(3—Dimensional Variational和4—Dimensional Variational，分别简称3DVAR和4DVAR)和滤波方法(Filtering)的原理、算法设计和实际应用进行系统地回顾，并对这些资料同化方法的优缺点进行分析和讨论。在滤波框架下，所有的现代资料同化方法都被统一了：()Ⅰ和3DVAR是不随时间变化的滤波器，4DVAR和卡曼滤波是线性滤波器，即非线性滤波的退化情形；而集合滤波能构建非线性的滤波器，因为集合在某种程度上体现了系统的非高斯信息。一个非线性滤波器的主要优点是能计算和应用随时间变化的各阶误差统计距，如误差协方差矩阵。将非线性滤波器计算的随时间变化的误差协方差矩阵引入到（）Ⅰ或4DVAR中，也许能实质性地改进这些传统方法。在实际应用中，方法的优劣可能取决于所选用的数值模式和可获得的计算资源，因此需针对不同的问题选取不同的资料同化方法。由于各种资料同化方法具有统一性，因此可建立测试系统来评价这些方法，从而对各种方法获得更深入的理解，改进现有的资料同化技术，并提高人们对海洋／大气环境的预测能力。     

集合变换卡尔曼滤波数据同化方法中的协方差局地化

在集合数据同化中,协方差局地化（covariance localization,CL）方法的使用存在限制。集合转换卡尔曼滤波（ensemble transform Kalman filter,ETKF）作为集合平方根滤波的变种方法,是一种应用较广、计算高效的数据同化方法。本文分析了CL方法应用于ETKF方法的困难,从而改进CL方法使其可以适用于ETKF方法。另外,结合浅水方程,利用Askey函数作为多元局地化函数,提出了一种适用于多元数值模型的CL方法。通过具体实验验证,得到了较好的分析结果。     

GNSS/INS组合导航系统中异常观测的探测与抑制

在隧道、城市峡谷、多径效应明显等恶劣环境中,GNSS定位结果易受异常值影响。GNSS异常值的存在使得观测噪声服从Gauss分布的假设不再成立,从而将严重影响GNSS/INS组合系统中Kalman滤波的性能。本文将抗差Kalman滤波用于GNSS/INS组合导航系统中异常观测的探测与抑制。根据设定的显著性水平对Kalman滤波的新息向量进行c2检验,原假设被拒绝时认为异常观测存在;然后引入尺度化因子对新息向量的协方差矩阵进行膨胀,以抑制异常观测对滤波结果的影响,并通过解析方法求解该因子。数值仿真的结果验证了算法的有效性。     

Microbial Biomass and Organic Nutrients in the Deep-Sea Sediments of the Central Indian Ocean Basin

In order to assess the impact of deep-sea mining on the in situ benthic life, we measured the microbial standing stock and concentration of organic nutrients in the deep-sea sediments of the Central Indian Ocean Basin in the Indian pioneer area. Sediments were collected using box core and grab samples during September 1996. The total bacterial numbers ranged from 10 10 -10 11 cells per g -1 dry weight sediment. There was a marginal decrease in the number of bacteria from surface to 30 cm depth, though the subsurface section registered a higher number than did the surface. The highest numbers were encountered at depths of 4-8 cm. The retrievable number of bacteria were two orders less in comparison with the direct total counts of bacteria. An almost homogeneous distribution of bacteria, total organic carbon, living biomass, and lipids throughout the depth of cores indicates active microbial and benthic processes in the deep sea sediments. On the other hand, a uniform distribution of total counts of bacteria, carbohydrates, and total organic carbon in all the cores indicates their stable nature and suggests that they can serve as useful parameters for long-term monitoring of the area after the benthic disturbance. Further studies on temporal variability in this region would not only verify the observed norms of distribution of these variables but would also help to understand restabilization processes after the simulated benthic disturbance.     

海上大直径钢管桩打桩过程中桩周土强度弱化研究

海上大直径钢管桩打桩过程中,桩周土体受到强烈扰动而发生强度弱化,掌握桩周土体强度弱化规律对于准确预测打桩过程、保证工程安全具有重要意义。为研究土体强度弱化规律,开展了环剪试验模拟打桩对桩周土体的扰动,测试土体强度随剪切速率的变化规律,建立了描述土体强度弱化规律的拟合公式,引入到打桩分析软件中。研究结果表明：土体的强度折减程度不仅与土体本身的性质有关还受到土体的埋深和剪切速率的影响,埋深越深土体强度折减程度越低,剪切速率越高土体强度折减越高,在打桩分析中可采用这里推荐的线性折减方法来模拟不同深度处土体强度的折减规律。     

Geoacoustic Inversion via Genetic Algorithm and Its Application to Manganese Sediment Identification

An acoustic inversion method using a wide-band signal and two near field receivers is proposed and applied to multiple layered seabed models including a manganese sediment. The inversion problem can be formulated into a probabilistic model comprised of signals, a forward model, and additive noise. The forward model simulates wide-band signals, such as chirp signals, and is chosen to be the source-waveletconvolution plane wave modeling method. The wavelet matching technique, using weighted least-squares fitting, estimates the sediment sound-speed and thickness on which determination of the possible numerical ranges for a priori uniform distribution is based. The genetic algorithm is applied to a global optimization problem to find a maximum a posteriori solution for determined a priori search space. Here the object function is defined by an L 2 norm of the difference between measured and modeled signals. Not only the marginal pdf but also its statistics are calculated by numerical evaluation of integrals using the samples selected during importance sampling process of the genetic algorithm.     

A Review of Some Environmental Issues Affecting Marine Mining

This article reviews information recently available from existing marine and coastal mining for responses to environmental issues affecting marine mining at different depths. It is particularly but not exclusively concerned with those issues affecting seabed biodiversity impact and recovery. Much information has been gathered in the past 10 years from shallow mining operations for construction aggregate, diamonds, and gold, from coastal mines discharging tailings to shallow and deep water, and from experimental deep mining tests. The responses to issues identified are summarized in a series of eight tables intended to facilitate site-specific consideration. Since impacts can spread widely in the surface mixing layer SML, and can affect the biologically productive euphotic zone, the main issues considered arise from the depth of mining relative to the SML of the sea. Where mining is below the SML, the issue is whether it is environmentally better to bring the extraction products to the surface vessel for processing (and waste discharge), or to process the extraction products as much as possible on the seabed. Responses to the issues need to be sitespecific, and dependent on adequate preoperational environmental impact and recovery prediction. For deep tailings disposal from a surface vessel, there are four important environmental unknowns: (1) the possible growth of "marine snow" (bacterial flocs) utilizing the enormous quantities of fine tailings particles (hundreds or thousands of metric tons per day) as nuclei for growth, (2) the possibility that local keystone plankton and nekton species may migrate diurnally down to and beyond the depth of deep discharge and hence be subjected to tailings impact at depth, (3) the burrow-up capability of deep benthos and their ability to survive high rates of tailings deposition, and (4) the pattern and rate of dispersion of a tailings density current through the deep water column from discharge point to seabed. Actions to obtain relevant information in general and site-specifically are suggested.     

Increased Particle Fluxes at the INDEX Site Attributable to Simulated Benthic Disturbance

Particle fluxes were measured 7 m above the sea bottom during the predisturbance, disturbance, and postdisturbance periods by using time series sediment traps attached to seven deep-sea moorings deployed in the INDEX experiment site in the Central Indian Basin. The predisturbance particle fluxes varied between 22.3 to 55.1 mg m -2 day -1 . Increased and variable particle fluxes were recorded by the sediment traps during the disturbance period. The increase observed was 0.5 to 4 times more than the background predisturbance fluxes. The increases in particle fluxes (~4 times) recorded by the sediment trap located in the southwestern direction (DMS-1) were the greatest, which could be the result of preferential movement of resuspended particles generated during the deep-sea benthic disturbance along the general current direction prevailing in this area during the experimental period. Also, the traps located closer to the disturbance area recorded greater fluxes than did the traps far away, across the Deep Sea Sediment Resuspension System path. This variability in recorded particle fluxes by the traps around the disturbance area clearly indicates that physical characteristics such as grain size and density of the resuspended particles produced during the disturbance had an important effect on particle movement. The postdisturbance measurements during ~5 days showed a reduction in particle fluxes of ~50%, indicating rapid particle settlement.     

Phylogenetic diversity and phenotypic characterization of cultivable bacterioplankton isolated from polar oceans

A set of 27 marine planktonic bacteria isolated from the polar regions was characterized by 16S rDNA sequencing and physiological and biochemical testing. More than half of these bacteria were positive for caseinase, gelatinase and 13-glucosidase, and could utilize glucose, maltose or malic acid as carbon source for cell growth. Twelve isolates expressed nitrate reduction activities. Except for one antarctic isolate BSwlO175 belonging to Actinobacteria phylum, these isolates were classified as γ-Proteobacteria, suggesting that γ-Proteobacteria dominated in cultivable marine bacterioplankton at both poles. Genus Pseudoalteromonas was the predominant group in the Chukchi Sea and the Bering Sea, and genus ShewaneUa dominated in cultivable bacterioplankton in the Prydz Bay. With sequence similarities above 97%, genus Psychrobacter was found at both poles. These 27 isolates were psychrotolerant, and significant 16S rDNA sequence similarities were found not only between arctic and antarctic marine bacteria （ 〉 99% ）, but also between polar marine bacteria and bacteria from other aquatic environments （ ≥ 98.8% ）, including temperate ocean, deep sea, pond and lake, suggesting that in the polar oceans less temperature-sensitive bacteria may be cosmopolitan and have a bipolar, even global, distribution at the species level.     

Determination of In Situ Sediment Shear Strength from Advanced Piston Corer Pullout Forces

The advanced piston corer (APC) has been used by the Ocean Drilling Program since 1985 for recovering soft sediments from the ocean floor. The pullout force measured on extracting the core barrel from the sediment is shown to correlate with the average shear strength of the sediment core measured in the ship's laboratory. A simple rule of thumb is derived relating the shear strength of the sediment to the pullout force. Multiple APC holes at individual sites allow the consistency of the pullout measurements to be assessed. The effects of different operational procedures during APC coring are also explored. Although generally applicable, the correlation between pullout force and laboratory measurements of shear strength breaks down for some APC holes, possibly because of the disturbance of some sediment types during the APC coring process. A better understanding of the physical process of APC coring, and its effect on the properties of the sediment both inside and immediately outside the core barrel, would indicate what confidence can be put on the measurement of pullout force as a way of evaluating the in situ shear strength of deep sea sediments.     

Deep-Sea Sediment Compression Curves: Some Controlling Factors, Spurious Overconsolidation, Predictions, and Geophysical Reproduction

Oedometer tests have been carried out on 70 undisturbed surficial clays (at approximately 250 mm below the mudline), mostly collected by free-fall corers from sites widely scattered throughout the deep-sea North Atlantic. Acoustic measurements were also made, initially on contiguous samples and ultimately on the same sample using a geophysically instrumented oedometer which also collected electrical resistivity data. Apart from those quiescent areas below the carbonate compensation depth, such as north of the West Indies where very fine clays exist, most of the samples are silty clays whose geotechnical-geophysical properties are dependent on the type of clay minerals present (and their ability to take in moisture), the sand-size fraction, and the quantity of carbonate present. Thus the pure clays have high compressibilities which decrease on the addition of coarse particles, while the converse is true for the acoustic parameters, these increasing with the sand fraction. Using the notion of the intrinsic compression line for all samples, and comparison to it of the measured compression curves, it is clear that, contrary to some previously held ideas, most deep-sea clays are normally consolidated; the addition of carbonate has the effect of creating an open, stronger sediment skeleton. Interestingly, where information is available, the variation with depth of a sample's acoustic velocity follows the void ratio pressure relationship of the compression curve. This allows the construction of an in-situ sediment compression curve using the in-situ geophysical observations.