Paleoceanographic records in the sedimentary cores from the middle Okinawa Trough

Two gravity piston cores (Cores 155 and 18) involved in this study were collected from the middle Okinawa Trough. Stratigraphy of the two cores was divided and classified based on the features of planktonic foraminifera oxygen isotope changes together with depositional sequence, millennium-scale climatic event comparison, carbonate cycles and AMS(14)C dating. Some paleoclimatic information contained in sediments of these cores was extracted to discuss the paleoclimatic change rules and the short-time scale events presented in interglacial period. Analysis on the variation of oxygen isotope values in stage two shows that the middle part of the Okinawa Trough may have been affected by fresh water from the Yellow River and the Yangtze River during the Last Glacial Maximum (LGM). The oxygen isotope value oscillating ranges of the cores have verified that the marginal sea has an amplifying effect on climate changes. The delta(13)C of benthic foraminifera Uvigerina was lighter in the glacial period than that in the interglacial period, which indicates that the Paleo-Kuroshio's main stream moved eastward and its influence area decreased. According to the temperature difference during the "YD" period existing in Core 180 and other data, we can reach the conclusion that the climatic changes in the middle Okinawa Trough area were controlled by global climatic changes, but some regional factors had also considerable influence on the climate changes. Some results in this paper support Fairbanks's point that the "YD" event was a brief stagnation of sea level rising during the global warming up procession. Moreover, the falling of sea level in the glacial period weakened the exchange between the bottom water of the Okinawa Trough and the deep water of the northwestern Pacific Ocean and resulted in low oxygen state of bottom water in this area. These procedures are the reasons for carbonate cycle in the Okinawa Trough area being consistent with the "Atlantic type" carbonate cycle.

Impact of difference accuracy on computational properties of vertical grids for a nonhydrostatic model

Starting from nonhydrostatic Boussinesq approximation equations, a general method is introduced to deduce the dispersion relationships. A comparative investigation is performed on inertia-gravity wave with horizontal lengths of 100, 10 and 1 km. These are examined using the second-order central difference scheme and the fourth-order compact difference scheme on vertical grids that are currently available from the perspectives of frequency, horizontal and vertical component of group velocity. These findings are compared to analytical solutions. The obtained results suggest that whether for the second-order central difference scheme or for the fourth-order compact difference scheme, Charny-Phillips and Lorenz ( L) grids are suitable for studying waves at the above-mentioned horizontal scales; the Lorenz time-staggered and Charny-Phillips time staggered (CPTS) grids are applicable only to the horizontal scales of less than 10 km, and N grid ( unstaggered grid) is unsuitable for simulating waves at any horizontal scale. Furthermore, by using fourth-order compact difference scheme with higher difference precision, the errors of frequency and group velocity in horizontal and vertical directions produced on all vertical grids in describing the waves with horizontal lengths of 1, 10 and 100 km cannot inevitably be decreased. So in developing a numerical model, the higher-order finite difference scheme, like fourth-order compact difference scheme, should be avoided as much as possible, typically on L and CPTS grids, since it will not only take many efforts to design program but also make the calculated group velocity in horizontal and vertical directions even worse in accuracy.

Characteristics of shallow gas hydrate in Okhotsk Sea

Multidisciplinary field investigations were carried out in Okhotsk Sea by R/V Akademik M.A. Lavrentyev (LV) of the Russian Academy of Sciences (RAS) in May 2006, supported by funding agencies from Korea, Russia, Japan and China. Geophysical data including echo-sounder, bottom profile, side-scan-sonar, and gravity core sample were obtained aimed to understand the characteristics and formation mechanism of shallow gas hydrates. Based on the geophysical data, we found that the methane flare detected by echo-sounder was the evidence of free gas in the sediment, while the dome structure detected by side-scan sonar and bottom profile was the root of gas venting. Gas hydrate retrieved from core on top of the dome structure which was interbedded as thin lamination or lenses with thickness varying from a few millimeters to 3 cm. Gas hydrate content in hydrate-bearing intervals visually amounted to 5%-30% of the sediment volume. This paper argued that gases in the sediment core were not all from gas hydrate decomposition during the gravity core lifting process, free gases must existed in the gas hydrate stability zone, and tectonic structure like dome structure in this paper was free gas central, gas hydrate formed only when gases over-saturated in this gas central, away from these structures, gas hydrate could not form due to low gas concentration.

Tectonics in the northwestern West Philippine basin

The West Philippine basin (WPB) is a currently inactive marginal basin belonging to Philippine Sea plate, which has a complex formation history and various crust structures. Based on gravity, magnetic and seismic data, the tectonics in West Philippine basin is characterized by amagma spreading stage and strike slip fractures. NNE trending Okinawa-Luzon fracture zone is a large fracture zone with apparent geomorphology and shows a right-handed movement. The results of joint gravity-magnetic-seismic inversion suggest that the Okinawa-Luzon fracture zone has intensive deformation and is a transform fault. Western existence of the NW trending fractures under Ryukyu Islands Arc is the main cause of the differences between south and north Okinawa Trough. The Urdaneta plateau is not a remained arc, but remnant of mantle plume although its lava chemistry is similar to oceanic island basalt (OIB).

Determination of fractional energy loss of waves in nearshore waters using an improved high-order Boussinesq-type model

Fractional energy losses of waves due to wave breaking when passing over a submerged bar are studied systematically using a modified numerical code that is based on the high-order Boussinesq-type equations. The model is first tested by the additional experimental data, and the model's capability of simulating the wave transformation over both gentle slope and steep slope is demonstrated. Then, the model's breaking index is replaced and tested. The new breaking index, which is optimized from the several breaking indices, is not sensitive to the spatial grid length and includes the bottom slopes. Numerical tests show that the modified model with the new breaking index is more stable and efficient for the shallow-water wave breaking. Finally, the modified model is used to study the fractional energy losses for the regular waves propagating and breaking over a submerged bar. Our results have revealed that how the nonlinearity and the dispersion of the incident waves as well as the dimensionless bar height (normalized by water depth) dominate the fractional energy losses. It is also found that the bar slope (limited to gentle slopes that less than 1:10) and the dimensionless bar length (normalized by incident wave length) have negligible effects on the fractional energy losses.

Impact of spatial resolution and spatial difference accuracy on the performance of Arakawa A-D grids

This paper alms at illustrating the impact of spatial difference scheme and spatial resolution on the performance of Arakawa A-D grids in physical space. Linear shallow water equations are discretized and forecasted on Arakawa A-D grids for 120-minute using the ordinary second-order (M and fourth-order (C4) finite difference schemes with the grid spacing being 100 km, 10 km and I km, respectively. Then the forecasted results are compared with the exact solution, the result indicates that when the grid spacing is I kin, the inertial gravity wave can be simulated on any grid with the same results from C2 scheme or C4 scheme, namely the impact of variable configuration is neglectable; while the inertial gravity wave is simulated with lengthened grid spacing, the effects of different variable configurations are different. However, whether for C2 scheme or for C4 scheme, the RMS is minimal (maximal) on C (D) grid. At the same time it is also shown that when the difference accuracy increases from C2 scheme to C4 scheme, the resulted forecasts do not uniformly decrease, which is validated by the change of the group A velocity relative error from C2 scheme to C4 scheme. Therefore, the impact of the grid spacing is more important than that of the difference accuracy on the performance of Arakawa A-D grid.

Strongly nonlinear long interfacial waves in uniform basic flows and their dispersion relationship

In this paper, long interfacial waves of finite amplitude in uniform basic flows are considered with the assumption that the aspect ratio between wavelength and water depth is small. A new model is derived using the velocities at arbitrary distances from the still water level as the velocity variables instead of the commonly used depth-averaged velocities. This significantly improves the dispersion properties and makes them applicable to a wider range of water depths. Since its derivation requires no assumption on wave amplitude, the model thus can be used to describe waves with arbitrary amplitude.

Ocean wave imaging mechanism by imaging radar

Analytical representations of the high frequency spectra of ocean wave and its variation due to the variation of ocean surface current are derived from the wave-number spectrum balance equation. The ocean surface imaging formulation of real aperture radar (RAR) is given using electromagnetic wave backscattering theory of ocean surface and the modulations of ocean surface winds, currents and their variations to RAR are described. A general representation of the phase modulation induced by the ocean surface motion is derived according to standard synthetic aperture radar (SAR) imaging theory. The detectability of ocean current and sea bottom topography by imaging radar is discussed. The results constitute the theoretical basis for detecting ocean wave fields, ocean surface winds, ocean surface current fields, sea bottom topography, internal wave and so on.

Effect of Stokes drift on upper ocean mixing

Stokes drift is the main source of vertical vorticity in the ocean mixed layer. In the ways of Coriolis - Stokes forcing and Langmuir circulations, Stokes drift can substantially affect the whole mixed layer. A modified Mellor-Yamada 2.5 level turbulence closure model is used to parameterize its effect on upper ocean mixing conventionally. Results show that comparing surface heating with wave breaking, Stokes drift plays the most important role in the entire ocean mixed layer, especially in the subsurface layer. As expected, Stokes drift elevates both the dissipation rate and the turbulence energy in the upper ocean mixing. Also, influence of the surface heating, wave breaking and wind speed on Stokes drift is investigated respectively. Research shows that it is significant and important to assessing the Stokes drift into ocean mixed layer studying. The laboratory observations are supporting numerical experiments quantitatively.

Mixed-layer water oscillations in tropical Pacific for ENSO cycle

The main modes of interannal variabilities of thermocline and sea surface wind stress in the tropical Pacific and their interactions are investigated, which show the following results. (1) The thermocline anomalies in the tropical Pacific have a zonal dipole pattern with 160 W as its axis and a meridional seesaw pattern with 6-8 degrees N as its transverse axis. The meridional oscillation has a phase lag of about 90 to the zonal oscillation, both oscillations get together to form the El Nino/La Nina cycle, which behaves as a mixed layer water oscillates anticlockwise within the tropical Pacific basin between equator and 12 degrees N. (2) There are two main patterns of wind stress anomalies in the tropical Pacific, of which the first component caused by trade wind anomaly is characterized by the zonal wind stress anomalies and its corresponding divergences field in the equatorial Pacific, and the abnormal cross- equatorial flow wind stress and its corresponding divergence field, which has a sign opposite to that of the equatorial region, in the off-equator of the tropical North Pacific, and the second component represents the wind stress anomalies and corresponding divergences caused by the ITCZ anomaly. (3) The trade winds anomaly plays a decisive role in the strength and phase transition of the ENSO cycle, which results in the sea level tilting, provides an initial potential energy to the mixed layer water oscillation, and causes the opposite thermocline displacement between the west side and east side of the equator and also between the equator and 12 degrees N of the North Pacific basin, therefore determines the amplitude and route for ENSO cycle. The ITCZ anomaly has some effects on the phase transition. (4) The thermal anomaly of the tropical western Pacific causes the wind stress anomaly and extends eastward along the equator accompanied with the mixed layer water oscillation in the equatorial Pacific, which causes the trade winds anomaly and produces the anomalous wind stress and the corresponding divergence in favor to conduce the oscillation, which in turn intensifies the oscillation. The coupled system of ocean-atmosphere interactions and the inertia gravity of the mixed layer water oscillation provide together a phase-switching mechanism and interannual memory for the ENSO cycle. In conclusion, the ENSO cycle essentially is an inertial oscillation of the mixed layer water induced by both the trade winds anomaly and the coupled ocean-atmosphere interaction in the tropical Pacific basin between the equator and 12 degrees N. When the force produced by the coupled ocean-atmosphere interaction is larger than or equal to the resistance caused by the mixed layer water oscillation, the oscillation will be stronger or maintain as it is, while when the force is less than the resistance, the oscillation will be weaker, even break.

海洋重力资料预处理软件的编制及应用

首先应用Visual Basic 6.0计算机语言，根据我们所熟知的海洋重力数据的处理步骤和海洋重力测量规范内容，开发出一套海洋重力资料整理软件，以便能使资料整理工作成为一体化，提高工作效率，保证数据处理的准确性。应用该软件对南海某工区的实际重力资料进行了处理，得到了满意的处理结果，证明该软件使用方便，计算准确可靠，能够快速准确的处理采集到的重力数据。对数据处理结果进行了分析和应用。 通过分析Eötvös校正值与重力仪读数之间的变化关系，可以对重力仪的动态性能进行分析和评价。得到2级海况条件下重力仪的实际阻尼延迟时间约为60s，比厂方推荐的理论值（110s）更接近实际，有助于减小数据处理过程中的定位误差，提高测量精确度。通过对重力仪的灵敏性进行分析，验证重力仪在测量过程中的准确性和可靠性，为获得精确的数据提供可靠的保障。 对该套规则重力测网数据，分别采用60s阻尼延迟时间和110s阻尼延迟时间进行阻尼延迟校正，然后进行交点差计算和平差处理，得到的结果显示：采用60s阻尼延迟时间对数据进行阻尼延迟校正后得到的测量精度为1.40mGal，经过平差后的测量精度为0.77mGal；而采用110s阻尼延迟时间对数据进行阻尼延迟校正后得到的测量精度为1.62mGal，经过平差后的测量精度为0.96mGal。此结果验证了利用船只机动转向法确定海洋重力仪的实际阻尼延迟时间的必要性和应用价值。 最后，对宋文尧等人的平差方法和刘保华等人的最小二乘平差方法进行了分析比较，从推导原理和计算结果两方面证明了两种方法在规则海洋重力测网平差作用上是等效的。对于不规则重力测网平差，最小二乘平差方法具有强大的优势，且计算机程序运算实现方便。

西菲律宾海盆西北部重磁资料处理及其地质意义

西菲律宾海盆是西太平洋边缘海盆中最古老的海盆之一，从属于菲律宾海板块，研究表明西菲律宾海盆形成过程复杂。 本文对包括KR03-04航次在内的众多航次的海洋重力及磁力调查数据处理，获得了该区域布格重力异常及磁异常等值线图、上下延拓图、导数异常图等常规成果图件。同时，本文还对本单位在该地区所做的多道地震调查中的三条剖面进行了处理和解释，并对其中的两条测线作了重磁震联合反演研究。通过这些工作，本文得出该地区的构造特征。 该区分为冲大东海岭，乌尔达内塔海台及冲绳－吕宋断裂带等次一级的构造单元，总的构造为NWW走向，也有NE向的大构造带。其中，冲绳－吕宋断裂带破碎严重，NE向走向，在琉球海沟附近(126°E)转为NNE向，为一右旋的深大断裂，可能是由于菲律宾海盆受到琉球岛弧阻拦形成的，也可能是一个转换断层；乌尔达内塔海台中部存在一较大的NW断裂。冲大东海岭和乌尔达内塔海台可能是残留的火山弧；结合ODP及DSDP的相关资料，本文确定了琉球海沟外缘俯冲挤压造山带的形成原因及形成年代。在西菲律宾海盆，伏在一层薄薄的沉积物之下有较厚的一层基岩破碎带，这有可能是西菲律宾海盆无岩浆扩张的有力证据之一，此次扩张的动力可能来源于俯冲于琉球岛弧下面板块的拖曳。

南黄海重磁资料处理与断裂系统研究

南黄海位于中国大陆和朝鲜半岛之间，北与北黄海相通，南与东海相连，是奠基于下扬子地台前震旦纪变质岩基底之上一个多旋回盆地的叠覆。 目前研究表明，南黄海兼跨中朝、扬子和华南3大构造单元，区域构造复杂，是我国近海唯一未获得工业性流油的海域，我国在该区与韩国尚存在海域及油气资源的划界争议问题。对南黄海开展研究，从开展理论研究、开发海洋资源、维护海洋权益等方面出发都有重大意义。由于我国在该区的勘探与研究成效不高，认识较浅，故应加大对该区的研究力度。 本文收集了南黄海重力与磁力资料，对这些资料进行了处理。通过对这些海洋重磁资料的常规处理，本文得到了该地区区域布格重力异常及磁异常的上下延拓图、导数异常图等常规成果图件，对研究区重磁异常进行了分区和特征解释；应用小波多尺度分解的方法对重磁数据进行了分析，并结合常规处理结果，对研究区的断裂系统进行了讨论；初步探讨了南黄海盆地的形成与演化。并获得如下几点结论： 1)研究区中断裂十分发育，优势走向为NE—NEE向，是太平洋板块挤压作用的结果，同时，研究区内不同方向（如NE、NW、EW和NS向等）的断裂相互交错，互有切割，有明显的多期活动的迹象； 2)证实了朝鲜半岛西缘断裂带的存在，其在靠近济州岛时表现的更为破碎；黄海中央断裂带为一条，走向为NNW，向南逐渐变为NW向延伸； 3)南黄海盆地是一个自晚元古代以来不断迁移叠合、并经多次改造而使原貌小完整的序列残留盆地，具有复杂的形成与演化过程。

输电线巡检机器人越障机理与试验

分析超高压输电线路架空地线上的障碍物类型以及跨越这些障碍的过程。可以发现,巡检机器人采用双臂交替跨越障碍,越障过程简单。但是由于受巡检机器人自身重力偏矩的影响和手臂长度尺寸的限制,当单臂悬架在架空地线上时,导致巡检机器人本体倾斜,另一手臂完成脱线和上线任务变得十分困难,有时甚至造成越障失败。为了解决上述问题,提出质量调节的控制方法。该方法通过调节巡检机器人的质心,使巡检机器人的本体保持水平状态。为了验证质量调节控制方法的正确性,采用Lagrange方法建立巡检机器人动力学模型,通过仿真试验、实验室模拟实际架空地线试验以及超高压实际现场试验说明了提出方法的可行性。

大型铰接并联六维测力平台摩擦建模

在已有制造工艺及标定技术基础上,为进一步改善大型铰接并联六维测力平台的测量精度,本文基于螺旋理论和影响系数原理,引入符号函数建立了Stewart结构大型铰接六维测力平台的摩擦模型。文中提出了关节摩擦对铰接并联六维测力平台测量精度的影响矩阵及I、H类误差表达式,绘制了在不同外载和关节摩擦系数条件下,六维测力平台的I、II类误差曲线,并总结丁关节摩擦和平台自重对测力平台测量精度的影响规律。为具有普通球形铰链人型Stewart平台六维测力下台精度的提高和改善提供了理论基础。