青藏高原东南缘构造应力场重构

采用有限元方法,针对青藏高原东南缘建立更细致、更精确的三维有限元弹性模型。选取9种不同的应力边界条件,分别进行优化分析后处理,将对应台站形变模拟值与GPS实测值进行误差分析,最终选取最佳方案作为古构造应力场。结果表明,青藏高原东南缘4 Ma BP的古应力场主要起源于中国大陆周围板块的相互作用,特别是印度板块NNE向强烈碰撞作用,成为中国大陆尤其是西南部青藏高原地区构造应力场最主要的动力来源,控制各个块体相互作用的方式和运动格局。青藏高原东南缘古应力场主要包括几个力源：西北部青藏高原侧向挤压造成的WE向应力约105 MPa;西南部直接来自于印度板块的NE-WS向应力约70 MPa;南部NS向作用力33 MPa;东南部扬子块体侧向NW-SE阻挡力56 MPa;北东部受扬子块体强烈EW向阻挡力90 MPa。这些力源共同作用于青藏高原东南缘,形成现今复杂应力场。
     

郯庐断裂带中南段重磁特征与深部构造

郯庐断裂带是中国东部地球物理场的分界线。它在磁场上呈现出一条NNE向线性正磁异常带,它的重力场特征表现为规模较大的布格重力异常梯度带。在断裂带的两侧,其区域磁场、区域重力场无论是场值大小,还是异常规模以及各异常的展布形态等诸方面都具有明显的区别。地球物理特征揭示郯庐断裂带深部构造,为研究重点矿产分布规律及地震灾害预报具有重要的实际意义。     

Particle trajectories of nonlinear gravity waves in deep water

The paper presents a numerical method for calculating the particle trajectories of nonlinear gravity waves in deep water. Particle trajectories, mass-transport velocity and Lagrangian wave period can be accurately determined by the proposed method. The high success rate of the proposed method is examined by comparing the present results with those of (a) Longuet-Higgins, M.S., 1986, 1987. Eulerian and Lagrangian aspects of surface waves. Journal of Fluid Mechanics 173, 683-707 and (b) Lagrangian moments and mass transport in Stokes waves. Journal of Fluid Mechanics 179, 547-555. It is shown that the dimensionless mass-transport velocity can exceed 10% for large waves, and the Lagrangian wave period is much larger than the Eulerian wave period for large waves.     

Radar cross section analysis of marine targets using a combining method of physical optics/geometric optics and a Monte-Carlo simulation

In this paper, the radar cross section of flat plates on ocean surfaces is statistically investigated. A combining method of physical optics and geometric optics is applied to establish an effective backscattering analysis procedure. This method is a high-frequency analysis method originally derived from a simplified Stratton-Chu integral equation, assuming that the radar is far away from the target so that Kirchhoff approximation is valid. A Monte-Carlo simulation method is adopted to statistically analyze the effects of undulated ocean surfaces. The ocean surfaces are randomly generated by Pierson-Moskowitz ocean wave spectrum and a directional distribution function. Numerical investigations are carried out for flat plates, with the same height and width but with different inclined angles, on ocean surfaces of various significant wave heights.     

Interaction of flexural gravity waves with shear current in shallow water

In the present study, the effect of shear current on the propagation of flexural gravity waves is analyzed under the assumptions of linearized shallow-water theory. Explicit expressions for the reflection and transmission coefficients associated with flexural gravity wave scattering by a step discontinuity in both water depth and current speed are derived. Further, trapping and scattering of flexural gravity waves by a jet-like shear current with a top-hat profile are examined and certain limiting conditions for the waves to exist are derived. The effects of change in water depth, current speed, incident wavelength and the angle of incidence on the group and phase velocities as well as on the reflection and transmission characteristics are analyzed through different numerical results.     

Interaction of water waves with vertical cylinders using null-field integral equations

The scattering of water waves by bottom-mounted vertical circular cylinders is solved by using the null-field integral equations in conjunction with degenerate kernels and Fourier series to avoid calculating the Cauchy and Hadamard principal values. In the implementation, the null-field point can be exactly located on the real boundary owing to the introduction of degenerate kernels for fundamental solutions. An adaptive observer system of polar coordinates is considered to fully employ the properties of degenerate kernels. For the hypersingular equation, vector decomposition for the radial and tangential gradients is carefully considered. This method can be seen as a semi-analytical approach since errors attribute from the truncation of Fourier series. Neither hypersingularity in the Burton and Miller approach nor the CHIEF concept was required to deal with the problem of irregular frequencies. Five advantages of free of calculating principal value, well-posed algebraic system, convergence rate of exponential order, meshfree and elimination of boundary-layer effect, are achieved by using the present approach. Numerical results are given for the forces and free-surface elevation around the circular boundaries. Also, the near-trapped behavior arisen from the physical resonance is detected. A general-purpose program for water wave impinging several circular cylinders with arbitrary number, radii, and positions was developed. Several examples of water wave structure interaction by vertical circular cylinders were demonstrated to see the validity of the present formulation.     

利用TOPEX卫星高度计资料分析台湾岛周边海域巨浪和大风特征

卫星高度计实现了对全球性或区域性的海洋参量的实时监测,TOPEX卫星高度计提供了迄今为止时间序列最长、数据质量最高的全球海面风速和有效波高的同步观测资料。利用TOPEX卫星高度计资料,分析了有效波高4m以上的巨浪在台湾岛周边海域的时空分布特征,结果表明台湾岛周边海域巨浪的分布具有明显的季节特征。每年平均有效波高最大值大多数出现在夏季,春季是1a中有效波高最小的季节,而秋季和冬季是巨浪出现频率最高的季节。波高大于6m的巨浪大都出现在台湾岛东部及东北部海域,在南部海域出现较少。     

Effects of wind waves of the Pacific westerly on the eastern Pacific wave transport

There exists a tongueshaped swelldominance pool known as Swell Pool (SP) in the Eastern Pacific region. The monthlymean wave transports (WT) for each month of 2000 is computed using the wave products of ECMWF reanalysis data. By comparing the 2000 monthlymean WT and monthlymean wind field from QUICKSCAT, large differences are found between the wave transport direction and the wind direction over the Eastern Pacific. This may serve as an evidence for the existence of the SP in this region. The work done in this study indicates that the sources of swell in the Tropical Eastern Pacific (TEP) are in the westerly regions of the Southern and Northern Pacific.     

Nonlinear numerical simulation on extreme-wave kinematics

A fully nonlinear numerical model based on a time-domain higher-order boundary element method (HOBEM) is founded to simulate the kinematics of extreme waves. In the model, the fully nonlinear free surface boundary conditions are satisfied and a semi-mixed Euler-Lagrange method is used to track free surface; a fourth-order Runga-Kutta technique is adopted to refresh the wave elevation and velocity potential on the free surface at each time step; an image Green function is used in the numerical wave tank so that the integrations on the lateral surfaces and bottom are excluded. The extreme waves are generated by the method of wave focusing. The physical experiments are carried out in a wave flume. On the horizontal velocity of the measured point, numerical solutions agree well with experimental results. The characteristics of the nonlinear extreme-wave kinematics and the velocity distribution are studied here.     

Estimation of the wave-related ripple characteristics and induced bed shear stress

A large data set on ripples was collected and examined. A set of new formulas for the prediction of the ripple characteristics is proposed with an emphasis on the disappearance of the ripples. The ripple wavelength was observed to be proportional to the bottom wave excursion but also to be a function of the grain-related Shields parameter and wave period parameter introduced by Mogridge et al. (1994). The ripple steepness was found to be nearly constant for orbital ripples, and with a sharp decrease for suborbital ripples. Two empirical functions are added including the effects of the critical Shields parameters (inception of transport and inception of sheet flow), i.e. giving the boundaries for the ripple existence's domain. The proposed formulas yield better prediction capabilities compared to the previously published formulas, especially when ripples are washed out. The effect of the ripple characteristics on the roughness height and the calculation of the bed shear stress is also discussed. It appeared that the bed shear stress calculation is more sensitive to the empirical coefficient a_r introduced in the estimation of the ripple-induced roughness height or to the limits of existence of the ripples than the ripple characteristics themselves.