Tsunami generation,propagation, and run-up with a high-order Boussinesq model

In this work we extend a high-order Boussinesq-type (finite difference) model, capable of simulating waves out to wavenumber times depth kh < 25, to include a moving sea-bed, for the simulation of earthquake- and landslide-induced tsunamis. The extension is straight forward, requiring only an additional term within the kinematic bottom condition. As first test cases we simulate linear and nonlinear surface waves generated from both positive and negative impulsive bottom movements. The computed results compare well against earlier theoretical, numerical, and experimental values. Additionally, we show that the long-time (fully nonlinear) evolution of waves resulting from an upthrusted bottom can eventually result in true solitary waves, consistent with theoretical predictions. It is stressed, however, that the nonlinearity used far exceeds that typical of geophysical tsunamis in the open ocean. The Boussinesq-type model is then used to simulate numerous tsunami-type events generated from submerged landslides, in both one and two horizontal dimensions. The results again compare well against previous experiments and/or numerical simulations. The new extension compliments recently developed run-up capabilities within this approach, and as demonstrated, the model can therefore treat tsunami events from their initial generation, through their later propagation, and final run-up phases. The developed model is shown to maintain reasonable computational efficiency, and is therefore attractive for the simulation of such events, especially in cases where dispersion is important.     

A temporally second-order accurate Godunov-type scheme for solving the extended Boussinesq equations

A numerical scheme for solving the class of extended Boussinesq equations is presented. Unlike previous schemes, where the governing equations are integrated through time using a fourth-order method, a second-order Godunov-type scheme is used thus saving storage and computational resources. The spatial derivatives are discretised using a combination of finite-volume and finite-difference methods. A fourth-order MUSCL reconstruction technique is used to compute the values at the cell interfaces for use in the local Riemann problems, whilst the bed source and dispersion terms are discretised using centred finite-differences of up to fourth-order accuracy. Numerical results show that the class of extended Boussinesq equations can be accurately solved without the need for a fourth-order time discretisation, thus improving the computational speed of Boussinesq-type numerical models. The numerical scheme has been applied to model a number of standard test cases for the extended Boussinesq equations and comparisons made to physical wave flume experiments.     

礁坪上波浪传播的数值模拟

本文基于具备间断捕捉能力的二阶全非线性Boussinesq数值模型,对规则波和随机波在礁坪地形上的传播变形进行了数值模拟。该模型采用高阶有限体积法和有限差分方法求解守恒格式的控制方程,将波浪破碎视为间断,同时采用静态重构技术处理了海岸动边界问题。重点针对礁坪上波浪传播过程中的波高空间分布和沿程衰减,礁坪上的平均水位变化,以及波浪能量频谱的移动和空间差异等典型水动力现象开展数值计算。将数值结果与实验结果对比,两者吻合情况良好,验证了模型具有良好的稳定性,具备模拟破碎波浪和海-岸动边界的能力,能较为准确地模拟波浪在礁坪地形上的传播过程中发生的各种水动力现象。     

Nonlinear dynamo effects for an inhomogeneously turbulent rotating fluid

Im Rahmen der Magnetohydrodynamik der mittleren Felder werden die turbulenzbedingte elektromotorische Kraft und der turbulenzbedingte Spannungstensor für eine leitende Flüssigkeit mit inhomogener Turbulenz auf einem rotierenden Körper berechnet. Dabei wird der Einflu$sZ des Magnetfeldes auf der Turbulenz berücksichtigt. So erscheinen nichtlineare Effekte bei den Vorgängen, die für das Vorhandensein von Magnetfeld und differentieller Rotation verantwortlich sind; die dafür ma$sZgebenden Koeffizienten sind Funktionen des Magnetfeldes. Die Nichtlinearitäten werden bis zur zweiten Ordnung im mittleren Magnetfeld (d. h. als schwache Nichtlinearitäten) erfa$sZt. Der $aL-Effekt beim $aL$OM-Dynamo erleidet sowohl bei langsamer als auch bei rascher Rotation eine Schwächung durch das Magnetfeld. Es werden einige astrophysikalische Anwendungen der Ergebnisse erörtert. Insbesondere wird auf die Möglichkeit hingewiesen, da$sZ die bei der Sonne beobachteten Torsionsschwingungen als Folge des im Aktivitätszyklus schwankenden magnetischen Einflusses auf den für die differentielle Rotation ma$sZgebenden A-Effekt erscheinen.     

地震波的场方程矩阵和能量的正定二次型及其意义

场方程是能够表述半空间地震波场的整体特征.波动方程、速度方程和能量方程.通过分析可知每个场方程都具有各自的“场方程矩阵”.能量方程能够对所有场方程矩阵进行综合和贯通,给出了能量方程以“弹性矩阵”为核心的普适性表达形式.最后,运用矩阵的正定二次型理论阐述了“能量矩阵与弹性矩阵”之间一致的对称性和正定性.能量矩阵蕴含的动态力的平衡关系、速度的时间_空间分布和能量的传播及变化的物理意义,能够从能量矩阵的正定二次型特性表述出来.本文研究分析问题的方法完全适用于复杂介质模型,相关的认识和结论可以拓展到均匀黏弹性各向同性介质、均匀弹性各向异性介质、均匀黏弹性各向异性介质以及比奥饱和流体介质.     

影响地下应力状态的构造作用

成岩成矿的外部条件中压力是很重要的因素。地下岩石所受的压力不仅来自上覆岩石的重力，而且来自构造力和热力及其它力。这些力合成一个三轴不等的应力状态，相当一个以最短轴为半径的球状应力迭加一个差应力。前者影响成岩成矿和变质作用，后者控制岩石变形和组构的变化。通过三维应变应力测量，这两部分力通过计算是可以区分开的，因而这将是建立地下深度测算新方法的一个基础。     

泉州湾赤潮藻类优势种细胞密度回归方程研究

于2006年5-11月对东海泉州湾赤潮监控区四个监测站位开展赤潮常规监测。根据监测结果,分别以各站位23项水质理化生物环境因子指标为自变量,相应赤潮藻类优势种的细胞密度为因变量,进行多元逐步回归分析,建立了各站位优势种中肋骨条藻(Skeletonema costatum)、太平洋海链藻(Thalassiosira pacifica)、微小原甲藻(Prorocentrum minimum)、尖刺拟菱形藻(Pseudo-nitzschia pungens)、丹麦细柱藻(Leptocylindrus danicus)和旋链角毛藻(Chaetoceros curvisetus)等的细胞密度多参量回归方程。结果表明,所有回归方程的复相关系数都接近于1,方差分析的结果均为回归极显著,表明所建立的回归方程可作为相应赤潮优势种细胞密度预报方程的高度有效性,将对今后泉州湾的赤潮预报提供良好的指导作用。     

Diel and seasonal variability in the fish community structure of a mud‐bottom estuarine habitat in the Gulf of Mexico

We hypothesized that temporal variation in fish species composition and community structure in a low complexity habitat in the Pueblo Viejo Lagoon, Mexico, is influenced by diel light/dark cycles and tidal stage, and by seasonal changes in salinity and temperature. We collected a total of 17,661 individuals during 2‐h interval sampling over six bi‐monthly 24‐h sampling cycles representing 53 species, of which 11 (～20%) were previously unknown in the system. Diel variation indicated that significantly higher numbers of individuals and species were caught at night, whereas diversity and evenness were higher during the day. Species richness was significantly higher in July and January, whereas diversity and evenness peaked around May; both were correlated with temperature. Diel variation in species composition was influenced primarily by the light/dark cycle. Cluster analyses of each diel cycle separated fish assemblages from midday samples from those of nocturnal samples, separated by an extended wide transition period as fish moved at dawn and during the late afternoon/dusk. Significant shifts (as determined by MANOVA) in assemblage structure occurred between months. Canonical correspondence analysis showed that temperature and day/night effects were the most important environmental variables structuring the fish community. This constrained ordination also defined species with specific habitat preferences as follows: (i) diurnal, warm temperature species (mainly planktivores) (Brevoortia gunteri, Cetengraulis edentulus, Diapterus auratus, and Membras martinica); (ii) nocturnal, warm temperature species (mainly predators) (Citharichthys spilopterus, Cathorops melanopus, and Bairdiella spp.); and (iii) low temperature, diurnal species (Brevoortia patronus and Mugil curema) or those with twilight and nocturnal distributions (Anchoa mitchilli, the most numerically abundant species). Our results indicate that diel and seasonal changes in fish community structure were mainly related to day/night cycles and temperature regimes.     

微扰Kepler系统轨道微分方程的近似Lie对称性与近似不变量

把极角θ视为独立变量,得到Kepler系统的轨道微分方程. 首先讨论Kepler系统轨道微分方程的Lie对称性和不变量,微扰Kepler系统轨道微分方程的精确Lie对称性和精确不变量,其次讨论微扰Kepler系统轨道微分方程的近似Lie对称性和近似不变量,并得到了微扰Kepler系统的9个一阶近似Lie对称性和6个一阶近似不变量,其中1个实为精确不变量,而其余5个分别等于微扰系数ε乘以Kepler系统相应的5个不变量.     

New finite-gap solutions for the coupled Burgers equations engendered by the Neumann systems

On the tangent bundle TS^N-1 of the unit sphere S^N-1, this paper reduces the coupled Burgers equations to two Neumann systems by using the nonlinearization of the Lax pair, whose Liouville integrability is displayed in the scheme of the r-matrix technique. Based on the Lax matrix of the Neumann systems, the Abel--Jacobi coordinates are appropriately chosen to straighten out the restricted Neumann flows on the complex torus, from which the new finite-gap solutions expressed by Riemann theta functions for the coupled Burgers equations are given in view of the Jacobi inversion.