Least Square Finite Element Method for Viscous Splitting of Unsteady Incompressible Navier–Stokes Equations

In order to solve unsteady incompressible Navier–Stokes(N–S) equations, a new stabilized finite element method,called the viscous-splitting least square FEM, is proposed. In the model, the N–S equations are split into diffusive and convective parts in each time step. The diffusive part is discretized by the backward difference method in time and discretized by the standard Galerkin method in space. The convective part is a first-order nonlinear equation.After the linearization of the nonlinear part by Newton's method, the convective part is also discretized by the backward difference method in time and discretized by least square scheme in space. C~0-type element can be used for interpolation of the velocity and pressure in the present model. Driven cavity flow and flow past a circular cylinder are conducted to validate the present model. Numerical results agree with previous numerical results, and the model has high accuracy and can be used to simulate problems with complex geometry.     

Penetration Depth of Torpedo Anchor in Two-Layered Cohesive Soil Bed by Free Fall

The penetration depth of torpedo anchor in two-layered soil bed was experimentally investigated. A total of 177 experimental data were obtained in laboratory by varying the undrained shear strength of the two-layered soil and the thickness of the top soil layer. The geometric parameters of the anchor and the soil properties (the liquid limit, plastic limit, specific gravity, undrained shear strength, density, and water content) were measured. Based on the energy analysis and present test data, an empirical formula to predict the penetration depth of torpedo anchor in two-layered soil bed was proposed. The proposed formula was extensively validated by laboratory and field data of previous researchers. The results were in good agreement with those obtained for two-layered and single-layered soil bed. Finally, a sensitivity analysis on the parameters in the formula was performed.     

粉质土海底边界层动态演化过程波浪水槽实验

The objectives of this study are carried out a series of controlled large wave flume experiments using fine-grained sediment from the Huanghe River Delta, exploring the complete sequence of sediment behavior in the bottom boundary layer(BBL) during wave-induced liquefaction. The results show that:(1) The BBL in silty seabed is exposed to a progressive wave, goes through a number of different stages including compaction before liquefaction, sediment liquefaction, and compaction after liquefaction, which determines the range and thickness of BBL.(2) With the introduction of waves, first, the sediment surface has settled by an amount S(S=1–2 cm) in the course of wave loadings with an insufficient accumulation of pore water pressure. And a thin high concentration layer formed the near-bed bottom.(3) Once the liquefaction sets in, the liquefied sediment with an ‘orbital motion' and the sub-liquefied sediment form a two-layer-sediment region. The range of BBL extends downwards and stopped at a certain depth, subsequently, develops upwards with the compaction process. Meanwhile, resuspended sediments diffuse to the upper water column.(4) During the dynamics process of the BBL beneath progressive waves, the re-suspended sediment increment ranked as sediment liquefaction erosion before liquefaction compaction after liquefaction.     

微型近红外光谱仪测量北极积雪覆盖海冰内部太阳辐照度可行性研究

The extremely low temperature, high humidity and limited power supply pose considerable challenges when using spectrometers within the Arctic sea ice. The feasibility of using a miniature low-power near-infrared spectrometer module to measure solar radiation in Arctic sea ice environments was investigated in this study.Temperature and integration time dependences of the spectrometer module were examined over the entire target operating range of –50℃ to 30℃, well below the specified operating range of this spectrometer. Using these observations, a dark output prediction model was developed to represent dark output as a function of temperature and integration time. Temperature-induced biases in the saturation output and linear operating range of the spectrometer were also determined. Temperature and integration time dependences of the signal output were evaluated. Two signal output correction models were developed and compared, to convert the signal output at any temperature within the operating temperature range and integration time to that measured at the reference temperature and integration time. The overall performance of the spectrometer was evaluated by integrating it into a refined fiber optic spectrometry system and measuring solar irradiance distribution in the ice cover with thickness of 1.85 m in the Arctic during the 9th Chinese National Arctic Research Expedition. The general shape of the measured solar irradiance above the snow surface agreed well with that measured by other commercial oceanographic spectroradiometers. The measured optical properties of the sea ice were generally comparable to those of similar ice measured using other instruments. This approach provides a general framework for assessing the feasibility of using spectrometers for applications in cold environments.     

安达曼海东部凹陷渐新世以来断裂–构造演化特征及其成因探讨

本文研究内容为印度洋东北部边缘海安达曼海的构造演化。利用安达曼海域东部大范围二维地震数据资料及钻井数据,结合区域地质概况以及前人研究成果,选取8条具有代表性的断层并将其划分为一级和二级断裂,运用生长指数法和古落差法对断层进行定量分析,再通过计算安达曼海东部凹陷4条主测线的构造沉降量,探讨构造演化过程。结果表明：选取的生长断层中3条属于一级断裂,跨度大,几乎切穿整个地层,属于控制安达曼海域地区沉降的大断裂;另外5条属于二级断裂,控制构造带的展布情况,属于构造带的分界线。渐新世时期,印度–澳大利亚板块与欧亚板块之间处于软碰撞阶段,断层发育缓慢,上下盘落差较小,生长指数与构造沉降量也处于低值;中新世时期,板块之间的耦合效应不断增强,断层发育速度加快,此时上下盘厚度最大,是形成多处断裂带以及多种断裂样式的关键时期,各地层生长指数和构造沉降量也达到峰值。上新世至今,安达曼海沟–弧–盆体系逐渐稳定,断层活动减弱,断裂上下盘厚度差基本一致,生长指数差异较小,构造沉降量基本稳定在 1 km 左右。     

Settlement Mode Analysis for An Immersed Tube Tunnel Considering A Nonuniform Foundation Under Tidal Load

Immersed tube tunnels are usually placed on soft soil layers in cross-sea tunnelling engineering. Owing to the influence of stratum conditions and slope design, the longitudinal distribution of substratum layers is generally uneven. Thus, the inhomogeneous deformation of the element-joint becomes the key factor in the failure of the immersed tube tunnel. Therefore, a corresponding calculation method for joint deformation is needed to explore the deformation law of immersed tube tunnels. By constructing a three-section immersed tube tunnel analysis model (TTM), the relationship between the two types of deformation of the immersed tube tunnel structure in a longitudinal nonuniform soft soil foundation is described, and the deformation characteristics of the immersed structure under different boundaries are discussed. Based on the mechanical behaviour of the joint and foundation, according to the Timoshenko beam on the Vlasov two-parameter foundation (VTM), considering the tidal cyclic load during the operation and maintenance period, an example analysis is given. Moreover, the deformation characteristics and development trend of the immersed tube tunnel under the influence of different soil layers are discussed. The obtained results have a certain guiding significance for the deformation calculation of immersed tube tunnels.

     

Semisupervised heterogeneous ensemble for ship target discrimination in synthetic aperture radar images

Ship detection using synthetic aperture radar (SAR) plays an important role in marine applications. The existing methods are capable of quickly obtaining many candidate targets, but numerous non-ship objects may be wrongly detected in complex backgrounds. These non-ship false alarms can be excluded by training discriminators, and the desired accuracy is obtained with enough verified samples. However, the reliable verification of targets in large-scene SAR images still inevitably requires manual interpretation, which is difficult and time consuming. To address this issue, a semisupervised heterogeneous ensemble ship target discrimination method based on a tri-training scheme is proposed to take advantage of the plentiful candidate targets. Specifically, various features commonly used in SAR image target discrimination are extracted, and several acknowledged classification models and their classic variants are investigated. Multiple discriminators are constructed by dividing these features into different groups and pairing them with each model. Then, the performance of all the discriminators is tested, and better discriminators are selected for implementing the semisupervised training process. These strategies enhance the diversity and reliability of the discriminators, and their heterogeneous ensemble makes more correct judgments on candidate targets, which facilitates further positive training. Experimental results demonstrate that the proposed method outperforms traditional tri-training.     

Automated multi-scale classification of the terrain units of the Jiaxie Guyots and their mineral resource characteristics

Given the advances in satellite altimetry and multibeam bathymetry, benthic terrain classification based on digital bathymetric models (DBMs) has been widely used for the mapping of benthic topographies. For instance, cobalt-rich crusts (CRCs) are important mineral resources found on seamounts and guyots in the western Pacific Ocean. Thick, plate-like CRCs are known to form on the summit and slopes of seamounts at the 1 000–3 000 m depth, while the relationship between seamount topography and spatial distribution of CRCs remains unclear. The benthic terrain classification of seamounts can solve this problem, thereby, facilitating the rapid exploration of seamount CRCs. Our study used an EM122 multibeam echosounder to retrieve high-resolution bathymetry data in the CRCs contract license area of China, i.e., the Jiaxie Guyots in 2015 and 2016. Based on the DBM construted by bathymetirc data, broad- and fine-scale bathymetric position indices were utilized for quantitative classification of the terrain units of the Jiaxie Guyots on multiple scales. The classification revealed four first-order terrain units (e.g., flat, crest, slope, and depression) and eleven second-order terrain units (e.g., local crests, depressions on crests, gentle slopes, crests on slopes, and local depressions, etc.). Furthermore, the classification of the terrain and geological analysis indicated that the Weijia Guyot has a large flat summit, with local crests at the southern summit, whereas most of the guyot flanks were covered by gentle slopes. “Radial” mountain ridges have developed on the eastern side, while large-scale gravitational landslides have developed on the western and southern flanks. Additionally, landslide masses can be observed at the bottom of these slopes. The coverage of local crests on the seamount is ~1 000 km2, and the local crests on the peak and flanks of the guyots may be the areas where thick and continuous plate-like CRCs are likely to occur.     

Reclamation-oriented spatiotemporal evolution of coastal wetland along Bohai Rim,China

Coastal wetlands are located in the ecotone of interaction between the land surface and sea, and anthropogenic activities extensively interfere with these wetlands through the reclamation of large tidal wetlands and destruction of the function of the ecosystems. In this study, we investigated the dynamic evolutionary characteristics of the Bohai Rim coastal area over the past 40 years using the Modified Normalized Difference Water Index, the fractal dimension, object-oriented classification, the land-use transfer trajectory, and regression analysis. Additionally, we quantified and monitored the evolution of reclamation and analyzed the correlation between reclamation and coastal wetlands based on 99 Landsat-2, -5, and -8 images (at 60 m and 30 m spatial resolution) over the period 1980–2019. The results are as follows. (1) The coastline of the Bohai Rim increased by 1 631.2 km from 1980 to 2019 with a zigzag variation. The artificial coastline increased by 2 946.1 km, whereas the natural coastline decreased by 90%. (2) The area of man-made wetlands increased by 3 736.9 km2, the area of construction land increased by 1 008.4 km2, and the natural wetland area decreased by 66%. The decrease of tidal flats is the main contributor to the decrease of natural wetland area (takes account for 91.1%). Coastal areas are affected by intense human disturbance, which was taken place across a large area of tidal flats and caused the landscape to fragment and be more heterogeneous. The coastal zone development activities were primarily concentrated in the southern Laizhou Bay, the Yellow River Delta, the Bohai Bay, the northern Liaodong Bay, and the Pulandian Bay. The solidified shorelines and increase in sea level have resulted in intertidal wetlands decreasing and impaired wetland ecology. (3) There is a good agreement between reclamation and the size of the coastal wetlands. Both land reclamation and the reduction in coastal wetland areas are significantly related to the population size, fishery output value, and urbanization rate. In summary, human activities, such as the construction of aquaculture ponds and salt pans, industrialization, and urbanization, are the primary forces that influence the environmental changes in the coastal region. This study is beneficial for establishing and improving the systems for the rational development and utilization of natural resources, and provides theoretical references for restoring wetland ecology and managing future reclamation activities in other coastal zone-related areas.     

A fast inversion method for ocean parameters based on dispersion curves with a single hydrophone

The dispersion characteristics of shallow water can be described by the dispersion curves, which contain substantial ocean parameter information. A fast ocean parameter inversion method based on dispersion curves with a single hydrophone is presented in this paper. The method is achieved through Bayesian theory. Several sets of dispersion curves extracted from measured data are used as the input function. The inversion is performed by matching a replica calculated with a dispersion formula. The bottom characteristics can be described by the bottom reflection phase shift parameter P. The propagation range and the depth can be inverted quickly when the seabed parameters are represented by on parameter P. The inversion results improve the inversion efficiency of the seabed parameters. Consequently, the inversion efficiency and accuracy are improved while the number of inversion parameters is decreased and the computational speed of replica is increased. The inversion results have lower error than the reference values, and the dispersion curves calculated with inversion parameters are also in good agreement with extracted curves from measured data; thus, the effectiveness of the inversion method is demonstrated.