Effect of a submerged bay-mouth breakwater on tsunami behavior analyzed by 2D/3D hybrid model simulation

The hybrid numerical model had been developed to simulate a complicated 3D flow around structures generated by tsunami. In the model, the conventional 2D model is adopted for the wide region far from structures and the 3D non-hydrostatic pressure model is used in the limited region adjacent to structures. The applicability of the model is shown by comparisons of the numerical results with the experimental results of the laboratory model tests and the numerical analysis results of the conventional whole 2D simulation. In addition, the effect of a submerged structure at the opening of a breakwater is discussed from the numerical simulations by the hybrid model. The submerged structure improves the stability of the rubble mound and reduces the tsunami inflow into the bay, while it increases the water surface velocity around the opening of the breakwater. The increase of surface velocity causes the increases of impulsive forces by collision with drafts and so on.     

Inferring the mass fraction of floc-deposited mud: application to fine-grained turbidites

Fine sediment deposition in the ocean is complicated by the cohesive nature of muds and their tendency to flocculate. The result is disaggregated inorganic grain size (DIGS) distributions of bottom sediment that are influenced by single‐grain and floc deposition. This study outlines a parametric model that characterizes bottom sediment DIGS distributions. Modelled parameters are then used to infer depositional conditions that account for the regional variation in the grain sizes deposited by turbidity currents on the Laurentian Fan–Sohm Abyssal Plain, offshore south‐eastern Canada. Results indicate that, on the channellized Laurentian Fan, the mass fraction of floc‐deposited mud increases only slightly downslope. The small evolution in this fraction arises because sediment concentration and turbulent energy are associated in turbidity currents. On the Sohm Abyssal Plain, however, the mass fraction of floc‐deposited mud decreases, probably as a result of lower sediment concentration at this source‐distal site. Estimates of the mass fraction of mud deposited as flocs suggest that floc deposition is the dominant mode by which sediment is lost from suspension, although single‐grain deposition contributes more to the depositional flux in proximal areas where high energy breaks flocs and in distal areas where low sediment concentration limits floc formation. It is concluded that, throughout the dispersal system, changes in the fraction of flocculated mud deposited from turbidity currents reflect changes in sediment concentration and energy downslope.     

Formation of a shelfbreak front by freshwater discharge

The circulation and transport of freshwater generated by an idealized buoyant source is studied using a three-dimensional primitive equation model. Freshwater enters the continental shelf, turns anticyclonically and moves downstream in the direction of Kelvin wave propagation. In the region close to the source, the flow reaches an equilibrium in the bottom boundary layer so that freshwater does not spread offshore any further. This offshore equilibrium distance increases as we move downstream until the freshwater is able to feel the presence of the shelfbreak. A shelfbreak front forms and the shelfbreak prevents any further offshore spreading of freshwater in the bottom boundary layer.Two complimentary mechanisms are responsible for the slow cross-shelf migration of freshwater and subsequent trapping of shelfbreak fronts: bottom stress and topographic changes. The shelfbreak creates an active, dynamic process preventing leakage from the continental shelf region to the slope region. However, the dynamical process that traps the front to the shelfbreak is still unclear.The location of the shelfbreak front depends on four dimensionless parameters: scaled inlet volume transport, scaled breadth, scaled “diffusivity” and scaled shelf width. We develop empirical relations for predicting the location of the frontal bottom intersection, given these parameters.     

抵偿高程面任意带高斯投影坐标系统的变形分析及应用

对投影变形问题进行了分析,结合实例,探讨了抵偿高程面任意带高斯投影对控制投影变形的实用性。     

噪声HVSR法在OBS探测中的应用——以珠江口外海OBS台阵为例

噪声谱比法(horizontal-to-vertical spectral ratio, HVSR)通常用于确定沉积层的厚度, 而应用于海底地震仪(OBS)以获取海底沉积厚度的研究较少.本文将HVSR法应用于珠江口外海域OBS台阵中, 验证了该方法用于OBS探测中的可靠性.通过对珠江口外海27台OBS收集的三分量地震噪声记录进行处理, 得到了所有台站的HVSR曲线.有23台OBS的HVSR曲线可获取与沉积层厚度相关的峰值频率, 通过经验公式可以计算得到沉积层厚度.结果显示, 珠江口外海域沉积层(新生代地层)厚度范围为100~3500 m, 研究区东南部沉积层远厚于西北部, 且沉积厚度受到NEE向滨海断裂带控制.本文获得的结果与以往依靠反射地震剖面和钻井等获得的沉积层厚度结果基本一致.进一步分析发现研究区新生代地层具有多层结构, 大部分HVSR曲线在5~12 Hz存在振幅最大的波峰, 认为是与第四系厚度相关的峰值频率, 继而确定第四系厚度为20~50 m.本研究拓展了HVSR法的应用领域, 为获取海底沉积层厚度提供了一种新的可靠方法.

     

西南印度洋中脊三维地震探测中炮点与海底地震仪的位置校正

在海上实施三维地震探测过程中,人工震源枪阵中心与船上GPS的距离及地震探测作业中的船行方向造成炮点实际位置与预设位置有一定偏差;自由落体投放的OBS由于海流的影响会偏离原定设计位置(投放点),因此,炮点与海底地震仪(OBS)的位置校正是三维地震结构研究中的基本环节.本文利用艏向信息校正了炮点位置;采用蒙特卡洛和最小二乘法方法对海底地震仪的位置进行了校正,并探讨了直达水波曲线特征.结果表明 OBS位置一般偏离设计点1 km左右,其误差范围在20 m以内,校正后的OBS记录剖面展示了真实的记录情况.该研究结果为下一步西南印度洋的三维层析成像研究提供了坚实数据基础,同时为今后南海的三维深部地壳结构探测提供经验与借鉴.     

考虑流固耦合的桥墩地震反应方法

跨海桥梁地震反应中,作用于桥墩上的动水压力具有明显的流固耦合特征。依据张量理论,推导时变区域散度变换关系及微分形式的几何守恒律;基于任意拉格朗日一欧拉描述,从采用欧拉描述的流体运动Navier—Stokes方程出发,推导时变区域的流体运动控制方程;给出流固耦合问题中的结构计算模型、耦合面接触条件、耦合场计算方法以及流场网格运动控制方法。以某跨海大桥为例说明桥墩地震反应方法,重点突出地震动输入、流场初始条件模拟等问题。计算结果表明：流固耦合理论能够模拟桥梁墩台地震反应中的流场和结构特性;流场初始条件的正确模拟可保证计算稳定性,并减少运算量;横向地震动激励下,桥墩基底剪力较大,纵向地震动激励下,结构运动剧烈;流固耦合系统中的线弹性结构在地震反应中具有明显的非线性特征。     

基于背景噪声互相关技术的雅浦俯冲带海底地震仪时间校正

由于电磁波无法在水中传播,导致海底地震仪(OBS)与全球定位系统(GPS)无法进行实时对钟,因此会产生一定的时钟漂移.而受船时、海况等影响,被动源OBS在回收后往往电量亏损,致使设备回收时并没有完成对钟,无法直接利用对钟前后的时钟漂移量进行线性校正.本文选取中国科学院海洋研究所2016—2017年间在雅浦俯冲带进行的长周期被动源OBS实验所获得的数据,利用背景噪声互相关方法对OBS进行了时间校正.结果发现,各OBS在经过一年左右时间的布放后,时钟漂移量在1～3 s之间,且基本呈线性趋势,将该结果平均分配到每天从而完成线性校正.另外,探讨了时钟的非线性漂移、误差源、参考台站的选取等因素对时钟准确度的影响,为OBS的数据处理进行了有益的探索.     

Tidal currents,energy flux and bottom boundary layer thickness in the Clyde Sea and North Channel of the Irish Sea

A high-resolution three-dimensional model of the Clyde Sea and the adjacent North Channel of the Irish Sea is used to compute the major diurnal and semidiurnal tides in the region, the associated energy fluxes and thickness of the bottom boundary layer. Initially, the accuracy of the model is assessed by performing a detailed comparison of computed tidal elevations and currents in the region, against an extensive database that exists for the M₂, S₂, N₂, K₁ and O₁ tides. Subsequently, the model is used to compute the tidal energy flux vectors in the region. These show that the major energy flux is confined to the North Channel region, with little energy flux into the Clyde Sea. Comparison with the observed energy flux in the North Channel shows that its across-channel distribution and its magnitude are particularly sensitive to the phase difference between elevation and current. Consequently, small changes in the computed values of these parameters due to slight changes of the order of the uncertainty in the open-boundary values to the model, can significantly influence the computed energy flux. The thickness of the bottom boundary layer in the region is computed using a number of formulations. Depending upon the definition adopted, the empirical coefficient C used to determine its thickness varies over the range 0.1 to 0.3, in good agreement with values found in the literature. In the North Channel, the boundary layer thickness occupies the whole water depth, and hence tidal turbulence produced at the sea bed keeps the region well mixed. In the Clyde Sea, the boundary layer thickness is a small fraction of the depth, and hence the region stratifies.Responsible Editor: Phil Dyke