Methods and procedures to determine the outer limits of the continental shelf beyond 200 nautical miles

This paper establishes techniques and methods to determine a variety of boundaries associated with 200 nautical miles beyond the continental shelf. The methods, based on topography, slope and second-deriv- ative profile integrated analysis, are now able to identify automatically the foot of the continental slope （FOS）. By analyzing the sedimentary profile, the points of 1% sediment thickness are recognized. Through the intersection, cut, deletion and mergence calculation of the extrapolated data set of fLxed-point series, the method succeeds in generating automatically the extrapolated boundaries, including the FOS＋60 M line, the 350 M line, and the 2500 m＋100 M line. In addition, based on the automatic analysis of the topographic profile, it can be applied to determine rapidly the points of maximum water depth. Taking the northern Okinawa Trough （OT） as an example, these methods are used to calculate and examine the boundaries in- cluded in the Submission by the People＇s Republic of China Concerning the Outer Limits of the Continental Shelf beyond 200 Nautical miles in Part of the East China Sea （ECS）; the boundaries thus derived have a solid scientific and rational basis.     

Underwater Terrain Positioning Method Based on Least Squares Estimation for AUV

To achieve accurate positioning of autonomous underwater vehicles, an appropriate underwater terrain database storage format for underwater terrain-matching positioning is established using multi-beam data as underwater terrainmatching data. An underwater terrain interpolation error compensation method based on fractional Brownian motion is proposed for defects of normal terrain interpolation, and an underwater terrain-matching positioning method based on least squares estimation(LSE) is proposed for correlation analysis of topographic features. The Fisher method is introduced as a secondary criterion for pseudo localization appearing in a topographic features flat area, effectively reducing the impact of pseudo positioning points on matching accuracy and improving the positioning accuracy of terrain flat areas. Simulation experiments based on electronic chart and multi-beam sea trial data show that drift errors of an inertial navigation system can be corrected effectively using the proposed method. The positioning accuracy and practicality are high, satisfying the requirement of underwater accurate positioning.     

Underwater Terrain Positioning Method Based on Least Squares Estimation for AUV

To achieve accurate positioning of autonomous underwater vehicles, an appropriate underwater terrain database storage format for underwater terrain-matching positioning is established using multi-beam data as underwater terrain-matching data. An underwater terrain interpolation error compensation method based on fractional Brownian motion is proposed for defects of normal terrain interpolation, and an underwater terrain-matching positioning method based on least squares estimation (LSE) is proposed for correlation analysis of topographic features. The Fisher method is introduced as a secondary criterion for pseudo localization appearing in a topographic features flat area, effectively reducing the impact of pseudo positioning points on matching accuracy and improving the positioning accuracy of terrain flat areas. Simulation experiments based on electronic chart and multi-beam sea trial data show that drift errors of an inertial navigation system can be corrected effectively using the proposed method. The positioning accuracy and practicality are high, satisfying the requirement of underwater accurate positioning.     

珠江三角洲思贤滘水文泥沙特征和现代演变

根据思贤滘的水文泥沙、地形和底质资料,分析了思贤滘的河床底质类型、分布规律、沉积物粒度组成、重矿物含量以及河床演变特征等,讨论了思贤滘的泥沙运移特点,指出西江和北江的推移质均已进入思贤滘,但未通过,因此,对思贤滘凹岸堤围的险段岁修工程不容忽视;并对思贤滘的整治规划提出了不设闸的刍议。     

面向海底成图基于DTM边界的等值线充填算法

提出一种基于DTM边界线的等值线矢量彩色充填图算法.算法在处理奇异格网点基础上跟踪并存储DTM边界线,通过边界点阵快速建立非封闭等值线与边界线间的双向联系,基于边界线建立非封闭等值线间的拓扑关系并顺序跟踪封闭多边形,通过多边形最小扩展框和转角法(或射线法)建立多边形间的嵌套关系.该算法已成功应用于海量多波束数据成图系统MBChart,解决了常用多波束后处理成图软件栅格充填与等值线之间的失配问题     

Reconstruction of the three-dimensional salinity and temperature fields of the Black Sea on the basis of satellite altimetry measurements

The methodology of joint processing of the satellite altimetry and occasional hydrological observations in the Black Sea for 1993–2012 is developed. The original technique for reconstruction of the 3D temperature and salinity fields in the deep-sea part is proposed and implemented. This technique makes it possible to identify the depths at which a contribution of adiabatic processes to the deformation of the temperature and salinity profiles of the sea is predominant. Daily-averaged 3D fields of the seawater temperature and salinity in a baroclinic layer on a regular grid are reconstructed. The evaluation of accuracy of the reconstructed temperature and salinity arrays is performed by comparing them with the data of hydrological exploration. Structures of the temperature and salinity fields are correlated naturally with topography of the altimetric level and clearly indicate the synoptic variability. Seasonal and interannual variabilities of the kinetic energy (averaged over horizons of the 63–400 m layer) of the geostrophic currents calculated using the dynamic method makes it possible to reveal a sharp increase in the kinetic energy of the currents in the winter season of 2002. A high correlation is found between the interannual variability of the ERA-Interim wind stress curl averaged over the surface of the deep sea part and the kinetic energy of the geostrophic currents in the 63–400 m layer.     

Modeling bistatic bottom scattering strength including a forwardscatter lobe

An upgrade to bistatic scattering strength modelling that is based on the authors' current understanding of bottom topographic scattering with an emphasis on modeling the `forward lobe' where Lambert's law fails quite significantly is reported. Low-frequency bottom scatter modeling is reviewed with particular emphasis on the issues of the forward scattered lobe. A specific model (a modified version of BISSM) is proposed, and the model's advantages and limitations are discussed. The requirement for certain high-resolution geomorphic data needed to support the model is discussed. Like the original BISSM, the version does not modify the accepted form for diffuse scattering, but it does modify the form of the forward lobe     

丁坝间潮滩地貌变化的经验正交函数分析

经验正交函 数(EOF)是用来处理原始数据、进行多元统计分析的一种方法。本文将原始距平矩阵分解为空间特征函数和相应的时间特征函数,并用前3个经验正交函数,即表示潮滩总体冲淤变化的第一特征函数,表示季节性冲淤变化的第二特征函数和表示偶然因素扰动引起冲淤变化的第三特征函数的线性组合来反映潮滩的变化,分析浙江瑞安淤泥质潮滩在建立丁坝后的演变规律。初步探讨了动力过程和地貌之间的关系,分析了淤泥质潮滩在丁坝工程影响下的时空变化特征与动力过程的关系;并对第二特征函数对应的动力过程进行详细阐述,结果显示波浪作用是影响浙江瑞安淤泥质潮滩季节性变化的主要因素。     

基于结构特征的遥感影像海岸码头分割方法

在分析水边线对码头提取的作用和海岸码头的结构特征的基础上,提出了一种基于结构特征的遥感影像海岸码头分割方法。该方法在遥感影像水边线提取的基础上,利用道格拉斯-普克算法通过压缩得到关键节点,接着在节点中选取满足要求的凹点,然后依据海岸结构特征通过对凹点的筛选分割出海岸规则突出物,最后通过对规则突出物的分析验证得到海岸码头。实验结果表明,该方法能够从水边线中将码头区域分割出来,并且受停靠船舶等噪声影响较小,为下一步海岸码头的精确提取奠定了坚实基础。     

围海工程前后瓯飞滩冲淤变化特征及动力机制分析

围海工程拓展出陆地空间的同时对工程海域地形变化产生较大的影响。本文以温州瑞安东部瓯飞滩海域为研究区,通过GIS处理水下地形资料,构建Delft3D数学模型,分析围海工程引起的瓯飞滩地形冲淤变化特征及其动力机制。研究表明,工程实施后研究区海域涨潮占优特性有所增强,近岸海域整体保持淤积态势,淤积速率略有增加;局部区域受影响较大,如围堤凸区附近受挑流作用表现为缓慢冲刷,而在围堤凹处受阻流效应影响表现为快速淤涨。研究对瓯飞滩及类似淤泥质海岸的开发利用与保护有重要的参考价值。     

Depth Contour Smoothing Based on the Fitting of Multi-Segment Bezier Curves

To overcome smoothness inconsistencies in depth contours interpolated from sounding data, this paper proposes a smoothing method based on the fitting of multi-segment Bezier curves considering navigational safety. The basic principle of this method is the division of a generated depth contour into a series of bends categorized as either convex bends or concave bends. Different smoothing strategies are designed for different bends while considering the safety constraint. Convexity-preserving Bezier curves are used to fit convex bends. For concave bends, the line segment of the bend is utilized as a boundary to construct a control polygon within the interior of the bend, after which a Bezier curve is constructed based on the polygon. Experiments with real data are implemented, and the results are compared with those of other methods. The results show that the smoothing procedure for depth contours with angularities is effective and respects the safety constraint; moreover, the shape characteristics are preserved, and high positional accuracy is achieved. In addition, the proposed smoothing method is free of topological inconsistencies.     

海底沙波特征线的最优方向剖面自动识别方法

海底沙波是发育在近海陆架上的一种常见海底地貌类型,海底沙波特征与运动规律的研究具有重要的科学意义与工程应用价值,沙波脊线与谷线是表征海底沙波的最基本特征,也是精确描述沙波运动的基本参量。本文提出了一种基于复合数字水深模型的沙波特征线自动识别方法——最优方向剖面法,基于水深曲面归算得到最优剖面方向,再依据最优剖面方向求导并判定极值,自动提取沙波形态特征点,最终形成沙波脊线和谷线。以台湾浅滩复合型沙波为例进行对比实验研究,结果表明,该方法能基于不同分辨率的数字水深模型自动准确地提取海底沙波脊线与谷线,勿需设置阈值,地形自动化识别程度得到进一步提升,具有重要的实际应用价值。     

Plane wave fitting of the internal tidal sea-surface height signatures extracted from the TOPEX/Poseidon altimeter data along satellite tracks

The sea-surface height (SSH) signatures of internal tides extracted from the TOPEX/Poseidon (T/P) altimeter data along satellite tracks are fitted with superposition of several plane waves which have different wavenumber vectors. The key problem of plane wave fitting with iterative method is how to determine the initial value of wavenumber of each plane wave. The previous solving method is to analyze the internal tidal SSH signatures along each track with wavenumber spectrum. But it is found that the problem cannot be solved completely with the wavenumber spectrum analysis method only. The method based on the combination of wavenumber spectrum analysis method and the exhaustive method is proposed to determine the initial values of wavenumbers for iteration. Numerical results indicate that the proposed method is not only reasonable and feasible but also better than the previous method. The proposed method is an improvement of the previous one, which is beneficial to improving the precision of plane wave fitting of the T/P internal tidal SSH signatures and deepening the understanding of the internal tides in ocean.     

Possible spreading centers in the Japan sea

The magnetic anomaly profiles across the Japan and the Yamato Basins are analyzed to find possible spreading centers. Magnetic anomalies are first transformed and then the correlation coefficient of symmetry of the pattern of profile is computed along each profile. If the coefficient of symmetry has a notable maximum at a point and such a point lines up for different profiles, the line of these points of maximum symmetry may be a spreading center. Three possible spreading centers are obtained. Two are situated at about the axis of the Japan and the Yamato Basins, and the other is at the topographic extension of the Yamato Ridge. The age of the opening has not been identified on the magnetic time scale yet, but the strike of the proposed spreading centers suggests a pre-Eocene time.     

Real-Time Adjustment Underwater Positioning System for Submarines

A rapid real-time adjustment scheme is proposed for improving the precision of the conventional short-base line (SBL) positioning fix system used by submarines and other underwater vehicles. In the proposed approach, an initial position estimate is obtained by solving the conventional SBL tracking equations of the submarine given the assumptions of a constant speed of sound in water and a straight-line propagation path. In the first stage of the real-time adjustment procedure, this initial estimate is corrected using an iterative computation scheme based on a 3D geometry model. The improved position estimate is then used to compute a new, more accurate value of the speed of sound in water. Finally, in the second stage of the real-time adjustment procedure, the corrected speed of sound in water and the discrepancy between the original and corrected position estimates obtained in the first adjustment procedure are applied to update the coordinates of the submarine based on the second signal received from the pinger. The numerical results show that the proposed real-time adjustment system yields a significant improvement in the accuracy of the positioning fix estimates compared to those obtained from the conventional SBL method or the SBL method with the first adjustment procedure only.     

基于Hough变换的DInSAR环状地表形变区自动提取

针对当前DInSAR地表形变区检测依赖于人工的现状,提出基于Hough变换的环状地表形变区自动提取方法。该方法在外部DEM的辅助下,由干涉数据获取去地形效应的差分干涉图,在对差分干涉图进行边缘提取的基础上,利用Hough变换实现环状地表形变区的自动提取。采用青海门源震区Sentinel-1A数据和山东济宁矿区ALOS PALSAR数据分别进行环状地表形变区自动提取实验,验证了本文方法的有效性。     

海上风机对潮间带地形变化影响的定量遥感监测研究

为科学分析潮间带风电场建设对潮间带滩涂稳定性的影响,文章综述国内外潮间带动态变化分析的研究进展,提出假想沙体表面判别法（HSSD）;从单时相DEM中定量估算海上风机建设引起的潮间带地形变化,分析海上风机建设及运行对近海沙体地形冲淤变化的贡献率和贡献方式;采用增强型水边线法（EWM）构建多时相潮间带DEM,开展潮间带滩涂冲淤变化影响海上风机安全性和稳定性的定量评估。实践应用证明,该方案可用于海上风机对潮间带地形变化影响的定量遥感监测。     

Multi-point tidal prediction using artificial neural network with tide-generating forces

This paper presents a neural network model of simulating tides at multi-points considering tide-generating forces. A comparison on the root mean square and correlation coefficient of three-year mixed tides at a single point computed with harmonic method, response–orthotide method, the NAO.99b model and the proposed model was made to show the prediction accuracy of each method. The proposed model is examined efficient as the harmonic method to estimate the tides at a single point. Extended application of the proposed model to predicting tides at some points neighboring to an original interest point identifies accurately simulating multi-point tides as the NAO.99b numerical model.     

Statistically average atmospheric bending angle model based on COSMIC experimental data

The retrieval of profiles of meteorological variables from radio occultation observations requires knowledge of bending angle profiles up to heights of no less than 60–70 km. Because of the residual error of the ionospheric correction, retrieved profiles become too noisy by a height of about 40 km. In order to invert the bending angle profiles, the statistical optimization is used. This makes it possible to construct an optimal linear combination of the a priori estimate of the average bending angle profile and a posteriori noisy estimate based on observations. The estimate of the average bending angle profile for the given coordinates and the time of year is usually based on the climatological atmospheric model. MSIS and CIRA models that have been used for this purpose are now obsolete and do not describe the global changes in the atmospheric state. The model of average bending angles BA–IAP (Bending Angle–Institute of Atmospheric Physics) is built based on the processing of the array of COSMIC radio occultation observations during 2006–2013. The proposed model is statistically validated based on the COSMIC database. It is shown that our model describes the average bending angle profiles more accurately than the MSIS model.     

Static and dynamic analyses of three-dimensional hollow concrete block revetments using polyhedral finite element method

In this study, we exploit the flexibility and advantages of polyhedral finite elements in modeling 3D hollow concrete block revetments. In the present method, a new concept is proposed, in which the polyhedral elements are subdivided into sub-tetrahedrons through virtual nodes located at the centroids of the faces and element. Piecewise linear shape functions of polyhedral elements are constructed for the sub-tetrahedrons. Then, the shape functions are used to formulate the stiffness and mass matrices. By this way, the present approach yields lower computational cost than that of the standard polyhedral finite element method based on rational basis functions. The presented method suits well general polyhedral meshes including concave elements, which are used to model 3D complicated structures in design of coastal defenses system protection such as interlocking revetments. Reliability and effectiveness of the present approach are validated by the well-known commercial ANSYS software through static and dynamic analysis.